Positive selection of synonymous mutations in vesicular stomatitis virus

Comparison of Endemic and Epidemic Vesicular Stomatitis Virus Lineages in Culicoides sonorensis Midges

Vesicular stomatitis virus (VSV) primarily infects livestock and is transmitted by direct contact and vectored by Culicoides midges (Diptera: Ceratopogonidae). Endemic to Central and South America, specific VSV lineages spread northward out of endemic regions of Mexico and into the U.S. sporadically every five to ten years. In 2012, a monophyletic epidemic lineage 1.1 successfully spread northward into the U.S. In contrast, the closest endemic ancestor, lineage 1.2, remained circulating exclusively in endemic regions in Mexico. It is not clear what roles virus-animal interactions and/or virus-vector interactions play in the ability of specific viral lineages to escape endemic regions in Mexico and successfully cause outbreaks in the U.S., nor the genetic basis for such incursions. Whole-genome sequencing of epidemic VSV 1.1 and endemic VSV 1.2 revealed significant differences in just seven amino acids. Previous studies in swine showed that VSV 1.1 was more virulent than VSV 1.2. Here, we compared the efficiency of these two viral lineages to infect the vector Culicoides sonorensis (Wirth and Jones) and disseminate to salivary glands for subsequent transmission. Our results showed that midges orally infected with the epidemic VSV 1.1 lineage had significantly higher infection dissemination rates compared to those infected with the endemic VSV 1.2 lineage. Thus, in addition to affecting virus-animal interactions, as seen with higher virulence in pigs, small genetic changes may also affect virus-vector interactions, contributing to the ability of specific viral lineages to escape endemic regions via vector-borne transmission.

Innate immune pathways act synergistically to constrain RNA virus evolution in Drosophila melanogaster

Host-pathogen interactions impose recurrent selective pressures that lead to constant adaptation and counter-adaptation in both competing species. Here, we sought to study this evolutionary arms-race and assessed the impact of the innate immune system on viral population diversity and evolution, using Drosophila melanogaster as model host and its natural pathogen Drosophila C virus (DCV). We isogenized eight fly genotypes generating animals defective for RNAi, Imd and Toll innate immune pathways as well as pathogen sensing and gut renewal pathways. Wild-type or mutant flies were then orally infected with DCV, and the virus was serially passaged ten times via reinfection in naïve flies. Viral population diversity was studied after each viral passage by high-throughput sequencing, and infection phenotypes were assessed at the beginning and at the end of the evolution experiment. We found that the absence of any of the various immune pathways studied increased viral genetic diversity while attenuating virulence. Strikingly, these effects were observed in a range of host factors described as having mainly antiviral or antibacterial functions. Together, our results indicate that the innate immune system as a whole, and not specific antiviral defense pathways in isolation, generally constrains viral diversity and evolution.

Mutation rate of SARS-CoV-2 and emergence of mutators during experimental evolution

Background and objectives

To understand how organisms evolve, it is fundamental to study how mutations emerge and establish. Here, we estimated the rate of mutation accumulation of SARS-CoV-2 in vitro and investigated the repeatability of its evolution when facing a new cell type but no immune or drug pressures.

Methodology

We performed experimental evolution with two strains of SARS-CoV-2, one carrying the originally described spike protein (CoV-2-D) and another carrying the D614G mutation that has spread worldwide (CoV-2-G). After 15 passages in Vero cells and whole genome sequencing, we characterized the spectrum and rate of the emerging mutations and looked for evidences of selection across the genomes of both strains.

Results

From the frequencies of the mutations accumulated, and excluding the genes with signals of selection, we estimate a spontaneous mutation rate of 1.3 × 10 ^-6 ± 0.2 × 10^-6 per-base per-infection cycle (mean across both lineages of SARS-CoV-2 ± 2SEM). We further show that mutation accumulation is larger in the CoV-2-D lineage and heterogeneous along the genome, consistent with the action of positive selection on the spike protein, which accumulated five times more mutations than the corresponding genomic average. We also observe the emergence of mutators in the CoV-2-G background, likely linked to mutations in the RNA-dependent RNA polymerase and/or in the error-correcting exonuclease protein.

Conclusions and implications

These results provide valuable information on how spontaneous mutations emerge in SARS-CoV-2 and on how selection can shape its genome toward adaptation to new environments. Lay Summary: Each time a virus replicates inside a cell, errors (mutations) occur. Here, via laboratory propagation in cells originally isolated from the kidney epithelium of African green monkeys, we estimated the rate at which the SARS-CoV-2 virus mutates-an important parameter for understanding how it can evolve within and across humans. We also confirm the potential of its Spike protein to adapt to a new environment and report the emergence of mutators-viral populations where mutations occur at a significantly faster rate.

SARS-CoV-2 genetic variations associated with COVID-19 pathogenicity

In this study, we performed genome-wide association analyses on SARS-CoV-2 genomes to identify genetic mutations associated with pre-symptomatic/asymptomatic COVID-19 cases. Various potential covariates and confounding factors of COVID-19 severity, including patient age, gender and country, as well as virus phylogenetic relatedness were adjusted for. In total, 3021 full-length genomes of SARS-CoV-2 generated from original clinical samples and whose patient status could be determined conclusively as either ‘pre-symptomatic/asymptomatic’ or ‘symptomatic’ were retrieved from the GISAID database. We found that the mutation 11 083G>T, located in the coding region of non-structural protein 6, is significantly associated with asymptomatic COVID-19. Patient age is positively correlated with symptomatic infection, while gender is not significantly correlated with the development of the disease. We also found that the effects of the mutation, patient age and gender do not vary significantly among countries, although each country appears to have varying baseline chances of COVID-19 symptom development.

Effects of synonymous mutations beyond codon bias: The evidence for adaptive synonymous substitutions from microbial evolution experiments

Synonymous mutations are often assumed to be neutral with respect to fitness because they do not alter the encoded amino acid and so cannot be 'seen' by natural selection. Yet a growing body of evidence suggests that synonymous mutations can have fitness effects that drive adaptive evolution through their impacts on gene expression and protein folding. Here, we review what microbial experiments have taught us about the contribution of synonymous mutations to adaptation. A survey of site-directed mutagenesis experiments reveals the distributions of fitness effects for nonsynonymous and synonymous mutations are more similar, especially for beneficial mutations, than expected if all synonymous mutations were neutral, suggesting they should drive adaptive evolution more often than is typically observed. A review of experimental evolution studies where synonymous mutations have contributed to adaptation shows they can impact fitness through a range of mechanisms including the creation of illicit RNA polymerase binding sites impacting transcription and changes to mRNA folding stability that modulate translation. We suggest that clonal interference in evolving microbial populations may be the reason synonymous mutations play a smaller role in adaptive evolution than expected based on their observed fitness effects. We finish by discussing the impacts of falsely assuming synonymous mutations are neutral and discuss directions for future work exploring the role of synonymous mutations in adaptive evolution.

