Previous Article | Next Article ![]()
Journal of Virology, May 2003, p. 5721-5730, Vol. 77, No. 10
0022-538X/03/$08.00+0 DOI: 10.1128/JVI.77.10.5721-5730.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Wilscott E. Naugler,1,
Nicole H. Tobin,1 and James I. Mullins2,3,6
Departments of Pediatrics,1 Laboratory Medicine,2 Microbiology,3 Biostatistics,4 Medicine, University of Washington, Seattle, Washington,6 Department of Pediatrics, Oregon Health and Science University, Portland, Oregon5
Received 2 December 2002/ Accepted 14 February 2003
|
|
|---|
0.9 log10, resulting in lower absolute DNA levels (P = 0.007). This study demonstrates the variable rates of viral replication when HAART has suppressed plasma HIV-1 RNA for years to a median of <50 copies/ml and that combinations of four or five antiretroviral drugs suppress viral replication even after short-term virologic failure of three-drug HAART and despite ongoing accumulation of drug-resistant mutants. Furthermore, the decrease of cellular HIV-1 DNA to low absolute levels in those without genetic evidence of viral replication suggests that monitoring viral DNA during HAART may gauge low-level replication. |
|
|---|
To gain insight into the emergence of drug resistant HIV-1 during apparently effective therapy, the viral population in the peripheral blood mononuclear cells (PBMC) of children was characterized when plasma HIV-1 RNA levels were below the limit of detection over a median of 5.1 years of HAART. Of particular interest were whether the evolution of resistant virus could be detected within PBMC prior to a rebound of plasma HIV-1 RNA to detectable levels, whether mutations accumulated over time with the eventual loss of suppression of viral replication, and whether simple identifiers were associated with low-level viral replication during apparently effective HAART.
|
|
|---|
Sampling of HIV-1 populations by limiting dilution PCR of PBMC DNA and quantification of DNA copy number in CD4+ cells. DNA for PCR was obtained from a known quantity of PBMC by lysis with a detergent-proteinase K solution (62) or by extraction using the IsoQuick nucleic acid extraction kit (Orca Research Inc, Bothell, Wash.). HIV-1 DNA was quantified by endpoint-dilution PCR (44). The viral DNA load per 106 CD4+ cells was then calculated. Specimens were diluted so that no more than 30% of the reactions were positive after nested PCR, providing 70 to 80% probability that a positive reaction was from a single viral template.
Nucleic acids encoding protease (PR), reverse transcriptase (RT), and the C2-V5 region of Env were coamplified in first-round PCRs using either two or three outer primer sets in the same reaction. The PR- and RT-encoding regions were amplified either individually, with PRL (5'-GGGACCAGCGGCTACACTAGAAGAAATGATGACAGCATGTCAGG-3') and PR2 (5'-GGAGTATTGTATGGATTTTCAGGCC-3') plus RT1 and RT2 (18), or amplified as one fragment with PRL and RT2. In either case, env was simultaneously amplified using ED31 and BH2 (17). Second-round PCRs were conducted separately, using pro primers PRC (5'-CTCCCCCTCAGAAGCAGGAGCCGATAGACAAGGAACTGTATCC-3') and PR4 (5'-GGGCCATCCATTCCTGGC-3'), rt primers RT4 and RT3 (18) (or pol primers PRC and RT3), and env primers ES7 and ES8 (9). Samples that did not amplify with the forward PR primers PRL and PRC were amplified using either pol-forward outer and pol-forward inner (19), PRF1 (5'-GAGCCAAGTAACAAATTCAGC-3') and PRF2 (5'-CACCAGAAGAGAGCTTCAGGT-3'), or PRA (5'-CCTAGGAAAAAGGGCTGTTGGAAATGTGG-3') and PRB (5'-ACTGAGAGACAGGCTAATTTTTTAGGGA-3'). PCR and direct DNA sequencing of PCR products were performed as previously described (45). The resulting sequences with evidence of polymorphic residues, thus indicating the presence of multiple templates, were excluded from our analyses.
Sequence analysis.
Sequences were assembled and error checked using Sequencher 3.0 (Gene Codes, Ann Arbor, Mich.). Those with substantial G
A mutational bias (as determined using HYPERMUT [28]), suggesting hypermutation (24, 57), were omitted from further analyses. Sequence alignments were obtained using ClustalW 1.7 (55) and manually edited as necessary. Regions of ambiguous alignment were removed from subsequent evolutionary analyses. Neighbor-joining phylogenetic trees of each subject's Env, PR, and RT coding sequences were constructed using PAUP* version 4.0b4 (54) with evolutionary models selected using the Akaike information criterion (1) under Modeltest 3.0 (42). Consequently, most env trees used the general time-reversible model with gamma-distributed heterogeneity of substitution rates (31). Model parameters are available from the authors upon request. Phylogenetic trees of the nucleotide sequences encoding HIV-1 Env, PR, and RT from each patient were constructed separately and rooted using HIV-1 sequences from his or her mother, when available, or several closely related HIV-1 sequences from GenBank. The reliability of clustering in phylograms was assessed by bootstrapping analyses (15). More complete information related to the phylogenetic analyses of virus from all subjects can be found at the website http://ubik.microbiol.washington.edu/HIV/Frenkel-2/index.html.
The most recent common ancestor (MRCA) sequence of a given infection was inferred as the sequence at the node that included all of the sequences for a given gene region (PR, RT, or env C2-V5) for a given patient. This sequence was obtained using maximum likelihood estimation using PAUP*, and the divergence rates from the MRCA were calculated using linear regression over time (30).
Statistical analyses. Logistic regression analysis was applied to the prevalence of genetic mutants over time, with tests to determine whether the coefficient for time was significantly different from zero. When analysis was conducted on more than one individual, generalized estimating equations were used to account for the correlations resulting from repeated measures from the same individual (10). The log10-transformed HIV-1 DNA levels before and during potent HAART were compared using Student's two-tailed t test.
Nucleotide sequence accession numbers. The GenBank accession numbers of the HIV-1 pol and env sequences derived in this study are AY075701 to AY077450.
|
|
|---|
![]() View larger version (49K): [in a new window] |
FIG. 1. HIV-1 RNA and CD4+ cell levels and antiretroviral treatments for eight subjects (data for subjects G1 and G2 are provided below in Fig. 4 and 5, respectively). The plasma HIV-1 RNA and CD4+ cell levels and antiretroviral drug regimens are shown for nucleoside analog-experienced children whose viral genetics were studied during a median of 5.1 years of HAART. Plasma HIV-1 RNA rebounded in G1, G2, H1, and R1 shortly after starting HAART. All subjects except H3 were ultimately treated with potent HAART, defined as combination therapy with three classes of drugs (NRTI, NNRTI, and PI) or with two NRTI and two PI, during which the median plasma viral RNA level was <50 copies/ml.
|
95% (missing zero to three doses per month) in 7 of the 10 children: in G2, H3, L1, R1, and B1 from the time therapy was initiated; in G1 after resolution of pneumonia that developed during the month he began three-drug HAART; and in H1 after the placement of a gastrostomy tube prior to starting a second potent HAART regimen. H2 had
95% adherence except during a brief time during his fourth year of therapy. In contrast, subjects M2 and S1 reported missing one to two doses per week. In addition, M2 missed multiple doses during 1 month in her fifth year of HAART due to problems within her family. Plasma HIV-1 RNA levels were determined a median of six times per year (range, two to nine) after falling below the limit of detection (Table 1; see also Fig. 1, 2, 4, and 5). Four children had a rebound of viral plasma RNA to >500 copies/ml during their initial two-class, three-drug HAART regimens (two NRTI plus one PI). Two of these, G2 (see Fig. 5) and R1 (see Fig. 1), reported complete adherence to HAART during this period; however, their regimens included nucleoside analogs that they had taken previously. Viral rebound in the other two was associated with difficulties in adherence, due to illness (G1 [see Fig. 4A]) and the disagreeable taste of ritonavir (H1 [see Fig. 1]). HAART potency was increased by the addition and switching of antiretroviral agents in these four and in a fifth child (L1 [see Fig. 1]) who did not have plasma viral rebound but began HAART with only two drugs. After these modifications in therapy, 9 of 10 children were receiving potent HAART. During potent HAART, plasma HIV-1 RNA levels were sustained at a median <50 copies/ml for a combined total of 40 patient years. Rebounds of >250 copies/ml occurred in only two during potent HAART, in H2 and M2, with 589 and 1,215 copies/ml, respectively, both during periods of nonadherence.
|
View this table: [in a new window] |
TABLE 1. Summary of subjects' characteristics and HIV-1 evolution
|
![]() View larger version (26K): [in a new window] |
FIG. 2. Summary of virologic studies. Subjects, identified by a letter and number, are grouped by intensity of therapy. "HAART" was treatment with NRTI and one PI, except for L1 who received a two-drug regimen of stavudine and ritonavir. "Potent HAART" was treatment with three classes of drugs, except for G1 who was prescribed two NRTI and two PI. The ranges of plasma HIV-1 RNA levels are indicated by vertical bars, as is the lower limit of detection of the assay used (shaded areas); a horizontal bar () indicates that all values were below the limit of detection. A decrease of 0.9 log10 in HIV-1 DNA in CD4+ cells within 1 to 2 years of potent HAART is indicated by an asterisk. Findings from sequence analyses are indicated in the boxes as follows: light gray, genetic distances of env from MRCA were stable or increased during HAART; dark gray, selection of new and/or increases in established HIV-1 drug-resistant mutations; black, evolution of HIV-1 in phylogram confirmed by distances from MRCA; white, decay of drug-resistant mutant prevalence or no phylogenetic evidence of viral evolution. na, not applicable due to no pre-HAART drug-resistant mutants, or not assessed. Viral replication and associated viral evolution were diminished in association with potent HAART.
|
![]() View larger version (43K): [in a new window] |
FIG. 4. Clinical history and viral sequence analyses of a subject with regressive viral evolution during potent HAART. (A) Clinical course and laboratory values, including drug history, plasma HIV-1 RNA levels, and CD4+ and CD8+ cell numbers. (B) The genetic distances of these sequences from the inferred MRCA for pol encoding PR are shown. (C) The neighbor-joining phylogenetic analysis for HIV-1 encoding PR is shown, with multiple iterations showing the virus detected at the indicated time points as filled circles and viruses from past time points as open circles. Pre-HAART specimens demonstrated time-ordered evolution, and within 5 months of beginning three-drug HAART in 1996 a drug-resistant mutant (N88S, within box) was selected. Subsequently, and in association with the intensification of HAART, there was a gradual diminution of viral sequences that grouped with more recent time points, except for a few drug-resistant variants (N88S) that persisted but did not evolve or increase in prevalence. After 5 years of HAART, the detection of PBMC sequences that grouped with those from earlier time points suggested die-off of more recently infected cells, also documented by the decrease in distance from the MRCA.
|
![]() View larger version (39K): [in a new window] |
FIG. 5. Clinical history and viral sequence analyses of a subject with evidence for significant viral replication and forward viral evolution during HAART with a median viral load of <50 copies/ml. (A) Clinical course and laboratory values, including drug history, plasma HIV-1 RNA levels, and CD4+ and CD8+ cell numbers. (B) Genetic distances of HIV-1 RNA sequences from the inferred MRCA for pol encoding PR. (C) Neighbor-joining phylogenetic analysis for HIV-1 encoding PR, with multiple iterations showing the virus detected at the indicated time points as filled circles and viruses from past time points as open circles. Three months after initiating HAART, when no plasma HIV-1 RNA was detected, mutations associated with resistance to ritonavir and other protease inhibitors were detected in PBMC. Viral rebound to 1,200 copies/ml occurred after 6 months of HAART. The regimen was intensified, after which the median plasma viral load was <50 copies/ml; however, during 3.8 years of potent HAART this patient had 9 plasma HIV-1 RNA determinations between 50 and 101 copies/ml, and on 10 occasions the patient had <50 copies/ml. Plasma viral variants (circled) grouped with PBMC-derived viral DNA from early in the course of infection (near the root of the 14-year-old tree, at left) and with highly mutated PBMC virus (in box) that included variants with new drug resistance mutations to all three classes of antiretrovirals. Low-level viremia ceased, and the HIV-1 DNA in CD4+ cells decreased by >1 log10 when lopinavir-ritonavir was substituted for ritonavir and saquinavir.
|
0.9 log10 between 1 and 2 years of potent HAART (Table 2). The pre-HAART viral DNA load in these four children was similar (P = 0.731) to that in those who had a <0.9 log10 decrease; however, their post-HAART loads were significantly lower (P = 0.007; Student's t test). |
View this table: [in a new window] |
TABLE 2. Marked decline and low absolute HIV-1 DNA load after 1 to 2 years of HAART in subjects without genetic evidence for viral replicationa
|
![]() View larger version (26K): [in a new window] |
FIG. 3. Frequency of drug-resistant mutations during HAART. The percentages of sequences with drug-resistant mutations in PBMC are shown during a median of 5.1 years of effective HAART. (A) Mutants selected by the HAART regimen increased slowly in relative frequency (P < 0.0001), indicating ongoing low-level viral replication, although plasma HIV-1 RNA levels were <50 copies/ml. (B) Over the same period of time, mutants not selected by the ongoing HAART decreased in frequency (P = 0.0013), either due to relative increases in selected viruses or to a die-off of cells with more recently selected virus and relative persistence of virus from early infection. (C) Shifts in the prevalence of mutants in subject H1 show selection of L74V during didanosine-containing HAART and diminution of M184V after cessation of lamivudine (linear regression of slope of L74V versus M184V frequencies after cessation of lamivudine, P = 0.018).
|
Although no sustained viral rebound was detected among subjects during potent HAART (Table 1 and Fig. 1 and 2), episodes of transient viremia were detected in subject G2 with plasma viral RNA levels between 50 and 101 copies/ml on 9 of 19 occasions during 3.8 years of three-class, five-drug HAART (Fig. 5). The prevalence of PI-resistant mutants increased in his PBMC (0% in 1996, 30% in 2000, and 55% in 2001; P = 0.0005) and in his plasma (58% in late 1999 to 100% in 2000). Notably, mutations conferring high-level resistance to all three classes of antiretrovirals predominated in his plasma and PBMC virus load (RT had T69SEA, K103N, V106A, Y181C, M184V, T215F/Y, and K219Q; PR had L10I, I54V, A71V/T, V82A/F, I84V, and L90M), yet viral rebound was not sustained. Subsequent to the substitution of lopinavir-ritonavir (Kaletra) for saquinavir and ritonavir, his plasma HIV-1 RNA was <50 copies/ml in all nine assays over 19 months and his viral DNA load decreased by >1 log10/106 CD4+ cells to 1.8 log10 copies/106 CD4+ PBMC.
Drug-resistant mutants detected in PBMC prior to the initiation of HAART generally increased in prevalence when selected for by three-drug and potent HAART regimens (P < 0.0001) (Fig. 3A), although in one subject, G1, the frequency of mutants (V75T) decreased over several years of HAART. Simultaneously, mutants selected by drugs in previous but not the current regimens decreased in prevalence (P = 0.0013) (Fig. 3B). An absence of relevant drug-resistant mutants prevented the evaluation of all subjects for these trends.
Phylogenetic and genetic distance analysis. Three viral phylogenetic patterns, forward, no, and regressive evolution, were determined by statistical analysis of genetic distances from the MRCA (Table 1). "Forward evolution" (a statistically significant increase in distance from the MRCA) occurred in G2 (Fig. 5B) and to a lesser degree in M2 (PR distance [P = 0.056]) during potent HAART (Table 1). "No evolution" was detected in H1, H2, H3, R1, and S1. Nonetheless, low levels of replication (Fig. 3A) were indicated in H1 and H3 by increases in the prevalence of certain drug-resistant mutants during potent and nonpotent HAART, respectively, and were suggested in R1 and S1 by persistence of genetic diversity during HAART (data not shown). However, there was no diversification of the viral population of H2, who began HAART during primary infection shortly after a course of peripartum zidovudine (data not shown). "Regressive evolution," with viral sequences becoming more similar to the MRCA of infection (i.e., to early virus) during HAART compared to pre-HAART sequences, was pronounced in three subjects, B1, G1, and L1, during potent HAART (Table 1). Notably, the three subjects with regressive evolution and H2 had pronounced decreases in viral DNA load (Table 2).
Another important observation in the phylogenetic analysis was that plasma viral RNA detected intermittently at low levels (50 to 101 copies/ml in G2 [Fig. 5C] and R1 [data not shown]) included genotypes typical of those found shortly after infection, indicating that viral genomes similar to those in long-lived cells were replication competent.
|
|
|---|
Viral replication during HAART when plasma viral RNA was below the level of detection was indicated in six children by several different analyses of viral genotypes. Phylogenetic analysis, as used in previous studies (60, 63), revealed viral replication. However, this was the least sensitive measure of replication, detecting replication in 1 of the 10 children studied. Ranking the methods utilized in this study according to increasingly more sensitive measurements of viral replication, after phylogenetic analysis, was as follows: increasing genetic distances from the MRCA; the selection of novel drug-resistant mutants; increases in the frequency of drug-resistant mutants antedating HAART; and stable genetic distances, or divergence, from the MRCA. The four children without evidence of replication by any of the aforementioned tests, including no viral diversification following primary infection or with genetic distances that approached the MRCA, were assumed to have insignificant viral replication.
HIV-1 DNA levels normalized to CD4+ PBMC appeared to be a sensitive indicator of residual viral replication in the children we studied. While HIV-1 DNA levels before HAART were similar among children that did and did not have evidence of replication during HAART, those without replication had a
0.9 log10 decrease in their HIV-1 DNA and significantly lower absolute HIV-1 DNA concentration in CD4+ PBMC after 1 to 2 years of HAART, compared to no change in levels among children with replication. These data suggest that HIV-1 DNA levels might be useful in gauging low-level replication among individuals with plasma HIV-1 RNA below the limits of detection. Data from others' studies similarly suggest that PBMC HIV-1 DNA loads correlate with markers of viral nucleic acid synthesis (2, 4, 20, 27, 37, 61), as have intracellular viral RNA levels (4) and the ratio of unspliced to multiply spliced HIV-1 mRNA (19). Decreases in total HIV-1 levels associated with viral replication well-suppressed by HAART appear primarily due to decreases in unintegrated HIV-1 DNA (27). HIV-1 DNA levels have also correlated closely with the level of plasma RNA, grouped as <3, 3 to 50, and 50 to 200 copies/ml (61), and a strong correlation was observed between viral DNA and intracellular viral RNA levels (P = 0.005; r = 0.69) among 18 subjects with plasma HIV-1 RNA levels of <200 copies/ml, including 15 with <20 copies/ml (4). In contrast to our study, these data correlating HIV-1 RNA and DNA levels do not prove that full cycles of viral replication with infection of new cells were ongoing. Viral particles and mRNA can be derived from provirus without resulting in the infection of additional cells. Thus, it would follow that individuals with higher levels of provirus would produce more viral RNA. Our study, by demonstrating increases in the proportion of viral DNA with drug resistance mutations, indicates full cycles of viral replication were ongoing, with HAART failing to completely suppress the infection of additional cells.
Genetic bottlenecks affecting viral replication capacity (8, 33; T. Wrin, A. Gamarnik, N. Whitehurst, J. Beauchaine, J. M. Whitcomb, N. S. Hellman, and C. J. Petropoulos, abstract from the 5th International Workshop on HIV Drug Resistance and Treatment Strategies 2001, Antivir. Ther. 6:20, 2001) and pharmacologic barriers (12, 40) imposed by HAART appeared critical to limiting viral replication among the children we studied. These effects were demonstrated by the suppression of viral replication by potent four- or five-drug HAART after virologic failure of three-drug regimens. The selection of increasingly drug-resistant but implicitly poorly fit viruses in one subject during nearly 4 years of potent HAART demonstrated that the genetic barrier posed by a therapeutic regimen could persist for a sustained period of time. The importance of pharmacologic barriers was observed in this same subject when episodes of transient viremias ceased following the substitution of lopinavir-ritonavir for saquinavir and ritonavir in his five-drug regimen. Furthermore, strengthening of the pharmacologic barrier by lopinavir-ritonavir was associated with a decrease in HIV-1 DNA in his CD4+ cells to an absolute level in the range of those in the children without viral genetic evidence of replication. While we suspect that the replication capacity of his drug-resistant mutants was impaired, this was not assessed in vitro (52; Wrin et al., Antivir. Ther. 6:20, 2001), and thus we did not estimate the relative contribution of the genetic and pharmacokinetic barriers of his HAART on the suppression of viral replication.
In the one case evaluated, new drug-resistant mutants were detected in PBMC DNA prior to rebound of plasma viral RNA. This finding suggests that monitoring virus in PBMC by sensitive assays (3, 13, 47, 53, 56, 58) may have clinical utility in forecasting viral rebound.
Among our subjects with 5 to 6 years of profound suppression of viral replication, the persisting PBMC viral genotypes were similar to virus detected early in the course of HIV-1 infection. The abundance of early genotypes in the persisting viral population most likely results from a relatively greater number of long-lived cells becoming infected during the period of acute infection compared to similar time intervals during later stages of disease. The persisting early viral population appeared to include replication-competent virus, since HIV-1 genotypes typical of early infection were detected in the plasma of two children we studied and by others studying latent cellular and low-level virus in the blood (23, 46). Importantly, PBMC containing drug-resistant mutants, while decreasing in relative frequency during HAART, persisted at low levels in all individuals for years following the cessation of the selecting drug, as others have observed (23, 46). The relative loss of more recently infected PBMC and the persistence of virus from early infection may explain why plasma viremia during effective HAART is often genetically similar to archival virus (23) and why the viral phenotype has been observed to revert from X4/syncytium- to R5/non-syncytium-inducing virus during HAART (14). However, immune reconstitution and return to a more healthful cytokine milieu could also limit the replication of X4/syncytium-inducing viruses.
In summary, analysis of viral population genetics over a median of 5.1 years of effective HAART detected ongoing replication in a significant subset of children, with the most sensitive indicators being a selection of new or an increasing frequency of preexisting drug-resistant mutants in PBMC and maintenance of the mean genetic distance from the MRCA of infection. Among children with ongoing selection of drug-resistant mutants documented, viral rebound appeared to be limited by potent four- and five-drug HAART regimens, presumably by selection of viruses with impaired replication capacity. Furthermore, subjects without genetic evidence of viral replication had a greater persistence of viral genotypes typical of early infection compared to recently evolved viral sequences. In addition, their HIV-1 DNA levels decreased markedly with HAART and to lower absolute levels compared to subjects with ongoing viral replication, providing a rationale for studies evaluating whether HIV-1 DNA levels predict subsequent virologic failure among subjects with plasma HIV-1 RNA levels of <50 copies/ml.
Present address: Department of Microbiology, University of Washington, Seattle, Wash. ![]()
Present address: Department of Medicine, Oregon Health and Science University, Portland, OR 97239. ![]()
|
|
|---|
A hypermutation of human immunodeficiency virus type I genome. J. Virol. 65:1779-1788.
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»