Previous Article | Next Article 
Journal of Virology, November 1999, p. 9404-9412, Vol. 73, No. 11
0022-538X/99/$04.00+0
Copyright © 1999, American Society for Microbiology. All rights reserved.
Evolution of Envelope Sequences of Human
Immunodeficiency Virus Type 1 in Cellular Reservoirs in the Setting of
Potent Antiviral Therapy
Huldrych F.
Günthard,1,
Simon D. W.
Frost,2
Andrew J.
Leigh-Brown,2
Caroline C.
Ignacio,1
Kristin
Kee,1
Alan S.
Perelson,3
Celsa A.
Spina,1,4
Diane V.
Havlir,1
Marjan
Hezareh,1
David J.
Looney,1,4
Douglas D.
Richman,1,4 and
Joseph K.
Wong1,4,*
University of California, San
Diego,1 and San Diego Veterans Affairs
Healthcare System,4 San Diego, California;
University of Edinburgh, Edinburgh,
Scotland2; and the Los Alamos
National Laboratory, Los Alamos, New Mexico3
Received 17 March 1999/Accepted 9 July 1999
 |
ABSTRACT |
In human immunodeficiency virus (HIV)-infected patients treated
with potent antiretroviral therapy, the persistence of latently infected cells may reflect the long decay half-life of this cellular reservoir or ongoing viral replication at low levels with continuous replenishment of the population or both. To address these
possibilities, sequences encompassing the C2 and V3 domains of HIV-1
env were analyzed from virus present in baseline plasma and
from viral isolates obtained after 2 years of suppressive therapy in
six patients. The presence of sequence changes consistent with
evolution was demonstrated for three subjects and correlated with less
complete suppression of viral replication, as indicated by the rapidity of the initial virus load decline or the intermittent reappearance of
even low levels of detectable viremia. Together, these results provide
evidence for ongoing replication. In the remaining three patients,
virus recovered after 2 years of therapy was either genotypically
contemporary with or ancestral to virus present in plasma 2 years
before, indicating that virus recovery had indeed resulted from
activation of latently infected cells.
 |
INTRODUCTION |
Treatment with combination
antiretroviral agents can reduce plasma and lymph node RNA levels below
the level of detection of sensitive assays for nucleic acid (2,
12, 31, 46). The persistence of cells harboring
replication-competent virus despite periods of up to 3 years of potent
suppressive therapy has raised questions about the feasibility of
eradication of infection with currently recommended treatments (3,
9, 10, 47, 49). Despite the apparent lack of development of drug
resistance in patients with undetectable virus (9, 47), the
inferred presence of unintegrated viral DNA in peripheral blood
mononuclear cells (PBMC) in one study (3) has raised
questions about the certainty that viral replication is controlled even
in this setting. Extrapolating from such observations, it can be argued
that the persistence of this reservoir of virus reflects ongoing viral replication at low levels and continuous replenishment rather than the
existence of long-lived, latently infected cells.
To address this question, we studied sequences of the human
immunodeficiency virus type 1 (HIV-1) envelope corresponding to the
relatively conserved C2 and the hypervariable V3 domains of env from HIV-1 present in baseline plasma and in virus
recovered from in vitro culture of PBMC following 2 years of therapy.
The V3 region is an important determinant of cell tropism and is a target of host immune responses. Consequently, it is among the most
variable regions in the HIV genome (15, 19, 21, 36). Sequence differences in this region have been used to compare the
relatedness of virus populations (20, 33, 42) and to study
viral evolution within individual hosts (1a, 6, 37-39, 45),
among patients (1a, 6, 25, 26, 37, 38, 45), and across
patient populations (21, 22, 27).
Samples were studied from six individuals from whom, despite prolonged
suppression of plasma viremia, virus had been previously recovered by
coculture of PBMC (47). We hypothesized that ongoing replication should be accompanied by changes in viral envelope sequence
and that phylogenetic analysis would permit the distinction of virus
produced by activation of latently infected cells from that produced as
a consequence of continued viral propagation during the course of therapy.
(Presented in part at the 2nd International Workshop on HIV Drug
Resistance and Treatment Strategies, Lake Maggiore, Italy, June 1998.)
 |
MATERIALS AND METHODS |
Patient selection.
Subjects were chosen from a
well-characterized cohort of HIV-infected patients undergoing treatment
with indinavir, lamivudine, and zidovudine (12, 47). Each
had apparent suppression of plasma and serum viral RNA to below the
50-copy/ml level with the exception of patient M, who had 50 copies of
HIV RNA/ml on the day of PBMC sampling for virus isolation (Fig.
1). Sequence comparisons were made between virus present in baseline plasma and
virus isolated while on drug therapy at year 2.

View larger version (10K):
[in this window]
[in a new window]
|
FIG. 1.
Log RNA copies/ml in serum or plasma plotted against
days under indinavir, zidovudine, and lamivudine combination therapy.
The baseline for the y axis is the limit of detection of the
Roche ultrasensitive assay (1.6 log copies/ml). The top panel
demonstrates the response of two patients (Pt.) with steep initial declines in viral load (A and B). The middle panel
shows two patients with steep initial declines followed by transient
reappearance of low levels of detectable viral RNA (K and L). The
bottom panel shows two patients with slow initial responses (C and M)
and, in one patient (M), intermittent transient reappearance of low
levels of viral RNA in the blood.
|
|
Viral isolates.
The materials and methods for viral
isolation have been previously published (47). Early culture
supernatants were chosen to minimize ex vivo culture effects on the
composition of the viral isolates (typically at day 14 of original
cultures and before peak p24 production).
RNA extraction.
RNA was extracted from cell-free viral
culture supernatants with the Qiagen HCV RNA extraction kit (Qiagen,
Chatsworth, Calif.) and from plasma by the AMPLICOR extraction
procedure (Roche Diagnostics, Branchburg, N.J.) according to the
manufacturers' protocols.
RT-PCR.
Following denaturation of the RNA template and
primer mix at 70°C for 3 min, reverse transcription (RT) was
performed at 45°C for 1 h in the following buffer: 100 mM Tris
(pH 8.3), 10 mM MgCl, 50 mM KCl, 10 mM dithiothreitol, 250 µM
dinucleoside triphosphate, 37 U of RNAgard (Pharmacia), 2.5 U of avian
myeloblastosis virus reverse transcriptase (Gibco BRL, Gaithersburg,
Md.), 12.5 pmol of V3Bout primer (5'-ATTACAGTAGAAAAATTCCCCT-3').
Nested PCRs were then performed in quadruplicate for each sample
to minimize sampling bias. RT product (cDNA; 10 µl) was carried over
to a first-step PCR employing the primers V3Bout and V3Fout
(5'-CAAAGGTATCCTTTGAGCCAAT-3') under the following reaction
and thermocycling conditions: Tris (pH 8.3), 3 mM MgCl2,
250 µM deoxynucleoside triphosphate, 20 pmol of primers, 1 mg of
bovine serum albumin/ml, 5 U of Taq polymerase (Boehringer
Mannheim) in a 50-µl reaction volume at 95°C for 2 min, 95°C for
30 s, 55°C for 30 s, 72°C for 1 min for 35 cycles; and
72°C for 1 min. One microliter of the first-step product was then
carried over to a second-step PCR with the primers V3Fin (5'-GAACAGGACCAGGATCCAATGTCAGCACAGTACAAT-3') and V3Bin
(5'-GCGTTAAAGCTTCTGGGTCCCCTCCTGAG-3') with reaction
conditions identical to those for the first-step PCR. Ten percent of
the PCR products were resolved on 1% agarose gels and visualized by
ethidium bromide staining. We have examined the rate of
misincorporation during amplification under these nested PCR conditions
by comparing the results from sequencing of unamplified pNL4-3
(1) with those from sequencing of amplified and subsequently
cloned products and have observed a misincorporation frequency of
approximately 1 to 3 per 3,500 bases.
Cloning.
After determination that each of quadruplicate
nested PCR products yielded comparable band intensities by ethidium
staining, 5 µl of each of the quadruplicate reactions were pooled
together and gel purified by agarose gel electrophoresis. DNA was
extracted from gel slices with the Qiaex II Gel Extraction kit (Qiagen) according to the manufacturer's protocol. Two microliters of the pooled, gel-purified PCR products was ligated into the plasmid vector
pCR 2.1 and cloned according to the instructions for the TOPO Ta
Cloning system (Invitrogen, Carlsbad, Calif.) under ampicillin selection.
Sequencing.
Ten to 13 individual clones were picked based on
blue-white screening and grown in 3-ml cultures. Plasmids were purified
with the Qiaprep Turbo-8 kit (Qiagen). Approximately 300 ng of each plasmid was sequenced with the sequencing primers ForTopoEcoR1 (5'-TGGATATCTGCAGAATTCG-3') and M13 Reverse
(5'-CAGGAAACAGCTATGAC-3') with a 373A automated sequencer
and Big-dye chain terminator chemistry (Applied Biosystems, Inc.,
Foster City, Calif.). The sequences were aligned and manually edited
with Sequence Navigator 1.1 (Applied Biosystems). In total, 115 evaluable sequences were obtained for the six patients.
Phylogenetic analysis.
Tree reconstructions were performed
by a neighbor-joining method implemented with DNADIST (Kimura
distances; transition transversion ratio, 2.0) and NEIGHBOR in Phylip
version 3.54 (7) with 100 bootstrapped data sets and by a
maximum-likelihood method with fastDNAML (7, 30), applying a
transition/transversion ratio of 2.0 and a uniform substitution rate.
Graphic representation of phylogenetic trees was accomplished with
Treeview version 1.5 (30a).
Determination of pairwise distances.
Pairwise distances from
the most recent common ancestor (MRCA; represented by the most distal
node which gives rise to all sequences from the same patient) to clonal
sequences on the maximum-likelihood phylogenetic tree were estimated
both by summing relevant branch lengths generated with fastDNAML
version 1.0 (random input order; transition/transversion ratio, 2.0;
single rate of substitution) (7, 30) and by direct
computation with DNADIST (PHYLIP version 3.54; maximum-likelihood
estimates; transition/transversion ratio, 2.0) (7) from the
inferred MRCA sequences obtained with DNAML (PHYLIP version 3.6)
(7, 8). Branch lengths from fastDNAML were chosen over
those provided by DNAML (PHYLIP version 3.6) because the latter
presently assigns a minimum distance of 0.00005 to every segment of
inferred maximum-likelihood trees, including branches connecting
identical sequences. However, the two programs produced trees with
comparable maximum-likelihood values. Duplicate analyses were
performed, with selection of the MRCA based on all viral isolate
sequences and the nearest plasma sequences (changed the selection of
the MRCA for patients A, C, and M). Means and standard errors for
distances were calculated with JMP version 3.1.5 (SAS, Cary, N.C.).
Mean pairwise distances between plasma virus and viral isolates were
calculated with the program SENDBS (Naoko Takeaki) with and without
correcting for within-sample diversity (29).
Estimation of residual replication.
Residual HIV replication
was approximated by calculation of the area under the curve (AUC)
generated by plotting log plasma RNA against time. The observed AUC
(AUCObs) represented the log RNA copy by time product which
results from the decay of free virus in plasma and virus produced by
cells infected prior to institution of therapy (productively infected
cells; a reservoir of "long-lived" infected cells; and release of
virus from the follicular dendritic network in lymphoid tissues and
activation of latently infected cells) as well as from any virus
produced as a consequence of ongoing viral replication during
combination therapy.
The theoretical AUC for each patient was estimated according to the
long-lived-infected-cell model of Perelson et al. (
31).
This
calculation assumes decay in the absence of residual viral
replication
under therapy and assumes that activation of latent
cells contributes
relatively little to the viral RNA measured
during first- or
second-phase decay of plasma virus. Log
10 viral
RNA copy
number (log
10V) as a function of time
(
t) can be described
as follows:
|
(1)
|
where
V0 is the viral load at baseline,

