Time-Sampled Population Sequencing Reveals the Interplay of Selection and Genetic Drift in Experimental Evolution of Potato Virus Y

Detection of single nucleotide polymorphisms in virus genomes assembled from high-throughput sequencing data: large-scale performance testing of sequence analysis strategies

Recent developments in high-throughput sequencing (HTS) technologies and bioinformatics have drastically changed research in virology, especially for virus discovery. Indeed, proper monitoring of the viral population requires information on the different isolates circulating in the studied area. For this purpose, HTS has greatly facilitated the sequencing of new genomes of detected viruses and their comparison. However, bioinformatics analyses allowing reconstruction of genome sequences and detection of single nucleotide polymorphisms (SNPs) can potentially create bias and has not been widely addressed so far. Therefore, more knowledge is required on the limitations of predicting SNPs based on HTS-generated sequence samples. To address this issue, we compared the ability of 14 plant virology laboratories, each employing a different bioinformatics pipeline, to detect 21 variants of pepino mosaic virus (PepMV) in three samples through large-scale performance testing (PT) using three artificially designed datasets. To evaluate the impact of bioinformatics analyses, they were divided into three key steps: reads pre-processing, virus-isolate identification, and variant calling. Each step was evaluated independently through an original, PT design including discussion and validation between participants at each step. Overall, this work underlines key parameters influencing SNPs detection and proposes recommendations for reliable variant calling for plant viruses. The identification of the closest reference, mapping parameters and manual validation of the detection were recognized as the most impactful analysis steps for the success of the SNPs detections. Strategies to improve the prediction of SNPs are also discussed.

Gathered from the Vine: A Survey of Seven Grapevine Viruses Within New England Vineyards

Vineyards in the Southeastern New England American Viticultural Area were surveyed for the incidence of seven major viruses: grapevine leafroll-associated viruses (GLRaV-1, GLRaV-2, GLRaV-3, and GLRaV-4), grapevine fanleaf virus (GFLV), tomato ringspot virus (ToRSV), and tobacco ringspot virus (TRSV). Viruses were detected by DAS-ELISA and confirmed by RT-PCR and Sanger sequencing. Multiple viruses were present in 19 out of the 25 vineyards surveyed between 2018 and 2020. GLRaV-3 (27.59%) was the most prevalent virus followed by GLRaV-4 (14.90%), GLRaV-1 (13.52%), GLRaV-2 (11.03%), ToRSV (6.34%), GFLV (5.24%), and TRSV (2.62%). Furthermore, phylogenetic analyses of the viral partial genome sequences acquired in this study revealed that the grapevine viruses present in this area are diverse, indicating that they may have been introduced from different sources. Our findings stress the need for improving the sanitary status of planting materials to avoid the introduction and dissemination of viruses to vineyards in this important wine-producing region of New England.

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.

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.

Determinants of Persistent Patterns of Pepino Mosaic Virus Mixed Infections

Pepino mosaic virus (PepMV) has become a pandemic virus in tomato crops, causing important economic losses worldwide. In Spain, isolates of the EU and CH2 strains co-circulate, with PepMV-EU predominantly found in mixed infections. Simultaneous in planta mixed infections result in an asymmetric antagonism against PepMV-CH2, but the outcome of over-infections has never been tested. PepMV-EU and PepMV-CH2 time-lagged inoculations were performed, and viral accumulation was measured 10 days after challenge inoculation. PepMV-EU had a protective effect over PepMV-CH2; in contrast, the accumulation of PepMV-EU increased in plants pre-inoculated with PepMV-CH2 as compared to single infections. We also studied the effect of the type of infection on viral transmission. Independently of the nature of the infection (single or mixed), we observed a strong positive correlation between virus accumulation in the source plant and transmission, excluding mixed infection effects different than modulating viral accumulation. Finally, in order to determine the genetic variability of PepMV strains in single and mixed infections, a 430 nucleotide region was RT-PCR amplified from samples from a serial passages experiment and deep-sequenced. No significant differences were found in the number of nucleotide substitutions between single and mixed infections for PepMV-EU; in contrast, significant differences were found for PepMV-CH2, which was more variable in single than in mixed infections. Comparing PepMV-EU with PepMV-CH2, a higher nucleotide diversity was found for PepMV-CH2. Collectively, our data strongly suggest that PepMV mixed infections can impact the virus epidemiology by modulating in planta virus strain accumulation and diversification.

