Slower evolution of human immunodeficiency virus type 1 quasispecies during progression to AIDS.

The evolution of human immunodeficiency virus type 1 (HIV-1) quasispecies at the envelope gene was studied from the time of infection in 11 men who experienced different rates of CD4+ cell count decline and 6 men with unknown dates of infection by using DNA heteroduplex mobility assays. Quasispecies were genetically homogeneous near the time of seroconversion. Subsequently, slower proviral genetic diversification and higher plasma viremia correlated with rapid CD4+ cell count decline. Except for the fastest progressors to AIDS, highly diverse quasispecies developed in all subjects within 3 to 4 years. High quasispecies diversity was then maintained for years until again becoming more homogeneous in a subset of late-stage AIDS patients. Individuals who maintained high CD4+ cell counts showed continuous genetic turnover of their complex proviral quasispecies, while more closely related sets of variants were found in longitudinal samples of severely immunocompromised patients. The limited number of variants that grew out in short-term PBMC cocultures were rare in the uncultured proviral quasispecies of healthy, long-term infected individuals but more common in vivo in patients with low CD4+ cell counts. The slower evolution of HIV-1 observed during rapid progression to AIDS and in advanced patients may reflect ineffective host-mediated selection pressures on replicating quasispecies.     

Fitness landscape of human immunodeficiency virus type 1 protease quasispecies

Here we show, at a high resolution (1%), the human immunodeficiency virus type 1 (HIV-1) protease gene quasispecies landscape from three infected na?ve individuals. A huge range of genetic configurations was found (67%, 71%, and 80% of the nucleotide clones from the three individuals, respectively, were different), and these configurations created a dense net that linked different parts of the viral population. Similarly, a vast diversity of different protease activities was also found. Importantly, 65% of the analyzed enzymes had detectable protease activity, and 11% of the minority individual variants showed similar or better fitness than the master (most abundant) enzyme, suggesting that the viral complexity in this genomic region does not exclusively depend on the enzyme's catalytic efficiency. Several high-fitness minority variants had only one substitution compared to the master sequence, supporting the possibility that the rugged HIV-1 protease quasispecies fitness landscape may be formed by a continuous network that can be traversed by single mutational steps without passing through defective or less-adapted proteins.     

Perinatal transmission of major, minor, and multiple maternal human immunodeficiency virus type 1 variants in utero and intrapartum

Previous studies have provided conflicting data on the presence of selective pressures in the transmission of a homogeneous maternal viral subpopulation to the infant. Therefore, the purpose of this study was to definitively characterize the human immunodeficiency virus type 1 (HIV-1) quasispecies transmitted in utero and intrapartum. HIV-1 envelope gene diversity from peripheral blood mononuclear cells and plasma was measured during gestation and at delivery in mothers who did and did not transmit HIV perinatally by using a DNA heteroduplex mobility assay. Children were defined as infected in utero or intrapartum based on the timing of the first detection of HIV. Untreated transmitting mothers (n = 19) had significantly lower HIV-1 quasispecies diversity at delivery than untreated nontransmittting mothers (n = 18) (median Shannon entropy, 0.711 [0.642 to 0.816] versus 0.853 [0.762 to 0.925], P = 0.005). Eight mothers transmitted a single major env variant to their infants in utero, and one mother transmitted a single major env variant intrapartum. Four mothers transmitted multiple HIV-1 env variants to their infants in utero, and two mothers transmitted multiple env variants intrapartum. The remaining six intrapartum- and two in utero-infected infants had a homogeneous HIV-1 env quasispecies which did not comigrate with their mothers' bands at their first positive time point. In conclusion, in utero transmitters were more likely to transmit single or multiple major maternal viral variants. In contrast, intrapartum transmitters were more likely to transmit minor HIV-1 variants. These data indicate that different selective pressures, depending on the timing of transmission, may be involved in determining the pattern of maternal HIV-1 variant transmission.     

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.     

Tracking hepatitis C virus quasispecies major and minor variants in symptomatic and asymptomatic liver transplant recipients.

