Replication at periodically changing multiplicity of infection promotes stable coexistence of competing viral populations

Effects of historical co-infection on host shift abilities of exploitative and competitive viruses

Rapid evolution contributes to frequent emergence of RNA viral pathogens on novel hosts. However, accurately predicting which viral genotypes will emerge has been elusive. Prior work with lytic RNA bacteriophage ɸ6 (family Cystoviridae) suggested that evolution under low multiplicity of infection (MOI; proportion of viruses to susceptible cells) selected for greater host exploitation, while evolution under high MOI selected for better intracellular competition against co-infecting viruses. We predicted that phage genotypes that had experienced 300 generations of low MOI ecological history would be relatively advantaged in initial growth on two novel hosts. We inferred viral growth through changes in host population density, specifically by analyzing five attributes of growth curves of infected bacteria. Despite equivalent growth of evolved viruses on the original host, low MOI evolved clones were generally advantaged relative to high MOI clones in exploiting novel hosts. However, the specific attributes of growth curves that supported their advantage differed by host, indicating interactions between both viral and host genotype. Although there will be host specificity in viral growth, we suggest based on infectivity differences of viruses from high versus low MOI histories that prior MOI selection can later affect emergence potential.

Competition between social cheater viruses is driven by mechanistically different cheating strategies

Cheater viruses, also known as defective interfering viruses, cannot replicate on their own yet replicate faster than the wild type upon coinfection. While there is growing interest in using cheaters as antiviral therapeutics, the mechanisms underlying cheating have been rarely explored. During experimental evolution of MS2 phage, we observed the parallel emergence of two independent cheater mutants. The first, a point deletion mutant, lacked polymerase activity but was advantageous in viral packaging. The second synonymous mutant cheater displayed a completely different cheating mechanism, involving an altered RNA structure. Continued evolution revealed the demise of the deletion cheater and rise of the synonymous cheater. A mathematical model inferred that while a single cheater is expected to reach an equilibrium with the wild type, cheater demise arises from antagonistic interactions between coinfecting cheaters. These findings highlight layers of parasitism: viruses parasitizing cells, cheaters parasitizing intact viruses, and cheaters may parasitize other cheaters.

Beneficial coinfection can promote within-host viral diversity

In many viral infections, a large number of different genetic variants can coexist within a host, leading to more virulent infections that are better able to evolve antiviral resistance and adapt to new hosts. But how is this diversity maintained? Why do faster-growing variants not outcompete slower-growing variants, and erode this diversity? One hypothesis is if there are mutually beneficial interactions between variants, with host cells infected by multiple different viral genomes producing more, or more effective, virions. We modelled this hypothesis with both mathematical models and simulations, and found that moderate levels of beneficial coinfection can maintain high levels of coexistence, even when coinfection is relatively rare, and when there are significant fitness differences between competing variants. Rare variants are more likely to be coinfecting with a different variant, and hence beneficial coinfection increases the relative fitness of rare variants through negative frequency dependence, and maintains diversity. We further find that coexisting variants sometimes reach unequal frequencies, depending on the extent to which different variants benefit from coinfection, and the ratio of variants which leads to the most productive infected cells. These factors could help drive the evolution of defective interfering particles, and help to explain why the different segments of multipartite viruses persist at different equilibrium frequencies.

Cooperation between distinct viral variants promotes growth of H3N2 influenza in cell culture

RNA viruses rapidly diversify into quasispecies of related genotypes. This genetic diversity has long been known to facilitate adaptation, but recent studies have suggested that cooperation between variants might also increase population fitness. Here, we demonstrate strong cooperation between two H3N2 influenza variants that differ by a single mutation at residue 151 in neuraminidase, which normally mediates viral exit from host cells. Residue 151 is often annotated as an ambiguous amino acid in sequenced isolates, indicating mixed viral populations. We show that mixed populations grow better than either variant alone in cell culture. Pure populations of either variant generate the other through mutation and then stably maintain a mix of the two genotypes. We suggest that cooperation arises because mixed populations combine one variant's proficiency at cell entry with the other's proficiency at cell exit. Our work demonstrates a specific cooperative interaction between defined variants in a viral quasispecies.

The Multiplicity of Cellular Infection Changes Depending on the Route of Cell Infection in a Plant Virus