Intra-Population Competition during Adaptation to Increased Temperature in an RNA Bacteriophage

Evolution of RNA bacteriophages of the family Leviviridae is governed by the high error rates of their RNA-dependent RNA polymerases. This fact, together with their large population sizes, leads to the generation of highly heterogeneous populations that adapt rapidly to most changes in the environment. Throughout adaptation, the different mutants that make up a viral population compete with each other in a non-trivial process in which their selective values change over time due to the generation of new mutations. In this work we have characterised the intra-population dynamics of a well-studied levivirus, Qβ, when it is propagated at a higher-than-optimal temperature. Our results show that adapting populations experienced rapid changes that involved the ascent of particular genotypes and the loss of some beneficial mutations of early generation. Artificially reconstructed populations, containing a fraction of the diversity present in actual populations, fixed mutations more rapidly, illustrating how population bottlenecks may guide the adaptive pathways. The conclusion is that, when the availability of beneficial mutations under a particular selective condition is elevated, the final outcome of adaptation depends more on the occasional occurrence of population bottlenecks and how mutations combine in genomes than on the selective value of particular mutations.

Increasing temperature alters the within-host competition of viral strains and influences virus genetic variability

Environmental conditions can affect viral accumulation, virulence and adaptation, which have implications in the disease outcomes and efficiency of control measures. Concurrently, mixed viral infections are relevant in plants, being their epidemiology shaped by within-host virus–virus interactions. However, the extent in which the combined effect of variations in abiotic components of the plant ecological niche and the prevalence of mixed infections affect the evolutionary dynamics of viral populations is not well understood. Here, we explore the interplay between ecological and evolutionary factors during viral infections and show that isolates of two strains of Pepino mosaic potexvirus coexisted in tomato plants in a temperature-dependent continuum between neutral and antagonistic interactions. After a long-term infection, the mutational analysis of the evolved viral genomes revealed strain-specific single-nucleotide polymorphisms that were modulated by the interaction between the type of infection and temperature. These results suggest that the temperature is an ecological driver of virus-virus interactions, with an effect on the genetic diversity of individual viruses that are co-infecting an individual host. This research provides insights into the effect that changes in host growth temperatures might have on the evolutionary dynamics of viral populations in mixed infections.

Next-Generation Sequencing and the CRISPR-Cas Nexus: A Molecular Plant Virology Perspective

In recent years, next-generation sequencing (NGS) and contemporary Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-CRISPR-associated (Cas) technologies have revolutionized the life sciences and the field of plant virology. Both these technologies offer an unparalleled platform for sequencing and deciphering viral metagenomes promptly. Over the past two decades, NGS technologies have improved enormously and have impacted plant virology. NGS has enabled the detection of plant viruses that were previously undetectable by conventional approaches, such as quarantine and archeological plant samples, and has helped to track the evolutionary footprints of viral pathogens. The CRISPR-Cas-based genome editing (GE) and detection techniques have enabled the development of effective approaches to virus resistance. Different versions of CRISPR-Cas have been employed to successfully confer resistance against diverse plant viruses by directly targeting the virus genome or indirectly editing certain host susceptibility factors. Applications of CRISPR-Cas systems include targeted insertion and/or deletion, site-directed mutagenesis, induction/expression/repression of the gene(s), epigenome re-modeling, and SNPs detection. The CRISPR-Cas toolbox has been equipped with precision GE tools to engineer the target genome with and without double-stranded (ds) breaks or donor templates. This technique has also enabled the generation of transgene-free genetically engineered plants, DNA repair, base substitution, prime editing, detection of small molecules, and biosensing in plant virology. This review discusses the utilities, advantages, applications, bottlenecks of NGS, and CRISPR-Cas in plant virology.

Population Genomics of Plant Viruses: The Ecology and Evolution of Virus Emergence

The genomics era has revolutionized studies of adaptive evolution by monitoring large numbers of loci throughout the genomes of many individuals. Ideally, the investigation of emergence in plant viruses requires examining the population dynamics of both virus and host, their interactions with each other, with other organisms and the abiotic environment. Genetic mechanisms that affect demographic processes are now being studied with high-throughput technologies, traditional genetics methods, and new computational tools for big-data. In this review, we discuss the utility of these approaches to monitor and detect changes in virus populations within cells and individuals, and over wider areas across species and communities of ecosystems. The advent of genomics in virology has fostered a multidisciplinary approach to tackling disease risk. The ability to make sense of the information now generated in this integrated setting is by far the most substantial obstacle to the ultimate goal of plant virology to minimize the threats to food security posed by disease. To achieve this goal, it is imperative to understand and forecast how populations respond to future changes in complex natural systems.