To evaluate the possibility that distinct viral quasispecies play a role in the pathogenesis of progressive hepatitis C virus (HCV) infection, we performed a detailed evaluation of HCV quasispecies before and after liver transplantation in five patients infected with HCV genotype 1, three of whom developed severe recurrent hepatitis C and two of whom developed asymptomatic posttransplant infections with high-titered viremia. HCV quasispecies were characterized by using a combination of nucleotide sequencing plus heteroduplex tracking assay of the second envelope gene hypervariable region (HVR). An average of 30 HVR clones were analyzed per specimen; an average of five specimens were analyzed per patient over a 6- to 24-month study period. The complexity of HCV quasispecies in pretransplant serum varied, ranging from one to nine genetically distinct variants for the five patients. However, in all five cases, relatively homogenous quasispecies variants emerged after liver transplantation. In the three patients who developed recurrent hepatitis, quasispecies major variants present in pretransplant serum were efficiently propagated immediately after liver transplantation and were propagated throughout the course of acute and chronic hepatitis. In contrast, in the two asymptomatic cases, we observed rapid depletion of pretransplant quasispecies major variants from posttransplant serum, followed by emergence of new quasispecies variants by posttransplant day 30. Genetic analysis suggested that in these cases, the new quasispecies variants were derived from minor variants present at relatively low clonal frequency (less than 5% of HVR clones) within the pretransplant quasispecies populations. These data demonstrate that quasispecies tracking patterns are associated with the rapidity and severity of HCV-associated liver disease after liver transplantation. Further characterization of HCV quasispecies in animal model systems is warranted.     

Sequencing-based detection of low-frequency human immunodeficiency virus type 1 drug-resistant mutants by an RNA/DNA heteroduplex generator-tracking assay

Drug-resistant viruses may be present as minority variants during early treatment failures or following discontinuation of failed antiretroviral regimens. A limitation of the traditional direct PCR population sequencing method is its inability to detect human immunodeficiency virus type 1 (HIV-1) variants present at frequencies lower than 20%. A drug resistance genotyping assay based on the isolation and DNA sequencing of minority HIV protease variants is presented here. A multiple-codon-specific heteroduplex generator probe was constructed to improve the separation of HIV protease genes varying in sequence at 12 codons associated with resistance to protease inhibitors. Using an RNA molecule as probe allowed the simple sequencing of protease variants isolated as RNA/DNA heteroduplexes with different electrophoretic mobilities. The protease gene RNA heteroduplex generator-tracking assay (RNA-HTA) was tested on plasma quasispecies from 21 HIV-1-infected persons in whom one or more protease resistance mutations emerged during therapy or following initiation of salvage regimens. In 11 of 21 cases, RNA-HTA testing of virus from the first episode of virologic failure identified protease resistance mutations not seen by population-based PCR sequencing. In 8 of these 11 cases, all of the low-frequency drug resistance mutations detected exclusively by RNA-HTA during the first episode became detectable by population-based PCR sequencing at the later time point. Distinct sets of protease mutations could be linked on different genomes in patients with high-frequency protease gene lineages. The enhanced detection of minority drug resistance variants using a sequencing-based assay may improve the efficacy of genotype-assisted salvage therapies.     

In-depth characterization of viral isolates from plasma and cells compared with plasma circulating quasispecies in early HIV-1 infection

Background

The use of in vitro models to unravel the phenotypic characteristics of circulating viral variants is key to understanding HIV-1 pathogenesis but limited by the availability of primary viral isolates from biological samples. However, overall in vivo genetic variability of HIV-1 within a subject may not be reflected in the viable viral population obtained after isolation. Although several studies have tried to determine whether viral populations expanded in vitro are representative of in vivo findings, the answer remains unclear due to the reduced number of clonal sequences analyzed or samples compared. In order to overcome previous experimental limitations, here we applied Deep Pyrosequencing (DPS) technology in combination with phenotypic experiments to analyze and compare with unprecedented detail the composition of viral isolates and in vivo quasispecies.

Methodology/Principal Findings

We amplified by DPS HIV-1 genomic regions covering gag, protease, integrase and env-V3 to characterize paired isolates from plasma and peripheral blood mononuclear cells and compare them with total plasma viral RNA in four recently HIV-1 infected subjects. Our study demonstrated the presence of unique haplotypes scattered between sample types with conservation of major variants. In addition, no differences in intra- and inter-population encoded protein variability were found between the different types of isolates or when these were compared to plasma viral RNA within subjects. Additionally, in vitro experiments demonstrated phenotypic similarities in terms of replicative capacity and co-receptor usage between viral isolates and plasma viral RNA.

Conclusion