The multiplicity of cellular infection (MOI) is the number of virus genomes of a given virus species that infect individual cells. This parameter chiefly impacts the severity of within-host population bottlenecks as well as the intensity of genetic exchange, competition, and complementation among viral genotypes. Only a few formal estimations of the MOI currently are available, and most theoretical reports have considered this parameter as constant within the infected host. Nevertheless, the colonization of a multicellular host is a complex process during which the MOI may dramatically change in different organs and at different stages of the infection. We have used both qualitative and quantitative approaches to analyze the MOI during the colonization of turnip plants by Turnip mosaic virus. Remarkably, different MOIs were observed at two phases of the systemic infection of a leaf. The MOI was very low in primary infections from virus circulating within the vasculature, generally leading to primary foci founded by a single genome. Each lineage then moved from cell to cell at a very high MOI. Despite this elevated MOI during cell-to-cell progression, coinfection of cells by lineages originating in different primary foci is severely limited by the rapid onset of a mechanism inhibiting secondary infection. Thus, our results unveil an intriguing colonization pattern where individual viral genomes initiate distinct lineages within a leaf. Kin genomes then massively coinfect cells, but coinfection by two distinct lineages is strictly limited.

IMPORTANCE

The MOI is the size of the viral population colonizing cells and defines major phenomena in virus evolution, like the intensity of genetic exchange and the size of within-host population bottlenecks. However, few studies have quantified the MOI, and most consider this parameter as constant during infection. Our results reveal that the MOI can depend largely on the route of cell infection in a systemically infected leaf. The MOI is usually one genome per cell when cells are infected from virus particles moving long distances in the vasculature, whereas it is much higher during subsequent cell-to-cell movement in mesophyll. However, a fast-acting superinfection exclusion prevents cell coinfection by merging populations originating from different primary foci within a leaf. This complex colonization pattern results in a situation where within-cell interactions are occurring almost exclusively among kin and explains the common but uncharacterized phenomenon of genotype spatial segregation in infected plants.

Congruent evolution of fitness and genetic robustness in vesicular stomatitis virus

Quasispecies theory is a case of mutation-selection balance for evolution at high mutation rates, such as those observed in RNA viruses. One of the main predictions of this model is the selection for robustness, defined as the ability of an organism to remain phenotypically unchanged in the face of mutation. We have used a collection of vesicular stomatitis virus strains that had been evolving either under positive selection or under random drift. We had previously shown that the former increase in fitness while the latter have overall fitness decreases (I. S. Novella, J. B. Presloid, T. Zhou, S. D. Smith-Tsurkan, B. E. Ebendick-Corpus, R. N. Dutta, K. L. Lust, and C. O. Wilke, J. Virol. 84:4960-4968, 2010). Here, we determined the robustness of these strains and demonstrated that strains under positive selection not only increase in fitness but also increase in robustness. In contrast, strains under drift not only decreased in fitness but also decreased in robustness. There was a good overall correlation between fitness and robustness. We also tested whether there was a correlation between fitness and thermostability, and we observed that the correlation was imperfect, indicating that the fitness effects of mutations are exerted in part at a level other than changing the resistance of the protein to temperature.

The role of environmental factors on the evolution of phenotypic diversity in vesicular stomatitis virus populations

Virus adaptation to an ever-changing environment requires the availability of variants with phenotypes that can fulfil new requirements for replication. High mutation rates result in the generation of these variants. The factors that contribute to the maintenance or elimination of this diversity, however, are not fully understood. This study used a collection of vesicular stomatitis virus strains generated under different conditions to measure the extent of variation within each population, and tested the effects of several environmental factors on diversity. It was found that the host-cell type used for selection sometimes had an effect on the extent of variation and that there may be different levels of variation over time. Persistent infections promoted higher levels of diversity than acute infections, presumably due to complementation. In contrast, environmental heterogeneity, host breadth and the cell type used for testing (as opposed to the cell type used for selection) did not seem to have an effect on the amount of phenotypic diversity observed.

Arenaviruses and lethal mutagenesis. Prospects for new ribavirin-based interventions

Lymphocytic choriomeningitis virus (LCMV) has contributed to unveil some of the molecular mechanisms of lethal mutagenesis, or loss of virus infectivity due to increased mutation rates. Here we review these developments, and provide additional evidence that ribavirin displays a dual mutagenic and inhibitory activity on LCMV that can be relevant to treatment designs. Using 5-fluorouracil as mutagenic agent and ribavirin either as inhibitor or mutagen, we document an advantage of a sequential inhibitor-mutagen administration over the corresponding combination treatment to achieve a low LCMV load in cell culture. This advantage is accentuated in the concentration range in which ribavirin acts mainly as an inhibitor, rather than as mutagen. This observation reinforces previous theoretical and experimental studies in supporting a sequential inhibitor-mutagen administration as a possible antiviral design. Given recent progress in the development of new inhibitors of arenavirus replication, our results suggest new options of ribavirin-based anti-arenavirus treatments.

