Exponential increases of RNA virus fitness during large population transmissions

Limits to selection on standing variation in an asexual population

We consider how a population of N haploid individuals responds to directional selection on standing variation, with no new variation from recombination or mutation. Individuals have trait values z1,…,zN, which are drawn from a distribution ψ; the fitness of individual i is proportional to [Formula: see text] . For illustration, we consider the Laplace and Gaussian distributions, which are parametrised only by the variance V0, and show that for large N, there is a scaling limit which depends on a single parameter NV0. When selection is weak relative to drift (NV0≪1), the variance decreases exponentially at rate 1/N, and the expected ultimate gain in log fitness (scaled by V0), is just NV0, which is the same as Robertson's (1960) prediction for a sexual population. In contrast, when selection is strong relative to drift (NV0≫1), the ultimate gain can be found by approximating the establishment of alleles by a branching process in which each allele competes independently with the population mean and the fittest allele to establish is certain to fix. Then, if the probability of survival to time t∼1/V0 of an allele with value z is P(z), with mean P¯, the winning allele is the fittest of NP¯ survivors drawn from a distribution ψP/P¯. The expected ultimate change is ∼2log(1.15NV0) for a Gaussian distribution, and ∼-12log0.36NV0-log-log0.36NV0 for a Laplace distribution. This approach also predicts the variability of the process, and its dynamics; we show that in the strong selection regime, the expected genetic variance decreases as ∼t-3 at large times. We discuss how these results may be related to selection on standing variation that is spread along a linear chromosome.

Differential immunophenotype of circulating monocytes from pregnant women in response to viral ligands

BACKGROUND

Viral infections during pregnancy can have deleterious effects on mothers and their offspring. Monocytes participate in the maternal host defense against invading viruses; however, whether pregnancy alters monocyte responses is still under investigation. Herein, we undertook a comprehensive in vitro study of peripheral monocytes to characterize the differences in phenotype and interferon release driven by viral ligands between pregnant and non-pregnant women.

METHODS

Peripheral blood was collected from third-trimester pregnant (n = 20) or non-pregnant (n = 20, controls) women. Peripheral blood mononuclear cells were isolated and exposed to R848 (TLR7/TLR8 agonist), Gardiquimod (TLR7 agonist), Poly(I:C) (HMW) VacciGrade™ (TLR3 agonist), Poly(I:C) (HMW) LyoVec™ (RIG-I/MDA-5 agonist), or ODN2216 (TLR9 agonist) for 24 h. Cells and supernatants were collected for monocyte phenotyping and immunoassays to detect specific interferons, respectively.

RESULTS

The proportions of classical (CD14^hiCD16^-), intermediate (CD14^hiCD16⁺), non-classical (CD14^loCD16⁺), and CD14^loCD16^- monocytes were differentially affected between pregnant and non-pregnant women in response to TLR3 stimulation. The proportions of pregnancy-derived monocytes expressing adhesion molecules (Basigin and PSGL-1) or the chemokine receptors CCR5 and CCR2 were diminished in response to TLR7/TLR8 stimulation, while the proportions of CCR5^- monocytes were increased. Such differences were found to be primarily driven by TLR8 signaling, rather than TLR7. Moreover, the proportions of monocytes expressing the chemokine receptor CXCR1 were increased during pregnancy in response to poly(I:C) stimulation through TLR3, but not RIG-I/MDA-5. By contrast, pregnancy-specific changes in the monocyte response to TLR9 stimulation were not observed. Notably, the soluble interferon response to viral stimulation by mononuclear cells was not diminished in pregnancy.

CONCLUSIONS

Our data provide insight into the differential responsiveness of pregnancy-derived monocytes to ssRNA and dsRNA, mainly driven by TLR8 and membrane-bound TLR3, which may help to explain the increased susceptibility of pregnant women to adverse outcomes resulting from viral infection as observed during recent and historic pandemics.

Mammarenavirus Genetic Diversity and Its Biological Implications

Members of the family Arenaviridae are classified into four genera: Antennavirus, Hartmanivirus, Mammarenavirus, and Reptarenavirus. Reptarenaviruses and hartmaniviruses infect (captive) snakes and have been shown to cause boid inclusion body disease (BIBD). Antennaviruses have genomes consisting of 3, rather than 2, segments, and were discovered in actinopterygian fish by next-generation sequencing but no biological isolate has been reported yet. The hosts of mammarenaviruses are mainly rodents and infections are generally asymptomatic. Current knowledge about the biology of reptarenaviruses, hartmaniviruses, and antennaviruses is very limited and their zoonotic potential is unknown. In contrast, some mammarenaviruses are associated with zoonotic events that pose a threat to human health. This review will focus on mammarenavirus genetic diversity and its biological implications. Some mammarenaviruses including lymphocytic choriomeningitis virus (LCMV) are excellent experimental model systems for the investigation of acute and persistent viral infections, whereas others including Lassa (LASV) and Junin (JUNV) viruses, the causative agents of Lassa fever (LF) and Argentine hemorrhagic fever (AHF), respectively, are important human pathogens. Mammarenaviruses were thought to have high degree of intra-and inter-species amino acid sequence identities, but recent evidence has revealed a high degree of mammarenavirus genetic diversity in the field. Moreover, closely related mammarenavirus can display dramatic phenotypic differences in vivo. These findings support a role of genetic variability in mammarenavirus adaptability and pathogenesis. Here, we will review the molecular biology of mammarenaviruses, phylogeny, and evolution, as well as the quasispecies dynamics of mammarenavirus populations and their biological implications.

Viral Fitness, Population Complexity, Host Interactions, and Resistance to Antiviral Agents

Fitness of viruses has become a standard parameter to quantify their adaptation to a biological environment. Fitness determinations for RNA viruses (and some highly variable DNA viruses) meet with several uncertainties. Of particular interest are those that arise from mutant spectrum complexity, absence of population equilibrium, and internal interactions among components of a mutant spectrum. Here, concepts, fitness measurements, limitations, and current views on experimental viral fitness landscapes are discussed. The effect of viral fitness on resistance to antiviral agents is covered in some detail since it constitutes a widespread problem in antiviral pharmacology, and a challenge for the design of effective antiviral treatments. Recent evidence with hepatitis C virus suggests the operation of mechanisms of antiviral resistance additional to the standard selection of drug-escape mutants. The possibility that high replicative fitness may be the driver of such alternative mechanisms is considered. New broad-spectrum antiviral designs that target viral fitness may curtail the impact of drug-escape mutants in treatment failures. We consider to what extent fitness-related concepts apply to coronaviruses and how they may affect strategies for COVID-19 prevention and treatment.

Virus Evolution on Fitness Landscapes