This study is the first in-depth comparison and characterization of viral isolates from different sources and plasma circulating quasispecies using DPS in recently HIV-1 infected subjects. Our data supports the use of primary isolates regardless of their plasma or cellular origin to define genetic variability and biological traits of circulating HIV-1 quasispecies.     

Human immunodeficiency virus type 1 evolution in vivo tracked by DNA heteroduplex mobility assays.

High mutation rates and strong selective pressures imposed on human immunodeficiency viruses in vivo result in the formation of pools of genetic variants known as quasispecies. DNA heteroduplex mobility and tracking analyses were used to monitor the generation of HIV sequence diversity, to estimate quasispecies complexity, and to assess the turnover of genetic variants to approach an understanding of the relationship between viral quasispecies evolution in vivo and disease progression. Proviral DNA pools were nearly homogeneous soon after sexual transmission. The emergence and clearance of individual variants then occurred at different rates in different individuals. High quasispecies complexity was found in long-term-infected, asymptomatic individuals, while rapid CD4+ cell decline and AIDS were often, but not always, associated with lower quasispecies complexity. Proviral genetic variation was often low following in vitro culture, because of the outgrowth of one or a few variants that often became more abundant only later as proviruses in peripheral blood mononuclear cells. These studies provide insight into the dynamics of human immunodeficiency virus sequence changes in vivo and illustrate the utility of heteroduplex analysis for the study of phenomena associated with rapid genetic changes.     

HIV-1 Tropism Dynamics and Phylogenetic Analysis from Longitudinal Ultra-Deep Sequencing Data of CCR5- and CXCR4-Using Variants

Objective

Coreceptor switch from CCR5 to CXCR4 is associated with HIV disease progression. The molecular and evolutionary mechanisms underlying the CCR5 to CXCR4 switch are the focus of intense recent research. We studied the HIV-1 tropism dynamics in relation to coreceptor usage, the nature of quasispecies from ultra deep sequencing (UDPS) data and their phylogenetic relationships.

Methods

Here, we characterized C2-V3-C3 sequences of HIV obtained from 19 patients followed up for 54 to 114 months using UDPS, with further genotyping and phylogenetic analysis for coreceptor usage. HIV quasispecies diversity and variability as well as HIV plasma viral load were measured longitudinally and their relationship with the HIV coreceptor usage was analyzed. The longitudinal UDPS data were submitted to phylogenetic analysis and sampling times and coreceptor usage were mapped onto the trees obtained.

Results

Although a temporal viral genetic structuring was evident, the persistence of several viral lineages evolving independently along the infection was statistically supported, indicating a complex scenario for the evolution of viral quasispecies. HIV X4-using variants were present in most of our patients, exhibiting a dissimilar inter- and intra-patient predominance as the component of quasispecies even on antiretroviral therapy. The viral populations from some of the patients studied displayed evidences of the evolution of X4 variants through fitness valleys, whereas for other patients the data favored a gradual mode of emergence.

Conclusions

CXCR4 usage can emerge independently, in multiple lineages, along the course of HIV infection. The mode of emergence, i.e. gradual or through fitness valleys seems to depend on both virus and patient factors. Furthermore, our analyses suggest that, besides becoming dominant after population-level switches, minor proportions of X4 viruses might exist along the infection, perhaps even at early stages of it. The fate of these minor variants might depend on both viral and host factors.     

Role of minority populations of human immunodeficiency virus type 1 in the evolution of viral resistance to protease inhibitors

Human immunodeficiency virus type 1 (HIV-1) drug resistance results from the accumulation of mutations in the viral genes targeted by the drugs. These genetic changes, however, are commonly detected and monitored by techniques that only take into account the dominant population of plasma virus. Because HIV-1-infected patients harbor a complex and diverse mixture of virus populations, the mechanisms underlying the emergence and the evolution of resistance are not fully elucidated. Using techniques that allow the quantification of resistance mutations in minority virus species, we have monitored the evolution of resistance in plasma virus populations from patients failing protease inhibitor treatment. Minority populations with distinct resistance genotypes were detected in all patients throughout the evolution of resistance. The emergence of new dominant genotypes followed two possible mechanisms: (i) emergence of a new mutation in a currently dominant genotype and (ii) emergence of a new genotype derived from a minority virus species. In most cases, these population changes were associated with an increase in resistance at the expense of a reduction in replication capacity. Our findings provide a preliminary indication that minority viral species, which evolve independently of the majority virus population, can eventually become dominant populations, thereby serving as a reservoir of diversity and possibly accelerating the development of drug resistance.     

