Share this article on:

HIV-1 Reverse Transcriptase Nucleotide Substitutions in Subtype C–Infected, Drug-Naive, and Treatment-Experienced Patients in South India

Saravanan, Shanmugam PhD*; Madhavan, Vidya PhD*; Balakrishnan, Pachamuthu PhD*; Solomon, Sunil S. PhD*; Kumarasamy, Nagalingeswaran PhD*; Mayer, Kenneth H. MD; Waldrop, Greer BA*; Solomon, Suniti MD*; Smith, Davey M. MD‡,§

JAIDS Journal of Acquired Immune Deficiency Syndromes: November 1st, 2011 - Volume 58 - Issue 3 - p e94–e95
doi: 10.1097/QAI.0b013e318232a13c
Letters to the Editor

*Y. R. Gaitonde Centre for AIDS Research and Education, Voluntary Health Services Hospital, Chennai, India

Department of Medicine, Brown University, Providence, RI

Division of Infectious Diseases, University of California San Diego, La Jolla, CA.

§Veterans Affairs Healthcare System, San Diego, CA

To the Editors:

HIV-1 can evolve to escape from antiretroviral therapy (ART) and host immune responses.1–3 The pattern of nucleotide changes in these evolutionary scenarios is dependent on the genetic background of the virus. Given the large genetic diversity of HIV-1 worldwide, the development of ART resistance is most likely different between subtypes. We studied the mutational patterns occurring during virologic failure of ART among patients infected with HIV-1 subtype C in India.

Two groups of individuals infected with HIV-1 subtype C in southern India were investigated. One group (n = 58) was ART naive with viral loads >3000 HIV RNA copies per milliliter (COBAS Amplicor HIV-1 Monitor, Version 1.5; Roche Molecular Systems Inc, Pleasanton, CA), and the other group (n = 524) was receiving first-line ART that consists of 2 nucleoside reverse transcriptase inhibitors and a nonnucleoside reverse transcriptase inhibitors, including combinations of zidovudine/stavudine plus lamivudine plus nevirapine/efavirenz and whose viral load was >3000 HIV RNA copies per milliliter. HIV RNA was extracted from blood plasma using Qiagen RNA MiniAmp kit (Qiagen, Valencia, CA), and bulk sequences of HIV-1 reverse transcriptase (630 base pairs) were generated, as described previously in a validated in-house genotyping assay.4 Sequences were aligned (ClustalX) to an Indian subtype C reference (C.IN.AFo67155) and examined for HIV-1 subtype (REGA v2),5 nucleotide diversity, mismatches, transition and transversion mutations (Highlighter HIV LANL),6 average pairwise distance (TN69 model in HyPHY),7 positive selection (REL in HyPHY),7 drug resistance–associated mutations (IAS-USA and Stanford HIV Resistance Database),8,9 and APOBEC signatures G → A hypermutation (Hypermut HIV LANL).10 Comparisons were made using 1-way analysis of variance and Tukey multiple posttests (Prism v5.03).

Hypermutation was not identified in either group, and overall nucleotide diversity was higher in the treated group (average pairwise distance 0.043 vs. 0.078). Compared with the ART-naive group, the treated group had higher mean number of G to A (14.6 vs. 15.9), A to G (13.91 vs. 16.02), and G to C (0.59 vs. 1.48) single nucleotide transitions and the dinucleotide GT to AT substitution (4.09 vs. 4.76) (all P < 0.001, Table 1). Overall, the A to G transition was the predominant pattern observed in resistance-associated mutations in the treated group. Specifically, the M184V (ATG to GTG) mutation was present in 78.7% of the treated group and was characterized by an A to G transition. Other observed resistance-associated mutations included 37% K103N (AAA to AAC), 33% Y181C (TAT to TGT), 33% D67N (GAC to AAC), 32% G190A (GGC to GCC), 30% with M41L (ATG to TTG), 27% T215Y (ACC to TAC), 23% K70R (AAG to AGG), 18% K101E (AAA to GAA), and 13% V106M (GTG to ATG) (Figure 1). These data are consistent with previous reports, which found that the most prevalent thymidine analogue mutations in HIV-1 subtype C from Africa included D67N and T215Y. Also, consistent with previous reports,3,11,12 we did not observe high rates of the M184I substitution, which is caused by a G to A hypermutation and commonly occurs before the selection and fixation of the M184V mutation. As expected, there was a trend for these resistance sites to demonstrate positive selection in the ART-treated vs. ART-naive groups (P = 0.06), and the only drug resistance–associated mutation in the ART-naive group was the K219Q and it was found in only 1 individual.