is the decay constant of productively infected cells,
µ
M is the decay
constant of long-lived cells,
c is the clearance constant for
virions,
T0 is the CD4 count at baseline,
N is
the viral burst
size, and
k is the infection constant
(
31). Perelson et al.
observed
NkT0
values between 2.84 and 2.92 for a
T0 of 205 to
500 cells/mm
3. We utilized a mean value for
NkT0 of 2.88 and the observed virus
clearance
constant,
c, of 3 day
1 in our analysis
(
32).
The serially measured viral RNA concentrations and those calculated
from equation (1) were plotted against time under treatment
with the
program PRISM 2.0 (Graphpad, San Diego, Calif.). AUCs
were calculated
by summation of trapezoids implemented in the
program PRISM 2.0. The
predicted AUC (AUC
Pred) was computed with
a

of 1.2 day
1 and a µ
M of 0.115, based on
the best-case decay constants
observed by Perelson et al.
(
31). Duplicate analyses with a

of 0.7 day
1 and a µ
M of 0.066 based on
mean decay constants (
31)
for the calculation of
AUC
Pred resulted in very similar results
(data not shown).
The difference between AUC
Obs and AUC
Pred,
then,
represented the proportion of the RNA concentration by time
product
that should be attributable to residual viral replication for
each patient while under therapy (AUC
Net).
In general, the use of AUC as an indicator of the difference in the
degree of ongoing viral replication between patients may
be best when
their baseline RNA concentrations are similar, while
the calculation of
AUC
Net from AUC
Obs and AUC
Pred
becomes most
necessary when baseline RNA concentrations vary greatly
among
the patients being compared. In the present study, the patients
had relatively uniform starting RNA
levels.
Our calculations assumed negligible contribution to total viremia from
sporadic activation of latently infected cells (a simplification
justified by observations of phase 1 and phase 2 decay by Perelson
et
al. [
31]) and may therefore tend to underestimate the
AUC
Pred,
leading to overestimates of AUC
Net.
Because of this and other
assumptions made for the derivation of these
terms, it would be
incorrect to treat the calculated values as
definitive for residual
viral replication. Rather, they are meant to
represent a relative
magnitude of residual replication (derived from
the best available
kinetic data and based on prevailing theory in the
published literature)
and are applied here to permit reasonable
interpatient
comparison.
Statistical analyses.
Correlations of AUCNet
with pairwise genetic distance and Wilcoxon rank sums for comparison of
pairwise distances were computed with the statistical analysis program
JMP 3.1.5. (SAS) and the program Prism 2.0.
Nucleotide sequence accession numbers.
All sequences
reported here have been deposited in GenBank and were given accession
no. AF185823 to AF185937.
 |
RESULTS |
Assessment of residual viral replication based on the response of
serum and plasma RNA to drug therapy permitted segregation of the
patients into three patterns (Fig. 1). Group 1 (patients A and B) had
subjects with the steepest initial declines in viral load and without
any detectable viral RNA for the duration of this study period. Group 2 patients (K and L) were subjects who had similarly rapid initial
declines in plasma virus concentrations upon initiation of therapy but
with sporadic reappearance of low but detectable levels of viremia
later. Finally, two patients (C and M) constituted a third group
exhibiting a protracted and gradual decline in plasma virus levels with
later sporadic reappearance of viremia throughout the first 2 years of
therapy. In the case of M, viral RNA was present in plasma at the same
time that samples were obtained for viral culture (47). The
amount of residual replication as reflected by AUCNet was
low to negligible for patients A, B, and L, intermediate for K and C,
and highest for patient M (Table 1).
Phylogenetic reconstruction.
In every case, sequences from
each patient clustered together and not with those from other patients
or with prototypic laboratory-adapted and "primary" HIV clones, as
shown by a reconstruction by maximum-likelihood analysis (Fig.
2) (24). Tree topology with a
neighbor-joining reconstruction was similar: sequences from each
patient formed individual clusters with highly significant bootstrap
values (>80 per 100 bootstrapped data sets examined) in the consensus
neighbor-joining tree (data not shown). In one case (patient B), year 2 isolates were indistinguishable from sequences present at baseline and were very homogeneous. For patients A and C, the viral isolates most
closely resembled one or a subset of sequences present in baseline
plasma, while for patients K, L, and M, no plasma sequences were
identified which clearly clustered with sequences from year 2 viral
isolates, although in the case of M, a "nearest plasma clone" could
be identified. For patients M and C, branch lengths from the closest
proximal common node (representing the most recent common ancestor) to
the year 2 isolates appeared longer than those to the baseline plasma
sequences.

