Previous Article | Next Article ![]()
Journal of Virology, August 2005, p. 9625-9634, Vol. 79, No. 15
0022-538X/05/$08.00+0 doi:10.1128/JVI.79.15.9625-9634.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Ann J. Melvin,1
Jennifer L. McKernan,1
Diane M. Pawluk,1,
Kathleen M. Mohan,1
Paul F. Lewis,4
James I. Mullins,2,5,6 and
Lisa M. Frenkel1,5*
Departments of Pediatrics,1 Microbiology,2 Biostatistics,3 Laboratory Medicine,5 Medicine, University of Washington, Seattle, Washington,6 Oregon Health Sciences University, Portland, Oregon4
Received 23 August 2004/ Accepted 19 April 2005
|
|
|---|
|
|
|---|
LLV appeared to be clinically insignificant in two studies, as subjects with and without LLV had similar rates of virologic failure of HAART (24, 45). Detailed studies of LLV have demonstrated viral sequences similar to that found early in the course of infection (25, 35, 38). These observations support the hypothesis that activation of latently infected cells is the source of LLV (7, 17, 25, 38). A modification of this hypothesis proposes that there is continuous low-level expression of virus from a stable proviral reservoir, with detectable virus occurring when random biological processes or assay variations increase viral levels above the limit of detection of the assays (12, 13, 35). Importantly, during effective HAART these virions would not perpetuate infection, as protease (PR) inhibitors interfere with the maturation of new virus, rendering it defective and noninfectious, and reverse transcriptase (RT) inhibitors would block infection of additional cells.
In other studies, ongoing viral replication during HAART was either suspected (12, 14, 20) or detected by phylogenetic analysis and by selection of new drug-resistant mutants (8, 21, 22). Virologic failure of HAART has also been associated with LLV, especially when repeated episodes were observed (21, 22).
To better understand their origin and significance, viral sequences derived from LLV plasma specimens were compared to sequences from peripheral blood mononuclear cells (PBMC) collected before and during HAART. The genetic distances of LLV plasma sequences from the inferred most recent common ancestor (MRCA) of infection was used to "date" the LLV virus. The presence of "ancestral" or "archived" virus (genetic distance relatively close to inferred MRCA of infection) in a subject's plasma was taken to suggest that LLV originated from production of viral particles from cells latently infected with virus prior to HAART, whereas "new" or "recently evolved" virus implied that plasma virus resulted from ongoing viral replication.
|
|
|---|
Plasma HIV-1 RNA levels. Plasma HIV-1 RNA values were measured at 1- to 4-month intervals after initiation of HAART as part of clinical care using the UltraSensitive Monitor 1.0 (Roche Diagnostics, Montclair, NJ) or, for subjects with non-B subtypes, the Versant HIV RNA 3.0 (Bayer, Emeryville, CA).
PBMC-associated HIV-1 DNA amplification. DNA was extracted from PBMC using the IsoQuick nucleic acid extraction kit (Orca Research Inc., Bothell, WA). Extracted DNA was then diluted so that 30 or fewer of 100 nested PCRs were positive, yielding a >70% probability that a single viral template was in each reaction mixture (42). Regions of pol and env were amplified in a first-round PCR with primers PRL (19), RT2 (18), ED31 (9), and BH2 (3), followed by separate second rounds of PCR to amplify pol and env. Second-round pol primers were PRC (19) and RT3 (18), and env primers were ES7 and ES8 (10). Occasionally, a three-region multiplex reaction was used, with pol regions encoding PR and RT, as well as env, amplified. In these cases, the first-round multiplex PCR used PRA (19), PR2 (19), RTA (5'-AACTTCTGTATGTCATTGACAGTCCA-3'), and RT1 (18) with the env primers listed above. The second-round PCR used PR4 (19), PRB (19) for the PR region of pol, RT4 (18) and RTB (5'-CATTTATCAGGATGGAGTTCATA-3') for the RT region of pol, and ES7 and ES8 for env. Alternative primers were then substituted when PCRs were negative using the aforementioned primers. These included ED5 and ED12 (10) for the first-round PCR of env in subject T2, PRF1 (19) and PR5 for the first-round PCR of pol encoding PR, and RTC and RTB for the first-round PCR and RTD and RT3 for the second-round PCR of pol encoding RT in subject G1. Cycling conditions for first and second rounds were 94°C for 5 min; 35 cycles of 94°C for 30 s, 55°C for 30 s, and 72°C for 1 min; followed by 72°C for 6 min.
Plasma RNA extraction and reverse transcription. RNA was extracted from 0.5 to 3 ml of plasma. A positive control for extraction was cell-free HIV-1RTMC (AIDS Research and Reference Reagent Program, Division of AIDS, National Institute of Allergy and Infectious Diseases, National Institutes of Health) (32), at a concentration of 100 viral c diluted into 1 ml of HIV-1-negative plasma.
Virions were pelleted by a 2-h centrifugation at 25,200 x g and 4°C. The pellet was suspended in 400 µl of Trizol (Life Technologies, Rockville, MD) for RNA extraction, using the manufacturer's protocol. Following a 5-min incubation, 1 µg of tRNA and 80 µl of chloroform were added. The solution was vortexed, incubated on ice for 10 min, and then centrifuged at 13,500 x g for 5 min. RNA was precipitated from the aqueous phase by addition of 250 µl of 100% isopropanol equilibrated at 20°C, collected by centrifugation at 13,500 x g at 4°C for 20 min, rinsed with 70% ethanol, and then pelleted at 16,000 x g for 5 min. The supernatant was discarded, and the RNA pellet was dried for 10 to 15 min at 56°C. The RNA was then immediately suspended in a small volume of nuclease-free water on ice with or without 0.5 U/µl of RNasin RNase inhibitor (Promega Co., Madison, WI).
The RNA from the subject's samples, the extraction control, and a control for reverse transcription (
1,000 c of HIV-1RTMC RNA) were transcribed into cDNA in a reaction with specific primers for three regions of the viral genome: pol region encoding PR with PR5 (5'-CCTGGCTTTAATTTTACTGGTACAGTTTCA-3'), pol region encoding RT with RTA, and env with BH2. The cDNA was then diluted to the end point and PCR amplified with nested primers. The first round was multiplexed, including PR5 and PRA, RTA and RTC (5'-GGCCCAAAAGTTAAACAATGG-3'), and ED31 and BH2. Three separate second-round PCRs were done, each using 2 µl of the first-round product as the template and primers ES7 and ES8 (env), PRB and UHGR (5'-CTAATGGGAAAATTTAAAGTGCA) for pol encoding PR, or RTB and RTD (5'-GGCCATTGACAGAAGAAAAAATAA-3') for pol encoding RT. PCR cycling conditions were as for the DNA amplification above, except that 45 cycles were conducted for each round.
Dideoxynucleotide sequencing. PCR products were purified using ExoSAP-IT (Amersham Biosciences, Piscataway, NJ) prior to direct bidirectional sequencing using the Big Dye Terminator Kit (Applied Biosystems, Foster City, CA) with the second-round PCR primers listed above. Amplicon lengths were 493 bp for the region of pol encoding PR (and including a portion of gag), 627 bp for the region of pol encoding RT, and 627 bp for the C2-V5 region of env. Sequencing reaction products were purified with DyeEx columns (QIAGEN, Bothell, WA) prior to analysis on a 377 or 310 automated Applied Biosystems sequencer.
Sequence analysis.
Sequences were assembled and error checked using Sequencher (Gene Codes, Ann Arbor, MI). Those with substantial G
A mutational bias as determined using HYPERMUT (31, 43), suggesting hypermutation (51), were omitted from further analyses. Sequence alignments were obtained using ClustalW 1.7 (49) and edited as necessary. Regions of ambiguous alignment were removed from subsequent evolutionary analyses. Neighbor-joining phylogenetic trees were constructed using PAUP* versions 4.0b4 through 4.0b10 (47) with evolutionary models selected using the Akaike information criterion (1) under Modeltest 3.06 (40). Model parameters were as follows: (i) for PR (GTR+I+
), equilibrium nucleotide frequencies fA = 0.3880, fC = 0.1653, fG = 0.2310, fT = 0.2175; shape parameter (
) of the
distribution reflecting site-to-site rate variability of variable sites,
= 0.5382; assumed proportion of invariable sites = 0.1638; R matrix values, RA
C = 4.2758, RA
G = 16.3460, RA
T = 1.1576, RC
G = 0.2308, RC
T = 20.2003, RG
T = 1; (ii) for RT (GTR+I+
), fA = 0.4642, fC = 0.1383, fG = 0.1755, fT = 0.2220;
= 0.7765; assumed proportion of invariable sites = 0.2908; R matrix values, RA
C = 1.0896, RA
G = 5.2510, RA
T = 0.2570, RC
G = 0.4460, RC
T = 7.7782, RG
T = 1; and (iii) for env (TVM+I+
), fA = 0.4734, fC = 0.1704, fG = 0.1603, fT = 0.1959;
= 0.8087; assumed proportion of invariable sites = 0.1198; R matrix values, RA
C = 1.5222, RA
G = RC
T = 3.9416, RA
T = 0.5910, RC
G = 0.6665, RG
T = 1. Supplemental data are available at our website (http://ubik.microbiol.washington.edu/HIV/Tobin-1/). Trees were constructed using all available sequences from the children and their mothers, along with outgroup HIV-1 sequences from GenBank. A sequence corresponding to the MRCA of infection was estimated at the basal node for each gene region of each subject using maximum-likelihood estimation and the evolutionary model determined for the respective gene region.
A plasma specimen was defined as having identical sequences when five or more identical pol-RT sequences or two or more identical env sequences were obtained by end point dilution of the plasma specimen. The region of pol encoding PR was not evaluated for identical sequences due to a high level of sequence conservation in this region.
Statistical analysis of viral genetic divergence and LLV frequency. PBMC viral sequences from each subject were evaluated for divergence from their MRCA over time, separately for each gene region, using a linear random-effect regression model with divergence as the outcome and time since the start of HAART as the covariate of interest. The regression model indicated ongoing replication when the estimated coefficient-of-time covariate was positive and significantly different than zero for any of three gene regions. Plasma sequences were not included in the random-effect model due to an imbalance in the quantity of data from the plasma and PBMC.
The frequency of LLV for each patient was analyzed based on an overdispersed negative binomial for the number of LLV that was offset for the number of years of follow-up. An indicator variable for replicator versus nonreplicator was used to test for the difference in the rate of LLV per year between these two groups.
Nucleotide sequence accession numbers. The GenBank accession numbers for the HIV-1 env and pol sequences derived or analyzed in this study are AY483287 to AY484389 and AY075701 to AY077450 (19).
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Characteristics of subjects with intermittent LLV
|
![]() View larger version (48K): [in a new window] |
FIG. 1. HIV-1 RNA and CD4 cell levels over time. Plasma HIV-1 RNA and CD4 cell levels of subjects are shown over time. The lower limits of detection of the plasma RNA assays used are shaded. Times at which HAART was initiated, intensified, or stopped are indicated with the symbols defined in the inset. The following very low LLV values of two subjects are not discernible: B2, 51 c/ml at 12.6 years of age; V1, 53 c/ml at 8.5 years of age. LPV/r, lopinavir-ritonavir combination.
|
Viral sequences from LLV plasma diverged from the inferred MRCA of infection during HAART in 3/11 subjects (Fig. 2, subjects F1, G2, and T2). This divergence continued over time, suggesting ongoing viral replication in two subjects. Viral divergence could not be confirmed in one case (subject T2) because only single plasma PR and RT sequences were amplified from his LLV specimen. Neither LLV plasma nor PBMC viral sequences from the remaining eight subjects diverged from the MRCA during HAART (e.g., subjects H2, H4, M2, and P1 in Fig. 2). The absence of diverging viral sequences in these cases suggests that LLV plasma virus was derived from provirus archived prior to HAART.
![]() View larger version (20K): [in a new window] |
FIG. 2. Genetic divergence from the most recent common ancestor of infection during HAART. Genetic divergence from the MRCA of pol encoding RT (left side) and env encoding C2-V5 (right side) is shown for representative subjects. The upper three panels (subjects F1, G2, and T2) show viral sequences diverging from the MRCA during HAART in both plasma (gray diamonds) and PBMC (black dots), suggesting ongoing viral replication despite a median plasma HIV-1 RNA level of <50 c/ml. In the lower panels (subjects H2, H4, M2, and P1), plasma and PBMC virus did not diverge from the MRCA, suggesting that viral replication was inhibited by HAART.
|
The divergence of PBMC-associated viral DNA from the MRCA, also indicative of viral replication, was assessed in a linear random-effect regression model. The 10 subjects included in the model had specimens analyzed from immediately pre-HAART and from one or more time points during HAART. Viral divergence during HAART was detected in the PBMC viral DNA of only subject G2 (Table 2) and in all three HIV-1 gene regions evaluated. While the rate of PBMC viral divergence post-HAART in F1 was similar to that in G2 (Fig. 3), the trend in F1 was not statistically significant in the logistic model, presumably due to a shorter period of observation and fewer available sequences.
|
View this table: [in a new window] |
TABLE 2. Linear random-effects model to evaluate divergence of PBMC-associated HIV-1 sequence from the MRCA during HAARTa
|
![]() View larger version (16K): [in a new window] |
FIG. 3. Divergence of PBMC-associated HIV-1 env sequences after initiation of HAART. Viral sequence divergence from the pre-HAART specimen was plotted (vertical axis) for each specimen analyzed during HAART (horizontal axis). The point of reference for calculating divergence was the mean distance of each subject's viral sequences from his or her specimen collected immediately prior to the initiation of HAART. While the divergence slopes for F1 and G2 are similarly positive, suggesting ongoing viral replication, only the slope of G2's data was statistically significant (P = 0.00002) in the linear random-effect model, likely due to the greater amount of G2 data and a longer period of observation.
|
Drug resistance mutations in plasma and PBMC viral sequences. Selection of new drug resistance mutations during HAART, a sensitive indicator of ongoing viral replication (19), was detected in the plasma of two subjects, F1 and G2 (Fig. 4). Mutations encoding resistance to nucleoside and nonnucleoside RT inhibitors also appeared to increase in prevalence in G2s PBMC over time (G2s RT sequences are shown in Fig. 5). Drug resistance mutations detected in LLV plasma viral sequences of other subjects (R1 and S1) (Table 1) were all initially selected prior to HAART. (All sequence data are available at http://ubik.microbiol.washington.edu/HIV/Tobin-1/.)
![]() View larger version (31K): [in a new window] |
FIG. 4. Neighbor-joining phylogram of pol encoding PR with drug-resistant genotypes selected during HAART shown for subject G2. Plasma (gray diamonds) and PBMC (black dots) viral sequence evolution of G2 is shown. Selection of drug-resistant mutants (with codon and encoded amino acid listed) contributed significantly to the evolution from the MRCA (circled). Sequences are rooted to subtype B sequences from GenBank and the percentage of substitutions per site indicated by the horizontal line.
|
![]() View larger version (48K): [in a new window] |
FIG. 5. Amino acid alignment of HIV-1 pol sequences encoding RT of subject G2. An alignment of a subset of G2 viral sequences is shown to demonstrate the increasing prevalence of mutants over time in PBMC. Sequences are shown in an abbreviated format to feature amino acids that confer resistance to antiretroviral drugs (nucleoside and nonnucleoside RT inhibitors). Dots indicate amino acid identity to the subject's MRCA (G2 Ancestor Seq). Sequences from PBMC are labeled with specimen dates and cell and sequence indicators, and plasma sequences are similarly labeled with dates and plasma and sequence indicators.
|
|
|
|---|
LLV plasma gene sequences from most subjects in our study were similar to those derived from their PBMC collected prior to the initiation of HAART. Neither viral evolution nor selection of drug-resistant mutants was detected among these subjects. Detection of low levels of virus with genetic sequences typical of that archived ancestral virus during effective HAART have been reported previously (25, 29, 35, 37, 39). However, our unique detection of identical viral sequences during LLV provides strong support for the hypothesis that LLV can result from production of virus following immune activation and clonal expansion of latently infected cells (36). The detection of multiple identical viral sequences within single plasma specimens (e.g., Fig. 2, 11/13 F1 RT sequences at 602 days of HAART and 9/12 H4 env sequences at 448 days) suggests that these virions were produced from clonal expansion of activated HIV-1-infected cells. Since derivation of all viral templates was by direct sequencing of the PCR product from plasma diluted to the end point, these identical sequences were not a methodological artifact (33). Importantly, full cycles of infection should be blocked by effective HAART, as PR inhibitors result in production of defective noninfectious virus and RT inhibitors block infection of additional cells. Thus, as observed in several studies (12, 24, 45), LLV during HAART that effectively suppresses viral replication should not increase the risk of virologic failure.
The genetic analyses of LLV sequences from two of our subjects (G2 and F1) demonstrated increases in viral divergence from their MRCA, indicating that viral replication was ongoing during HAART even though their median plasma HIV-1 RNA level was maintained at <50 c/ml. The detection of drug-resistant mutants in all G2s LLV sequences (28/28 RT sequences) but only one of F1s sequences (1/33 RT sequences) suggests that selective pressure from antiretrovirals (4) was greater in G2. Prolonged "nonsuppressive" antiretroviral therapy in G2 prior to the initiation of HAART selected for drug-resistant mutants that conferred partial resistance to his later HAART regimen. Partial resistance to drugs in his HAART probably allowed low-level replication to occur despite unfailing adherence to HAART. In contrast, F1 had periods of nonadherence to HAART. While a sustained viral rebound was not detected, her suspension of HAART presumably allowed viral replication and evolution without selection of drug-resistant mutants. We considered the possibility that F1's median plasma HIV-1 RNA may not have actually been <50 c/ml, due to inaccurate quantification of her subtype D virus by the Amplicor 1.0 RT-PCR assay (28, 50). However, this seemed unlikely given that her HIV-1 RNA was also undetectable in the branched-chain DNA assay sensitive to 75 c/ml plasma, an assay that has been found to adequately quantify HIV-1 subtype D (5).
Two or more serial LLV have been associated with a greater risk of virologic failure of HAART compared to a single LLV (5.8-fold versus 2.0-fold) (21). In our study, serial LLV were detected in the two subjects with evidence of ongoing low-level viral replication (F1 and G2) and in none of the subjects that appeared to express only archived virus. When HAART was intensified in G2 by switching from ritonavir-saquinavir to lopinavir-ritonavir with continuation of the same RT inhibitors, his LLV ceased (0/10 assays with detectable HIV-1 RNA), precluding determination of whether virologic failure would have eventually ensued under his initial HAART regimen.
Selection of mutant virus during HAART has not always resulted in virologic failure (8, 19). Virologic failure is most likely modulated by multiple interacting parameters, including the genetic and pharmacologic barriers that each regimen imposes on the selection of "fit" drug-resistant mutants. Of note, our subjects with multiple and sequential episodes of LLV, F1 and G2, were treated with four to six antiretroviral agents that should have provided greater genetic and pharmacologic barriers compared to commonly prescribed three-drug HAART regimens. These relatively high genetic and pharmacologic barriers may have contributed to G2 not experiencing virologic failure even though several new primary drug resistance mutations (27) were selected, as detected in his plasma and PBMC during HAART.
LLV could not be confirmed in one-quarter of our subjects' LLV plasma specimens. The LLV in these instances could have been due to false-positive HIV-1 RNA quantification results, usually near the lower assay limits of detection, or incorrect reporting of results by the laboratory. The former appears likely given that the reproducibility of detection of LLV has varied in others studies, with a range of 41 to 97% (26, 41). If a large percentage of LLV are false positives, then clinical studies would be biased to underestimate the risk of virologic failure associated with LLV.
An important caveat to our and others' studies is that following long-term effective HAART an evaluation of sequence divergence could theoretically fail to identify low-level viral replication. The viral genotypes that persist in PBMC after long-term suppression of viral replication date mostly to early infection, with more recently evolved viruses persisting to a relatively lesser degree (19, 39, 44, 46). Virions produced by these archived viral sequences, and infecting additional cells, may not demonstrate sufficient divergence from the MRCA to be recognized as evolving virus (19). Indeed, full cycles of viral replication, even when plasma viremia is undetectable by current assays, may account for the persistence of the latent PBMC reservoir (29, 39). We have shown that sequence divergence and phylogenetic analysis of PBMC-associated virus were less sensitive in gauging low-level viral replication compared to measuring shifts in the size of drug-resistant viral populations in PBMC (19). Prior studies suggest that viral replication may be more readily detected by divergence of plasma viral RNA compared to PBMC-associated viral DNA (19, 39, 44, 46). Here, we did not directly compare the sensitivity of increases in plasma HIV-1 genetic diversity to increases in the size of drug-resistant viral populations in PBMC to gauge low-level viral replication during effective HAART.
Production of identical viral sequences in plasma provided strong evidence that LLV resulted from activation and proliferation of one or more cells infected concomitantly in the past, whereas LLV with ancestral but nonidentical sequences could have resulted from activation of multiple persistently infected PBMC or similar reservoirs of archived virus (6, 15, 30). However, for the reasons mentioned above, we cannot exclude the possibility that the latter instances of LLV resulted from low-level viral replication that was below the sensitivity of our genetic analyses.
The long-term efficacy of HAART would most likely differ, depending on which of the two processes shown here predominateproduction of virus from a subset of long-lived infected cells or ongoing viral evolution resulting from new rounds of infection. Furthermore, the inconsistent association of LLV with subsequent virologic failure of HAART in studies of large populations (21, 22, 24, 45) may be due to a greater propensity for one or the other of these processes. Assays that differentiate LLV resulting from expression of ancestral virus from those derived from ongoing productive cycles of viral replication may therefore be of clinical value.
This work was supported by an NIH training grant to the University of Washington and grants to L.M.F. and J.I.M. from the U.S. Public Health Service and to the University of Washington Center for AIDS Research.
Present address: Department of Biostatistics, University of Washington School of Medicine, Seattle, WA 98195-4806. ![]()
Present address: P.O. Box 416, Kirkland, WA 98083. ![]()
|
|
|---|
A hypermutation. Bioinformatics 16:400-401.
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»