HIV-1 evolution under pressure of protease inhibitors: Climbing the stairs of viral fitness

The human immunodeficiency virus (HIV-1) has evolved into a viral quasispecies with a high replication capacity or fitness. Antiretroviral drugs potently inhibit replication of the wild-type virus, but HIV-1 responds by selection of drug-resistant variants. Here we review, in brief, the evolution of resistance to protease inhibitors that is characterized by severe fitness losses and an abundance of subsequent repair strategies. The possibility to restrict HIV-1 fitness is discussed in relation to the control of HIV-1 pathogenesis.     

Population Genetic Analysis of the Protease Locus of Human Immunodeficiency Virus Type 1 Quasispecies Undergoing Drug Selection, Using a Denaturing Gradient-Heteroduplex Tracking Assay

Monitoring the evolution of human immunodeficiency virus type 1 (HIV-1) drug resistance requires measuring the frequency of closely related genetic variants making up the complex viral quasispecies found in vivo. In order to resolve both major and minor (>/=2%) protease gene variants differing by one or more nucleotide substitutions, we analyzed PCR products derived from plasma viral quasispecies by using a combination of denaturing gradient gel electrophoresis and DNA heteroduplex tracking assays. Correct population sampling of the high level of genetic diversity present within viral quasispecies could be documented by parallel analysis of duplicate, independently generated PCR products. The composition of genetically complex protease gene quasispecies remained constant over short periods of time in the absence of treatment and while plasma viremia fell >100-fold following the initiation of protease inhibitor ritonavir monotherapy. Within a month of initiating therapy, a strong reduction in the genetic diversity of plasma viral populations at the selected protease locus was associated with rising plasma viremia and the emergence of drug resistance. The high levels of protease genetic diversity seen before treatment reemerged only months later. In one patient, reduction in genetic diversity at the protease gene was observed concomitantly with an increase in diversity at the envelope gene (E. L. Delwart, P. Heng, A. Neumann, and M. Markowitz, J. Virol. 72:2416-2421, 1998), indicating that opposite population genetic changes can take place in different HIV-1 loci. The rapid emergence of drug-resistant HIV-1 was therefore associated with a strong, although only transient, reduction in genetic diversity at the selected locus. The denaturing gradient-heteroduplex tracking assay is a simple method for the separation and quantitation of very closely related, low-frequency, genetic variants within complex viral populations.     

Can modulation of viral fitness represent a target for anti-HIV-1 strategies?

Prolonged use of anti-retroviral compounds in human immunodeficiency virus type 1 (HIV-1) infection selects for drug-resistant and often mutidrug-resistant viral variants. Drug-resistance mutations may also affect viral fitness. Interestingly, recent research has indicated that some of the unfit drug-resistant variants may be less pathogenic, suggesting that decreased viral fitness is beneficial for the host and may be driven by specific treatments during anti-HIV-1 infection. A second potential antiviral strategy starting with profound inteference with viral fitness aims at forcing viruses towards lethal mutagenesis (the so-called "error catastrophe"). This review summarizes the methods for addressing HIV-1 fitness in vitro and ex vivo, the current understanding of clinical implications of reduced HIV-1 fitness, and the potential use of anti-HIV-1 strategies aiming at modulating viral fitness. Finally, it is emphasized how the peculiar features of HIV-1 quasispecies (displaying two different forms of memory, a replicative and a non-replicative form) may sharply influence the design of future diagnostic methodologies for fitness analysis.     

Heterosexual Transmission of Subtype C HIV-1 Selects Consensus-Like Variants without Increased Replicative Capacity or Interferon-α Resistance

