A trade-off between neutrality and adaptability limits the optimization of viral quasispecies

Theoretical studies of quasispecies usually focus on two properties of those populations at the mutation-selection equilibrium, namely asymptotic growth rate and population diversity. It has been postulated that, as a consequence of the high error rate of quasispecies replication, an increase of neutrality facilitates population optimization by reducing the amount of mutations with a deleterious effect on fitness. In this study we analyse how the optimization of equilibrium properties is affected when a quasispecies evolves in an environment perturbed through frequent bottleneck events. By means of a simple model we demonstrate that high neutrality may be detrimental when the population has to overcome repeated reductions in the population size, and that the property to be optimized in this situation is the time required to regenerate the quasispecies, i.e. its adaptability. In the scenario described, neutrality and adaptability cannot be simultaneously optimized. When fitness is equated with long-term survivability, high neutrality is the appropriate strategy in constant environments, while populations evolving in fluctuating environments are fitter when their neutrality is low, such that they can respond faster to perturbations. Our results might be relevant to better comprehend how a minority virus could displace the circulating quasispecies, a fact observed in natural infections and essential in viral evolution.     

Quantitation of relative fitness and great adaptability of clonal populations of RNA viruses.

We describe a sensitive, internally controlled method for comparing the genetic adaptability and relative fitness of virus populations in constant or changing host environments. Certain monoclonal antibody-resistant mutants of vesicular stomatitis virus can compete equally during serial passages in mixtures with the parental wild-type clone from which they were derived. These genetically marked "surrogate wild-type" neutral mutants, when mixed with wild-type virus, allow reliable measurement of changes in virus fitness and of virus adaptation to different host environments. Quantitative fitness vector plots demonstrate graphically that even clones of an RNA virus are composed of complex variant populations (quasispecies). Variants of greater fitness (competitive replication ability) were selected within very few passages of virus clones in new host cells or animals. Even clones which were well adapted to BHK21 cells gained further fitness during repeated passages in BHK21 cells.     

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.     

Genetic bottlenecks and population passages cause profound fitness differences in RNA viruses.

Repeated clone-to-clone (genetic bottleneck) passages of an RNA phage and vesicular stomatitis virus have been shown previously to result in loss of fitness due to Muller's ratchet. We now demonstrate that Muller's ratchet also operates when genetic bottleneck passages are carried out at 37 rather than 32 degrees C. Thus, these fitness losses do not depend on growth of temperature-sensitive (ts) mutants at lowered temperatures. We also demonstrate that during repeated genetic bottleneck passages, accumulation of deleterious mutations does occur in a stepwise (ratchet-like) manner as originally proposed by Muller. One selected clone which had undergone significant loss of fitness after only 20 genetic bottleneck passages was passaged again in clone-to-clone series. Additional large losses of fitness were observed in five of nine independent bottleneck series; the relative fitnesses of the other four series remained close to the starting fitness. In sharp contrast, when the same selected clone was transferred 20 more times as large populations (10(5) to 10(6) PFU transferred at each passage), significant increases in fitness were observed in all eight passage series. Finally, we selected several clones which had undergone extreme losses of fitness during 20 bottleneck passages. When these low-fitness clones were passaged many times as large virus populations, they always regained very high relative fitness. We conclude that transfer of large populations of RNA viruses regularly selects those genomes within the quasispecies population which have the highest relative fitness, whereas bottleneck transfers have a high probability of leading to loss of fitness by random isolation of genomes carrying debilitating mutations. Both phenomena arise from, and underscore, the extreme mutability and variability of RNA viruses.     

Fitness analyses of vesicular stomatitis strains with rearranged genomes reveal replicative disadvantages

Gene expression of the nonsegmented negative-strand RNA viruses is determined by the position of each gene relative to that the single 3' promoter. The general order of genes among all of the viruses of the order Mononegavirales is highly conserved. In previous work we generated recombinant viruses in which the order of the three central genes of the prototypical rhabdovirus, vesicular stomatitis virus, was rearranged to all six possible permutations. While some of these viruses replicated less well than the wild type when assayed by single-step growth analyses in BSC-1 cells, others replicated as well or slightly better. In the work reported here, we used competition assays to compare the fitness of the viruses with alternative gene orders to that of the wild-type (wt) virus. We found that the relative fitness of these recombinant viruses depended on the multiplicity of infection (MOI) but not on the population size. However, during competitions at low MOI, when complementation cannot compensate for the defects of the populations with rearranged genomes, the virus with the wt gene order was always the most fit.     

Mutagenesis-induced, large fitness variations with an invariant arenavirus consensus genomic nucleotide sequence