Adaptation and Codon-Usage Preference of Apple and Pear-Infecting Apple Stem Grooving Viruses

Apple stem grooving virus (ASGV; genus Capillovirus) is an economically important virus. It has an approx. 6.5 kb, monopartite, linear, positive-sense, single-stranded RNA genome. The present study includes identification of 24 isolates—13 isolates from apple (Pyrus malus L.) and 11 isolates from pear (Pyrus communis L.)—from different agricultural fields in South Korea. The coat protein (CP) gene of the corresponding 23 isolates were amplified, sequenced, and analyzed. The CP sequences showed phylogenetic separation based on their host species, and not on the geography, indicating host adaptation. Further analysis showed that the ASGV isolated in this study followed host adaptation influenced and preferred by the host codon-usage.

Experimental virus evolution in cancer cell monolayers, spheroids, and tissue explants

Viral laboratory evolution has been used for different applications, such as modeling viral emergence, drug-resistance prediction, and therapeutic virus optimization. However, these studies have been mainly performed in cell monolayers, a highly simplified environment, raising concerns about their applicability and relevance. To address this, we compared the evolution of a model virus in monolayers, spheroids, and tissue explants. We performed this analysis in the context of cancer virotherapy by performing serial transfers of an oncolytic vesicular stomatitis virus (VSV-Δ51) in 4T1 mouse mammary tumor cells. We found that VSV-Δ51 gained fitness in each of these three culture systems, and that adaptation to the more complex environments (spheroids or explants) correlated with increased fitness in monolayers. Most evolved lines improved their ability to suppress β-interferon secretion compared to the VSV-Δ51 founder, suggesting that the selective pressure exerted by antiviral innate immunity was important in the three systems. However, system-specific patterns were also found. First, viruses evolved in monolayers remained more oncoselective that those evolved in spheroids, since the latter showed concomitant adaptation to non-tumoral mouse cells. Second, deep sequencing indicated that viral populations evolved in monolayers or explants tended to be more genetically diverse than those evolved in spheroids. Finally, we found highly variable outcomes among independent evolutionary lines propagated in explants. We conclude that experimental evolution in monolayers tends to be more reproducible than in spheroids or explants, and better preserves oncoselectivity. Our results also suggest that monolayers capture at least some relevant selective pressures present in more complex systems.

Dengue Virus Serotype 2 Intrahost Diversity in Patients with Different Clinical Outcomes

Intrahost genetic diversity is thought to facilitate arbovirus adaptation to changing environments and hosts, and it might also be linked to viral pathogenesis. Dengue virus serotype 2 (DENV-2) has circulated in Brazil since 1990 and is associated with severe disease and explosive outbreaks. Intending to shed light on the viral determinants for severe dengue pathogenesis, we sought to analyze the DENV-2 intrahost genetic diversity in 68 patient cases clinically classified as dengue fever (n = 31), dengue with warning signs (n = 19), and severe dengue (n = 18). Unlike previous DENV intrahost diversity studies whose approaches employed PCR, here we performed viral whole-genome deep sequencing from clinical samples with an amplicon-free approach, representing the real intrahost diversity scenario. Striking differences were detected in the viral population structure between the three clinical categories, which appear to be driven mainly by different infection times and selection pressures, rather than being linked with the clinical outcome itself. Diversity in the NS2B gene, however, showed to be constrained, irrespective of clinical outcome and infection time. Finally, 385 non-synonymous intrahost single-nucleotide variants located along the viral polyprotein, plus variants located in the untranslated regions, were consistently identified among the samples. Of them, 124 were exclusively or highly detected among cases with warning signs and among severe cases. However, there was no variant that by itself appeared to characterize the cases of greater severity, either due to its low intrahost frequency or the conservative effect on amino acid substitution. Although further studies are necessary to determine their real effect on viral proteins, this heightens the possibility of epistatic interactions. The present analysis represents an initial effort to correlate DENV-2 genetic diversity to its pathogenic potential and thus contribute to understanding the virus’s dynamics within its human host.

Identifying Candidate Genetic Markers of CDV Cross-Species Pathogenicity in African Lions

Canine distemper virus (CDV) is a multi-host pathogen with variable clinical outcomes of infection across and within species. We used whole-genome sequencing (WGS) to search for viral markers correlated with clinical distemper in African lions. To identify candidate markers, we first documented single-nucleotide polymorphisms (SNPs) differentiating CDV strains associated with different clinical outcomes in lions in East Africa. We then conducted evolutionary analyses on WGS from all global CDV lineages to identify loci subject to selection. SNPs that both differentiated East African strains and were under selection were mapped to a phylogenetic tree representing global CDV diversity to assess if candidate markers correlated with documented outbreaks of clinical distemper in lions (n = 3). Of 54 SNPs differentiating East African strains, ten were under positive or episodic diversifying selection and 20 occurred in the clinical strain despite strong purifying selection at those loci. Candidate markers were in functional domains of the RNP complex (n = 19), the matrix protein (n = 4), on CDV glycoproteins (n = 5), and on the V protein (n = 1). We found mutations at two loci in common between sequences from three CDV outbreaks of clinical distemper in African lions; one in the signaling lymphocytic activation molecule receptor (SLAM)-binding region of the hemagglutinin protein and another in the catalytic center of phosphodiester bond formation on the large polymerase protein. These results suggest convergent evolution at these sites may have a functional role in clinical distemper outbreaks in African lions and uncover potential novel barriers to pathogenicity in this species.

Positive Selection of ORF1ab, ORF3a, and ORF8 Genes Drives the Early Evolutionary Trends of SARS-CoV-2 During the 2020 COVID-19 Pandemic

