|Year : 2009 | Volume
| Issue : 3 | Page : 231-236
Drug resistant mutations detected by genotypic drug resistance testing in patients failing therapy in clade C HIV-1 infected individuals from India
AJ Kandathil1, R Kannangai1, VP Verghese2, SA Pulimood3, P Rupali4, G Sridharan1, P Grant5, D Pillay5, OC Abraham4
1 Department of Clinical Virology, Christian Medical College, Vellore 632 004, India
2 Department of Child Health, Christian Medical College, Vellore 632 004, India
3 Department of Dermatology, Christian Medical College, Vellore 632 004, India
4 Department of Internal Medicine, Christian Medical College, Vellore 632 004, India
5 Centre of Virology, University College London, London, United Kingdom
|Date of Submission||22-Jan-2009|
|Date of Acceptance||07-May-2009|
|Date of Web Publication||4-Jul-2009|
O C Abraham
Department of Internal Medicine, Christian Medical College, Vellore 632 004
Source of Support: None, Conflict of Interest: None
Purpose: There has been an increase in the number of individuals administered antiretroviral therapy (ART) in India but treatment outcome is hampered by increasing development of drug resistance. Previous reports from India have shown M184V as the commonest mutation in treated individuals. However, there is no evidence for any protease mutations in these reports. This study was done to observe the common/unique mutational patterns observed in reverse transcriptase (RT) and protease (Pr) genes of clade C HIV-1 strains from individuals showing treatment failure in India. Materials and Methods: The assay was done by sequencing the Pr and RT genes of the HIV-1 strains from 18 individuals failing ART. Analysis was carried out using Stanford HIV drug resistance database (SHDB). The sequences were also submitted to the calibrated population resistance tool of SHDB and Rega HIV-1 sub typing tool. Phylogenetic analysis and quality control were performed with Mega 4. Results: Among the 20 strains, 19 showed resistance to both nucleoside reverse transcriptase inhibitors (NRTIs) and non-nucleoside reverse transcriptase inhibitors (NNRTIs), one strain to NNRTIs and five strains showed protease inhibitors (PI) resistance and 3-class resistance. The most common mutation conferring NRTI resistance was M184V (90%) while K103N (45%) was the most common mutation conferring NNRTI resistance. The M46I mutation was seen in 20% of the Pr sequences. Conclusion: Resistance testing to check the prevalence of drug resistance mutations that arise following failure of the first line regimen to establish guidelines for second line regimens in India is a must. Studies are needed to confirm if mutation patterns that arise among clade C following failure of ART are the same as for clade B strains.
Keywords: Clade C, drug resistance, HIV-1, India, treatment failure
|How to cite this article:|
Kandathil A J, Kannangai R, Verghese V P, Pulimood S A, Rupali P, Sridharan G, Grant P, Pillay D, Abraham O C. Drug resistant mutations detected by genotypic drug resistance testing in patients failing therapy in clade C HIV-1 infected individuals from India. Indian J Med Microbiol 2009;27:231-6
|How to cite this URL:|
Kandathil A J, Kannangai R, Verghese V P, Pulimood S A, Rupali P, Sridharan G, Grant P, Pillay D, Abraham O C. Drug resistant mutations detected by genotypic drug resistance testing in patients failing therapy in clade C HIV-1 infected individuals from India. Indian J Med Microbiol [serial online] 2009 [cited 2015 Aug 29];27:231-6. Available from: http://www.ijmm.org/text.asp?2009/27/3/231/53205
| ~ Introduction|| |
There has been an increase in the number of individuals administered antiretroviral therapy (ART) in India under the auspices of the National AIDS Control Organization (NACO), India. ART is known to reduce morbidity and mortality associated with HIV infection but there are a number of factors that can hamper its benefits, with adherence being one of them.  Adherence is a dynamic process and can vary not only among HIV infected individuals but also in a given individual from time to time.  Individuals with a high level of adherence have a greater chance of viral suppression than individuals who are moderately adherent or non-adherent.  Incomplete suppression can lead to emergence of drug-resistant strains.  Previous reports from India have reported M184V as the most commonly observed mutation in the reverse transcriptase (RT) of HIV-1 strains in patients failing ART. ,
This is a preliminary report based on findings of a study on drug resistant mutations that emerge in the RT and protease (Pr) genes of HIV-1 strains in infected individuals on ART in India. This report would contribute to the global data on the common mutational patterns observed in clade C HIV-1 strains from India available in the drug resistance database like the Stanford HIV database.
| ~ Materials and Methods|| |
Blood samples were collected from 18 HIV-1 infected individuals at the clinical virology department of a tertiary care center in South India. These individuals were referred from the infectious disease clinic of the hospital for CD4+ T cell estimation and HIV-1 viral load testing.
For this study, only individuals with detectable HIV-1 viral loads more than eight months after initiation of ART during the year 2005-2007 were recruited. Informed consent was obtained from all the individuals before collection of the samples. This study was approved by the Institution Review Board (IRB). Among the 18 treatment experienced patients, follow up information were available for two patients. The Guava® Easy CD4™ System (Guava Technologies, Hayward, CA, USA), as described earlier, was used for CD4 estimation.  Viral load estimation was done using Artus RealArt™ HIV-1 kits (Qiagen Hamburg Gmbh, Germany) using Rotor- Gene™ 3000. 
The genotypic resistance assay was carried out by sequencing the Pr and RT genes. HIV-1 RNA was extracted from plasma using QIAamp® viral RNA extraction kit (Qiagen GmbH, Hilden, Germany). Amplification and sequencing of the Pr and RT genes and was carried out as previously described. 
All the sequences were also submitted to the calibrated population resistance tool (CPR) available at the Stanford HIV drug resistance database. Information on the frequency of mutations observed in our strains was generated based on the drug resistance mutation list of the IAS-USA major (2007).
Phylogenetic analysis and quality control of sequences were performed for the 20 sequences with Mega 4 software using a minimum evolution method with Kimura two-parameter and bootstrap value of 500 replicates.  These were compared to the HIV-1 pol sequences obtained from the Los Alamos HIV sequence database. The sequences included were consensus sequences for M group, A1,A2,B,C,D,AE,F1,F2,G,H; Indian subtype C sequences (AF286232, AF067155, AF067157, AY713414, AY049708, DQ826595 and DQ826669), non-Indian subtype C sequences from China (AY967806), Ethiopia (AY255823), Myanmar (AB097871) and Zambia (AF544009 and DQ093604). The M group consensus sequence was used to root the tree. The sequences were also submitted to the Rega HIV-1 sub typing tool. 
| ~ Results|| |
The mean HIV-1 viral load in this group was log 10 5.54 copies/ml with a median value of log 10 4.72 copies/ml. The CD4 counts available on 16 patients ranged from 26 - 1155 cells/mm 3 , with a median value of 306 cells/mm 3 and mean of 313 cells/mm 3 . Among the 18 patients who had an age of range of 11-47 years there were 15 (83%) males and three (17%) females. The study population consisted of individuals from South (n=16, 89%), East (n=1, 5.5%) and North (n=1, 5.5%) India.
In the phylogenetic tree [Figure 1] the study sequences were individually distinct and seen to segregate with consensus C at the major node indicating robust sequencing information. Similar results were obtained using the REGA HIV-1 Sub typing Tool - v2.0.
Among the 20 strains, 19 showed resistance to both nucleoside reverse transcriptase inhibitors (NRTIs) and non-nucleoside reverse transcriptase inhibitors (NNRTIs) while one strain was resistant only to the NNRTIs. Among the 20 strains, five strains were resistant to protease inhibitors (PIs). These five strains showed three-class resistance i.e., they were resistant to NRTIs, NNRTIs and PIs.
The RT gene showed major mutations at 20 sites with 24 mutation patterns while the Pr gene showed 4 mutation patterns at 4 sites. The combination of mutations observed in the RT and PR sequences along with the ART combinations are shown in [Table 1]. Amplification and sequencing of the reverse transcriptase gene revealed eighteen (90%) strains with M184V substitution, a known 3TC resistance-conferring mutation. Among the observed substitutions, K103N, Y181C, Y188L and G190A have implications for NNRTI resistance while M41L, D67N, K70R, T215F and T215Y have been implicated in NRTI resistance.
The most common substitution responsible for resistance to PIs was M46I (n=4) in these strains. The other substitutions observed were D30N and L90M (n=2, each). The prevalence of major amino acid substitutions seen in RT and Pr known to affect susceptibility to the different classes of drugs used in ART are shown in [Table 2]. The difference in the prevalence of the mutations between the study sequences and the Stanford database is shown in [Table 2].
| ~ Discussion|| |
The RT amino acid sequence was found to have more drug resistance-conferring mutations than the Pr amino acid sequence. This could be explained as the outcome of all patients having received at least one regimen containing NRTIs plus NNRTIs. PI mutations were seen in only 4 of 18 patients, as PI use is still limited in India.
Eighteen (90%) of the 19 strains sequenced in this group had the M184V mutation. This mutation is selected for by lamividune (3TC), emtricitabine (FTC) or abacavir (ABC) use and confers high-level resistance to 3TC.  Patients receiving incompletely suppressive 3TC regimens usually develop M184V as their first mutation. 
The other mutations observed included nucleotide excision mutations (NEM) 41, 67, 70, 210, 215 and 219.  The most commonly occurring NEM was at position 215. The T215Y mutation was observed in 25% of the sequences while T215F was seen in 20%. These NEMs cause drastic reduction in susceptibility to azidothymidine (AZT). They do not, however, drastically affect the action of 3TC although they can affect the synergism of 3TC along with AZT.  Only 5% of sequences showed the Q151M mutation. The strain showing the Q151M substitution did not show any of the thymidine analogue mutations (TAMs) except for K70R; this strain also showed secondary mutations at positions 75, 77 and 116. The Q151M mutation is seen in patients receiving didanosine (ddI) in combination with AZT and stavudine (d4T) and is rarely seen in 3TC-containing NRTI regimens. , As most first line ART regimens in India contain 3TC this may explain the low prevalence of Q151M in the study sequences. Previous publications from India have highlighted the absence or low prevalence of the Q151M mutation. , Based on data from the West from previously published studies, M184V was the most commonly observed mutation in patients on NRTIs.  We have very limited data but the mutational patterns observed appear similar to published data. ,, One of the worrying factors is a recent report from western India which documented M41L, D67N, M184V, and A98G in the RT gene of treatment naοve individuals. 
Mutations conferring resistance to NNRTIs were also seen, with one strain showing resistance only to NNRTIs. The NNRTI resistance-conferring mutations included K103N (45%), Y181C (15%) and Y188L (20%). Additionally, Y181I mutation conferring resistance to NNRTIs was seen in 5%. The minor mutations observed in the RT conferring resistance to NNRTIs that were observed included V90I, K101E and V106I. We were unable to discern any identifiable pattern of association between the major and minor mutations. Furthermore, the RT mutations were not restricted to any particular site domain of the RT molecule relative to its function.
Four patients showed resistance to PIs. With the exception of L90M, the other mutations were observed near or in the flap region of the protease as deduced from 3D modeling. The flaps must be flexible in order to allow movement of the substrate and product.  D30N was observed in one of the patients showing failure to PIs. D30N confers resistance only to NFV with no cross-resistance to other PIs.  The patient who had the D30N mutation also had M46I mutation that was, however, not observed in the sequence of the sample sequenced at a second visit two years later. This could be indicative of reduced viral fitness of the strain, indirectly indicated by a reduced viral load with a log 10 0.43 reductions at the patient's second visit. An alternative explanation could be amplification of a different quasispecies in the second sample. The M46I mutation contributes to reduced susceptibility to nelfinavir and fosamprenavir.  The N88D mutation was seen at both visits.
L90M mutation was observed in two sequences and was present along with M46I. L90M has been reported from patients treated with saquinavir, nelfinavir, indinavir and ritonavir.  L90M has been implicated in clinical cross-resistance to each of the PIs.  In India, PIs are not normally part of first line antiretroviral regimens because of their cost. The low prevalence of protease mutations are similar to previously published literature from India among treatment experienced patients. , The L89V mutation was seen in our treatment naοve population. In a recent study from western India among treatment naοve individuals has documented V82A mutation in the protease gene.  Based on the Stanford drug resistance database the mutation contributes to resistance to darunavir (DRV), which is however not used in ART regimens in India. This may possibly explain the lack of selection of this mutation in our treatment-experienced group.
Based on data obtained from the Stanford HIV genotypic drug resistance database all the submitted sequences belonged to clade C. This was also supplemented with phenotypic analysis and results obtained from the REGA HIV-1 sub typing tool. This is consistent with genotypic data for strains from different regions of India including the south that have shown the predominant subtype to be clade C.  Prevalence (%) of major mutations observed in Reverse Transcriptase (1-300 aa positions) and Protease (1-99 aa positions) sequences of ART naοve (n=101) and ART experienced (n=20) HIV-1 infected individuals were harmonized for those observed based on SDRM (2008) and IAS-USA (2007) respectively. The RT amino acid sequence in the treatment-experienced group in our study had a higher prevalence of D67N, K103N, M184V and G190A than those seen in the HIV drug resistance Stanford database. This may be an outcome of the first-line ART regimen (two NRTIs and one NNRTI) followed in India.
Based on our findings, we suggest the use of resistance testing to check the prevalence of the drug resistance mutation that arises following failure of the first line regimen. This will help in establishing guidelines for second line regimens in India. Further studies are needed to confirm if mutation patterns that arise following failure of therapy for clade C are the same for clade B strains.
| ~ Acknowledgement|| |
We would like to thank the Christian Medical College Fluid Research for the partial funding for this study.
| ~ References|| |
|1.||Nischal KC, Khopkar U, Saple DG. Improving adherence to antiretroviral therapy. Indian J Dermatol Venereol Leprol 2005;71:316-20. [PUBMED] |
|2.||Carrieri P, Cailleton V, Le Moing V, Spire B, Dellamonica P, Bouvet E, et al . The dynamic of adherence to highly active antiretroviral therapy: results from the French National APROCO cohort. J Acquir Immune Defic Syndr 2001;28:232-9. |
|3.||Clavel F, Hance AJ. HIV drug resistance. N Engl J Med 2004;350:1023-35. [PUBMED] [FULLTEXT]|
|4.||Deshpande A, Jauvin V, Magnin N, Pinson P, Faure M, Masquelier B, et al . Resistance mutations in subtype C HIV type 1 isolates from Indian patients of Mumbai receiving NRTIs plus NNRTIs and experiencing a treatment failure: resistance to AR. AIDS Res Hum Retroviruses 2007;23:335-40. |
|5.||Sen S, Tripathy SP, Chimanpure VM, Patil AA, Bagul RD, Paranjape RS. Human immunodeficiency virus type 1 drug resistance mutations in peripheral blood mononuclear cell proviral DNA among antiretroviral treatment-naive and treatment-experienced patients from Pune, India. AIDS Res Hum Retroviruses 2007;23:489-97. [PUBMED] [FULLTEXT]|
|6.||Kandathil AJ, Kannangai R, David S, Nithyanandam G, Solomon S, Balakrishnan P, et al . Comparison of Microcapillary Cytometry Technology and Flow Cytometry for CD4+ and CD8+ T-Cell Estimation. Clin Diagn Lab Immunol 2005;12:1006-9. [PUBMED] [FULLTEXT]|
|7.||Kannangai R, Kandathil AJ, Ebenezer DL, Nithyanandam G, Samuel P, Abraham OC, et al . Evidence for lower CD4+ T cell and higher viral load in asymptomatic HIV-1 infected individuals of India (south): Implications for therapy initiation. Indian J Med Microbiol 2008;26:217-21. [PUBMED] |
|8.||Kandathil AJ, Kannangai R, Abraham OC, Priscilla R, Pulimood SA, Verghese VP, et al . The frequency of HIV-I drug resistance mutations among treatment naοve individuals at a tertiary care centre in south India. Int J STD and AIDS- (in press) |
|9.||Tamura K, Dudley J, Nei M, Kumar S. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 2007;24:1596-9. |
|10.||de Oliveira T, Deforche K, Cassol S, Salminen M, Paraskevis D, Seebregts C, et al . An automated genotyping system for analysis of HIV-1 and other microbial sequences. Bioinformatics 2005;21:3797-800. [PUBMED] [FULLTEXT]|
|11.||Geretti AM. Clinical implications of HIV drug resistance to nucleoside and nucleotide reverse transcriptase inhibitors. AIDS Rev 2006;8:210-20. [PUBMED] |
|12.||Shafer RW. Genotypic testing for human immunodeficiency virus type 1 drug resistance. Clin Microbiol Rev 2002;15:247-77. [PUBMED] [FULLTEXT]|
|13.||Van Vaerenbergh K, Van Laethem K, Albert J, Boucher CA, Clotet B, Floridia M, et al . Prevalence and characteristics of multinucleoside-resistant human immunodeficiency virus type 1 among European patients receiving combinations of nucleoside analogues. Antimicrob Agents Chemother 2000;44:2109-17. [PUBMED] [FULLTEXT]|
|14.||Coakley EP, Gillis JM, Hammer SM. Phenotypic and genotypic resistance patterns of HIV-1 isolates derived from individuals treated with didanosine and stavudine. AIDS 2000;14:F9-15. [PUBMED] [FULLTEXT]|
|15.||Lall M, Gupta RM, Sen S, Kapila K, Tripathy SP, Paranjape RS. Profile of primary resistance in HIV-1-infected treatment-naive individuals from Western India. AIDS Res Hum Retroviruses 2008;24:987-90. [PUBMED] [FULLTEXT]|
|16.||Shao W, Everitt L, Manchester M, Loeb DD, Hutchison CA 3rd, Swanstrom R. Sequence requirements of the HIV-1 protease flap region determined by saturation mutagenesis and kinetic analysis of flap mutants. Proc Natl Acad Sci U S A 1997;94:2243-8. [PUBMED] [FULLTEXT]|
|17.||Martinez-Cajas JL, Wainberg MA. Protease inhibitor resistance in HIV-infected patients: Molecular and clinical perspectives. Antiviral Res 2007;76:203-21. [PUBMED] [FULLTEXT]|
|18.||Sen S, Tripathy SP, Patil AA, Chimanpure VM, Paranjape RS. High prevalence of human immunodeficiency virus type 1 drug resistance mutations in antiretroviral treatment-experienced patients from Pune, India. AIDS Res Hum Retroviruses 2007;23:1303-8. [PUBMED] [FULLTEXT]|
|19.||Kandathil AJ, Kannangai R, Abraham OC, Sudarsanam TD, Pulimood SA, Sridharan G. Genotypic resistance profile of HIV-1 protease gene: A preliminary report from Vellore, south India. Indian J Med Microbiol 2008;26:151-4. [PUBMED] |
|20.||Kandathil AJ, Ramalingam S, Kannangai R, David S, Sridharan G. Molecular epidemiology of HIV. Indian J Med Res 2005;121:333-44. [PUBMED] [FULLTEXT]|
[Table 1], [Table 2]
|This article has been cited by|
||Viremia and HIV-1 Drug Resistance Mutations Among Patients Receiving Second-Line Highly Active Antiretroviral Therapy in Chennai, Southern India
| || |
| ||Clinical Infectious Diseases. 2012; |
||Selection of nonnucleoside reverse transcriptase inhibitor-associated mutations in HIV-1 subtype C: evidence of etravirine cross-resistance :
| ||Ujjwal Neogi, Anita Shet, Ranjani Shamsundar, Maria L Ekstrand |
| ||AIDS. 2011; 25(8): 1123 |
||Drug resistance mutations to protease and reverse transcriptase inhibitors in treatment naive HIV-1 clade c infected individuals from Mumbai, India
| ||Potdar, P.D., Daswani, B.S., Rane, N.J. |
| || International Journal of Virology. 2011; 7(1): 13-23 |
||Evolutionary Dynamics of Complex Networks of HIV Drug-Resistant Strains: The Case of San Francisco
| ||R. J. Smith, J. T. Okano, J. S. Kahn, E. N. Bodine, S. Blower |
| ||Science. 2010; 327(5966): 697-701 |
||Analysis of RT Sequences of Subtype C HIV-Type 1 Isolates from Indian Patients at Failure of a First-Line Treatment According to Clinical and/or Immunological WHO Guidelines
| ||Alaka Deshpande,Anne Cécile Jeannot,Marie Hélène Schrive,Linda Wittkop,Patricia Pinson,Hervé J. Fleury |
| ||AIDS Research and Human Retroviruses. 2010; 26(3): 343 |
||Viremia and drug resistance among HIV-1 patients on antiretroviral treatment: a cross-sectional study in Soweto, South Africa :
| ||Ziad El-Khatib, Anna Mia Ekstrom, Johanna Ledwaba, Lerato Mohapi, Fatima Laher, Alan Karstaedt, Salome Charalambous, Max Petzold, David Katzenstein, Lynn Morris |
| ||AIDS. 2010; 24(11): 1679 |
||Analysis of RT sequences of subtype C HIV-type 1 isolates from indian patients at failure of a first-line treatment according to clinical and/or immunological WHO guidelines
| ||Deshpande, A., Jeannot, A.C., Schrive, M.H., Wittkop, L., Pinson, P., Fleury, H.J. |
| ||AIDS Research and Human Retroviruses. 2010; 26(3): 343-350 |