Incongruent fitness landscapes, not tradeoffs, dominate the adaptation of vesicular stomatitis virus to novel host types

Host radiation refers to the ability of parasites to adapt to new environments and expand or change their niches. Adaptation to one specific environment may involve a loss in adaptation to a second environment. Thus, fitness costs may impose limits to niche expansion and constitute the cost of specialization. Several reports have addressed the cost of host radiation in vesicular stomatitis virus (VSV), but in some cases the experimental setup may have resulted in the overestimation of fitness costs. To clarify this issue, experiments were carried out in which a reference strain of VSV was allowed to adapt to HeLa, MDCK and BHK-21 cells, and to a regime of alternation between HeLa and Madin–Darby canine kidney (MDCK) cells. Measurement of viral fitness on each cell type showed that most virus populations behaved as generalists, and increased in fitness in all environments. Tradeoffs, where a fitness increase in one environment led to a fitness decrease in another environment, were rare. These results highlight the importance of using appropriate methods to measure fitness in evolved virus populations, and provide further support to a model of evolutionary dynamics in which costs due to incongruent landscapes provided by different environments are more common than tradeoffs.

Rapid adaptive amplification of preexisting variation in an RNA virus

The amount and nature of preexisting variation in a population of RNA viruses is an important determinant of the virus's ability to adapt rapidly to a changed environment. However, direct quantification of this preexisting variation may be cumbersome, because potentially beneficial alleles are typically rare, and isolation of a large number of subclones is required. Here, we propose a simpler method. We infer the initial population structure of vesicular stomatitis virus (VSV) by fitting a mathematical model of asexual evolution to an extensive set of measurements of VSV fitness dynamics under various conditions, including new and previously published data. The inferred variation of fitness in the initial population agrees very well with the results of direct experiments with subclone fitness quantification. From the same procedure, we also estimate the mean fitness effect of beneficial mutations (selection coefficient s), the percentage of sites in the genome that are under moderate positive or negative selection, and the percentage of sites where beneficial mutations may potentially occur. For VSV strain MARM U evolving in BHK-21 cells, the three parameters have values of 0.39, 9%, and 0.06%, respectively. The method can be generalized and applied easily to other rapidly evolving microbes, including both asexual microorganisms and those with recombination.

Viral Quasispecies: Dynamics, Interactions, and Pathogenesis*

Quasispecies theory is providing a solid, evolving conceptual framework for insights into virus population dynamics, adaptive potential, and response to lethal mutagenesis. The complexity of mutant spectra can influence disease progression and viral pathogenesis, as demonstrated using virus variants selected for increased replicative fidelity. Complementation and interference exerted among components of a viral quasispecies can either reinforce or limit the replicative capacity and disease potential of the ensemble. In particular, a progressive enrichment of a replicating mutant spectrum with interfering mutant genomes prompted by enhanced mutagenesis may be a key event in the sharp transition of virus populations into error catastrophe that leads to virus extinction. Fitness variations are influenced by the passage regimes to which viral populations are subjected, notably average fitness decreases upon repeated bottleneck events and fitness gains upon competitive optimization of large viral populations. Evolving viral quasispecies respond to selective constraints by replication of subpopulations of variant genomes that display higher fitness than the parental population in the presence of the selective constraint. This has been profusely documented with fitness effects of mutations associated with resistance of pathogenic viruses to antiviral agents. In particular, selection of HIV-1 mutants resistant to one or multiple antiretroviral inhibitors, and the compensatory effect of mutations in the same genome, offers a compendium of the molecular intricacies that a virus can exploit for its survival. This chapter reviews the basic principles of quasispecies dynamics as they can serve to explain the behavior of viruses.

Frequency-Dependent Selection in a Periodic Environment

We examine the action of natural selection in a periodically changing environment where two competing strains are specialists respectively for each environmental state. When the relative fitness of the strains is subject to a very general class of frequency-dependent selection, we show that coexistence rather than extinction is the likely outcome. This coexistence may be a stable periodic equilibrium, stable limit cycles of varying lengths, or be deterministically chaotic. Our model is applicable to the population dynamics commonly found in many types of viruses.