In this study, we analyzed full-length SARS-CoV-2 genomes from multiple countries to determine early trends in the evolutionary dynamics of the novel COVID-19 pandemic. Results indicated SARS-CoV-2 evolved early into at least three phylogenetic groups, characterized by positive selection at specific residues of the accessory proteins ORF3a and ORF8. Also, we are reporting potential relevant sites under positive selection at specific sites of non-structural proteins nsp6 and helicase. Our analysis of co-evolution showed evidence of epistatic interactions among sites in the genome that may be important in the generation of variants adapted to humans. These observations might impact not only public health but also suggest that more studies are needed to understand the genetic mechanisms that may affect the development of therapeutic and preventive tools, like antivirals and vaccines. Collectively, our results highlight the identification of ongoing selection even in a scenario of conserved sequences collected over the first 3 months of this pandemic.

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.

Molecular basis of genetic variation of viruses

Genetic variation is a necessity of all biological systems. Viruses use all known mechanisms of variation; mutation, several forms of recombination, and segment reassortment in the case of viruses with a segmented genome. These processes are intimately connected with the replicative machineries of viruses, as well as with fundamental physical-chemical properties of nucleotides when acting as template or substrate residues. Recombination has been viewed as a means to rescue viable genomes from unfit parents or to produce large modifications for the exploration of phenotypic novelty. All types of genetic variation can act conjointly as blind processes to provide the raw materials for adaptation to the changing environments in which viruses must replicate. A distinction is made between mechanistically unavoidable and evolutionarily relevant mutation and recombination.

Increased Virulence of an Epidemic Strain of Vesicular Stomatitis Virus Is Associated With Interference of the Innate Response in Pigs

Vesicular stomatitis virus (VSV) causes sporadic outbreaks of vesicular disease in the southwestern United States. The intrinsic characteristics of epidemic strains associated with these outbreaks are poorly understood. In this study, we report the distinctive genomic and biological characteristics of an epidemic (NJ0612NME6) strain of VSV compared with an endemic (NJ0806VCB) strain. Genomic comparisons between the two strains revealed a total of 111 nucleotide differences (23 non-synonymous) with potentially relevant replacements located in the P, G, and L proteins. When tested in experimentally infected pigs, a natural host of VSV, the epidemic strain caused higher fever and an increased number of vesicular lesions compared to pigs infected with the endemic strain. Pigs infected with the epidemic strain showed decreased systemic antiviral activity (type I - IFN), lower antibody levels, higher levels of interleukin 6, and lower levels of tumor necrosis factor during the acute phase of disease compared to pigs infected with the endemic strain. Furthermore, we document the existence of an RNAemia phase in pigs experimentally infected with VSV and explored the cause for the lack of recovery of infectious virus from blood. Finally, the epidemic strain was shown to be more efficient in down-regulating transcription of IRF-7 in primary porcine macrophages. Collectively, the data shows that the epidemic strain of VSV we tested has an enhanced ability to modulate the innate immune response of the vertebrate host. Further studies are needed to examine other epidemic strains and what contributions a phenotype of increased virulence might have on the transmission of VSV during epizootics.

Introduction of Viral Hemorrhagic Septicemia Virus into Freshwater Cultured Rainbow Trout Is Followed by Bursts of Adaptive Evolution

Viral hemorrhagic septicemia virus (VHSV), a rhabdovirus infecting teleost fish, has repeatedly crossed the boundary from marine fish species to freshwater cultured rainbow trout. These naturally replicated cross-species transmission events permit the study of general and repeatable evolutionary events occurring in connection with viral emergence in a novel host species. The purpose of the present study was to investigate the adaptive molecular evolution of the VHSV glycoprotein, one of the key virus proteins involved in viral emergence, following emergence from marine species into freshwater cultured rainbow trout. A comprehensive phylogenetic reconstruction of the complete coding region of the VHSV glycoprotein was conducted, and adaptive molecular evolution was investigated using a maximum likelihood approach to compare different codon substitution models allowing for heterogeneous substitution rate ratios among amino acid sites. Evidence of positive selection was detected at six amino acid sites of the VHSV glycoprotein, within the signal peptide, the confirmation-dependent major neutralizing epitope, and the intracellular tail. Evidence of positive selection was found exclusively in rainbow trout-adapted virus isolates, and amino acid combinations found at the six sites under positive selection pressure differentiated rainbow trout- from non-rainbow trout-adapted isolates. Furthermore, four adaptive sites revealed signs of recurring identical changes across phylogenetic groups of rainbow trout-adapted isolates, suggesting that repeated VHSV emergence in freshwater cultured rainbow trout was established through convergent routes of evolution that are associated with immune escape.IMPORTANCE This study is the first to demonstrate that VHSV emergence from marine species into freshwater cultured rainbow trout has been accompanied by bursts of adaptive evolution in the VHSV glycoprotein. Furthermore, repeated detection of the same adaptive amino acid sites across phylogenetic groups of rainbow trout-adapted isolates indicates that adaptation to rainbow trout was established through parallel evolution. In addition, signals of convergent evolution toward the maintenance of genetic variation were detected in the conformation-dependent neutralizing epitope or in close proximity to disulfide bonds involved in the structural conformation of the neutralizing epitope, indicating adaptation to immune response-related genetic variation across freshwater cultured rainbow trout.

The Evolution of Molecular Compatibility between Bacteriophage ΦX174 and its Host

Viruses rely upon their hosts for biosynthesis of viral RNA, DNA and protein. This dependency frequently engenders strong selection for virus genome compatibility with potential hosts, appropriate gene regulation and expression necessary for a successful infection. While bioinformatic studies have shown strong correlations between codon usage in viral and host genomes, the selective factors by which this compatibility evolves remain a matter of conjecture. Engineered to include codons with a lesser usage and/or tRNA abundance within the host, three different attenuated strains of the bacterial virus ФX174 were created and propagated via serial transfers. Molecular sequence data indicate that biosynthetic compatibility was recovered rapidly. Extensive computational simulations were performed to assess the role of mutational biases as well as selection for translational efficiency in the engineered phage. Using bacteriophage as a model system, we can begin to unravel the evolutionary processes shaping codon compatibility between viruses and their host.

Significant differences in terms of codon usage bias between bacteriophage early and late genes: a comparative genomics analysis

Background

Viruses undergo extensive evolutionary selection for efficient replication which effects, among others, their codon distribution. In the current study, we aimed at understanding the way evolution shapes the codon distribution in early vs. late viral genes in terms of their expression during different stages in the viral replication cycle. To this end we analyzed 14 bacteriophages and 11 human viruses with available information about the expression phases of their genes.

Results

