- Short report
- Open Access
Highly diversified multiply drug-resistant HIV-1 quasispecies in PBMCs: a case report
Retrovirologyvolume 5, Article number: 43 (2008)
Although drug resistance is a major challenge in HIV therapy, the effect of drug resistance mutations on HIV evolution in vivo is not well understood. We have now investigated genetic heterogeneity in HIV-1 by performing drug resistance genotyping of the PR-RT regions of viruses derived from plasma and peripheral blood mononuclear cells (PBMCs) of a single patient who had failed multiple regimens of anti-retroviral therapy.
Patterns of drug resistance mutations showed that the viral populations in PBMCs were more heterogeneous than in plasma. Extensive analysis of HIV from infected PBMCs in this patient showed that high-level diversity existed among 109 cloned PR-RT sequences and that the majority of mutations were related to drug resistance. Moreover, the PBMCs included archival species that reflected the treatment history of the patient while those in plasma were mainly related to the most recent treatment. Some of the proviral clones contained single or multiple mutations in various combinations. Approximately eighteen percent of the proviral clones derived from infected PBMCs were defective, i.e. 5.5% contained single nucleotide deletions (frameshift mutations) and 12.8% encoded in-frame stop codons (nonsense mutations). Amino acid substitutions in PR and the polymerase region of RT occurred in 12–15% of cases but were much less frequent in the RNase H region of RT, which might not have been under drug selection pressure.
Selective drug pressure can yield multiple drug-resistant quasispecies that include archival and replication-incompetent species in PBMC reservoirs.
HIV quasispecies within infected individuals can rapidly adapt to hosts [1–7] due, in part, to variations in replicative fitness that enable some viruses to grow faster than others[3, 8]. This is of obvious clinical relevance, since viral genetic changes can result in alterations in receptor usage, escape from drug and host immune pressure, and can impact on viral pathogenesis. HIV-1 may also evolve separately in different physiological compartments, e.g., peripheral blood mononuclear cells (PBMCs) vs. the central nervous system. Here, we report on an individual who failed multiple antiviral therapies (ART), including use of nucleoside and non-nucleoside RT inhibitors (NRTIs and NNRTIs) and protease inhibitors (PIs). After initiating therapy elsewhere with undisclosed regimens, the patient was treated in 1999 at the Jewish General Hospital, Montreal, Canada, with zidovudine (AZT)/lamivudine (3TC)/efavirenz (EFV) plus unboosted indinavir (IDV) and nelfinavir (NFV) for 9 months and 3 months, respectively, and was switched to stavudine (d4T)/3TC/amprenavir (APV) for 12 months, at which time viral samples were obtained for resistance testing. Viral RNA from plasma and proviral DNA from PBMCs were purified using commercial kits (Qiagen, Mississauga, ON, Canada). Initial HIV-1 genotyping was performed using Trugene HIV-1 genotyping kits (Siemens Diagnostics Inc., Toronto, Canada). All studied were performed with approval of the Ethics Review Committee, Jewish General Hospital.
The degree of quasispecies heterogeneity was higher in PBMCs than in plasma
Mutations in PR and RT associated with drug resistance were compared in plasma vs PBMCs. Both types of samples contained viruses with multiple primary (M46I/L, G48V, I54V, V82A or L90M) and secondary resistance mutations (e.g. L10I) in PR as well as multiple mutations in RT (M41L, E44A, T69N, V118I, M184V, L210W, T215Y, K219R for NRTIs) (A98G, K101E, V1081, Y181C and G190A for NNRTIs) (Table 1). Both the plasma and PBMC samples contained mixtures of mutations, although some mutational motifs were only detected in the PBMCs. For example, mixtures of 41K/R, 54I/V, 64I/V, 82V/A, 90M/I in PR and 181Y/C, 190A/G, 219K/R in RT were identified in PBMCs but not in plasma. Conversely, 35D and 69N in PR and 108I in RT were detected only in plasma but not PBMCs, as determined by genotyping. These results were confirmed by clonal sequencing of PBMC DNA. In general, viruses harbouring the unboosted protease motif including L90M were exclusively present in PBMCs. This is consistent with the fact that genotyping often fails to detect minority species that are represented at levels <10 to 35% in a given population .
Resistance-associated mutations in PR-RT clones reveal heterogeneous viral populations within infected PBMCs
Viral genetic diversity in the infected PBMCs was analyzed by randomly selecting and sequencing 109 clones of two independent cloning efforts. Nested PCR was performed to amplify the entire PR-RT region. One pair of primers, forward 5'-ACTGAGAGACAGGCTAATTTTTTAGG and backward 5'-TTGGGCCTTATCTATTTCCAT (Bio S&T, Montreal, Canada) was used for the first round of PCR using Taq polymerase (Invitrogen, Burlington, ON, Canada) with 30 cycles of annealing at 55°C for 1 min and extension at 72°C for 3 min. A second round of 25 cycle PCR with primers (forward 5'-ACTATCCATGGTCCCTCAGATCACTCTTTGG and backward 5'-ACTAATTTGTCGACTTGTTCATTTCCTCC (Bio S&T)) was used to generate a 2.1 kb DNA fragment spanning the PR and RT genes. The sample was diluted 100 fold in the second round PCR. The fragment was cloned into the Nco 1-Sal I sites of a modified version of vector pTWIN2 (New England Biolabs, Toronto, Canada) by standard molecular cloning methods. Positive clones were amplified and sequenced using ABI fluorescence sequencing kits (Applied Biosystems, Foster City, CA, USA). Analysis of DNA sequences and protein translation were performed with GeneTool software (BioTools Inc., Alberta, Canada). All viral nucleotide positions were assigned based on differences from the HxB2 reference consensus sequence.
No wild-type (wt) PR or RT sequences were detected in any of these clones, each of which contained a minimum of 8 and a maximum of 12 mutations in PR and a minimum of 11 or a maximum of 19 mutations in RT. A total of 15 mutations were present at resistance-associated sites in PR among the 109 clones tested, while mixtures existed at 9 of these sites. In regard to RT, a total of 19 resistance mutations were detected, of which 11 were also present within mixtures. The frequency of resistance motifs in PBMCs are shown in Table 2. Thymidine analog mutations (TAMs, M41L, L210W, T215Y) were present in all clones while 184V was detected in 81% of sequences. Two reasons may explain the difference. First, the patient had used thymidine analogues for a longer period than 3TC. Second, viruses carrying TAMs are more fit than those with M184V, and the latter mutation may have been under selection pressure even in the presence of 3TC.
ART apparently influenced the emergence of multiple resistance-associated quasispecies. Several PR secondary resistance mutations were found in all of the quasispecies, e.g. 10L, 16E, 37T, 48V, 60E, 63P, 77I. Other mutations were found in only a portion of the PR sequences analyzed, i.e. 35D, 41K, 46I/L, 54V 69Y, 82A and 90M, compared with the HxB2 reference sequence, while the other viruses sequenced were wt in regard to these positions. The 82A and 90M mutations in PR were not found together in any of the clones; neither were 35D and 41K in PR, reflecting the fact that the switch from unboosted to boosted PR inhibitors in the regimen selected for the 82A mutational pathway as opposed to 90M, that represented an archival species present in PBMCs rather than plasma. Resistance-associated substitutions in RT were also found in multiple combinations. 41L, 44A, 98G, 118I, 208I, 210W, 211R, 214F, 215Y, and 283I were found in all the clones (Tables 2 &3). K101E/V108I and Y181C/G190A were mostly present in paired fashion. The frequencies of resistance-associated single substitutions that were only detected in some quasispecies within RT were 74V, 108I, 181C, 184V, and 219R. Several sites contained two or three different amino acid substitutions, i.e., 37T/37S, 46I/46L/46S, 63P/63H within PR and 44A/44T, 69N/69A, 10lE/101Q, 102Q/102R, 190A/190T, 210W/210R/210C, 211D/211E within RT (Table 2). These mixtures may reflect a change from AZT-to d4T-containing regimens. PR had the highest variations of all the sequences analyzed. The average pairwise nucleotide distances of the viral quasispecies within the PBMCs of this patient were 0.022(0.024) for PR, 0.015 (0.012) for the N-terminal 220 amino acids of RT, 0.011(0.007) for the second part of RT (amino acids 221–440) and 0.015(0.013) for RNase H by Mega package (version 4.0) . Numbers in brackets were obtained using the first and second nucleotides of each codon.
Non resistance-associated substitutions in PR-RT
Variations also occurred at 4 sites within PR that are not recognized as conferring drug resistance, i.e. P39S, I62V, I64V and I72V. Of these, I64V might be pre-supposed to be of biological significance as a secondary or compensatory mutation, since it was found to cluster together with the 46L, 54V and 82A mutations associated with APV but not 46I and 90M associated with the unboosted IDV and NFV regimens, that were present in 45.8% of the 109 clones analyzed. In the case of RT, prevalent substitutions included V35I, K43E, E122K, S162D, K166R, D177E, E203K, K223Q, L228R, I293V, P294Q, I326V, Q334E, P345Q, A355G, R356K I375V, T376A, T377Q and T400A. A158P, R206S and G359S also appeared at relatively low frequency (i.e. 6%–17%) and in a dispersed pattern, i.e. not as a cluster. By way of comparison, we also analyzed the RNase H region of RT and found that variability was slight, with 3 substitutions present among 90 amino acids at its N terminus in all clones, i.e. L452S, N460D, and N519S (519S is, in fact, the consensus amino acid among subtype B viruses although N519 is present within HxB2). Various frequencies were detected for V466M (24.8%), V467I (20.2%), Y483H (54.1%) and N514S (22%); these substitutions, in fact, are naturally present among various HIV-l subtypes and were probably not selected in the patient. V466M was found to be dispersed yet clustered with Y483H. V4671, Y483H and N514S did not co-exist within any of the clones that were studied. Thus, only 3 to 6 substitutions were identified within RNase H, suggesting that viral evolution had occurred more slowly in this region than in the similarly sized PR region, reflecting differences in drug selection pressure. Although RT inhibitors have been reported to affect RNase H, no mutations in RNase H were identified in our recent report .
Frameshift and nonsense mutations
HIV-1 RT may commonly introduce addition/deletion and/or stop codons into the HIV-1 genome that result, in turn, in an accumulation of defective quasispecies. Our results suggest that defective viral particles may have represented a substantial proportion of viruses within the patient whose samples were analyzed. Of the 109 clones, 6 (5.5%) contained single nucleotide deletions; all of these were located at or near the poly A and C tracts at which RT is believed to be especially error-prone. In-frame nonsense mutations were found at greater frequency, i.e. in 3.7% of cases in PR and in 9.2% of cases in RT. About 18% of the HIV-1 clones analyzed appeared to be defective in PR and/or RT and would be expected to yield non-infectious particles. Defective quasispecies constitute a high proportion of HIV-l particles in infected individuals. Previous reports showed that sequences in uncultured vs cultured PBMCs can be different. 18% of the sequences analyzed in our study were defective in either PR or RT, similar to other reports that examined the vif gene. Since the size of the HIV-1 genome is greater than PR-RT, we would expect that the proportion of defective particles produced within an infected individual might be much higher than 18%. Previous results had shown that only one of 10 proviral clones derived from an infected brain sample had the ability to replicate in tissue culture. The significance of defective viruses in HIV pathogenesis should be further investigated.
Evolutionary relationships among clones
The 109 cloned sequences from PBMCs were genetically related but existed as distinct quasispecies. Only three identical pairs were detected among the 109 clones. Phylogenetic tree analysis of PR sequences suggested that the majority of PBMC quasispecies belonged to 2 distinct groups. Considering the fact that most mutations were synonymous, we further examined protein sequences. In general, the clones could be assigned into 4 groups and 31 subgroups (Table 3). The 82A, 90M, 46I/L, 54V and 64V mutations within PR were important in assignment of the clones to groups, i.e. A-D. Genotyping of plasma viral RNA only detected 82A (Table 1) while proviral DNA clones from PBMCs revealed that about half of them contained either 82A (54.2%) or 90 M (45.8%)(Table 2). These findings indicate the existence of two distinct protease resistance pathways. Table 3 shows that clusters existed among PR46I/90M for group A, PR46L/54V/64V/82A for group B, PR46L/90M for group C, and PR46I/82A for group D, and that 89.9% of the samples belonged to group A (40.4%) or B (49.5%), while groups C and D represented 5.5% and 4.6% of samples, respectively.
Phylogenetic tree analysis of RT suggested that the PBMC quasispecies were closely related, different from HxB2 (not shown). An analysis of drug resistance mutations showed that M184V might exist alone while others mutations usually co-existed with other mutations in the quasispecies (Table 3). The major quasispecies within the PBMCs were subgroups containing PR4IK/46I/90M/RT184V and PR41K/46L/54V/82A/RT184V, and these were represented at 18.3% and 10.1%, respectively (Table 3).
No apparent relationship existed between resistance-associated mutations in the PR vs. RT genes (Table 3), reflecting the fact that they emerged independently during therapy. In contrast, relationships did exist in regard to mutations within either PR or RT. For example, all the PR sequences contained either 82A or 90 M but not both, while 46L/54V/64V clustered mostly with 82A and 46I clustered mostly with 90 M, i.e., groups A and B. Either 41K or 35D, but not both, occurred alongside 90 M or 82A. In the case of RT, 184V or 181C/190A/219R might have emerged independently since both were detected. When combined with the 4 patterns of mutations in PR, i.e. 35D/46I/90M, 41K/46I/90M, 35D/46L/54V/64V/82A and 41K/46L/54V/64V/82A, it appears as though 8 distinct lineages might have existed in this subject from which quasispecies might have been derived. For example, RT181C/184V/190A/219R might have emerged from 181C/190A/219R while RT101E/108I/181C/184V/190A/219R might have emerged from RT181C/184V/190A/219R. In this regard, there was a 3 month interruption in the use of efavirenz (March to June, 1998) that might have contributed to the selective development of the two distinct NNRTI resistance pathways.
The low frequencies of groups C and D suggest that they might have been generated as a result of recombination within the PR region with the cross-over site potentially being between positions 46 and 82. To reduce the possibility of recombination during PCR, we diluted our samples during the second PCR and limited the number of cycles. A control reaction, that employed a mixture of known DNAs, containing wt or mutated PR-RT mutations, revealed that the recombination rate was lower than 3 percent in our experiments (not shown). Several quasispecies re RT might also have been generated through recombination, e.g. RT69N/101E/108I/181C/184V/190A/219R. Variant PR41K/46L/54V/64V/90M could have emerged from PR4IK/46L/54V/64V/82A and 41K/46I/90M crossed over between positions 64 and 82. The RT variant, RT69N/101E/108I/181C/184V/190A/219R could have been generated by the crossing over of RT181C/184V/190A/219R and RT69N/101E/108I/184V between positions 108 and 181.
Infection by HIV is a highly dynamic process, with viral turnover rates being as high as 1010 particles per day. Although PBMCs are probably not the most important site of viral replication in the body, they are important in viral spread and represent a source of virus particles that are found in plasma. Hence, PBMCs carry information about the evolutionary process of HIV infection while not necessarily reflecting the majority of viral species at any point in time. Our results are consistent in that we were able to detect certain resistance-associated mutations in plasma, e.g. RT69N, even though this mutation was represented at a level of only 10.2% in cloned sequences and was not detectable in PBMCs by genotyping. The emergence of RT69N was related to the addition of d4T to the treatment regimen. Interestingly, quasispecies carrying PR46I and PR46L were detected at almost equal frequencies in PBMCs (Table 2), and quasispecies containing PR90M and PR46I constituted about half of the viral population in PBMCs but were undetectable in plasma by genotyping. This may be attributed to archival quasispecies identified in PBMCs, as compared to the more recent plasma circulating species.
Dominant quasispecies in plasma represented only a small proportion of those detected in PBMCs, and might reflect the type of drug employed, e.g. RT69N emerged following use of d4T and 3TC (subgroup 16–21 of Table 3). Dominance of PR46L/54V/64V/82A in plasma and disappearance of PR46I/90M from plasma might be related to a switch in treatment from unboosted IDV and NFV to APV, since PR46I and PR90M are primary mutations for IDV and NFV, respectively, but not for APV. The crystal structure of a multiply mutated HIV-1 PR including 46I, 82A, 90M reveals an expanded active site cavity compared to that of wild type. Quasispecies carrying PR46I/90M combined with other mutations, i.e. 10I/16E/37T/48V/63P/77I, are probably not as replication-competent as PR46L/48V/54V/64V/82A under APV pressure.
In our study, wt sequences in regard to PR or RT were not detected in any of 109 clones. However besides 3 identical pairs, 98% (106/109) of these sequences were related but different. All could be divided into about 31 subgroups, with resistance mutations existing in different combinations. The existence of large quasispecies can serve as reservoirs for the emergence of viral species resistant to different drugs, since infected cells might survive for several years [20, 21]. Our results on uncultured, infected PBMCs are representative of only a tiny fraction of in vivo infection but indicate an astonishingly high degree of heterogeneity in regard to resistance.
Yamaguchi Y, Gojobori T: Evolutionary mechanisms and population dynamics of the third variable envelope region of HIV within single hosts. Proc Natl Acad Sci U S A. 1997, 94 (4): 1264-1269. 10.1073/pnas.94.4.1264.
Zhang L, Diaz RS, Ho DD, Mosley JW, Busch MP, Mayer A: Host-specific driving force in human immunodeficiency virus type 1 evolution in vivo. J Virol. 1997, 71 (3): 2555-2561.
Domingo E, Menendez-Arias L, Holland JJ: RNA virus fitness. Rev Med Virol. 1997, 7 (2): 87-96. 10.1002/(SICI)1099-1654(199707)7:2<87::AID-RMV188>3.0.CO;2-0.
Goodenow M, Huet T, Saurin W, Kwok S, Sninsky J, Wain-Hobson S: HIV-1 isolates are rapidly evolving quasispecies: evidence for viral mixtures and preferred nucleotide substitutions. J Acquir Immune Defic Syndr. 1989, 2 (4): 344-352.
Wolinsky SM, Korber BT, Neumann AU, Daniels M, Kunstman KJ, Whetsell AJ, Furtado MR, Cao Y, Ho DD, Safrit JT: Adaptive evolution of human immunodeficiency virus-type 1 during the natural course of infection. Science. 1996, 272 (5261): 537-542. 10.1126/science.272.5261.537.
Vignuzzi M, Stone JK, Arnold JJ, Cameron CE, Andino R: Quasispecies diversity determines pathogenesis through cooperative interactions in a viral population. Nature. 2006, 439 (7074): 344-348. 10.1038/nature04388.
Lemey P, Rambaut A, Pybus OG: HIV evolutionary dynamics within and among hosts. AIDS Rev. 2006, 8 (3): 125-140.
Tebit DM, Nankya I, Arts EJ, Gao Y: HIV diversity, recombination and disease progression: how does fitness "fit" into the puzzle?. AIDS Rev. 2007, 9 (2): 75-87.
Liner KJ, Hall CD, Robertson KR: Impact of human immunodeficiency virus (HIV) subtypes on HIV-associated neurological disease. J Neurovirol. 2007, 13 (4): 291-304. 10.1080/13550280701422383.
Stingele K, Haas J, Zimmermann T, Stingele R, Hubsch-Muller C, Freitag M, Storch-Hagenlocher B, Hartmann M, Wildemann B: Independent HIV replication in paired CSF and blood viral isolates during antiretroviral therapy. Neurology. 2001, 56 (3): 355-361.
Palmer S, Kearney M, Maldarelli F, Halvas EK, Bixby CJ, Bazmi H, Rock D, Falloon J, Davey RT, Dewar RL, Metcalf JA, Hammer S, Mellors JW, Coffin JM: Multiple, linked human immunodeficiency virus type 1 drug resistance mutations in treatment-experienced patients are missed by standard genotype analysis. J Clin Microbiol. 2005, 43 (1): 406-413. 10.1128/JCM.43.1.406-413.2005.
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.
Ntemgwa M, Wainberg MA, Oliveira M, Moisi D, Lalonde R, Micheli V, Brenner BG: Variations in reverse transcriptase and RNase H domain mutations in human immunodeficiency virus type 1 clinical isolates are associated with divergent phenotypic resistance to zidovudine. Antimicrob Agents Chemother. 2007, 51 (11): 3861-3869. 10.1128/AAC.00646-07.
Sanchez G, Xu X, Chermann JC, Hirsch I: Accumulation of defective viral genomes in peripheral blood mononuclear cells of human immunodeficiency virus type 1-infected individuals. J Virol. 1997, 71 (3): 2233-2240.
Tominaga K, Kato S, Negishi M, Takano T: A high frequency of defective vif genes in peripheral blood mononuclear cells from HIV type 1-infected individuals. AIDS Res Hum Retroviruses. 1996, 12 (16): 1543-1549.
Li Y, Kappes JC, Conway JA, Price RW, Shaw GM, Hahn BH: Molecular characterization of human immunodeficiency virus type 1 cloned directly from uncultured human brain tissue: identification of replication-competent and -defective viral genomes. J Virol. 1991, 65 (8): 3973-3985.
Perelson AS, Neumann AU, Markowitz M, Leonard JM, Ho DD: HIV-1 dynamics in vivo: virion clearance rate, infected cell life-span, and viral generation time. Science. 1996, 271 (5255): 1582-1586. 10.1126/science.271.5255.1582.
Michael NL, Chang G, Ehrenberg PK, Vahey MT, Redfield RR: HIV-1 proviral genotypes from the peripheral blood mononuclear cells of an infected patient are differentially represented in expressed sequences. J Acquir Immune Defic Syndr. 1993, 6 (10): 1073-1085.
Logsdon BC VJF: Crystal structures of a multidrug-resistant human immunodeficiency virus type 1 protease reveal an expanded active-site cavity. J Virol. 2004, 78 (6): 3123-3132. 10.1128/JVI.78.6.3123-3132.2004.
Ramratnam B MJE: The decay of the latent reservoir of replication-competent HIV-1 is inversely correlated with the extent of residual viral replication during prolonged anti-retroviral therapy. Nat Med. 2000, 6: 82-85. 10.1038/71577.
Palmer S, Maldarelli F, Wiegand A, Bernstein B, Hanna GJ, Brun SC, Kempf DJ, Mellors JW, Coffin JM, King MS: Low-level viremia persists for at least 7 years in patients on suppressive antiretroviral therapy. Proc Natl Acad Sci U S A. 2008, 105 (10): 3879-3884. 10.1073/pnas.0800050105.
This research was supported by the Canadian Institutes for Health Research. We thank Maureen Oliveira and Daniela Moisi for technical assistance and Monica Romero for preparation of the manuscript.
YQ did the experiments, analyzed data and wrote the first version of the manuscript. BGB and MAW supervised the research and revised the manuscript. AD provided all the clinical samples for this work, contributed intellectually to the design of the experiments, and assisted with manuscript revision. All authors read and approved the final manuscript.