The rate and spectrum of spontaneous mutations in a plant RNA virus

High-resolution mapping reveals the mechanism and contribution of genome insertions and deletions to RNA virus evolution

RNA viruses rapidly adapt to selective conditions due to the high intrinsic mutation rates of their RNA-dependent RNA polymerases (RdRps). Insertions and deletions (indels) in viral genomes are major contributors to both deleterious mutational load and evolutionary novelty, but remain understudied. To characterize the mechanistic details of their formation and evolutionary dynamics during infection, we developed a hybrid experimental-bioinformatic approach. This approach, called MultiMatch, extracts insertions and deletions from ultradeep sequencing experiments, including those occurring at extremely low frequencies, allowing us to map their genomic distribution and quantify the rates at which they occur. Mapping indel mutations in adapting poliovirus and dengue virus populations, we determine the rates of indel generation and identify mechanistic and functional constraints shaping indel diversity. Using poliovirus RdRp variants of distinct fidelity and genome recombination rates, we demonstrate tradeoffs between fidelity and Indel generation. Additionally, we show that maintaining translation frame and viral RNA structures constrain the Indel landscape and that, due to these significant fitness effects, Indels exert a significant deleterious load on adapting viral populations. Conversely, we uncover positively selected Indels that modulate RNA structure, generate protein variants, and produce defective interfering genomes in viral populations. Together, our analyses establish the kinetic and mechanistic tradeoffs between misincorporation, recombination, and Indel rates and reveal functional principles defining the central role of Indels in virus evolution, emergence, and the regulation of viral infection.

Single-cell mutation rate of turnip crinkle virus (-)-strand replication intermediates

Viruses with single-stranded, positive-sense (+) RNA genomes incur high numbers of errors during replication, thereby creating diversified genome populations from which new, better adapted viral variants can emerge. However, a definitive error rate is known for a relatively few (+) RNA plant viruses, due to challenges to account for perturbations caused by natural selection and/or experimental set-ups. To address these challenges, we developed a new approach that exclusively profiled errors in the (-)-strand replication intermediates of turnip crinkle virus (TCV), in singly infected cells. A series of controls and safeguards were devised to ensure errors inherent to the experimental process were accounted for. This approach permitted the estimation of a TCV error rate of 8.47 X 10-5 substitution per nucleotide site per cell infection. Importantly, the characteristic error distribution pattern among the 50 copies of 2,363-base-pair cDNA fragments predicted that nearly all TCV (-) strands were products of one replication cycle per cell. Furthermore, some of the errors probably elevated error frequencies by lowering the fidelity of TCV RNA-dependent RNA polymerase, and/or permitting occasional re-replication of progeny genomes. In summary, by profiling errors in TCV (-)-strand intermediates incurred during replication in single cells, this study provided strong support for a stamping machine mode of replication employed by a (+) RNA virus.

Mutation rate, selection, and epistasis inferred from RNA virus haplotypes via neural posterior estimation

RNA viruses are particularly notorious for their high levels of genetic diversity, which is generated through the forces of mutation and natural selection. However, disentangling these two forces is a considerable challenge, and this may lead to widely divergent estimates of viral mutation rates, as well as difficulties in inferring the fitness effects of mutations. Here, we develop, test, and apply an approach aimed at inferring the mutation rate and key parameters that govern natural selection, from haplotype sequences covering full-length genomes of an evolving virus population. Our approach employs neural posterior estimation, a computational technique that applies simulation-based inference with neural networks to jointly infer multiple model parameters. We first tested our approach on synthetic data simulated using different mutation rates and selection parameters while accounting for sequencing errors. Reassuringly, the inferred parameter estimates were accurate and unbiased. We then applied our approach to haplotype sequencing data from a serial passaging experiment with the MS2 bacteriophage, a virus that parasites Escherichia coli. We estimated that the mutation rate of this phage is around 0.2 mutations per genome per replication cycle (95% highest density interval: 0.051-0.56). We validated this finding with two different approaches based on single-locus models that gave similar estimates but with much broader posterior distributions. Furthermore, we found evidence for reciprocal sign epistasis between four strongly beneficial mutations that all reside in an RNA stem loop that controls the expression of the viral lysis protein, responsible for lysing host cells and viral egress. We surmise that there is a fine balance between over- and underexpression of lysis that leads to this pattern of epistasis. To recap, we have developed an approach for joint inference of the mutation rate and selection parameters from full haplotype data with sequencing errors and used it to reveal features governing MS2 evolution.

Insect allies-Assessment of a viral approach to plant genome editing

The Insect Allies program of the Defense Advanced Research Projects Agency has already sparked scientific debate concerning technology assessment-related issues, among which the most prevalent is that of dual use. Apart from the issues concerning peaceful applications, the technology also provides the blueprint for a potential bioweapon. However, the combination of a virus-induced genetic modification of crop plants in the field using genetically modified insect vectors poses a greater risk than the hitherto existing use of genetically modified organisms. The technology's great depth of intervention allows a number of sources for hazard and a tendency towards high exposure, but it is also encumbered with notable deficits in knowledge. These issues call for a thorough technology assessment. This article aims to provide an initial characterization from a technology assessment perspective, focusing on potential sources of risk for this novel invasive environmental biotechnology at an early stage of research and development. Integr Environ Assess Manag 2022;18:1488-1499. © 2022 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Analysis of the Contribution of Intrinsic Disorder in Shaping Potyvirus Genetic Diversity

Intrinsically disordered regions (IDRs) are abundant in the proteome of RNA viruses. The multifunctional properties of these regions are widely documented and their structural flexibility is associated with the low constraint in their amino acid positions. Therefore, from an evolutionary stand point, these regions could have a greater propensity to accumulate non-synonymous mutations (NS) than highly structured regions (ORs, or 'ordered regions'). To address this hypothesis, we compared the distribution of non-synonymous mutations (NS), which we relate here to mutational robustness, in IDRs and ORs in the genome of potyviruses, a major genus of plant viruses. For this purpose, a simulation model was built and used to distinguish a possible selection phenomenon in the biological datasets from randomly generated mutations. We analyzed several short-term experimental evolution datasets. An analysis was also performed on the natural diversity of three different species of potyviruses reflecting their long-term evolution. We observed that the mutational robustness of IDRs is significantly higher than that of ORs. Moreover, the substitutions in the ORs are very constrained by the conservation of the physico-chemical properties of the amino acids. This feature is not found in the IDRs where the substitutions tend to be more random. This reflects the weak structural constraints in these regions, wherein an amino acid polymorphism is naturally conserved. In the course of evolution, potyvirus IDRs and ORs follow different evolutive paths with respect to their mutational robustness. These results have forced the authors to consider the hypothesis that IDRs and their associated amino acid polymorphism could constitute a potential adaptive reservoir.

Empirical estimates of the mutation rate for an alphabaculovirus

Mutation rates are of key importance for understanding evolutionary processes and predicting their outcomes. Empirical mutation rate estimates are available for a number of RNA viruses, but few are available for DNA viruses, which tend to have larger genomes. Whilst some viruses have very high mutation rates, lower mutation rates are expected for viruses with large genomes to ensure genome integrity. Alphabaculoviruses are insect viruses with large genomes and often have high levels of polymorphism, suggesting high mutation rates despite evidence of proofreading activity by the replication machinery. Here, we report an empirical estimate of the mutation rate per base per strand copying (s/n/r) of Autographa californica multiple nucleopolyhedrovirus (AcMNPV). To avoid biases due to selection, we analyzed mutations that occurred in a stable, non-functional genomic insert after five serial passages in Spodoptera exigua larvae. Our results highlight that viral demography and the stringency of mutation calling affect mutation rate estimates, and that using a population genetic simulation model to make inferences can mitigate the impact of these processes on estimates of mutation rate. We estimated a mutation rate of μ = 1×10⁻⁷ s/n/r when applying the most stringent criteria for mutation calling, and estimates of up to μ = 5×10⁻⁷ s/n/r when relaxing these criteria. The rates at which different classes of mutations accumulate provide good evidence for neutrality of mutations occurring within the inserted region. We therefore present a robust approach for mutation rate estimation for viruses with stable genomes, and strong evidence of a much lower alphabaculovirus mutation rate than supposed based on the high levels of polymorphism observed.