The landscape paradigm is revisited in the light of evolution in simple systems. A brief overview of different classes of fitness landscapes is followed by a more detailed discussion of the RNA model, which is currently the only evolutionary model that allows for a comprehensive molecular analysis of a fitness landscape. Neutral networks of genotypes are indispensable for the success of evolution. Important insights into the evolutionary mechanism are gained by considering the topology of sequence and shape spaces. The dynamic concept of molecular quasispecies is viewed in the light of the landscape paradigm. The distribution of fitness values in state space is mirrored by the population structures of mutant distributions. Two classes of thresholds for replication error or mutations are important: (i) the-conventional-genotypic error threshold, which separates ordered replication from random drift on neutral networks, and (ii) a phenotypic error threshold above which the molecular phenotype is lost. Empirical landscapes are reviewed and finally, the implications of the landscape concept for virus evolution are discussed.

A Two-Level, Intramutant Spectrum Haplotype Profile of Hepatitis C Virus Revealed by Self-Organized Maps

RNA viruses replicate as complex mutant spectra termed viral quasispecies. The frequency of each individual genome in a mutant spectrum depends on its rate of generation and its relative fitness in the replicating population ensemble. The advent of deep sequencing methodologies allows for the first-time quantification of haplotype abundances within mutant spectra. There is no information on the haplotype profile of the resident genomes and how the landscape evolves when a virus replicates in a controlled cell culture environment. Here, we report the construction of intramutant spectrum haplotype landscapes of three amplicons of the NS5A-NS5B coding region of hepatitis C virus (HCV). Two-dimensional (2D) neural networks were constructed for 44 related HCV populations derived from a common clonal ancestor that was passaged up to 210 times in human hepatoma Huh-7.5 cells in the absence of external selective pressures. The haplotype profiles consisted of an extended dense basal platform, from which a lower number of protruding higher peaks emerged. As HCV increased its adaptation to the cells, the number of haplotype peaks within each mutant spectrum expanded, and their distribution shifted in the 2D network. The results show that extensive HCV replication in a monotonous cell culture environment does not limit HCV exploration of sequence space through haplotype peak movements. The landscapes reflect dynamic variation in the intramutant spectrum haplotype profile and may serve as a reference to interpret the modifications produced by external selective pressures or to compare with the landscapes of mutant spectra in complex in vivo environments. IMPORTANCE The study provides for the first time the haplotype profile and its variation in the course of virus adaptation to a cell culture environment in the absence of external selective constraints. The deep sequencing-based self-organized maps document a two-layer haplotype distribution with an ample basal platform and a lower number of protruding peaks. The results suggest an inferred intramutant spectrum fitness landscape structure that offers potential benefits for virus resilience to mutational inputs.

Tradeoffs for a viral mutant with enhanced replication speed

RNA viruses exist as genetically heterogeneous populations due to high mutation rates, and many of these mutations reduce fitness and/or replication speed. However, it is unknown whether mutations can increase replication speed of a virus already well adapted to replication in cultured cells. By sequentially passaging coxsackievirus B3 in cultured cells and collecting the very earliest progeny, we selected for increased replication speed. We found that a single mutation in a viral capsid protein, VP1-F106L, was sufficient for the fast-replication phenotype. Characterization of this mutant revealed quicker genome release during entry compared to wild-type virus, highlighting a previously unappreciated infection barrier. However, this mutation also reduced capsid stability in vitro and reduced replication and pathogenesis in mice. These results reveal a tradeoff between overall replication speed and fitness. Importantly, this approach-selecting for the earliest viral progeny-could be applied to a variety of viral systems and has the potential to reveal unanticipated inefficiencies in viral replication cycles.

From RNA World to SARS-CoV-2: The Edited Story of RNA Viral Evolution

The current SARS-CoV-2 pandemic underscores the importance of understanding the evolution of RNA genomes. While RNA is subject to the formation of similar lesions as DNA, the evolutionary and physiological impacts RNA lesions have on viral genomes are yet to be characterized. Lesions that may drive the evolution of RNA genomes can induce breaks that are repaired by recombination or can cause base substitution mutagenesis, also known as base editing. Over the past decade or so, base editing mutagenesis of DNA genomes has been subject to many studies, revealing that exposure of ssDNA is subject to hypermutation that is involved in the etiology of cancer. However, base editing of RNA genomes has not been studied to the same extent. Recently hypermutation of single-stranded RNA viral genomes have also been documented though its role in evolution and population dynamics. Here, we will summarize the current knowledge of key mechanisms and causes of RNA genome instability covering areas from the RNA world theory to the SARS-CoV-2 pandemic of today. We will also highlight the key questions that remain as it pertains to RNA genome instability, mutations accumulation, and experimental strategies for addressing these questions.

Animal and human RNA viruses: genetic variability and ability to overcome vaccines

RNA viruses, in general, exhibit high mutation rates; this is mainly due to the low fidelity displayed by the RNA-dependent polymerases required for their replication that lack the proofreading machinery to correct misincorporated nucleotides and produce high mutation rates. This lack of replication fidelity, together with the fact that RNA viruses can undergo spontaneous mutations, results in genetic variants displaying different viral morphogenesis, as well as variation on their surface glycoproteins that affect viral antigenicity. This diverse viral population, routinely containing a variety of mutants, is known as a viral ‘quasispecies’. The mutability of their virions allows for fast evolution of RNA viruses that develop antiviral resistance and overcome vaccines much more rapidly than DNA viruses. This also translates into the fact that pathogenic RNA viruses, that cause many diseases and deaths in humans, represent the major viral group involved in zoonotic disease transmission, and are responsible for worldwide pandemics.

A new implication of quasispecies dynamics: Broad virus diversification in absence of external perturbations

RNA genetic elements include many important animal and plant pathogens. They share high mutability, a trait that has multiple implications for the interactions with their host organisms. Here we review evidence of a new adaptive feature of RNA viruses that we term "broadly diversifying selection". It constitutes a new type of positive selection without participation of any external selective agent, and which is built upon a progressive increase of the number of different genomes that dominate the population. The evidence was provided by analyses of mutant spectrum composition of two important viral pathogens, foot-and-mouth disease virus (FMDV) and hepatitis C virus (HCV) after prolonged replication in their respective cell culture environment. Despite being fueled by mutations that arise randomly and in absence of an external guiding selective force, this type of selection prepares the viral population for a response to selective forces still to occur. Since current evidence suggests that broadly diversifying selection is favored by elevated mutation rates and population sizes, it may constitute a more general behavior, relevant also to the adaptive dynamics of microbial populations and cancer cells.

Viral fitness: history and relevance for viral pathogenesis and antiviral interventions