Heterosexual transmission of HIV-1 is characterized by a genetic bottleneck that selects a single viral variant, the transmitted/founder (TF), during most transmission events. To assess viral characteristics influencing HIV-1 transmission, we sequenced 167 near full-length viral genomes and generated 40 infectious molecular clones (IMC) including TF variants and multiple non-transmitted (NT) HIV-1 subtype C variants from six linked heterosexual transmission pairs near the time of transmission. Consensus-like genomes sensitive to donor antibodies were selected for during transmission in these six transmission pairs. However, TF variants did not demonstrate increased viral fitness in terms of particle infectivity or viral replicative capacity in activated peripheral blood mononuclear cells (PBMC) and monocyte-derived dendritic cells (MDDC). In addition, resistance of the TF variant to the antiviral effects of interferon-α (IFN-α) was not significantly different from that of non-transmitted variants from the same transmission pair. Thus neither in vitro viral replicative capacity nor IFN-α resistance discriminated the transmission potential of viruses in the quasispecies of these chronically infected individuals. However, our findings support the hypothesis that within-host evolution of HIV-1 in response to adaptive immune responses reduces viral transmission potential.     

Resistance-Associated Loss of Viral Fitness in Human Immunodeficiency Virus Type 1: Phenotypic Analysis of Protease and gag Coevolution in Protease Inhibitor-Treated Patients

We have studied the phenotypic impact of adaptative Gag cleavage site mutations in patient-derived human immunodeficiency virus type 1 (HIV-1) variants having developed resistance to the protease inhibitor ritonavir or saquinavir. We found that Gag mutations occurred in a minority of resistant viruses, regardless of the duration of the treatment and of the protease mutation profile. Gag mutations exerted only a partial corrective effect on resistance-associated loss of viral fitness. Reconstructed viruses with resistant proteases displayed multiple Gag cleavage defects, and in spite of Gag adaptation, several of these defects remained, explaining the limited corrective effect of cleavage site mutations on fitness. Our data provide clear evidence of the interplay between resistance and fitness in HIV-1 evolution in patients treated with protease inhibitors.     

Viral quasispecies profiles as the result of the interplay of competition and cooperation

Background  

Viral quasispecies can be regarded as a swarm of genetically related mutants. A common approach employed to describe viral quasispecies is by means of the quasispecies equation (QE). However, a main criticism of QE is its lack of frequency-dependent selection. This can be overcome by an alternative formulation for the evolutionary dynamics: the replicator-mutator equation (RME). In turn, a problem with the RME is how to quantify the interaction coefficients between viral variants. Here, this is addressed by adopting an ecological perspective and resorting to the niche theory of competing communities, which assumes that the utilization of resources primarily determines ecological segregation between competing individuals (the different viral variants that constitute the quasispecies). This provides a theoretical framework to estimate quantitatively the fitness landscape.     

Origin of human immunodeficiency virus type 1 quasispecies emerging after antiretroviral treatment interruption in patients with therapeutic failure

The emergence of antiretroviral (ARV) drug-resistant human immunodeficiency virus type 1 (HIV-1) quasispecies is a major cause of treatment failure. These variants are usually replaced by drug-sensitive ones when the selective pressure of the drugs is removed, as the former have reduced fitness in a drug-free environment. This was the rationale for the design of structured ARV treatment interruption (STI) studies for the management of HIV-1 patients with treatment failure. We have studied the origin of drug-sensitive HIV-1 quasispecies emerging after STI in patients with treatment failure due to ARV drug resistance. Plasma and peripheral blood mononuclear cell samples were obtained the day of treatment interruption (day 0) and 30 and 60 days afterwards. HIV-1 pol and env were partially amplified, cloned, and sequenced. At day 60 drug-resistant variants were replaced by completely or partially sensitive quasispecies. Phylogenetic analyses of pol revealed that drug-sensitive variants emerging after STI were not related to their immediate temporal ancestors but formed a separate cluster, demonstrating that STI leads to the recrudescence and reemergence of a sequestrated viral population rather than leading to the back mutation of drug-resistant forms. No evidence for concomitant changes in viral tropism was seen, as deduced from env sequences. This study demonstrates the important role that the reemergence of quasispecies plays in HIV-1 population dynamics and points out the difficulties that may be found when recycling ARV therapies with patients with treatment failure.     

