Experimental Rugged Fitness Landscape in Protein Sequence Space

The fitness landscape in sequence space determines the process of biomolecular evolution. To plot the fitness landscape of protein function, we carried out in vitro molecular evolution beginning with a defective fd phage carrying a random polypeptide of 139 amino acids in place of the g3p minor coat protein D2 domain, which is essential for phage infection. After 20 cycles of random substitution at sites 12–130 of the initial random polypeptide and selection for infectivity, the selected phage showed a 1.7×10(4)-fold increase in infectivity, defined as the number of infected cells per ml of phage suspension. Fitness was defined as the logarithm of infectivity, and we analyzed (1) the dependence of stationary fitness on library size, which increased gradually, and (2) the time course of changes in fitness in transitional phases, based on an original theory regarding the evolutionary dynamics in Kauffman's n-k fitness landscape model. In the landscape model, single mutations at single sites among n sites affect the contribution of k other sites to fitness. Based on the results of these analyses, k was estimated to be 18–24. According to the estimated parameters, the landscape was plotted as a smooth surface up to a relative fitness of 0.4 of the global peak, whereas the landscape had a highly rugged surface with many local peaks above this relative fitness value. Based on the landscapes of these two different surfaces, it appears possible for adaptive walks with only random substitutions to climb with relative ease up to the middle region of the fitness landscape from any primordial or random sequence, whereas an enormous range of sequence diversity is required to climb further up the rugged surface above the middle region

Mosquitoes Put the Brake on Arbovirus Evolution: Experimental Evolution Reveals Slower Mutation Accumulation in Mosquito Than Vertebrate Cells

Like other arthropod-borne viruses (arboviruses), mosquito-borne dengue virus (DENV) is maintained in an alternating cycle of replication in arthropod and vertebrate hosts. The trade-off hypothesis suggests that this alternation constrains DENV evolution because a fitness increase in one host usually diminishes fitness in the other. Moreover, the hypothesis predicts that releasing DENV from host alternation should facilitate adaptation. To test this prediction, DENV was serially passaged in either a single human cell line (Huh-7), a single mosquito cell line (C6/36), or in alternating passages between Huh-7 and C6/36 cells. After 10 passages, consensus mutations were identified and fitness was assayed by evaluating replication kinetics in both cell types as well as in a novel cell type (Vero) that was not utilized in any of the passage series. Viruses allowed to specialize in single host cell types exhibited fitness gains in the cell type in which they were passaged, but fitness losses in the bypassed cell type, and most alternating passages, exhibited fitness gains in both cell types. Interestingly, fitness gains were observed in the alternately passaged, cloned viruses, an observation that may be attributed to the acquisition of both host cell–specific and amphi-cell-specific adaptations or to recovery from the fitness losses due to the genetic bottleneck of biological cloning. Amino acid changes common to both passage series suggested convergent evolution to replication in cell culture via positive selection. However, intriguingly, mutations accumulated more rapidly in viruses passed in Huh-7 cells than in those passed in C6/36 cells or in alternation. These results support the hypothesis that releasing DENV from host alternation facilitates adaptation, but there is limited support for the hypothesis that such alternation necessitates a fitness trade-off. Moreover, these findings suggest that patterns of genetic evolution may differ between viruses replicating in mammalian and mosquito cells

Experimental Passage of St. Louis Encephalitis Virus In Vivo in Mosquitoes and Chickens Reveals Evolutionarily Significant Virus Characteristics

St. Louis encephalitis virus (SLEV; Flaviviridae, flavivirus) was the major cause of epidemic flaviviral encephalitis in the U.S. prior to the introduction of West Nile virus (WNV) in 1999. However, outbreaks of SLEV have been significantly more limited then WNV in terms of levels of activity and geographic dispersal. One possible explanation for these variable levels of activity is that differences in the potential for each virus to adapt to its host cycle exist. The need for arboviruses to replicate in disparate hosts is thought to result in constraints on both evolution and host-specific adaptation. If cycling is the cause of genetic stability observed in nature and arboviruses lack host specialization, then sequential passage should result in both the accumulation of mutations and specialized viruses better suited for replication in that host. Previous studies suggest that WNV and SLEV differ in capacity for both genetic change and host specialization, and in the costs each accrues from specializing. In an attempt to clarify how selective pressures contribute to epidemiological patterns of WNV and SLEV, we evaluated mutant spectra size, consensus genetic change, and phenotypic changes for SLEV in vivo following 20 sequential passages via inoculation in either Culex pipiens mosquitoes or chickens. Results demonstrate that the capacity for genetic change is large for SLEV and that the size of the mutant spectrum is host-dependent using our passage methodology. Despite this, a general lack of consensus change resulted from passage in either host, a result that contrasts with the idea that constraints on evolution in nature result from host cycling alone. Results also suggest that a high level of adaptation to both hosts already exists, despite host cycling. A strain significantly more infectious in chickens did emerge from one lineage of chicken passage, yet other lineages and all mosquito passage strains did not display measurable host-specific fitness gains. In addition, increased infectivity in chickens did not decrease infectivity in mosquitoes, which further contrasts the concept of fitness trade-offs for arboviruses

