Four principles to establish a universal virus taxonomy.

A universal taxonomy of viruses is essential for a comprehensive view of the virus world and for communicating the complicated evolutionary relationships among viruses. However, there are major differences in the conceptualisation and approaches to virus classification and nomenclature among virologists, clinicians, agronomists, and other interested parties. Here, we provide recommendations to guide the construction of a coherent and comprehensive virus taxonomy, based on expert scientific consensus. Firstly, assignments of viruses should be congruent with the best attainable reconstruction of their evolutionary histories, i.e., taxa should be monophyletic. This fundamental principle for classification of viruses is currently included in the International Committee on Taxonomy of Viruses (ICTV) code only for the rank of species. Secondly, phenotypic and ecological properties of viruses may inform, but not override, evolutionary relatedness in the placement of ranks. Thirdly, alternative classifications that consider phenotypic attributes, such as being vector-borne (e.g., "arboviruses"), infecting a certain type of host (e.g., "mycoviruses," "bacteriophages") or displaying specific pathogenicity (e.g., "human immunodeficiency viruses"), may serve important clinical and regulatory purposes but often create polyphyletic categories that do not reflect evolutionary relationships. Nevertheless, such classifications ought to be maintained if they serve the needs of specific communities or play a practical clinical or regulatory role. However, they should not be considered or called taxonomies. Finally, while an evolution-based framework enables viruses discovered by metagenomics to be incorporated into the ICTV taxonomy, there are essential requirements for quality control of the sequence data used for these assignments. Combined, these four principles will enable future development and expansion of virus taxonomy as the true evolutionary diversity of viruses becomes apparent.

Emergence of two distinct spatial folds in a pair of plant virus proteins encoded by nested genes.

Virus genomes may encode overlapping or nested open reading frames that increase their coding capacity. It is not known whether the constraints on spatial structures of the two encoded proteins limit the evolvability of nested genes. We examine the evolution of a pair of proteins, p22 and p19, encoded by nested genes in plant viruses from the genus Tombusvirus. The known structure of p19, a suppressor of RNA silencing, belongs to the RAGNYA fold from the alpha+beta class. The structure of p22, the cell-to-cell movement protein from the 30K family widespread in plant viruses, is predicted with the AlphaFold approach, suggesting a single jelly-roll fold core from the all-beta class, structurally similar to capsid proteins from plant and animal viruses. The nucleotide and codon preferences impose modest constraints on the types of secondary structures encoded in the alternative reading frames, nonetheless allowing for compact, well-ordered folds from different structural classes in two similarly-sized nested proteins. Tombusvirus p22 emerged through radiation of the widespread 30K family, which evolved by duplication of a virus capsid protein early in the evolution of plant viruses, whereas lineage-specific p19 may have emerged by a stepwise increase in the length of the overprinted gene and incremental acquisition of functionally active secondary structure elements by the protein product. This evolution of p19 toward the RAGNYA fold represents one of the first documented examples of protein structure convergence in naturally occurring proteins.

Identification of epigenetically regulated genes involved in plant-virus interaction and their role in virus-triggered induced resistance.

BACKGROUND: Plant responses to a wide range of stresses are known to be regulated by epigenetic mechanisms. Pathogen-related investigations, particularly against RNA viruses, are however scarce. It has been demonstrated that Arabidopsis thaliana plants defective in some members of the RNA-directed DNA methylation (RdDM) or histone modification pathways presented differential susceptibility to the turnip mosaic virus. In order to identify genes directly targeted by the RdDM-related RNA Polymerase V (POLV) complex and the histone demethylase protein JUMONJI14 (JMJ14) during infection, the transcriptomes of infected mutant and control plants were obtained and integrated with available chromatin occupancy data for various epigenetic proteins and marks. RESULTS: A comprehensive list of virus-responsive gene candidates to be regulated by the two proteins was obtained. Twelve genes were selected for further characterization, confirming their dynamic regulation during the course of infection. Several epigenetic marks on their promoter sequences were found using in silico data, raising confidence that the identified genes are actually regulated by epigenetic mechanisms. The altered expression of six of these genes in mutants of the methyltransferase gene CURLY LEAF and the histone deacetylase gene HISTONE DEACETYLASE 19 suggests that some virus-responsive genes may be regulated by multiple coordinated epigenetic complexes. A temporally separated multiple plant virus infection experiment in which plants were transiently infected with one virus and then infected by a second one was designed to investigate the possible roles of the identified POLV- and JMJ14-regulated genes in wild-type (WT) plants. Plants that had previously been stimulated with viruses were found to be more resistant to subsequent virus challenge than control plants. Several POLV- and JMJ14-regulated genes were found to be regulated in virus induced resistance in WT plants, with some of them poisoned to be expressed in early infection stages. CONCLUSIONS: A set of confident candidate genes directly regulated by the POLV and JMJ14 proteins during virus infection was identified, with indications that some of them may be regulated by multiple epigenetic modules. A subset of these genes may also play a role in the tolerance of WT plants to repeated, intermittent virus infections.