We demonstrated evidence of selection for distinct composition of synonymous codons in early and late viral genes in 50% of the analyzed bacteriophages. Among others, this phenomenon may be related to the time specific adaptation of the viral genes to the translation efficiency factors involved at different bacteriophage developmental stages. Specifically, we showed that the differences in codon composition in different temporal gene groups cannot be explained only by phylogenetic proximities between the analyzed bacteriophages, and can be partially explained by differences in the adaptation to the host tRNA pool, nucleotide bias, GC content and more.

In contrast, no difference in temporal regulation of synonymous codon usage was observed in human viruses, possibly because of a stronger selection pressure due to a larger effective population size in bacteriophages and their bacterial hosts.

Conclusions

The codon distribution in large fractions of bacteriophage genomes tend to be different in early and late genes. This phenomenon seems to be related to various aspects of the viral life cycle, and to various intracellular processes. We believe that the reported results should contribute towards better understanding of viral evolution and may promote the development of relevant procedures in synthetic virology.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4248-7) contains supplementary material, which is available to authorized users.

Human Rhinovirus Diversity and Evolution: How Strange the Change from Major to Minor

Rhinoviruses are the most common causes of the common cold. Their many distinct lineages fall into "major" and "minor" groups that use different cell surface receptors to enter host cells. Minor-group rhinoviruses are more immunogenic in laboratory studies, although their patterns of transmission and their cold symptoms are broadly similar to those of the major group. Here we present evolutionary evidence that minor-group viruses are also more immunogenic in humans. A key finding is that rates of amino acid substitutions at exposed sites in the capsid proteins VP2, VP3, and VP1 tend to be elevated in minor-group relative to major-group viruses, while rates at buried sites show no consistent differences. A reanalysis of historical virus watch data also indicates a higher immunogenicity of minor-group viruses, consistent with our findings about evolutionary rates at amino acid positions most directly exposed to immune surveillance. The increased immunogenicity and speed of evolution in minor-group lineages may contribute to the very large numbers of rhinovirus serotypes that coexist while differing in virulence.IMPORTANCE Most colds are caused by rhinoviruses (RVs). Those caused by a subset known as the minor-group members of rhinovirus species A (RV-A) are correlated with the inception and aggravation of asthma in at-risk populations. Genetically, minor-group viruses are similar to major-group RV-A, from which they were derived, although they tend to elicit stronger immune responses. Differences in their rates and patterns of molecular evolution should be highly relevant to their epidemiology. All RV-A strains show high rates of amino acid substitutions in the capsid proteins at exposed sites not previously identified as being immunogenic, and this increase is significantly greater in minor-group viruses. These findings will inform future studies of the recently discovered RV-C, which also appears to exacerbate asthma in adults and children. In addition, these findings draw attention to the difficult problem of explaining the long-term coexistence of many serotypes of major- and minor-group RVs.

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.

Increasing Clinical Severity during a Dengue Virus Type 3 Cuban Epidemic: Deep Sequencing of Evolving Viral Populations

During the dengue virus type 3 (DENV-3) epidemic that occurred in Havana in 2001 to 2002, severe disease was associated with the infection sequence DENV-1 followed by DENV-3 (DENV-1/DENV-3), while the sequence DENV-2/DENV-3 was associated with mild/asymptomatic infections. To determine the role of the virus in the increasing severity demonstrated during the epidemic, serum samples collected at different time points were studied. A total of 22 full-length sequences were obtained using a deep-sequencing approach. Bayesian phylogenetic analysis of consensus sequences revealed that two DENV-3 lineages were circulating in Havana at that time, both grouped within genotype III. The predominant lineage is closely related to Peruvian and Ecuadorian strains, while the minor lineage is related to Venezuelan strains. According to consensus sequences, relatively few nonsynonymous mutations were observed; only one was fixed during the epidemic at position 4380 in the NS2B gene. Intrahost genetic analysis indicated that a significant minor population was selected and became predominant toward the end of the epidemic. In conclusion, greater variability was detected during the epidemic's progression in terms of significant minority variants, particularly in the nonstructural genes. An increasing trend of genetic diversity toward the end of the epidemic was observed only for synonymous variant allele rates, with higher variability in secondary cases. Remarkably, significant intrahost genetic variation was demonstrated within the same patient during the course of secondary infection with DENV-1/DENV-3, including changes in the structural proteins premembrane (PrM) and envelope (E). Therefore, the dynamic of evolving viral populations in the context of heterotypic antibodies could be related to the increasing clinical severity observed during the epidemic.

IMPORTANCE

Based on the evidence that DENV fitness is context dependent, our research has focused on the study of viral factors associated with intraepidemic increasing severity in a unique epidemiological setting. Here, we investigated the intrahost genetic diversity in acute human samples collected at different time points during the DENV-3 epidemic that occurred in Cuba in 2001 to 2002 using a deep-sequencing approach. We concluded that greater variability in significant minor populations occurred as the epidemic progressed, particularly in the nonstructural genes, with higher variability observed in secondary infection cases. Remarkably, for the first time significant intrahost genetic variation was demonstrated within the same patient during the course of secondary infection with DENV-1/DENV-3, including changes in structural proteins. These findings indicate that high-resolution approaches are needed to unravel molecular mechanisms involved in dengue pathogenesis.

Large-Effect Beneficial Synonymous Mutations Mediate Rapid and Parallel Adaptation in a Bacterium

Contrary to previous understanding, recent evidence indicates that synonymous codon changes may sometimes face strong selection. However, it remains difficult to generalize the nature, strength, and mechanism(s) of such selection. Previously, we showed that synonymous variants of a key enzyme-coding gene (fae) of Methylobacterium extorquens AM1 decreased enzyme production and reduced fitness dramatically. We now show that during laboratory evolution, these variants rapidly regained fitness via parallel yet variant-specific, highly beneficial point mutations in the N-terminal region of fae. These mutations (including four synonymous mutations) had weak but consistently positive impacts on transcript levels, enzyme production, or enzyme activity. However, none of the proposed mechanisms (including internal ribosome pause sites or mRNA structure) predicted the fitness impact of evolved or additional, engineered point mutations. This study shows that synonymous mutations can be fixed through strong positive selection, but the mechanism for their benefit varies depending on the local sequence context.

Molecular Basis of Genetic Variation of Viruses

Genetic variation is a necessity of all biological systems. Viruses use all known mechanisms of variation: mutation, several forms of recombination, and segment reassortment in the case of viruses with a segmented genome. These processes are intimately connected with the replicative machineries of viruses, as well as with fundamental physico-chemical properties of nucleotides when acting as template or substrate residues. Recombination has been viewed as a means to rescue viable genomes from unfit parents, or to produce large modifications for the exploration of phenotypic novelty. All types of genetic variation can act conjointly as blind processes to provide the raw materials for adaptation to the changing environments in which viruses must replicate.