Virus populations can evolve rapidly, driven by the large number of mutations that occur during virus replication. It is challenging to measure mutation rates because selection will affect which mutations are observed: beneficial mutations are overrepresented in virus populations, while deleterious mutations are selected against and therefore underrepresented. Few mutation rates have been estimated for viruses with large DNA genomes, and there are no estimates for any insect virus. Here, we estimate the mutation rate for an alphabaculovirus, a virus that infects caterpillars and has a large, 134 kilobase pair DNA genome. To ensure that selection did not bias our estimate of mutation rate, we studied which mutations occurred in a large artificial region inserted into the virus genome, where mutations did not affect viral fitness. We deep sequenced evolved virus populations, and compared the distribution of observed mutants to predictions from a simulation model to estimate mutation rate. We found evidence for a relatively low mutation rate, of one mutation in every 10 million bases replicated. This estimate is in line with expectations for a DNA virus with self-correcting replication machinery and a large genome.

OUP accepted manuscript

Potyviridae, the largest family of known RNA viruses (realm Riboviria), belongs to the picorna-like supergroup and has important agricultural and ecological impacts. Potyvirid genomes are translated into polyproteins, which are in turn hydrolyzed to release mature products. Recent sequencing efforts revealed an unprecedented number of potyvirids with a rich variability in gene content and genomic layouts. Here, we review the heterogeneity of non-core modules that expand the structural and functional diversity of the potyvirid proteomes. We provide a family-wide classification of P1 proteinases into the functional Types A and B, and discuss pretty interesting sweet potato potyviral ORF (PISPO), putative zinc fingers, and alkylation B (AlkB)—non-core modules found within P1 cistrons. The atypical inosine triphosphate pyrophosphatase (ITPase/HAM1), as well as the pseudo tobacco mosaic virus-like coat protein (TMV-like CP) are discussed alongside homologs of unrelated virus taxa. Family-wide abundance of the multitasking helper component proteinase (HC-pro) is revised. Functional connections between non-core modules are highlighted to support host niche adaptation and immune evasion as main drivers of the Potyviridae evolutionary radiation. Potential biotechnological and synthetic biology applications of potyvirid leader proteinases and non-core modules are finally explored.

RNA Viral Vectors for Accelerating Plant Synthetic Biology

The tools of synthetic biology have enormous potential to help us uncover the fundamental mechanisms controlling development and metabolism in plants. However, their effective utilization typically requires transgenesis, which is plagued by long timescales and high costs. In this review we explore how transgenesis can be minimized by delivering foreign genetic material to plants with systemically mobile and persistent vectors based on RNA viruses. We examine the progress that has been made thus far and highlight the hurdles that need to be overcome and some potential strategies to do so. We conclude with a discussion of biocontainment mechanisms to ensure these vectors can be used safely as well as how these vectors might expand the accessibility of plant synthetic biology techniques. RNA vectors stand poised to revolutionize plant synthetic biology by making genetic manipulation of plants cheaper and easier to deploy, as well as by accelerating experimental timescales from years to weeks.

Loss-of-function of Nicotiana tabacum L. eukaryotic translation initiation factors eIF4E1-S and eIF(iso)4E-T synergistically confers high-level resistance to both Potato virus Y (PVY) and resistance-breaking PVY

The plant eukaryotic translation-initiation factors eIF4E and eIF(iso)4E play key roles in infection by plant RNA viruses, especially potyviruses. Mutations in the genes that encode these factors reduce susceptibility to the viruses. In the amphidiploid plant tobacco (Nicotiana tabacum L.), eIF4E1-S deletion mutants resist Potato virus Y (PVY), but resistance-breaking strains (RB-PVY) have appeared. In an earlier study, we demonstrated that the loss-of-function of eIF(iso)4E-T reduces susceptibility to RB-PVY. Here, we show that simultaneous inhibition of eIF4E1-S and eIF(iso)4E-T synergistically confers enhanced resistance to both PVY and RB-PVY without host growth or development defects. PVY symptoms and accumulation in a tobacco line lacking eIF4E1-S were detected at 14 days post-inoculation (dpi) and RB-PVY symptoms in lines without functional eIF(iso)4E-T were observed at 24 dpi. RB-PVY emerged in a PVY-infected tobacco line lacking eIF4E1-S. In contrast, lines without functional eIF4E1-S and eIF(iso)4E-T were nearly immune to PVY and RB-PVY, and little accumulation of either virus was detected even at 56 dpi. Thus, the lines will be promising for PVY-resistance breeding. This study provides a novel strategy to develop tobacco highly resistant to PVY and RB-PVY, and insights into the mechanisms responsible for high-level resistance.

Genetic Variability and Molecular Evolution of Maize Yellow Mosaic Virus Populations from Different Geographic Origins

Maize yellow mosaic virus (MaYMV) hosted in various gramineous plants was assigned to the genus Polerovirus (family Luteoviridae) in 2018. However, little is known about its genetic diversity and population structure. In this study, 509 sugarcane leaf samples with mosaic symptoms were collected in 2017 to 2019 from eight sugarcane-growing provinces in China. Reverse-transcription PCR results revealed that four positive-sense RNA viruses were found to infect sugarcane, and the incidence of MaYMV among samples from Fujian, Sichuan, and Guangxi Provinces was 52.1, 9.8, and 2.5%, respectively. Based on 82 partial MaYMV sequences and 46 whole-genome sequences from different host plants, phylogenetic analysis revealed that MaYMV populations are very closely associated with their source geographical regions (China, Africa, and South America). Pairwise identity analysis showed significant variability in genome sequences among MaYMV isolates with genomic nucleotide identities of 91.1 to 99.9%. In addition to codon mutations, insertions or deletions also contributed to genetic variability in individual coding regions, especially in the readthrough protein (P3-P5 fusion protein). Low gene flow and significant genetic differentiation of MaYMV were observed among the three geographical populations, suggesting that environmental adaptation is an important evolutionary force that shapes the genetic structure of MaYMV. Genes in the MaYMV genome were subject to strong negative or purification selection during evolution, except for the movement protein (MP), which was under positive selection pressure. This finding suggests that the MP may play an important role in MaYMV evolution. Taken together, our findings provide basic information for the development of an integrated disease management strategy against MaYMV.

Regulatory barriers to improving global food security

Crop agriculture and food production constantly face climactic challenges to the supply of safe, nutritious food. These challenges highlight the importance of innovation resulting in improved crop technologies, capable of providing consistently increasing yields in the face of abiotic and biotic stresses. This article addresses the challenge that regulatory barriers are, and can, have on the adoption of innovative crop and food technologies that improve food security. Evidence of increased crop yield and the potential for increased yields, are presented from innovative plant breeding technologies, especially gene editing. Recent advances from the use of gene editing in the pharmaceutical field may offer opportunities to reduce regulatory burdens.

•

The EU regulatory framework for GM crops is in gridlock with no signs of ending.

•

The CJEU 2018 gene editing ruling has resulted in a decline in EU R&D investment.

•

Gene editing regulation in the EU will have spillover effects in Africa and Asia.

•

Covid-19 vaccines may require regulation as equivalent to GMOs, delaying approval.

Analysis of codon usage bias in potato virus Y non-recombinant strains