Plant Virus Adaptation to New Hosts: A Multi-scale Approach.

Viruses are studied at each level of biological complexity: from within-cells to ecosystems. The same basic evolutionary forces and principles operate at each level: mutation and recombination, selection, genetic drift, migration, and adaptive trade-offs. Great efforts have been put into understanding each level in great detail, hoping to predict the dynamics of viral population, prevent virus emergence, and manage their spread and virulence. Unfortunately, we are still far from this. To achieve these ambitious goals, we advocate for an integrative perspective of virus evolution. Focusing in plant viruses, we illustrate the pervasiveness of the above-mentioned principles. Beginning at the within-cell level, we describe replication modes, infection bottlenecks, and cellular contagion rates. Next, we move up to the colonization of distal tissues, discussing the fundamental role of random events. Then, we jump beyond the individual host and discuss the link between transmission mode and virulence. Finally, at the community level, we discuss properties of virus-plant infection networks. To close this review we propose the multilayer network theory, in which elements at different layers are connected and submit to their own dynamics that feed across layers, resulting in new emerging properties, as a way to integrate information from the different levels.

Plant virus evolution under strong drought conditions results in a transition from parasitism to mutualism.

Environmental conditions are an important factor driving pathogens' evolution. Here, we explore the effects of drought stress in plant virus evolution. We evolved turnip mosaic potyvirus in well-watered and drought conditions in Arabidopsis thaliana accessions that differ in their response to virus infection. Virus adaptation occurred in all accessions independently of watering status. Drought-evolved viruses conferred a significantly higher drought tolerance to infected plants. By contrast, nonsignificant increases in tolerance were observed in plants infected with viruses evolved under standard watering. The magnitude of this effect was dependent on the plant accessions. Differences in tolerance were correlated to alterations in the expression of host genes, some involved in regulation of the circadian clock, as well as in deep changes in the balance of phytohormones regulating defense and growth signaling pathways. Our results show that viruses can promote host survival in situations of abiotic stress, with the magnitude of such benefit being a selectable trait.

An Evolved 5' Untranslated Region of Alfalfa Mosaic Virus Allows the RNA Transport of Movement-Defective Variants.