Genetic characterization of rebounding human immunodeficiency virus type 1 in plasma during multiple interruptions of highly active antiretroviral therapy

Various strategies of interrupting highly active antiretroviral therapy (HAART) are being investigated for the treatment of human immunodeficiency virus (HIV) infection. Interruptions of greater than 2 weeks frequently result in rebound of plasma HIV RNA. In order to discern changes in the viral population that might occur during cycles of treatment interruption, we evaluated the homology of HIV-1 envelope gene sequences over time in 12 patients who received four to seven cycles of 4 weeks off HAART followed by 8 weeks on HAART by using the heteroduplex tracking assay and novel statistical tools. HIV populations in 9 of 12 patients diverged from those found in the first cycle in at least one subsequent cycle. The substantial genetic changes noted in HIV env did not correlate with increased or decreased log changes in levels of plasma HIV RNA (P > 0.5). Thus, genetic changes in HIV env itself did not contribute in a systematic way to changes in levels of plasma viremia from cycle to cycle of treatment interruption. In addition, the data suggest that there may be multiple compartments contributing to the rebound of plasma viremia and to viral diversity from cycle to cycle of intermittent therapy.     

Multiple V1/V2 env variants are frequently present during primary infection with human immunodeficiency virus type 1

Human immunodeficiency virus type 1 (HIV-1) exists as a complex population of multiple genotypic variants in persons with chronic infection. However, acute HIV-1 infection via sexual transmission is a low-probability event in which there is thought to be low genetic complexity in the initial inoculum. In order to assess the viral complexity present during primary HIV-1 infection, the V1/V2 and V3 variable regions of the env gene were examined by using a heteroduplex tracking assay (HTA) capable of resolving these genotypic variants. Blood plasma samples from 26 primary HIV-1-infected subjects were analyzed for their level of diversity. Half of the subjects had more than one V1/V2 viral variant during primary infection, indicating the frequent transmission of multiple variants. This observation is inconsistent with the idea of infrequent transmission based on a small transmitting inoculum of cell-free virus. In chronically infected subjects, the complexity of the viral populations was even greater in both the V1/V2 and the V3 regions than in acutely infected subjects, indicating that in spite of the presence of multiple variants in acute infection, the virus does pass through a genetic bottleneck during transmission. We also examined how well the infecting virus penetrated different anatomical compartments by using the HTA. Viral variants detected in blood plasma were compared to those detected in seminal plasma and/or cerebral spinal fluid of six individuals. The virus in each of these compartments was to a large extent identical to virus in blood plasma, a finding consistent with rapid penetration of the infecting variant(s). The low-probability transmission of multiple variants could be the result of transient periods of hyperinfectiousness or hypersusceptibility. Alternatively, the inefficient transfer of a multiply infected cell could account for both the low probability of transmission and the transfer of multiple variants.     

Genetic Variation in a Human Immunodeficiency Virus Type 2 Live-Virus Macaca nemestrina Vaccine Model

Four pigtailed macaques were inoculated with an infectious, apathogenic human immunodeficiency virus type 2 (HIV-2) molecular clone (HIV-2_KR) and subsequently challenged with a highly pathogenic strain, HIV-2₂₈₇, together with two naive control animals. After challenge, two animals inoculated with a high dose of the immunizing strain were protected from CD4 decline and immunodeficiency. To examine the role of genetic heterogeneity in protection, fragments of the env gene were amplified from peripheral blood mononuclear cell DNA and plasma RNA of challenged animals by PCR, examined by using a heteroduplex tracking assay (HTA), and sequenced. By HTA, variation was detected principally within the V1 and V2 regions of envelope. Extent of variation in viral DNA clones as assessed by HTA correlated with inoculum size, as did the degree of variation in sequences of clones derived from viral DNA. Conversely, a rapid reduction in the number of plasma viral RNA variants was noted by HTA at 8 weeks postinfection in protected animals; this reduction was not present in naive or unprotected macaques. Sequences derived from plasma viral RNA were found to be more closely related than corresponding viral DNA sequences, and protection correlated with a significant reduction in variation in plasma RNA sequences in animals given the identical inocula of HIV-2₂₈₇. Nonsynonymous mutations were significantly less prevalent in the protected animals. An additional potential glycosylation site was predicted to be present in the V2 region in all but one clone, and amino acid signatures related to protection were identified in viral DNA and RNA clones within both the V1 and V2 regions. Examination of the role of viral variation in this HIV-2 live-virus vaccine model may provide valuable insights into immunopathogenesis.     