Experimental evolution of an RNA virus in cells with innate immunity defects

Experimental evolution studies have shown that RNA viruses respond rapidly to directional selection and thus can adapt efficiently to changes in host cell tropism, antiviral drugs, or other imposed selective pressures. However, the evolution of RNA viruses under relaxed selection has been less extensively explored. Here, we evolved vesicular stomatitis virus in mouse embryonic fibroblasts knocked-out for PKR, a protein with a central role in antiviral innate immunity. Vesicular stomatitis virus adapted to PKR-negative mouse embryonic fibroblasts in a gene-specific manner, since the evolved viruses exhibited little or no fitness improvement in PKR-positive cells. Full-length sequencing revealed the presence of multiple parallel nucleotide substitutions arising in independent evolution lines. However, site-directed mutagenesis showed that the effects of these substitutions were not PKR dependent. In contrast, we found evidence for sign epistasis, such that a given substitution which was positively selected was strongly deleterious when tested as a single mutation. Our results suggest that virus evolution in cells with specific innate immunity defects may drive viral specialization. However, this process is not deterministic at the molecular level, probably because the fixation of mutations which are tolerated under a relaxed selection regime is governed mainly by random genetic drift.

Temporal dynamics of intrahost molecular evolution for a plant RNA virus

Populations of plant RNA viruses are highly polymorphic in infected plants, which may allow rapid within-host evolution. To understand tobacco etch potyvirus (TEV) evolution, longitudinal samples from experimentally evolved populations in the natural host tobacco and from the alternative host pepper were phenotypically characterized and genetically analyzed. Temporal and compartmental variabilities of TEV populations were quantified using high throughput Illumina sequencing and population genetic approaches. Of the two viral phenotypic traits measured, virulence increased in the novel host but decreased in the original one, and viral load decreased in both hosts, though to a lesser extent in the novel one. Dynamics of population genetic diversity were also markedly different among hosts. Population heterozygosity increased in the ancestral host, with a dominance of synonymous mutations fixed, whereas it did not change or even decreased in the new host, with an excess of nonsynonymous mutations. All together, these observations suggest that directional selection is the dominant evolutionary force in TEV populations evolving in a novel host whereas either diversifying selection or random genetic drift may play a fundamental role in the natural host. To better understand these evolutionary dynamics, we developed a computer simulation model that incorporates the effects of mutation, selection, and drift. Upon parameterization with empirical data from previous studies, model predictions matched the observed patterns, thus reinforcing our idea that the empirical patterns of mutation accumulation represent adaptive evolution.

Evolution of oncolytic viruses

Owing to their replicative capacity, oncolytic viruses (OVs) can evolve under the action of natural selection. Reversion to virulence and recombination with wild-type strains may compromise OV safety, therefore requiring evolutionary risk assessment studies. On the other hand, evolution can be directed in the laboratory to create more potent and safer OVs. Previous work in the experimental evolution field provides a background for OV directed evolution, and has identified interesting exploitable features. While genetic engineering has greatly advanced the field of oncolytic virotherapy, this approach is sometimes curtailed by the complexity and diversity of virus-host interactions. Directed evolution provides an alternative approach that may help to obtain new OVs without prejudice toward the underlying molecular mechanisms involved.

Evidence for stabilizing selection on codon usage in chromosomal rearrangements of Drosophila pseudoobscura

There has been a renewed interest in investigating the role of stabilizing selection acting on genome-wide traits such as codon usage bias. Codon bias, when synonymous codons are used at unequal frequencies, occurs in a wide variety of taxa. Standard evolutionary models explain the maintenance of codon bias through a balance of genetic drift, mutation and weak purifying selection. The efficacy of selection is expected to be reduced in regions of suppressed recombination. Contrary to observations in Drosophila melanogaster, some recent studies have failed to detect a relationship between the recombination rate, intensity of selection acting at synonymous sites, and the magnitude of codon bias as predicted under these standard models. Here, we examined codon bias in 2798 protein coding loci on the third chromosome of D. pseudoobscura using whole-genome sequences of 47 individuals, representing five common third chromosome gene arrangements. Fine-scale recombination maps were constructed using more than 1 million segregating sites. As expected, recombination was demonstrated to be significantly suppressed between chromosome arrangements, allowing for a direct examination of the relationship between recombination, selection, and codon bias. As with other Drosophila species, we observe a strong mutational bias away from the most frequently used codons. We find the rate of synonymous and nonsynonymous polymorphism is variable between different amino acids. However, we do not observe a reduction in codon bias or the strength of selection in regions of suppressed recombination as expected. Instead, we find that the interaction between weak stabilizing selection and mutational bias likely plays a role in shaping the composition of synonymous codons across the third chromosome in D. pseudoobscura.

Contribution of silent mutations to thermal adaptation of RNA bacteriophage Qβ

Changes in protein function and other biological properties, such as RNA structure, are crucial for adaptation of organisms to novel or inhibitory environments. To investigate how mutations that do not alter amino acid sequence may be positively selected, we performed a thermal adaptation experiment using the single-stranded RNA bacteriophage Qβ in which the culture temperature was increased from 37.2°C to 41.2°C and finally to an inhibitory temperature of 43.6°C in a stepwise manner in three independent lines. Whole-genome analysis revealed 31 mutations, including 14 mutations that did not result in amino acid sequence alterations, in this thermal adaptation. Eight of the 31 mutations were observed in all three lines. Reconstruction and fitness analyses of Qβ strains containing only mutations observed in all three lines indicated that five mutations that did not result in amino acid sequence changes but increased the amplification ratio appeared in the course of adaptation to growth at 41.2°C. Moreover, these mutations provided a suitable genetic background for subsequent mutations, altering the fitness contribution from deleterious to beneficial. These results clearly showed that mutations that do not alter the amino acid sequence play important roles in adaptation of this single-stranded RNA virus to elevated temperature.

IMPORTANCE