Although the coat protein (CP) has a relevant role in the long-distance movement of alfalfa mosaic virus (AMV) and brome mosaic virus (BMV), its precise function is not fully understood. Previous results showed that a specific interaction between the C termini of the movement protein (MP) and the cognate CP is required for systemic transport. Thus, we have performed a compensatory evolution experiment using an AMV RNA3 derivative defective in long-distance transport that carries a BMV MP lacking the C-terminal 48 residues and unable to interact with the AMV CP. After several passages, five independent evolution lineages were able to move long distance. The analysis of the viral RNA of these lineages showed the presence of three different modifications located exclusively at the 5' untranslated region (5' UTR). The three evolved 5' UTR variants accumulated comparable levels of viral RNA and CP but reduced the accumulation of virus particles and the affinity between the 5' UTR and the AMV CP. In addition, the evolved 5' UTR increased cell-to-cell transport for both the AMV RNA3 carrying the BMV MP and that carrying the AMV MP. Finally, the evolved 5' UTRs allowed the systemic transport of an AMV RNA3 carrying a CP mutant defective in virus particles and increased the systemic transport of several AMV RNA3 derivatives carrying different viral MPs associated with the 30K superfamily. Altogether, our findings indicate that virus particles are not required for the systemic transport of AMV but also that BMV MP is competent for the short- and long-distance transport without the interaction with the CP. IMPORTANCE The results obtained in the present work could challenge the view of the role of the virus particle in the systemic transport of plant viruses. In this sense, we show that two different MPs are competent to systemically transport the AMV genome without the requirement of the virus particles, as reported for viruses lacking a CP (e.g., Umbravirus). The incapability of the viral MP to interact with the CP triggered virus variants that evolved to reduce the formation of virus particles, probably to increase the accessibility of the MP to the viral progeny. Our results point to the idea that virus particles would not be necessary for the viral systemic transport but would be necessary for vector virus transmission. This idea is reinforced by the observation that heterologous MPs also increased the systemic transport of the AMV constructs that have reduced encapsidation capabilities.

Virus classification based on in-depth sequence analyses and development of demarcation criteria using the Betaflexiviridae as a case study.

Currently, many viruses are classified based on their genome organization and nucleotide/amino acid sequence identities of their capsid and replication-associated proteins. Although biological traits such as vector specificities and host range are also considered, this later information is scarce for the majority of recently identified viruses, characterized only from genomic sequences. Accordingly, genomic sequences and derived information are being frequently used as the major, if not only, criteria for virus classification and this calls for a full review of the process. Herein, we critically addressed current issues concerning classification of viruses in the family Betaflexiviridae in the era of high-throughput sequencing and propose an updated set of demarcation criteria based on a process involving pairwise identity analyses and phylogenetics. The proposed framework has been designed to solve the majority of current conundrums in taxonomy and to facilitate future virus classification. Finally, the analyses performed herein, alongside the proposed approaches, could be used as a blueprint for virus classification at-large.

Transmission modes affect the population structure of potato virus Y in potato.

Transmission is a crucial part of a viral life cycle and transmission mode can have an important impact on virus biology. It was demonstrated that transmission mode can influence the virulence and evolution of a virus; however, few empirical data are available to describe the direct underlying changes in virus population structure dynamics within the host. Potato virus Y (PVY) is an RNA virus and one of the most damaging pathogens of potato. It comprises several genetically variable strains that are transmitted between plants via different transmission modes. To investigate how transmission modes affect the within-plant viral population structure, we have used a deep sequencing approach to examine the changes in the genetic structure of populations (in leaves and tubers) of three PVY strains after successive passages by horizontal (aphid and mechanical) and vertical (via tubers) transmission modes. Nucleotide diversities of viral populations were significantly influenced by transmission modes; lineages transmitted by aphids were the least diverse, whereas lineages transmitted by tubers were the most diverse. Differences in nucleotide diversities of viral populations between leaves and tubers were transmission mode-dependent, with higher diversities in tubers than in leaves for aphid and mechanically transmitted lineages. Furthermore, aphid and tuber transmissions were shown to impose stronger genetic bottlenecks than mechanical transmission. To better understand the structure of virus populations within the host, transmission mode, movement of the virus within the host, and the number of replication cycles after transmission event need to be considered. Collectively, our results suggest a significant impact of virus transmission modes on the within-plant diversity of virus populations and provide quantitative fundamental data for understanding how transmission can shape virus diversity in the natural ecosystems, where different transmission modes are expected to affect virus population structure and consequently its evolution.

Why are viral genomes so fragile? The bottleneck hypothesis.

If they undergo new mutations at each replication cycle, why are RNA viral genomes so fragile, with most mutations being either strongly deleterious or lethal? Here we provide theoretical and numerical evidence for the hypothesis that genetic fragility is partly an evolutionary response to the multiple population bottlenecks experienced by viral populations at various stages of their life cycles. Modelling within-host viral populations as multi-type branching processes, we show that mutational fragility lowers the rate at which Muller's ratchet clicks and increases the survival probability through multiple bottlenecks. In the context of a susceptible-exposed-infectious-recovered epidemiological model, we find that the attack rate of fragile viral strains can exceed that of more robust strains, particularly at low infectivities and high mutation rates. Our findings highlight the importance of demographic events such as transmission bottlenecks in shaping the genetic architecture of viral pathogens.

