|Daily Sampling of an HIV-1 Patient with Slowly Progressing Disease Displays Persistence of Multiple env Subpopulations Consistent with Neutrality.|
|Jump to Full Text|
|PMID: 21829600 Owner: NLM Status: In-Data-Review|
|The molecular evolution of HIV-1 is characterized by frequent substitutions, indels and recombination events. In addition, a HIV-1 population may adapt through frequency changes of its variants. To reveal such population dynamics we analyzed HIV-1 subpopulation frequencies in an untreated patient with stable, low plasma HIV-1 RNA levels and close to normal CD4+ T-cell levels. The patient was intensively sampled during a 32-day period as well as approximately 1.5 years before and after this period (days -664, 1, 2, 3, 11, 18, 25, 32 and 522). 77 sequences of HIV-1 env (approximately 3100 nucleotides) were obtained from plasma by limiting dilution with 7-11 sequences per time point, except day -664. Phylogenetic analysis using maximum likelihood methods showed that the sequences clustered in six distinct subpopulations. We devised a method that took into account the relatively coarse sampling of the population. Data from days 1 through 32 were consistent with constant within-patient subpopulation frequencies. However, over longer time periods, i.e. between days 1…32 and 522, there were significant changes in subpopulation frequencies, which were consistent with evolutionarily neutral fluctuations. We found no clear signal of natural selection within the subpopulations over the study period, but positive selection was evident on the long branches that connected the subpopulations, which corresponds to >3 years as the subpopulations already were established when we started the study. Thus, selective forces may have been involved when the subpopulations were established. Genetic drift within subpopulations caused by de novo substitutions could be resolved after approximately one month. Overall, we conclude that subpopulation frequencies within this patient changed significantly over a time period of 1.5 years, but that this does not imply directional or balancing selection. We show that the short-term evolution we study here is likely representative for many patients of slow and normal disease progression.|
|Helena Skar; Ryan N Gutenkunst; Karin Wilbe Ramsay; Annette Alaeus; Jan Albert; Thomas Leitner|
|Type: Journal Article Date: 2011-08-02|
|Title: PloS one Volume: 6 ISSN: 1932-6203 ISO Abbreviation: PLoS ONE Publication Date: 2011|
|Created Date: 2011-08-10 Completed Date: - Revised Date: -|
Medline Journal Info:
|Nlm Unique ID: 101285081 Medline TA: PLoS One Country: United States|
|Languages: eng Pagination: e21747 Citation Subset: IM|
|Department of Virology, Swedish Institute for Infectious Disease Control, Solna, Sweden.|
|APA/MLA Format Download EndNote Download BibTex|
Journal ID (nlm-ta): PLoS One
Journal ID (publisher-id): plos
Journal ID (pmc): plosone
Publisher: Public Library of Science, San Francisco, USA
This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Received Day: 4 Month: 4 Year: 2011
Accepted Day: 6 Month: 6 Year: 2011
collection publication date: Year: 2011
Electronic publication date: Day: 2 Month: 8 Year: 2011
Volume: 6 Issue: 8
E-location ID: e21747
PubMed Id: 21829600
Publisher Id: PONE-D-11-06212
|Daily Sampling of an HIV-1 Patient with Slowly Progressing Disease Displays Persistence of Multiple env Subpopulations Consistent with Neutrality Alternate Title:Neutral Subpopulation Fluctuations in HIV-1|
|Ryan N. Gutenkunst4|
|Karin Wilbe Ramsay12|
1Department of Virology, Swedish Institute for Infectious Disease Control, Solna, Sweden
2Department of Microbiology, Tumor and Cell Biology, Karolinska Institute, Stockholm, Sweden
3Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, New Mexico, United States of America
4Department of Molecular and Cellular Biology, University of Arizona, Tucson, Arizona, United States of America
5Department of Medicine, Karolinska Institute, Stockholm, Sweden
|Correspondence: * E-mail: email@example.com
Contributed by footnote: Conceived and designed the experiments: TL HS KWR JA AA. Performed the experiments: KWR AA. Analyzed the data: HS RNG TL. Wrote the paper: TL HS RNG JA. Conceived and designed the population genetics modeling: RNG TL.
The HIV-1 envelope gene (env) displays the largest genetic diversity in the HIV-1 genome. The evolutionary rate (nucleotide substitution rate) of env is affected by the strength of the pressure of the immune system ,  so that both the immune pressure and the evolutionary rate are higher during the chronic, asymptomatic phase than during end-stage disease. Similarly, the immune pressure in long-term non-progressors lasts longer and is often stronger than in typical patients. Thus, HIV-1 genetic evolution in env during the chronic disease stage has been characterized by positive selection for escape mutants due to continuous immune surveillance , , , . However, other studies have found HIV-1 evolution during chronic infection to be consistent with a neutral model of evolution, characterized by small effective population sizes (Ne) strongly influenced by random genetic drift , . Whether the mutation process will be deterministic or stochastic is generally believed to be dependent on the population size. Deterministic models assume an infinite population size, which given the large amount of HIV-1 particles produced daily in an infected individual (1010 virions/day) is not unreasonable . However, it has been proposed that the Ne of HIV-1 during chronic infection is several orders of magnitude lower , , , , which would suggest that stochastic processes could influence HIV-1 evolution. To date, a few models have tried to unify the estimated small Ne and the strong positive selection believed to act on HIV during chronic infection. A meta-population model, where a large collection of small subpopulations is subject to frequent migration, extinction, and recolonization, was shown to agree with the low effective population sizes seen in chronic HIV infection . Another example is a combination of both directional and neutral forces acting on the HIV population, where random genetic drift of neutral mutations predominates combined with brief episodes of directional selection . A combination of the two, where the meta-population model and selective sweeps both are factors that act together to reduce the intra-host effective population size of HIV-1 has been proposed to be the most likely explanation of the reduced Ne. Thus, it is still unclear how HIV diversity is affected by selection in an infected individual, and furthermore on which time scale selection operates.
Here we compare short-term (days, weeks, months) and long-term (years) HIV-1 evolution in a treatment naïve, asymptomatic patient with low plasma HIV-1 RNA levels (viral load) and fluctuating, often close to normal CD4+ T-lymphocyte (CD4) counts. In patients like this the immune system generally puts a strong pressure on the virus for a longer time than in typical patients that, in the absence of antiretroviral drugs, develop AIDS quicker. We find that multiple distinct subpopulations persist over years, but that their frequencies fluctuate over time. The fluctuations during the time period of days to months showed no significant signature of variable selection across sequence sites, and the fluctuations were consistent with a neutral model of evolution. Hence, we find no need for balancing selection to explain the persistence of the subpopulations over these time intervals. However, over the period of years, we could detect a signal of positive selection, especially at potential N-linked glycosylation sites (PNGS), which may have shaped the subpopulation structure. Finally, we show that it is important to correctly handle subpopulation fluctuations when using genetic distances to estimate the number of de novo mutations.
Seventy-seven individual virus sequences of approximately 3100 nucleotides covering vpu, env and the first half of nef were analyzed. The sequences were sampled by limiting dilution from plasma samples obtained at 9 different time points spanning a period of 3 years (Table 1). The limiting dilution sequencing methodology (aka. SGA and SGS) applied here ensures that PCR and sequencing artifacts are virtually absent in the sequences , . The majority of sequences were sampled during a time period of 32 days, where the first time point was called day 1. In addition, two samples were drawn approximately 1.5 years before (day −664) and 1.5 years after (day 522) the main sampling period. At each time point 7 to 11 sequences were generated with exception for the earliest time point from which only 3 sequences could be amplified. As this patient was a slow progressor with low virus load it was difficult to obtain additional sequences (see Methods). Only 2 sequences were identical (at day 18), whereas all other sequences were unique.
In total, 18 unique deleterious mutations were observed, including 9 nucleotide substitutions that caused stop codons and 9 that caused frame shifts. Because deleterious mutations are unlikely to survive to the next generation, this suggests a minimum rate of 7.67×10−5 deleterious substitutions per site per generation, i.e., in the same order of magnitude as other point mutations and recombination occur , , . In addition to point mutations, sequence −664.2 had two large deletions, one of 149 nucleotides (nts) in the beginning of gp120 and another of 435 nts in the end of gp41. Sequence 1.10 had two large deletions of 54 and 60 nts in the middle of gp120, and 522.1 had a large deletion of 48 nts in the middle of gp120. In addition, several smaller, mostly in-frame indels, were present, which thus resulted in amino acid insertions or deletions.
Putative recombinants were identified using the PHI-NNet test. Two sequences (s4.−664.3 and s5.522.10) were classified as putative recombinants within subpopulations and six sequences were identified as putative recombinants with ancestors derived from at least two subpopulations (s2. −664.2, s4. −664.1, s4.18.7, s4.32.9, s4.522.9, s5.1.7) (Figure S1). If these sequences were removed no recombination signal remained in the dataset according to the c-AIC criteria in a GARD Single Breakpoint analysis. To confirm that the identified putative recombinants were robustly identified, we performed 100 iterations of removal of 8 random sequences. None of these iterations rendered the dataset free from recombination signal according to the PHI-NNet test.
The general time reversible model with variable rates among sites and a proportion of invariable sites (GTR+G+I) was the best substitution model for our data according to a Modeltest analysis. This model was used to infer a maximum likelihood (ML) tree of the HIV population (Figure 1). The tree displayed six phylogenetic clades, designated as subpopulations s1 through s6, which were all supported by ML bootstrap values 61–100%. Independent of the inferred tree, and thus less affected by any remaining recombination signal, Hudson's population subdivision test supported all subpopulations except s2 at p[K*s]≤0.0005 (s2 p[K*s] = 0.0668). A majority of the subpopulations (4 of 5) persisted over the entire study period, if we excluded day −664 that was insufficiently sampled. Thus, at the last time point (day 522), representatives of subpopulations s3, s4, s5 and s6 were still present.
To test if the different subpopulations were under different potential selection pressures, we investigated dN/dS ratios on all sites in the sequences and on all branches in the phylogeny. Four dN/dS categories (0, 0.42, 1.24, 10000) were found to best explain the data according to the c-AIC criterion and the GAbranch model available in Hyphy to test lineage specific selection on branches (Figure S2). A majority (79%) of the branches in the tree fell into dN/dS categories 0.42 and 1.24 and of these 72% suggested positive selection, but there was no clear pattern of where in the tree they occurred. The deep branches that connected the subpopulations displayed selection in either direction, i.e., 0.42≤dN/dS≤1.24. Branches displaying either no synonymous or non-synonymous mutations (dN/dS categories 0 and 10000) occurred exclusively within the subpopulations, where the total number of mutations on most of the branches was very low.
In agreement with the branch analyses, codon models could not identify any sites under selection within the subpopulations, suggesting neutral evolution over the time of this study (days 1…522). Furthermore, there was clear evidence of variable selection pressure over sites when we analyzed all subpopulations together in a single phylogenetic tree (p<0.01, likelihood ratio tests with M0 vs. M3: df = 4, and M1a vs. M2a: df = 2). This indicates that individual sites may have been under selection when the subpopulations were established, i.e., over a time much longer than 522 days as the subpopulations already existed and were defined by relatively long branches at day 1 (see further below for an evaluation of how much branches grew over the study period). Interestingly, potential N-linked glycosylation sites (PNGS) were significantly overrepresented among positively selected sites in analyses with the variable selection model M3 (p<0.001, Chi-square test). Further, while both amino acid substitutions and PNGS replacements correlated well with positive selection strength, the response to positive selection was stronger on branches separating the subpopulations than on branches within subpopulations (Figure S3).
Overall, these results suggest that natural selection has little or no impact over short time periods (≤1.5 years), but over long time periods (>3 years) positively selected sites, especially those that involve PNGS, may contribute to the subpopulation structure.
Figure 2 shows the count at which the different subpopulations were observed at each time point. We were interested in understanding the persistence of subpopulations over time. Could these experimental data be expected under a neutral process, or would the data be better explained by selection, and in particular balancing selection?
We first asked whether the experimental data provided evidence that the frequencies of the different subpopulations were changing within the patient. To do so, we used a χ2 statistic to assess whether the frequencies observed at each day J were consistent with the frequencies observed for the prior days 1…J-1. To account for the relatively few samples per time point we accounted for the number of sequences per sample (to limit stochastic sampling effects) and pooled days together when no frequency differences were observed (to increase the power of our analysis). Hence, statistical significance was assessed by simulating both the multinomial sampling of day J's observation, and the inference of the within-patient frequencies for days 1…J-1. The resulting p-values for constant within-patient subpopulation frequencies are shown in Figure 2. From this analysis, we observed no statistically significant changes in within-patient subpopulation frequencies over the first 32 days, but at day 522 the frequency fluctuations became statistically significant (p<0.05). As detailed in Supporting Information, repeated analyses that excluded the putative recombinant sequences yielded similar results (Table S1), although the fluctuations at day 522 had a p-value of 0.064.
Under a neutral model, the Ne controls the strength of fluctuations in within-patient subpopulation frequencies. Large Ne results in small fluctuations, while small Ne results in large fluctuations. Eventually, if no new subpopulations arise (as was observed in days 1…522), under a neutral model one subpopulation would eventually take over the entire population, eliminating all subpopulation diversity. Thus, we now ask whether the observed diversity at day 522 is consistent with a neutral model, given realistic values for Ne.
Since no significant frequency changes occurred during days 1…32, we pooled those sequences together (n = 64) to derive more accurate subpopulation frequencies (Table 2). Given these inferred frequencies, we compared expected and observed subpopulation frequencies on day 522. We noted a number of potentially unlikely events under a neutral model. Some of these events indicated small frequency fluctuations, and thus large Ne. One was that the observed frequency of subpopulation s6 on day 522 was 5, exactly the expected value. Also, we still observed subpopulations s3 and s4 on day 522, at frequencies similar to that expected from constant within-patient frequencies. Additionally, we observed 4 subpopulations present on day 522, indicating that not much diversity had been lost. Thus, we asked which values of Ne are large enough to be consistent with these observations. On the other hand, other aspects of the data suggested large frequency fluctuations. In particular, subpopulation s1 was not observed on day 522, while our expectation was 2 observations. Hence, we also asked which values of Ne are small enough to be consistent with this observation. Finally, subpopulation s5 was observed at much higher frequency than expected at day 522. For small Ne, we would have expected subpopulation s5 to go extinct, whereas for very large Ne it would have been unlikely to rise to high frequency. Thus we also asked whether any value of Ne makes the observed count of subpopulation s5 plausible. Nearly identical results were observed when putative recombinants were removed (Table S2). Again, these calculations accounted for the stochastic effects of the sampling size at day 522.
Independently, the genetic diversity (θ) and Ne were estimated using three methods (Table 1). The results were similar and unaffected by the exclusion of putative recombinants. The estimates based on Fluctuate showed that θ ranged from 0.015 to 0.051 substitutions/site during the study period corresponding to = 512 with 2σ range ±162 (pink region in Figure 3).
Figure 3 plots the likelihood versus Ne of the scenarios considered above, under a neutral model. While we investigated a large range of possible Ne's (range 5–50,000), all the aspects of the data we have tested are not significantly unlikely under a neutral scenario within the Ne range consistent with our other analyses (Table 1). In particular, all the likelihoods for these individual aspects of the data are larger than 0.05 for Ne≈800. Thus, we cannot reject a neutral model for these data, even though some uncertainty may remain because of the low likelihood (p<0.1) of observing ≥3 taxa of subpopulation s5 at day 522. The modeling results were robust to whether potential recombinants were included or excluded (Figure S4).
Because the subpopulations were present at different frequencies over time we were interested in the potential impact of such fluctuations on the measured evolutionary rate. Clearly, the apparent substitution rate varied greatly over time (Figure S5). Thus, naively measuring the genetic difference between time points may mislead the estimation of the de novo substitution rate. However, the fluctuations may be another mechanism that HIV-1 uses to adapt and evolve its population structure. Hence, to accurately estimate the substitution rate one must take the phylogeny into account.
We used a Bayesian coalescent method to infer the de novo substitution rate within each subpopulation. To account for the frequency variation of each subpopulation we used the Bayesian skyline demographic model, which allows Ne to vary over time in a non-parametric way. The evolutionary rate was inferred as a hyper-parameter using separate, independent trees to describe each of the subpopulations. A relaxed clock model was used to infer a hyper-parameter with individual distributions for each subpopulation. The mean estimated within-subpopulation substitution rate was 2.33×10−3 substitutions site−1 year−1, with a 95% highest posterior density (HPD) interval of 0.94–3.74×10−3 substitutions site−1 year−1. In agreement with our frequency analysis above (Figure 2), no significant deviation from a constant Ne could be observed as the Bayesian skyline could contain a constant Ne within the 95% HPD.
The genetic divergence (nucleotide substitution rate) was further analyzed in subpopulation s6, which constituted the largest group of sequences (Figure 4). The mean pairwise distance (MPD) between clones sampled at the same time, i.e. the population diversity (0 days), was then compared to the divergence at later sampling times. Interestingly, we found that the MPD of sequences separated by about one month's interval (and longer) differed significantly from the intra-sample diversity (p<0.01, Wilcoxon rank sum test), but no significant divergence was seen in shorter time intervals. Hence, this HIV-1 subpopulation had moved significantly in sequence space after about a month.
HIV-1 evolves by introducing mutations (substitutions, indels, recombination) through a “sloppy” replication mechanism, mainly due to the unfaithful replication by the viral reverse transcriptase. These mutations are often deleterious  or otherwise detrimental to virus fitness , . However, some mutants have an advantage as they may allow escape from immune surveillance ,  or more effective infection of certain tissue compartments or cell types, such as cells in the brain or the genital tract , ,  or naïve CD4+ T-cells, which express CXCR4 , . Here we show that in addition to the mutational processes, HIV-1 can alter its population structure by frequency shifts among subpopulations. Because we analyzed a relatively small number of sequences per time point, we were careful to include the sampling into our analysis method. Over short time (days, weeks, months) these fluctuations were consistent with a constant population size, and most mutations that occurred at this time scale were neutral or only weakly selected. On longer time scales we noticed that the fluctuations became significant movements.
Here, we focused on short-term evolutionary processes (days, weeks, and months), whereas earlier studies as well as the generation of, and escape from, neutralizing antibody responses involve time frames of months to years , , , . Clinically, our patient was classified as a slow disease progressor. Genetically, the virus population in our patient was described by co-existing subpopulations. Thus, it is interesting to compare the HIV population genetics of our patient to previously published patients with normal and slow disease progression. In a study by Shankarappa et al, 5 patients had slow disease progression (p2, p3, p7, p9, p11) and 5 had normal progression (p1, p5, p6, p8) . These patients were followed over many years, but interestingly over a sampling period equivalent to ours (522 days, but with fewer samples), patients in both clinical groups showed subpopulation structure qualitatively similar to our patient (Figure S6). Thus, the short-term evolution we study here is likely representative for many patients regardless of disease progression rate.
One might have expected that the persisting subpopulations found in this patient were controlled by balancing selection , , . Directional selection would have favored the fittest of the subpopulations and it would have been unexpected to see them coexist for so long, let alone to have several well separated subpopulations, which implies that they have existed for longer than the study period. Hence, some type of frequency-dependent selection, where the fitness of a variant/subpopulation is dependent on its relative frequency, would be the alternative hypothesis to neutral drift. Here we show that although the immune system partly controls virus replication during the chronic phase of the disease, particularly well in a slow progressor, and where one would expect escape mutants to dominate in env, the genetic evolution is consistent with a neutral process, at least over the time period studied here. In agreement with this, it was recently shown that genetic drift was a main contributor to HIV evolution in culture . Similarly, in several other virus systems with large population sizes and high mutation rates, where deterministic processes are expected, genetic drift was shown to have a larger effect than expected , , , . Furthermore, stochastic evolution during drug treatment of HIV-1 has previously been demonstrated . Thus, also in vivo evolution of HIV-1 during the chronic phase may be largely described by neutral and stochastic processes. We speculate that this might be due to that the immune system “hits” all subpopulations with near equal efficiency.
Our tests for neutrality of the subpopulation frequency fluctuations are of necessity informal. A more formal procedure would assess the likelihood of the data under neutral models with varying Ne and compare with models that additionally include either balancing or directional selection. There exist methods for estimating Ne from multi-allele temporal data (e.g. ), as well as methods for inferring directional selection from two-allele temporal data (e.g. ). However, we are not aware of likelihood methods that include balancing selection and multiple alleles, and their development is beyond the scope of the present study. Hence, we have relied on a less formal method that may not have optimal power, but nevertheless is informative. In addition, our Ne estimate from sequence data were in the order of previously estimated Ne of HIV-1 in chronic infection , , , however, subpopulation structure or non-neutral evolution may bias these estimates, therefore we included a large range of plausible Ne in our test of neutrality (Ne = 5–50,000).
We have sampled free HIV-1 viral particles in plasma but we do not know where these virions were produced. The degree of compartmentalization of HIV-1 replication is uncertain; some researchers have found evidence of compartmentalization whereas others have not , , , , , , . However, in untreated patients most virus in plasma is produced by short-lived activated CD4+ T-lymphocytes ,  and there is no or limited compartmentalization between virus in plasma and lymphocytes , , . Thus, the plasma virus population should be competing for the same resources, which would justify our analysis of whether balancing selection exists. However, we cannot exclude that the frequency fluctuations we see may be due to differential production from different compartments. The subpopulations were present in actively replicating virus since 1) the subpopulations were detected over time at high frequencies (i.e., detected in 7–11 single molecules per time point), 2) the molecules sequenced must represent virions which were replication competent at least in the previous generation, and 3) subpopulation s6 evolved at a measurable evolutionary rate.
It was interesting to note that PNGS were significantly over-represented among positively selected sites. Glycosylation and movement of glycans have been suggested to be an important immune escape mechanism of HIV-1 , , , . Our data are compatible with such a scenario, which suggests that immune escape and positive selection on PNGS may have contributed to the evolution of the genetic subpopulations in our patient. Previously, it was demonstrated quantitatively that a wide range in strength of the autologous neutralizing antibody response between patients and corresponding differences in the impact on the viral population . The fact that we observed positive selection on PNGS, and earlier studies have shown consecutive replacement HIV-1 env sequences ,  and continuous neutralization escape , , do not contradict our observation that evolution in our patient overall was consistent with a neutral model of evolution.
The estimation of the evolutionary rate may be misled if phylogenetic relationships are not accounted for. This may occur in a naïve analysis where genetic distances between two (or more) time points are directly compared. Without accounting for the phylogeny, especially when the population is divided into clear subpopulations, frequency shifts between existing variants will masquerade as de novo mutations. We show that in env, which has the highest evolutionary rate in the HIV-1 genome, the number of de novo mutations accumulated a significant distance after about one month within a subpopulation. Hence, this suggests that sampling more frequently than this may not be useful to estimate the evolutionary rate in patients during the chronic phase. This is in good agreement with previous estimates of significant temporal changes at 22 months based on genomic sequences of similar length (∼1100 nt) in gag-pol, which evolves much slower than env. Note, however, that selection during drug treatment may potentially act upon existing variants/subpopulations at a much faster rate , , but, as we show here, during chronic, asymptomatic, untreated viral infection evolution proceeds mostly by neutral drift over shorter time frames.
In conclusion, we have performed high-frequency sampling of HIV-1 evolution in a chronically infected, untreated patient with slowly progressing disease. We shown that multiple well-separated subpopulations may persist for years, and over weeks and possibly months their frequencies remained constant. Over the time period of years, however, their fluctuations became significant, but were still consistent with a neutral model of evolution. However, sequence-based methods showed that individual sites had experienced positive selection, possibly as the subpopulations were being formed over several years. While the subpopulation frequencies fluctuated consistent with neutrality, the divergence within a subpopulation showed a temporal trend that was resolved at about one month's time.
The patient was a treatment naive, asymptomatic man that had been HIV-1 infected for approximately 7 years at the start of the study (day 1). The plasma viral load had been stable and relatively low for several years and ranged from 450 to 1220 copies per ml during the main study period. The CD4 count was around 600 at the time point for the first sample, but with previous CD4 fluctuations including values below 500. Thus, the patient did not fulfill the definition of a long-term non-progressor (CD4 counts >500 for more than ten years without antiretroviral), and instead we classify the patient as a slow progressor . In support of this, the patient was put on treatment after more than 15 years of infection.
Blood samples were collected each morning for 12 consecutive days (day 1 through day 12) and then once every week for 3 weeks (day 18, 25 and 32). One later sample, that was collected 1.5 years after the first sample (day 522) was also analyzed as well as an earlier sample that was collected 1.5 years before the start of the study (day −664). Plasma was prepared by centrifugation at 2000 rpm for 10 min at room temperature and stored at −70°C until analysis. Viral RNA was extracted from plasma using the Nuclisense RNA extraction kit (NASBA Diagnostics, Organon Teknika, Boxtel, The Netherlands) according to the manufacturer's instructions, and cDNA synthesis was performed using the First-Strand cDNA synthesis kit (Amersham Pharmacia Biotech, Uppsala, Sweden).
The patient gave written consent and the study was approved by the regional ethics committee (Karolinska Sjukhuset, Lokal forskningsetikkommitte Nord) in Stockholm (Dnr: 98–336).
Single viral molecules were obtained by limiting dilution of the cDNA . The method was selected to minimize the influence of PCR errors in the sequences and to allow sequencing of the entire env gene. According to the Poisson distribution, the likelihood that a positive PCR reaction originates from a single molecule is 0.95 if the fraction of positive reactions is 1∶3. After a dilution series, we determined the template load for each PCR and diluted our template accordingly. Hence, positive PCR samples from dilutions containing less than 1∶3 of positive reactions were sequenced and analyzed. The single molecule status was confirmed by screening for mixed nucleotide positions in the final sequence chromatograms and sequences with mixed positions were excluded. Hence, this procedure will identify PCR errors after the cDNA synthesis as they would be seen in the chromatograms at frequencies ≤25%. In addition, bidirectional sequencing was performed. In one sequence only one mixed position (at 50% in overlapping sequence fragments) was detected (A/G) and in this case both possible sequences were included. A 3.1-kb region covering vpu, env, and one-half of nef was amplified and sequenced as previously described . A nested amplification was used with outer primers JL86 ([gene: 5′-CCGTCTAGATGCTGTTTATTCATTTCAGAATTGG-3′]) and JL89 ([gene: 5′-TCCAGTCCCCCCTTTTCTTTTAAAAA-3′]), and inner primers ED3 ([gene: 5′-TTAGGCATCTCCTATGGCAGGAAGAAGCGG-3′]) and JL88 ([gene: 5′-TAAGTCATTGGTCTTAAAGGTACCTG-3′]). The Expand Long Template kit (Boehringer Mannheim, Indianapolis, IN) was used according to the manufacturer's recommendations and a hot start was achieved by separating the primers and the template from the enzymes (Tgo DNA polymerase and Taq DNA polymerase) with a wax layer (DynaWax; Finnzymes, Espoo, Finland). The PCR program was 94°C for 10 sec, 55°C for 30 sec, and 68°C for 4 min for a total of 30 cycles. Concentrations of 0.4 µM primer and 0.2 mM total dNTP in a final volume of 50 µl were used, and 2 µl of first-round product was transferred to the second-round reaction. Positive reactions were purified using the GFX purification kit (Amersham Biosciences Corp, Piscataway, NJ) and directly sequenced with a walking primer approach using standard dideoxy-terminator fluorescent automated sequencing methodology (Applied Biosystems, Foster City, CA) on ABI 310 or 3100 sequencing machines. Sequencing primers were designed so that each nucleotide of the PCR fragment was detected by at least two separate primers. Hence, all nucleotide calls were made based on at least two sequencing reactions, ensuring high base-call accuracy. The sequences were evaluated and assembled into contigs using Sequencher software (Genecodes Inc, Ann Arbor, USA). The sequences are deposited in Genbank under the accession numbers: JN251812–JN251888.
Sequences were manually aligned to HIV-1 reference sequences using the Se-Al software . The program Modeltest v 3.7  was used to search for the substitution model that best described the evolution of the dataset. ML trees were inferred using PhyML 3.0  using 5 random starting trees with SPR and NNI tree search algorithms. Substitution model parameters were estimated from the data. Topological uncertainty was estimated using maximum likelihood evaluated non-parametric bootstrap analysis with 1000 replicates. Whether the sequence data generally supported a neutral model of evolution was tested using Tajima's D-test, Fu and Li's D*-test and Fay and Wu's H-test , , . Coalescent estimation of Ne was made with the coalescent-likelihood programs Recombine and Fluctuate implemented in the Lamarc 2.1.3 package  as well as by calculating the mean-pairwise distance (MPD) using the program PAUP* .
To exclude possible laboratory contamination and sample mix-up, a phylogenetic tree was constructed where other subtype B env sequences from the HIV sequence database  were included together with the current dataset. This analysis showed that all our sequences formed a monophyletic cluster (not shown).
In order to assess the extent of recombination in our dataset, and possibly identify the recombinants, we applied a procedure that has been shown to be able to identify intra-host recombination . Conflicting phylogenetic signals in the dataset are visualized using the Neighbor Net (NNet) algorithm  implemented in SplitsTree version 4.10  and the presence for recombination signal is then specifically tested with the pairwise homoplasy index (PHI) statistic . The PHI statistic measures the similarity between closely linked sites and the significance of the observed test statistic is obtained using a permutation test. If there is no recombination in the data the genealogical correlation of adjacent sites is invariant to permutation . But in the presence of finite recombination, the order of the sites is important, and distant sites will tend to have less genealogical correlation than adjacent sites. ,  Subpopulations were screened one at a time by the PHI-NNet test. Intra-subpopulation recombinants were removed before screening for putative inter-subpopulation recombinants. As the identification process of putative recombinants may be subjective we wanted to control for human bias in selecting putative recombinants. We therefore randomly removed an equal number of sequences as were determined recombinant and calculated the PHI p-value. This randomized reduction was performed a hundred times.
To verify that the removal of the putative recombinants as determined by the PHI-NNET analysis rendered the dataset free from recombination signal we tested the two alternative datasets with the single breakpoint analysis available at www.datamonkey.org/GARD/.
Recombinant sequences, as determined by the PHI-NNet test, were removed and the alignment stripped so that only single-frame coding regions were present, i.e., only env without vpu/rev. A few spurious stop-codons were conservatively changed to the weighted nucleotide in the corresponding column of the alignment. This will reduce diversity and will not lead to false positive selection detection. We tested whether the identified subpopulations had evolved under different selection pressures by using GAbranch  available at the www.datamonkey.org website . GAbranch automatically partitions all branches in the tree into several selective regimes and performs multi-model inference enabling us to infer dN/dS rates for each branch in the tree without subjectively choosing which branches to test for differential selection. In addition, we tested the subpopulations for site-specific selection or variation using Nielsen & Yang's hierarchical model-pairs (M0, M1a, M2a, M3,) in HyPhy , . Individual amino acid changes were identified over the ML tree within or between subpopulations using MacClade .
Putative subpopulations were identified by high bootstrap values as above. To test whether these subpopulations were statistically significant, we conducted a test for population subdivision originally developed by Hudson et al. and further developed to test HIV-1 intra-patient evolution by Achaz et al. , . The method calculates matrices of pairwise sequence differences for the putative subpopulations as well as for the whole dataset. To assess the significance of the structure the sequences are randomly relabeled into new subsets of populations (keeping n1 and n2 constant), which generates a p-value for the probability that the structure observed was due simply to chance. The test does not rely upon a common genealogy for all sites, which makes it robust to the presence of recombination .
Our significance test includes stochastic effects due to limited sample sizes. To determine whether subpopulation frequencies were significantly fluctuating within the patient, we asked whether the sample frequencies observed on each day J were consistent with the within-patient frequencies inferred from days 1…J-1. On each day J we have NJ total observed sequences, within which subpopulation i has frequency count fiJ. If we assume the within-patient frequencies φi are constant, then given the observations from days 1…J-1 the maximum-likelihood estimate of the within-patient frequency of subpopulation i is simply the fraction of all previously observed sequences that are from subpopulation i. To assess whether the observed frequencies fiJ on a given day J are consistent with the inferred patient frequencies, we use Pearson's χ2 statistic:.
In assessing the significance of the observed χ2 values, we account for our uncertainty in estimating the within-patient frequencies from the day 1…J-1 data using simulations of the estimation process. To do so, we begin with our maximum-likelihood estimate of the within-patient frequencies from the observed data. We then generate simulated data sets using constant within-patient frequencies of , matching the number of samples NJ from each day. For each of these simulations, we calculate χ2 for day J using population frequencies estimated from simulated days 1…J-1. One subtlety is that χ2 will be divergence if a subpopulation is first observed on day J. To avoid this, when calculating χ2 for both the real and simulated data, we introduce one addition count (a pseudocount) for each subpopulation on the first day. Our p-values are then the proportion of simulated data sets that yield a larger χ2 than the real data.
Our neutral simulations begin at day 1 with a population of Ne individuals, divided amongst the 5 subpopulations based on their inferred within-population frequencies from days 1…32. In each generation a new population of Ne individuals is created by sampling with replacement from the prior generation, i.e., a Wright-Fisher model of reproduction. The generation time for HIV-1 is assumed to be 2 days , so that 260 generations separate our samples at day 1 and day 522. For each simulation, based on the final subpopulation frequencies after 260 generations, we calculate the likelihood of observing particular aspects of the real day 522 data. The likelihood of sampling a given number of sequences from population i can be calculated directly from the multinomial probabilities. The likelihood of observing n subpopulations out of T possible is calculated by summing the likelihood of obtaining 0 observations for T-n subpopulations, carefully accounting for the number of ways to do this. Our overall likelihoods represent the average over 104 population simulations. Also, our results are unchanged if we incorporate our uncertainty in estimating by initializing each simulation with different within-patient frequencies consistent with our observations from days 1…32.
The programs TreeRate ,  and BEAST v1.5.1  were used to infer evolutionary rates. TreeRate optimizes the root and the evolutionary rate for a given tree by minimizing tip-height variances at two specified sampling times. The given tree was inferred by PhyML 3.0 as described above, thereby not preconditioned on a molecular clock. In addition we used Bayesian analysis (BEAST) assuming a relaxed molecular clock (uncorrelated lognormal) and a non-parametric population growth model (Bayesian skyline).
Neighbor-Net diagrams showing the evolutionary relationships in the viral population including incompatible signals. Clones from the different time points are indicated with different symbols and colors as shown. The subpopulations are labeled with letters s1–s6 A. All 77 taxa with 8 putative recombinants as determined by the PHI test and indicated with a star. B. The resulting network when these 8 putative recombinants were removed.
Cladogram of the non-recombinant sequences where the dN/dS values inferred through GAbranch are shown. Taxa labels are colored according to genetic subpopulation.
Correlation of the normalized frequency of amino acid (aa) substitutions (A) and normalized potential N-linked glycosylation site (PNGS) replacements (B) to the probability of strong positive selection pressure. The frequencies of aa and PNGS replacements were normalized by the number of tree branches between subpopulations (N = 8) or within subpopulations (N = 130). The selection pressures were partitioned into 3 rate classes (dN/dS = 3.92, dN/dS = 0.55, dN/dS = 0.13), optimized using the Nielsen-Yang model M3 in HyPhy , . The probability of each site of belonging to the dN/dS = 3.92 class was used to measure the strong positive selection pressure. The correlations to strong positive selection were R = 0.78 (between subpopulation aa substitutions), R = 0.69 (within subpopulation aa substitutions), R = 0.40 (between subpopulation PNGS replacements), and R = 0.38 (within subpopulation PNGS replacements). The response to strong positive selection, as measured by OLS regression slopes, was 23 times stronger to between than within subpopulation aa substitutions and 25 times stronger to between than within subpopulation PNGS replacements (p<<0.001, F-statistic, in both cases).
Click here for additional data file (pone.0021747.s003.pdf)
Likelihood under a neutral model of aspects of data with putative recombinants excluded.
The evolutionary rates given in percent substitutions per site and year, measured between all time points with TreeRate. Each arrow begins at the first time-point, and the end of each shaded area represents the second time point. The impact of inclusion or exclusion of putative recombinant is shown; the upper part of each arrow represents exclusion of recombinants (Figure S1), and lower part of arrow represents the results when recombinants were included.
Genetic diversity and divergence over time in our patient and previously published patient data. Our patient (Study patient) was sampled over approximately 3 years, with most samples days, weeks and months apart up to 522 days (∼1.5 years). Comparing our results to an equivalent sampling period of patients in a study by Shankarappa et al , shows that regardless of disease progression rate similar subpopulation structure as in the study patient occurs in at least 5/9 Shankarappa patients (3, 5, 7, 8, 9). Shankarappa patients 2, 3, 7, 9, 11 had slow disease progression, as our patient, and the others normal disease progression. All trees are on the same scale (see scale bar) and the sampling time intervals are also on the same scale, 120 evenly divided colors over 12 years (12 colors shown in legend).
Click here for additional data file (pone.0021747.s006.jpg)
Subpopulation frequency fluctuations. Excluding recombinant sequences from our analysis, we obtain the results shown in Table S1 for the significance of within-patient frequency fluctuations.
Click here for additional data file (pone.0021747.s007.docx)
Subpopulation frequency fluctuations are consistent with neutral drift.Table S2 shows, excluding recombinants, the subpopulation frequencies inferred from our data, along with the expected and observed counts in day 522. Again, we observe several aspects of the data that are informative about potential deviations from neutrality. We test for the likelihood of 1) s1 not being observed, 2) s3 being observed, 3) s4 not being observed, 4) s5 being observed at frequency 2 or greater, 5) s6 being observed at frequency between 3 and 5 (indicating a fluctuation of less than 1 from expected), and 6) observing 3 or more populations. Figure S1 shows the results that none of these aspects of the data are significantly unlikely (p<0.05) under a neutral model.
Click here for additional data file (pone.0021747.s008.docx)
Competing Interests: The authors have declared that no competing interests exist.
Funding: This work was supported by grants from the National Institutes of Health (NIH) [grant 1R01AI087520-01A1], the Swedish Research Council, and the Swedish International Development Cooperation Agency [grant no. SWE-2006-018]. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
The authors would like to thank Sergei Kosakowski Pond for help with the GAbranch analysis.
|1.||Lee HY,Perelson AS,Park SC,Leitner T. Year: 2008Dynamic correlation between intrahost HIV-1 quasispecies evolution and disease progression.PLoS Comput Biol4e100024019079613|
|2.||Shankarappa R,Margolick JB,Gange SJ,Rodrigo AG,Upchurch D,et al. Year: 1999Consistent viral evolutionary changes associated with the progression of human immunodeficiency virus type 1 infection.J Virol73104891050210559367|
|3.||Lemey P,Kosakovsky Pond SL,Drummond AJ,Pybus OG,Shapiro B,et al. Year: 2007Synonymous substitution rates predict HIV disease progression as a result of underlying replication dynamics.PLoS Comput Biol3e2917305421|
|4.||Richman DD,Wrin T,Little SJ,Petropoulos CJ. Year: 2003Rapid evolution of the neutralizing antibody response to HIV type 1 infection.Proc Natl Acad Sci U S A1004144414912644702|
|5.||Lee SK,Xu Z,Lieberman J,Shankar P. Year: 2002The functional CD8 T cell response to HIV becomes type-specific in progressive disease.J Clin Invest1101339134712417573|
|6.||Nielsen R,Yang Z. Year: 1998Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene.Genetics1489299369539414|
|7.||Shriner D,Shankarappa R,Jensen MA,Nickle DC,Mittler JE,et al. Year: 2004Influence of random genetic drift on human immunodeficiency virus type 1 env evolution during chronic infection.Genetics1661155116415082537|
|8.||Brown AJ. Year: 1997Analysis of HIV-1 env gene sequences reveals evidence for a low effective number in the viral population.Proc Natl Acad Sci U S A94186218659050870|
|9.||Perelson AS,Neumann AU,Markowitz M,Leonard JM,Ho DD. Year: 1996HIV-1 dynamics in vivo: virion clearance rate, infected cell life-span, and viral generation time.Science271158215868599114|
|10.||Kouyos RD,Althaus CL,Bonhoeffer S. Year: 2006Stochastic or deterministic: what is the effective population size of HIV-1?Trends Microbiol1450751117049239|
|11.||Nijhuis M,Boucher CA,Schipper P,Leitner T,Schuurman R,et al. Year: 1998Stochastic processes strongly influence HIV-1 evolution during suboptimal protease-inhibitor therapy.Proc Natl Acad Sci U S A9514441144469826719|
|12.||Shriner D,Liu Y,Nickle DC,Mullins JI. Year: 2006Evolution of intrahost HIV-1 genetic diversity during chronic infection.Evolution601165117616892967|
|13.||Achaz G,Palmer S,Kearney M,Maldarelli F,Mellors JW,et al. Year: 2004A robust measure of HIV-1 population turnover within chronically infected individuals.Mol Biol Evol211902191215215321|
|14.||Simmonds P,Balfe P,Peutherer JF,Ludlam CA,Bishop JO,et al. Year: 1990Human immunodeficiency virus-infected individuals contain provirus in small numbers of peripheral mononuclear cells and at low copy numbers.J Virol648648722296085|
|15.||Salazar-Gonzalez JF,Bailes E,Pham KT,Salazar MG,Guffey MB,et al. Year: 2008Deciphering human immunodeficiency virus type 1 transmission and early envelope diversification by single-genome amplification and sequencing.J Virol823952397018256145|
|16.||Gao F,Chen Y,Levy DN,Conway JA,Kepler TB,et al. Year: 2004Unselected mutations in the human immunodeficiency virus type 1 genome are mostly nonsynonymous and often deleterious.J Virol782426243314963138|
|17.||Mansky LM,Temin HM. Year: 1995Lower in vivo mutation rate of human immunodeficiency virus type 1 than that predicted from the fidelity of purified reverse transcriptase.J Virol69508750947541846|
|18.||Neher RA,Leitner T. Year: 2010Recombination rate and selection strength in HIV intra-patient evolution.PLoS Comput Biol6e100066020126527|
|19.||Tajima F. Year: 1989Statistical method for testing the neutral mutation hypothesis by DNA polymorphism.Genetics1235855952513255|
|20.||Fu YX,Li WH. Year: 1993Statistical tests of neutrality of mutations.Genetics1336937098454210|
|21.||Fay JC,Wu CI. Year: 2000Hitchhiking under positive Darwinian selection.Genetics1551405141310880498|
|22.||Edwards CT,Holmes EC,Pybus OG,Wilson DJ,Viscidi RP,et al. Year: 2006Evolution of the human immunodeficiency virus envelope gene is dominated by purifying selection.Genetics1741441145316951087|
|23.||Troyer RM,McNevin J,Liu Y,Zhang SC,Krizan RW,et al. Year: 2009Variable fitness impact of HIV-1 escape mutations to cytotoxic T lymphocyte (CTL) response.PLoS Pathog5e100036519343217|
|24.||Frost SD,Dumaurier MJ,Wain-Hobson S,Brown AJ. Year: 2001Genetic drift and within-host metapopulation dynamics of HIV-1 infection.Proc Natl Acad Sci U S A986975698011381143|
|25.||Allen TM,Altfeld M,Geer SC,Kalife ET,Moore C,et al. Year: 2005Selective escape from CD8+ T-cell responses represents a major driving force of human immunodeficiency virus type 1 (HIV-1) sequence diversity and reveals constraints on HIV-1 evolution.J Virol79132391324916227247|
|26.||Schnell G,Price RW,Swanstrom R,Spudich S. Year: 2010Compartmentalization and clonal amplification of HIV-1 variants in the cerebrospinal fluid during primary infection.J Virol842395240720015984|
|27.||Anderson JA,Ping LH,Dibben O,Jabara CB,Arney L,et al. Year: 2010HIV-1 Populations in Semen Arise through Multiple Mechanisms.PLoS Pathog6|
|28.||Zhang L,Rowe L,He T,Chung C,Yu J,et al. Year: 2002Compartmentalization of surface envelope glycoprotein of human immunodeficiency virus type 1 during acute and chronic infection.J Virol769465947312186928|
|29.||Heeregrave EJ,Geels MJ,Brenchley JM,Baan E,Ambrozak DR,et al. Year: 2009Lack of in vivo compartmentalization among HIV-1 infected naive and memory CD4+ T cell subsets.Virology393243219698967|
|30.||van Rij RP,Blaak H,Visser JA,Brouwer M,Rientsma R,et al. Year: 2000Differential coreceptor expression allows for independent evolution of non-syncytium-inducing and syncytium-inducing HIV-1.J Clin Invest1061039105211032864|
|31.||Wei X,Decker JM,Wang S,Hui H,Kappes JC,et al. Year: 2003Antibody neutralization and escape by HIV-1.Nature42230731212646921|
|32.||Bunnik EM,Pisas L,van Nuenen AC,Schuitemaker H. Year: 2008Autologous neutralizing humoral immunity and evolution of the viral envelope in the course of subtype B human immunodeficiency virus type 1 infection.J Virol827932794118524815|
|33.||Albert J,Abrahamsson B,Nagy K,Aurelius E,Gaines H,et al. Year: 1990Rapid development of isolate-specific neutralizing antibodies after primary HIV-1 infection and consequent emergence of virus variants which resist neutralization by autologous sera.AIDS41071122328092|
|34.||Frost SD,Wrin T,Smith DM,Kosakovsky Pond SL,Liu Y,et al. Year: 2005Neutralizing antibody responses drive the evolution of human immunodeficiency virus type 1 envelope during recent HIV infection.Proc Natl Acad Sci U S A102185141851916339909|
|35.||Trotter MV,Spencer HG. Year: 2008The generation and maintenance of genetic variation by frequency-dependent selection: constructing polymorphisms under the pairwise interaction model.Genetics1801547155718791254|
|36.||Yuste E,Moya A,Lopez-Galindez C. Year: 2002Frequency-dependent selection in human immunodeficiency virus type 1.J Gen Virol8310310611752706|
|37.||Forster R,Wilke CO. Year: 2007Frequency-Dependent Selection in a Periodic Environment.Physica A38125526417940581|
|38.||Voronin Y,Holte S,Overbaugh J,Emerman M. Year: 2009Genetic drift of HIV populations in culture.PLoS Genet5e100043119300501|
|39.||Reuter G,Pankovics P,Szucs G. Year: 2008Genetic drift of norovirus genotype GII-4 in seven consecutive epidemic seasons in Hungary.J Clin Virol4213514018420454|
|40.||Wang WK,Sung TL,Lee CN,Lin TY,King CC. Year: 2002Sequence diversity of the capsid gene and the nonstructural gene NS2B of dengue-3 virus in vivo.Virology30318119112482670|
|41.||Nelson MI,Simonsen L,Viboud C,Miller MA,Taylor J,et al. Year: 2006Stochastic processes are key determinants of short-term evolution in influenza a virus.PLoS Pathog2e12517140286|
|42.||Allain JP,Dong Y,Vandamme AM,Moulton V,Salemi M. Year: 2000Evolutionary rate and genetic drift of hepatitis C virus are not correlated with the host immune response: studies of infected donor-recipient clusters.J Virol742541254910684268|
|43.||Seo TK,Thorne JL,Hasegawa M,Kishino H. Year: 2002Estimation of effective population size of HIV-1 within a host: a pseudomaximum-likelihood approach.Genetics1601283129311973287|
|44.||Edo-Matas D,van Gils MJ,Bowles EJ,Navis M,Rachinger A,et al. Year: 2010Genetic composition of replication competent clonal HIV-1 variants isolated from peripheral blood mononuclear cells (PBMC), HIV-1 proviral DNA from PBMC and HIV-1 RNA in serum in the course of HIV-1 infection.Virology40549250420638697|
|45.||Gantt S,Carlsson J,Heath L,Bull ME,Shetty AK,et al. Year: 2010Genetic analyses of HIV-1 env sequences demonstrate limited compartmentalization in breast milk and suggest viral replication within the breast that increases with mastitis.J Virol84108121081920660189|
|46.||Harrington PR,Schnell G,Letendre SL,Ritola K,Robertson K,et al. Year: 2009Cross-sectional characterization of HIV-1 env compartmentalization in cerebrospinal fluid over the full disease course.AIDS2390791519414991|
|47.||Ince WL,Harrington PR,Schnell GL,Patel-Chhabra M,Burch CL,et al. Year: 2009Major coexisting human immunodeficiency virus type 1 env gene subpopulations in the peripheral blood are produced by cells with similar turnover rates and show little evidence of genetic compartmentalization.J Virol834068408019211740|
|48.||Kemal KS,Foley B,Burger H,Anastos K,Minkoff H,et al. Year: 2003HIV-1 in genital tract and plasma of women: compartmentalization of viral sequences, coreceptor usage, and glycosylation.Proc Natl Acad Sci U S A100129721297714557540|
|49.||Poles MA,Elliott J,Vingerhoets J,Michiels L,Scholliers A,et al. Year: 2001Despite high concordance, distinct mutational and phenotypic drug resistance profiles in human immunodeficiency virus type 1 RNA are observed in gastrointestinal mucosal biopsy specimens and peripheral blood mononuclear cells compared with plasma.J Infect Dis18314314811087199|
|50.||Ho DD,Neumann AU,Perelson AS,Chen W,Leonard JM,et al. Year: 1995Rapid turnover of plasma virions and CD4 lymphocytes in HIV-1 infection.Nature3731231267816094|
|51.||Wei X,Ghosh SK,Taylor ME,Johnson VA,Emini EA,et al. Year: 1995Viral dynamics in human immunodeficiency virus type 1 infection.Nature3731171227529365|
|52.||Heath L,Fox A,McClure J,Diem K,van 't Wout AB,et al. Year: 2009Evidence for limited genetic compartmentalization of HIV-1 between lung and blood.PLoS One4e694919759830|
|53.||Simmonds P,Zhang LQ,McOmish F,Balfe P,Ludlam CA,et al. Year: 1991Discontinuous 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 Virol65626662761920632|
|54.||Sagar M,Wu X,Lee S,Overbaugh J. Year: 2006Human immunodeficiency virus type 1 V1–V2 envelope loop sequences expand and add glycosylation sites over the course of infection, and these modifications affect antibody neutralization sensitivity.J Virol809586959816973562|
|55.||Ganeshan S,Dickover RE,Korber BT,Bryson YJ,Wolinsky SM. Year: 1997Human immunodeficiency virus type 1 genetic evolution in children with different rates of development of disease.J Virol716636778985398|
|56.||Richman DD,Havlir D,Corbeil J,Looney D,Ignacio C,et al. Year: 1994Nevirapine resistance mutations of human immunodeficiency virus type 1 selected during therapy.J Virol68166016667509000|
|57.||Chen R,Yokoyama M,Sato H,Reilly C,Mansky LM. Year: 2005Human immunodeficiency virus mutagenesis during antiviral therapy: impact of drug-resistant reverse transcriptase and nucleoside and nonnucleoside reverse transcriptase inhibitors on human immunodeficiency virus type 1 mutation frequencies.J Virol79120451205716140780|
|58.||Klein MR,Miedema F. Year: 1995Long-term survivors of HIV-1 infection.Trends Microbiol33863918564357|
|59.||Liitsola K,Holm K,Bobkov A,Pokrovsky V,Smolskaya T,et al. Year: 2000An AB recombinant and its parental HIV type 1 strains in the area of the former Soviet Union: low requirements for sequence identity in recombination. UNAIDS Virus Isolation Network.AIDS Res Hum Retroviruses161047105310933619|
|60.||Rambaut A. Year: 1996 (1996) Se-Al: Sequence Alignment Editor, Available from http://tree.bio.ed.ac.uk/software/seal/.|
|61.||Posada D,Crandall KA. Year: 1998MODELTEST: testing the model of DNA substitution.Bioinformatics148178189918953|
|62.||Guindon S,Gascuel O. Year: 2003A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood.Syst Biol5269670414530136|
|63.||Kuhner MK. Year: 2006LAMARC 2.0: maximum likelihood and Bayesian estimation of population parameters.Bioinformatics2276877016410317|
|64.||Swofford DL. Year: 2002PAUP* Phylogenetic Analysis Using Parsimony (*and Other Methods) 4.0b10Sunderland, MAEd. Sinauer Associates|
|65.||Kuiken C,Foley B,Leitner T,Apetrei C,Hahn BH,et al. Year: 2010HIV Sequence Compendium 2010Los Alamos, New MexicoLos Alamos National Laboratory, Theoretical Biology and Biophysics, LA-UR|
|66.||Salemi M,Gray RR,Goodenow MM. Year: 2008An exploratory algorithm to identify intra-host recombinant viral sequences.Mol Phylogenet Evol4961862818801446|
|67.||Bryant D,Moulton V. Year: 2004Neighbor-net: an agglomerative method for the construction of phylogenetic networks.Mol Biol Evol2125526514660700|
|68.||Huson DH,Bryant D. Year: 2006Application of phylogenetic networks in evolutionary studies.Mol Biol Evol2325426716221896|
|69.||Bruen TC,Philippe H,Bryant D. Year: 2006A simple and robust statistical test for detecting the presence of recombination.Genetics1722665268116489234|
|70.||Hudson RR,Kaplan NL. Year: 1985Statistical properties of the number of recombination events in the history of a sample of DNA sequences.Genetics1111471644029609|
|71.||Kosakovsky Pond SL,Posada D,Gravenor MB,Woelk CH,Frost SD. Year: 2006GARD: a genetic algorithm for recombination detection.Bioinformatics223096309817110367|
|72.||Pond SL,Frost SD. Year: 2005A genetic algorithm approach to detecting lineage-specific variation in selection pressure.Mol Biol Evol2247848515509724|
|73.||Pond SL,Frost SD. Year: 2005Datamonkey: rapid detection of selective pressure on individual sites of codon alignments.Bioinformatics212531253315713735|
|74.||Pond SL,Frost SD,Muse SV. Year: 2005HyPhy: hypothesis testing using phylogenies.Bioinformatics2167667915509596|
|75.||Maddison DR,Maddison WP. Year: 2003MacClade 4: Analysis of phylogeny and character evolution. Version 4.06.Sunderland, MassachusettsSinauer Associates|
|76.||Maljkovic Berry I,Athreya G,Kothari M,Daniels M,Bruno WJ,et al. Year: 2009The evolutionary rate dynamically tracks changes in HIV-1 epidemics: Application of a simple method for optimizing the evolutionary rate in phylogenetic trees with longitudinal data.Epidemics In Press, Corrected Proof.|
|77.||Maljkovic Berry I,Ribeiro R,Kothari M,Athreya G,Daniels M,et al. Year: 2007Unequal Evolutionary Rates in the Human Immunodeficiency Virus Type 1 (HIV-1) Pandemic: the Evolutionary Rate of HIV-1 Slows Down When the Epidemic Rate Increases.J Virol81106251063517634235|
|78.||Drummond AJ,Rambaut A. Year: 2007BEAST: Bayesian evolutionary analysis by sampling trees.BMC Evol Biol721417996036|
[Figure ID: pone-0021747-g001]
Figure 1 Maximum likelihood tree of the phylogenetic relationships of the viral subpopulations.
Sequences from the different time points (in days from day 1) are indicated with different symbols and colors as shown. The subpopulations are labeled with letters s1–s6 and the corresponding bootstrap values are shown as ratios of 1000 replicates.
[Figure ID: pone-0021747-g002]
Figure 2 Bar chart showing the observed within-patient frequency fluctuations of the genetic subpopulations during the study period.
Subpopulations as defined in Figure 1 are shown in respective colour and recombinant sequences are marked with diagonal stripes. P-values for tests of constant within-patient subpopulation frequencies are shown above the histogram for each day. Thus for each day J, subpopulation frequencies φi of days 1…J-1 are compared to the φi frequencies of day J. See text for details.
[Figure ID: pone-0021747-g003]
Figure 3 Likelihoods of various aspects of the data under neutral evolution.
The pink shaded region denotes the 2σ range of Ne (512±162) inferred using Fluctuate (Table 1) and the dotted line denotes a cut-off at p = 0.05.
[Figure ID: pone-0021747-g004]
Figure 4 Genetic divergence in subpopulation s6.
Mean pairwise distances were calculated between sequences sampled with different time intervals. At an interval of one month or more, the genetic distances were significantly greater than the intra-sample diversity (0 days interval) (p<0.01, Wilcoxon rank sum test). Sampling intervals of 1–2 days and 3–4 weeks were estimated together and are named 1 day, and 1 month, respectively.
Table 1 Sequence data.
|Day||No. clones||RNA load||MPDa||Fluctuate||Recombine|
aMean Pairwise Distance, as measured by PAUP* using a GTR substitution model.
bGenetic diversity (substitutions/site).
cExponential growth rate.
dEffective population size determined from θ = 2Neμ with μ = 3.4×10−5 substitutions site−1 generation−1.
fNot analyzed because of the small sample size.
gRecombination rate, C/μ, where C is the rate of recombination per inter-site link per generation, and μ is the substitution rate per site per generation.
Table 2 Subpopulation frequencies: inferred, expected, and observed.
|Subpopulation||φi,32±2σ||Expected fi,522 (out of 10)||Observed fi,522|
Previous Document: Identification of Novel Immunoregulatory Molecules in Human Thymic Regulatory CD4CD25 T Cells by Pha...
Next Document: Transforming Growth Factor ? Signaling Pathway Associated Gene Polymorphisms May Explain Lower Breas...