Transmission of single and multiple viral variants in primary HIV-1 subtype C infection

To address whether sequences of viral gag and env quasispecies collected during the early post-acute period can be utilized to determine multiplicity of transmitted HIV's, recently developed approaches for analysis of viral evolution in acute HIV-1 infection [1,2] were applied. Specifically, phylogenetic reconstruction, inter- and intra-patient distribution of maximum and mean genetic distances, analysis of Poisson fitness, shape of highlighter plots, recombination analysis, and estimation of time to the most recent common ancestor (tMRCA) were utilized for resolving multiplicity of HIV-1 transmission in a set of viral quasispecies collected within 50 days post-seroconversion (p/s) in 25 HIV-infected individuals with estimated time of seroconversion. The decision on multiplicity of HIV infection was made based on the model's fit with, or failure to explain, the observed extent of viral sequence heterogeneity. The initial analysis was based on phylogeny, inter-patient distribution of maximum and mean distances, and Poisson fitness, and was able to resolve multiplicity of HIV transmission in 20 of 25 (80%) cases. Additional analysis involved distribution of individual viral distances, highlighter plots, recombination analysis, and estimation of tMRCA, and resolved 4 of the 5 remaining cases. Overall, transmission of a single viral variant was identified in 16 of 25 (64%) cases, and transmission of multiple variants was evident in 8 of 25 (32%) cases. In one case multiplicity of HIV-1 transmission could not be determined. In primary HIV-1 subtype C infection, samples collected within 50 days p/s and analyzed by a single-genome amplification/sequencing technique can provide reliable identification of transmission multiplicity in 24 of 25 (96%) cases. Observed transmission frequency of a single viral variant and multiple viral variants were within the ranges of 64% to 68%, and 32% to 36%, respectively.     

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.     

Initial Fitness Recovery of HIV-1 Is Associated with Quasispecies Heterogeneity and Can Occur without Modifications in the Consensus Sequence

Background

Fitness recovery of HIV-1 “in vitro” was studied using viral clones that had their fitness decreased as a result of plaque-to-plaque passages.

Principal Findings

After ten large population passages, the viral populations showed an average increase of fitness, although with wide variations among clones. While 5 clones showed significant fitness increases, 3 clones showed increases that were only marginally significant (p<0.1), and 4 clones did not show any change. Fitness recovery was not accompanied by an increase in p24 production, but was associated with an increase in viral titer. Few mutations (an average of 2 mutations per genome) were detected in the consensus nucleotide sequence of the entire genome in all viral populations. Five of the populations did not fix any mutation, and three of them displayed marginally significant fitness increases, illustrating that fitness recovery can occur without detectable alterations of the consensus genomic sequence. The investigation of other possible viral factors associated with the initial steps of fitness recovery, showed that viral quasispecies heterogeneity increased between the initial clones and the passaged populations. A direct statistical correlation between viral heterogeneity and viral fitness was obtained.

Conclusions

Thus, the initial fitness recovery of debilitated HIV-1 clones was mediated by an increase in quasispecies heterogeneity. This observation, together with the invariance of the consensus sequence despite fitness increases demonstrates the relevance of quasispecies heterogeneity in the evolution of HIV-1 in cell culture.     

Fitness increase of memory genomes in a viral quasispecies

