Distribution of mutational fitness effects and of epistasis in the 5' untranslated region of a plant RNA virus

Experimental evolution of cowpea mild mottle virus reveals recombination-driven reduction in virulence accompanied by increases in diversity and viral fitness

Major themes in pathogen evolution are emergence, evolution of virulence, host adaptation and the processes that underlie them. RNA viruses are of particular interest due to their rapid evolution. The in vivo molecular evolution of an RNA plant virus was demonstrated here using a necrotic isolate of cowpea mild mottle virus (CPMMV) and a susceptible soybean genotype submitted to serial inoculations. We show that the virus lost the capacity to cause necrosis after six passages through the host plant. When a severe bottleneck was imposed, virulence reduction occurred in the second passage. The change to milder symptoms had fitness benefits for the virus (higher RNA accumulation) and for its vector, the whitefly Bemisia tabaci. Genetic polymorphisms were highest in ORF1 (viral replicase) and were independent of the symptom pattern. Recombination was a major contributor to this diversity - even with the strong genetic bottleneck, recombination events and hot spots were detected within ORF1. Virulence reduction was associated with different sites in ORF1 associated to recombination events in both experiments. Overall, the results demonstrate that the reduction in virulence was a consequence of the emergence of new variants, driven by recombination. Besides providing details of the evolutionary mechanisms behind a reduction in virulence and its effect under viral and vector fitness, we propose that this recombination-driven switch in virulence allows the pathogen to rapidly adapt to a new host and, potentially, switch back.

Recent insights into the genotype-phenotype relationship from massively parallel genetic assays

With the molecular revolution in Biology, a mechanistic understanding of the genotype-phenotype relationship became possible. Recently, advances in DNA synthesis and sequencing have enabled the development of deep mutational scanning assays, capable of scoring comprehensive libraries of genotypes for fitness and a variety of phenotypes in massively parallel fashion. The resulting empirical genotype-fitness maps pave the way to predictive models, potentially accelerating our ability to anticipate the behaviour of pathogen and cancerous cell populations from sequencing data. Besides from cellular fitness, phenotypes of direct application in industry (e.g. enzyme activity) and medicine (e.g. antibody binding) can be quantified and even selected directly by these assays. This review discusses the technological basis of and recent developments in massively parallel genetics, along with the trends it is uncovering in the genotype-phenotype relationship (distribution of mutation effects, epistasis), their possible mechanistic bases and future directions for advancing towards the goal of predictive genetics.

Developmental structuring of phenotypic variation: A case study with a cellular automata model of ontogeny

Developmental mechanisms not only produce an organismal phenotype, but they also structure the way genetic variation maps to phenotypic variation. Here, we revisit a computational model for the evolution of ontogeny based on cellular automata, in which evolution regularly discovered two alternative mechanisms for achieving a selected phenotype, one showing high modularity, the other showing morphological integration. We measure a primary variational property of the systems, their distribution of fitness effects of mutation. We find that the modular ontogeny shows the evolution of mutational robustness and ontogenic simplification, while the integrated ontogeny does not. We discuss the wider use of this methodology on other computational models of development as well as real organisms.

Evolutionary Dynamics in the RNA Bacteriophage Qβ Depends on the Pattern of Change in Selective Pressures

The rate of change in selective pressures is one of the main factors that determines the likelihood that populations can adapt to stress conditions. Generally, the reduction in the population size that accompanies abrupt environmental changes makes it difficult to generate and select adaptive mutations. However, in systems with high genetic diversity, as happens in RNA viruses, mutations with beneficial effects under new conditions can already be present in the population, facilitating adaptation. In this work, we have propagated an RNA bacteriophage (Qβ) at temperatures higher than the optimum, following different patterns of change. We have determined the fitness values and the consensus sequences of all lineages throughout the evolutionary process in order to establish correspondences between fitness variations and adaptive pathways. Our results show that populations subjected to a sudden temperature change gain fitness and fix mutations faster than those subjected to gradual changes, differing also in the particular selected mutations. The life-history of populations prior to the environmental change has great importance in the dynamics of adaptation. The conclusion is that in the bacteriophage Qβ, the standing genetic diversity together with the rate of temperature change determine both the rapidity of adaptation and the followed evolutionary pathways.

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.

Genomic Loss and Silencing on the Y Chromosomes of Rumex

Across many unrelated lineages of plants and animals, Y chromosomes show a recurrent pattern of gene degeneration and loss, but the relative importance of inefficient selection, adaptive gene silencing, and neutral genetic drift in causing degeneration remain poorly understood. Here, we use next-generation genome and transcriptome sequencing to investigate patterns of ongoing Y chromosome degeneration in two annual plant species of Rumex (Polygonaceae) differing in their degree of degeneration and sex chromosome heteromorphism. We find evidence for both gene loss as well as silencing in these young plant sex chromosomes. Our analyses revealed significantly more gene deletion relative to silencing in R. rothschildianus, which has had a larger nonrecombining region for a longer period than R. hastatulus, consistent with this system being at a more advanced stage of degeneration. Intra- and interspecific comparisons of genomic coverage and heterozygosity indicated that loss of expression precedes gene deletion, implying that the final stages of mutation accumulation and gene loss may often occur neutrally. We found no evidence for adaptive silencing of genes that have lost expression. Our results suggest that the initial spread of deleterious regulatory variants and/or epigenetic silencing may be an important driver of early degeneration of Y chromosomes.