The quasispecies dynamics of viral populations (continuous generation of variant genomes and competition among them) has as one of its frequent consequences variations in overall multiplication capacity, a major component of viral fitness. This parameter has multiple implications for viral pathogenesis and viral disease control, some of them unveiled thanks to deep sequencing of viral populations. Darwinian fitness is an old concept whose quantification dates back to the early developments of population genetics. It was later applied to viruses (mainly to RNA viruses) to quantify relative multiplication capacities of individual mutant clones or complex populations. The present article reviews the fitness concept and its relevance for the understanding of the adaptive dynamics of viruses in constant and changing environments. Many studies have addressed the fitness cost of escape mutations (to antibodies, cytotoxic T cells or inhibitors) as an influence on the efficacy of antiviral interventions. Here, we summarize the evidence that the basal fitness level can be a determinant of inhibitor resistance.

Estimating viral bottleneck sizes for FMDV transmission within and between hosts and implications for the rate of viral evolution

RNA viruses exist as populations of closely related genomes, characterized by a high diversity of low-frequency variants. As viral genomes from one population disperse to establish new sites of replication, the fate of these low-frequency variants depends to a large extent on the size of the founding population. Focusing on foot-and-mouth disease virus (FMDV) we conjecture that variants are more likely to be transmitted through wide bottlenecks, but more likely to approach fixation in new populations following narrow bottlenecks; therefore, the longer-term rate of accumulation of 'nearly neutral' variants at high frequencies is likely to be inversely related to the bottleneck size. We examine this conjecture in vivo by estimating bottleneck sizes relating 'parent' and 'daughter' populations observed at different scales ranging from within host to between host (within the same herd, and in different herds) using a previously established method. Within hosts, we find bottleneck sizes to range from 5 to 20 viral genomes between populations transmitted from the pharynx to the serum, and from 4 to 54 between serum and lesion populations. Between hosts, we find bottleneck sizes to range from 2 to 39, suggesting inter-host bottlenecks are of a similar size to intra-host bottlenecks. We establish a statistically significant negative relationship between the probability of genomic consensus level change and bottleneck size, and present a simple sampling model that captures this empirical relationship. We also present a novel in vitro experiment to investigate the impact of bottleneck size on the frequency of mutations within FMDV populations, demonstrate that variant frequency in a population increases more rapidly during small population passages, and provide evidence for positive selection during the passage of large populations.

Rare haplotype load as marker for lethal mutagenesis

RNA viruses replicate with a template-copying fidelity, which lies close to an extinction threshold. Increases of mutation rate by nucleotide analogues can drive viruses towards extinction. This transition is the basis of an antiviral strategy termed lethal mutagenesis. We have introduced a new diversity index, the rare haplotype load (RHL), to describe NS5B (polymerase) mutant spectra of hepatitis C virus (HCV) populations passaged in absence or presence of the mutagenic agents favipiravir or ribavirin. The increase in RHL is more prominent in mutant spectra whose expansions were due to nucleotide analogues than to multiple passages in absence of mutagens. Statistical tests for paired mutagenized versus non-mutagenized samples with 14 diversity indices show that RHL provides consistently the highest standardized effect of mutagenic treatment difference for ribavirin and favipiravir. The results indicate that the enrichment of viral quasispecies in very low frequency minority genomes can serve as a robust marker for lethal mutagenesis. The diagnostic value of RHL from deep sequencing data is relevant to experimental studies on enhanced mutagenesis of viruses, and to pharmacological evaluations of inhibitors suspected to have a mutagenic activity.

Viral Fitness Correlates with the Magnitude and Direction of the Perturbation Induced in the Host's Transcriptome: The Tobacco Etch Potyvirus-Tobacco Case Study

Determining the fitness of viral genotypes has become a standard practice in virology as it is essential to evaluate their evolutionary potential. Darwinian fitness, defined as the advantage of a given genotype with respect to a reference one, is a complex property that captures, in a single figure, differences in performance at every stage of viral infection. To what extent does viral fitness result from specific molecular interactions with host factors and regulatory networks during infection? Can we identify host genes in functional classes whose expression depends on viral fitness? Here, we compared the transcriptomes of tobacco plants infected with seven genotypes of tobacco etch potyvirus that differ in fitness. We found that the larger the fitness differences among genotypes, the more dissimilar the transcriptomic profiles are. Consistently, two different mutations, one in the viral RNA polymerase and another in the viral suppressor of RNA silencing, resulted in significantly similar gene expression profiles. Moreover, we identified host genes whose expression showed a significant correlation, positive or negative, with the virus' fitness. Differentially expressed genes which were positively correlated with viral fitness activate hormone- and RNA silencing-mediated pathways of plant defense. In contrast, those that were negatively correlated with fitness affect metabolism, reducing growth, and development. Overall, these results reveal the high information content of viral fitness and suggest its potential use to predict differences in genomic profiles of infected hosts.

Quasispecies and virus

Quasispecies theory has been instrumental in the understanding of RNA virus population dynamics because it considered for the first time mutation as an integral part of the replication process. The key influences of quasispecies theory on experimental virology have been: (1) to disclose the mutant spectrum nature of viral populations and to evaluate its consequences; (2) to unveil collective properties of genome ensembles that can render a mutant spectrum a unit of selection; and (3) to identify new vulnerability points of pathogenic RNA viruses on three fronts: the need to apply multiple selective constraints (in the form of drug combinations) to minimize selection of treatment-escape variants, to translate the error threshold concept into antiviral designs, and to construct attenuated vaccine viruses through alterations of viral polymerase copying fidelity or through displacements of viral genomes towards unfavorable regions of sequence space. These three major influences on the understanding of viral pathogens preceded extensions of quasispecies to non-viral systems such as bacterial and tumor cell collectivities and prions. These developments are summarized here.

Bacteria Facilitate Enteric Virus Co-infection of Mammalian Cells and Promote Genetic Recombination

RNA viruses exist in genetically diverse populations due to high levels of mutations, many of which reduce viral fitness. Interestingly, intestinal bacteria can promote infection of several mammalian enteric RNA viruses, but the mechanisms and consequences are unclear. We screened a panel of 41 bacterial strains as a platform to determine how different bacteria impact infection of poliovirus, a model enteric virus. Most bacterial strains, including those extracted from cecal contents of mice, bound poliovirus, with each bacterium binding multiple virions. Certain bacterial strains increased viral co-infection of mammalian cells even at a low virus-to-host cell ratio. Bacteria-mediated viral co-infection correlated with bacterial adherence to cells. Importantly, bacterial strains that induced viral co-infection facilitated genetic recombination between two different viruses, thereby removing deleterious mutations and restoring viral fitness. Thus, bacteria-virus interactions may increase viral fitness through viral recombination at initial sites of infection, potentially limiting abortive infections.

Differences in adaptive dynamics determine the success of virus variants that propagate together

