Human immunodeficiency virus type 1 genetic evolution in children with different rates of development of disease.

The rate of development of disease varies considerably among human immunodeficiency virus type 1 (HIV-1)-infected children. The reasons for these observed differences are not clearly understood but most probably depend on the dynamic interplay between the HIV-1 quasispecies virus population and the immune constraints imposed by the host. To study the relationship between disease progression and genetic diversity, we analyzed the evolution of viral sequences within six perinatally infected children by examining proviral sequences spanning the C2 through V5 regions of the viral envelope gene by PCR of blood samples obtained at sequential visits. PCR product DNAs from four sample time points per child were cloned, and 10 to 13 clones from each sample were sequenced. Greater genetic distances relative to the time of infection were found for children with low virion-associated RNA burdens and slow progression to disease relative to those found for children with high virion-associated RNA burdens and rapid progression to disease. The greater branch lengths observed in the phylogenetic reconstructions correlated with a higher accumulation rate of nonsynonymous base substitutions per potential nonsynonymous site, consistent with positive selection for change rather than a difference in replication kinetics. Viral sequences from children with slow progression to disease also showed a tendency to form clusters that associated with different sampling times. These progressive shifts in the viral population were not found in viral sequences from children with rapid progression to disease. Therefore, despite the HIV-1 quasispecies being a diverse, rapidly evolving, and competing population of genetic variants, different rates of genetic evolution could be found under different selective constraints. These data suggest that the evolutionary dynamics exhibited by the HIV-1 quasispecies virus populations are compatible with a Darwinian system evolving under the constraints of natural selection.     

The Genealogical Population Dynamics of HIV-1 in a Large Transmission Chain: Bridging within and among Host Evolutionary Rates

Transmission lies at the interface of human immunodeficiency virus type 1 (HIV-1) evolution within and among hosts and separates distinct selective pressures that impose differences in both the mode of diversification and the tempo of evolution. In the absence of comprehensive direct comparative analyses of the evolutionary processes at different biological scales, our understanding of how fast within-host HIV-1 evolutionary rates translate to lower rates at the between host level remains incomplete. Here, we address this by analyzing pol and env data from a large HIV-1 subtype C transmission chain for which both the timing and the direction is known for most transmission events. To this purpose, we develop a new transmission model in a Bayesian genealogical inference framework and demonstrate how to constrain the viral evolutionary history to be compatible with the transmission history while simultaneously inferring the within-host evolutionary and population dynamics. We show that accommodating a transmission bottleneck affords the best fit our data, but the sparse within-host HIV-1 sampling prevents accurate quantification of the concomitant loss in genetic diversity. We draw inference under the transmission model to estimate HIV-1 evolutionary rates among epidemiologically-related patients and demonstrate that they lie in between fast intra-host rates and lower rates among epidemiologically unrelated individuals infected with HIV subtype C. Using a new molecular clock approach, we quantify and find support for a lower evolutionary rate along branches that accommodate a transmission event or branches that represent the entire backbone of transmitted lineages in our transmission history. Finally, we recover the rate differences at the different biological scales for both synonymous and non-synonymous substitution rates, which is only compatible with the ‘store and retrieve’ hypothesis positing that viruses stored early in latently infected cells preferentially transmit or establish new infections upon reactivation.     

High-avidity, high-IFNγ-producing CD8 T-cell responses following immune selection during HIV-1 infection

HIV-1 mutations, which reduce or abolish CTL responses against virus-infected cells, are frequently selected in acute and chronic HIV infection. Among population HIV-1 sequences, immune selection is evident as human leukocyte antigen (HLA) allele-associated substitutions of amino acids within or near CD8 T-cell epitopes. In these cases, the non-adapted epitope is susceptible to immune recognition until an escape mutation renders the epitope less immunogenic. However, several population-based studies have independently identified HLA-associated viral changes, which lead to the formation of a new T-cell epitope, suggesting that the immune responses that these variants or 'neo-epitopes' elicit provide an evolutionary advantage to the virus rather than the host. Here, we examined the functional characteristics of eight CD8 T-cell responses that result from viral adaptation in 125 HLA-genotyped individuals with chronic HIV-1 infection. Neo-epitopes included well-characterized immunodominant epitopes restricted by common HLA alleles, and in most cases the T-cell responses against the neo-epitope showed significantly greater functional avidity and higher IFNγ production than T cells for non-adapted epitopes, but were not more cytotoxic. Neo-epitope formation and emergence of cognate T-cell response coincident with a rise in viral load was then observed in vivo in an acutely infected individual. These findings show that HIV-1 adaptation not only abrogates the immune recognition of early targeted epitopes, but may also increase immune recognition to other epitopes, which elicit immunodominant but non-protective T-cell responses. These data have implications for immunodominance associated with polyvalent vaccines based on the diversity of chronic HIV-1 sequences.     

