Rate of novel host invasion affects adaptability of evolving RNA virus lineages

Although differing rates of environmental turnover should be consequential for the dynamics of adaptive change, this idea has been rarely examined outside of theory. In particular, the importance of RNA viruses in disease emergence warrants experiments testing how differing rates of novel host invasion may impact the ability of viruses to adaptively shift onto a novel host. To test whether the rate of environmental turnover influences adaptation, we experimentally evolved 144 Sindbis virus lineages in replicated tissue-culture environments, which transitioned from being dominated by a permissive host cell type to a novel host cell type. The rate at which the novel host ‘invaded’ the environment varied by treatment. The fitness (growth rate) of evolved virus populations was measured on each host type, and molecular substitutions were mapped via whole genome consensus sequencing. Results showed that virus populations more consistently reached high fitness levels on the novel host when the novel host ‘invaded’ the environment more gradually, and gradual invasion resulted in less variable genomic outcomes. Moreover, virus populations that experienced a rapid shift onto the novel host converged upon different genotypes than populations that experienced a gradual shift onto the novel host, suggesting a strong effect of historical contingency.     

Clonal interference is alleviated by high mutation rates in large populations

When a beneficial mutation is fixed in a population that lacks recombination, the genetic background linked to that mutation is fixed. As a result, beneficial mutations on different backgrounds experience competition, or "clonal interference," that can cause asexual populations to evolve more slowly than their sexual counterparts. Factors such as a large population size (N) and high mutation rates (mu) increase the number of competing beneficial mutations, and hence are expected to increase the intensity of clonal interference. However, recent theory suggests that, with very large values of Nmu, the severity of clonal interference may instead decline. The reason is that, with large Nmu, genomes including both beneficial mutations are rapidly created by recurrent mutation, obviating the need for recombination. Here, we analyze data from experimentally evolved asexual populations of a bacteriophage and find that, in these nonrecombining populations with very large Nmu, recurrent mutation does appear to ameliorate this cost of asexuality.     

Fitness recovery and compensatory evolution in natural mutant lines of C. elegans

Deleterious mutation accumulation plays a central role in evolutionary genetics, conservation biology, human health, and evolutionary medicine (e.g., methods of viral attenuation for live vaccines). It is therefore important to understand whether and how quickly populations with accumulated deleterious mutational loads can recover fitness through adaptive evolution. We used laboratory experimental evolution with four long-term mutation-accumulation (MA) lines of Caenorhabditis elegans nematodes to study the dynamics of such fitness evolution. We previously showed that when homozygous mutant populations are evolved in large population sizes, they can rapidly achieve wild-type fitness through the accumulation of new beneficial or compensatory epistatic mutations. Here, we expand this approach to demonstrate that when replicate lineages are initiated from the same mutant genotype, phenotypic evolution is only sometimes repeatable. MA genotypes that recovered ancestral fitness in the previous experiment did not always do so here. Further, the pattern of adaptive evolution in independently evolved replicates was contingent upon the MA genotype and varied among fitness-related traits. Our findings suggest that new beneficial mutations can drive rapid fitness evolution, but that the adaptive process is rendered somewhat unpredictable by its susceptibility to chance events and sensitivity to the evolutionary history of the starting population.     

Prior evolution in stochastic versus constant temperatures affects RNA virus evolvability at a thermal extreme

It is unclear how historical adaptation versus maladaptation in a prior environment affects population evolvability in a novel habitat. Prior work showed that vesicular stomatitis virus (VSV) populations evolved at constant 37°C improved in cellular infection at both 29°C and 37°C; in contrast, those evolved under random changing temperatures between 29°C and 37°C failed to improve. Here, we tested whether prior evolution affected the rate of adaptation at the thermal‐niche edge: 40°C. After 40 virus generations in the new environment, we observed that populations historically evolved at random temperatures showed greater adaptability. Deep sequencing revealed that most of the newly evolved mutations were de novo. Also, two novel evolved mutations in the VSV glycoprotein and replicase genes tended to co‐occur in the populations previously evolved at constant 37°C, whereas this parallelism was not seen in populations with prior random temperature evolution. These results suggest that prior adaptation under constant versus random temperatures constrained the mutation landscape that could improve fitness in the novel 40°C environment, perhaps owing to differing epistatic effects of new mutations entering genetic architectures that earlier diverged. We concluded that RNA viruses maladapted to their previous environment could “leapfrog” over counterparts of higher fitness, to achieve faster adaptability in a novel environment.     