Genetic Diversity of Tomato Black Ring Virus Satellite RNAs and Their Impact on Virus Replication.

Viral satellite RNAs (satRNAs) are small subviral particles that are associated with the genomic RNA of a helper virus (HV). Their replication, encapsidation, and movement depend on the HV. In this paper, we performed a global analysis of the satRNAs associated with different isolates of tomato black ring virus (TBRV). We checked the presence of satRNAs in 42 samples infected with TBRV, performed recombination and genetic diversity analyses, and examined the selective pressure affecting the satRNAs population. We identified 18 satRNAs in total that differed in length and the presence of point mutations. Moreover, we observed a strong effect of selection operating upon the satRNA population. We also constructed infectious cDNA clones of satRNA and examined the viral load of different TBRV isolates in the presence and absence of satRNAs, as well as the accumulation of satRNA molecules on infected plants. Our data provide evidence that the presence of satRNAs significantly affects viral load; however, the magnitude of this effect differs among viral isolates and plant hosts. We also showed a positive correlation between the number of viral genomic RNAs (gRNAs) and satRNAs for two analysed TBRV isolates.

Phase transitions in virology.

Viruses have established relationships with almost every other living organism on Earth and at all levels of biological organization: from other viruses up to entire ecosystems. In most cases, they peacefully coexist with their hosts, but in most relevant cases, they parasitize them and induce diseases and pandemics, such as the AIDS and the most recent avian influenza and COVID-19 pandemic events, causing a huge impact on health, society, and economy. Viruses play an essential role in shaping the eco-evolutionary dynamics of their hosts, and have been also involved in some of the major evolutionary innovations either by working as vectors of genetic information or by being themselves coopted by the host into their genomes. Viruses can be studied at different levels of biological organization, from the molecular mechanisms of genome replication, gene expression and encapsidation, to global pandemics. All these levels are different and yet connected through the presence of threshold conditions allowing for the formation of a capsid, the loss of genetic information or epidemic spreading. These thresholds, as occurs with temperature separating phases in a liquid, define sharp qualitative types of behaviour. Thesephase transitionsare very well known in physics. They have been studied by means of simple, but powerful models able to capture their essential properties, allowing us to better understand them. Can the physics of phase transitions be an inspiration for our understanding of viral dynamics at different scales? Here we review well-known mathematical models of transition phenomena in virology. We suggest that the advantages of abstract, simplified pictures used in physics are also the key to properly understanding the origins and evolution of complexity in viruses. By means of several examples, we explore this multilevel landscape and how minimal models provide deep insights into a diverse array of problems. The relevance of these transitions in connecting dynamical patterns across scales and their evolutionary and clinical implications are outlined.

Viral Fitness Determines the Magnitude of Transcriptomic and Epigenomic Reprograming of Defense Responses in Plants.

Although epigenetic factors may influence the expression of defense genes in plants, their role in antiviral responses and the impact of viral adaptation and evolution in shaping these interactions are still poorly explored. We used two isolates of turnip mosaic potyvirus with varying degrees of adaptation to Arabidopsis thaliana to address these issues. One of the isolates was experimentally evolved in the plant and presented increased load and virulence relative to the ancestral isolate. The magnitude of the transcriptomic responses was larger for the evolved isolate and indicated a role of innate immunity systems triggered by molecular patterns and effectors in the infection process. Several transposable elements located in different chromatin contexts and epigenetic-related genes were also affected. Correspondingly, mutant plants having loss or gain of repressive marks were, respectively, more tolerant and susceptible to turnip mosaic potyvirus, with a more efficient response against the ancestral isolate. In wild-type plants, both isolates induced similar levels of cytosine methylation changes, including in and around transposable elements and stress-related genes. Results collectively suggested that apart from RNA silencing and basal immunity systems, DNA methylation and histone modification pathways may also be required for mounting proper antiviral defenses and that the effectiveness of this type of regulation strongly depends on the degree of viral adaptation to the host.