Infection with multiple human immunodeficiency virus type 1 variants is associated with faster disease progression

Human immunodeficiency virus type 1 (HIV-1)-infected individuals develop a genetically diverse virus population over time, but often only a limited number of viral variants are transmitted from a chronic carrier to a newly infected person. Interestingly, many women but few men are infected by multiple HIV-1 variants from a single partner. To determine whether the complexity of the infecting virus population influences clinical outcome, we examined viral diversity in the HIV-1 envelope sequences present at primary infection in 156 women from Kenya for whom we had follow-up data on viral RNA levels and CD4 T-cell counts. Eighty-nine women had multiple viral genotypes, while 67 women had a single genotype at primary infection. Women who acquired multiple viral genotypes had a significantly higher viral load (median, 4.84 versus 4.64 log(10) copies/ml, P = 0.04) and a significantly lower CD4(+)-T-cell count (median, 416 versus 617 cells/mm(3), P = 0.01) 4 to 24 months after infection compared to women who were infected with a single viral genotype. These studies suggest that early HIV-1 genetic diversity is linked to faster disease progression.     

Subtypes of human immunodeficiency virus type 1 and disease stage among women in Nairobi, Kenya

In sub-Saharan Africa, where the effects of human immunodeficiency virus type 1 (HIV-1) have been most devastating, there are multiple subtypes of this virus. The distribution of different subtypes within African populations is generally not linked to particular risk behaviors. Thus, Africa is an ideal setting in which to examine the diversity and mixing of viruses from different subtypes on a population basis. In this setting, it is also possible to address whether infection with a particular subtype is associated with differences in disease stage. To address these questions, we analyzed the HIV-1 subtype, plasma viral loads, and CD4 lymphocyte levels in 320 women from Nairobi, Kenya. Subtype was determined by a combination of heteroduplex mobility assays and sequence analyses of envelope genes, using geographically diverse subtype reference sequences as well as envelope sequences of known subtype from Kenya. The distribution of subtypes in this population was as follows: subtype A, 225 (70.3%); subtype D, 65 (20.5%); subtype C, 22 (6.9%); and subtype G, 1 (0.3%). Intersubtype recombinant envelope genes were detected in 2.2% of the sequences analyzed. Given that the sequences analyzed represented only a small fraction of the proviral genome, this suggests that intersubtype recombinant viral genomes may be very common in Kenya and in other parts of Africa where there are multiple subtypes. The plasma viral RNA levels were highest in women infected with subtype C virus, and women infected with subtype C virus had significantly lower CD4 lymphocyte levels than women infected with the other subtypes. Together, these data suggest that women in Kenya who are infected with subtype C viruses are at more advanced stages of immunosuppression than women infected with subtype A or D. There are at least two models to explain the data from this cross-sectional study; one is that infection with subtype C is associated with a more rapid disease progression, and the second is that subtype C represents an older epidemic in Kenya. Discriminating between these possibilities in a longitudinal study will be important for increasing our understanding of the role of specific subtypes in the transmission and pathogenesis of HIV-1.     

Unequal 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

HIV-1 sequences in intravenous drug user (IDU) networks are highly homogenous even after several years, while this is not observed in most sexual epidemics. To address this disparity, we examined the human immunodeficiency virus type 1 (HIV-1) evolutionary rate on the population level for IDU and heterosexual transmissions. All available HIV-1 env V3 sequences from IDU outbreaks and heterosexual epidemics with known sampling dates were collected from the Los Alamos HIV sequence database. Evolutionary rates were calculated using phylogenetic trees with a t test root optimization of dated samples. The evolutionary rate of HIV-1 subtype A1 was found to be 8.4 times lower in fast spread among IDUs in the former Soviet Union (FSU) than in slow spread among heterosexual individuals in Africa. Mixed epidemics (IDU and heterosexual) showed intermediate evolutionary rates, indicating a combination of fast- and slow-spread patterns. Hence, if transmissions occur repeatedly during the initial stage of host infection, before selective pressures of the immune system have much impact, the rate of HIV-1 evolution on the population level will decrease. Conversely, in slow spread, where HIV-1 evolves under the pressure of the immune system before a donor infects a recipient, the virus evolution at the population level will increase. Epidemiological modeling confirmed that the evolutionary rate of HIV-1 depends on the rate of spread and predicted that the HIV-1 evolutionary rate in a fast-spreading epidemic, e.g., for IDUs in the FSU, will increase as the population becomes saturated with infections and the virus starts to spread to other risk groups.     