Potato virus Y (PVY) is a member of the genus Potyvirus, family Potyviridae, is considered one of the most devastating pest affecting economically important crops, such as potato, tobacco, tomato and pepper, representing a serious threat due to high incidence and worldwide distribution. Its economic significance as well as it biological and molecular complexities have aroused great attention, thus several studies have explore it genetic characteristics. However, little is known about PVY codon usage. To shed light on the relation of codon usage among viruses and their hosts is extremely important to understand virus survival, fitness and evolution. In this study, we performed a comprehensive analysis of codon usage and composition of PVY non-recombinant strains (PVY^N-NA, PVY^Eu-N, PVY^O, PVY^O5, PVY^C) based on 130 complete open reading frame sequences extracted from public databases. Furthermore, similarities between the synonymous codon usage of PVY and its main hosts were investigated. The results obtained in the current study suggest that the overall codon usage among PVY genotypes is similar and slightly biased. PVY codon usage is strongly influenced by mutational bias, but also by G + C compositional constraint and dinucleotide composition. Furthermore, similarities among codon usage preferences between PVY strains and analyzed hosts were observed.

Identification of Loci Associated with Enhanced Virulence in Spodoptera litura Nucleopolyhedrovirus Isolates Using Deep Sequencing

Spodoptera litura is an emerging pest insect in cotton and arable crops in Central Asia. To explore the possibility of using baculoviruses as biological control agents instead of chemical pesticides, in a previous study we characterized a number of S. litura nucleopolyhedrovirus (SpltNPV) isolates from Pakistan. We found significant differences in speed of kill, an important property of a biological control agent. Here we set out to understand the genetic basis of these differences in speed of kill, by comparing the genome of the fast-killing SpltNPV-Pak-TAX1 isolate with that of the slow-killing SpltNPV-Pak-BNG isolate. These two isolates and the SpltNPV-G2 reference strain from China were deep sequenced with Illumina. As expected, the two Pakistani isolates were closely related with >99% sequence identity, whereas the Chinese isolate was more distantly related. We identified two loci that may be associated with the fast action of the SpltNPV-Pak-TAX1 isolate. First, an analysis of rates of synonymous and non-synonymous mutations identified neutral to positive selection on open reading frame (ORF) 122, encoding a viral fibroblast growth factor (vFGF) that is known to affect virulence in other baculoviruses. Second, the homologous repeat region hr17, a putative enhancer of transcription and origin of replication, is absent in SpltNPV-Pak-TAX1 suggesting it may also affect virulence. Additionally, we found there is little genetic variation within both Pakistani isolates, and we identified four genes under positive selection in both isolates that may have played a role in adaptation of SpltNPV to conditions in Central Asia. Our results contribute to the understanding of the enhanced activity of SpltNPV-Pak-TAX1, and may help to select better SpltNPV isolates for the control of S. litura in Pakistan and elsewhere.

On the stability of sequences inserted into viral genomes

Viruses are widely used as vectors for heterologous gene expression in cultured cells or natural hosts, and therefore a large number of viruses with exogenous sequences inserted into their genomes have been engineered. Many of these engineered viruses are viable and express heterologous proteins at high levels, but the inserted sequences often prove to be unstable over time and are rapidly lost, limiting heterologous protein expression. Although virologists are aware that inserted sequences can be unstable, processes leading to insert instability are rarely considered from an evolutionary perspective. Here, we review experimental work on the stability of inserted sequences over a broad range of viruses, and we present some theoretical considerations concerning insert stability. Different virus genome organizations strongly impact insert stability, and factors such as the position of insertion can have a strong effect. In addition, we argue that insert stability not only depends on the characteristics of a particular genome, but that it will also depend on the host environment and the demography of a virus population. The interplay between all factors affecting stability is complex, which makes it challenging to develop a general model to predict the stability of genomic insertions. We highlight key questions and future directions, finding that insert stability is a surprisingly complex problem and that there is need for mechanism-based, predictive models. Combining theoretical models with experimental tests for stability under varying conditions can lead to improved engineering of viral modified genomes, which is a valuable tool for understanding genome evolution as well as for biotechnological applications, such as gene therapy.

Stochastic Modeling and Simulation of Viral Evolution

RNA viruses comprise vast populations of closely related, but highly genetically diverse, entities known as quasispecies. Understanding the mechanisms by which this extreme diversity is generated and maintained is fundamental when approaching viral persistence and pathobiology in infected hosts. In this paper, we access quasispecies theory through a mathematical model based on the theory of multitype branching processes, to better understand the roles of mechanisms resulting in viral diversity, persistence and extinction. We accomplish this understanding by a combination of computational simulations and the theoretical analysis of the model. In order to perform the simulations, we have implemented the mathematical model into a computational platform capable of running simulations and presenting the results in a graphical format in real time. Among other things, we show that the establishment of virus populations may display four distinct regimes from its introduction into new hosts until achieving equilibrium or undergoing extinction. Also, we were able to simulate different fitness distributions representing distinct environments within a host which could either be favorable or hostile to the viral success. We addressed the most used mechanisms for explaining the extinction of RNA virus populations called lethal mutagenesis and mutational meltdown. We were able to demonstrate a correspondence between these two mechanisms implying the existence of a unifying principle leading to the extinction of RNA viruses.

Impact of genetic drift, selection and accumulation level on virus adaptation to its host plants

The efficiency of plant major resistance genes is limited by the emergence and spread of resistance-breaking mutants. Modulation of the evolutionary forces acting on pathogen populations constitutes a promising way to increase the durability of these genes. We studied the effect of four plant traits affecting these evolutionary forces on the rate of resistance breakdown (RB) by a virus. Two of these traits correspond to virus effective population sizes (N_e ) at either plant inoculation or during infection. The third trait corresponds to differential selection exerted by the plant on the virus population. Finally, the fourth trait corresponds to within-plant virus accumulation (VA). These traits were measured experimentally on Potato virus Y (PVY) inoculated to a set of 84 pepper doubled-haploid lines, all carrying the same pvr2³ resistance gene, but having contrasting genetic backgrounds. The lines showed extensive variation for the rate of pvr2³ RB by PVY and for the four other traits of interest. A generalized linear model showed that three of these four traits, with the exception of N_e at inoculation, and several pairwise interactions between them had significant effects on RB. RB increased with increasing values of N_e during plant infection or VA. The effect of differential selection was more complex because of a strong interaction with VA. When VA was high, RB increased as the differential selection increased. An opposite relationship between RB and differential selection was observed when VA was low. This study provides a framework to select plants with appropriate virus evolution-related traits to avoid or delay RB.

Tobamoviruses as Models for the Study of Virus Evolution

The study of tobacco mosaic virus and other tobamovirus species has greatly contributed to the development of all areas of virology, including virus evolution. Research with tobamoviruses has been pioneer, or particularly significant, in all major areas of research in this field, including: the characterization of the genetic diversity of virus populations, the mechanisms and rates of generation of genetic diversity, the analysis of the genetic structure of virus populations and of the factors that shape it, the adaptation of viruses to hosts and the evolution of host range, and the evolution of virus taxa and of virus-host interactions. Many of these continue to be hot topics in evolutionary biology, or have been identified recently as such, including (i) host-range evolution, (ii) predicting the overcoming of resistance in crops, (iii) trade-offs between virus life-history traits in virus evolution, and (iv) the codivergence of viruses and hosts at different taxonomical and spatial scales. Tobamoviruses may be particularly appropriate to address these topics with plant viruses, as they provide convenient experimental systems, and as the detailed knowledge on their molecular and structural biology allows the analysis of the mechanisms behind evolutionary processes. Also, the extensive information on parameters related to infection dynamics and population structure may facilitate the development of realistic models to predict virus evolution. Certainly, tobamoviruses will continue to be favorite system for the study of virus evolution.

The diversity, evolution and epidemiology of plant viruses: A phylogenetic view

During the past four decades, the scientific community has seen an exponential advance in the number, sophistication, and quality of molecular techniques and bioinformatics tools for the genetic characterization of plant virus populations. Predating these advances, the field of Phylogenetics has significantly contributed to understand important aspects of plant virus evolution. This review aims at summarizing the impact of Phylogenetics in the current knowledge on three major aspects of plant virus evolution that have benefited from the development of phylogenetic inference: (1) The identification and classification of plant virus diversity. (2) The mechanisms and forces shaping the evolution of plant virus populations. (3) The understanding of the interaction between plant virus evolution, epidemiology and ecology. The work discussed here highlights the important role of phylogenetic approaches in the study of the dynamics of plant virus populations.