Virus fitness is a complex parameter that results from the interaction of virus-specific characters (e.g. intracellular growth rate, adsorption rate, virion extracellular stability, and tolerance to mutations) with others that depend on the underlying fitness landscape and the internal structure of the whole population. Individual mutants usually have lower fitness values than the complex population from which they come from. When they are propagated and allowed to attain large population sizes for a sufficiently long time, they approach mutation-selection equilibrium with the concomitant fitness gains. The optimization process follows dynamics that vary among viruses, likely due to differences in any of the parameters that determine fitness values. As a consequence, when different mutants spread together, the number of generations experienced by each of them prior to co-propagation may determine its particular fate. In this work we attempt a clarification of the effect of different levels of population diversity in the outcome of competition dynamics. To this end, we analyze the behavior of two mutants of the RNA bacteriophage Qβ that co-propagate with the wild-type virus. When both competitor viruses are clonal, the mutants rapidly outcompete the wild type. However, the outcome in competitions performed with partially optimized virus populations depends on the distance of the competitors to their clonal origin. We also implement a theoretical population dynamics model that describes the evolution of a heterogeneous population of individuals, each characterized by a fitness value, subjected to subsequent cycles of replication and mutation. The experimental results are explained in the framework of our theoretical model under two non-excluding, likely complementary assumptions: (1) The relative advantage of both competitors changes as populations approach mutation-selection equilibrium, as a consequence of differences in their growth rates and (2) one of the competitors is more robust to mutations than the other. The main conclusion is that the nearness of an RNA virus population to mutation-selection equilibrium is a key factor determining the fate of particular mutants arising during replication.

Dynamical trade-offs arise from antagonistic coevolution and decrease intraspecific diversity

Trade-offs play an important role in evolution. Without trade-offs, evolution would maximize fitness of all traits leading to a "master of all traits". The shape of trade-offs has been shown to determine evolutionary trajectories and is often assumed to be static and independent of the actual evolutionary process. Here we propose that coevolution leads to a dynamical trade-off. We test this hypothesis in a microbial predator-prey system and show that the bacterial growth-defense trade-off changes from concave to convex, i.e., defense is effective and cheap initially, but gets costly when predators coevolve. We further explore the impact of such dynamical trade-offs by a novel mathematical model incorporating de novo mutations for both species. Predator and prey populations diversify rapidly leading to higher prey diversity when the trade-off is concave (cheap). Coevolution results in more convex (costly) trade-offs and lower prey diversity compared to the scenario where only the prey evolves.

Limiting fitness distributions in evolutionary dynamics

Natural selection works on variation in fitness, but how should we measure "variation" to predict the rate of future evolution? Fisher's fundamental theorem of natural selection provides the short-run answer: the instantaneous rate of growth of a population's mean fitness is its variance in fitness. This identity captures an important feature of the evolutionary process, but, because it does not specify how the variance itself evolves in time, it cannot be used to predict evolutionary dynamics in the long run. In this paper we reconsider the problem of computing evolutionary trajectories from limited statistical information. We identify the feature of fitness distributions which controls their late-time evolution: their (suitably defined) tail indices. We show that the location, scale and shape of the fitness distribution can be predicted far into the future from the measurement of this tail index at some initial time. Unlike the "fitness waves" studied in the literature, this pattern encompasses both positive and negative selection and is not restricted to rapidly adapting populations. Our results are well supported by numerical simulations, both from the Wright-Fisher model and from a less structured genetic algorithm.

Antigenic diversification is correlated with increased thermostability in a mammalian virus

The theory of plastogenetic congruence posits that ultimately, the pressure to maintain function in the face of biomolecular destabilization produces robustness. As temperature goes up so does destabilization. Thus, genetic robustness, defined as phenotypic constancy despite mutation, should correlate with survival during thermal challenge. We tested this hypothesis using vesicular stomatitis virus (VSV). We produced two sets of evolved strains after selection for higher thermostability by either preincubation at 37°C or by incubation at 40°C during infection. These VSV populations became more thermostable and also more fit in the absence of thermal selection, demonstrating an absence of tradeoffs. Eleven out of 12 evolved populations had a fixed, nonsynonymous substitution in the nucleocapsid (N) open reading frame. There was a partial correlation between thermostability and mutational robustness that was observed when the former was measured at 42°C, but not at 37°C. These results are consistent with our earlier work and suggest that the relationship between robustness and thermostability is complex. Surprisingly, many of the thermostable strains also showed increased resistance to monoclonal antibody and polyclonal sera, including sera from natural hosts. These data suggest that evolved thermostability may lead to antigenic diversification and an increased ability to escape immune surveillance in febrile hosts, and potentially to an improved robustness. These relationships have important implications not only in terms of viral pathogenesis, but also for the development of vaccine vectors and oncolytic agents.

Genetic variation in fitness within a clonal population of a plant RNA virus

A long-standing observation in evolutionary virology is that RNA virus populations are highly polymorphic, composed by a mixture of genotypes whose abundances in the population depend on complex interaction between fitness differences, mutational coupling and genetic drift. It was shown long ago, though in cell cultures, that most of these genotypes had lower fitness than the population they belong, an observation that explained why single-virion passages turned on Muller’s ratchet while very large population passages resulted in fitness increases in novel environments. Here we report the results of an experiment specifically designed to evaluate in vivo the fitness differences among the subclonal components of a clonal population of the plant RNA virus tobacco etch potyvirus (TEV). Over 100 individual biological subclones from a TEV clonal population well adapted to the natural tobacco host were obtained by infectivity assays on a local lesion host. The replicative fitness of these subclones was then evaluated during infection of tobacco relative to the fitness of large random samples taken from the starting clonal population. Fitness was evaluated at increasing number of days post-inoculation. We found that at early days, the average fitness of subclones was significantly lower than the fitness of the clonal population, thus confirming previous observations that most subclones contained deleterious mutations. However, as the number of days of viral replication increases, population size expands exponentially, more beneficial and compensatory mutations are produced, and selection becomes more effective in optimizing fitness, the differences between subclones and the population disappeared.

Arenavirus Quasispecies and Their Biological Implications