Limits to adaptation in asexual populations

In asexual populations, the rate of adaptation is basically limited by the frequency and properties of spontaneous beneficial mutations. Hence, knowledge of these mutational properties and how they are affected by particular evolutionary conditions is a precondition for understanding the process of adaptation. Here, we address how the rate of adaptation of asexual populations is limited by its population size and mutation rate, as well as by two factors affecting the fraction of mutations that confer a benefit, i.e. the initial adaptedness of the population and the variability of the environment. These factors both influence which mutations are likely to occur, as well as the probability that they will ultimately contribute to adaptation. We attempt to separate the consequences of these basic population features in terms of their effect on the rate of adaptation by using results from evolution experiments with microorganisms.     

The Impact of Macroscopic Epistasis on Long-Term Evolutionary Dynamics

Genetic interactions can strongly influence the fitness effects of individual mutations, yet the impact of these epistatic interactions on evolutionary dynamics remains poorly understood. Here we investigate the evolutionary role of epistasis over 50,000 generations in a well-studied laboratory evolution experiment in Escherichia coli. The extensive duration of this experiment provides a unique window into the effects of epistasis during long-term adaptation to a constant environment. Guided by analytical results in the weak-mutation limit, we develop a computational framework to assess the compatibility of a given epistatic model with the observed patterns of fitness gain and mutation accumulation through time. We find that a decelerating fitness trajectory alone provides little power to distinguish between competing models, including those that lack any direct epistatic interactions between mutations. However, when combined with the mutation trajectory, these observables place strong constraints on the set of possible models of epistasis, ruling out many existing explanations of the data. Instead, we find that the data are consistent with a “two-epoch” model of adaptation, in which an initial burst of diminishing-returns epistasis is followed by a steady accumulation of mutations under a constant distribution of fitness effects. Our results highlight the need for additional DNA sequencing of these populations, as well as for more sophisticated models of epistasis that are compatible with all of the experimental data.     

Magnitude and sign epistasis among deleterious mutations in a positive-sense plant RNA virus

How epistatic interactions between mutations determine the genetic architecture of fitness is of central importance in evolution. The study of epistasis is particularly interesting for RNA viruses because of their genomic compactness, lack of genetic redundancy, and apparent low complexity. Moreover, interactions between mutations in viral genomes determine traits such as resistance to antiviral drugs, virulence and host range. In this study we generated 53 Tobacco etch potyvirus genotypes carrying pairs of single-nucleotide substitutions and measured their separated and combined deleterious fitness effects. We found that up to 38% of pairs had significant epistasis for fitness, including both positive and negative deviations from the null hypothesis of multiplicative effects. Interestingly, the sign of epistasis was correlated with viral protein-protein interactions in a model network, being predominantly positive between linked pairs of proteins and negative between unlinked ones. Furthermore, 55% of significant interactions were cases of reciprocal sign epistasis (RSE), indicating that adaptive landscapes for RNA viruses maybe highly rugged. Finally, we found that the magnitude of epistasis correlated negatively with the average effect of mutations. Overall, our results are in good agreement to those previously reported for other viruses and further consolidate the view that positive epistasis is the norm for small and compact genomes that lack genetic robustness.     

RNA editing by the host ADAR system affects the molecular evolution of the Zika virus