View larger version (17K):
[in this window]
[in a new window]
|
FIG. 2.
Maximum-likelihood phylogenetic tree displaying
relationships of sequences of clones derived from plasma collected at
the start of therapy (time zero) and from viral isolates obtained
approximately 2 years into therapy. The letters adjacent to individual
sequence clusters identify the patient source of the sequences. In
every case, the sequences from an individual patient clustered with
sequences from the same patient and with no others. Sequences of
prototypic T-lymphocytotrophic and macrophage-tropic viruses are
denoted by the respective reference names (NL4-3, SF162, etc.).
Sequences derived from baseline plasma are represented by solid
circles. Sequences derived from viral isolates are indicated by open
circles. The arrows with solid arrowheads identify nodes representing
the MRCA identified for all plasma and viral isolate sequences from
each patient. The arrows with open arrowheads identify nodes
representing the MRCA for the viral isolate cluster and the plasma
viral sequence(s) nearest it. Horizontal branch lengths are scaled and
correspond approximately to the percent nucleotide differences between
sequences and nodes. Sequences whose branch lengths were not
significantly different from zero are shown with collapsed branches.
|
|
Protein sequence comparisons.
In four patients, inferred
protein sequences demonstrated little difference between baseline and
year 2 isolates in the C2 region (Fig.
3). In the V3 region,
baseline sequences differed from year 2 viral isolates at multiple
positions in the case of patient M, while few differences were noted
for the other five subjects. For patient M, an S-to-R substitution was
seen at position 306 (according to HXB2 sequence by convention of the
Los Alamos National Laboratories) in the V3 loop which has been
associated with a switch from M- to T-cell tropism (4, 11),
but cell tropism and presumably coreceptor usage are context dependent (40). In this patient there was an accompanying R-to-S
replacement at residue 322. Cultivation of the year 2 virus in an MT-2
cell assay (34) showed that it was non-syncytium inducing
(data not shown). In the case of patient K, multiple replacements were
noted in the region downstream from the V3 loop, towards V4.

View larger version (45K):
[in this window]
[in a new window]
|
FIG. 3.
Predicted protein sequences. All letter designations for
amino acids (AA) conform to the standard International Union of Pure
and Applied Chemistry code. The HXB2 reference sequence (Los Alamos
National Laboratory) is given at the top with the region corresponding
to the V3 loop marked by a heavy line and residues 306 and 322 identified with arrows. Subject and source are identified in the first
column, with the consensus sequence in plasma at baseline given in its
entirety (based on the most frequently occurring AA at each position,
or when two AA occurred with equal frequency from the same source, the
AA corresponding to the HXB2 reference at the same position). Shown
next are clonal sequences from year 0 plasma, followed by the consensus
sequence from the year 2 viral isolates and the clonal sequences from
year 2 viral isolates. These sequences are shown with the following
abbreviations in reference to the baseline plasma consensus sequence
for the same patient: period, identity with the plasma consensus; dash,
gap inserted to maintain alignment of all sequences shown; question
mark, unresolved AA position; asterisk, nonsense mutation. Replacements
are indicated by the appropriate code letter. The frequency with which
a particular clonal sequence was encountered is shown in the first
column adjacent to each sequence. For patient M, residues 306 and 322 have been boxed to highlight the substitutions noted.
|
|
Pairwise genetic distances between baseline plasma virus and year 2 viral isolates.
Pairwise genetic distances (both corrected and
uncorrected for overall diversity of clones from each source)
demonstrated a continuum of small but measurable distances between
baseline plasma-derived sequences and 2-year viral isolates. Corrected distances demonstrated considerably greater variability than
uncorrected ones due to substantial variation among patients in the
level of diversity in plasma samples (Table
2). A more conservative calculation of
net divergence was also made by using the diversity of the time zero
plasma population as an estimate of the true diversity of virus in the
population of viably infected cells at year 2. The rank orders of
divergence in these analyses were similar, with M and K showing the
greatest divergence (data not shown).
Genetic distance to MRCA.
Based on coalescence theory
(17, 18, 35, 41), each sequence can be traced back to the
most recent common ancestral sequence (the MRCA, represented by the
most distal node common to all sequences from each patient on the
phylogenetic tree [Fig. 2]), with branch lengths proportional to the
number of elapsed generations. Here we assume that the mutation rate µ is relatively constant among individuals so that sequence
divergence is determined by gt, where g is the
number of generations (replication cycles) times time
1
and t is the elapsed time. Although mutations within the
reverse transcriptase have been reported to affect the fidelity of
HIV-1 RT (43), such mutations were not observed in these
viral isolates (47).
The difference between mean genetic distances from year 2 viral
sequences and those from year zero plasma sequences to the
MRCAs for
each patient was calculated to provide a relative measure
of
generations elapsed since the common ancestor. Larger positive
differences would be expected when year 2 isolates are more divergent
from the MRCA than baseline plasma virus, and increasingly more
negative values would be expected when the viral isolates are
genealogically older than baseline plasma virus (Table
3). When
MRCAs for each patient were
identified based on all available
sequences, the distances from the
MRCA to year 2 viral isolates
exceeded the distances from the MRCA to
the baseline plasma sequences
derived by summation of branch lengths
for three patients (M [
P = 0.0006 for Wilcoxon rank
sum], C [
P = 0.001], and K [
P = 0.0208]),
were equivalent for one (A [
P = 0.667]) and were less for two
(B [
P = 0.0023]
and L [
P = 0.0002]) subjects. The results were
qualitatively similar for distances generated by direct calculation
of
distance from inferred and observed sequences, although the
differences
only approached borderline significance for K and
were not significant
for L and C (Table
3). Analysis based on
MRCAs selected for the year 2 viral isolate clusters and the single
plasma clone nearest each year 2 cluster gave comparable results
(data not shown).
Both the summation of branch length method and the direct estimation of
distances were based on results from phylogenetic
reconstruction. As
such, they reflect the accuracies and limitations
of these methods
(
7,
16). In the case of analysis of branch
lengths, the
individual branch lengths given by fastDNAML were
typically
associated with fairly wide confidence intervals. The
confidence
intervals for the segments connecting the root of the
year 2 viral
isolate clusters and either the nearest plasma sequence
(patients A, C,
and M) or the root of plasma sequence clusters
(B, K, and L) to the
MRCA are shown in Table
4.
Correlation of evolutionary distance with residual replication
represented by AUCNet.
Mean pairwise distances between
plasma virus and virus isolated at year 2 (both corrected and
uncorrected for intracompartment diversity) correlated with
AUCNet, with r2 values of 0.55 and
0.6, respectively, but P values only approached statistical
significance (0.09 and 0.067, respectively). However, differences in
mean distance by analysis of branch lengths from MRCA to all virus
isolates and to all plasma sequences from the same patient showed a
significant correlation with AUCNet:
r2 = 0.889 (P = 0.0048). Using a
direct calculation from the inferred ancestral sequences, the
r2 value was 0.68 (P = 0.043)
with AUCNet (Fig. 4). When
MRCAs were selected based on all viral isolates but only the nearest
plasma sequence(s) for each patient, similar correlations were noted (data not shown).

View larger version (13K):
[in this window]
[in a new window]
|
FIG. 4.
Regression analysis of the AUCNet and the
differences in mean distances from plasma virus and virus isolates to
the MRCA with inclusion of all plasma sequences calculated with
summation of branch lengths (a) or by direct calculation of distance
from sequences (b). The solid line represents the line of fit. The
broken lines represent confidence curves for fit.
|
|
 |