The family Arenaviridae currently comprises over 20 viral species, each of them associated with a main rodent species as the natural reservoir and in one case possibly phyllostomid bats. Moreover, recent findings have documented a divergent group of arenaviruses in captive alethinophidian snakes. Human infections occur through mucosal exposure to aerosols or by direct contact of abraded skin with infectious materials. Arenaviruses merit interest both as highly tractable experimental model systems to study acute and persistent infections and as clinically important human pathogens including Lassa (LASV) and Junin (JUNV) viruses, the causative agents of Lassa and Argentine hemorrhagic fevers (AHFs), respectively, for which there are no FDA-licensed vaccines, and current therapy is limited to an off-label use of ribavirin (Rib) that has significant limitations. Arenaviruses are enveloped viruses with a bi-segmented negative strand (NS) RNA genome. Each genome segment, L (ca 7.3 kb) and S (ca 3.5 kb), uses an ambisense coding strategy to direct the synthesis of two polypeptides in opposite orientation, separated by a noncoding intergenic region (IGR). The S genomic RNA encodes the virus nucleoprotein (NP) and the precursor (GPC) of the virus surface glycoprotein that mediates virus receptor recognition and cell entry via endocytosis. The L genome RNA encodes the viral RNA-dependent RNA polymerase (RdRp, or L polymerase) and the small (ca 11 kDa) RING finger protein Z that has functions of a bona fide matrix protein including directing virus budding. Arenaviruses were thought to be relatively stable genetically with intra- and interspecies amino acid sequence identities of 90-95 % and 44-63 %, respectively. However, recent evidence has documented extensive arenavirus genetic variability in the field. Moreover, dramatic phenotypic differences have been documented among closely related LCMV isolates. These data provide strong evidence of viral quasispecies involvement in arenavirus adaptability and pathogenesis. Here, we will review several aspects of the molecular biology of arenaviruses, phylogeny and evolution, and quasispecies dynamics of arenavirus populations for a better understanding of arenavirus pathogenesis, as well as for the development of novel antiviral strategies to combat arenavirus infections.

Statistical Mechanics and Thermodynamics of Viral Evolution

This paper uses methods drawn from physics to study the life cycle of viruses. The paper analyzes a model of viral infection and evolution using the "grand canonical ensemble" and formalisms from statistical mechanics and thermodynamics. Using this approach we enumerate all possible genetic states of a model virus and host as a function of two independent pressures-immune response and system temperature. We prove the system has a real thermodynamic temperature, and discover a new phase transition between a positive temperature regime of normal replication and a negative temperature "disordered" phase of the virus. We distinguish this from previous observations of a phase transition that arises as a function of mutation rate. From an evolutionary biology point of view, at steady state the viruses naturally evolve to distinct quasispecies. This paper also reveals a universal relationship that relates the order parameter (as a measure of mutational robustness) to evolvability in agreement with recent experimental and theoretical work. Given that real viruses have finite length RNA segments that encode proteins which determine virus fitness, the approach used here could be refined to apply to real biological systems, perhaps providing insight into immune escape, the emergence of novel pathogens and other results of viral evolution.

Influence of mutation and recombination on HIV-1 in vitro fitness recovery

The understanding of the evolutionary processes underlying HIV-1 fitness recovery is fundamental for HIV-1 pathogenesis, antiretroviral treatment and vaccine design. It is known that HIV-1 can present very high mutation and recombination rates, however the specific contribution of these evolutionary forces in the "in vitro" viral fitness recovery has not been simultaneously quantified. To this aim, we analyzed substitution, recombination and molecular adaptation rates in a variety of HIV-1 biological clones derived from a viral isolate after severe population bottlenecks and a number of large population cell culture passages. These clones presented an overall but uneven fitness gain, mean of 3-fold, respect to the initial passage values. We found a significant relationship between the fitness increase and the appearance and fixation of mutations. In addition, these fixed mutations presented molecular signatures of positive selection through the accumulation of non-synonymous substitutions. Interestingly, viral recombination correlated with fitness recovery in most of studied viral quasispecies. The genetic diversity generated by these evolutionary processes was positively correlated with the viral fitness. We conclude that HIV-1 fitness recovery can be derived from the genetic heterogeneity generated through both mutation and recombination, and under diversifying molecular adaptation. The findings also suggest nonrandom evolutionary pathways for in vitro fitness recovery.

Sustainable deployment of QTLs conferring quantitative resistance to crops: first lessons from a stochastic model

Quantitative plant disease resistance is believed to be more durable than qualitative resistance, since it exerts less selective pressure on the pathogens. However, the process of progressive pathogen adaptation to quantitative resistance is poorly understood, which makes it difficult to predict its durability or to derive principles for its sustainable deployment. Here, we study the dynamics of pathogen adaptation in response to quantitative plant resistance affecting pathogen reproduction rate and its colonizing capacity. We developed a stochastic model for the continuous evolution of a pathogen population within a quantitatively resistant host. We assumed that pathogen can adapt to a host by the progressive restoration of reproduction rate or of colonizing capacity, or of both. Our model suggests that a combination of quantitative trait loci (QTLs) affecting distinct pathogen traits was more durable if the evolution of repressed traits was antagonistic. Otherwise, quantitative resistance that depressed only pathogen reproduction was more durable. In order to decelerate the progressive pathogen adaptation, QTLs that decrease the pathogen's maximum capacity to colonize must be combined with QTLs that decrease the spore production per lesion or the infection efficiency or that increase the latent period. Our theoretical framework can help breeders to develop principles for sustainable deployment of QTLs.

Increased replicative fitness can lead to decreased drug sensitivity of hepatitis C virus

Passage of hepatitis C virus (HCV) in human hepatoma cells resulted in populations that displayed partial resistance to alpha interferon (IFN-α), telaprevir, daclatasvir, cyclosporine, and ribavirin, despite no prior exposure to these drugs. Mutant spectrum analyses and kinetics of virus production in the absence and presence of drugs indicate that resistance is not due to the presence of drug resistance mutations in the mutant spectrum of the initial or passaged populations but to increased replicative fitness acquired during passage. Fitness increases did not alter host factors that lead to shutoff of general host cell protein synthesis and preferential translation of HCV RNA. The results imply that viral replicative fitness is a mechanism of multidrug resistance in HCV. Importance: Viral drug resistance is usually attributed to the presence of amino acid substitutions in the protein targeted by the drug. In the present study with HCV, we show that high viral replicative fitness can confer a general drug resistance phenotype to the virus. The results exclude the possibility that genomes with drug resistance mutations are responsible for the observed phenotype. The fact that replicative fitness can be a determinant of multidrug resistance may explain why the virus is less sensitive to drug treatments in prolonged chronic HCV infections that favor increases in replicative fitness.

Neutral null models for diversity in serial transfer evolution experiments

Evolution experiments with microorganisms coupled with genome-wide sequencing now allow for the systematic study of population genetic processes under a wide range of conditions. In learning about these processes in natural, sexual populations, neutral models that describe the behavior of diversity and divergence summaries have played a pivotal role. It is therefore natural to ask whether neutral models, suitably modified, could be useful in the context of evolution experiments. Here, we introduce coalescent models for polymorphism and divergence under the most common experimental evolution assay, a serial transfer experiment. This relatively simple setting allows us to address several issues that could affect diversity patterns in evolution experiments, whether selection is operating or not: the transient behavior of neutral polymorphism in an experiment beginning from a single clone, the effects of randomness in the timing of cell division and noisiness in population size in the dilution stage. In our analyses and discussion, we emphasize the implications for experiments aimed at measuring diversity patterns and making inferences about population genetic processes based on these measurements.