Zika virus (ZIKV) is a mosquito‐transmitted flavivirus, linked to microcephaly and fetal death in humans. Here, we investigate whether host‐mediated RNA editing of adenosines (ADAR) plays a role in the molecular evolution of ZIKV. Using complete coding sequences for the ZIKV polyprotein, we show that potential ADAR substitutions are underrepresented at the ADAR‐resistant GA dinucleotides of both the positive and negative strands, that these changes are spatially and temporally clustered (as expected of ADAR editing) for three branches of the viral phylogeny, and that ADAR mutagenesis can be linked to its codon usage. Furthermore, resistant GA dinucleotides are enriched on the positive (but not negative) strand, indicating that the former is under stronger purifying selection than the latter. ADAR editing also affects the evolution of the rhabdovirus sigma. Our study now documents that host ADAR editing is a mutation and evolutionary force of positive‐ as well as negative‐strand RNA viruses.     

Phylogenomics and molecular evolution of foot-and-mouth disease virus

This report describes the use of Bayesian methods to analyze polyprotein coding region sequences (n = 217) obtained from GenBank to define the genome-wide phylogeny of foot and mouth disease virus (FMDV). The results strongly supported the monophyly of five FMDV serotypes, O, A, Asia 1, C, and SAT 3, while sequences for the two remaining FMDV serotypes, SAT 1 and SAT 2 did not separate into entirely distinct clades. The phylogenomic tree revealed three sister-group relationships, serotype O + Asia 1, A + C, and SAT 1 + 3 + 2, with a new branching pattern: {[(O, Asia 1), (A, C)], (SAT 1, 2, 3)}. Within each serotype, there was no apparent periodic, geographic, or host species influence on the evolution of global FMDVs. Analysis of the polyprotein coding region of these sequences provided evidence for the influence of purifying selection on the evolution of FMDV. Using a Bayesian coalescent approach, the evolutionary rate of FMDV isolates that circulated during the years 1932-2007 was estimated to be 1.46 × 10(-3) substitutions/site/year, and the most recent common ancestor of the virus existed approximately 481 years ago. Bayesian skyline plot revealed a population expansion in the early 20(th) century that was followed by a rapid decline in population size from the late 20(th) century to the present day. These findings provide new insights into the mechanisms that impact on the evolution of this important livestock pathogen.     

Delayed transmission selects for increased survival of vesicular stomatitis virus

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

A phylogenetic method for detecting positive epistasis in gene sequences and its application to RNA virus evolution

RNA virus genomes are compact, often containing multiple overlapping reading frames and functional secondary structure. Consequently, it is thought that evolutionary interactions between nucleotide sites are commonplace in the genomes of these infectious agents. However, the role of epistasis in natural populations of RNA viruses remains unclear. To investigate the pervasiveness of epistasis in RNA viruses, we used a parsimony-based computational method to identify pairs of co-occurring mutations along phylogenies of 177 RNA virus genes. This analysis revealed widespread evidence for positive epistatic interactions at both synonymous and nonsynonymous nucleotide sites and in both clonal and recombining viruses, with the majority of these interactions spanning very short sequence regions. These findings have important implications for understanding the key aspects of RNA virus evolution, including the dynamics of adaptation. Additionally, many comparative analyses that utilize the phylogenetic relationships among gene sequences assume that mutations represent independent, uncorrelated events. Our results show that this assumption may often be invalid.     

The impact of high‐order epistasis in the within‐host fitness of a positive‐sense plant RNA virus

RNA viruses are the main source of emerging infectious diseases because of the evolutionary potential bestowed by their fast replication, large population sizes and high mutation and recombination rates. However, an equally important property, which is usually neglected, is the topography of the fitness landscape. How many fitness maxima exist and how well they are connected is especially interesting, as this determines the number of accessible evolutionary pathways. To address this question, we have reconstructed a region of the fitness landscape of tobacco etch potyvirus constituted by mutations observed during the experimental adaptation of the virus to the novel host Arabidopsis thaliana. Fitness was measured for many genotypes and showed the existence of multiple peaks and holes in the landscape. We found prevailing epistatic effects between mutations, with cases of reciprocal sign epistasis being common among pairs of mutations. We also found that high‐order epistasis was as important as pairwise epistasis in their contribution to fitness. Therefore, results suggest that the landscape was rugged due to the existence of holes caused by lethal genotypes, that a very limited number of potential neutral paths exist and that it contained a single adaptive peak.     