DISCUSSION |
Latently infected resting T cells harboring replication-competent
virus persist in treated patients despite the absence of detectable
free virus in the blood (3, 9, 47). Whether this reservoir
of virus is maintained as a result of a low level of ongoing viral
replication or as a result of a slow decay rate carries important
implications for the risk of development of drug resistance and
treatment failure in such patients. This question is central to the
success or failure of any eradication strategy based on currently
available antiviral drugs.
By conventional measures, all subjects studied had excellent responses
to therapy, with viral loads rapidly falling to below 200 copies/ml of
serum or plasma upon initiation of antiviral therapy. These responses
were sustained at levels below the threshold of detection of a more
sensitive assay (<50 copies/ml) for periods approaching 2 years at the
time virus was isolated (12, 13). Previous studies of these
subjects had demonstrated the lack of development of new genotypic drug
resistance over this period (47). However, a detailed
examination of the serial HIV RNA levels in blood from these patients
revealed varying patterns based on rapidity of decline of the viral
load and in some cases on the intermittent reappearance of low-level
viremia during the course of therapy. These qualitative differences
were reflected quantitatively in the calculated AUCNet
values representing residual viremia from ongoing viral replication
(and/or, to a small extent, from activation of latently infected cells
in vivo [see Materials and Methods]). It is not known why the
antiviral drug effect may vary among these patients in spite of
identical therapies, similar prestudy drug treatment histories, and
similar incidence of preexisting drug resistance (13, 46,
47). Such differences could be attributable to individual
differences in drug absorption or drug metabolism or to variable adherence.
Pairwise genetic distances from year 2 viral isolates to corresponding
baseline plasma virus and pairwise distances from year 2 viral isolates
and baseline plasma virus to the ancestral sequence inferred from a
maximum-likelihood phylogenetic reconstruction were examined to
estimate viral evolution in these subjects. Use of simple pairwise
distances from year 2 viral isolates to baseline plasma virus would be
expected to be an accurate reflection of generations elapsed between
time points if year 2 viral isolates were, in every case, descendent
from the virus present in plasma at baseline. However, if viral
isolates at year 2 arise from latently infected cells that might have
been ancestral to those present in plasma at baseline or which might
have descended independently from a common ancestor, simple pairwise
distances and generations elapsed need not be directly proportional.
So, while those three patients (C, K, and M) having the largest
AUCNet (e.g., those with slower initial decline of viremia
and those with sporadic reappearance of viremia) had the highest mean
pairwise distances and the three subjects (A, B, and L) with the lowest
AUCNet had the lowest pairwise distances, the correlation
of these parameters failed to reach statistical significance.
To quantitatively analyze genetic relationships of virus isolated at
year 2 and virus present in baseline plasma (time zero) without having
to assume that year 2 virus descended directly from baseline plasma
virus, we examined pairwise distances from each cloned sequence to the
corresponding MRCAs. The approach of summing relevant branch segments
benefits from the incorporation of the structure of the inferred
phylogenies into the distance measurements but does not readily yield
an appropriate confidence interval. Therefore, we have also
reconstructed the MRCA sequences with a maximum-likelihood approach
(implemented in PHYLIP version 3.6) (7, 8) and directly
estimated evolutionary distance between these and the baseline (plasma)
and year 2 (virus isolate) sequences. The estimates resulting from both
of these approaches demonstrated significant correlations of
differences in distance to the MRCA with AUCNet (Fig. 4).
Findings consistent with sequence evolution in cultured virus from some
patients but not from others could be explained by a sampling artifact
(isolation of virus from few latently infected cells) or by a selection
artifact from in vitro coculture conditions with some patients but not
others (5, 23, 28). We have analyzed sequences of virus
cultured at time zero from patient A (virus isolates from time zero
were not available from the other patients). This baseline isolate
clustered more closely with virus in plasma at time zero than with
viral isolates from the year 2 time point (data not shown), indicating
greater effects from the time of sampling than from in vitro culturing
in this setting. Moreover, the strong correlation between sequence
evolution (measured by pairwise distances to the MRCA) and residual
replication inferred by the magnitude of AUCNet supports in
vivo biological differences over random sampling artifacts or artifacts
from in vitro culture to explain these differences. Additional studies
of the nature of the sequence changes observed have been performed
which show that in the cases of M and K, but not the other four
subjects, the changes in fact reflect selective differences between
year 0 and year 2 (11a).
In the case of patient M and possibly patient K, the intermittent
appearance of detectable plasma viral RNA into the second year of
therapy and the finding of evolution in the envelope sequences would be
consistent with continuous or intermittent viral replication. Alternatively, but not mutually exclusively, in the case of patient M,
where nonsynonymous substitutions in V3 were so frequent, sequences of
year 2 viral isolates might reflect differences of a viral population
maintained in an isolated anatomic compartment or a sanctuary site
(20, 33, 48, 49). For patient C, reappearance of detectable
free virus did not occur in the second year, although it was seen late
into the first year, and the initial decline of plasma viremia was much
more gradual than that seen for patient A, B, or L. In this case, the
intermediate level of sequence evolution observed might be attributable
either to viral replication occurring soon after initiation of therapy
or to low levels of ongoing replication throughout the course of treatment.
Finally, in the cases of A, B, and L, treatment was associated with
both the rapid decay of plasma virus levels and the sustained suppression of plasma viremia to undetectable levels. In these subjects, corrected pairwise genetic distances between sequences from
baseline plasma and virus isolated later in therapy were smaller. In
these three cases, mean distances from sequences of viral isolates at
year 2 were equivalent to or even shorter than mean distances from
baseline plasma to the MRCA, findings consistent with recovery of virus
from cells infected at the time of or before the baseline plasma
sequences. However, the failure to demonstrate sequence evolution for
these patients in this particular viral reservoir does not rule out the
possibility that very low levels of replication or replication in other
compartments and in sanctuary sites not accessible to drug therapy persist.
These findings have several important implications for the
understanding of HIV latency and the evaluation of antiviral therapies. First, they confirm that, in subjects who have sustained suppression of
plasma viremia documented by the frequent application of the most
sensitive assays (such as patients A, B, and L), a population of
latently infected cells likely persists into the third year of therapy
and that the maintenance of this reservoir of cells appears to be
largely independent of discernible viral replication. From the
available data, one would surmise that in such patients, the likelihood
for emergence of drug resistance should be small. Second, even among
patients who have apparent "suppression" of plasma viremia based on
less sensitive assays (thresholds between 50 and 200 copies/ml), a
subset will not have achieved complete suppression of viral
replication. In these cases, the reservoir of replication-competent
virus (whether in latently or productively infected cells) may in fact
be replenished by continued viral propagation, a conclusion consistent
with the demonstration of spliced RNA message and, in one case, changes
seen in pol sequences in subjects previously studied
(14).
Frequent testing for viremia with the most sensitive assays may permit
the identification of such individuals with incomplete suppression of
viral replication. The magnitude of risk for development of drug
resistance under these circumstances is not known. Consequently, it is
also not known whether, in the absence of demonstrable drug resistance,
additions to or substitutions in antiviral regimens are warranted.
However, these observations suggest that a less rapid initial decline
in viral load and the sporadic reappearance of low levels of detectable
viremia during therapy may signify incomplete suppression of viral
replication and that such assessments may be meaningful in evaluating
and comparing the potency of antiviral therapies, as has been proposed
by others (44). Finally, these findings indicate that any
measurement of the "clearance rate" of the reservoir of latently
infected T cells can only be accurately performed on the very best
characterized patients, employing, at a minimum, both sensitive and
frequent assays of residual RNA in blood to screen for incomplete
suppression of viral replication.
 |
ADDENDUM |
Following the submission of our manuscript, Zhang et al.
(50) reported the demonstration of viral replication in two
of eight patients with primary HIV infection treated with antiviral therapy who had achieved undetectable plasma RNA levels. In the same
issue, Furtado and colleagues (11b) reported the persistence of viral transcriptional activity in PBMC from treated, chronically infected patients with undetectable viral RNA, suggesting either ongoing replication or a transcriptionally active form of latent infection.
 |
ACKNOWLEDGMENTS |
We gratefully acknowledge the technical assistance of Linda Terry
and Nancy Keating; the administrative assistance of Sharon Wilcox,
Darica Smith, and Mark Biedermann; and the participation of the
patients from the San Diego Cohort of M035. We further acknowledge
helpful comments and constructive criticism from Steven Wolinsky, John
Guatelli, James Mullins, Allen Rodrigo, Gerald Learn, Jon Condra, and
Lisa Frenkel and the cooperation of Emilio Emini, Leslie Jonas, and
Robin Isaacs.
This work was supported by NIH grants AI 01361 and AI 43752 and a VA
Career Development Award (J.K.W.); NIH grants AI 27670, AI 38858, and
AI 29164 (D.D.R.); NIH grant AI 36214 (Center for AIDS Research)
(D.D.R. and D.J.L.); NIH grant RR 06555 (A.S.P.); Swiss National
Science Foundation Grant 84AD-046176 (H.F.G.); a Medical Research
Council grant (S.D.W.F.); an unrestricted educational grant from the
Merck Research Laboratories; and a grant from the Research Center for
AIDS and HIV Infection of the San Diego Veterans Affairs Medical Center.
 |
FOOTNOTES |
*
Corresponding author. Mailing address: Stein Clinical
Research Building No. 326, University of California, San Diego, 9500 Gilman Dr., La Jolla, CA 92093-0679. Phone: (619) 552-8585, ext. 7193. Fax: (619) 552-7445. E-mail: j2wong{at}ucsd.edu.
Present address: University of Zürich, Zürich, Switzerland.
 |