High genetic diversity and adaptive potential of two simian hemorrhagic fever viruses in a wild primate population

Key biological properties such as high genetic diversity and high evolutionary rate enhance the potential of certain RNA viruses to adapt and emerge. Identifying viruses with these properties in their natural hosts could dramatically improve disease forecasting and surveillance. Recently, we discovered two novel members of the viral family Arteriviridae: simian hemorrhagic fever virus (SHFV)-krc1 and SHFV-krc2, infecting a single wild red colobus (Procolobus rufomitratus tephrosceles) in Kibale National Park, Uganda. Nearly nothing is known about the biological properties of SHFVs in nature, although the SHFV type strain, SHFV-LVR, has caused devastating outbreaks of viral hemorrhagic fever in captive macaques. Here we detected SHFV-krc1 and SHFV-krc2 in 40% and 47% of 60 wild red colobus tested, respectively. We found viral loads in excess of 10(6)-10(7) RNA copies per milliliter of blood plasma for each of these viruses. SHFV-krc1 and SHFV-krc2 also showed high genetic diversity at both the inter- and intra-host levels. Analyses of synonymous and non-synonymous nucleotide diversity across viral genomes revealed patterns suggestive of positive selection in SHFV open reading frames (ORF) 5 (SHFV-krc2 only) and 7 (SHFV-krc1 and SHFV-krc2). Thus, these viruses share several important properties with some of the most rapidly evolving, emergent RNA viruses.

Experimental selection reveals a trade-off between fecundity and lifespan in the coliphage Qß

Understanding virus evolution is key for improving ways to counteract virus-borne diseases. Results from comparative analyses have previously suggested a trade-off between fecundity and lifespan for viruses that infect the bacterium Escherichia coli (i.e. for coliphages), which, if confirmed, would define a particular constraint on the evolution of virus fecundity. Here, the occurrence of such a trade-off is investigated through a selection experiment using the coliphage Qß. Selection was applied for increased fecundity in three independent wild-type Qß populations, and the ability of the virions to remain viable outside the host was determined. The Qß life-history traits involved in the evolution of fecundity and the genetic changes associated with this evolution were also investigated. The results reveal that short-term evolution of increased fecundity in Qß was associated with decreased viability of phage virions. This trade-off apparently arose because fecundity increased at the expense of reducing the amount of resources (mainly time) invested per produced virion. Thus, the results also indicate that Qß fecundity may be enhanced through increases in the rates of adsorption to the host and progeny production. Finally, genomic sequencing of the evolved populations pinpointed sequences likely to be involved in the evolution of Qß fecundity.

Genetic diversity among pandemic 2009 influenza viruses isolated from a transmission chain

Background

Influenza viruses such as swine-origin influenza A(H1N1) virus (A(H1N1)pdm09) generate genetic diversity due to the high error rate of their RNA polymerase, often resulting in mixed genotype populations (intra-host variants) within a single infection. This variation helps influenza to rapidly respond to selection pressures, such as those imposed by the immunological host response and antiviral therapy. We have applied deep sequencing to characterize influenza intra-host variation in a transmission chain consisting of three cases due to oseltamivir-sensitive viruses, and one derived oseltamivir-resistant case.

Methods

Following detection of the A(H1N1)pdm09 infections, we deep-sequenced the complete NA gene from two of the oseltamivir-sensitive virus-infected cases, and all eight gene segments of the viruses causing the remaining two cases.

Results

No evidence for the resistance-causing mutation (resulting in NA H275Y substitution) was observed in the oseltamivir-sensitive cases. Furthermore, deep sequencing revealed a subpopulation of oseltamivir-sensitive viruses in the case carrying resistant viruses. We detected higher levels of intra-host variation in the case carrying oseltamivir-resistant viruses than in those infected with oseltamivir-sensitive viruses.

Conclusions

Oseltamivir-resistance was only detected after prophylaxis with oseltamivir, suggesting that the mutation was selected for as a result of antiviral intervention. The persisting oseltamivir-sensitive virus population in the case carrying resistant viruses suggests either that a small proportion survive the treatment, or that the oseltamivir-sensitive virus rapidly re-establishes itself in the virus population after the bottleneck. Moreover, the increased intra-host variation in the oseltamivir-resistant case is consistent with the hypothesis that the population diversity of a RNA virus can increase rapidly following a population bottleneck.

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 importance of intergenic recombination in norovirus GII.3 evolution

Norovirus genotype II.3 (GII.3) strains are a major cause of sporadic gastroenteritis. Intergenic recombination between the capsid and RNA-dependent RNA polymerase (RdRp) genes is common and results in the acquisition of an alternative RdRp genotype. This study aimed to explore the evolution of the GII.3 capsid gene, focusing on the influence of intergenic recombination. The capsid genes from six GII.3 norovirus strains, collected from Australian children between 2001 and 2010, were sequenced and aligned with 66 GII.3 capsid sequences from GenBank, spanning 1975 to 2010. The GII.3 capsid gene evolved at a rate of 4.16 × 10(-3) to 6.97 × 10(-3) nucleotide substitutions/site/year from 1975 to 2010 and clustered into five temporally sequential lineages. Clustering of the GII.3 capsid gene sequences was associated with intergenic recombination and switches between RdRp genotypes GII.3, GII.a, GII.b, GII.12, and an undefined ancestral RdRp. Comparison of the substitution rate of the GII.3 and GII.b RdRps suggested that RdRp switching allows a higher evolutionary rate, leading to increased genetic diversity and adaptability. Alignment of GII.3 capsid sequences revealed 36 lineage-specific conserved amino acid substitutions, four of which were under positive selection. Many conserved substitutions were within predicted antibody binding regions and close to host attachment factor binding sites. In conclusion, evolution of GII.3 noroviruses was primarily driven by intergenic recombination. The acquisition of new RdRps may lead to a faster mutation rate and increased genetic diversity, improving overall GII.3 fitness.

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.

A single amino acid change resulting in loss of fluorescence of eGFP in a viral fusion protein confers fitness and growth advantage to the recombinant vesicular stomatitis virus