Coinfection and superinfection in RNA virus populations: a selection-mutation model

In this paper, we present a general selection-mutation model of evolution on a one-dimensional continuous fitness space. The formulation of our model includes both the classical diffusion approach to mutation process as well as an alternative approach based on an integral operator with a mutation kernel. We show that both approaches produce fundamentally equivalent results. To illustrate the suitability of our model, we focus its analytical study into its application to recent experimental studies of in vitro viral evolution. More specifically, these experiments were designed to test previous theoretical predictions regarding the effects of multiple infection dynamics (i.e., coinfection and superinfection) on the virulence of evolving viral populations. The results of these experiments, however, did not match with previous theory. By contrast, the model we present here helps to understand the underlying viral dynamics on these experiments and makes new testable predictions about the role of parameters such the time between successive infections and the growth rates of resident and invading populations.     

Rapid evolution of adaptive niche construction in experimental microbial populations

Many species engage in adaptive niche construction: modification of the local environment that increases the modifying organism's competitive fitness. Adaptive niche construction provides an alternative pathway to higher fitness, shaping the environment rather than conforming to it. Yet, experimental evidence for the evolutionary emergence of adaptive niche construction is lacking, leaving its role in evolution uncertain. Here we report a direct observation of the de novo evolution of adaptive niche construction in populations of the bacteria Pseudomonas fluorescens. In a laboratory experiment, we allowed several bacterial populations to adapt to a novel environment and assessed whether niche construction evolved over time. We found that adaptive niche construction emerged rapidly, within approximately 100 generations, and became ubiquitous after approximately 400 generations. The large fitness effect of this niche construction was dominated by the low fitness of evolved strains in the ancestrally modified environment: evolved niche constructors were highly dependent on their specific environmental modifications. Populations were subjected to frequent resetting of environmental conditions and severe reduction of spatial habitat structure, both of which are thought to make adaptive niche construction difficult to evolve. Our finding that adaptive niche construction nevertheless evolved repeatably suggests that it may play a more important role in evolution than generally thought.     

Synchronous Waves of Failed Soft Sweeps in the Laboratory: Remarkably Rampant Clonal Interference of Alleles at a Single Locus

It has increasingly been recognized that adapting populations of microbes contain not one, but many lineages continually arising and competing at once. This process, termed “clonal interference,” alters the rate and dynamics of adaptation and biases winning mutations toward those with the largest selective effect. Here we uncovered a dramatic example of clonal interference between multiple similar mutations occurring at the same locus within replicate populations of Methylobacterium extorquens AM1. Because these mutational events involved the transposition of an insertion sequence into a narrow window of a single gene, they were both readily detectable at low frequencies and could be distinguished due to differences in insertion sites. This allowed us to detect up to 17 beneficial alleles of this type coexisting in a single population. Despite conferring a large selective benefit, the majority of these alleles rose and then fell in frequency due to other lineages emerging that were more fit. By comparing allele-frequency dynamics to the trajectories of fitness gains by these populations, we estimated the fitness values of the genotypes that contained these mutations. Collectively across all populations, these alleles arose upon backgrounds with a wide range of fitness values. Within any single population, however, multiple alleles tended to rise and fall synchronously during a single wave of multiple genotypes with nearly identical fitness values. These results suggest that alleles of large benefit arose repeatedly in failed “soft sweeps” during narrow windows of adaptation due to the combined effects of epistasis and clonal interference.     

Patterns of epistasis in RNA viruses: a review of the evidence from vaccine design

Epistasis results when the fitness effects of a mutation change depending on the presence or absence of other mutations in the genome. The predictions of many influential evolutionary hypotheses are determined by the existence and form of epistasis. One rich source of data on the interactions among deleterious mutations that has gone untapped by evolutionary biologists is the literature on the design of live, attenuated vaccine viruses. Rational vaccine design depends upon the measurement of individual and combined effects of deleterious mutations. In the current study, we have reviewed data from 29 vaccine-oriented studies using 14 different RNA viruses. Our analyses indicate that (1) no consistent tendency towards a particular form of epistasis exists across RNA viruses and (2) significant interactions among groups of mutations within individual viruses occur but are not common. RNA viruses are significant pathogens of human disease, and are tractable model systems for evolutionary studies--we discuss the relevance of our findings in both contexts.     