Enhanced mutagenesis may result in RNA virus extinction, but the molecular events underlying this process are not well understood. Here we show that 5-fluorouracil (FU)-induced mutagenesis of the arenavirus lymphocytic choriomeningitis virus (LCMV) resulted in preextinction populations whose consensus genomic nucleotide sequence remained unaltered. Furthermore, fitness recovery passages in the absence of FU, or alternate virus passages in the presence and absence of FU, led to profound differences in the capacity of LCMV to produce progeny, without modification of the consensus genomic sequence. Molecular genetic analysis failed to produce evidence of hypermutated LCMV genomes. The results suggest that low-level mutagenesis to enrich the viral population with defector, interfering genomes harboring limited numbers of mutations may mediate the loss of infectivity that accompanies viral extinction.     

Subclonal components of consensus fitness in an RNA virus clone.

Most RNA virus populations exhibit extremely high mutation frequencies which generate complex, genetically heterogeneous populations referred to as quasi-species. Previous work has shown that when a large spectrum of the quasi-species is transferred, natural selection operates, leading to elimination of noncompetitive (inferior) genomes and rapid gains in fitness. However, whenever the population is repeatedly reduced to a single virion, variable declines in fitness occur as predicted by the Muller's ratchet hypothesis. Here, we quantitated the fitness of 98 subclones isolated from an RNA virus clonal population. We found a normal distribution around a lower fitness, with the average subclone being less fit than the parental clonal population. This finding demonstrates the phenotypic diversity in RNA virus populations and shows that, as expected, a large fraction of mutations generated during virus replication is deleterious. This clarifies the operation of Muller's ratchet and illustrates why a large number of virions must be transferred for rapid fitness gains to occur. We also found that repeated genetic bottleneck passages can cause irregular stochastic declines in fitness, emphasizing again the phenotypic heterogeneity present in RNA virus populations. Finally, we found that following only 60 h of selection (15 passages in which virus yields were harvested after 4 h), RNA virus populations can undergo a 250% average increase in fitness, even on a host cell type to which they were already well adapted. This is a remarkable ability; in population biology, even a much lower fitness gain (e.g., 1 to 2%) can represent a highly significant reproductive advantage. We discuss the biological implications of these findings for the natural transmission and pathogenesis of RNA viruses.     

The order of events on Avian coronavirus life cycle shapes the order of quasispecies evolution during host switches

To understand the population genetics events during coronavirus host switches, the Beaudette strain of Avian coronavirus (AvCoV) adapted to BHK-21 cells was passaged 15 times in VERO cells, the virus load and the variants at each passage being determined by RT-qPCR and genome-length deep sequencing. From BHK-21 P2 to VERO P3, a trend for the extinction of variants was followed by stability up to VERO P11 and both the emergence and the rise in frequency in some variants, while the virus loads were stable up to VERO P12. At the spillover from BHK-21 to VERO cells, variants that both emerged, showed a rise in frequency or were extinguished were detected on the spike, while variants at the M gene showed the same pattern only at VERO passage 13. Furthermore, nsps 3-5, 9 and 15 variants were detected at lower passages compared to the consensus sequences, with those at nsp3 being detected in the spectra also at higher passages. This suggests that quasispecies coronavirus evolution in spillovers follows the virus life cycle, starting with the evolution of the receptor binding proteins, followed by the replicase and then proteins involved in virion assembly, keeping the general fitness of the mutant spectrum stable.Keywords: Coronavirus, quasispecies, cell culture, spillover     

Unusual distribution of mutations associated with serial bottleneck passages of human immunodeficiency virus type 1

Repeated bottleneck passages result in fitness losses of RNA viruses. In the case of human immunodeficiency virus type 1 (HIV-1), decreases in fitness after a limited number of plaque-to-plaque transfers in MT-4 cells were very drastic. Here we report an analysis of entire genomic nucleotide sequences of four HIV-1 clones derived from the same HIV-1 isolate and their low-fitness progeny following 7 to 15 plaque-to-plaque passages. Clones accumulated 4 to 28 mutations per genome, with dominance of A --> G and G --> A transitions (57% of all mutations) and 49% nonsynonymous replacements. One clone-but not three sibling clones-showed an overabundance of G --> A transitions, evidencing the highly stochastic nature of some types of mutational bias. The distribution of mutations along the genome was very unusual in that mutation frequencies in gag were threefold higher than in env. Particularly striking was the complete absence of replacements in the V3 loop of gp120, confirmed with partial nucleotide sequences of additional HIV-1 clones subjected to repeated bottleneck passages. The analyses revealed several amino acid replacements that have not been previously recorded among natural HIV-1 isolates and illustrate how evolution of an RNA virus genome, with regard to constant and variable regions, can be profoundly modified by alterations in population dynamics.     

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.     

