BRIEF REPORT
Detection of Low-Level K65R Variants in Nucleoside Reverse Transcriptase Inhibitor–Naive Chronic and Acute HIV-1 Subtype C Infections Jin-fen Li,1 Jonathan T. Lipscomb,1 Xierong Wei,1 Neil A. Martinson,2 Lynn Morris,3 Walid Heneine,1 and Jeffrey A. Johnson1 1Division of HIV/AIDS Prevention, Centers for Disease Control and Prevention, Atlanta, Georgia; 2Perinatal HIV Research Unit, University of the Witwatersrand and Johns Hopkins University School of Medicine, Baltimore Maryland; and 3AIDS Unit, National Institute for Communicable Diseases, Johannesburg, South Africa
To substantiate reports of greater emergence of the K65R nucleoside reverse transcriptase inhibitor (NRTI) mutation in human immunodeficiency virus type 1 (HIV-1) subtype C, we examined natural low-level K65R expression in subtype C relative to subtypes B and AE. We used allele-specific polymerase chain reaction to screen HIV-1 amplified by reversetranscription high-fidelity polymerase chain reaction from subtype C–infected South African women and infants and CRF01(subtype AE) from Thailand; all subjects were NRTI naive. We found low-level K65R of unknown clinical significance in NRTI-naive subtype C–infected women and infants at frequencies above the natural occurrence in subtypes B and AE. The frequent appearance of subtype C frameshift deletions at codon 65 supports a propensity for transcription error in this region. A global scale-up of antiretroviral therapy (ART) and increasing interest in antiretroviral drugs for use in prophylaxis to prevent human immunodeficiency virus (HIV) transmission has raised
Received 23 August 2010; accepted 1 November 2010. Potential conflicts of interest: J.A.J. and W.H. are authors on a patent application for the real-time polymerase chain reaction assays. All other authors report no conflicts of interest. Disclaimer: The findings and conclusions in this manuscript are those of the authors and do not necessarily represent the official view of the Centers for Disease Control and Prevention (CDC). The use of trade names is for information purposes only and does not constitute CDC endorsement. Presented in part: Conference on Retroviruses and Opportunistic Infections in San Francisco, California, 16-19 February, 2010. Reprints or correspondence: Jeffrey A. Johnson, MD, 1600 Clifton Rd NE, MS G45, Atlanta, Georgia 30333 (
[email protected]). The Journal of Infectious Diseases 2011;203:798–802 Published by Oxford University Press on behalf of the Infectious Diseases Society of America 2011. 1537-6613/2011/2036-0001$15.00 DOI: 10.1093/infdis/jiq126
798
d
JID 2011:203 (15 March)
d
Li et al.
concern over whether subtype variation may differentially select for drug resistance. A few reports have claimed that subtype C, the predominant HIV-1 clade worldwide, may select resistance to nonnucleoside reverse transcriptase (RT) inhibitors [1] and nucleoside/nucleotide RT inhibitors (NRTIs) (reviewed in [2]) more easily than other HIV subtypes. The drug resistance mutation K65R has received particular attention, because it confers resistance to tenofovir and moderate resistance to lamivudine/ emtricitabine, NRTIs much favored for their high potency and tolerability. In the absence of drug pressure, for a region of the HIV genome to yield a frequency of error higher than what would be expected from random transcriptional inaccuracy, an element of the virus genome may itself compel introduction of the mutations. Possible genetic features that might explain this include nucleotide sequences or folding that causes the polymerase to slip and pause on the template [3]. If indeed particular HIV-1 subtype sequences can more frequently generate RT missteps at sites associated with drug resistance, then drug resistance selection might be greater in ART-experienced areas where those subtypes are prevalent. We previously evaluated subtype B clinical specimens from the era before the introduction of ART and found that the K65R mutation was not detectable at frequencies ..3% [4]. To obtain evidence of natural K65R quasispecies expression in HIV-1 subtype C at frequencies ..3%, we used sensitive real-time polymerase chain reaction (PCR) assays to screen subtype C–infected NRTI-naive women and their infants who became infected. Mother-infant pairs were further assessed for evidence of K65R vertical transmission. To compare subtype C with an additional subtype, we also screened subtype AE specimens for K65R above the frequency seen with subtype B. Our findings shed light on the differential natural expression of K65R in subtype C HIV-1. METHODS Study Populations
To assess the natural frequency of the K65R mutation, we evaluated masked plasma RNA samples from (1) 99 South African ART-naive, chronically subtype C–infected pregnant women before they received intrapartum single-dose nevirapine and (2) 20 8-week-old infected infants born to transmitting mothers and only given single-dose nevirapine at birth [5, 6]. The women were participants in mother-to-child prevention studies in 2001–2005. In 5 infants, a plasma sample was also collected 10 days after birth (n 5 25 infant samples). The
subtype C results were compared with those from 112 randomized subtype AE–infected, ART-naive participants enrolled in a Thai open-label trial in 1999–2000 (EPV40001; baseline viral loads, 104–106 copies/mL) and with our historical data on 138 pre–ART era subtype B clinical samples collected in 1982–1985 [4]. The Thai specimens were previously available samples kindly provided by ViiV/GlaxoSmithKline for drug resistance assay evaluation. The South African and Thai studies were approved by local institutional review boards, and the subjects provided informed consent for storage of blood for future HIV testing. The institutional review board of the Centers for Disease Control and Prevention approved this study after determining that the recoded, archived specimens were not linked to personally identifiable subjects.
CT values of ,9 and ,11 PCR cycles for the subtype C and subtype AE assays, respectively, corresponded to ..3% K65R. Real-time PCR amplicons from reactions that were positive for K65R were sequenced to confirm the presence of intact mutant sequences, as described elsewhere [4, 7]. In-house bulk genotyping was also performed to determine whether K65R was detectable by conventional testing. Clonal Confirmation
Clonal sequencing was subsequently performed on specimens that screened positive for K65R by real-time PCR to estimate mutant frequencies; 90–135 clonal sequences of 190 base pairs (RT codons 56–119) per positive specimen were compared with the sample bulk genotype and the K65R-positive real-time amplicon sequences to elucidate mutant subpopulations.
Reverse-Transcription PCR
HIV genomic RNA was extracted from patient blood plasma using the Qiagen UltraSens Virus Kit or Qiagen MagAttract Viral RNA M48 kit (Qiagen M48 BioRobot). The templates for real-time PCR testing were generated by reverse-transcription PCR with murine leukemia virus RT and GeneAmp High Fidelity PCR polymerase (Applied Biosystems), using primers that spanned codons 51–236 of the HIV-B/AE RT [4] and codons 20–272 of the HIV-C RT, as detailed elsewhere [5]. K65R-Specific Real-Time PCR Screening
All samples were screened using the Bio-Rad iCycler plus iQ4 optical unit and real-time PCR–based drug resistance assays [4] for the K65R mutation. The assay discriminates the second position of RT codon 65 for a guanosine transition (A/G), which is indicative of the arginine mutation. The oligonucleotide sequences for the subtype AE K65R assay were as follows: forward primer, HIV-AE_K65R.1F (5#-ATAYAATACTCCARTATTTGCTATAAACAG); reverse primer, HIV-B comREV [4]; FAMlabeled probes, AE Com 3.1P (5#-FAM-TCAGTAACAG‘‘T’’ACTAGATGTGGGAGATGCATAT, 80%) and AE Com 3.2P (5#FAM-TCAGTAACAG‘‘T’’ACTGGATGTGGGGGATGCATAT, 20%). The K65R assay oligonucleotides for subtype C were as follows: forward primers, 65R.6F (5#-CAATACTCCAGTATT TGTCATACCAAG, 80%) and HIV-C 65R.5.1F (5#-AACACTCCARTATTTGCYATACCAAG, 20%); reverse primer, HIV-C 65.1REV (5#-TYTTTAACCCTGMTG GGTGTGGTAT); FAMlabeled probes, HIV-C 65.1P (5-FAM-TCAGGGARC‘‘T’’C AATAAAAGAACTCAAGACTTYTGGGA, 80%) and HIV-C 65.2P (5-FAM-TCAGG GAAC‘‘T’’YAAYAAAAGAACTCAAGACTTYTGGGA, 20%). Quotation marks indicate the nucleotide to which the Black Hole Quencher is attached. K65R-specific and total copy amplification cycles that cross the fluorescence thresholds were compared (DCT) to extrapolate the mutant frequency, as previously described. K65R in subtypes AE and C was screened for above the .3% threshold we had previously established as the quasispecies background level for subtype B [4]. Delta
RESULTS Subtype C K65R Assay Screening
K65R screening of the 99 South African antepartum samples yielded positive results for 6 women (6.0%) (Table 1). The K65Rpositive real-time amplicon sequences revealed that for 2 women (SA1-m40 and SA2-m23) the K65R mutation was in frame, as evidenced by the primer 3’-terminus unambiguously recognizing a guanosine in the second position of codon 65 (Supplemental Figure S1). For 2 other women (SA1-m15 and SA1-m32), the mutation-specific amplicon sequence revealed that the primer was correctly binding a second-position guanosine at codon 65, but also a third-position guanosine in the second position, the latter a result of an adenosine missing in codon 65 (21 frameshift product). For the remaining 2 women with positive screening results for K65R (SA2-m16 and SA2-m19), the amplicon sequence revealed that the assay detected only thirdbase guanosine frameshifts in codon 65, suggesting that no intact K65R was present. Of the 99 women evaluated, 20 had HIV-1–infected infants who received only single-dose nevirapine at birth [7]. Screening of the infant plasma obtained at day 10 and week 8 after birth identified 3 infants (SA2-i6, SA2-i19, and SA2i21) who were positive for K65R at preliminary screening (Table 1). K65R-specific amplicon sequencing predicted that all 3 infants had intact K65R variants as well as the 21 frameshift, suggested by mixed sequences at codon 65. The median DCTs for mothers and infants were 11.1 and 10.4 cycles, respectively. Subtype AE K65R Assay Screening. Sensitive screening of the 112 ART-naive Thai subjects identified 1 individual (.9%) who was positive for K65R, with a mutation frequency estimated at .5% (DCT, 10.5 cycles). The mutation amplicon sequence demonstrated an intact mutation and showed no downstream resistance mutations linked to the K65R. Low-Level K65R in Subtype C HIV-1
d
JID 2011:203 (15 March)
d
799
Table 1. Results of Clonal Analysis in Samples Positive for K65R at Preliminary Real-Time Polymerase Chain Reaction Screening
Screening DCT, PCR cyclesa
No. of Intact K65R Clones (%)
SA1-m15
8.7
1/135 (.7)
SA1-m32
8.8
3/90 (3.3)
SA1-m40 SA2-m16
8.8 6.5
1/90 (1.1) 0/135 (0)b
SA2-m19
8.4
0/135 (0)b
SA2-m23
8.9
Subject Subtype C–infected mothers
1/90 (1.1)
Subtype C–infected infants SA2-i6 (day 10)
8.4
6/90 (6.7)
SA2-i19 (week 8)
8.4
3/90 (3.3)
SA2-i21 (day 10)
8.7
2/90 (2.2)
SA2-i21 (week 8) SA2-i23 (day 10)
9.0 9.8c
1/90 (1.1) 0/135 (0)
Subtype AE–infected subject Th_28
10.5
1/90 (1.1)
a
NOTE. The K65R DCT cutoff was ,9 cycles for subtype C and ,11 cycles for subtype AE. b
Two of the 3 infants with intact K65R, SA2-i6 and SA2-i21, also had the nonnucleoside RT inhibitor mutations, Y181C and K103N 1 Y181C, respectively, from the single-dose nevirapine exposure, a finding not known to be connected to the appearance of K65R, as indicated in the Stanford HIV Drug Resistance Database RT Treatment Profile for NVP (http://hivdb. stanford. edu/cgi-bin/RTMutSummary.cgi). Clonal analysis of the subtype AE–infected subject who had a positive screening result for K65R identified the mutation intact in 1 of 90 clones (1.1%), with no other mutations linked in this K65R clone. We further screened this individual with sensitive testing for other thymidine analog mutations (M41L, K70R, T215 F/Y), and randomly sequenced 45 additional clones from this subject, and could not find evidence of other mutations to indicate transmitted drug resistance or any frameshifts at codon 65. Phylogenetic analysis of K65R variants to exclude low-level cross-contamination revealed that intact K65R clones from infants formed discrete branches with their source bulk sequences and were distinct from the variants identified in mothers (MEGA 4.1 software; data not shown).
Only frameshift variants were identified.
c
Above the cutoff but screened to verify absence of K65R vertical transmission from the infant’s mother (SA2-m23).
Bulk sequencing revealed no detectable mutations in the subtype AE and C samples that were positive using real-time PCR screening. Clonal Screening
K65R clonal screening confirmed intact K65R mutants in 4 (4%) of 99 ART-naive, subtype C–infected South African women (SA1-m15, SA1-m32, SA1-m40, and SA2-m23) at frequencies of .7%–3.3% (Table 1). Clonal sequences surrounding codon 65 are shown in Figure 1. Two women who had intact K65R (SA1-m15 and SA1-m32) also had clones with an adenosine missing from codon 65. The other 2 women (SA2-m16 and m19) who were initially positive at screening had only codon 65 frameshift clones, which were consistent with their screening assay amplicon sequences. For all 3 HIV-1–infected infants who were screening assay positive (SA2-i6, SA2-i19, and SA2-i21), clonal analysis confirmed intact K65R mutants at frequencies of 2.2%–6.7% (3/20 infants overall; 15%). One infant, SA2-i21, had intact K65R at both the 10-day and 8-week time points. Furthermore, the 3 K65R-positive infants all had codon 65 frameshift mutants. One clone from SA2-i6 had 2 nucleotides missing from codon 65 and 2 missing from codon 66 (24 frameshift). In no matched pair was intact K65R identified in both mother and infant. The difference in prevalence of K65R between infants (15%) and mothers (4%) did not achieve statistical significance (P 5 .057; Fishers exact test). 800
d
JID 2011:203 (15 March)
d
Li et al.
DISCUSSION We describe the natural appearance of low-level K65R variants in South African NRTI-naive adults and infants with subtype C infections at frequencies higher than those detected for subtype B and AE quasispecies. The absence of evidence for motherinfant K65R transmission supports a spontaneous emergence of K65R that was nearly 4-fold more prevalent in acutely infected infants than in chronically infected women. However, incomplete infant viral load data prevented us from ascertaining whether greater viremia, typically associated with rampant, acute infection, could have contributed to increased K65R expression. The most striking finding with subtype C viruses was the frequent appearance of frameshift mutations at codon 65, which suggested that the region is a mutagenesis ‘‘hot spot.’’ This transcription error has recently been demonstrated in vitro for subtype C and was suggested to be the possible result of sequence-induced RT slippage and pausing, leading to nucleotide dislocation [8, 9]. An earlier evaluation of the same mothers using an M184Vspecific assay found that the M184V quasispecies frequency was very similar to that of subtype B (data not shown), further indicating that the extent of mutations observed around codon 65 was not a result of indiscriminate mutagenesis in subtype C. The only K65R variant identified in the subtype AE population could not be explained by the same dislocation mechanism reported for subtype C [8, 9], because there was no adjacent guanosine in codon 65, and, furthermore, the absence of codon 65 frameshifts also implied a less-volatile genomic region.
Figure 1. Subtype C (SA1 and SA2) and AE (Thai) sample sequences surrounding reverse transcriptase codon 65; arrow identifies the second position of codon 65 where A/G (boxed) indicates a K65R mutation (AGG, AGA). Bulk sequences are followed by clones for each subject. Dashes represent nucleotide deletions; def, defective frameshift clones that have a downstream G shifted into the second position of codon 65. A subtype C reference sequence (ref1) is shown at the top, with the amino acid frames translated and codon 65 in boldface.
Although additional screening did not identify intact K65R clones in mothers of infants with K65R, we did not screen the number of clones that could have compensated for random selection error. Therefore, it is possible that intact variants were present above the target frequency in those mothers but were missed in the random clone analyses. Furthermore, the years sampled for each subtype varied, and unrecognized thymidine analog resistance may have been
transmitted in the contemporary subtype C and AE populations. However, the absence of other mutations provided no indication that our results may have been biased by cases of transmitted drug resistance. Although recent reports have indicated that PCR generates spurious low-level K65R with subtype C [10], we found intact variants in infants at frequencies above the reported artifact frequencies. Moreover, if the observed codon 65 errors are indeed template Low-Level K65R in Subtype C HIV-1
d
JID 2011:203 (15 March)
d
801
driven, they would be expected to occur with lower-fidelity HIV RT in vivo. The data suggest that subtype C generates a higher quasispecies frequency of codon 65 mutations. This finding may help explain what has been reported as greater selection of K65R in subtype C–infected subjects receiving ART, mostly from the use of stavudine [2, 11–13]. However, longitudinal studies of K65R at baseline and during therapy are needed to ascribe clinical significance to low-level quasispecies variants. The detection of subtype-specific spontaneous frameshift variants may confound resistance screening; thus, data from point mutation testing must be analyzed carefully. Funding This work was supported entirely by intramural funds from the Centers for Disease Control and Prevention.
References 1. Eshleman SH, Hoover DR, Chen S, et al. Nevirapine (NVP) resistance in women with HIV-1 subtype C, compared with subtypes A and D, after the administration of single-dose NVP. J Infect Dis 2005; 192:30–6. 2. Martinez-Cajas JL, Pai NP, Klein MB, Wainberg MA. Differences in resistance mutations among HIV-1 non-subtype B infections: a systematic review of evidence (1996-2008). J Int AIDS Soc 2009; 12:11. 3. Coustinos D, Invernizzi CF, Moisi D, et al. The K65R mutation in subtype B and C HIV-1: rates of development and the implications of template specific dislocation mutagenesis. [abstract]. Antivir Ther 2010; (Suppl 15):A13.
802
d
JID 2011:203 (15 March)
d
Li et al.
4. Johnson JA, Li J-F, Wei X, et al. Simple PCR assays improve the sensitivity of HIV-1 subtype B drug resistance testing and allow linking of resistance mutations. PLoS One 2007; 2:e638. 5. Johnson JA, Li J-F, Morris L, et al. Emergence of drug-resistant HIV-1 after intrapartum administration of single-dose nevirapine is substantially underestimated. J Infect Dis 2005; 192:16–23. 6. Martinson NA, Morris L, Johnson J, et al. Women exposed to singledose nevirapine in successive pregnancies: effectiveness and nonnucleoside reverse transcriptase inhibitor resistance. AIDS 2009; 23:809–16. 7. Johnson JA, Li J-F, Wei X, et al. Minority HIV-1 drug resistance mutations are present in antiretroviral treatment-naı¨ve populations and associate with reduced treatment efficacy. PLoS Med 2008; 5:e158. 8. Coutsinos D, Invernizzi CF, Xu H, et al. 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 2009; 83:2029–33. 9. Coutsinos D, Invernizzi CF, Xu H, Brenner BG, Wainberg MA. Factors affecting template usage in the development of K65R resistance in subtype C variants of HIV type-1. Antivir Chem Chemother 2010; 20:117–31. 10. Varghese V, Wang E, Babrzadeh F, et al. Nucleic acid template and the risk of a PCR-induced HIV-1 drug resistance mutation. PLoS One 2010; 5:e10992. 11. Turner D, Shahar E, Katchman E, et al. Prevalence of the K65R resistance reverse transcriptase mutation in different HIV-1 subtypes in Israel. J Med Virol 2009; 81:1509–12. 12. Pillay V, Pillay C, Kantor R, et al. HIV type 1 subtype C drug resistance among pediatric and adult South African patients failing, C drug resistance among pediatric and adult South African patients failing antiretroviral therapy. AIDS Res Hum Retroviruses 2008; 24:1449–54. 13. Hosseinipour MC, van Oosterhout JJ, Weigel R, et al. The public health approach to identify antiretroviral therapy failure: high-level nucleoside reverse transcriptase inhibitor resistance among Malawians failing first-line antiretroviral therapy. AIDS 2009; 23:1127–34.