Evaluating the within-host fitness effects of mutations fixed during virus adaptation to different ecotypes of a new host

The existence of genetic variation for resistance in host populations is assumed to be essential to the spread of an emerging virus. Models predict that the rate of spread slows down with the increasing frequency and higher diversity of resistance alleles in the host population. We have been using the experimental pathosystem Arabidopsis thaliana—tobacco etch potyvirus (TEV) to explore the interplay between genetic variation in host''s susceptibility and virus diversity. We have recently shown that TEV populations evolving in A. thaliana ecotypes that differ in susceptibility to infection gained within-host fitness, virulence and infectivity in a manner compatible with a gene-for-gene model of host–parasite interactions: hard-to-infect ecotypes were infected by generalist viruses, whereas easy-to-infect ecotypes were infected by every virus. We characterized the genomes of the evolved viruses and found cases of host-driven convergent mutations. To gain further insights in the mechanistic basis of this gene-for-gene model, we have generated all viral mutations individually as well as in specific combinations and tested their within-host fitness effects across ecotypes. Most of these mutations were deleterious or neutral in their local ecotype and only a very reduced number had a host-specific beneficial effect. We conclude that most of the mutations fixed during the evolution experiment were so by drift or by selective sweeps along with the selected driver mutation. In addition, we evaluated the ruggedness of the underlying adaptive fitness landscape and found that mutational effects were mostly multiplicative, with few cases of significant epistasis.     

The ecology and genetics of fitness in Chlamydomonas. XIII. Fitness of long-term sexual and asexual populations in benign environments

We measured the mean fitness of populations of Chlamydomonas reinhardtii maintained in the laboratory as obligately sexual or asexual populations for about 100 sexual cycles and about 1000 asexual generations. Sexuality (random gamete fusion followed by meiosis) is expected to reduce mutational load and increase mean fitness by combining deleterious mutations from different lines of descent. We found no evidence for this process of mutation clearance: the mean fitness of sexual populations did not exceed that of asexual populations, whether measured through competition or in pure culture. We found instead that sexual progeny suffer an immediate loss in fitness, and that sexual lines maintain genetic variance for fitness. We suggest that sexual populations at equilibrium with selection in a benign environment may be mixtures of several or many epistatic genotypes with nearly equal fitness. Recombination between these genotypes reduces mean fitness and creates genetic variance for fitness. This may provide fuel for continued selection should the environment change.     

Exploring the evolution of multicellularity in Saccharomyces cerevisiae under bacteria environment: An experimental phylogenetics approach

There have been over 25 independent unicellular to multicellular evolutionary transitions, which have been transformational in the complexity of life. All of these transitions likely occurred in communities numerically dominated by unicellular organisms, mostly bacteria. Hence, it is reasonable to expect that bacteria were involved in generating the ecological conditions that promoted the stability and proliferation of the first multicellular forms as protective units. In this study, we addressed this problem by analyzing the occurrence of multicellularity in an experimental phylogeny of yeasts (Sacharomyces cerevisiae) a model organism that is unicellular but can generate multicellular clusters under some conditions. We exposed a single ancestral population to periodic divergences, coevolving with a cocktail of environmental bacteria that were inoculated to the environment of the ancestor, and compared to a control (no bacteria). We quantified culturable microorganisms to the level of genera, finding up to 20 taxa (all bacteria) that competed with the yeasts during diversification. After 600 generations of coevolution, the yeasts produced two types of multicellular clusters: clonal and aggregative. Whereas clonal clusters were present in both treatments, aggregative clusters were only present under the bacteria treatment and showed significant phylogenetic signal. However, clonal clusters showed different properties if bacteria were present as follows: They were more abundant and significantly smaller than in the control. These results indicate that bacteria are important modulators of the occurrence of multicellularity, providing support to the idea that they generated the ecological conditions‐promoting multicellularity.