A comparison of parallel pyrosequencing and sanger clone-based sequencing and its impact on the characterization of the genetic diversity of HIV-1

Background

Pyrosequencing technology has the potential to rapidly sequence HIV-1 viral quasispecies without requiring the traditional approach of cloning. In this study, we investigated the utility of ultra-deep pyrosequencing to characterize genetic diversity of the HIV-1 gag quasispecies and assessed the possible contribution of pyrosequencing technology in studying HIV-1 biology and evolution.

Methodology/Principal Findings

HIV-1 gag gene was amplified from 96 patients using nested PCR. The PCR products were cloned and sequenced using capillary based Sanger fluorescent dideoxy termination sequencing. The same PCR products were also directly sequenced using the 454 pyrosequencing technology. The two sequencing methods were evaluated for their ability to characterize quasispecies variation, and to reveal sites under host immune pressure for their putative functional significance. A total of 14,034 variations were identified by 454 pyrosequencing versus 3,632 variations by Sanger clone-based (SCB) sequencing. 11,050 of these variations were detected only by pyrosequencing. These undetected variations were located in the HIV-1 Gag region which is known to contain putative cytotoxic T lymphocyte (CTL) and neutralizing antibody epitopes, and sites related to virus assembly and packaging. Analysis of the positively selected sites derived by the two sequencing methods identified several differences. All of them were located within the CTL epitope regions.

Conclusions/Significance

Ultra-deep pyrosequencing has proven to be a powerful tool for characterization of HIV-1 genetic diversity with enhanced sensitivity, efficiency, and accuracy. It also improved reliability of downstream evolutionary and functional analysis of HIV-1 quasispecies.     

HIV-1 Dual Infected LTNP-EC Patients Developed an Unexpected Antibody Cross-Neutralizing Activity

This study evaluated the neutralization breadth in dually infected (DI) HIV-1 long-term non-progressor elite controller patients (LTNP-EC) using a representative minipanel of 6 viruses from 5 different subtypes. Our results showed an improved neutralization breadth in DI LTNP-EC patients when compared with matched LTNP single-infected patients. The role of viral diversity in neutralization was estimated with the Shannon Entropy and the p-distance in viral quasispecies. We found a positive correlation between neutralization breadth and diversity within the viral quasispecies. This correlation could explain why a group of LTNP-EC patients developed a broad neutralizing response despite having undetectable levels of viremia.     

Whole Genome Deep Sequencing of HIV-1 Reveals the Impact of Early Minor Variants Upon Immune Recognition During Acute Infection

Deep sequencing technologies have the potential to transform the study of highly variable viral pathogens by providing a rapid and cost-effective approach to sensitively characterize rapidly evolving viral quasispecies. Here, we report on a high-throughput whole HIV-1 genome deep sequencing platform that combines 454 pyrosequencing with novel assembly and variant detection algorithms. In one subject we combined these genetic data with detailed immunological analyses to comprehensively evaluate viral evolution and immune escape during the acute phase of HIV-1 infection. The majority of early, low frequency mutations represented viral adaptation to host CD8+ T cell responses, evidence of strong immune selection pressure occurring during the early decline from peak viremia. CD8+ T cell responses capable of recognizing these low frequency escape variants coincided with the selection and evolution of more effective secondary HLA-anchor escape mutations. Frequent, and in some cases rapid, reversion of transmitted mutations was also observed across the viral genome. When located within restricted CD8 epitopes these low frequency reverting mutations were sufficient to prime de novo responses to these epitopes, again illustrating the capacity of the immune response to recognize and respond to low frequency variants. More importantly, rapid viral escape from the most immunodominant CD8+ T cell responses coincided with plateauing of the initial viral load decline in this subject, suggestive of a potential link between maintenance of effective, dominant CD8 responses and the degree of early viremia reduction. We conclude that the early control of HIV-1 replication by immunodominant CD8+ T cell responses may be substantially influenced by rapid, low frequency viral adaptations not detected by conventional sequencing approaches, which warrants further investigation. These data support the critical need for vaccine-induced CD8+ T cell responses to target more highly constrained regions of the virus in order to ensure the maintenance of immunodominant CD8 responses and the sustained decline of early viremia.