Recent studies using whole-genome analysis technology suggested the importance of mutations that do not alter the amino acid sequence for adaptation of organisms to novel environmental conditions. It is necessary to investigate how these mutations may be positively selected and to determine to what degree such mutations that do not alter amino acid sequences contribute to adaptive evolution. Here, we report the roles of these silent mutations in thermal adaptation of RNA bacteriophage Qβ based on experimental evolution during which Qβ showed adaptation to growth at an inhibitory temperature. Intriguingly, four synonymous mutations and one mutation in the untranslated region that spread widely in the Qβ population during the adaptation process at moderately high temperature provided a suitable genetic background to alter the fitness contribution of subsequent mutations from deleterious to beneficial at a higher temperature.

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.

Experimental evolution of an emerging plant virus in host genotypes that differ in their susceptibility to infection

This study evaluates the extent to which genetic differences among host individuals from the same species condition the evolution of a plant RNA virus. We performed a threefold replicated evolution experiment in which Tobacco etch potyvirus isolate At17b (TEV-At17b), adapted to Arabidopsis thaliana ecotype Ler-0, was serially passaged in five genetically heterogeneous ecotypes of A. thaliana. After 15 passages we found that evolved viruses improved their fitness, showed higher infectivity and stronger virulence in their local host ecotypes. The genome of evolved lineages was sequenced and putative adaptive mutations identified. Host-driven convergent mutations have been identified. Evidences supported selection for increased translational efficiency. Next, we sought for the specificity of virus adaptation by infecting all five ecotypes with all 15 evolved virus populations. We found that some ecotypes were more permissive to infection than others, and that some evolved virus isolates were more specialist/generalist than others. The bipartite network linking ecotypes with evolved viruses was significantly nested but not modular, suggesting that hard-to-infect ecotypes were infected by generalist viruses whereas easy-to-infect ecotypes were infected by all viruses, as predicted by a gene-for-gene model of infection.

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.

Molecular signature of high yield (growth) influenza a virus reassortants prepared as candidate vaccine seeds

Background

Human influenza virus isolates generally grow poorly in embryonated chicken eggs. Hence, gene reassortment of influenza A wild type (wt) viruses is performed with a highly egg adapted donor virus, A/Puerto Rico/8/1934 (PR8), to provide the high yield reassortant (HYR) viral ‘seeds’ for vaccine production. HYR must contain the hemagglutinin (HA) and neuraminidase (NA) genes of wt virus and one to six ‘internal’ genes from PR8. Most studies of influenza wt and HYRs have focused on the HA gene. The main objective of this study is the identification of the molecular signature in all eight gene segments of influenza A HYR candidate vaccine seeds associated with high growth in ovo.

Methodology

The genomes of 14 wt parental viruses, 23 HYRs (5 H1N1; 2, 1976 H1N1-SOIV; 2, 2009 H1N1pdm; 2 H2N2 and 12 H3N2) and PR8 were sequenced using the high-throughput sequencing pipeline with big dye terminator chemistry.

Results

Silent and coding mutations were found in all internal genes derived from PR8 with the exception of the M gene. The M gene derived from PR8 was invariant in all 23 HYRs underlining the critical role of PR8 M in high yield phenotype. None of the wt virus derived internal genes had any silent change(s) except the PB1 gene in X-157. The highest number of recurrent silent and coding mutations was found in NS. With respect to the surface antigens, the majority of HYRs had coding mutations in HA; only 2 HYRs had coding mutations in NA.

Significance

In the era of application of reverse genetics to alter influenza A virus genomes, the mutations identified in the HYR gene segments associated with high growth in ovo may be of great practical benefit to modify PR8 and/or wt virus gene sequences for improved growth of vaccine ‘seed’ viruses.

Comparative genomic and phylogenetic analysis of the first Usutu virus isolate from a human patient presenting with neurological symptoms

Usutu virus (USUV) is a mosquito-borne flavivirus, belonging to the Japanese encephalitis antigenic complex, that circulates among mosquitoes and birds. We describe and analyze the complete genome sequence of the first USUV strain isolated from an immunocompromised patient with neuroinvasive disease. This USUV isolate showed an overall nucleotide identity of 99% and 96%, respectively, with the genomes of isolates from Europe and Africa. Comparison of the human USUV complete polyprotein sequence with bird-derived strains, showed two unique amino acid substitutions. In particular, one substitution (S595G) was situated in the DIII domain of the viral Envelope protein that is recognized by flavivirus neutralizing antibodies. An additional amino acid substitution (D3425E) was identified in the RNA-dependent RNA polymerase (RdRp) domain of the NS5 protein. This substitution is remarkable since E3425 is highly conserved among the other USUV isolates that were not associated with human infection. However, a similar substitution was observed in Japanese encephalitis and in West Nile viruses isolated from humans. Phylogenetic analysis of the human USUV strain revealed a close relationship with an Italian strain isolated in 2009. Analysis of synonymous nucleotide substitutions (SNSs) among the different USUV genomes showed a specific evolutionary divergence among different countries. In addition, 15 SNSs were identified as unique in the human isolate. We also identified four specific nucleotide substitutions in the 5' and 3' untranslated regions (UTRs) in the human isolate that were not present in the other USUV sequences. Our analyses provide the basis for further experimental studies aimed at defining the effective role of these mutations in the USUV genome, their potential role in the development of viral variants pathogenic for humans and their evolution and dispersal out of Africa.

Human microRNA hsa-miR-296-5p suppresses enterovirus 71 replication by targeting the viral genome

Enterovirus 71 (EV71) has emerged as a major cause of neurological disease following the near eradication of poliovirus. Accumulating evidence suggests that mammalian microRNAs (miRNAs), a class of noncoding RNAs of 18 to 23 nucleotides (nt) with important regulatory roles in many cellular processes, participate in host antiviral defenses. However, the roles of miRNAs in EV71 infection and pathogenesis are still unclear. Here, hsa-miR-296-5p expression was significantly increased in EV71-infected human cells. As determined by virus titration, quantitative real-time PCR (qRT-PCR), and Western blotting, overexpression of hsa-miR-296-5p inhibited, while inhibition of endogenous hsa-miR-296-5p facilitated, EV71 infection. Additionally, two potential hsa-miR-296-5p targets (nt 2115 to 2135 and nt 2896 to 2920) located in the EV71 genome (strain BrCr) were bioinformatically predicted and validated by luciferase reporter assays and Western blotting. Genomic alignment of various EV71 strains revealed synonymous mutations in hsa-miR-296-5p target sequences. Furthermore, the introduction of synonymous mutations into the EV71 BrCr genome by site-directed mutagenesis impaired the viral inhibitory effects of hsa-miR-296-5p and facilitated mutant virus infection. Meanwhile, compensatory mutations in corresponding hsa-miR-296-5p target sequences of the EV71 HeN strain (GenBank accession number JN256064) restored the inhibitory effects of the miRNA. These results indicate that hsa-miR-296-5p inhibits EV71 replication by targeting the viral genome. Our findings support the notion that cellular miRNAs can inhibit virus infection and that the virus mutates to escape suppression by cellular miRNAs.