A Path toward SARS-CoV-2 Attenuation: Metabolic Pressure on CTP Synthesis Rules the Virus Evolution

In the context of the COVID-19 pandemic, we describe here the singular metabolic background that constrains enveloped RNA viruses to evolve toward likely attenuation in the long term, possibly after a step of increased pathogenicity. Cytidine triphosphate (CTP) is at the crossroad of the processes allowing SARS-CoV-2 to multiply, because CTP is in demand for four essential metabolic steps. It is a building block of the virus genome, it is required for synthesis of the cytosine-based liponucleotide precursors of the viral envelope, it is a critical building block of the host transfer RNAs synthesis and it is required for synthesis of dolichol-phosphate, a precursor of viral protein glycosylation. The CCA 3'-end of all the transfer RNAs required to translate the RNA genome and further transcripts into the proteins used to build active virus copies is not coded in the human genome. It must be synthesized de novo from CTP and ATP. Furthermore, intermediary metabolism is built on compulsory steps of synthesis and salvage of cytosine-based metabolites via uridine triphosphate that keep limiting CTP availability. As a consequence, accidental replication errors tend to replace cytosine by uracil in the genome, unless recombination events allow the sequence to return to its ancestral sequences. We document some of the consequences of this situation in the function of viral proteins. This unique metabolic setup allowed us to highlight and provide a raison d'être to viperin, an enzyme of innate antiviral immunity, which synthesizes 3'-deoxy-3',4'-didehydro-CTP as an extremely efficient antiviral nucleotide.

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.

Experimental evolution for niche breadth in bacteriophage T4 highlights the importance of structural genes

Ecologists have long studied the evolution of niche breadth, including how variability in environments can drive the evolution of specialism and generalism. This concept is of particular interest in viruses, where niche breadth evolution may explain viral disease emergence, or underlie the potential for therapeutic measures like phage therapy. Despite the significance and potential applications of virus-host interactions, the genetic determinants of niche breadth evolution remain underexplored in many bacteriophages. In this study, we present the results of an evolution experiment with a model bacteriophage system, Escherichia virus T4, in several host environments: exposure to Escherichia coli C, exposure to E. coli K-12, and exposure to both E. coli C and E. coli K-12. This experimental framework allowed us to investigate the phenotypic and molecular manifestations of niche breadth evolution. First, we show that selection on different hosts led to measurable changes in phage productivity in all experimental populations. Second, whole-genome sequencing of experimental populations revealed signatures of selection. Finally, clear and consistent patterns emerged across the host environments, especially the presence of new mutations in phage structural genes-genes encoding proteins that provide morphological and biophysical integrity to a virus. A comparison of mutations found across functional gene categories revealed that structural genes acquired significantly more mutations than other categories. Our findings suggest that structural genes are central determinants in bacteriophage niche breadth.

Theoretical conditions for the coexistence of viral strains with differences in phenotypic traits: a bifurcation analysis

We investigate the dynamics of a wild-type viral strain which generates mutant strains differing in phenotypic properties for infectivity, virulence and mutation rates. We study, by means of a mathematical model and bifurcation analysis, conditions under which the wild-type and mutant viruses, which compete for the same host cells, can coexist. The coexistence conditions are formulated in terms of the basic reproductive numbers of the strains, a maximum value of the mutation rate and the virulence of the pathogens. The analysis reveals that parameter space can be divided into five regions, each with distinct dynamics, that are organized around degenerate Bogdanov-Takens and zero-Hopf bifurcations, the latter of which gives rise to a curve of transcritical bifurcations of periodic orbits. These results provide new insights into the conditions by which viral populations may contain multiple coexisting strains in a stable manner.

Identification and Molecular Characterization of Recombinant Potato Virus Y (PVY) in Potato from South Korea, PVYNTN Strain