Dynamics of deletion genotypes in an experimental insect virus population

Defective viruses, that are deficient in certain essential genes, are maintained in the population by trans-complementation, exploiting the gene products of complete genotypes in co-infected cells. This process becomes prevalent only when cells are frequently infected by several virus particles, and only then will the fitness of defective viruses be subjected to frequency-dependent selection. Deletion variants that are not infectious per os are present in a multicapsid nucleopolyhedrovirus (SfMNPV, Baculoviridae) that infects the fall army worm, Spodoptera frugiperda. These variants enhance the pathogenicity and, therefore, the likelihood of transmission of the virus when co-infecting cells with complete genotypes, resulting in occlusion bodies (OBs) that may contain both genotypes co-occluded. Mixtures of complete (B) and defective (C) variants in ratios of 90% B+10% C, 50% B+50% C and 10% B+90% C were used to inoculate by injection S. frugiperda larvae. Viral OBs extracted from diseased insects were subjected to four or five successive rounds of per os infection. Following successive passages, genotype frequencies in all three experimental populations converged to a single equilibrium frequency comprising approximately 20% of deletion genotype C and approximately 80% of complete genotype B. This mirrors the relative proportions of deletion (22%) and complete (78%) genotypes observed in the wild-type SfMNPV population. The pathogenicity of experimental populations at the final passage was not significantly different from that of the wild-type isolate. In contrast, OBs of all genotype mixtures were significantly more pathogenic than OBs of genotype B alone. A population genetics model, in which virus populations were assigned linear frequency-dependent transmissibility values, was in remarkably close agreement to empirical data. Clearly, non-infectious deletion variants can profoundly affect the likelihood of transmission and the genetic structure and stability of virus populations.

Mechanisms of genetic robustness in RNA viruses

Two key features of RNA viruses are their compacted genomes and their high mutation rate. Accordingly, deleterious mutations are common and have an enormous impact on viral fitness. In their multicellular hosts, robustness can be achieved by genomic redundancy, including gene duplication, diploidy, alternative metabolic pathways and biochemical buffering mechanisms. However, here we review evidence suggesting that during RNA virus evolution, alternative robustness mechanisms may have been selected. After briefly describing how genetic robustness can be quantified, we discuss mechanisms of intrinsic robustness arising as consequences of RNA-genome architecture, replication peculiarities and quasi-species population dynamics. These intrinsic robustness mechanisms operate efficiently at the population level, despite the mutational sensitivity shown by individual genomes. Finally, we discuss the possibility that viruses might exploit cellular buffering mechanisms for their own benefit, producing a sort of extrinsic robustness.

Viruses as quasispecies: biological implications

During viral infections, the complex and dynamic distributions of variants, termed viral quasispecies, play a key role in the adaptability of viruses to changing environments and the fate of the population as a whole. Mutant spectra are continuously and avoidably generated during RNA genome replication, and they are not just a by-product of error-prone replication, devoid of biological relevance. On the contrary, current evidence indicates that mutant spectra contribute to viral pathogenesis, can modulate the expression of phenotypic traits by subpopulations of viruses, can include memory genomes that reflect the past evolutionary history of the viral lineage, and, furthermore, can participate in viral extinction through lethal mutagenesis. Also, mutant spectra are the target on which selection and random drift act to shape the long-term evolution of viruses. The biological relevance of mutant spectra is the central topic of this chapter.

Tradeoff between short-term and long-term adaptation in a changing environment

We investigate the competition dynamics of two microbial or viral strains that live in an environment that switches periodically between two states. One of the strains is adapted to the long-term environment, but pays a short-term cost, while the other is adapted to the short-term environment and pays a cost in the long term. We explore the tradeoff between these alternative strategies in extensive numerical simulations and present a simple analytic model that can predict the outcome of these competitions as a function of the mutation rate and the time scale of the environmental changes. Our model is relevant for arboviruses, which alternate between different host species on a regular basis.

Quasispecies theory in the context of population genetics

BACKGROUND

A number of recent papers have cast doubt on the applicability of the quasispecies concept to virus evolution, and have argued that population genetics is a more appropriate framework to describe virus evolution than quasispecies theory.

RESULTS

I review the pertinent literature, and demonstrate for a number of cases that the quasispecies concept is equivalent to the concept of mutation-selection balance developed in population genetics, and that there is no disagreement between the population genetics of haploid, asexually-replicating organisms and quasispecies theory.