Going, going, gone: predicting the fate of genomic insertions in plant RNA viruses

Horizontal gene transfer is common among viruses, while they also have highly compact genomes and tend to lose artificial genomic insertions rapidly. Understanding the stability of genomic insertions in viral genomes is therefore relevant for explaining and predicting their evolutionary patterns. Here, we revisit a large body of experimental research on a plant RNA virus, tobacco etch potyvirus (TEV), to identify the patterns underlying the stability of a range of homologous and heterologous insertions in the viral genome. We obtained a wide range of estimates for the recombination rate-the rate at which deletions removing the insertion occur-and these appeared to be independent of the type of insertion and its location. Of the factors we considered, recombination rate was the best predictor of insertion stability, although we could not identify the specific sequence characteristics that would help predict insertion instability. We also considered experimentally the possibility that functional insertions lead to higher mutational robustness through increased redundancy. However, our observations suggest that both functional and non-functional increases in genome size decreased the mutational robustness. Our results therefore demonstrate the importance of recombination rates for predicting the long-term stability and evolution of viral RNA genomes and suggest that there are unexpected drawbacks to increases in genome size for mutational robustness.

Dynamics of PVY strains in field grown potato: Impact of strain competition and ability to overcome host resistance mechanisms

Potato virus Y (PVY) is the most important viral pathogen affecting potato crops worldwide. PVY can be transmitted non-persistently by aphids that do not colonize the host plant, resulting in a rapid acquisition and transmission of the virus between plants. PVY exists as a complex of strains that can be distinguished according to their pathogenicity, serology and genomic analysis. While virus incidence remains low in Scottish seed potato crops, PVY has become the increasingly prevalent virus. The monitoring of PVY^N and PVY^O serotypes has revealed a recent shift towards PVY^N which now accounts for more than 90% of all PVY cases. A survey of the molecular diversity of PVY^N isolates indicated that 80%-90% belong to the recombinant European (EU)-NTN group, with North-American (NA)-NTN and non-recombinant EU-N variants accounting for the remainder. The shift from non-recombinant to recombinant PVY isolates is a common trend observed worldwide. Surveys of a range of PVY isolates representing the main strain and phylogenetic groups suggest that PVY has the ability to overcome hypersensitive response-mediated resistance with significant differences between isolates of the same strain group. Contrastingly, genes mediating extreme resistance (Ry_adg, Ry_sto) provide efficient resistance to PVY transmission to progeny tubers. Transmission experiments in field conditions of PVY isolates representing the three main molecular groups (PVY^O, PVY^EU-NTN, PVY^NA-NTN) indicate that PVY^EU-NTN has the highest transmission rate. Our results suggest that PVY^EU-NTN isolate has a competitive advantage over PVY^O and PVY^NA-NTN isolates which is likely to be an important factor in shaping the evolution of viruses and the population dynamics of PVY.

A novel twelve class fluctuation test reveals higher than expected mutation rates for influenza A viruses

Influenza virus’ low replicative fidelity contributes to its capacity for rapid evolution. Clonal sequencing and fluctuation tests have suggested that the influenza virus mutation rate is 2.7 × 10^–6 - 3.0 × 10^–5 substitutions per nucleotide per strand copied (s/n/r). However, sequencing assays are biased toward mutations with minimal fitness impacts and fluctuation tests typically investigate only a subset of all possible single nucleotide mutations. We developed a fluctuation test based on reversion to fluorescence in a set of virally encoded mutant green fluorescent proteins, which allowed us to measure the rates of selectively neutral mutations representative of the twelve different mutation types. We measured an overall mutation rate of 1.8 × 10^–4 s/n/r for PR8 (H1N1) and 2.5 × 10^–4 s/n/r for Hong Kong 2014 (H3N2) and a transitional bias of 2.7–3.6. Our data suggest that each replicated genome will have an average of 2–3 mutations and highlight the importance of mutational load in influenza virus evolution.

DOI:http://dx.doi.org/10.7554/eLife.26437.001

Efficient Transmission and Propagation of Tomato Chlorosis Virus by Simple Single-Leaflet Grafting

Tomato chlorosis virus (ToCV), a member of the genus Crinivirus, has caused an epidemic disease in tomato worldwide. ToCV is phloem-limited and transmitted by whiteflies in a semi-persistent manner, but not by mechanical inoculation. Experimental propagation of ToCV has been performed primarily by using whitefly-mediated inoculation. To develop a simple and convenient method for transmission of ToCV, we investigated grafting single-leaflets from tomato plants infected with ToCV to recipient tomato seedlings. Forty-one of 46 tomato seedlings tested were grafted successfully with single-leaflets infected with ToCV. Among them, 36 seedlings (87.8%) were systemically infected with ToCV and developed typical symptoms. Our results demonstrated that single-leaflet grafting could provide a sufficient amount of inoculum for the transmission of ToCV to the grafted seedlings.

Mechanisms of viral mutation

The remarkable capacity of some viruses to adapt to new hosts and environments is highly dependent on their ability to generate de novo diversity in a short period of time. Rates of spontaneous mutation vary amply among viruses. RNA viruses mutate faster than DNA viruses, single-stranded viruses mutate faster than double-strand virus, and genome size appears to correlate negatively with mutation rate. Viral mutation rates are modulated at different levels, including polymerase fidelity, sequence context, template secondary structure, cellular microenvironment, replication mechanisms, proofreading, and access to post-replicative repair. Additionally, massive numbers of mutations can be introduced by some virus-encoded diversity-generating elements, as well as by host-encoded cytidine/adenine deaminases. Our current knowledge of viral mutation rates indicates that viral genetic diversity is determined by multiple virus- and host-dependent processes, and that viral mutation rates can evolve in response to specific selective pressures.

Effect of Host Species on Topography of the Fitness Landscape for a Plant RNA Virus

Adaptive fitness landscapes are a fundamental concept in evolutionary biology that relate the genotypes of individuals to their fitness. In the end, the evolutionary fate of evolving populations depends on the topography of the landscape, that is, the numbers of accessible mutational pathways and possible fitness peaks (i.e., adaptive solutions). For a long time, fitness landscapes were only theoretical constructions due to a lack of precise information on the mapping between genotypes and phenotypes. In recent years, however, efforts have been devoted to characterizing the properties of empirical fitness landscapes for individual proteins or for microbes adapting to artificial environments. In a previous study, we characterized the properties of the empirical fitness landscape defined by the first five mutations fixed during adaptation of tobacco etch potyvirus (TEV) to a new experimental host, Arabidopsis thaliana Here we evaluate the topography of this landscape in the ancestral host Nicotiana tabacum By comparing the topographies of the landscapes for the two hosts, we found that some features remained similar, such as the existence of fitness holes and the prevalence of epistasis, including cases of sign and reciprocal sign epistasis that created rugged, uncorrelated, and highly random topographies. However, we also observed significant differences in the fine-grained details between the two landscapes due to changes in the fitness and epistatic interactions of some genotypes. Our results support the idea that not only fitness tradeoffs between hosts but also topographical incongruences among fitness landscapes in alternative hosts may contribute to virus specialization. IMPORTANCE Despite its importance for understanding virus evolutionary dynamics, very little is known about the topography of virus adaptive fitness landscapes, and even less is known about the effects that different host species and environmental conditions may have on this topography. To bridge this gap, we evaluated the topography of a small fitness landscape formed by all genotypes that result from every possible combination of the first five mutations fixed during adaptation of TEV to the novel host A. thaliana To assess the effect that host species may have on this topography, we evaluated the fitness of every genotype in both the ancestral and novel hosts. We found that both landscapes share some macroscopic properties, such as the existence of holes and being highly rugged and uncorrelated, yet they differ in microscopic details due to changes in the magnitude and sign of fitness and epistatic effects.