Viral quasispecies may contain a subset of minority genomes that reflect those genomic sequences that were dominant at an early phase of quasispecies evolution. Such minority genomes are referred to as memory in viral quasispecies. A memory marker previously characterized in foot-and-mouth disease virus (FMDV) is an internal oligoadenylate tract of variable length that became dominant upon serial plaque-to-plaque transfers of FMDV clones. During large population passages, genomes with internal oligoadenylate were outcompeted by wild-type revertants but remained in the mutant spectra as memory genomes. Here, we report a quantification of relative fitness of several FMDV clones, harboring internal oligoadenylate tracts of different length, and that were retrieved at early or late times (passage number) after implementation of memory. The results show that for any given length range of the oligoadenylate, maintenance in memory resulted in an increase in relative fitness, comparable to the increase undergone by the entire population. The fitness increase is in agreement with the Red Queen hypothesis, and implies a replicative memory mechanism. Thus, permanence of memory genomes may be a source of high fitness variants despite their initial low fitness, and despite having remained hidden in mutant spectra. This reinforces the interest of diagnosing minority genomes during chronic human and animal viral infections.     

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.     

Impaired fitness of human immunodeficiency virus type 1 variants with high-level resistance to protease inhibitors.

One hope to maintain the benefits of antiviral therapy against the human immunodeficiency virus type 1 (HIV-1), despite the development of resistance, is the possibility that resistant variants will show decreased viral fitness. To study this possibility, HIV-1 variants showing high-level resistance (up to 1,500-fold) to the substrate analog protease inhibitors BILA 1906 BS and BILA 2185 BS have been characterized. Active-site mutations V32I and I84V/A were consistently observed in the protease of highly resistant viruses, along with up to six other mutations. In vitro studies with recombinant mutant proteases demonstrated that these mutations resulted in up to 10(4)-fold increases in the Ki values toward BILA 1906 BS and BILA 2185 BS and a concomitant 2,200-fold decrease in catalytic efficiency of the enzymes toward a synthetic substrate. When introduced into viral molecular clones, the protease mutations impaired polyprotein processing, consistent with a decrease in enzyme activity in virions. Despite these observations, however, most mutations had little effect on viral replication except when the active-site mutations V32I and I84V/A were coexpressed in the protease. The latter combinations not only conferred a significant growth reduction of viral clones on peripheral blood mononuclear cells but also caused the complete disappearance of mutated clones when cocultured with wild-type virus on T-cell lines. Furthermore, the double nucleotide mutation I84A rapidly reverted to I84V upon drug removal, confirming its impact on viral fitness. Therefore, high-level resistance to protease inhibitors can be associated with impaired viral fitness, suggesting that antiviral therapies with such inhibitors may maintain some clinical benefits.     

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.     

Genetic characterization of human immunodeficiency virus type 1 in blood and genital secretions: evidence for viral compartmentalization and selection during sexual transmission.

To explore the mechanism of sexual transmission of human immunodeficiency virus type 1 (HIV-1), we compared HIV-1 gp120 sequences in longitudinal samples from five acute seroconvertors with those from their corresponding sexual partners (transmitters). We used a quantitative homoduplex tracking assay to compare the overall genetic composition of HIV-1 quasispecies in each transmission pair and to track the transmitted viruses during the acute and asymptomatic stages of HIV-1 infection. In the chronically infected transmitters, HIV-1 variants in genital secretions differed from those in blood and variants in cells differed from those in cell-free plasma, indicating remarkable sequence heterogeneity in these subjects as well as compartmentalization of the virus in different bodily sites. Conversely, two of five seroconvertors had only a few related variants and three of five harbored only one viral population, indicating that in these subjects the transmitted viruses were typically homogeneous. Transmitted viruses were evident in the donor's seminal plasma (one of five cases) and even more so in their seminal cells (three of five cases), suggesting that both cell-associated and cell-free viruses can be transmitted. In every pair studied, the transmitted variant(s) represents only a minor population in the semen of the corresponding transmitter, thereby providing evidence that HIV-1 selection indeed occurs during sexual transmission.     

Sequential immunization with a subtype B HIV-1 envelope quasispecies partially mimics the in vivo development of neutralizing antibodies