Using a recombinant vesicular stomatitis virus encoding eGFP fused in-frame with an essential viral replication protein, the phosphoprotein P, we show that during passage in culture, the virus mutates the nucleotide C289 within eGFP of the fusion protein PeGFP to A or T, resulting in R97S/C amino acid substitution and loss of fluorescence. The resultant non-fluorescent virus exhibits increased fitness and growth advantage over its fluorescent counterpart. The growth advantage of the non-fluorescent virus appears to be due to increased transcription and replication activities of the PeGFP protein carrying the R97S/C substitution. Further, our results show that the R97S/C mutation occurs prior to accumulation of mutations that can result in loss of expression of the gene inserted at the G-L gene junction. These results suggest that fitness gain is more important for the recombinant virus than elimination of expression of the heterologous gene.

Viral quasispecies evolution

Evolution of RNA viruses occurs through disequilibria of collections of closely related mutant spectra or mutant clouds termed viral quasispecies. Here we review the origin of the quasispecies concept and some biological implications of quasispecies dynamics. Two main aspects are addressed: (i) mutant clouds as reservoirs of phenotypic variants for virus adaptability and (ii) the internal interactions that are established within mutant spectra that render a virus ensemble the unit of selection. The understanding of viruses as quasispecies has led to new antiviral designs, such as lethal mutagenesis, whose aim is to drive viruses toward low fitness values with limited chances of fitness recovery. The impact of quasispecies for three salient human pathogens, human immunodeficiency virus and the hepatitis B and C viruses, is reviewed, with emphasis on antiviral treatment strategies. Finally, extensions of quasispecies to nonviral systems are briefly mentioned to emphasize the broad applicability of quasispecies theory.

Kinetics of virus production from single cells

The production of virus by infected cells is an essential process for the spread and persistence of viral diseases, the effectiveness of live-viral vaccines, and the manufacture of viruses for diverse applications. Yet despite its importance, methods to precisely measure virus production from cells are lacking. Most methods test infected-cell populations, masking how individual cells behave. Here we measured the kinetics of virus production from single cells. We combined simple steps of liquid-phase infection, serial dilution, centrifugation, and harvesting, without specialized equipment, to track the production of virus particles from BHK cells infected with vesicular stomatitis virus. Remarkably, cell-to-cell differences in latent times to virus release were within a factor of two, while production rates and virus yields spanned over 300-fold, highlighting an extreme diversity in virus production for cells from the same population. These findings have fundamental and technological implications for health and disease.

Specific and nonspecific host adaptation during arboviral experimental evolution

During the past decade or so, there has been a substantial body of work to dissect arboviral evolution and to develop models of adaptation during host switching. Regardless of what species serve as host or vectors, and of the geographic distribution and the mechanisms of replication, arboviruses tend to have slow evolutionary rates in nature. The hypothesis that this is the result of replication in the disparate environments provided by host and vector did not receive solid experimental support in any of the many viral species tested. Instead, it seems that from the virus's point of view, either the two environments are sufficiently similar or one of the environments so dominates viral evolution that there is tolerance for suboptimal adaptation to the other environment. Replication in alternating environments has an unexpected cost in that there is decreased genetic variance that translates into a compromised adaptability for bypassed environments. Arboviruses under strong and continuous positive selection may have unusual patterns of genomic changes, with few or no mutations accumulated in the consensus sequence or with dN/dS values typically consistent with random drift in DNA-based organisms.

Stickbreaking: a novel fitness landscape model that harbors epistasis and is consistent with commonly observed patterns of adaptive evolution

In relating genotypes to fitness, models of adaptation need to both be computationally tractable and qualitatively match observed data. One reason that tractability is not a trivial problem comes from a combinatoric problem whereby no matter in what order a set of mutations occurs, it must yield the same fitness. We refer to this as the bookkeeping problem. Because of their commutative property, the simple additive and multiplicative models naturally solve the bookkeeping problem. However, the fitness trajectories and epistatic patterns they predict are inconsistent with the patterns commonly observed in experimental evolution. This motivates us to propose a new and equally simple model that we call stickbreaking. Under the stickbreaking model, the intrinsic fitness effects of mutations scale by the distance of the current background to a hypothesized boundary. We use simulations and theoretical analyses to explore the basic properties of the stickbreaking model such as fitness trajectories, the distribution of fitness achieved, and epistasis. Stickbreaking is compared to the additive and multiplicative models. We conclude that the stickbreaking model is qualitatively consistent with several commonly observed patterns of adaptive evolution.

Quasispecies as a matter of fact: viruses and beyond

We review the origins of the quasispecies concept and its relevance for RNA virus evolution, viral pathogenesis and antiviral treatment strategies. We emphasize a critical point of quasispecies that refers to genome collectivities as the unit of selection, and establish parallels between RNA viruses and some cellular systems such as bacteria and tumor cells. We refer also to tantalizing new observations that suggest quasispecies behavior in prions, perhaps as a result of the same quantum-mechanical indeterminations that underlie protein conformation and error-prone replication in genetic systems. If substantiated, these observations with prions could lead to new research on the structure-function relationship of non-nucleic acid biological molecules.

Two optimal mutation rates in obligate pathogens subject to deleterious mutation

Pathogen species with high mutation rates are likely to accumulate deleterious mutations that reduce their reproductive potential within the host. By altering the within-host growth rate of the pathogen, the deleterious mutation load has the potential to affect epidemiological properties such as prevalence, mean pathogen load, and the mean duration of infections. Here, I examine an epidemiological model that allows for multiple segregating mutations that affect within-host replication efficiency. The model demonstrates a complex range of outcomes depending on pathogen mutation rate, including two distinct, widely separated mutation rates associated with high pathogen prevalence. The low mutation rate prevalence peak is associated with small amounts of genetic diversity within the pathogen population, relatively stable prevalence and infection dynamics, and genetic variation partitioned between hosts. The high mutation rate peak is characterized by considerable genetic diversity both within and between hosts, relatively frequent invasions by more virulent types, and is qualitatively similar to an RNA virus quasispecies. The two prevalence peaks are separated by a valley where natural selection favors evolution toward the optimal within-host state, which is associated with high virulence and relatively rapid host mortality. Both chronic and acute infections are examined using stochastic forward simulations.

Profound differences in virus population genetics correspond to protection from CD4 decline resulting from feline lentivirus coinfection

CD4 decline is a hallmark of disease onset in individuals infected with Feline Immunodeficiency Virus (FIV) or Human Immunodeficiency Virus type 1 (HIV-1). Cats that are infected with a poorly replicating, apathogenic FIV (PLV) prior to exposure to a virulent FIV strain (FIVC) maintain CD4 numbers by mechanisms that are not correlated with a measurable adaptive immune response or reduction in circulating viral load. We employed population genetic approaches based on the 3′ portion of the viral genome to estimate the population structure of FIVC from single and dual infected cats. In dual infected cats, FIVC effective population size was decreased during the initial viral expansion phase, and after three weeks of infection, the population declined sharply. The FIVC population recovered to pre-bottleneck levels approximately seven weeks post-FIVC infection. However, the population emerging from the bottleneck in dual infected cats was distinct based on estimates of temporal population structure and substitution profiles. The transition to transversion rate ratio (κ) increased from early to late phases in dual infected cats due primarily to a decrease in transversions whereas in single infected cats, κ declined over time. Although one clone with extensive G to A substitutions, indicative of host cytidine deaminase editing, was recovered from a dual infected cat during the bottleneck, the post bottleneck population had an overall reduction in G to A substitutions. These data are consistent with a model of PLV-induced host restriction, putatively involving host DNA editing, that alters the dynamics of FIVC throughout the course of infection leading to disease attenuation.