Integrated analyses using RNA-Seq data reveal viral genomes, single nucleotide variations, the phylogenetic relationship, and recombination for Apple stem grooving virus

Background

Next-generation sequencing (NGS) provides many possibilities for plant virology research. In this study, we performed integrated analyses using plant transcriptome data for plant virus identification using Apple stem grooving virus (ASGV) as an exemplar virus. We used 15 publicly available transcriptome libraries from three different studies, two mRNA-Seq studies and a small RNA-Seq study.

Results

We de novo assembled nearly complete genomes of ASGV isolates Fuji and Cuiguan from apple and pear transcriptomes, respectively, and identified single nucleotide variations (SNVs) of ASGV within the transcriptomes. We demonstrated the application of NGS raw data to confirm viral infections in the plant transcriptomes. In addition, we compared the usability of two de novo assemblers, Trinity and Velvet, for virus identification and genome assembly. A phylogenetic tree revealed that ASGV and Citrus tatter leaf virus (CTLV) are the same virus, which was divided into two clades. Recombination analyses identified six recombination events from 21 viral genomes.

Conclusions

Taken together, our in silico analyses using NGS data provide a successful application of plant transcriptomes to reveal extensive information associated with viral genome assembly, SNVs, phylogenetic relationships, and genetic recombination.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2994-6) contains supplementary material, which is available to authorized users.

Genome Sequence Analysis of New Isolates of the Winona Strain of Plum pox virus and the First Definitive Evidence of Intrastrain Recombination Events

Plum pox virus (PPV) is genetically diverse with nine different strains identified. Mutations, indel events, and interstrain recombination events are known to contribute to the genetic diversity of PPV. This is the first report of intrastrain recombination events that contribute to PPV's genetic diversity. Fourteen isolates of the PPV strain Winona (W) were analyzed including nine new strain W isolates sequenced completely in this study. Isolates of other strains of PPV with more than one isolate with the complete genome sequence available in GenBank were included also in this study for comparison and analysis. Five intrastrain recombination events were detected among the PPV W isolates, one among PPV C strain isolates, and one among PPV M strain isolates. Four (29%) of the PPV W isolates analyzed are recombinants; one of which (P2-1) is a mosaic, with three recombination events identified. A new interstrain recombinant event was identified between a strain M isolate and a strain Rec isolate, a known recombinant. In silico recombination studies and pairwise distance analyses of PPV strain D isolates indicate that a threshold of genetic diversity exists for the detectability of recombination events, in the range of approximately 0.78×10(-2) to 1.33×10(-2) mean pairwise distance. RDP4 analyses indicate that in the case of PPV Rec isolates there may be a recombinant breakpoint distinct from the obvious transition point of strain sequences. Evidence was obtained that indicates that the frequency of PPV recombination is underestimated, which may be true for other RNA viruses where low genetic diversity exists.

The Multiplicity of Cellular Infection Changes Depending on the Route of Cell Infection in a Plant Virus

The multiplicity of cellular infection (MOI) is the number of virus genomes of a given virus species that infect individual cells. This parameter chiefly impacts the severity of within-host population bottlenecks as well as the intensity of genetic exchange, competition, and complementation among viral genotypes. Only a few formal estimations of the MOI currently are available, and most theoretical reports have considered this parameter as constant within the infected host. Nevertheless, the colonization of a multicellular host is a complex process during which the MOI may dramatically change in different organs and at different stages of the infection. We have used both qualitative and quantitative approaches to analyze the MOI during the colonization of turnip plants by Turnip mosaic virus. Remarkably, different MOIs were observed at two phases of the systemic infection of a leaf. The MOI was very low in primary infections from virus circulating within the vasculature, generally leading to primary foci founded by a single genome. Each lineage then moved from cell to cell at a very high MOI. Despite this elevated MOI during cell-to-cell progression, coinfection of cells by lineages originating in different primary foci is severely limited by the rapid onset of a mechanism inhibiting secondary infection. Thus, our results unveil an intriguing colonization pattern where individual viral genomes initiate distinct lineages within a leaf. Kin genomes then massively coinfect cells, but coinfection by two distinct lineages is strictly limited.

IMPORTANCE

The MOI is the size of the viral population colonizing cells and defines major phenomena in virus evolution, like the intensity of genetic exchange and the size of within-host population bottlenecks. However, few studies have quantified the MOI, and most consider this parameter as constant during infection. Our results reveal that the MOI can depend largely on the route of cell infection in a systemically infected leaf. The MOI is usually one genome per cell when cells are infected from virus particles moving long distances in the vasculature, whereas it is much higher during subsequent cell-to-cell movement in mesophyll. However, a fast-acting superinfection exclusion prevents cell coinfection by merging populations originating from different primary foci within a leaf. This complex colonization pattern results in a situation where within-cell interactions are occurring almost exclusively among kin and explains the common but uncharacterized phenomenon of genotype spatial segregation in infected plants.

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.

Viruses roll the dice: the stochastic behavior of viral genome molecules accelerates viral adaptation at the cell and tissue levels

Recent studies on evolutionarily distant viral groups have shown that the number of viral genomes that establish cell infection after cell-to-cell transmission is unexpectedly small (1-20 genomes). This aspect of viral infection appears to be important for the adaptation and survival of viruses. To clarify how the number of viral genomes that establish cell infection is determined, we developed a simulation model of cell infection for tomato mosaic virus (ToMV), a positive-strand RNA virus. The model showed that stochastic processes that govern the replication or degradation of individual genomes result in the infection by a small number of genomes, while a large number of infectious genomes are introduced in the cell. It also predicted two interesting characteristics regarding cell infection patterns: stochastic variation among cells in the number of viral genomes that establish infection and stochastic inequality in the accumulation of their progenies in each cell. Both characteristics were validated experimentally by inoculating tobacco cells with a library of nucleotide sequence-tagged ToMV and analyzing the viral genomes that accumulated in each cell using a high-throughput sequencer. An additional simulation model revealed that these two characteristics enhance selection during tissue infection. The cell infection model also predicted a mechanism that enhances selection at the cellular level: a small difference in the replication abilities of coinfected variants results in a large difference in individual accumulation via the multiple-round formation of the replication complex (i.e., the replication machinery). Importantly, this predicted effect was observed in vivo. The cell infection model was robust to changes in the parameter values, suggesting that other viruses could adopt similar adaptation mechanisms. Taken together, these data reveal a comprehensive picture of viral infection processes including replication, cell-to-cell transmission, and evolution, which are based on the stochastic behavior of the viral genome molecules in each cell.

Molecular biology of potyviruses

Potyvirus is the largest genus of plant viruses causing significant losses in a wide range of crops. Potyviruses are aphid transmitted in a nonpersistent manner and some of them are also seed transmitted. As important pathogens, potyviruses are much more studied than other plant viruses belonging to other genera and their study covers many aspects of plant virology, such as functional characterization of viral proteins, molecular interaction with hosts and vectors, structure, taxonomy, evolution, epidemiology, and diagnosis. Biotechnological applications of potyviruses are also being explored. During this last decade, substantial advances have been made in the understanding of the molecular biology of these viruses and the functions of their various proteins. After a general presentation on the family Potyviridae and the potyviral proteins, we present an update of the knowledge on potyvirus multiplication, movement, and transmission and on potyvirus/plant compatible interactions including pathogenicity and symptom determinants. We end the review providing information on biotechnological applications of potyviruses.

Comparative mutational analyses of influenza A viruses

