Figures
Abstract
Background
We and others have shown that subtype C HIV-1 isolates from patients failing on a regimen containing stavudine (d4T) or zidovudine (AZT) exhibit thymidine-associated mutations (TAMs) and K65R which can impair the efficacy of Tenofovir (TDF) at second line. Depending on the various studies, the prevalence of K65R substitution as determined by the Sanger method ranges from 4 to 30%. Our aim was to determine whether ultra-deep pyrosequencing (UDPS) could provide more information than the Sanger method about selection of K65R in this population of patients.
Methods
27 subtype C HIV-1 isolates from treated patients failing on a regimen with d4T or AZT plus lamivudine (3TC) plus nevirapine (NVP) or efavirenz (EFV) and who had been sequenced by Sanger were investigated by UDPS at codon 65 of the reverse transcriptase (RT). 18 isolates from naïve patients and dilutions of a control K65R plasmid were analysed by Sanger plus UDPS.
Results
Analysis of Sanger sequences of subtype C HIV-1 isolates from naïve patients exhibited expected polymorphic substitutions compared to subtype B but no drug resistance mutations (DRMs). Quantitation of K65R variants by UDPS ranged from <0.4% to 3.08%. Sanger sequences of viral isolates from patients at failure of d4T or AZT plus 3TC plus NVP or EFV showed numerous DRMs to nucleoside reverse transcriptase inhibitors (NRTIs) including M184V, thymidine-associated mutations (TAMs) plus DRMs to non- nucleoside reverse transcriptase inhibitors (NNRTIs). Two K65R were observed by Sanger in this series of 27 samples with UDPS percentages of 27 and 87%. Other samples without K65R by Sanger exhibited quantities of K65R variants ranging from <0.4% to 0.80%, which were below the values observed in isolates from naïve patients.
Citation: Recordon-Pinson P, Papuchon J, Reigadas S, Deshpande A, Fleury H (2012) K65R in Subtype C HIV-1 Isolates from Patients Failing on a First-Line Regimen Including d4T or AZT: Comparison of Sanger and UDP Sequencing Data. PLoS ONE 7(5): e36549. https://doi.org/10.1371/journal.pone.0036549
Editor: Luis Menéndez-Arias, Centro de Biología Molecular Severo Ochoa (CSIC-UAM), Spain
Received: February 3, 2012; Accepted: April 9, 2012; Published: May 16, 2012
Copyright: © 2012 Recordon-Pinson et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This study was funded by Merck Sharp and Dohme-Chibret (grant 3Z22204938). JP and SR were supported by Agence Nationale de Recherches sur le Sida et les Hépatites Virales.The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have received funding from Merck Sharp and Dohme-Chibret for this study. There are no patents, products in development or marketed products to declare. This does not alter the authors' adherence to all the PLoS ONE policies on sharing data and materials.
Introduction
We and others [1]–[2] have shown that subtype C HIV-1 isolates from Indian patients who fail on first-line HAART composed of stavudine (d4T) or zidovudine (AZT) plus lamivudine (3TC) plus nevirapine (NVP) or efavirenz (EFV) and according to WHO clinical and/or immunological criteria exhibit numerous drug resistance mutations (DRMs) to nucleoside reverse transcriptase inhibitors (NRTIs) and to non-nucleoside reverse transcriptase inhibitors (NNRTIs) in the reverse transcriptase (RT) part of the viral genome, including thymidine-associated mutations (TAMs) and K65R (the prevalence of which was around 8% in our series). When the RT sequences were introduced into the ANRS and Stanford algorithms, both algorithms showed that the DRMs of the first line induce a decreased susceptibility to tenofovir (TDF), an NRTI drug that is still used as second line in some southern countries. This has implications for public health because patients who fail with a first-line regimen including d4T or AZT plus 3TC plus NVP or EFV and who switch to 3TC plus TDF plus ritonavir-boosted lopinavir (LPV/RTV) will in fact not be fully susceptible to TDF and therefore to the second-line regimen. TAMs and K65R are known to induce partial or full resistance to TDF [3]. Regarding K65R, similar studies carried out in the same context of failure on a first-line regimen including d4T, AZT or dideoxyinosine (ddI) showed a prevalence of 4% in a South African population where subtype C was predominant [4], 10.9% in CRF02_AG plus G viruses in Nigeria [5], 14% in subtype C isolates from South Africa [6], 24% in Malawi with subtype C viruses [7] and up to 30% in Botswana [8].
From a molecular point of view, it has been demonstrated that the RT KKK nucleotide motif at codons 64, 65, 66 in reverse transcriptase of subtype C HIV-1 appears to lead to template pausing that facilitates the selection of K65R, even in isolates from untreated patients [9]–[12]. Moreover, it has also been shown that the KKK motif in this subtype can lead to PCR-induced K65R [13].
The aim of the present study was to clarify the prevalence of K65R in these subtype C isolates from patients failing on a first line including d4T or AZT. Since we had the prevalence of K65R by the Sanger sequencing method, we investigated K65R variants by ultradeep pyrosequencing (UDPS) in the same samples as those already used for bulk sequencing and sought whether UDPS provided additional knowledge to the Sanger results.
Methods
Subtype C HIV-1 isolates from Indian patients failing (according to WHO clinical and/or immunological criteria) on a first-line treatment including d4T or AZT plus 3TC plus NVP or EFV were sequenced on RT by the Sanger method, and the sequences were recorded in the Los Alamos database (GenBank JF895621–JF895673). Among these samples, 27 were randomly selected for investigation by UDPS using the Roche GS Junior equipment. RNA extracted previously from the samples was used to amplify a short region of RT with primers including specific sequences for the GS Junior system (Table 1). The reverse transcription used SuperScriptIII RTPCR enzyme (Invitrogen, Carlsbad, CA) with 10 ul RNA, one cDNA synthesis cycle at 50°C for 30 min and 40 cycles of PCR amplification. The nested PCR used FastStart HiFi (Roche) with 2 ul of RTPCR product, 40 cycles of PCR amplification. Amplicons were purified by AMPure kit (Agencourt Biosciences), quantified using Quant-iT Picogreen (Invitrogen) and pooled at equimolar concentrations. Clonal amplification on beads (EmPCR) was performed using the 454 Life Science reagents that enable bidirectional sequencing, composed of a 30 cycle PCR amplification. DNA containing beads were recovered and UDPS was performed on the GS Junior sequencer (454 Life Sciences). Most of the samples had HIV viral loads >100,000 copies/mL (mean: 379,753 copies/mL; IQR: 11209–5,817,977 copies/mL) and at least 1,000 clonal sequencing reads were used for the analysis, allowing a 0.4% accuracy in the quantitation [14]. Positions studied were codon 65 and codon 70, which was chosen as a TAM position in NRTI-treated patients leading to K70R, because it has been frequently observed in the studied isolates and because 70R is considered to be a DRM and not at all a substitution potentially related to polymorphism.
Since the viral isolates obtained at initiation of first-line therapy were not available, the data obtained by Sanger and UDPS in isolates from patients at failure were compared to those of subtype C isolates from 18 naïve patients. Some of their bulk sequences have been previously published by our group [15]. All codons were analysed by Sanger while potential polymorphism at codons 65 and 70 was investigated by both Sanger and UDPS. As a control, we used two subtype C MJ4 plasmids, one wild type and one bearing K65R (both provided by Mark Wainberg's group in Montreal). The UDPS results of the study are available in GenBank under accession number SRA 050640.
Results
Table 2 shows the Sanger results of viral isolates from naïve patients and UDPS results of codons 65 plus 70 in these isolates, together with UDPS results of control plasmids for codons 65 and 70. Analysis of Sanger results for naïve patients showed an extensive polymorphism compared to subtype B without involvement of substitutions 65R and 70R. K70R was <0.4% by UDPS in all isolates from naïve patients. K65R ranged from <0.4% to 3.08% (mean 0.66±0.76 standard deviation, SD).
Regarding the 27 isolates from treated patients at failure, the drug resistance mutations (DRMs) according to the French ANRS algorithm and following bulk DNA sequencing (Sanger) are listed in Table 3. Most of them exhibited the M184V mutations to 3TC of the regimen plus TAMs to d4T and/or AZT and DRMs to NNRTIs. Only two samples (455 and 493) bore a K65R mutation, one (455) cumulating K65R and the Q151M nucleoside analog mutation (NAM). Table 4 compares the Sanger and UDPS data of these isolates at codons 65 and 70. Quantitation of K65R by UDPS ranged from <0.4% to 87%. In the two isolates with K65R by Sanger, percentages of K65R by UDPS were 27% and 87%. If we only consider positions 65 found not to have K65R with Sanger (25 samples), the quantities of K65R ranged from <0.4% to 0.80% (0.162±0.22 SD). K70R with UDPS ranged from <0.4% to 100%. Eight samples exhibited K70R with Sanger and the corresponding UDPS values ranged from 34.50 to 100%. Regarding positions 70 found not to have the K70R mutation with Sanger (19), all of them were <0.4% for 70R by UDPS except two samples (470 and 489 with values of 1.80 and 4.40% respectively).
Discussion
Sanger and UDPS results in isolates from naïve patients
Although UDPS has limitations particularly with regard to polymerization and pyrosequencing errors [13], [16], recent studies with different methods (UDPS, allele specific PCR) have shown that K65R is identified more frequently in subtype C HIV-1 from naïve patients [14], [17]. In our naïve patients, there was a clear difference between K70R (mean 0%) and K65R (mean 0.66%) (Table 2). Our data on position 65 are in agreement with those of Kozal et al [14], ranging from <0.4% to 1.33% apart from one sample (1131) at 3.08%. We were not expecting selection of K65R in this population of patients who were quite distant from primary infection and therefore from potential transmission of K65R mutants. For us, this naïve population was the basis of an evaluation of natural variation at codon 65 in subtype C HIV-1.
Sanger and UDPS results of isolates from treated patients at failure
The Sanger results of isolates from treated patients were as expected with a predominance of M184V, numerous TAMs of pathway1 (M41L, D67N, K70R, L210W, T215Y/F) and DRMs to NNRTIs (mainly K101E, K103N, V106M, Y181C, G190A). As mentioned above, 2 isolates of the series exhibited a K65R substitution.
With regard to UDPS results at codon 70, the quantitative data were different from those recorded in naïve patients: 8 isolates exhibiting K70R with Sanger had UDPS K70R values above 34.60%, while 2 samples without K70R with Sanger (470 and 489) had K70R variants at quantities above the <0.4% background observed in naïve patients. We hypothesize that these isolates are undergoing a process of selecting K70R mutations. Regarding the K65R values apart the two samples with K65R by Sanger, the UDPS quantities of K65R variants were low and below those of isolates from naïve patients. Our results are not in accordance with those obtained by another group [18] using an allele- specific PCR which exhibited minority variants of K65R in four subtype C HIV-1 isolates out of 30 patients having received NRTIs at first line; it must be pointed out that this technique uses an intercalating dye and high-melt resolution assay which can be difficult to interpret due to genomic variability in the flanking region of codon 65.
K65R substitutions are generated in subtype C isolates from naïve patients due to the 64–65–66 motif. There are some constraints in experienced patients failing on a suboptimal regimen with d4T or AZT plus 3TC plus NVP or EFV. We first hypothesize that 184V, which was the most prevalent mutation observed in our series of treated patients at failure, has dampened the emergence of 65R as noted by others [3]. Second, there is an antagonism between TAMs and K65R, while the latter can be found in association with NAMs (Q151M) and is considered to be increasingly selected in the presence of DRMs to NNRTIs and particularly Y181C and G190A. As noted above, the prevalence of K65R in this clinical context of failure ranges from 4 to 30%. In our series, we estimate this prevalence to be 8% and the UDPS data did not reveal any process of K65R selection that cannot be assessed by Sanger sequencing.
Acknowledgments
We thank M. Wainberg's group in McGill University, Montreal, for providing the control K65R plasmid. We thank Roche for providing the UDPS GS Junior equipment. The authors are indebted to Ray Cooke for linguistic assistance.
Author Contributions
Conceived and designed the experiments: PP HF. Analyzed the data: PP SR HF. Contributed reagents/materials/analysis tools: PP JP AD. Wrote the paper: PP SR HF.
References
- 1. Deshpande A, Jeannot AC, Schrive MH, Wittkop L, Pinson P, et al. (2010) 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. AIDS Res Hum Retroviruses 3: 343–350.
- 2. Kumarasamy N, Madhavan K, Venkatesh K, Saravanan S, Kantor R, et al. (2009) High frequency of clinically significant mutations after first-line generic highly active antiretroviral therapy failure: implications for second-line options in resource-limited settings. Clin Infect Dis 49: 306–309.
- 3. Brenner BG, Coutsinos D (2009) The K65R mutation in HIV-1 reverse transcriptase: genetic barriers, resistance profile and clinical implications. HIV Ther 3: 583–594.
- 4. van Zyl GU, van der Merwe L, Claassen M, Zeier M, Preiser W (2011) Antiretroviral resistance patterns and factors associated with resistance in adult patients failing NNRTI-based regimens in the Western Cape, South Africa. J Med Virol 83: 1764–1769.
- 5. Hawkins CA, Chaplin B, Idoko J, Ekong E, Adewole I, et al. (2009) Clinical and genotypic findings in HIV-infected patients with the K65R mutation failing first-line antiretroviral therapy in Nigeria. J Acquir Immune Defic Syndr 52: 228–234.
- 6. Pillay V, Pillay C, Kantor R, Venter F, Levin L, et al. (2008) HIV type 1 subtype C drug resistance among pediatric and adult South African patients failing antiretroviral therapy. AIDS Res Hum Retroviruses 11: 1449–1454.
- 7. Hosseinipour MC, van Oosterhout JJ, Weigel R, Phiri S, et al. (2009) The public health approach to identify antiretroviral therapy failure: high-level nucleoside reverse transcriptase inhibitor resistance among Malawians failing first-line antiretroviral therapy. J AIDS Jun 1;23(9): 1127–34.
- 8. Avalos A, Brenner B, Gaolathe T, Mine M, Gaseitsiwe S, et al. (2006) High prevalence of the K65R mutation in human immunodeficiency virus type 1 subtype C isolates from infected patients in Botswana treated with didanosine-based regimens. Antimicrob Agents Chemother 50: 4182–4185.
- 9. Coutsinos D, Invernizzi CF, Xu H, Moisi D, Oliveira M, et al. (2009) Template usage is responsible for the preferential acquisition of the K65R reverse transcriptase mutation in subtype C variants of human immunodeficiency virus type 1. J Virol 83: 2029–2033.
- 10. Coutsinos D, Invernizzi CF, Xu H, Brenner BG, Wainberg MA (2010) Factors affecting template usage in the development of K65R resistance in subtype C variants of HIV type-1. Antivir Chem Chemother 20: 117–131.
- 11. Coutsinos D, Invernizzi CF, Moisi D, Oliveira M, Martinez-Cajas JL, et al. (2011) A template-dependent dislocation mechanism potentiates K65R reverse transcriptase mutation development in subtype C variants of HIV-1. PLoS One 6: e20208.
- 12. Invernizzi CF, Coutsinos D, Oliveira M, Moisi D, Brenner BG, et al. (2009) Signature nucleotide polymorphisms at positions 64 and 65 in reverse transcriptase favor the selection of the K65R resistance mutation in HIV-1 subtype C. J Infect Dis 200: 1202–1206.
- 13. Varghese V, Wang E, Babrzadeh F, Bachmann MH, Shahriar R, et al. (2010) Nucleic acid template and the risk of a PCR-Induced HIV-1 drug resistance mutation. PLoS One 5: e10992.
- 14. Kozal MJ, Chiarella J, St John EP, Moreno EA, Simen BB, et al. (2011) Prevalence of low-level HIV-1 variants with reverse transcriptase mutation K65R and the effect of antiretroviral drug exposure on variant levels. Antivir Ther 16: 925–929.
- 15. Deshpande A, Karki S, Recordon-Pinson P, Fleury HJ (2011) Drug Resistance Mutations in HIV Type 1 Isolates from Naive Patients Eligible for First Line Antiretroviral Therapy in JJ Hospital, Mumbai, India. AIDS Res Hum Retroviruses. AIDS Res Hum Retroviruses 27: 1345–1347.
- 16. Shafer RW (2009) Low-abundance drug-resistant HIV-1 variants: finding significance in an era of abundant diagnostic and therapeutic options. J Infect Dis 199: 610–612.
- 17. Li JF, Lipscomb JT, Wei X, Martinson NA, Morris L, et al. (2011) Detection of low-level K65R variants in nucleoside reverse transcriptase inhibitor-naive chronic and acute HIV-1 subtype C infections. J Infect Dis 203: 798–802.
- 18. Toni TA, Brenner BG, Asahchop EL, Ntemgwa M, Moisi D, et al. (2010) Development of an allele-specific PCR for detection of the K65R resistance mutation in patients infected with subtype C human immunodeficiency virus type 1. Antimicrob Agents Chemother 54: 907–911.