Tracking viral evolution during a disease outbreak: the rapid and complete selective sweep of a circovirus in the endangered Echo parakeet

Circoviruses are among the smallest and simplest of all viruses, but they are relatively poorly characterized. Here, we intensively sampled two sympatric parrot populations from Mauritius over a period of 11 years and screened for the circovirus Beak and feather disease virus (BFDV). During the sampling period, a severe outbreak of psittacine beak and feather disease, which is caused by BFDV, occurred in Echo parakeets. Consequently, this data set presents an ideal system for studying the evolution of a pathogen in a natural population and to understand the adaptive changes that cause outbreaks. Unexpectedly, we discovered that the outbreak was most likely caused by changes in functionally important regions of the normally conserved replication-associated protein gene and not the immunogenic capsid. Moreover, these mutations were completely fixed in the Echo parakeet host population very shortly after the outbreak. Several capsid alleles were linked to the replication-associated protein outbreak allele, suggesting that whereas the key changes occurred in the latter, the scope of the outbreak and the selective sweep may have been influenced by positive selection in the capsid. We found evidence for viral transmission between the two host populations though evidence for the invasive species as the source of the outbreak was equivocal. Finally, the high evolutionary rate that we estimated shows how rapidly new variation can arise in BFDV and is consistent with recent results from other small single-stranded DNA viruses.

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.

POSITIVE DARWINIAN SELECTION: DOES THE COMPARATIVE METHOD RULE?

To detect positive Darwinian selection it is thought essential to compare two sequences. Despite its defects, "the comparative method rules." However, genes evolving rapidly under positive selection conflict more with internal forces (the genome phenotype) than genes evolving slowly under negative selection. In particular, there is conflict with stem-loop potential. The conflict between protein-encoding potential (primary information) and stem-loop potential (secondary information) permits detection of positive selection in a single sequence. The degree to which secondary information is compromised provides a measure of the speed of transmission of primary information. Thus, the sovereignty of the comparative method is challenged not only by its own defects, but also by the availability of a single-sequence method. However, while of limited utility for positive selection, the comparative method casts new light on Darwin's great question — the origin of species. Comparison of rates of synonymous and non-synonymous mutation suggests that branching into new species begins with synonymous mutations.

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.

The fitness effects of synonymous mutations in DNA and RNA viruses

Despite being silent with respect to protein sequence, synonymous nucleotide substitutions can be targeted by natural selection directly at the DNA or RNA level. However, there has been no systematic assessment of how frequent this type of selection is. Here, we have constructed 53 single random synonymous substitution mutants of the bacteriophages Qβ and ΦX174 by site-directed mutagenesis and assayed their fitness. Analysis of this mutant collection and of previous studies undertaken with a variety of single-stranded (ss) viruses demonstrates that selection at synonymous sites is stronger in RNA viruses than in DNA viruses. We estimate that this type of selection contributes approximately 18% of the overall mutational fitness effects in ssRNA viruses under our assay conditions and that random synonymous substitutions have a 5% chance of being lethal to the virus, whereas in ssDNA viruses, these figures drop to 1.4% and 0%, respectively. In contrast, the effects of nonsynonymous substitutions appear to be similar in ssRNA and ssDNA viruses.

Identifying genetic markers of adaptation for surveillance of viral host jumps

Adaptation is often thought to affect the likelihood that a virus will be able to successfully emerge in a new host species. If so, surveillance for genetic markers of adaptation could help to predict the risk of disease emergence. However, adaptation is difficult to distinguish conclusively from the other processes that generate genetic change. In this Review we survey the research on the host jumps of influenza A, severe acute respiratory syndrome-coronavirus, canine parvovirus and Venezuelan equine encephalitis virus to illustrate the insights that can arise from combining genetic surveillance with microbiological experimentation in the context of epidemiological data. We argue that using a multidisciplinary approach for surveillance will provide a better understanding of when adaptations are required for host jumps and thus when predictive genetic markers may be present.

Adaptive evolutionary walks require neutral intermediates in RNA fitness landscapes

In RNA fitness landscapes with interconnected networks of neutral mutations, neutral precursor mutations can play an important role in facilitating the accessibility of epistatic adaptive mutant combinations. I use an exhaustively surveyed fitness landscape model based on short sequence RNA genotypes (and their secondary structure phenotypes) to calculate the minimum rate at which mutants initially appearing as neutral are incorporated into an adaptive evolutionary walk. I show first, that incorporating neutral mutations significantly increases the number of point mutations in a given evolutionary walk when compared to estimates from previous adaptive walk models. Second, that incorporating neutral mutants into such a walk significantly increases the final fitness encountered on that walk - indeed evolutionary walks including neutral steps often reach the global optimum in this model. Third, and perhaps most importantly, evolutionary paths of this kind are often extremely winding in their nature and have the potential to undergo multiple mutations at a given sequence position within a single walk; the potential of these winding paths to mislead phylogenetic reconstruction is briefly considered.

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.

Accumulation of point mutations and reassortment of genomic RNA segments are involved in the microevolution of Puumala hantavirus in a bank vole (Myodes glareolus) population

The genetic diversity of Puumala hantavirus (PUUV) was studied in a local population of its natural host, the bank vole (Myodes glareolus). The trapping area (2.5 x 2.5 km) at Konnevesi, Central Finland, included 14 trapping sites, at least 500 m apart; altogether, 147 voles were captured during May and October 2005. Partial sequences of the S, M and L viral genome segments were recovered from 40 animals. Seven, 12 and 17 variants were detected for the S, M and L sequences, respectively; these represent new wild-type PUUV strains that belong to the Finnish genetic lineage. The genetic diversity of PUUV strains from Konnevesi was 0.2-4.9 % for the S segment, 0.2-4.8 % for the M segment and 0.2-9.7 % for the L segment. Most nucleotide substitutions were synonymous and most deduced amino acid substitutions were conservative, probably due to strong stabilizing selection operating at the protein level. Based on both sequence markers and phylogenetic clustering, the S, M and L sequences could be assigned to two groups, 'A' and 'B'. Notably, not all bank voles carried S, M and L sequences belonging to the same group, i.e. S(A)M(A)L(A) or S(B)M(B)L(B). A substantial proportion (8/40, 20 %) of the newly characterized PUUV strains possessed reassortant genomes such as S(B)M(A)L(A), S(A)M(B)L(B) or S(B)M(A)L(B). These results suggest that at least some of the PUUV reassortants are viable and can survive in the presence of their parental strains.