Strain-dependent mutational effects for Pepino mosaic virus in a natural host

Background

Pepino mosaic virus (PepMV) is an emerging plant pathogen that infects tomatoes worldwide. Understanding the factors that influence its evolutionary success is essential for developing new control strategies that may be more robust against the evolution of new viral strains. One of these evolutionary factors is the distribution of mutational fitness effect (DMFE), that is, the fraction of mutations that are lethal, deleterious, neutral, and beneficial on a given viral strain and host species. The goal of this study was to characterize the DMFE of introduced nonsynonymous mutations on a mild isolate of PepMV from the Chilean 2 strain (PepMV-P22). Additionally, we also explored whether the fitness effect of a given mutation depends on the gene where it appears or on epistatic interactions with the genetic background. To address this latter possibility, a subset of mutations were also introduced in a mild isolate of the European strain (PepMV-P11) and the fitness of the resulting clones measured.

Results

A collection of 25 PepMV clones each containing a single nucleotide nonsynonymous substitution was created by site-directed mutagenesis and the fitness of each mutant was determined. PepMV-P22 genome showed a high degree of robustness against point mutations, with 80% of mutations being either neutral or even beneficial and only 20% being deleterious or lethal. We found that the effect of mutations strongly depended on the gene in which they were introduced. Mutations with the largest average beneficial effects were those affecting the RdRp gene, in contrast to mutations affecting TGB1 and CP genes, for which the average effects were deleterious. Moreover, significant epistatic interactions were observed between nonsynonymous mutations and the genetic background, meaning that the effect of a given nucleotide substitution on a particular genomic context cannot be predicted by knowing its effect in a different one.

Conclusions

Our results indicated that PepMV genome has a surprisingly high robustness against mutations. We also found that fitness consequences of a given mutation differ between the two strains analyzed. This discovery suggests that the strength of selection, and thus the rates of evolution, vary among PepMV strains.

Evolutionary ecology of virus emergence

The cross-species transmission of viruses into new host populations, termed virus emergence, is a significant issue in public health, agriculture, wildlife management, and related fields. Virus emergence requires overlap between host populations, alterations in virus genetics to permit infection of new hosts, and adaptation to novel hosts such that between-host transmission is sustainable, all of which are the purview of the fields of ecology and evolution. A firm understanding of the ecology of viruses and how they evolve is required for understanding how and why viruses emerge. In this paper, I address the evolutionary mechanisms of virus emergence and how they relate to virus ecology. I argue that, while virus acquisition of the ability to infect new hosts is not difficult, limited evolutionary trajectories to sustained virus between-host transmission and the combined effects of mutational meltdown, bottlenecking, demographic stochasticity, density dependence, and genetic erosion in ecological sinks limit most emergence events to dead-end spillover infections. Despite the relative rarity of pandemic emerging viruses, the potential of viruses to search evolutionary space and find means to spread epidemically and the consequences of pandemic viruses that do emerge necessitate sustained attention to virus research, surveillance, prophylaxis, and treatment.

On the importance of skewed offspring distributions and background selection in virus population genetics

Many features of virus populations make them excellent candidates for population genetic study, including a very high rate of mutation, high levels of nucleotide diversity, exceptionally large census population sizes, and frequent positive selection. However, these attributes also mean that special care must be taken in population genetic inference. For example, highly skewed offspring distributions, frequent and severe population bottleneck events associated with infection and compartmentalization, and strong purifying selection all affect the distribution of genetic variation but are often not taken into account. Here, we draw particular attention to multiple-merger coalescent events and background selection, discuss potential misinference associated with these processes, and highlight potential avenues for better incorporating them into future population genetic analyses.

A Statistical Guide to the Design of Deep Mutational Scanning Experiments

The characterization of the distribution of mutational effects is a key goal in evolutionary biology. Recently developed deep-sequencing approaches allow for accurate and simultaneous estimation of the fitness effects of hundreds of engineered mutations by monitoring their relative abundance across time points in a single bulk competition. Naturally, the achievable resolution of the estimated fitness effects depends on the specific experimental setup, the organism and type of mutations studied, and the sequencing technology utilized, among other factors. By means of analytical approximations and simulations, we provide guidelines for optimizing time-sampled deep-sequencing bulk competition experiments, focusing on the number of mutants, the sequencing depth, and the number of sampled time points. Our analytical results show that sampling more time points together with extending the duration of the experiment improves the achievable precision disproportionately compared with increasing the sequencing depth or reducing the number of competing mutants. Even if the duration of the experiment is fixed, sampling more time points and clustering these at the beginning and the end of the experiment increase experimental power and allow for efficient and precise assessment of the entire range of selection coefficients. Finally, we provide a formula for calculating the 95%-confidence interval for the measurement error estimate, which we implement as an interactive web tool. This allows for quantification of the maximum expected a priori precision of the experimental setup, as well as for a statistical threshold for determining deviations from neutrality for specific selection coefficient estimates.