Revisiting Orthotospovirus phylogeny using full-genome data and testing the contribution of selection, recombination and segment reassortment in the origin of members of new species.

The family Tospoviridae of the order Bunyavirales is constituted of tri-segmented negative-sense single-stranded RNA viruses that infect plants and are also able to replicate in their insect vectors in a persistent manner. The family is composed of a single genus, Orthotospovirus, whose type species is Tomato spotted wilt orthotospovirus. Previous studies assessing the phylogenetic relationships within this genus were based on partial genomic sequences, resulting in unresolved clades and a poor assessment of the roles of recombination and segment reassortment during mixed infections. Full genome sequences of members of recognized Orthotospovirus species are now available at NCBI. In this study, we examined 67 complete genome sequences from members of 22 species. Our study confirms the existence of four phylogroups (A to D), grouped in two major clades (A-B and C-D) within the genus. We found strong evidence that within-segment recombination events and reassortment of segments during mixed infections have been involved in the origin of new orthotospoviruses. Also, selection pressures were analyzed for each gene, and evidence of positive selection was found in all genes.

Host-virus evolutionary dynamics with specialist and generalist infection strategies: Bifurcations, bistability, and chaos.

In this work, we have investigated the evolutionary dynamics of a generalist pathogen, e.g., a virus population, that evolves toward specialization in an environment with multiple host types. We have particularly explored under which conditions generalist viral strains may rise in frequency and coexist with specialist strains or even dominate the population. By means of a nonlinear mathematical model and bifurcation analysis, we have determined the theoretical conditions for stability of nine identified equilibria and provided biological interpretation in terms of the infection rates for the viral specialist and generalist strains. By means of a stability diagram, we identified stable fixed points and stable periodic orbits, as well as regions of bistability. For arbitrary biologically feasible initial population sizes, the probability of evolving toward stable solutions is obtained for each point of the analyzed parameter space. This probability map shows combinations of infection rates of the generalist and specialist strains that might lead to equal chances for each type becoming the dominant strategy. Furthermore, we have identified infection rates for which the model predicts the onset of chaotic dynamics. Several degenerate Bogdanov-Takens and zero-Hopf bifurcations are detected along with generalized Hopf and zero-Hopf bifurcations. This manuscript provides additional insights into the dynamical complexity of host-pathogen evolution toward different infection strategies.

Viral Fitness Correlates with the Magnitude and Direction of the Perturbation Induced in the Host's Transcriptome: The Tobacco Etch Potyvirus-Tobacco Case Study.

Determining the fitness of viral genotypes has become a standard practice in virology as it is essential to evaluate their evolutionary potential. Darwinian fitness, defined as the advantage of a given genotype with respect to a reference one, is a complex property that captures, in a single figure, differences in performance at every stage of viral infection. To what extent does viral fitness result from specific molecular interactions with host factors and regulatory networks during infection? Can we identify host genes in functional classes whose expression depends on viral fitness? Here, we compared the transcriptomes of tobacco plants infected with seven genotypes of tobacco etch potyvirus that differ in fitness. We found that the larger the fitness differences among genotypes, the more dissimilar the transcriptomic profiles are. Consistently, two different mutations, one in the viral RNA polymerase and another in the viral suppressor of RNA silencing, resulted in significantly similar gene expression profiles. Moreover, we identified host genes whose expression showed a significant correlation, positive or negative, with the virus' fitness. Differentially expressed genes which were positively correlated with viral fitness activate hormone- and RNA silencing-mediated pathways of plant defense. In contrast, those that were negatively correlated with fitness affect metabolism, reducing growth, and development. Overall, these results reveal the high information content of viral fitness and suggest its potential use to predict differences in genomic profiles of infected hosts.

Evolving by deleting: patterns of molecular evolution of Apple stem pitting virus isolates from Poland.