Using human immunodeficiency virus type 1 sequences to infer historical features of the acquired immune deficiency syndrome epidemic and human immunodeficiency virus evolution

In earlier work, human immunodeficiency virus type 1 (HIV-1) sequences were analysed to estimate the timing of the ancestral sequence of the main group of HIV-1, the virus that is responsible for the acquired immune deficiency syndrome pandemic, yielding a best estimate of 1931 (95% confidence interval of 1915-1941). That work will be briefly reviewed, outlining how phylogenetic tools were extended to incorporate improved evolutionary models, how the molecular clock model was adapted to incorporate variable periods of latency, and how the approach was validated by correctly estimating the timing of two historically documented dates. The advantages, limitations, and assumptions of the approach will be summarized, with particular consideration of the implications of branch length uncertainty and recombination. We have recently undertaken new phylogenetic analysis of an extremely diverse set of human immunodeficiency virus envelope sequences from the Democratic Republic of the Congo (the DRC, formerly Zaire). This analysis both corroborates and extends the conclusions of our original study. Coalescent methods were used to infer the demographic history of the HIV-1 epidemic in the DRC, and the results suggest an increase in the exponential growth rate of the infected population through time.     

Evolution of human immunodeficiency virus type 1 coreceptor usage during antiretroviral Therapy: a Bayesian approach

There is substantial evidence for ongoing replication and evolution of human immunodeficiency virus type 1 (HIV-1), even in individuals receiving highly active antiretroviral therapy. Viral evolution in the presence of antiviral therapy needs to be considered when developing new therapeutic strategies. Phylogenetic analyses of HIV-1 sequences can be used for this purpose but may give rise to misleading results if rates of intrapatient evolution differ significantly. To improve analyses of HIV-1 evolution relevant to studies of pathogenesis and treatment, we developed a Bayesian hierarchical model that incorporates all available sequence data while simultaneously allowing the phylogenetic parameters of each patient to vary. We used this method to examine evolutionary changes in HIV-1 coreceptor usage in response to treatment. We examined patients whose viral populations exhibited a shift in coreceptor utilization in response to therapy. CXCR4 (X4) strains emerged in each patient but were suppressed following initiation of new antiretroviral regimens, so that CCR5-utilizing (R5) strains predominated. By phylogenetically reconstructing the evolutionary relationship of HIV-1 obtained longitudinally from each patient, it was possible to examine the origin of the reemergent R5 virus. Using our Bayesian hierarchical approach, we found that the reemergent R5 virus detectable after therapy was more closely related to the predecessor R5 virus than to the X4 strains. The Bayesian hierarchical approach, unlike more traditional methods, makes it possible to evaluate competing hypotheses across patients. This model is not limited to analyses of HIV-1 but can be used to elucidate evolutionary processes for other organisms as well.     

Exceptional Heterogeneity in Viral Evolutionary Dynamics Characterises Chronic Hepatitis C Virus Infection

The treatment of HCV infection has seen significant progress, particularly since the approval of new direct-acting antiviral drugs. However these clinical achievements have been made despite an incomplete understanding of HCV replication and within-host evolution, especially compared with HIV-1. Here, we undertake a comprehensive analysis of HCV within-host evolution during chronic infection by investigating over 4000 viral sequences sampled longitudinally from 15 HCV-infected patients. We compare our HCV results to those from a well-studied HIV-1 cohort, revealing key differences in the evolutionary behaviour of these two chronic-infecting pathogens. Notably, we find an exceptional level of heterogeneity in the molecular evolution of HCV, both within and among infected individuals. Furthermore, these patterns are associated with the long-term maintenance of viral lineages within patients, which fluctuate in relative frequency in peripheral blood. Together, our findings demonstrate that HCV replication behavior is complex and likely comprises multiple viral subpopulations with distinct evolutionary dynamics. The presence of a structured viral population can explain apparent paradoxes in chronic HCV infection, such as rapid fluctuations in viral diversity and the reappearance of viral strains years after their initial detection.     

Predicting demographic group structures based on DNA sequence data