Antagonistic pleiotropy involving promoter sequences in a virus

Selection of specialist genotypes, that is, populations with limited niche width, promotes the maintenance of diversity. Specialization to a particular environment may have a cost in other environments, including fitness tradeoffs. When the tradeoffs are the result of mutations that have a beneficial effect in the selective environment but a deleterious effect in other environments, we have antagonistic pleiotropy. Alternatively, tradeoffs can result from the fixation of mutations that are neutral in the selective environment but have a negative effect in other environments, and thus the tradeoff is due to mutation accumulation. We tested the mechanisms underlying the fitness tradeoffs observed during adaptation to persistent infection of vesicular stomatitis virus in insect cells by sequencing the full-length genomes of 12 strains with a history of replication in a single niche (acute mammalian infection or persistent insect infection) or in temporally heterogeneous niches and correlated genetic and fitness changes. Ecological theory predicts a correlation between the selective environment and the niche width of the evolved populations, such that adaptation to single niches should lead to the selection of specialists and niche cycling should result in the selection of generalists. Contrary to this expectation, adaptation to one of the single niches resulted in a generalist and adaptation to a heterogeneous environment led to the selection of a specialist. Only one-third of the mutations that accumulated during persistent infection had a fitness cost that could be explained in all cases by antagonistic pleiotropy. Mutations involved in fitness tradeoffs included changes in regulatory sequences, particularly at the 3' termini of the genomes, which contain the single promoter that controls viral transcription and replication.

Quasispecies dynamics and fixation of a synonymous mutation in hantavirus transmission

RNA-dependent RNA polymerases, the key enzymes in replication of RNA viruses, have a low fidelity; thus, these viruses replicate as a swarm of mutants termed viral quasispecies. Constant generation of new mutations allows RNA viruses to adapt swiftly to a novel environment through selection of both pre-existing and de novo-generated genetic variants. Here, quasispecies dynamics were studied in vivo in controlled hantavirus transmission from experimentally infected to naïve rodents through infested cage bedding. An elementary step of virus microevolution was apparent, as one synonymous mutation (A759G) repeatedly became fixed in the viral RNA quasispecies populations in the recipient animals.

Rapid adaptive amplification of preexisting variation in an RNA virus

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

Evolutionary genomics of host adaptation in vesicular stomatitis virus

Populations experiencing similar selection pressures can sometimes diverge in the genetic architectures underlying evolved complex traits. We used RNA virus populations of large size and high mutation rate to study the impact of historical environment on genome evolution, thus increasing our ability to detect repeatable patterns in the evolution of genetic architecture. Experimental vesicular stomatitis virus populations were evolved on HeLa cells, on MDCK cells, or on alternating hosts. Turner and Elena (2000. Cost of host radiation in an RNA virus. Genetics. 156:1465-1470.) previously showed that virus populations evolved in single-host environments achieved high fitness on their selected hosts but failed to increase in fitness relative to their ancestor on the unselected host and that alternating-host-evolved populations had high fitness on both hosts. Here we determined the complete consensus sequence for each evolved population after 95 generations to gauge whether the parallel phenotypic changes were associated with parallel genomic changes. We also analyzed the patterns of allele substitutions to discern whether differences in fitness across hosts arose through true pleiotropy or the presence of not only a mutation that is beneficial in both hosts but also 1 or more mutations at other loci that are costly in the unselected environment (mutation accumulation [MA]). We found that ecological history may influence to what extent pleiotropy and MA contribute to fitness asymmetries across environments. We discuss the degree to which current genetic architecture is expected to constrain future evolution of complex traits, such as host use by RNA viruses.

Theme and variations in the evolutionary pathways to virulence of an RNA plant virus species

The diversity of a highly variable RNA plant virus was considered to determine the range of virulence substitutions, the evolutionary pathways to virulence, and whether intraspecific diversity modulates virulence pathways and propensity. In all, 114 isolates representative of the genetic and geographic diversity of Rice yellow mottle virus (RYMV) in Africa were inoculated to several cultivars with eIF(iso)4G-mediated Rymv1-2 resistance. Altogether, 41 virulent variants generated from ten wild isolates were analyzed. Nonconservative amino acid replacements at five positions located within a stretch of 15 codons in the central region of the 79-aa-long protein VPg were associated with virulence. Virulence substitutions were fixed predominantly at codon 48 in most strains, whatever the host genetic background or the experimental conditions. There were one major and two isolate-specific mutational pathways conferring virulence at codon 48. In the prevalent mutational pathway I, arginine (AGA) was successively displaced by glycine (GGA) and glutamic acid (GAA). Substitutions in the other virulence codons were displaced when E48 was fixed. In the isolate-specific mutational pathway II, isoleucine (ATA) emerged and often later coexisted with valine (GTA). In mutational pathway III, arginine, with the specific S2/S3 strain codon usage AGG, was displaced by tryptophane (TGG). Mutational pathway I never arose in the widely spread West African S2/S3 strain because G48 was not infectious in the S2/S3 genetic context. Strain S2/S3 least frequently overcame resistance, whereas two geographically localized variants of the strain S4 had a high propensity to virulence. Codons 49 and 26 of the VPg, under diversifying selection, are candidate positions in modulating the genetic barriers to virulence. The theme and variations in the evolutionary pathways to virulence of RYMV illustrates the extent of parallel evolution within a highly variable RNA plant virus species.

Quantifying antagonistic epistasis in a multifunctional RNA secondary structure of the Rous sarcoma virus

Recent studies have suggested that antagonistic epistasis (i.e. mutations having smaller effects in combination than alone) may be common among RNA viruses, in contrast to other biological systems. Here, by re-analysing previously published data from a random viral library, selection and epistasis coefficients were estimated in the U5-IR stem and loop of theRous sarcoma virus, a region that adopts a conserved secondary structure and is involved in various essential steps of viral infection. The estimated mutational fitness effects are extremely high and genetic interactions are antagonistic on average. This pattern might be representative of RNA virus genomes, which show high compaction and frequent secondary structures. The implications for RNA virus adaptability are explored.