A major goal of human immunodeficiency virus type 1 (HIV-1) vaccine efforts is the design of Envelope (Env)-based immunogens effective at eliciting heterologous or broad neutralizing antibodies (NAbs). We hypothesized that programming the B-cell response could be achieved by sequentially exposing the host to a collection of env variants representing the viral quasispecies members isolated from an individual that developed broad NAbs over time. This ordered vaccine approach (sequential) was compared to exposure to a cocktail of env clones (mixture) and to a single env variant (clonal). The three strategies induced comparable levels of the autologous and heterologous neutralization of tier 1 pseudoviruses. Sequential and mixture exposure to quasispecies led to epitope targeting similar to that observed in the simian-human immunodeficiency virus (SHIV)-infected animal from which the env variants were cloned, while clonal and sequential exposure led to greater antibody maturation than the mixture. Therefore, the sequential vaccine approach best replicated the features of the NAb response observed in that animal. This study is the first to explore the use of a collection of HIV-1 env quasispecies variants as immunogens and to present evidence that it is possible to educate the B-cell response by sequential exposure to native HIV-1 quasispecies env variants derived from an individual with a broadened NAb response.     

Quasispecies spatial models for RNA viruses with different replication modes and infection strategies

Empirical observations and theoretical studies suggest that viruses may use different replication strategies to amplify their genomes, which impact the dynamics of mutation accumulation in viral populations and therefore, their fitness and virulence. Similarly, during natural infections, viruses replicate and infect cells that are rarely in suspension but spatially organized. Surprisingly, most quasispecies models of virus replication have ignored these two phenomena. In order to study these two viral characteristics, we have developed stochastic cellular automata models that simulate two different modes of replication (geometric vs stamping machine) for quasispecies replicating and spreading on a two-dimensional space. Furthermore, we explored these two replication models considering epistatic fitness landscapes (antagonistic vs synergistic) and different scenarios for cell-to-cell spread, one with free superinfection and another with superinfection inhibition. We found that the master sequences for populations replicating geometrically and with antagonistic fitness effects vanished at low critical mutation rates. By contrast, the highest critical mutation rate was observed for populations replicating geometrically but with a synergistic fitness landscape. Our simulations also showed that for stamping machine replication and antagonistic epistasis, a combination that appears to be common among plant viruses, populations further increased their robustness by inhibiting superinfection. We have also shown that the mode of replication strongly influenced the linkage between viral loci, which rapidly reached linkage equilibrium at increasing mutations for geometric replication. We also found that the strategy that minimized the time required to spread over the whole space was the stamping machine with antagonistic epistasis among mutations. Finally, our simulations revealed that the multiplicity of infection fluctuated but generically increased along time.     

Virus genome dynamics under different propagation pressures: reconstruction of whole genome haplotypes of west nile viruses from NGS data

Background

Extensive focus is placed on the comparative analyses of consensus genotypes in the study of West Nile virus (WNV) emergence. Few studies account for genetic change in the underlying WNV quasispecies population variants. These variants are not discernable in the consensus genome at the time of emergence, and the maintenance of mutation-selection equilibria of population variants is greatly underestimated. The emergence of lineage 1 WNV strains has been studied extensively, but recent epidemics caused by lineage 2 WNV strains in Hungary, Austria, Greece and Italy emphasizes the increasing importance of this lineage to public health. In this study we explored the quasispecies dynamics of minority variants that contribute to cell-tropism and host determination, i.e. the ability to infect different cell types or cells from different species from Next Generation Sequencing (NGS) data of a historic lineage 2 WNV strain.

Results

Minority variants contributing to host cell membrane association persist in the viral population without contributing to the genetic change in the consensus genome. Minority variants are shown to maintain a stable mutation-selection equilibrium under positive selection, particularly in the capsid gene region.

Conclusions

This study is the first to infer positive selection and the persistence of WNV haplotype variants that contribute to viral fitness without accompanying genetic change in the consensus genotype, documented solely from NGS sequence data. The approach used in this study streamlines the experimental design seeking viral minority variants accurately from NGS data whilst minimizing the influence of associated sequence error.