Host Alternation Is Necessary to Maintain the Genome Stability of Rift Valley Fever Virus

Arthropod-borne viruses are transmitted among vertebrate hosts by insect vectors. Unusually, Rift Valley fever virus (RVFV) can also be transmitted by direct contacts of animals/humans with infectious tissues. What are the molecular mechanisms and evolutionary events leading to adopt one mode of transmission rather than the other? Viral replication is implied to be different in a vertebrate host and an invertebrate host. The alternating host cycle tends to limit virus evolution by adopting a compromise fitness level for replication in both hosts. To test this hypothesis, we used a cell culture model system to study the evolution of RVFV. We found that freeing RVFV from alternating replication in mammalian and mosquito cells led to large deletions in the NSs gene carrying the virulence factor. Resulting NSs-truncated viruses were able to protect mice from a challenge with a virulent RVFV. Thus, in nature, virulence is likely maintained by continuous alternating passages between vertebrates and insects. Thereby, depending on the mode of transmission adopted, the evolution of RVFV will be of major importance to predict the outcome of outbreaks

Effect of Host Species on the Distribution of Mutational Fitness Effects for an RNA Virus

Knowledge about the distribution of mutational fitness effects (DMFE) is essential for many evolutionary models. In recent years, the properties of the DMFE have been carefully described for some microorganisms. In most cases, however, this information has been obtained only for a single environment, and very few studies have explored the effect that environmental variation may have on the DMFE. Environmental effects are particularly relevant for the evolution of multi-host parasites and thus for the emergence of new pathogens. Here we characterize the DMFE for a collection of twenty single-nucleotide substitution mutants of Tobacco etch potyvirus (TEV) across a set of eight host environments. Five of these host species were naturally infected by TEV, all belonging to family Solanaceae, whereas the other three were partially susceptible hosts belonging to three other plant families. First, we found a significant virus genotype-by-host species interaction, which was sustained by differences in genetic variance for fitness and the pleiotropic effect of mutations among hosts. Second, we found that the DMFEs were markedly different between Solanaceae and non-Solanaceae hosts. Exposure of TEV genotypes to non-Solanaceae hosts led to a large reduction of mean viral fitness, while the variance remained constant and skewness increased towards the right tail. Within the Solanaceae hosts, the distribution contained an excess of deleterious mutations, whereas for the non-Solanaceae the fraction of beneficial mutations was significantly larger. All together, this result suggests that TEV may easily broaden its host range and improve fitness in new hosts, and that knowledge about the DMFE in the natural host does not allow for making predictions about its properties in an alternative host

Experimental Evolution of an Oncolytic Vesicular Stomatitis Virus with Increased Selectivity for p53-Deficient Cells

Experimental evolution has been used for various biotechnological applications including protein and microbial cell engineering, but less commonly in the field of oncolytic virotherapy. Here, we sought to adapt a rapidly evolving RNA virus to cells deficient for the tumor suppressor gene p53, a hallmark of cancer cells. To achieve this goal, we established four independent evolution lines of the vesicular stomatitis virus (VSV) in p53-knockout mouse embryonic fibroblasts (p53−/− MEFs) under conditions favoring the action of natural selection. We found that some evolved viruses showed increased fitness and cytotoxicity in p53−/− cells but not in isogenic p53+/+ cells, indicating gene-specific adaptation. However, full-length sequencing revealed no obvious or previously described genetic changes associated with oncolytic activity. Half-maximal effective dose (EC50) assays in mouse p53-positive colon cancer (CT26) and p53-deficient breast cancer (4T1) cells indicated that the evolved viruses were more effective against 4T1 cells than the parental virus or a reference oncolytic VSV (MΔ51), but showed no increased efficacy against CT26 cells. In vivo assays using 4T1 syngeneic tumor models showed that one of the evolved lines significantly delayed tumor growth compared to mice treated with the parental virus or untreated controls, and was able to induce transient tumor suppression. Our results show that RNA viruses can be specifically adapted typical cancer features such as p53 inactivation, and illustrate the usefulness of experimental evolution for oncolytic virotherapy

Variation in RNA Virus Mutation Rates across Host Cells

It is well established that RNA viruses exhibit higher rates of spontaneous mutation than DNA viruses and microorganisms. However, their mutation rates vary amply, from 10−6 to 10−4 substitutions per nucleotide per round of copying (s/n/r) and the causes of this variability remain poorly understood. In addition to differences in intrinsic fidelity or error correction capability, viral mutation rates may be dependent on host factors. Here, we assessed the effect of the cellular environment on the rate of spontaneous mutation of the vesicular stomatitis virus (VSV), which has a broad host range and cell tropism. Luria-Delbrück fluctuation tests and sequencing showed that VSV mutated similarly in baby hamster kidney, murine embryonic fibroblasts, colon cancer, and neuroblastoma cells (approx. 10−5 s/n/r). Cell immortalization through p53 inactivation and oxygen levels (1–21%) did not have a significant impact on viral replication fidelity. This shows that previously published mutation rates can be considered reliable despite being based on a narrow and artificial set of laboratory conditions. Interestingly, we also found that VSV mutated approximately four times more slowly in various insect cells compared with mammalian cells. This may contribute to explaining the relatively slow evolution of VSV and other arthropod-borne viruses in nature