REFERENCES |
| 1.
|
Adachi, A.,
H. E. Gendelman,
S. Koenig,
T. Folks,
R. Willey,
A. Rabson, and M. A. Martin.
1986.
Production of acquired immunodeficiency syndrome-associated retrovirus in human and nonhuman cells transfected with an infectious molecular clone.
J. Virol.
59:284-291[Abstract/Free Full Text].
|
| 1a.
|
Balfe, P.,
P. Simmonds,
C. A. Ludlam,
J. O. Bishop, and A. J. Leigh-Brown.
1990.
Concurrent evolution of human immunodeficiency virus type 1 in patients infected from the same source: rate of sequence change and low frequency of inactivating mutations.
J. Virol.
64:6221-6233[Abstract/Free Full Text].
|
| 2.
|
Cavert, W.,
D. W. Notermans,
K. Staskus,
S. Wietgrefe,
M. Zupancic,
K. Gebhard,
K. Henry,
Z. Q. Zhang,
R. Mills,
H. McDade,
J. Goudsmit,
S. A. Danner, and A. T. Haase.
1997.
Kinetics of response in lymphoid tissues to antiretroviral therapy of HIV-1 infection.
Science
276:960-964[Abstract/Free Full Text].
|
| 3.
|
Chun, T. W.,
L. Stuyver,
S. B. Mizell,
L. A. Ehler,
J. M. Mican,
M. Baseler,
A. L. Lloyd,
M. A. Nowak, and A. S. Fauci.
1997.
Presence of an inducible HIV latent reservoir during highly active antiretroviral therapy.
Proc. Natl. Acad. Sci. USA
94:13193-13197[Abstract/Free Full Text].
|
| 4.
|
De Jong, J.-J.,
J. Goudsmit,
W. Keulen,
B. Klaver,
W. Krone,
M. Tersmette, and A. de Ronde.
1992.
Human immunodeficiency virus type 1 clones chimeric for the envelope V3 domain differ in syncytium formation and replication capacity.
J. Virol.
66:757-765[Abstract/Free Full Text].
|
| 5.
|
Delassus, S.,
R. Cheynier, and S. Wain-Hobson.
1991.
Evolution of human immunodeficiency virus type 1 nef and long terminal repeat sequences over 4 years in vivo and in vitro.
J. Virol.
65:225-231[Abstract/Free Full Text].
|
| 6.
|
Delwart, E. L.,
P. Heng,
H. W. Sheppard,
D. Wolpert,
A. U. Neumann,
B. Korber, and J. I. Mullins.
1997.
Slower evolution of human immunodeficiency virus type 1 quasispecies during progression to AIDS.
J. Virol.
71:7498-7508[Abstract].
|
| 7.
|
Felsenstein, J.
1993.
Phylogenetic inference package (Phylip) 3.5.
University of Washington, Seattle.
|
| 8.
| Felsenstein, J. Personal communication.
|
| 9.
|
Finzi, D.,
M. Hermankova,
T. Pierson,
L. Carruth,
C. Buck,
R. Chaisson,
T. Quinn,
K. Chadwick,
J. Margolick,
R. Brookmeyer,
J. Gallant,
M. Markowitz,
D. D. Ho,
D. Richman, and R. Siliciano.
1997.
Identification of a reservoir for HIV in patients on highly active antiretroviral therapy.
Science
278:1295-1300[Abstract/Free Full Text].
|
| 10.
|
Finzi, D.,
J. Blankson,
J. D. Siliciano,
J. B. Margolick,
K. Chadwick,
T. Pierson,
K. Smith,
J. Lisziewicz,
F. Lori,
C. Flexner,
T. C. Quinn,
R. E. Chaisson,
E. Rosenberg,
B. Walker,
S. Gange,
J. Gallant, and R. F. Siliciano.
1999.
Latent infection of CD4 T cells provides a mechanism for lifelong persistence of HIV-1, even in patients on effective combination therapy.
Nat. Med.
5:512-517[Medline].
|
| 11.
|
Fouchier, R. A. M.,
M. Groenink,
N. A. Kootstra,
M. Tersmette,
H. G. Huisman,
F. Miedema, and H. Schuitemaker.
1992.
Phenotype-associated sequence variation in the third variable domain of the human immunodeficiency virus type 1 gp120 molecule.
J. Virol.
66:3183-3187[Abstract/Free Full Text].
|
| 11a.
| Frost, S. D. W., et al. Unpublished data.
|
| 11b.
|
Furtado, M. R.,
D. S. Callaway,
J. P. Phair,
K. Kunstman,
J. L. Stanton,
C. A. Macken,
A. S. Perelson, and S. M. Wolinsky.
1999.
Persistence of HIV-1 transcription in peripheral blood mononuclear cells in patients receiving potent antiretroviral therapy.
N. Engl. J. Med.
340:1614-1622[Abstract/Free Full Text].
|
| 12.
|
Gulick, R.,
J. W. Mellors,
D. Havlir,
J. Eron,
C. Gonzalez,
D. McMahon,
L. Jonas,
A. Meibohm,
D. Holder,
W. A. Schleif,
J. H. Condra,
E. Emini,
R. Isaacs,
J. Chodakewitz, and D. Richman.
1998.
Simultaneous vs. sequential initiation of therapy with indinavir, zidovudine, and lamivudine for HIV-1 infection: 100 week follow-up.
JAMA
280:35-41[Abstract/Free Full Text].
|
| 13.
|
Gulick, R. M.,
J. W. Mellors,
D. Havlir,
J. J. Eron,
C. Gonzalez,
D. McMahon,
D. D. Richman,
F. T. Valentine,
P. Deutsch,
A. Meibohm,
D. Holder,
W. A. Schleif,
J. H. Condra,
E. A. Emini, and J. A. Chodakewitz.
1997.
Treatment with a combination of indinavir, zidovudine and lamivudine in HIV-infected adults with prior antiretroviral use.
N. Engl. J. Med.
337:734-739[Abstract/Free Full Text].
|
| 14.
|
Günthard, H. F.,
J. K. Wong,
C. C. Ignacio,
J. C. Guatelli,
N. L. Riggs,
D. Havlir, and D. D. Richman.
1998.
Human immunodeficiency virus replication and genotypic resistance in blood and lymph nodes after a year of potent antiretroviral therapy.
J. Virol.
72:2422-2428[Abstract/Free Full Text].
|
| 15.
|
Hahn, B. H.,
G. M. Shaw,
M. E. Taylor,
R. R. Redfield,
P. D. Markham,
S. Z. Salahuddin,
F. Wong-Staal,
R. C. Gallo,
E. S. Parks, and W. P. Parks.
1986.
Genetic variation in HTLV-III/LAV over time in patients with AIDS or at risk for AIDS.
Science
232:1548-1553[Abstract/Free Full Text].
|
| 16.
|
Hillis, D.,
J. Huelsenbeck, and C. Cunningham.
1994.
Application and accuracy of molecular phylogenies.
Science
264:671-675[Abstract/Free Full Text].
|
| 17.
|
Hudson, R. R.
1990.
Gene genealogies and the coalescent process, p. 1-44.
In
D. Futuyma, and J. Antonovics (ed.), Oxford surveys in evolutionary biology. Oxford University Press, Oxford, United Kingdom.
|
| 18.
|
Kingman, J. F. C.
1982.
On the genealogy of large populations, p. 27-43.
In
J. Gani, and E. J. Hannan (ed.), Applied probability trust.Sheffield, England.
|
| 19.
|
Korber, B. T.,
E. E. Allen,
A. D. Farmer, and G. Myers.
1995.
Heterogeneity of HIV-1 and HIV-2.
AIDS
9(Suppl. A):S5-S18.
|
| 20.
|
Korber, B. T. M.,
K. J. Kunstman,
B. K. Patterson,
M. Furtado,
M. M. McEvilly,
R. Levy, and S. M. Wolinsky.
1994.
Genetic differences between blood- and brain-derived viral sequences from human immunodeficiency virus type 1-infected patients: evidence of conserved elements in the V3 region of the envelope protein of brain-derived sequences.
J. Virol.
68:7467-7481[Abstract/Free Full Text].
|
| 21.
|
Korber, B. T. M.,
K. MacInnes,
R. F. Smith, and G. Myers.
1994.
Mutational trends in V3 loop protein sequences observed in different genetic lineages of human immunodeficiency virus type 1.
J. Virol.
68:6730-6744[Abstract/Free Full Text].
|
| 22.
|
Kuiken, C.,
G. Zwart,
E. Baan,
R. Coutinho,
J. van den Hoek, and J. Goudsmit.
1993.
Increasing antigenic and genetic diversity of the V3 variable domain of the HIV-1 envelope protein in the course of the AIDS epidemic.
Proc. Natl. Acad. Sci. USA
90:9061-9065[Abstract/Free Full Text].
|
| 23.
|
Kusumi, K.,
B. Conway,
S. Cunningham,
A. Berson,
C. Evans,
A. K. N. Iversen,
D. Colvin,
M. V. Gallo,
S. Coutre,
E. G. Shpaer,
D. V. Faulkner,
A. deRonde,
S. Volkman,
C. Williams,
M. S. Hirsch, and J. I. Mullins.
1992.
Human immunodeficiency virus type 1 envelope gene structure and diversity in vivo and after cocultivation in vitro.
J. Virol.
66:875-885[Abstract/Free Full Text].
|
| 24.
|
Learn, G. H.,
B. T. Korber,
B. Foley,
B. Hahn,
S. M. Wolinsky, and J. I. Mullins.
1996.
Maintaining the integrity of the HIV sequence databases.
J. Virol.
70:5720-5730[Abstract/Free Full Text].
|
| 25.
|
Leitner, T.,
D. Escanilla,
C. Franzen,
M. Uhlen, and J. Albert.
1996.
Accurate reconstruction of a known HIV-1 transmission history by phylogenetic tree analysis.
Proc. Natl. Acad. Sci. USA
93:10864-10869[Abstract/Free Full Text].
|
| 26.
|
Leitner, T.,
S. Kumar, and J. Albert.
1997.
Tempo and mode of nucleotide substitutions in gag and env gene fragments in human immunodeficiency virus type 1 populations with known transmission history.
J. Virol.
71:4761-4770[Abstract].
|
| 27.
|
Lukashov, V. V.,
C. Kuiken,
D. Vlahov,
R. Coutinho, and J. Goudsmit.
1996.
Evidence for HIV-1 strains of U.S. intravenous drug users as founders of AIDS epidemic among intravenous drug users in northern Europe.
AIDS Res. Hum. Retroviruses
12:1179-1183[Medline].
|
| 28.
|
Michael, N. L.,
G. Chang,
P. K. Ehrenberg,
M. T. Vahey, and R. R. Redfield.
1993.
HIV-1 proviral genotypes from the peripheral blood mononuclear cells of an infected patient are differentially represented in expressed sequences.
J. AIDS
6:1073-1085.
|
| 29.
|
Nei, M., and L. Jin.
1989.
Variances of the average numbers of nucleotide substitutions within and between populations.
Mol. Biol. Evol.
6:290-300[Abstract].
|
| 30.
|
Olsen, G. J.,
H. Matsuda,
R. Hagstrom, and R. Overbeek.
1994.
FastDNAML 1.0.
University of Illinois, Urbana, and Argonne National Laboratory, Argonne, Ill.
|
| 30a.
|
Page, R. D. M.
1996.
Treeview: and application to display phylogenetic trees on personal computers.
Comput. Appl. Biosci.
12:357-358[Free Full Text].
|
| 31.
|
Perelson, A. S.,
P. Essunger,
Y. Cao,
M. Vesanen,
A. Hurley,
K. Saksela,
M. Markowitz, and D. D. Ho.
1997.
Decay characteristics of HIV-1 infected compartments during combination therapy.
Nature
387:188-191[Medline].
|
| 32.
|
Perelson, A. S.,
A. U. Neumann,
M. Markowitz,
J. M. Leonard, and D. D. Ho.
1996.
HIV-1 dynamics in vivo: virion clearance rate, infected cell lifetime, and viral generation time.
Science
271:1582-1586[Abstract].
|
| 33.
|
Power, C.,
J. C. McArthur,
R. T. Johnson,
D. E. Griffin,
J. D. Glass,
S. Perryman, and B. Chesebro.
1994.
Demented and nondemented patients with AIDS differ in brain-derived human immunodeficiency virus type 1 envelope sequences.
J. Virol.
68:4643-4649[Abstract/Free Full Text].
|
| 34.
|
Richman, D. D., and S. A. Bozzette.
1994.
The impact of syncytium-inducing phenotype of human immunodeficiency virus on disease progression.
J. Infect. Dis.
169:968-974[Medline].
|
| 35.
|
Rodrigo, A.
1998.
Coalescent approaches to HIV population genetics, p. 54.
In
Proceedings of the 5th Annual HIV Dynamics and Evolution Conference.
|
| 36.
|
Saag, M. S.,
B. H. Hahn,
J. Gibbons,
Y. Li,
E. S. Parks,
W. P. Parks, and G. M. Shaw.
1988.
Extensive variation of human immunodeficiency virus type-1 in vivo.
Nature
334:440-444[Medline].
|
| 37.
|
Shankarappa, R.,
P. Gupta,
G. Learn,
A. Rodrigo,
C. R. Rinaldo, Jr.,
M. C. Gorry,
J. I. Mullins,
P. L. Nara, and G. D. Ehrlich.
1998.
Evolution of HIV-1 envelope sequences in infected individuals with differing disease progression profiles.
Virology
241:251-259[Medline].
|
| 38.
|
Simmonds, P.,
P. Balfe,
C. A. Ludlam,
J. O. Bishop, and A. J. Leigh-Brown.
1990.
Analysis of sequence diversity in hypervariable regions of the external glycoprotein of human immunodeficiency virus type 1.
J. Virol.
64:5840-5850[Abstract/Free Full Text].
|
| 39.
|
Simmonds, P.,
L. Q. Zhang,
F. McOmish,
P. Balfe,
C. A. Ludlam, and A. J. Leigh-Brown.
1991.
Discontinuous sequence change of human immunodeficiency virus (HIV) type 1 env sequences in plasma viral and lymphocyte-associated proviral populations in vivo: implications for models of HIV pathogenesis.
J. Virol.
65:6266-6276[Abstract/Free Full Text].
|
| 40.
|
Speck, R. F.,
K. Wehrly,
E. J. Platt,
R. E. Atchison,
I. F. Charo,
D. Kabat,
B. Chesebro, and M. A. Goldsmith.
1997.
Selective employment of chemokine receptors as human immunodeficiency virus type 1 coreceptors determined by individual amino acids in the envelope V3 loop.
J. Virol.
71:7136-7139[Abstract].
|
| 41.
|
Tajima, F.
1983.
Evolutionary relationship of DNA sequences in finite populations.
Genetics
105:437-460[Abstract/Free Full Text].
|
| 42.
|
van't Wout, A.,
L. Ran,
C. Kuiken,
N. Kootstra,
S. Pals, and H. Schuitemaker.
1998.
Analysis of the temporal relationship between human immunodeficiency virus type 1 quasispecies in sequential blood samples and various organs obtained at autopsy.
J. Virol.
72:488-496[Abstract/Free Full Text].
|
| 43.
|
Wainberg, M. A.,
W. C. Drosopoulous,
H. Salomon,
M. Hsu,
G. Borkow,
M. Parniak,
Z. Gu,
Q. Song,
J. Manne, and S. Islam.
1996.
Enhanced fidelity of 3TC-selected mutant HIV-1 reverse transcriptase.
Science
271:1282-1285[Abstract].
|
| 44.
|
Weverling, G. J.,
J. M. A. Lange,
S. Jurriaans,
J. M. Prins,
V. Lukashov,
D. W. Notermans,
M. Roos,
H. Schuitemaker,
R. M. W. Hoetelmans,
S. V. Danner,
J. Goudsmit, and F. de Wolf.
1998.
Alternative multidrug regimen provides improved suppression of HIV-1 replication over triple therapy.
AIDS
12:F117-F122[Medline].
|
| 45.
|
Wolinsky, S. M.,
B. T. Korber,
A. U. Neumann,
M. Daniels,
K. J. Kuntsman,
A. J. Whetsell,
M. R. Furtado,
Y. Cao,
D. D. Ho,
J. T. Safrit, et al.
1996.
Adaptive evolution of human immunodeficiency virus-type 1 during the natural course of infection.
Science
272:537-542[Abstract].
|
| 46.
|
Wong, J. K.,
H. F. Gunthard,
D. V. Havlir,
Z. Q. Zhang,
A. T. Haase,
C. C. Ignacio,
S. Kwok,
E. A. Emini, and D. D. Richman.
1997.
Reduction of HIV-1 in blood and lymph nodes following potent anti-retroviral therapy and the virologic correlates of treatment failure.
Proc. Natl. Acad. Sci. USA
94:12574-12579[Abstract/Free Full Text].
|
| 47.
|
Wong, J. K.,
M. Hezareh,
H. F. Gunthard,
D. Havlir,
C. C. Ignacio,
C. Spina, and D. D. Richman.
1997.
Recovery of replication competent HIV despite prolonged suppression of plasma viremia.
Science
278:1291-1295[Abstract/Free Full Text].
|
| 48.
|
Wong, J. K.,
C. C. Ignacio,
F. Torriani,
D. Havlir,
N. J. S. Fitch, and D. D. Richman.
1997.
In vivo compartmentalization of HIV: evidence from the examination of pol sequences from autopsy tissues.
J. Virol.
70:2059-2071[Abstract].
|
| 49.
|
Zhang, H.,
G. Dornadula,
M. Beumont,
L. Livornese, Jr.,
B. Van Uitert,
K. Henning, and R. J. Pomerantz.
1998.
Human Immunodeficiency Virus type 1 in the semen of men receiving highly active antiretroviral therapy.
N. Engl. J. Med.
339:1803-1809[Abstract/Free Full Text].
|
| 50.
|
Zhang, L.,
B. Ramratnam,
K. Tenner-Racz,
Y. He,
M. Vesanen,
S. Lewin,
A. Talal,
P. Racz,
A. S. Perelson,
B. T. Korber,
M. Markowitz, and D. D. Ho.
1999.
Quantifying residual HIV-1 replication in patients receiving combination antiretroviral therapy.
N. Engl. J. Med.
340:1605-1613[Abstract/Free Full Text].
|
Journal of Virology, November 1999, p. 9404-9412, Vol. 73, No. 11
0022-538X/99/$04.00+0
Copyright © 1999, American Society for Microbiology. All rights reserved.
This article has been cited by other articles:
-
Di Mascio, M., Srinivasula, S., Bhattacharjee, A., Cheng, L., Martiniova, L., Herscovitch, P., Lertora, J., Kiesewetter, D.
(2009). Antiretroviral Tissue Kinetics: In Vivo Imaging Using Positron Emission Tomography. Antimicrob. Agents Chemother.
53: 4086-4095
[Abstract]
[Full Text]
-
Richman, D. D., Margolis, D. M., Delaney, M., Greene, W. C., Hazuda, D., Pomerantz, R. J.
(2009). The Challenge of Finding a Cure for HIV Infection. Science
323: 1304-1307
[Abstract]
[Full Text]
-
Joos, B., Fischer, M., Kuster, H., Pillai, S. K., Wong, J. K., Boni, J., Hirschel, B., Weber, R., Trkola, A., Gunthard, H. F., The Swiss HIV Cohort Study,
(2008). HIV rebounds from latently infected cells, rather than from continuing low-level replication. Proc. Natl. Acad. Sci. USA
105: 16725-16730
[Abstract]
[Full Text]
-
Joos, B., Trkola, A., Kuster, H., Aceto, L., Fischer, M., Stiegler, G., Armbruster, C., Vcelar, B., Katinger, H., Gunthard, H. F.
(2006). Long-Term Multiple-Dose Pharmacokinetics of Human Monoclonal Antibodies (MAbs) against Human Immunodeficiency Virus Type 1 Envelope gp120 (MAb 2G12) and gp41 (MAbs 4E10 and 2F5).. Antimicrob. Agents Chemother.
50: 1773-1779
[Abstract]
[Full Text]
-
Tobin, N. H., Learn, G. H., Holte, S. E., Wang, Y., Melvin, A. J., McKernan, J. L., Pawluk, D. M., Mohan, K. M., Lewis, P. F., Mullins, J. I., Frenkel, L. M.
(2005). Evidence that Low-Level Viremias during Effective Highly Active Antiretroviral Therapy Result from Two Processes: Expression of Archival Virus and Replication of Virus. J. Virol.
79: 9625-9634
[Abstract]
[Full Text]
-
Joos, B., Trkola, A., Fischer, M., Kuster, H., Rusert, P., Leemann, C., Boni, J., Oxenius, A., Price, D. A., Phillips, R. E., Wong, J. K., Hirschel, B., Weber, R., Gunthard, H. F., the Swiss HIV Cohort Study,
(2005). Low Human Immunodeficiency Virus Envelope Diversity Correlates with Low In Vitro Replication Capacity and Predicts Spontaneous Control of Plasma Viremia after Treatment Interruptions. J. Virol.
79: 9026-9037
[Abstract]
[Full Text]
-
Monie, D., Simmons, R. P., Nettles, R. E., Kieffer, T. L., Zhou, Y., Zhang, H., Karmon, S., Ingersoll, R., Chadwick, K., Zhang, H., Margolick, J. B., Quinn, T. C., Ray, S. C., Wind-Rotolo, M., Miller, M., Persaud, D., Siliciano, R. F.
(2005). A Novel Assay Allows Genotyping of the Latent Reservoir for Human Immunodeficiency Virus Type 1 in the Resting CD4+ T Cells of Viremic Patients. J. Virol.
79: 5185-5202
[Abstract]
[Full Text]
-
Sharkey, M., Triques, K., Kuritzkes, D. R., Stevenson, M.
(2005). In Vivo Evidence for Instability of Episomal Human Immunodeficiency Virus Type 1 cDNA. J. Virol.
79: 5203-5210
[Abstract]
[Full Text]
-
Casazza, J. P., Betts, M. R., Hill, B. J., Brenchley, J. M., Price, D. A., Douek, D. C., Koup, R. A.
(2005). Immunologic Pressure within Class I-Restricted Cognate Human Immunodeficiency Virus Epitopes during Highly Active Antiretroviral Therapy. J. Virol.
79: 3653-3663
[Abstract]
[Full Text]
-
Nettles, R. E., Kieffer, T. L., Kwon, P., Monie, D., Han, Y., Parsons, T., Cofrancesco, J. Jr, Gallant, J. E., Quinn, T. C., Jackson, B., Flexner, C., Carson, K., Ray, S., Persaud, D., Siliciano, R. F.
(2005). Intermittent HIV-1 Viremia (Blips) and Drug Resistance in Patients Receiving HAART. JAMA
293: 817-829
[Abstract]
[Full Text]
-
Strain, M. C., Letendre, S., Pillai, S. K., Russell, T., Ignacio, C. C., Gunthard, H. F., Good, B., Smith, D. M., Wolinsky, S. M., Furtado, M., Marquie-Beck, J., Durelle, J., Grant, I., Richman, D. D., Marcotte, T., McCutchan, J. A., Ellis, R. J., Wong, J. K.
(2005). Genetic Composition of Human Immunodeficiency Virus Type 1 in Cerebrospinal Fluid and Blood without Treatment and during Failing Antiretroviral Therapy. J. Virol.
79: 1772-1788
[Abstract]
[Full Text]
-
Di Mascio, M., Markowitz, M., Louie, M., Hurley, A., Hogan, C., Simon, V., Follmann, D., Ho, D. D., Perelson, A. S.
(2004). Dynamics of Intermittent Viremia during Highly Active Antiretroviral Therapy in Patients Who Initiate Therapy during Chronic versus Acute and Early Human Immunodeficiency Virus Type 1 Infection. J. Virol.
78: 10566-10573
[Abstract]
[Full Text]
-
Wightman, F., Walters, T., Ayres, A., Bowden, S., Bartholomeusz, A., Lau, D., Locarnini, S., Lewin, S. R.
(2004). Comparison of Sequence Analysis and a Novel Discriminatory Real-Time PCR Assay for Detection and Quantification of Lamivudine-Resistant Hepatitis B Virus Strains. J. Clin. Microbiol.
42: 3809-3812
[Abstract]
[Full Text]
-
Scott, P., Arnold, E., Evans, B., Pozniak, A., Moyle, G., Shahmenesh, M., White, D., Shirley, J., Cane, P., Pillay, D.
(2004). Surveillance of HIV antiretroviral drug resistance in treated individuals in England: 1998-2000. J Antimicrob Chemother
53: 469-473
[Abstract]
[Full Text]
-
Bello, G., Casado, C., Garcia, S., Rodriguez, C., del Romero, J., Lopez-Galindez, C.
(2004). Co-existence of recent and ancestral nucleotide sequences in viral quasispecies of human immunodeficiency virus type 1 patients. J. Gen. Virol.
85: 399-407
[Abstract]
[Full Text]
-
Persaud, D., Siberry, G. K., Ahonkhai, A., Kajdas, J., Monie, D., Hutton, N., Watson, D. C., Quinn, T. C., Ray, S. C., Siliciano, R. F.
(2004). Continued Production of Drug-Sensitive Human Immunodeficiency Virus Type 1 in Children on Combination Antiretroviral Therapy Who Have Undetectable Viral Loads. J. Virol.
78: 968-979
[Abstract]
[Full Text]
-
Rodrigo, A. G., Goode, M., Forsberg, R., Ross, H. A., Drummond, A.
(2003). Inferring Evolutionary Rates Using Serially Sampled Sequences from Several Populations. Mol Biol Evol
20: 2010-2018
[Abstract]
[Full Text]
-
Di Mascio, M., Markowitz, M., Louie, M., Hogan, C., Hurley, A., Chung, C., Ho, D. D., Perelson, A. S.
(2003). Viral Blip Dynamics during Highly Active Antiretroviral Therapy. J. Virol.
77: 12165-12172
[Abstract]
[Full Text]
-
Havlir, D. V., Strain, M. C., Clerici, M., Ignacio, C., Trabattoni, D., Ferrante, P., Wong, J. K.
(2003). Productive Infection Maintains a Dynamic Steady State of Residual Viremia in Human Immunodeficiency Virus Type 1-Infected Persons Treated with Suppressive Antiretroviral Therapy for Five Years. J. Virol.
77: 11212-11219
[Abstract]
[Full Text]
-
Frenkel, L. M., Wang, Y., Learn, G. H., McKernan, J. L., Ellis, G. M., Mohan, K. M., Holte, S. E., De Vange, S. M., Pawluk, D. M., Melvin, A. J., Lewis, P. F., Heath, L. M., Beck, I. A., Mahalanabis, M., Naugler, W. E., Tobin, N. H., Mullins, J. I.
(2003). Multiple Viral Genetic Analyses Detect Low-Level Human Immunodeficiency Virus Type 1 Replication during Effective Highly Active Antiretroviral Therapy. J. Virol.
77: 5721-5730
[Abstract]
[Full Text]
-
Strain, M. C., Gunthard, H. F., Havlir, D. V., Ignacio, C. C., Smith, D. M., Leigh-Brown, A. J., Macaranas, T. R., Lam, R. Y., Daly, O. A., Fischer, M., Opravil, M., Levine, H., Bacheler, L., Spina, C. A., Richman, D. D., Wong, J. K.
(2003). Heterogeneous clearance rates of long-lived lymphocytes infected with HIV: Intrinsic stability predicts lifelong persistence. Proc. Natl. Acad. Sci. USA
100: 4819-4824
[Abstract]
[Full Text]
-
Persaud, D., Zhou, Y., Siliciano, J. M., Siliciano, R. F.
(2003). Latency in Human Immunodeficiency Virus Type 1 Infection: No Easy Answers. J. Virol.
77: 1659-1665
[Full Text]
-
Ruff, C. T., Ray, S. C., Kwon, P., Zinn, R., Pendleton, A., Hutton, N., Ashworth, R., Gange, S., Quinn, T. C., Siliciano, R. F., Persaud, D.
(2002). Persistence of Wild-Type Virus and Lack of Temporal Structure in the Latent Reservoir for Human Immunodeficiency Virus Type 1 in Pediatric Patients with Extensive Antiretroviral Exposure. J. Virol.
76: 9481-9492
[Abstract]
[Full Text]
-
van Rij, R. P., Visser, J. A., van Praag, R. M. E., Rientsma, R., Prins, J. M., Lange, J. M. A., Schuitemaker, H.
(2002). Both R5 and X4 Human Immunodeficiency Virus Type 1 Variants Persist during Prolonged Therapy with Five Antiretroviral Drugs. J. Virol.
76: 3054-3058
[Abstract]
[Full Text]
-
Brooks, D. G., Zack, J. A.
(2002). Effect of Latent Human Immunodeficiency Virus Infection on Cell Surface Phenotype. J. Virol.
76: 1673-1681
[Abstract]
[Full Text]
-
Zhu, T., Muthui, D., Holte, S., Nickle, D., Feng, F., Brodie, S., Hwangbo, Y., Mullins, J. I., Corey, L.
(2002). Evidence for Human Immunodeficiency Virus Type 1 Replication In Vivo in CD14+ Monocytes and Its Potential Role as a Source of Virus in Patients on Highly Active Antiretroviral Therapy. J. Virol.
76: 707-716
[Abstract]
[Full Text]
-
Ortiz, G. M., Hu, J., Goldwitz, J. A., Chandwani, R., Larsson, M., Bhardwaj, N., Bonhoeffer, S., Ramratnam, B., Zhang, L., Markowitz, M. M., Nixon, D. F.
(2002). Residual Viral Replication during Antiretroviral Therapy Boosts Human Immunodeficiency Virus Type 1-Specific CD8+ T-Cell Responses in Subjects Treated Early after Infection. J. Virol.
76: 411-415
[Abstract]
[Full Text]
-
Lawn, S. D., Butera, S. T., Folks, T. M.
(2001). Contribution of Immune Activation to the Pathogenesis and Transmission of Human Immunodeficiency Virus Type 1 Infection. Clin. Microbiol. Rev.
14: 753-777
[Abstract]
[Full Text]
-
Havlir, D. V., Bassett, R., Levitan, D., Gilbert, P., Tebas, P., Collier, A. C., Hirsch, M. S., Ignacio, C., Condra, J., Gunthard, H. F., Richman, D. D., Wong, J. K.
(2001). Prevalence and Predictive Value of Intermittent Viremia With Combination HIV Therapy. JAMA
286: 171-179
[Abstract]
[Full Text]
-
Hermankova, M., Ray, S. C., Ruff, C., Powell-Davis, M., Ingersoll, R., D'Aquila, R. T., Quinn, T. C., Siliciano, J. D., Siliciano, R. F., Persaud, D.
(2001). HIV-1 Drug Resistance Profiles in Children and Adults With Viral Load of <50 Copies/mL Receiving Combination Therapy. JAMA
286: 196-207
[Abstract]
[Full Text]
-
Deeks, S. G.
(2001). Durable HIV Treatment Benefit Despite Low-Level Viremia: Reassessing Definitions of Success or Failure. JAMA
286: 224-226
[Full Text]
-
Martinez-Picado, J., DePasquale, M. P., Kartsonis, N., Hanna, G. J., Wong, J., Finzi, D., Rosenberg, E., Gunthard, H. F., Sutton, L., Savara, A., Petropoulos, C. J., Hellmann, N., Walker, B. D., Richman, D. D., Siliciano, R., D'Aquila, R. T.
(2000). Antiretroviral resistance during successful therapy of HIV type 1 infection. Proc. Natl. Acad. Sci. USA
97: 10948-10953
[Abstract]
[Full Text]
-
Crowe, S. M., Sonza, S.
(2000). HIV-1 can be recovered from a variety of cells including peripheral blood monocytes of patients receiving highly active antiretroviral therapy: a further obstacle to eradication. J. Leukoc. Biol.
68: 345-350
[Abstract]
[Full Text]
-
Korber, B., Muldoon, M., Theiler, J., Gao, F., Gupta, R., Lapedes, A., Hahn, B. H., Wolinsky, S., Bhattacharya, T.
(2000). Timing the Ancestor of the HIV-1 Pandemic Strains. Science
288: 1789-1796
[Abstract]
[Full Text]
-
Kuster, H., Opravil, M., Ott, P., Schlaepfer, E., Fischer, M., Gunthard, H. F., Luthy, R., Weber, R., Cone, R. W.
(2000). Treatment-Induced Decline of Human Immunodeficiency Virus-1 p24 and HIV-1 RNA in Lymphoid Tissue of Patients with Early Human Immunodeficiency Virus-1 Infection. Am. J. Pathol.
156: 1973-1986
[Abstract]
[Full Text]
-
Carpenter, C. C. J., Cooper, D. A., Fischl, M. A., Gatell, J. M., Gazzard, B. G., Hammer, S. M., Hirsch, M. S., Jacobsen, D. M., Katzenstein, D. A., Montaner, J. S. G., Richman, D. D., Saag, M. S., Schechter, M., Schooley, R. T., Thompson, M. A., Vella, S., Yeni, P. G., Volberding, P. A.
(2000). Antiretroviral Therapy in Adults: Updated Recommendations of the International AIDS Society-USA Panel. JAMA
283: 381-390
[Abstract]
[Full Text]