The error-prone RNA-dependent RNA polymerase (RdRP) and external selective pressures are the driving forces for RNA viral diversity. When confounded by selective pressures, it is difficult to assess if influenza A viruses (IAV) that have a wide host range possess comparable or distinct spontaneous mutational frequency in their RdRPs. We used in-depth bioinformatics analyses to assess the spontaneous mutational frequencies of two RdRPs derived from human seasonal (A/Wuhan/359/95; Wuhan) and H5N1 (A/Vietnam/1203/04; VN1203) viruses using the mini-genome system with a common firefly luciferase reporter serving as the template. High-fidelity reverse transcriptase was applied to generate high-quality mutational spectra which allowed us to assess and compare the mutational frequencies and mutable motifs along a target sequence of the two RdRPs of two different subtypes. We observed correlated mutational spectra (τ correlation P < 0.0001), comparable mutational frequencies (H3N2:5.8 ± 0.9; H5N1:6.0 ± 0.5), and discovered a highly mutable motif "(A)AAG" for both Wuhan and VN1203 RdRPs. Results were then confirmed with two recombinant A/Puerto Rico/8/34 (PR8) viruses that possess RdRP derived from Wuhan or VN1203 (RG-PR8×Wuhan(PB2, PB1, PA, NP) and RG-PR8×VN1203(PB2, PB1, PA, NP)). Applying novel bioinformatics analysis on influenza mutational spectra, we provide a platform for a comprehensive analysis of the spontaneous mutation spectra for an RNA virus.

Direct versus indirect radiation action in irradiated vegetal embryos

Maize is one of model plants useful for genetic investigations and also very important for its agrotechnical utilizations. Here the genotoxic effects of low dose X-rays and accelerated electrons in maize caryopses was carried out with focus on the influence of water content at the moment of seed irradiation. X-ray photon beam as well as accelerated electrons were provided with 2.40 Gy min−1 dose rate. Pre-soaked and dry maize caryopses were irradiated with 0.5–3.0–6.0 Gy. Cytogenetic investigations were carried out based on microscope observations of chromosomes stained by Feulgen method. The mitotic index was found diminished in hydrated samples indicating the negative influence of indirect effects of water radicals. As known the water radiolysis release free radicals that attack biomolecules in addition to the directly absorbed radiation impact. Slight positive influence of 0.5 Gy radiation dose on cell division was evidenced. Chromosomal aberrations were identified like: vagrand chromosomes, C-metaphases, picnotic chromosomes, chromatide bridges. General tendency of aberrant mitoses enhancing was recorded in watered samples — with up to the twice increase for 6.0 Gy radiation dose. The results evidenced the hydration role in monitoring cytogenetic effects of low dose radiations in plant systems -with possible biotechnological applications.

Variability in the mutation rates of RNA viruses

ABSTRACT: 

It is well established that RNA viruses show extremely high mutation rates, but less attention has been paid to the fact that their mutation rates also vary strongly, from 10-6 to 10-4 substitutions per nucleotide per cell infection. The causes explaining this variability are still poorly understood, but candidate factors are the viral genome size and polarity, host-specific gene expression patterns, or the intracellular environment. Differences between animal and plant viruses, or between arthropod-borne and directly transmitted viruses have also been postulated. Finally, RNA viruses may be able to regulate the rate at which new mutations spread in the population by modifying features of the viral infection cycle, such as lysis time.

Shrinkage of genome size in a plant RNA virus upon transfer of an essential viral gene into the host genome

Nonretroviral integrated RNA viruses (NIRVs) are genes of nonretroviral RNA viruses found in the genomes of many eukaryotic organisms. NIRVs are thought to sometimes confer virus resistance, meaning that they could impact spread of the virus in the host population. However, a NIRV that is expressed may also impact the evolution of virus populations within host organisms. Here, we experimentally addressed the evolution of a virus in a host expressing a NIRV using Tobacco etch virus (TEV), a plant RNA virus, and transgenic tobacco plants expressing its replicase, NIb. We found that a virus missing the NIb gene, TEV-ΔNIb, which is incapable of autonomous replication in wild-type plants, had a higher fitness than the full-length TEV in the transgenic plants. Moreover, when the full-length TEV was evolved by serial passages in transgenic plants, we observed genomic deletions within NIb--and in some cases the adjacent cistrons--starting from the first passage. When we passaged TEV and TEV-ΔNIb in transgenic plants, we found mutations in proteolytic sites, but these only occurred in TEV-ΔNIb lineages, suggesting the adaptation of polyprotein processing to altered NIb expression. These results raise the possibility that NIRV expression can indeed induce the deletion of the corresponding genes in the viral genome, resulting in the formation of viruses that are replication defective in hosts that do not express the same NIRV. Moreover, virus genome evolution was contingent upon the deletion of the viral replicase, suggesting NIRV expression could also alter patterns of virus evolution.

Replication of Plant Viruses

Viruses replicate using both their own genetic information and host cell components and machinery. The different genome types have different replication pathways which contain controls on linking the process with translation and movement around the cell as well as not compromising the infected cell. This chapter discusses the replication mechanisms, faults in replication and replication of viruses co-infecting cells.

Viruses replicate using both their own genetic information and host cell components and machinery. The different genome types have different replication pathways which contain controls on linking the process with translation and movement around the cell as well as not compromising the infected cell. This chapter discusses the replication mechanisms, faults in replication and replication of viruses coinfecting cells.

Estimation of the in vivo recombination rate for a plant RNA virus

Phylogenomic evidence suggested that recombination is an important evolutionary force for potyviruses, one of the larger families of plant RNA viruses. However, mixed-genotype potyvirus infections are marked by low levels of cellular coinfection, precluding template switching and recombination events between virus genotypes during genomic RNA replication. To reconcile these conflicting observations, we evaluated the in vivo recombination rate (rg) of Tobacco etch virus (TEV; genus Potyvirus, family Potyviridae) by coinfecting plants with pairs of genotypes marked with engineered restriction sites as neutral markers. The recombination rate was then estimated using two different approaches: (i) a classical approach that assumed recombination between marked genotypes can occur in the whole virus population, rendering an estimate of rg = 7.762 × 10(-8) recombination events per nucleotide site per generation, and (ii) an alternative method that assumed recombination between marked genotypes can occur only in coinfected cells, rendering a much higher estimate of rg = 3.427 × 10(-5) recombination events per nucleotide site per generation. This last estimate is similar to the TEV mutation rate, suggesting that recombination should be at least as important as point mutation in creating variability. Finally, we compared our mutation and recombination rate estimates to those reported for animal RNA viruses. Our analysis suggested that high recombination rates may be an unavoidable consequence of selection for fast replication at the cost of low fidelity.

Stability and fitness impact of the visually discernible Rosea1 marker in the Tobacco etch virus genome

Antirrhinum majus Rosea1 (Ros1) is an MYB-related transcription factor that induces anthocyanin biosynthesis in plant tissues, and has been shown to be suitable for visual tracking of virus infection in plants. However, activation of anthocyanin biosynthesis has far reaching effects on plant physiology and could consequently have negative effects on viral replication. Therefore, viruses carrying the Ros1 marker might have a low fitness and consequently rapidly lose the marker. To compare the stability of the Ros1 marker, we generated Tobacco etch virus (TEV) based constructs containing either Ros1 or the enhanced green fluorescent protein (eGFP) between the NIb and CP cistrons (TEV-Ros1 and TEV-eGFP, respectively). We measured the within-host competitive fitness of both viruses by direct competitions with a common competitor during infection of Nicotiana tabacum. The fitness of TEV-Ros1 was significantly lower than that of TEV-eGFP, and both recombinant viruses had a significantly lower fitness than the wild-type virus. Nevertheless, after seven weeks of infection in N. tabacum, similar levels of marker gene instability where found for both viruses. Despite lower fitness of the marked virus, Ros1 is therefore a viable alternative marker for tracking viral infection in plants.

Genotypic but not phenotypic historical contingency revealed by viral experimental evolution

BACKGROUND

The importance of historical contingency in determining the potential of viral populations to evolve has been largely unappreciated. Identifying the constraints imposed by past adaptations is, however, of importance for understanding many questions in evolutionary biology, such as the evolution of host usage dynamics by multi-host viruses or the emergence of escape mutants that persist in the absence of antiviral treatments. To address this issue, we undertook an experimental approach in which sixty lineages of Tobacco etch potyvirus that differ in their past evolutionary history and degree of adaptation to Nicotiana tabacum were allowed to adapt to this host for 15 rounds of within host multiplication and transfer. We thereafter evaluated the degree of adaptation to the new host as well as to the original ones and characterized the consensus sequence of each lineage.