Mutation Bias Favors Protein Folding Stability in the Evolution of Small Populations

Mutation bias in prokaryotes varies from extreme adenine and thymine (AT) in obligatory endosymbiotic or parasitic bacteria to extreme guanine and cytosine (GC), for instance in actinobacteria. GC mutation bias deeply influences the folding stability of proteins, making proteins on the average less hydrophobic and therefore less stable with respect to unfolding but also less susceptible to misfolding and aggregation. We study a model where proteins evolve subject to selection for folding stability under given mutation bias, population size, and neutrality. We find a non-neutral regime where, for any given population size, there is an optimal mutation bias that maximizes fitness. Interestingly, this optimal GC usage is small for small populations, large for intermediate populations and around 50% for large populations. This result is robust with respect to the definition of the fitness function and to the protein structures studied. Our model suggests that small populations evolving with small GC usage eventually accumulate a significant selective advantage over populations evolving without this bias. This provides a possible explanation to the observation that most species adopting obligatory intracellular lifestyles with a consequent reduction of effective population size shifted their mutation spectrum towards AT. The model also predicts that large GC usage is optimal for intermediate population size. To test these predictions we estimated the effective population sizes of bacterial species using the optimal codon usage coefficients computed by dos Reis et al. and the synonymous to non-synonymous substitution ratio computed by Daubin and Moran. We found that the population sizes estimated in these ways are significantly smaller for species with small and large GC usage compared to species with no bias, which supports our prediction

Quasispecies Spatial Models for RNA Viruses with Different Replication Modes and Infection Strategies

Empirical observations and theoretical studies suggest that viruses may use different replication strategies to amplify their genomes, which impact the dynamics of mutation accumulation in viral populations and therefore, their fitness and virulence. Similarly, during natural infections, viruses replicate and infect cells that are rarely in suspension but spatially organized. Surprisingly, most quasispecies models of virus replication have ignored these two phenomena. In order to study these two viral characteristics, we have developed stochastic cellular automata models that simulate two different modes of replication (geometric vs stamping machine) for quasispecies replicating and spreading on a two-dimensional space. Furthermore, we explored these two replication models considering epistatic fitness landscapes (antagonistic vs synergistic) and different scenarios for cell-to-cell spread, one with free superinfection and another with superinfection inhibition. We found that the master sequences for populations replicating geometrically and with antagonistic fitness effects vanished at low critical mutation rates. By contrast, the highest critical mutation rate was observed for populations replicating geometrically but with a synergistic fitness landscape. Our simulations also showed that for stamping machine replication and antagonistic epistasis, a combination that appears to be common among plant viruses, populations further increased their robustness by inhibiting superinfection. We have also shown that the mode of replication strongly influenced the linkage between viral loci, which rapidly reached linkage equilibrium at increasing mutations for geometric replication. We also found that the strategy that minimized the time required to spread over the whole space was the stamping machine with antagonistic epistasis among mutations. Finally, our simulations revealed that the multiplicity of infection fluctuated but generically increased along time

Infidelity of SARS-CoV Nsp14-Exonuclease Mutant Virus Replication Is Revealed by Complete Genome Sequencing

Most RNA viruses lack the mechanisms to recognize and correct mutations that arise during genome replication, resulting in quasispecies diversity that is required for pathogenesis and adaptation. However, it is not known how viruses encoding large viral RNA genomes such as the Coronaviridae (26 to 32 kb) balance the requirements for genome stability and quasispecies diversity. Further, the limits of replication infidelity during replication of large RNA genomes and how decreased fidelity impacts virus fitness over time are not known. Our previous work demonstrated that genetic inactivation of the coronavirus exoribonuclease (ExoN) in nonstructural protein 14 (nsp14) of murine hepatitis virus results in a 15-fold decrease in replication fidelity. However, it is not known whether nsp14-ExoN is required for replication fidelity of all coronaviruses, nor the impact of decreased fidelity on genome diversity and fitness during replication and passage. We report here the engineering and recovery of nsp14-ExoN mutant viruses of severe acute respiratory syndrome coronavirus (SARS-CoV) that have stable growth defects and demonstrate a 21-fold increase in mutation frequency during replication in culture. Analysis of complete genome sequences from SARS-ExoN mutant viral clones revealed unique mutation sets in every genome examined from the same round of replication and a total of 100 unique mutations across the genome. Using novel bioinformatic tools and deep sequencing across the full-length genome following 10 population passages in vitro, we demonstrate retention of ExoN mutations and continued increased diversity and mutational load compared to wild-type SARS-CoV. The results define a novel genetic and bioinformatics model for introduction and identification of multi-allelic mutations in replication competent viruses that will be powerful tools for testing the effects of decreased fidelity and increased quasispecies diversity on viral replication, pathogenesis, and evolution