Potato is an important source of food in South Korea, and viruses represent a significant threat to sustainable and profitable potato production. However, information about viruses affecting the potato crop in South Korea is limited. In 2017, potato plants of five cultivars exhibiting foliar mosaic, crinkling, and mottle were collected in two seed potato production areas, in Gangwon-do and Jeollabuk-do Provinces, and subjected to virus testing and characterization. Potato virus Y (PVY) was found associated with mosaic symptoms, and samples were characterized using reverse transcription polymerase chain reaction (RT-PCR) and whole genome sequencing. All analyzed PVY-positive samples were found to represent the same recombinant PVY strain: PVY^NTN. Three PVY isolates were subjected to whole genome sequencing using overlapping RT-PCR fragments and Sanger methodology, and all three were confirmed to represent strain PVY^NTNa after a recombination analysis of the complete genomes. In phylogenetic analysis, the three South Korean isolates were placed most closely to several PVY^NTNa isolates reported from Japan and Vietnam, suggesting a common source of infection. This is the first report and complete molecular characterization of a PVY^NTN strain present in the country, and because this strain induces tuber necrotic ringspot disease in susceptible cultivars of potato, appropriate management tools need to be implemented to mitigate potential tuber quality losses.

High-Throughput Sequencing Facilitates Characterization of a "Forgotten" Plant Virus: The Case of a Henbane Mosaic Virus Infecting Tomato

High-throughput sequencing has dramatically broadened the possibilities for plant virus research and diagnostics, enabling discovery of new or obscure viruses, and virus strains and rapid sequencing of their genomes. In this research, we employed high-throughput sequencing to discover a new virus infecting tomato, Henbane mosaic virus (Potyvirus, Potyviridae), which was first discovered at the beginning of 20th century in the United Kingdom in cultivated henbane. A field tomato plant with severe necrotic symptoms of unknown etiology was sampled in Slovenia and high-throughput sequencing analysis using small RNA and ribosomal RNA depleted total RNA approaches revealed a mixed infection with Potato virus M (Carlavirus, Betaflexiviridae), Southern tomato virus (Amalgavirus, Amalgamaviridae) and henbane mosaic virus in the sample. The complete genomic sequence of henbane mosaic virus was assembled from the sequencing reads. By re-inoculation of the infected material on selected test plants, henbane mosaic virus was isolated and a host range analysis was performed, demonstrating the virus was pathogenic on several plant species. Due to limited metadata in public repositories, the taxonomic identification of the virus isolate was initially putative. Thus, in the next step, we used small RNA sequencing to determine genomic sequences of four historic isolates of the virus, obtained from different virus collections. Phylogenetic analyses performed using this new sequence information enabled us to taxonomically position Henbane mosaic virus as a member of the Potyvirus genus within the chili veinal mottle virus phylogenetic cluster and define the relationship of the new tomato isolate with the historic ones, indicating the existence of at least four putative strains of the virus. The first detection of henbane mosaic virus in tomato and demonstration of its pathogenicity on this host is important for plant protection and commercial tomato production. Since the virus was initially present in a mixed infection, and its whole genome was not sequenced, it has probably been overlooked in routine diagnostics. This study confirms the applicability of a combination of high-throughput sequencing and classic plant virus characterization methods for identification and phylogenetic classification of obscure viruses and historical viral isolates, for which no or limited genome sequence data is available.

Next Generation Sequencing for Detection and Discovery of Plant Viruses and Viroids: Comparison of Two Approaches

Next generation sequencing (NGS) technologies are becoming routinely employed in different fields of virus research. Different sequencing platforms and sample preparation approaches, in the laboratories worldwide, contributed to a revolution in detection and discovery of plant viruses and viroids. In this work, we are presenting the comparison of two RNA sequence inputs (small RNAs vs. ribosomal RNA depleted total RNA) for the detection of plant viruses by Illumina sequencing. This comparison includes several viruses, which differ in genome organization and viroids from both known families. The results demonstrate the ability for detection and identification of a wide array of known plant viruses/viroids in the tested samples by both approaches. In general, yield of viral sequences was dependent on viral genome organization and the amount of viral reads in the data. A putative novel Cytorhabdovirus, discovered in this study, was only detected by analysing the data generated from ribosomal RNA depleted total RNA and not from the small RNA dataset, due to the low number of short reads in the latter. On the other hand, for the viruses/viroids under study, the results showed higher yields of viral sequences in small RNA pool for viroids and viruses with no RNA replicative intermediates (single stranded DNA viruses).