RESULTS

We found that past evolutionary history did not determine the phenotypic outcome of this common host evolution phase, and that the signal of local adaptation to past hosts had largely disappeared. By contrast, evolutionary history left footprints at the genotypic level, since the majority of host-specific mutations present at the beginning of this experiment were retained in the end-point populations and may have affected which new mutations were consequently fixed. This resulted in further divergence between the sequences despite a shared selective environment.

CONCLUSIONS

The present experiment reinforces the idea that the answer to the question "How important is historical contingency in evolution?" strongly depends on the level of integration of the traits studied. A strong historical contingency was found for TEV genotype, whereas a weak effect of on phenotypic evolution was revealed. In an applied context, our results imply that viruses are not easily trapped into suboptimal phenotypes and that (re)emergence is not evolutionarily constrained.

Ultradeep sequencing analysis of population dynamics of virus escape mutants in RNAi-mediated resistant plants

Plant artificial micro-RNAs (amiRs) have been engineered to target viral genomes and induce their degradation. However, the exceptional evolutionary plasticity of RNA viruses threatens the durability of the resistance conferred by these amiRs. It has recently been shown that viral populations not experiencing strong selective pressure from an antiviral amiR may already contain enough genetic variability in the target sequence to escape plant resistance in an almost deterministic manner. Furthermore, it has also been shown that viral populations exposed to subinhibitory concentrations of the antiviral amiR speed up this process. In this article, we have characterized the molecular evolutionary dynamics of an amiR target sequence in a viral genome under both conditions. The use of Illumina ultradeep sequencing has allowed us to identify virus sequence variants at frequencies as low as 2 × 10(-6) and to track their variation in time before and after the viral population was able of successfully infecting plants fully resistant to the ancestral virus. We found that every site in the amiR-target sequence of the viral genome presented variation and that the variant that eventually broke resistance was sampled among the many coexisting ones. In this system, viral evolution in fully susceptible plants results from an equilibrium between mutation and genetic drift, whereas evolution in partially resistant plants originates from more complex dynamics involving mutation, selection, and drift.

Effects of potyvirus effective population size in inoculated leaves on viral accumulation and the onset of symptoms

Effective population size (N(e)) is a key parameter for understanding evolutionary processes, but it is generally not considered in epidemiological studies or in studying infections of individual hosts. Whether N(e) has an effect on the onset of symptoms and viral accumulation in Tobacco etch virus (TEV) infection of Nicotiana tabacum plants is considered here. Using mixtures of TEV variants carrying fluorescent markers, the dose dependence of N(e) was confirmed, and the inoculation procedure was found to be the main source of variation in these experiments. Whereas the onset of symptoms was independent of N(e), there was less and more variable accumulation at 6 days postinoculation for small N(e) values (N(e) < 5). The observed variation in accumulation was not heritable, however, suggesting that this variation was not due to the fixation of deleterious mutations in the small founder populations. On the other hand, virus-induced fluorescence and accumulation in the inoculated leaf were strongly N(e) dependent. Systemic accumulation was independent of N(e), although removal of the inoculated leaf led to a small reduction in systemic accumulation for small N(e) values. For whole plants, N(e)-dependent effects on accumulation were no longer observed at 9 days postinoculation. Therefore, the effects of N(e) on accumulation are due mainly to limited expansion in the inoculated leaf and are transient. In this system, N(e)-dependent effects will be strongest at low doses and early in infection. We conclude that N(e) can have implications for epidemiology and infection at the individual host level, beyond determining the rate of mixed-genotype infection.

The three faces of riboviral spontaneous mutation: spectrum, mode of genome replication, and mutation rate

Riboviruses (RNA viruses without DNA replication intermediates) are the most abundant pathogens infecting animals and plants. Only a few riboviral infections can be controlled with antiviral drugs, mainly because of the rapid appearance of resistance mutations. Little reliable information is available concerning i) kinds and relative frequencies of mutations (the mutational spectrum), ii) mode of genome replication and mutation accumulation, and iii) rates of spontaneous mutation. To illuminate these issues, we developed a model in vivo system based on phage Qß infecting its natural host, Escherichia coli. The Qß RT gene encoding the Read-Through protein was used as a mutation reporter. To reduce uncertainties in mutation frequencies due to selection, the experimental Qß populations were established after a single cycle of infection and selection against RT⁻ mutants during phage growth was ameliorated by plasmid-based RT complementation in trans. The dynamics of Qß genome replication were confirmed to reflect the linear process of iterative copying (the stamping-machine mode). A total of 32 RT mutants were detected among 7,517 Qß isolates. Sequencing analysis of 45 RT mutations revealed a spectrum dominated by 39 transitions, plus 4 transversions and 2 indels. A clear template•primer mismatch bias was observed: A•C>C•A>U•G>G•U> transversion mismatches. The average mutation rate per base replication was ≈9.1×10⁻⁶ for base substitutions and ≈2.3×10⁻⁷ for indels. The estimated mutation rate per genome replication, μ_g, was ≈0.04 (or, per phage generation, ≈0.08), although secondary RT mutations arose during the growth of some RT mutants at a rate about 7-fold higher, signaling the possible impact of transitory bouts of hypermutation. These results are contrasted with those previously reported for other riboviruses to depict the current state of the art in riboviral mutagenesis.

Viral disease is a subject of major concern in public health. Diseases produced by riboviruses (RNA viruses sensu stricto) represent a special urgency, because these viruses display an exceptional capability to generate resistance mutations against antiviral drugs. Unfortunately, little is known about the rate and nature of spontaneous mutation in riboviruses. Thus, characterization of their mutation process may be helpful in the development of improved ways to counteract riboviral diseases. In this study, we investigated the mutation process in vivo of a model ribovirus, the bacteriophage Qß, focusing on three key aspects: i) the kinds and relative frequencies of mutations, ii) the mode of genome replication, and iii) the rate of spontaneous mutation. Our results, combined with other information about riboviral mutagenesis, depict a ribovirus mutation spectrum largely dominated by transitions, a predominantly linear mode of genome replication, and a mutation rate per genome replication on the order of 0.04 for bacteriophages and plant viruses but perhaps an order of magnitude higher for mammalian riboviruses.

Biomedical implications of viral mutation and evolution

Mutation rates vary hugely across viruses and strongly determine their evolution. In addition, viral mutation and evolution are biomedically relevant because they can determine pathogenesis, vaccine efficacy and antiviral resistance. We review experimental methods for estimating viral mutation rates and how these estimates vary across viral groups, paying special attention to the more general trends. Recent advances positing a direct association between viral mutation rates and virulence, or the use of high-fidelity variants as attenuated vaccines, are also discussed. Finally, we review the implications of viral mutation and evolution for the design of rational antiviral therapies and for efficient epidemiological surveillance.

Luria-delbruck estimation of turnip mosaic virus mutation rate in vivo

A potential drawback of recent antiviral therapies based on the transgenic expression of artificial microRNAs is the ease with which viruses may generate escape mutations. Using a variation of the classic Luria-Delbrück fluctuation assay, we estimated that the spontaneous mutation rate in the artificial microRNA (amiR) target of a plant virus was ca. 6 × 10(-5) per replication event.

Simultaneous mutations in multi-viral proteins are required for soybean mosaic virus to gain virulence on soybean genotypes carrying different R genes

BACKGROUND

Genetic resistance is the most effective and sustainable approach to the control of plant pathogens that are a major constraint to agriculture worldwide. In soybean, three dominant R genes, i.e., Rsv1, Rsv3 and Rsv4, have been identified and deployed against Soybean mosaic virus (SMV) with strain-specificities. Molecular identification of virulent determinants of SMV on these resistance genes will provide essential information for the proper utilization of these resistance genes to protect soybean against SMV, and advance knowledge of virus-host interactions in general.