Viral fitness can influence the repertoire of virus variants selected by antibodies

Minority genomes in the mutant spectra of viral quasispecies may differ in relative fitness. Here, we report experiments designed to evaluate the contribution of relative fitness to selection by a neutralizing monoclonal antibody (mAb). We have reconstructed a foot-and-mouth disease virus (FMDV) quasispecies, with two matched pairs of distinguishable mAb-escape mutants as minority genomes of the mutant spectrum. Each mutant of a pair differs from the other by 11-fold or 33-fold in relative fitness. Analysis of the mutant spectra of virus populations selected with different concentrations of antibody in infections in liquid culture medium has documented a dominance of the high fitness counterpart in the selected population. Plaque development as a function of increasing concentration of the antibody has shown that each mutant of a matched pair yielded the same number of plaques, although the high fitness mutant required less time for plaque formation, and attained a larger plaque size at any given time-point. This result documents equal intrinsic resistance to the antibody of each mutant of a matched pair, confirming previous biochemical, structural, and genetic studies, which indicated that the epitopes of each mutant pair were indistinguishable regarding reactivity with the monoclonal antibody. Thus, relative viral fitness can influence in a significant way the repertoire of viral mutants selected from a viral quasispecies by a neutralizing antibody. We discuss the significance of these results in relation to antibody selection, and to other selective forces likely encountered by viral quasispecies in vivo.     

Evolution of a persistent aphthovirus in cytolytic infections: partial reversion of phenotypic traits accompanied by genetic diversification.

Foot-and-mouth disease virus (FMDV) shows a dual potential to be cytolytic or to establish persistent infections in cell culture. FMDV R100, a virus rescued after 100 passages of carrier BHK-21 cells persistently infected with FMDV clone C-S8c1, showed multiple genetic and phenotypic alterations relative to the parental clone C-S8c1. Several FMDV R100 populations have been subjected to 100 serial cytolytic infections in BHK-21 cells, and the reversion of phenotypic and genetic alterations has been analyzed. An extreme temperature sensitivity of R100 reverted totally or partially in some passage series but not in others. The small-plaque morphology reverted to normal size in all cases. The hypervirulence for BHK-21 cells did not revert, and even showed an increase, upon cytolytic passage. Most of the mutations that had been fixed in the R100 genome during persistence did not revert in the course of cytolytic passages, but the extended polyribocytidylate tract of R100 (about 460 residues, versus 290 in C-S8c1) decreased dramatically in length, to the range of 220 to 260 residues in all passage series examined. In passages involving very large viral populations, a variant with two amino acid substitutions (L-144-->V and A-145-->P) next to the highly conserved Arg-Gly-Asp (RGD motif; positions 141 to 143) within the G-H loop of capsid protein VP1 became dominant. A clonal analysis allowed isolation of a mutant with the single replacement A-145-->P. Viral production and growth competition experiments showed the two variants to have a fitness very close to that of the parental virus. The results provide evidence that the repertoire of variants that could potentially become dominant in viral quasispecies may be influenced by the population size of the evolving virus. The net results of a series of persistent-infection passages followed by a series of cytolytic passages was progressive genomic diversification despite reversion or stasis of phenotypic traits. Implications for the evolution of RNA viruses are discussed.     

Duration and fitness dependence of quasispecies memory.

The duration and fitness dependence of memory in viral quasispecies evolving in cell culture have been investigated using two genetic markers of foot-and-mouth disease virus (FMDV). In lineages of antigenic variant FMDV RED, which reverted to FMDV RGD, memory FMDV RED genomes were detected after 50 infectious cycles, and memory level was fitness dependent. In growth-competition experiments between a reference FMDV RGD and two different FMDV RED populations, a 7.6-fold higher fitness of the initial FMDV RED population resulted in 30 to 100-fold higher memory level. In lineages of low-fitness clones containing an elongated internal polyadenylate tract, revertants lacking excess adenylate residues became dominant by passage 20. However, genomes including a larger number of adenylate residues were detected as memory genomes after at least 150 infectious cycles. Thus, quasispecies memory can be durable and is fitness dependent, as predicted from the growth competition of two mutant forms of a genome. An understanding of factors influencing quasispecies memory levels and duration may have implications for the extended diagnosis of viruses based on the quantification of minority genomes.     

Adaptation of Puumala hantavirus to cell culture is associated with point mutations in the coding region of the L segment and in the noncoding regions of the S segment