The ability to infer relationships between groups of sequences, either by searching for their evolutionary history or by comparing their sequence similarity, can be a crucial step in hypothesis testing. Interpreting relationships of human immunodeficiency virus type 1 (HIV-1) sequences can be challenging because of their rapidly evolving genomes, but it may also lead to a better understanding of the underlying biology. Several studies have focused on the evolution of HIV-1, but there is little information to link sequence similarities and evolutionary histories of HIV-1 to the epidemiological information of the infected individual. Our goal was to correlate patterns of HIV-1 genetic diversity with epidemiological information, including risk and demographic factors. These correlations were then used to predict epidemiological information through analyzing short stretches of HIV-1 sequence. Using standard phylogenetic and phenetic techniques on 100 HIV-1 subtype B sequences, we were able to show some correlation between the viral sequences and the geographic area of infection and the risk of men who engage in sex with men. To help identify more subtle relationships between the viral sequences, the method of multidimensional scaling (MDS) was performed. That method identified statistically significant correlations between the viral sequences and the risk factors of men who engage in sex with men and individuals who engage in sex with injection drug users or use injection drugs themselves. Using tree construction, MDS, and newly developed likelihood assignment methods on the original 100 samples we sequenced, and also on a set of blinded samples, we were able to predict demographic/risk group membership at a rate statistically better than by chance alone. Such methods may make it possible to identify viral variants belonging to specific demographic groups by examining only a small portion of the HIV-1 genome. Such predictions of demographic epidemiology based on sequence information may become valuable in assigning different treatment regimens to infected individuals.     

Coestimation of recombination,substitution and molecular adaptation rates by approximate Bayesian computation

The estimation of parameters in molecular evolution may be biased when some processes are not considered. For example, the estimation of selection at the molecular level using codon-substitution models can have an upward bias when recombination is ignored. Here we address the joint estimation of recombination, molecular adaptation and substitution rates from coding sequences using approximate Bayesian computation (ABC). We describe the implementation of a regression-based strategy for choosing subsets of summary statistics for coding data, and show that this approach can accurately infer recombination allowing for intracodon recombination breakpoints, molecular adaptation and codon substitution rates. We demonstrate that our ABC approach can outperform other analytical methods under a variety of evolutionary scenarios. We also show that although the choice of the codon-substitution model is important, our inferences are robust to a moderate degree of model misspecification. In addition, we demonstrate that our approach can accurately choose the evolutionary model that best fits the data, providing an alternative for when the use of full-likelihood methods is impracticable. Finally, we applied our ABC method to co-estimate recombination, substitution and molecular adaptation rates from 24 published human immunodeficiency virus 1 coding data sets.     

Long-term experimental evolution of HIV-1 reveals effects of environment and mutational history

An often-returning question for not only HIV-1, but also other organisms, is how predictable evolutionary paths are. The environment, mutational history, and random processes can all impact the exact evolutionary paths, but to which extent these factors contribute to the evolutionary dynamics of a particular system is an open question. Especially in a virus like HIV-1, with a large mutation rate and large population sizes, evolution is expected to be highly predictable if the impact of environment and history is low, and evolution is not neutral. We investigated the effect of environment and mutational history by analyzing sequences from a long-term evolution experiment, in which HIV-1 was passaged on 2 different cell types in 8 independent evolutionary lines and 8 derived lines, 4 of which involved a switch of the environment. The experiments lasted for 240–300 passages, corresponding to approximately 400–600 generations or almost 3 years. The sequences show signs of extensive parallel evolution—the majority of mutations that are shared between independent lines appear in both cell types, but we also find that both environment and mutational history significantly impact the evolutionary paths. We conclude that HIV-1 evolution is robust to small changes in the environment, similar to a transmission event in the absence of an immune response or drug pressure. We also find that the fitness landscape of HIV-1 is largely smooth, although we find some evidence for both positive and negative epistatic interactions between mutations.

Analysis of the longest evolutionary experiment with HIV-1 to-date reveals continuous viral adaptation over several years. The authors quantify the environment-specific mutations that arise and determine the fraction of mutations that co-occur with significantly different frequencies than expected by chance.     

Intra-Host Diversity and Emergence of Unique GBV-C Viral Lineages in HIV Infected Subjects in Central China