In this study, 267 coat protein gene (CP) sequences from 48 Polish isolates of Apple stem pitting virus (ASPV) were determined. The genetic structure of the virus population was analysed and possible mechanisms of molecular evolution explored. We found evidence of recombination within the ASPV population and the presence of 17 ASPV molecular variants that differ in the length, number and arrangement of deletions in the CP. Population genetic analyses showed significant variation among isolates from pear and apple trees, between isolates from the same host species and, more interestingly, within isolates, supporting the existence of significant levels of variability within individual hosts, as expected by a quasispecies population structure. In addition, different tests support that selection might have been an important force driving diversification within isolates: positive selection was found acting upon certain amino acids. Phylogenetic analyses also showed that isolates did not classify according to the host species (pear or apple trees) but according to the pattern of deletions, suggesting a possible role for deletions during clade diversification.

Viroids: survivors from the RNA world?

Because RNA can be a carrier of genetic information and a biocatalyst, there is a consensus that it emerged before DNA and proteins, which eventually assumed these roles and relegated RNA to intermediate functions. If such a scenario--the so-called RNA world--existed, we might hope to find its relics in our present world. The properties of viroids that make them candidates for being survivors of the RNA world include those expected for primitive RNA replicons: (a) small size imposed by error-prone replication, (b) high G + C content to increase replication fidelity, (c) circular structure for assuring complete replication without genomic tags, (d) structural periodicity for modular assembly into enlarged genomes, (e) lack of protein-coding ability consistent with a ribosome-free habitat, and (f) replication mediated in some by ribozymes, the fingerprint of the RNA world. With the advent of DNA and proteins, those protoviroids lost some abilities and became the plant parasites we now know.

Viral replication modes in single-peak fitness landscapes: A dynamical systems analysis.

Positive-sense, single-stranded RNA viruses are important pathogens infecting almost all types of organisms. Experimental evidence from distributions of mutations and from viral RNA amplification suggest that these pathogens may follow different RNA replication modes, ranging from the stamping machine replication (SMR) to the geometric replication (GR) mode. Although previous theoretical work has focused on the evolutionary dynamics of RNA viruses amplifying their genomes with different strategies, little is known in terms of the bifurcations and transitions involving the so-called error threshold (mutation-induced dominance of mutants) and lethal mutagenesis (extinction of all sequences due to mutation accumulation and demographic stochasticity). Here we analyze a dynamical system describing the intracellular amplification of viral RNA genomes evolving on a single-peak fitness landscape focusing on three cases considering neutral, deleterious, and lethal mutants. We analytically derive the critical mutation rates causing lethal mutagenesis and error threshold, governed by transcritical bifurcations that depend on parameters α (parameter introducing the mode of replication), replicative fitness of mutants (k1), and on the spontaneous degradation rates of the sequences (ϵ). Our results relate the error catastrophe with lethal mutagenesis in a model with continuous populations of viral genomes. The former case involves dominance of the mutant sequences, while the latter, a deterministic extinction of the viral RNAs during replication due to increased mutation. For the lethal case the critical mutation rate involving lethal mutagenesis is µc=1-ɛ/α. Here, the SMR involves lower critical mutation rates, being the system more robust to lethal mutagenesis replicating closer to the GR mode. This result is also found for the neutral and deleterious cases, but for these later cases lethal mutagenesis can shift to the error threshold once the replication mode surpasses a threshold given by α=ϵ/k1.

Complete nucleotide sequence of a novel mycovirus from Trichoderma harzianum in China.

A new mycovirus was identified in Trichoderma harzianum strain 137 isolated in Xinjiang province, China. The whole genome sequence of the mycovirus was determined by metagenomic sequencing, RT-PCR, and RACE cloning. The mycovirus contained two genomic segments. The first segment was 2088 bp long and contained a single ORF (ORF1) encoding the RNA-dependent RNA polymerase (RdRP) (72.26 kDa). The second segment was 1634 bp long and also contained a single ORF (ORF2) encoding a hypothetical protein of 37.472 kDa. We named this novel mycovirus "Trichoderma harzianum bipartite mycovirus 1" (ThBMV1). Phylogenetic analysis showed that ThBMV1 clusters with other unclassified dsRNA mycoviruses.