The diversity of HIV-1 subtypes is a challenge for prevention and treatment programs worldwide. The differences among HIV-1 subtypes may have a profound impact on clinical management and surveillance of drug resistance, particularly as ART is expanded to patients infected with non–subtype B HIV-1. As expected, in this large cohort of patients infected with HIV-1 subtype C, virologic failure during ART was associated with high rates of resistance. Although this study does not allow for a direct comparison between subtype C vs. B rate of drug resistance development, this is an important issue that should be investigated in well-designed clinical trials that allow for uniform monitoring of ART failure. Interestingly, virologic failure during ART was also associated with greater numbers of nucleotide changes, but many patterns of these changes were not specific to the observed resistance-associated mutations. These patterns likely represent the mutational predilection of the replicating viral population while under ART pressure. The clinical ramifications of these predilections remain unclear and should be investigated in the future.

Back to Top | Article Outline


1. Preston BD, Dougherty JP. Mechanisms of retroviral mutation. Trends Microbiol. 1996;4:16–21
2. Quinones-Mateu M, Weber J, Rangel H, Charkaborty B. HIV-1 fitness and antiretroviral drug resistance. AIDS Rev. 2001;3:223–242
3. Rezende LF, Drosopoulos WC, Prasad VR. The influence of 3TC resistance mutation M184I on the fidelity and error specificity of human immunodeficiency virus type 1 reverse transcriptase. Nucleic Acids Res. 1998;26:3066–3072
4. Saravanan S, Vidya M, Balakrishanan P, et al. Evaluation of two human immunodeficiency virus-1 genotyping systems: ViroSeq 2.0 and an in-house method. J Virol Methods. 2009;159:211–216
5. de Oliveira T, Deforche K, Cassol S, et al. An automated genotyping system for analysis of HIV-1 and other microbial sequences. Bioinfomatics. 2005;21:3797–3800
6. . Los Alamos HIV Databases & HIV Database Tools Available at: December 5, 2010
7. Pond SK, Muse SV. Site-to-site variation of synonymous substitution rates. Mol Biol Evol. 2005;22:2375–2385
8. Johnson VA, Brun-Vézinet F, Clotet B, et al. Update of the drug resistance mutations in HIV-1: December 2010. Top HIV Med. 2010;18:156–163
9. Liu TF, Shafer RW. Web resources for HIV type 1 genotypic-resistance test interpretation. Clin Infect Dis. 2006;42:1608–1618
10. Rose PP, Korber BT. Detecting hypermutations in viral sequences with an emphasis on G → A hypermutation. Bioinformatics. 2000;16:400–401
11. Keulen W, Boucher C, Berkhout B. Nucleotide substitution patterns can predict the requirements for drug-resistance of HIV-1 proteins. Antiviral Res. 1996;31:45–57
12. Novitsky V, Wester CW, DeGruttola V, et al. The reverse transcriptase 67N 70R 215Y genotype is the predominant TAM pathway associated with virologic failure among HIV type 1C-infected adults treated with ZDV/ddI-containing HAART in southern Africa. AIDS Res Hum Retrovir. 2007;23:868–878
© 2011 Lippincott Williams & Wilkins, Inc.