We previously developed a model for studies on hantavirus host adaptation and initiated genetic analysis of Puumala virus variants passaged in colonized bank voles and in cultured Vero E6 cells. With the data presented in this paper, the sequence comparison of the wild-type and Vero E6-adapted variants of Puumala virus, strain Kazan, has been completed. The only amino acid substitution that distinguished the two virus variants was found in the L protein, Ser versus Phe at position 2053. Another mutation found in the L segment, the silent transition C1053U, could result from the selection of a variant with altered L RNA folding. Nucleotide substitutions observed in individual L cDNA clones, most of them A-->G and U-->C transitions, suggested that the population of L RNA molecules is represented by quasispecies. The mutation frequency in the L segment quasispecies appeared to be similar to the corresponding values for the S and M quasispecies. Analysis of the cDNA clones with the complete S segment sequences from passage 20 confirmed our earlier conclusion that the cell-adapted genotype of the virus is represented mostly by variants with mutated S segment noncoding regions. However, the spectrum of the S segment quasispecies appeared to be changing, suggesting that, after the initial adaptation (passages 1 to 11), the viral population is still being driven by selection for variants with higher fitness.     

Red queen dynamics, competition and critical points in a model of RNA virus quasispecies.

RNA viruses offer a unique opportunity for the study of evolution at the molecular level. Recent experiments involving clonal populations of RNA viruses have shown that competition among virus strains of approximately equal relative fitness can result in the eventual competitive exclusion of one of the species. As competition proceeds in time, both the winners and the losers exhibited absolute gains in fitness, consistent with the "Red Queen" hypothesis of evolution. Further experiments involving closely related evolving quasispecies revealed a highly predictable nonlinear behavior suggesting a deterministic component in the underlying quasispecies dynamics. This is apparently in contradiction with the standard view of RNA virus evolution as a highly unpredictable process. In this paper we present a simple model which allows previous hypothesis to be tested and provides an interpretation for the observed experimental results.     

The quasispecies (extremely heterogeneous) nature of viral RNA genome populations: biological relevance--a review

We review evidence that cloned (or uncloned) populations of most RNA viruses do not consist of a single genome species of defined sequence, but rather of heterogeneous mixtures of related genomes (quasispecies). Due to very high mutation rates, genomes of a quasispecies virus population share a consensus sequence but differ from each other and from the consensus sequence by one, several, or many mutations. Viral genome analyses by sequencing, fingerprinting, cDNA cloning etc. indicate that most viral RNA populations (quasispecies) contain all possible single and double genomic site mutations and varying proportions of triple, quadruple, etc. site mutations. This quasispecies structure of RNA virus populations has many important theoretical and practical implications because mutations at only one or a few sites may alter the phenotype of an RNA virus.     

EVOLUTIONARY DYNAMICS OF FITNESS RECOVERY FROM THE DEBILITATING EFFECTS OF MULLER'S RATCHET

The great adaptability shown by RNA viruses is a consequence of their high mutation rates. The evolution of fitness in a severely debilitated, clonal population of the nonsegmented ribovirus vesicular stomatitis virus (VSV) has been compared under five different demographic regimes, ranging from severe serial bottleneck passages (one virion) to large population passages (10⁵ virions or more) under similar environmental conditions (cell culture type and temperature). No matter how small the bottleneck, the fitness of the evolved populations was always higher than the fitness of the starting population; this result is clearly different from that previously reported for viruses with higher fitness. The reattainment of fitness under a regime of serial population passages showed two main characteristics: (1) the rate of adaptation was higher during early passages; and (2) a maximum fitness value was reached after a large number of passages. The maximum fitness reached by this initially debilitated clone was similar to the fitness of wild-type virus. The practical implications of these findings in the design of vaccines using attenuated viruses are also discussed.     

Delayed transmission selects for increased survival of vesicular stomatitis virus

Life-history theory predicts that traits for survival and reproduction cannot be simultaneously maximized in evolving populations. For this reason, in obligate parasites such as infectious viruses, selection for improved between-host survival during transmission may lead to evolution of decreased within-host reproduction. We tested this idea using experimental evolution of RNA virus populations, passaged under differing transmission times in the laboratory. A single ancestral genotype of vesicular stomatitis virus (VSV), a negative-sense RNA Rhabdovirus, was used to found multiple virus lineages evolved in either ordinary 24-h cell-culture passage, or in delayed passages of 48 h. After 30 passages (120 generations of viral evolution), we observed that delayed transmission selected for improved extracellular survival, which traded-off with lowered viral fecundity (slower exponential population growth and smaller mean plaque size). To further examine the confirmed evolutionary trade-off, we obtained consensus whole-genome sequences of evolved virus populations, to infer phenotype–genotype associations. Results implied that increased virus survival did not occur via convergence; rather, improved virion stability was gained via independent mutations in various VSV structural proteins. Our study suggests that RNA viruses can evolve different molecular solutions for enhanced survival despite their limited genetic architecture, but suffer generalized reproductive trade-offs that limit overall fitness gains.