GB virus C (GBV-C), which is highly prevalent among HIV/AIDS, seemed to slow the HIV disease progression. The HIV/GBV-C co-infected individuals may represent an interesting model for the investigation of the role played by HIV infection and/or the immune system in driving the evolution of the GBV-C viral populations. The present study investigated the prevalence and population dynamics of GB virus C in HIV infected individuals representing 13 geographic regions of Hubei Province of China. Approximately 37% of HIV-1 infected individuals were infected with GBV-C and genotype 3 is appeared to be predominant. Utilizing the 196 complete E2 nucleotide sequence data from 10 HIV/GBV-C infected individuals and employing coalescence based phylogenetic approaches; the present study has investigated the intra-host dynamics of GBV-C. The results revealed patient-specific unique GBV-C viral lineages and each viral lineage showed the evidence of rapid population expansion in respective HIV-1 infected patients, thus suggesting HIV-1 was unlikely to have been inhibiting effect on the GBV-C viral replication. GBV-C in all patients has experienced intense purifying selection, suggesting the GBV-C viral invasion and subsequent expansion within the HIV-1 infected hosts without any modification of the functional epitopes at their membrane protein. The finding of within host GBV-C recombinant sequences indicated recombination was one of the significant forces in the evolution and divergence of GBV-C.     

Episomal viral cDNAs identify a reservoir that fuels viral rebound after treatment interruption and that contributes to treatment failure

Viral reservoirs that persist in HIV-1 infected individuals on antiretroviral therapy (ART) are the major obstacle to viral eradication. The identification and definition of viral reservoirs in patients on ART is needed in order to understand viral persistence and achieve the goal of viral eradication. We examined whether analysis of episomal HIV-1 genomes provided the means to characterize virus that persists during ART and whether it could reveal the virus that contributes to treatment failure in patients on ART. For six individuals in which virus replication was highly suppressed for at least 20 months, proviral and episomal genomes present just prior to rebound were phylogenetically compared to RNA genomes of rebounding virus after therapy interruption. Episomal envelope sequences, but not proviral envelope sequences, were highly similar to sequences in rebounding virus. Since episomes are products of recent infections, the phylogenetic relationships support the conclusion that viral rebound originated from a cryptic viral reservoir. To evaluate whether the reservoir revealed by episomal sequence analysis was of clinical relevance, we examined whether episomal sequences define a viral population that contributes to virologic failure in individuals receiving the CCR5 antagonist, Vicriviroc. Episomal envelope sequences at or near baseline predicted treatment failure due to the presence of X4 or D/M (dual/mixed) viral variants. In patients that did not harbor X4 or D/M viruses, the basis for Vicriviroc treatment failure was indeterminate. Although these samples were obtained from viremic patients, the assay would be applicable to a large percentage of aviremic patients, based on previous studies. Summarily, the results support the use of episomal HIV-1 as an additional or alternative approach to traditional assays to characterize virus that is maintained during long-term, suppressive ART.     

Selection for human immunodeficiency virus type 1 envelope glycosylation variants with shorter V1-V2 loop sequences occurs during transmission of certain genetic subtypes and may impact viral RNA levels

Designing an effective human immunodeficiency virus type 1 (HIV-1) vaccine will rely on understanding which variants, from among the myriad of circulating HIV-1 strains, are most commonly transmitted and determining whether such variants have an Achilles heel. Here we show that heterosexually acquired subtype A HIV-1 envelopes have signature sequences that include shorter V1-V2 loop sequences and fewer predicted N-linked glycosylation sites relative to the overall population of circulating variants. In contrast, recently transmitted subtype B variants did not, and this was true for cases where the major risk factor was homosexual contact, as well as for cases where it was heterosexual contact. This suggests that selection during HIV-1 transmission may vary depending on the infecting subtype. There was evidence from 23 subtype A-infected women for whom there was longitudinal data that those who were infected with viruses with fewer potential N-linked glycosylation sites in V1-V2 had lower viral set point levels. Thus, our study also suggests that the extent of glycosylation in the infecting virus could impact disease progression.     

Memory in viral quasispecies

Biological adaptive systems share some common features: variation among their constituent elements and continuity of core information. Some of them, such as the immune system, are endowed with memory of past events. In this study we provide direct evidence that evolving viral quasispecies possess a molecular memory in the form of minority components that populate their mutant spectra. The experiments have involved foot-and-mouth disease virus populations with known evolutionary histories. The composition and behavior of the viral population in response to a selective constraint were influenced by past evolutionary history in a way that could not be predicted from examination of consensus nucleotide sequences of the viral populations. The molecular memory of the viral quasispecies influenced both the nature and the intensity of the response of the virus to a selective constraint.