METHODOLOGY/PRINCIPAL FINDINGS

To study the gain and loss of SMV virulence on all the three resistance loci, SMV strains G7 and two G2 isolates L and LRB were used as parental viruses. SMV chimeras and mutants were created by partial genome swapping and point mutagenesis and then assessed for virulence on soybean cultivars PI96983 (Rsv1), L-29 (Rsv3), V94-5152 (Rsv4) and Williams 82 (rsv). It was found that P3 played an essential role in virulence determination on all three resistance loci and CI was required for virulence on Rsv1- and Rsv3-genotype soybeans. In addition, essential mutations in HC-Pro were also required for the gain of virulence on Rsv1-genotype soybean. To our best knowledge, this is the first report that CI and P3 are involved in virulence on Rsv1- and Rsv3-mediated resistance, respectively.

CONCLUSIONS/SIGNIFICANCE

Multiple viral proteins, i.e., HC-Pro, P3 and CI, are involved in virulence on the three resistance loci and simultaneous mutations at essential positions of different viral proteins are required for an avirulent SMV strain to gain virulence on all three resistance loci. The likelihood of such mutations occurring naturally and concurrently on multiple viral proteins is low. Thus, incorporation of all three resistance genes in a soybean cultivar through gene pyramiding may provide durable resistance to SMV.

Advances in plant virus evolution: translating evolutionary insights into better disease management

Recent studies in plant virus evolution are revealing that genetic structure and behavior of virus and viroid populations can explain important pathogenic properties of these agents, such as host resistance breakdown, disease severity, and host shifting, among others. Genetic variation is essential for the survival of organisms. The exploration of how these subcellular parasites generate and maintain a certain frequency of mutations at the intra- and inter-host levels is revealing novel molecular virus-plant interactions. They emphasize the role of host environment in the dynamic genetic composition of virus populations. Functional genomics has identified host factors that are transcriptionally altered after virus infections. The analyses of these data by means of systems biology approaches are uncovering critical plant genes specifically targeted by viruses during host adaptation. Also, a next-generation resequencing approach of a whole virus genome is opening new avenues to study virus recombination and the relationships between intra-host virus composition and pathogenesis. Altogether, the analyzed data indicate that systematic disruption of some specific parameters of evolving virus populations could lead to more efficient ways of disease prevention, eradication, or tolerable virus-plant coexistence.

Tempo and mode of plant RNA virus escape from RNA interference-mediated resistance

A biotechnological application of artificial microRNAs (amiRs) is the generation of plants that are resistant to virus infection. This resistance has proven to be highly effective and sequence specific. However, before these transgenic plants can be deployed in the field, it is important to evaluate the likelihood of the emergence of resistance-breaking mutants. Two issues are of particular interest: (i) whether such mutants can arise in nontransgenic plants that may act as reservoirs and (ii) whether a suboptimal expression level of the transgene, resulting in subinhibitory concentrations of the amiR, would favor the emergence of escape mutants. To address the first issue, we experimentally evolved independent lineages of Turnip mosaic virus (TuMV) (family Potyviridae) in fully susceptible wild-type Arabidopsis thaliana plants and then simulated the spillover of the evolving virus to fully resistant A. thaliana transgenic plants. To address the second issue, the evolution phase took place with transgenic plants that expressed the amiR at subinhibitory concentrations. Our results show that TuMV populations replicating in susceptible hosts accumulated resistance-breaking alleles that resulted in the overcoming of the resistance of fully resistant plants. The rate at which resistance was broken was 7 times higher for TuMV populations that experienced subinhibitory concentrations of the antiviral amiR. A molecular characterization of escape alleles showed that they all contained at least one nucleotide substitution in the target sequence, generally a transition of the G-to-A and C-to-U types, with many instances of convergent molecular evolution. To better understand the viral population dynamics taking place within each host, as well as to evaluate relevant population genetic parameters, we performed in silico simulations of the experiments. Together, our results contribute to the rational management of amiR-based antiviral resistance in plants.

Dynamics of a plant RNA virus intracellular accumulation: stamping machine vs. geometric replication

The tremendous evolutionary potential of RNA viruses allows them to thrive despite host defense mechanisms and endows them with properties such as emergence, host switching, and virulence. The frequency of mutant viruses after an infectious process results from the interplay between the error rate of the viral replicase, from purifying mechanisms acting after transcription on aberrant RNAs, and from the amplification dynamics of virus RNA positive (+) and negative (-) strands. Two extreme scenarios describing viral RNA amplification are the geometric growth, in which each RNA strand serves as template for the synthesis of complementary strands with the same efficiency, and the stamping machine, where a strand is reiteratively used as template to synthesize multiple copies of the complementary. The resulting mutation frequencies are completely different, being geometric growth largely more mutagenic than stamping machine. In this work we evaluate the contribution of geometric growth and stamping machine to the overall genome amplification of the plant (+)-strand RNA virus turnip mosaic virus (TuMV). By means of transfection experiments of Nicotiana benthamiana protoplasts with a TuMV cDNA infectious clone and by using strand-specific quantitative real-time PCR, we determined the amplification dynamics of viral (+) and (-) RNA during a single-cell infectious process. A mathematical model describing the amplification of each viral strand was fitted to the data. Analyses of the model parameters showed that TuMV (+) and (-) RNA amplification occurs through a mixed strategy with ∼93% of genomes produced via stamping machine and only ∼7% resulting from geometric growth.

Tridimensional model structure and patterns of molecular evolution of Pepino mosaic virus TGBp3 protein

Background

Pepino mosaic virus (PepMV) is considered one of the most dangerous pathogens infecting tomatoes worldwide. The virus is highly diverse and four distinct genotypes, as well as inter-strain recombinants, have already been described. The isolates display a wide range on symptoms on infected plant species, ranging from mild mosaic to severe necrosis. However, little is known about the mechanisms and pattern of PepMV molecular evolution and about the role of individual proteins in host-pathogen interactions.

Methods

The nucleotide sequences of the triple gene block 3 (TGB3) from PepMV isolates varying in symptomatology and geographic origin have been analyzed. The modes and patterns of molecular evolution of the TGBp3 protein were investigated by evaluating the selective constraints to which particular amino acid residues have been subjected during the course of diversification. The tridimensional structure of TGBp3 protein has been modeled de novo using the Rosetta algorithm. The correlation between symptoms development and location of specific amino acids residues was analyzed.

Results

The results have shown that TGBp3 has been evolving mainly under the action of purifying selection operating on several amino acid sites, thus highlighting its functional role during PepMV infection. Interestingly, amino acid 67, which has been previously shown to be a necrosis determinant, was found to be under positive selection.

Conclusions

Identification of diverse selection events in TGB3p3 will help unraveling its biological functions and is essential to an understanding of the evolutionary constraints exerted on the Potexvirus genome. The estimated tridimensional structure of TGBp3 will serve as a platform for further sequence, structural and function analysis and will stimulate new experimental advances.

Experimental evolution of RNA versus DNA viruses

Based on their extremely high mutation rates, RNA viruses have been traditionally considered as the fastest evolving entities in nature. However, recent work has revealed that, despite their greater replication fidelity, single-stranded (ss) DNA viruses can evolve fast in a similar way. To further investigate this issue, we have compared the rates of adaptation and molecular evolution of ssRNA and ssDNA viruses under highly controlled laboratory conditions using the bacteriophages ΦX174, G4, f1, Qβ, SP, and MS2 as model systems. Our results indicate that ssRNA phages evolve faster than ssDNA phages under strong selective pressure, and that their extremely high mutation rates appear to be optimal for this kind of scenario. However, their performance becomes similar to that of ssDNA phages over the longer term or when the population is moderately well-adapted. Interestingly, the roughly 100-fold difference between the mutation rates of ssRNA and ssDNA phages yields less than a fivefold difference in adaptation and nucleotide substitution rates. The results are therefore consistent with the observation that, despite their lower mutation rates, ssDNA viruses can sometimes match the evolvability of RNA viruses.