Dynamics of in vitro fitness recovery of HIV-1

The study on the evolutionary consequences of an RNA viral population's fluctuations can be approached by in vitro experiments. This work describes the fitness recovery of HIV-1 after 20 large-population passages in 10 debilitated clones. The serial passages promoted an increase in viral fitness. In addition, we detected a significant number of mutations fixed in the complete genome consensus sequence of the final viral populations. Among the mutations, events of convergent evolution with important phenotypic characteristics occurred in several independent clones. One common change, V35I, in the nuclear localization signal of the p17 protein appeared in four viruses of three different lineages. Other common alterations mapped in position E196K of the reverse transcriptase or in position S316K of the V3 loop of the gp120 residue that is associated with the X4/R5 phenotype. Together with this mutational analysis, we studied the quasispecies heterogeneity of the initial and final viruses, revealing that fitness increase correlated with an augmentation in the genetic heterogeneity of viral quasispecies. However, while heterogeneity was mostly composed of synonymous (dS) mutations in the first 10 passages performed, at passage 21 it switched to nonsynonymous (dN) substitutions, with significant differences in dN - dS values between passages 11 and 21. In summary, the HIV-1 in vitro fitness recovery depicts a multiphase process occurring first by generation of mutations followed by fixation of the beneficial ones, depicting a classical Darwinian process.

Mechanisms of viral emergence

A number of virologic and environmental factors are involved in the emergence and re-emergence of viral disease. Viruses do not conservatively occupy a single and permanent ecological niche. Rather, due to their intrinsic capacity for genetic change, and to the evolvability of fitness levels, viruses display a potential to parasitize alternative host species. Mutation, recombination and genome segment reassortment, and combination of these molecular events, produce complex and phenotypically diverse populations of viruses, which constitute the raw material on which selection acts. The majority of emerging viral diseases of humans have a zoonotic origin. Sociologic and ecologic factors produce diverse and changing environments in which viral subpopulations have ample opportunities to be selected from intrinsically heterogeneous viral populations, particularly in the case of RNA viruses. In this manner, new human, animal and plant viruses have emerged periodically and, from all evidence, will continue to emerge. This article reviews some of the mechanisms that have been identified in viral emergence, with a focus on the importance of genetic variation of viruses, and on the general concept of biological complexity.

Diminishing returns from beneficial mutations and pervasive epistasis shape the fitness landscape for rifampicin resistance in Pseudomonas aeruginosa

Because adaptation depends upon the fixation of novel beneficial mutations, the fitness effects of beneficial mutations that are substituted by selection are key to our understanding of the process of adaptation. In this study, we experimentally investigated the fitness effects of beneficial mutations that are substituted when populations of the pathogenic bacterium Pseudomonas aeruginosa adapt to the antibiotic rifampicin. Specifically, we isolated the first beneficial mutation to be fixed by selection when 96 populations of three different genotypes of P. aeruginosa that vary considerably in fitness in the presence of rifampicin were challenged with adapting to a high dose of this antibiotic. The simple genetics of rifampicin resistance allowed us to determine the genetic basis of adaptation in the majority of our populations. We show that the average fitness effects of fixed beneficial mutations show a simple and clear pattern of diminishing returns, such that selection tends to fix mutations with progressively smaller effects as populations approach a peak on the adaptive landscape. The fitness effects of individual mutations, on the other hand, are highly idiosyncratic across genetic backgrounds, revealing pervasive epistasis. In spite of this complexity of genetic interactions in this system, there is an overall tendency toward diminishing-returns epistasis. We argue that a simple overall pattern of diminishing-returns adaptation emerges, despite pervasive epistasis between beneficial mutations, because many beneficial mutations are available, and while the fitness landscape is rugged at the fine scale, it is smooth and regular when we consider the average over possible routes to adaptation. In the context of antibiotic resistance, these results show that acquiring mutations that confer low levels of antibiotic resistance does not impose any constraint on the ability to evolve high levels of resistance.

Virulence correlates with fitness in vivo for two M group genotypes of Infectious hematopoietic necrosis virus (IHNV)

The nature of the association between viral fitness and virulence remains elusive in vertebrate virus systems, partly due to a lack of in vivo experiments using statistically sufficient numbers of replicate hosts. We examined the relationship between virulence and fitness in Infectious hematopoietic necrosis virus (IHNV), in vivo, in intact living rainbow trout. Trout were infected with a high or low virulence genotype of M genogroup IHNV, or a mixture of the two genotypes, so as to calculate relative fitness and the effect of a competition environment on fitness. Fitness was measured as total viral load in the host at time of peak viral density, quantified by genotype-specific quantitative RT-PCR (qRT-PCR). The more virulent IHNV genotype reached higher densities in both single and mixed infections. There was no effect of competition on the performance of either genotype. Our results suggest a positive link between IHNV genotype fitness and virulence.

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.

Genomic evolution of vesicular stomatitis virus strains with differences in adaptability

Virus strains with a history of repeated genetic bottlenecks frequently show a diminished ability to adapt compared to strains that do not have such a history. These differences in adaptability suggest differences in either the rate at which beneficial mutations are produced, the effects of beneficial mutations, or both. We tested these possibilities by subjecting four populations (two controls and two mutants with lower adaptabilities) to multiple replicas of a regimen of positive selection and then determining the fitnesses of the progeny through time and the changes in the consensus, full-length sequences of 56 genomes. We observed that at a given number of passages, the overall fitness gains observed for control populations were larger than fitness gains in mutant populations. However, these changes did not correlate with differences in the numbers of mutations accumulated in the two types of genomes. This result is consistent with beneficial mutations having a lower beneficial effect on mutant strains. Despite the overall fitness differences, some replicas of one mutant strain at passage 50 showed fitness increases similar to those observed for the wild type. We hypothesized that these evolved, high-fitness mutants may have a lower robustness than evolved, high-fitness controls. Robustness is the ability of a virus to avoid phenotypic changes in the face of mutation. We confirmed our hypothesis in mutation-accumulation experiments that showed a normalized fitness loss that was significantly larger in mutant bottlenecked populations than in control populations.

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.