CONCLUSION

Since quasispecies theory and mutation-selection balance are two sides of the same medal, the discussion about which is more appropriate to describe virus evolution is moot. In future work on virus evolution, we would do good to focus on the important questions, such as whether we can develop accurate, quantitative models of virus evolution, and to leave aside discussions about the relative merits of perfectly equivalent concepts.

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.

A segmented form of foot-and-mouth disease virus interferes with standard virus: a link between interference and competitive fitness

Serial passage of foot-and-mouth disease virus (FMDV) in BHK-21 cells at high multiplicity of infection resulted in dominance of particles containing defective RNAs that were infectious by complementation in the absence of standard viral RNA. In the present study, we show that the defective FMDV particles interfere with replication of the cognate standard virus. Coinfections of defective FMDV with standard FMDV mutants that differ up to 151-fold in relative fitness have documented that the degree of interference is higher for low fitness than for high fitness standard virus. These comparisons suggest a likely overlap between those mechanisms of intracellular competition that underlie viral interference and those expressed as fitness differences between two viruses when they coinfect the same cells. Interference may contribute to the selective pressures that help maintain dominance of segmented defective RNAs over the standard FMDV genome.

Suppression of viral infectivity through lethal defection

RNA viruses replicate with a very high error rate and give rise to heterogeneous, highly plastic populations able to adapt very rapidly to changing environments. Viral diseases are thus difficult to control because of the appearance of drug-resistant mutants, and it becomes essential to seek mechanisms able to force the extinction of the quasispecies before adaptation emerges. An alternative to the use of conventional drugs consists in increasing the replication error rate through the use of mutagens. Here, we report about persistent infections of lymphocytic choriomeningitis virus treated with fluorouracil, where a progressive debilitation of infectivity leading to eventual extinction occurs. The transition to extinction is accompanied by the production of large amounts of RNA, indicating that the replicative ability of the quasispecies is not strongly impaired by the mutagen. By means of experimental and theoretical approaches, we propose that a fraction of the RNA molecules synthesized can behave as a defective subpopulation able to drive the viable class extinct. Our results lead to the identification of two extinction pathways, one at high amounts of mutagen, where the quasispecies completely loses its ability to infect and replicate, and a second one, at lower amounts of mutagen, where replication continues while the infective class gets extinct because of the action of defectors. The results bear on a potential application of increased mutagenesis as an antiviral strategy in that low doses of a mutagenic agent may suffice to drive persistent virus to extinction.

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.

Density-dependent selection in vesicular stomatitis virus

We used vesicular stomatitis virus to test the effect of complementation on the relative fitness of a deleterious mutant, monoclonal antibody-resistant mutant (MARM) N, in competition with its wild-type ancestor. We carried out competitions of MARM N and wild-type populations at different multiplicities of infection (MOIs) and initial ratios of the wild type to the mutant and found that the fitness of MARM N relative to that of the wild type is very sensitive to changes in the MOI (i.e., the degree of complementation) but depends little, if at all, on the initial frequencies of MARM N and the wild type. Further, we developed a mathematical model under the assumption that during coinfection both viruses contribute to a common pool of protein products in the infected cell and that they both exploit this common pool equally. Under such conditions, the fitness of all virions that coinfect a cell is the average fitness in the absence of coinfection of that group of virions. In the absence of coinfection, complementation cannot take place and the relative fitness of each competitor is only determined by the selective value of its own products. We found good agreement between our experimental results and the model predictions, which suggests that the wild type and MARM N freely share all of their gene products under coinfection.

Negative effect of genetic bottlenecks on the adaptability of vesicular stomatitis virus

Muller's ratchet is a principle of evolutionary genetics describing mutant accumulation in populations that are repeatedly subjected to genetic bottleneck. The immediate effect of Muller's ratchet, overall loss of fitness, has been confirmed in several viral systems belonging to different groups. This report shows that in addition to fitness loss, genetic bottlenecks also have longer-term effects, namely changes in the capacity of viral populations to adapt. Thus, vesicular stomatitis virus strains with a history of genetic bottleneck have lower adaptability than strains maintained at relatively large population sizes. This lower adaptability is illustrated by their reduced ability to regain fitness and by their inability to outcompete wild-type populations in situations where the initial fitness of the bottlenecked mutant is the same or even higher than the initial fitness of the wild-type.