Deep sequencing of virus-derived small interfering RNAs and RNA from viral particles shows highly similar mutational landscapes of a plant virus population

Diagnosis and Characterization of Plant Viruses Using HTS to Support Virus Management and Tomato Breeding

Viral diseases pose a significant threat to tomato crops (Solanum lycopersicum L.), one of the world's most economically important vegetable crops. The limited genetic diversity of cultivated tomatoes contributes to their high susceptibility to viral infections. To address this challenge, tomato breeding programs must harness the genetic resources found in native populations and wild relatives. Breeding efforts may aim to develop broad-spectrum resistance against the virome. To identify the viruses naturally infecting 19 advanced lines, derived from native tomatoes, high-throughput sequencing (HTS) of small RNAs and confirmation with PCR and RT-PCR were used. Single and mixed infections with tomato mosaic virus (ToMV), tomato golden mosaic virus (ToGMoV), and pepper huasteco yellow vein virus (PHYVV) were detected. The complete consensus genomes of three variants of Mexican ToMV isolates were reconstructed, potentially forming a new ToMV clade with a distinct 3' UTR. The absence of reported mutations associated with resistance-breaking to ToMV suggests that the Tm-1, Tm-2, and Tm-22 genes could theoretically be used to confer resistance. However, the high mutation rates and a 63 nucleotide insertion in the 3' UTR, as well as amino acid mutations in the ORFs encoding 126 KDa, 183 KDa, and MP of Mexican ToMV isolates, suggest that it is necessary to evaluate the capacity of these variants to overcome Tm-1, Tm-2, and Tm-22 resistance genes. This evaluation, along with the characterization of advanced lines using molecular markers linked to these resistant genes, will be addressed in future studies as part of the breeding strategy. This study emphasizes the importance of using HTS for accurate identification and characterization of plant viruses that naturally infect tomato germplasm based on the consensus genome sequences. This study provides crucial insights to select appropriate disease management strategies and resistance genes and guide breeding efforts toward the development of virus-resistant tomato varieties.

Application of High-Throughput Sequencing for Comprehensive Virome Profiling in Grapevines Shows Yellows in Iran

A comprehensive study on the whole spectrum of viruses and viroids in five Iranian grapevine cultivars was carried out using sRNA libraries prepared from phloem tissue. A comparison of two approaches to virus detection from sRNAome data indicated a significant difference in the results and performance of the aligners in viral genome reconstruction. The results showed a complex virome in terms of viral composition, abundance, and richness. Thirteen viruses and viroids were identified in five Iranian grapevine cultivars, among which the grapevine red blotch virus and grapevine satellite virus were detected for the first time in Iranian vineyards. Grapevine leafroll-associated virus 1 (GLRaV1) and grapevine fanleaf virus (GFLV) were highly dominant in the virome. However, their frequency and abundance were somewhat different among grapevine cultivars. The results revealed a mixed infection of GLRaV1/grapevine yellow speckle viroid 1 (GYSVd1) and GFLV/GYSVd1 in grapevines that exhibited yellows and vein banding. We also propose a threshold of 14% of complete reconstruction as an appropriate threshold for detection of grapevine viruses that can be used as indicators for reliable grapevine virome profiling or in quarantine stations and certification programs.

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.

Generation and Functional Analysis of Defective Viral Genomes during SARS-CoV-2 Infection

Defective viral genomes (DVGs) have been identified in many RNA viruses as a major factor influencing antiviral immune response and viral pathogenesis. However, the generation and function of DVGs in SARS-CoV-2 infection are less known. In this study, we elucidated DVG generation in SARS-CoV-2 and its relationship with host antiviral immune response. We observed DVGs ubiquitously from transcriptome sequencing (RNA-seq) data sets of in vitro infections and autopsy lung tissues of COVID-19 patients. Four genomic hot spots were identified for DVG recombination, and RNA secondary structures were suggested to mediate DVG formation. Functionally, bulk and single-cell RNA-seq analysis indicated the interferon (IFN) stimulation of SARS-CoV-2 DVGs. We further applied our criteria to the next-generation sequencing (NGS) data set from a published cohort study and observed a significantly higher amount and frequency of DVG in symptomatic patients than those in asymptomatic patients. Finally, we observed exceptionally diverse DVG populations in one immunosuppressive patient up to 140 days after the first positive test of COVID-19, suggesting for the first time an association between DVGs and persistent viral infections in SARS-CoV-2. Together, our findings strongly suggest a critical role of DVGs in modulating host IFN responses and symptom development, calling for further inquiry into the mechanisms of DVG generation and into how DVGs modulate host responses and infection outcome during SARS-CoV-2 infection. IMPORTANCE Defective viral genomes (DVGs) are generated ubiquitously in many RNA viruses, including SARS-CoV-2. Their interference activity to full-length viruses and IFN stimulation provide the potential for them to be used in novel antiviral therapies and vaccine development. SARS-CoV-2 DVGs are generated through the recombination of two discontinuous genomic fragments by viral polymerase complex, and this recombination is also one of the major mechanisms for the emergence of new coronaviruses. Focusing on the generation and function of SARS-CoV-2 DVGs, these studies identify new hot spots for nonhomologous recombination and strongly suggest that the secondary structures within viral genomes mediate the recombination. Furthermore, these studies provide the first evidence for IFN stimulation activity of de novo DVGs during natural SARS-CoV-2 infection. These findings set up the foundation for further mechanism studies of SARS-CoV-2 recombination and provide evidence to harness the immunostimulatory potential of DVGs in the development of a vaccine and antivirals for SARS-CoV-2.

ViReMa: a virus recombination mapper of next-generation sequencing data characterizes diverse recombinant viral nucleic acids

BACKGROUND

Genetic recombination is a tremendous source of intrahost diversity in viruses and is critical for their ability to rapidly adapt to new environments or fitness challenges. While viruses are routinely characterized using high-throughput sequencing techniques, characterizing the genetic products of recombination in next-generation sequencing data remains a challenge. Viral recombination events can be highly diverse and variable in nature, including simple duplications and deletions, or more complex events such as copy/snap-back recombination, intervirus or intersegment recombination, and insertions of host nucleic acids. Due to the variable mechanisms driving virus recombination and the different selection pressures acting on the progeny, recombination junctions rarely adhere to simple canonical sites or sequences. Furthermore, numerous different events may be present simultaneously in a viral population, yielding a complex mutational landscape.

FINDINGS

We have previously developed an algorithm called ViReMa (Virus Recombination Mapper) that bootstraps the bowtie short-read aligner to capture and annotate a wide range of recombinant species found within virus populations. Here, we have updated ViReMa to provide an "error density" function designed to accurately detect recombination events in the longer reads now routinely generated by the Illumina platforms and provide output reports for multiple types of recombinant species using standardized formats. We demonstrate the utility and flexibility of ViReMa in different settings to report deletion events in simulated data from Flock House virus, copy-back RNA species in Sendai viruses, short duplication events in HIV, and virus-to-host recombination in an archaeal DNA virus.

Wild and weedy Hesperis matronalis hosts turnip mosaic virus across heterogeneous landscapes in upstate New York

Turnip mosaic virus (TuMV) is a widespread and economically important pathogen in agricultural crops and has the widest known host range in the virus family Potyviridae. While management of the virus and its aphid vectors in agricultural fields decreases virus incidence, many alternative wild hosts for TuMV may serve as source populations for crop infection through spillover. Over thirty years ago, research demonstrated that the introduced brassica, Dame's Rocket (Hesperis matronalis) hosts several viruses, including TuMV. Here, we use both enzyme-linked immunosorbent assays (ELISA) and next generation sequencing to document the frequent infection by TuMV of Dame's Rocket, which is common and widespread in disturbed areas around crop fields in upstate New York. Deep sequencing of multiple tissue types of symptomatic hosts indicate that the infection is systemic and causes diagnostic, visible symptoms. In a common garden experiment using host populations from across upstate New York, we found evidence for genetic tolerance to TuMV infection in H. matronalis. Field surveys show that TuMV prevalence varies across populations, but is generally higher in agricultural areas. Examining disease dynamics in this and other common alternative hosts will enhance our understanding of TuMV epidemiology and, more broadly, virus distribution in wild plants.

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.

Insight into aphid mediated Potato Virus Y transmission: A molecular to bioinformatics prospective

Potato, the world's most popular crop is reported to provide a food source for nearly a billion people. It is prone to a number of biotic stressors that affect yield and quality, out of which Potato Virus Y (PVY) occupies the top position. PVY can be transmitted mechanically and by sap-feeding aphid vectors. The application of insecticide causes an increase in the resistant vector population along with detrimental effects on the environment; genetic resistance and vector-virus control are the two core components for controlling the deadly PVY. Using transcriptomic tools together with differential gene expression and gene discovery, several loci and genes associated with PVY resistance have been widely identified. To combat this virus we must increase our understanding on the molecular response of the PVY-potato plant-aphid interaction and knowledge of genome organization, as well as the function of PVY encoded proteins, genetic diversity, the molecular aspects of PVY transmission by aphids, and transcriptome profiling of PVY infected potato cultivars. Techniques such as molecular and bioinformatics tools can identify and monitor virus transmission. Several studies have been conducted to understand the molecular basis of PVY resistance/susceptibility interactions and their impact on PVY epidemiology by studying the interrelationship between the virus, its vector, and the host plant. This review presents current knowledge of PVY transmission, epidemiology, genome organization, molecular to bioinformatics responses, and its effective management.

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.

Comparison of High Throughput Sequencing to Standard Protocols for Virus Detection in Berry Crops

We completed a comprehensive study comparing virus detection between high throughput sequencing (HTS) and standard protocols in 30 berry selections (12 Fragaria, 10 Vaccinium, and eight Rubus) with known virus profiles. The study examined temporal detection of viruses at four sampling times encompassing two growing seasons. Within the standard protocols, reverse transcription (RT) PCR proved better than biological indexing. Detection of known viruses by HTS and RT-PCR nearly mirrored each other. HTS provided superior detection compared with RT-PCR on a wide spectrum of variants and discovery of novel viruses. More importantly, in most cases in which the two protocols showed parallel virus detection, 11 viruses in 16 selections were not consistently detected by both methods at all sampling points. Based on these data, we propose a testing requirement of four sampling times over two growing seasons for berry and potentially other crops, to ensure that no virus remains undetected independent of titer, distribution, or other virus-virus or virus-host interactions.

Use of High-Throughput Sequencing and Two RNA Input Methods to Identify Viruses Infecting Tomato Crops

We used high-throughput sequencing to identify viruses on tomato samples showing virus-like symptoms. Samples were collected from crops in the Iberian Peninsula. Either total RNA or double-stranded RNA (dsRNA) were used as starting material to build the cDNA libraries. In total, seven virus species were identified, with pepino mosaic virus being the most abundant one. The dsRNA input provided better coverage and read depth but missed one virus species compared with the total RNA input. By performing in silico analyses, we determined a minimum sequencing depth per sample of 0.2 and 1.5 million reads for dsRNA and rRNA-depleted total RNA inputs, respectively, to detect even the less abundant viruses. Primers and TaqMan probes targeting conserved regions in the viral genomes were designed and/or used for virus detection; all viruses were detected by qRT-PCR/RT-PCR in individual samples, with all except one sample showing mixed infections. Three virus species (Olive latent virus 1, Lettuce ring necrosis virus and Tomato fruit blotch virus) are herein reported for the first time in tomato crops in Spain.

A Primer on the Analysis of High-Throughput Sequencing Data for Detection of Plant Viruses

High-throughput sequencing (HTS) technologies have become indispensable tools assisting plant virus diagnostics and research thanks to their ability to detect any plant virus in a sample without prior knowledge. As HTS technologies are heavily relying on bioinformatics analysis of the huge amount of generated sequences, it is of utmost importance that researchers can rely on efficient and reliable bioinformatic tools and can understand the principles, advantages, and disadvantages of the tools used. Here, we present a critical overview of the steps involved in HTS as employed for plant virus detection and virome characterization. We start from sample preparation and nucleic acid extraction as appropriate to the chosen HTS strategy, which is followed by basic data analysis requirements, an extensive overview of the in-depth data processing options, and taxonomic classification of viral sequences detected. By presenting the bioinformatic tools and a detailed overview of the consecutive steps that can be used to implement a well-structured HTS data analysis in an easy and accessible way, this paper is targeted at both beginners and expert scientists engaging in HTS plant virome projects.

Effect of RNA silencing suppression activity of chrysanthemum virus B p12 protein on small RNA species

Chrysanthemum virus B encodes a multifunctional p12 protein that acts as a transcriptional activator in the nucleus and as a suppressor of RNA silencing in the cytoplasm. Here, we investigated the impact of p12 on accumulation of major classes of small RNAs (sRNAs). The results show dramatic changes in the sRNA profiles characterised by an overall reduction in sRNA accumulation, changes in the pattern of size distribution of canonical siRNAs and in the ratio between sense and antisense strands, lower abundance of siRNAs with a U residue at the 5'-terminus, and changes in the expression of certain miRNAs, most of which were downregulated.

Mutations in virus-derived small RNAs

RNA viruses exist as populations of genome variants. Virus-infected plants accumulate 21–24 nucleotide small interfering RNAs (siRNAs) derived from viral RNA (virus-derived siRNAs) through gene silencing. This paper describes the profile of mutations in virus-derived siRNAs for three members of the family Potyviridae: Turnip mosaic virus (TuMV), Papaya ringspot virus (PRSV) and Wheat streak mosaic virus (WSMV). For TuMV in Arabidopsis thaliana, profiles were obtained for mechanically inoculated rosette leaves and systemically infected cauline leaves and inflorescence. Results are consistent with selection pressure on the viral genome imposed by local and systemic movement. By genetically removing gene silencing in the plant and silencing suppression in the virus, our results showed that antiviral gene silencing imposes selection in viral populations. Mutations in siRNAs derived from a PRSV coat protein transgene in the absence of virus replication showed the contribution of cellular RNA-dependent RNA polymerases to the generation of mutations in virus-derived siRNAs. Collectively, results are consistent with two sources of mutations in virus-derived siRNAs: viral RNA-dependent RNA polymerases responsible for virus replication and cellular RNA-dependent RNA polymerases responsible for gene silencing amplification.

Characterization of Potato Virus Y Isolates and Assessment of Nanopore Sequencing to Detect and Genotype Potato Viruses

Potato virus Y (PVY) is the most economically important virus infecting cultivated potato (Solanum tuberosum L.). Accurate diagnosis is crucial to regulate the trade of tubers and for the sanitary selection of plant material for propagation. However, high genetic diversity of PVY represents a challenge for the detection and classification of isolates. Here, the diversity of Irish PVY isolates from a germplasm collection and commercial sites was investigated using conventional molecular and serological techniques. Recombinant PVY isolates were prevalent, with PVYNTNa being the predominant genotype. In addition, we evaluated Nanopore sequencing to detect and reconstruct the whole genome sequence of four viruses (PVY, PVX, PVS, PLRV) and five PVY genotypes in a subset of eight potato plants. De novo assembly of Nanopore sequencing reads produced single contigs covering greater than 90% of the viral genome and sharing greater than 99.5% identity to the consensus sequences obtained with Illumina sequencing. Interestingly, single near full genome contigs were obtained for different isolates of PVY co-infecting the same plant. Mapping reads to available reference viral genomes enabled us to generate near complete genome sequences sharing greater than 99.90% identity to the Illumina-derived consensus. This is the first report describing the use of Oxford Nanopore's MinION to detect and genotype potato viruses. We reconstructed the genome of PVY and other RNA viruses; indicating the technologies potential for virus detection in potato production systems, and for the study of genetic diversity of highly heterogeneous viruses such as PVY.

The Potyviruses: An Evolutionary Synthesis Is Emerging

In this review, encouraged by the dictum of Theodosius Dobzhansky that "Nothing in biology makes sense except in the light of evolution", we outline the likely evolutionary pathways that have resulted in the observed similarities and differences of the extant molecules, biology, distribution, etc. of the potyvirids and, especially, its largest genus, the potyviruses. The potyvirids are a family of plant-infecting RNA-genome viruses. They had a single polyphyletic origin, and all share at least three of their genes (i.e., the helicase region of their CI protein, the RdRp region of their NIb protein and their coat protein) with other viruses which are otherwise unrelated. Potyvirids fall into 11 genera of which the potyviruses, the largest, include more than 150 distinct viruses found worldwide. The first potyvirus probably originated 15,000-30,000 years ago, in a Eurasian grass host, by acquiring crucial changes to its coat protein and HC-Pro protein, which enabled it to be transmitted by migrating host-seeking aphids. All potyviruses are aphid-borne and, in nature, infect discreet sets of monocotyledonous or eudicotyledonous angiosperms. All potyvirus genomes are under negative selection; the HC-Pro, CP, Nia, and NIb genes are most strongly selected, and the PIPO gene least, but there are overriding virus specific differences; for example, all turnip mosaic virus genes are more strongly conserved than those of potato virus Y. Estimates of dN/dS (ω) indicate whether potyvirus populations have been evolving as one or more subpopulations and could be used to help define species boundaries. Recombinants are common in many potyvirus populations (20%-64% in five examined), but recombination seems to be an uncommon speciation mechanism as, of 149 distinct potyviruses, only two were clear recombinants. Human activities, especially trade and farming, have fostered and spread both potyviruses and their aphid vectors throughout the world, especially over the past five centuries. The world distribution of potyviruses, especially those found on islands, indicates that potyviruses may be more frequently or effectively transmitted by seed than experimental tests suggest. Only two meta-genomic potyviruses have been recorded from animal samples, and both are probably contaminants.

Phytovirome Analysis of Wild Plant Populations: Comparison of Double-Stranded RNA and Virion-Associated Nucleic Acid Metagenomic Approaches

Metagenomic studies have indicated that the diversity of plant viruses was until recently far underestimated. As important components of ecosystems, there is a need to explore the diversity and richness of the viruses associated with plant populations and to understand the drivers shaping their diversity in space and time. Two viral sequence enrichment approaches, double-stranded RNA (dsRNA) and virion-associated nucleic acids (VANA), have been used and compared here for the description of the virome of complex plant pools representative of the most prevalent plant species in unmanaged and cultivated ecosystems. A novel bioinformatics strategy was used to assess viral richness not only at the family level but also by determining operational taxonomic units (OTU) following the clustering of conserved viral domains. A large viral diversity dominated by novel dsRNA viruses was detected in all sites, while a large between-site variability limited the ability to draw a clear conclusion on the impact of cultivation. A trend for a higher diversity of dsRNA viruses was nevertheless detected in unmanaged sites (118 versus 77 unique OTUs). The dsRNA-based approach consistently revealed a broader and more comprehensive diversity for RNA viruses than the VANA approach, whatever the assessment criterion. In addition, dissimilarity analyses indicated both approaches to be largely reproducible but not necessarily convergent. These findings illustrate features of phytoviromes in various ecosystems and a novel strategy for precise virus richness estimation. These results allow us to reason methodological choices in phytovirome studies and likely in other virome studies where RNA viruses are the focal taxa.IMPORTANCE There are today significant knowledge gaps on phytovirus populations and on the drivers impacting them but also on the comparative performance-methodological approaches for their study. We used and compared two viral sequence enrichment approaches, double-stranded RNAs (dsRNA) and virion-associated nucleic acids (VANA), for phytovirome description in complex pools representative of the most prevalent plant species in unmanaged and cultivated ecosystems. Viral richness was assessed by determining operational taxonomic units (OTU) following the clustering of conserved viral domains. There is some limited evidence of an impact of cultivation on viral populations. These results provide data allowing us to reason the methodological choices in virome studies. For researchers primarily interested in RNA viruses, the dsRNA approach is recommended because it consistently provided a more comprehensive description of the analyzed phytoviromes, but it understandably underrepresented DNA viruses and bacteriophages.

Multiomics analysis of tolerant interaction of potato with potato virus Y

Potato virus Y (PVY) is the most economically important viral pathogen of potato worldwide. Different potato cultivars react to the pathogen differently, resulting in resistant, tolerant or disease outcome of the interaction. Here we focus on tolerant interaction between potato cv. Désirée and PVYNTN. To capture the response in its full complexity, we analyzed the dynamic changes on multiple molecular levels, including transcriptomics, sRNAomics, degradomics, proteomics and hormonomics. The analysis was complemented by the measurements of viral accumulation, photosynthetic activity and phenotypisation of the symptoms. Besides cv. Désirée we also studied its transgenic counterpart depleted for the accumulation of salicylic acid (NahG-Désirée). This multiomics analysis provides better insights into the mechanisms leading to tolerant response of potato to viral infection and can be used as a base in further studies of plant immunity regulation.

From current knowledge to best practice: A primer on viral diagnostics using deep sequencing of virus-derived small interfering RNAs (vsiRNAs) in infected plants

Plants have evolved many defense strategies for combating viral infections. One major surveillance strategy adopted by them is manipulating viral sequences to generate distinct small RNA products via Dicer-like enzymes (DCL), and thereby restricting virus multiplication through the RNA interference (RNAi) mechanism. The power of high-throughput sequencing technologies, with diverse computational tools to handle small RNA sequencing (sRNA-Seq) data, bestows unprecedented opportunities to answer fundamental questions in plant virology. Here, we present some basic concepts of virus-derived, small interfering RNA (vsiRNA) biogenesis in plants, optimization strategies, caveats, and best practices for efficient discovery and diagnosis of known as well as novel plant viruses/viroids using deep sequencing of small RNA (sRNA) pools.

Potato virus Y; the Andean connection

Potato virus Y (PVY) causes disease in potatoes and other solanaceous crops. The appearance of its necrogenic strains in the 1980s made it the most economically important virus of potatoes. We report the isolation and genomic sequences of 32 Peruvian isolates of PVY which, together with 428 published PVY genomic sequences, gave an alignment of 460 sequences. Of these 190 (41%) were non-recombinant, and 162 of these provided a dated phylogeny, that corresponds well with the likely history of PVY, and show that PVY originated in South America which is where potatoes were first domesticated. The most basal divergences of the PVY population produced the N and C: O phylogroups; the origin of the N phylogroup is clearly Andean, but that of the O and C phylogroups is unknown, although they may have been first to establish in European crops. The current PVY population originated around 156 CE. PVY was probably first taken from South America to Europe in the 16th century in tubers. Most of the present PVY diversity emerged in the second half of the 19th century, after the Phytophthora infestans epidemics of the mid-19th century destroyed the European crop and stimulated potato breeding. Imported breeding lines were shared, and there was no quarantine. The early O population was joined later by N phylogroup isolates and their recombinants generated the R1 and R2 populations of damaging necrogenic strains. Our dating study has confirmed that human activity has dominated the phylodynamics of PVY for the last two millennia.

Immunocapture of virions with virus-specific antibodies prior to high-throughput sequencing effectively enriches for virus-specific sequences

Virus discovery based on high-throughput sequencing relies on enrichment for virus sequences prior to library preparation to achieve a sufficient number of viral reads. In general, preparations of double-stranded RNA or total RNA preparations treated to remove rRNA are used for sequence enrichment. We used virus-specific antibodies to immunocapture virions from plant sap to conduct cDNA synthesis, followed by library preparation and HTS. For the four potato viruses PLRV, PVY, PVA and PYV, template preparation by virion immunocapture provided a simpler and less expensive method than the enrichment of total RNA by ribosomal depletion. Specific enrichment of viral sequences without an intermediate amplification step was achieved, and this high coverage of sequences across the viral genomes was important to identify rare sequence variations. Using this approach, the first complete genome sequence of a potato yellowing virus isolate (PYV, DSMZ PV-0706) was determined in this study. PYV can be confidently assigned as a distinct species in the genus Ilarvirus.

High Throughput Sequencing For Plant Virus Detection and Discovery

Over the last decade, virologists have discovered an unprecedented number of viruses using high throughput sequencing (HTS), which led to the advancement of our knowledge on the diversity of viruses in nature, particularly unraveling the virome of many agricultural crops. However, these new virus discoveries have often widened the gaps in our understanding of virus biology; the forefront of which is the actual role of a new virus in disease, if any. Yet, when used critically in etiological studies, HTS is a powerful tool to establish disease causality between the virus and its host. Conversely, with globalization, movement of plant material is increasingly more common and often a point of dispute between countries. HTS could potentially resolve these issues given its capacity to detect and discover. Although many pipelines are available for plant virus discovery, all share a common backbone. A description of the process of plant virus detection and discovery from HTS data are presented, providing a summary of the different pipelines available for scientists' utility in their research.

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.

Small RNA-Omics for Plant Virus Identification, Virome Reconstruction, and Antiviral Defense Characterization

RNA interference (RNAi)-based antiviral defense generates small interfering RNAs that represent the entire genome sequences of both RNA and DNA viruses as well as viroids and viral satellites. Therefore, deep sequencing and bioinformatics analysis of small RNA population (small RNA-ome) allows not only for universal virus detection and genome reconstruction but also for complete virome reconstruction in mixed infections. Viral infections (like other stress factors) can also perturb the RNAi and gene silencing pathways regulating endogenous gene expression and repressing transposons and host genome-integrated endogenous viral elements which can potentially be released from the genome and contribute to disease. This review describes the application of small RNA-omics for virus detection, virome reconstruction and antiviral defense characterization in cultivated and non-cultivated plants. Reviewing available evidence from a large and ever growing number of studies of naturally or experimentally infected hosts revealed that all families of land plant viruses, their satellites and viroids spawn characteristic small RNAs which can be assembled into contigs of sufficient length for virus, satellite or viroid identification and for exhaustive reconstruction of complex viromes. Moreover, the small RNA size, polarity and hotspot profiles reflect virome interactions with the plant RNAi machinery and allow to distinguish between silent endogenous viral elements and their replicating episomal counterparts. Models for the biogenesis and functions of small interfering RNAs derived from all types of RNA and DNA viruses, satellites and viroids as well as endogenous viral elements are presented and discussed.

Characterization of maize translational responses to sugarcane mosaic virus infection

Sugarcane mosaic virus (SCMV) frequently causes dramatic losses in maize production as the main pathogen of maize dwarf mosaic disease. It is important to understand the translational responses in maize to SCMV infection since viruses have to recruit host translation apparatus to express their proteins. However, due to technical limitations, research on virus translation lags far behind that on transcription. Here, we studied the relationship between systemic symptom expression and virus accumulation and found that both SCMV RNA and proteins accumulated rapidly during the systemic infection process in which varying degrees of chlorosis to mosaic symptoms developed on non-inoculated leaves. In addition, we applied ribosome profiling, which couples polysomal mRNA isolation with high-throughput sequencing, on the symptomatic leaves infected with SCMV to unravel the translational responses of maize to viral infection on a genome-wide scale. The results showed that only the genomic positive-stranded RNA of SCMV was involved in translation, and SCMV only occupied a small amount of translational resources of host plant at the early stage of infection. Further analyses on a global gene expression and gene ontology (GO) enrichment revealed that photosynthesis and metabolism were dramatically repressed at both transcriptional and translational levels. Altogether, our results laid a foundation for dissecting the molecular mechanism of plant translational responses to viral infection.

Small RNA-Omics for Virome Reconstruction and Antiviral Defense Characterization in Mixed Infections of Cultivated Solanum Plants

In plants, RNA silencing-based antiviral defense generates viral small RNAs (sRNAs) faithfully representing the viral genomes. We employed sRNA sequencing and bioinformatics (sRNA-omics) to characterize antiviral defense and to reconstruct the full genomic sequences and their variants in the evolving viral quasispecies in cultivated solanaceous plants carrying mixed infections. In naturally infected Solanum tuberosum (potato), one case study revealed a virome comprising Potato virus Y (genus Potyvirus) and Potato virus X (genus Potexvirus), which was reconstructed by de novo-assembling separate genome-size sRNA contigs. Another case study revealed a virome comprising NTN and O strains of Potato virus Y, whose sRNAs assembled in chimeric contigs, which could be disentangled on the basis of reference genome sequences. Both viromes were stable in vegetative potato progeny. In a cross-protection trial of Solanum lycopersicum (tomato), the supposedly protective mild strain CH2 of Pepino mosaic virus (genus Potexvirus) was tested for protection against strain LP of the same virus. Reciprocal mechanical inoculations eventually resulted in co-infection of all individual plants with CH2 and LP strains, reconstructed as separate sRNA contigs. LP invasions into CH2-preinfected plants and vice versa were accompanied by alterations of consensus genome sequences in viral quasispecies, indicating a potential risk of cross-protection measures. Additionally, the study also revealed, by reconstruction from sRNAs, the presence of the mechanically nontransmissible Southern tomato virus (genus Amalgavirus) in some plants. Our in-depth analysis of sRNA sizes, 5'-nucleotide frequencies and hotspot maps revealed similarities in sRNA-generating mechanisms in potato and tomato, differential silencing responses to virome components and potential for sRNA-directed cross-targeting between viral strains which could not, however, prevent the formation of stable viromes.

Monolith Chromatography as Sample Preparation Step in Virome Studies of Water Samples

Viruses exist in aquatic media and many of them use this media as transmission route. Next-generation sequencing (NGS) technologies have opened new doors in virus research, allowing also to reveal a hidden diversity of viral species in aquatic environments. Not surprisingly, many of the newly discovered viruses are found in environmental fresh and marine waters. One of the problems in virome research can be the low amount of viral nucleic acids present in the sample in contrast to the background ones (host, eukaryotic, prokaryotic, environmental). Therefore, virus enrichment prior to NGS is necessary in many cases. In water samples, an added problem resides in the low concentration of viruses typically present in aquatic media. Different concentration strategies have been used to overcome such limitations. CIM monoliths are a new generation of chromatographic supports that due to their particular structural characteristics are very efficient in concentration and purification of viruses. In this chapter, we describe the use of CIM monolithic chromatography for sample preparation step in NGS studies targeting viruses in fresh or marine water. The step-by-step protocol will include a case study where CIM concentration was used to study the virome of a wastewater sample using NGS.

Salicylic Acid Perturbs sRNA-Gibberellin Regulatory Network in Immune Response of Potato to Potato virus Y Infection

Potato virus Y is the most economically important potato viral pathogen. We aimed at unraveling the roles of small RNAs (sRNAs) in the complex immune signaling network controlling the establishment of tolerant response of potato cv. Désirée to the virus. We constructed a sRNA regulatory network connecting sRNAs and their targets to link sRNA level responses to physiological processes. We discovered an interesting novel sRNAs-gibberellin regulatory circuit being activated as early as 3 days post inoculation (dpi) before viral multiplication can be detected. Two endogenous sRNAs, miR167 and phasiRNA931 were predicted to regulate gibberellin biosynthesis genes GA20-oxidase and GA3-oxidase. The increased expression of phasiRNA931 was also reflected in decreased levels of GA3-oxidase transcripts. Moreover, decreased concentration of gibberellin confirmed this regulation. The functional relation between lower activity of gibberellin signaling and reduced disease severity was previously confirmed in Arabidopsis-virus interaction using knockout mutants. We further showed that this regulation is salicylic acid-dependent as the response of sRNA network was attenuated in salicylic acid-depleted transgenic counterpart NahG-Désirée expressing severe disease symptoms. Besides downregulation of gibberellin signaling, regulation of immune receptor transcripts by miR6022 as well as upregulation of miR164, miR167, miR169, miR171, miR319, miR390, and miR393 in tolerant Désirée, revealed striking similarities to responses observed in mutualistic symbiotic interactions. The intertwining of different regulatory networks revealed, shows how developmental signaling, disease symptom development, and stress signaling can be balanced.

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).

Differential Characteristics of Viral siRNAs between Leaves and Roots of Wheat Plants Naturally Infected with Wheat Yellow Mosaic Virus, a Soil-Borne Virus

RNA silencing is an important innate antiviral defense in plants. Soil-borne plant viruses naturally infect roots via soil-inhabiting vectors, but it is unclear how antiviral RNA silencing responds to virus infection in this particular tissue. In this study, viral small interfering RNA (siRNA) profiles from leaves and roots of wheat plants naturally infected with a soil-borne virus, wheat yellow mosaic virus (WYMV, genus Bymovirus), were analyzed by deep sequencing. WYMV siRNAs were much more abundant in roots than leaves, which was positively correlated with the accumulation of viral RNA. WYMV siRNAs in leaves and roots were predominantly 21- and 22-nt long and equally derived from the positive- and negative-strands of the viral genome. WYMV siRNAs from leaves and roots differed in distribution pattern along the viral genome. Interestingly, compared to siRNAs from leaves (and most other reports), those from roots obviously had a lower A/U bias at the 5'-terminal nucleotide. Moreover, the expression of Dicer-like genes upon WYMV infection were differently regulated between leaves and roots. Our data suggest that RNA silencing in roots may operate differently than in leaves against soil-borne virus invasion.

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

RNA viruses are one of the fastest-evolving biological entities. Within their hosts, they exist as genetically diverse populations (i.e., viral mutant swarms), which are sculpted by different evolutionary mechanisms, such as mutation, natural selection, and genetic drift, and also the interactions between genetic variants within the mutant swarms. To elucidate the mechanisms that modulate the population diversity of an important plant-pathogenic virus, we performed evolution experiments with Potato virus Y (PVY) in potato genotypes that differ in their defense response against the virus. Using deep sequencing of small RNAs, we followed the temporal dynamics of standing and newly generated variations in the evolving viral lineages. A time-sampled approach allowed us to (i) reconstruct theoretical haplotypes in the starting population by using clustering of single nucleotide polymorphisms' trajectories and (ii) use quantitative population genetics approaches to estimate the contribution of selection and genetic drift, and their interplay, to the evolution of the virus. We detected imprints of strong selective sweeps and narrow genetic bottlenecks, followed by the shift in frequency of selected haplotypes. Comparison of patterns of viral evolution in differently susceptible host genotypes indicated possible diversifying evolution of PVY in the less-susceptible host (efficient in the accumulation of salicylic acid).IMPORTANCE High diversity of within-host populations of RNA viruses is an important aspect of their biology, since they represent a reservoir of genetic variants, which can enable quick adaptation of viruses to a changing environment. This study focuses on an important plant virus, Potato virus Y, and describes, at high resolution, temporal changes in the structure of viral populations within different potato genotypes. A novel and easy-to-implement computational approach was established to cluster single nucleotide polymorphisms into viral haplotypes from very short sequencing reads. During the experiment, a shift in the frequency of selected viral haplotypes was observed after a narrow genetic bottleneck, indicating an important role of the genetic drift in the evolution of the virus. On the other hand, a possible case of diversifying selection of the virus was observed in less susceptible host genotypes.

The pepper virome: natural co-infection of diverse viruses and their quasispecies

BACKGROUND

The co-infection of diverse viruses in a host plant is common; however, little is known about viral populations and their quasispecies in the host.

RESULTS

Here, we report the first pepper viromes that were co-infected by different types of viral genomes. The pepper viromes are dominated by geminivirus DNA-A followed by a novel carlavirus referred to as Pepper virus A. The two pepper cultivars share similar viral populations and replications. However, the quasispecies for double-stranded RNA virus and two satellite DNAs were heterogeneous and homogenous in susceptible and resistant cultivars, respectively, indicating the quasispecies of an individual virus depends on the host.

CONCLUSIONS

Taken together, we provide the first evidence that the host plant resistant to viruses has an unrevealed antiviral system, affecting viral quasispecies, not replication.

Diversity, Distribution, and Evolution of Tomato Viruses in China Uncovered by Small RNA Sequencing

Tomato is a major vegetable crop that has tremendous popularity. However, viral disease is still a major factor limiting tomato production. Here, we report the tomato virome identified through sequencing small RNAs of 170 field-grown samples collected in China. A total of 22 viruses were identified, including both well-documented and newly detected viruses. The tomato viral community is dominated by a few species, and they exhibit polymorphisms and recombination in the genomes with cold spots and hot spots. Most samples were coinfected by multiple viruses, and the majority of identified viruses are positive-sense single-stranded RNA viruses. Evolutionary analysis of one of the most dominant tomato viruses, Tomato yellow leaf curl virus (TYLCV), predicts its origin and the time back to its most recent common ancestor. The broadly sampled data have enabled us to identify several unreported viruses in tomato, including a completely new virus, which has a genome of ∼13.4 kb and groups with aphid-transmitted viruses in the genus Cytorhabdovirus Although both DNA and RNA viruses can trigger the biogenesis of virus-derived small interfering RNAs (vsiRNAs), we show that features such as length distribution, paired distance, and base selection bias of vsiRNA sequences reflect different plant Dicer-like proteins and Argonautes involved in vsiRNA biogenesis. Collectively, this study offers insights into host-virus interaction in tomato and provides valuable information to facilitate the management of viral diseases.IMPORTANCE Tomato is an important source of micronutrients in the human diet and is extensively consumed around the world. Virus is among the major constraints on tomato production. Categorizing virus species that are capable of infecting tomato and understanding their diversity and evolution are challenging due to difficulties in detecting such fast-evolving biological entities. Here, we report the landscape of the tomato virome in China, the leading country in tomato production. We identified dozens of viruses present in tomato, including both well-documented and completely new viruses. Some newly emerged viruses in tomato were found to spread fast, and therefore, prompt attention is needed to control them. Moreover, we show that the virus genomes exhibit considerable degree of polymorphisms and recombination, and the virus-derived small interfering RNA (vsiRNA) sequences indicate distinct vsiRNA biogenesis mechanisms for different viruses. The Chinese tomato virome that we developed provides valuable information to facilitate the management of tomato viral diseases.

Pathogenic seedborne viruses are rare but Phaseolus vulgaris endornaviruses are common in bean varieties grown in Nicaragua and Tanzania

Common bean (Phaseolus vulgaris) is an annual grain legume that was domesticated in Mesoamerica (Central America) and the Andes. It is currently grown widely also on other continents including Africa. We surveyed seedborne viruses in new common bean varieties introduced to Nicaragua (Central America) and in landraces and improved varieties grown in Tanzania (eastern Africa). Bean seeds, harvested from Nicaragua and Tanzania, were grown in insect-controlled greenhouse or screenhouse, respectively, to obtain leaf material for virus testing. Equal amounts of total RNA from different samples were pooled (30-36 samples per pool), and small RNAs were deep-sequenced (Illumina). Assembly of the reads (21-24 nt) to contiguous sequences and searches for homologous viral sequences in databases revealed Phaseolus vulgaris endornavirus 1 (PvEV-1) and PvEV-2 in the bean varieties in Nicaragua and Tanzania. These viruses are not known to cause symptoms in common bean and are considered non-pathogenic. The small-RNA reads from each pool of samples were mapped to the previously characterized complete PvEV-1 and PvEV-2 sequences (genome lengths ca. 14 kb and 15 kb, respectively). Coverage of the viral genomes was 87.9-99.9%, depending on the pool. Coverage per nucleotide ranged from 5 to 471, confirming virus identification. PvEV-1 and PvEV-2 are known to occur in Phaseolus spp. in Central America, but there is little previous information about their occurrence in Nicaragua, and no information about occurrence in Africa. Aside from Cowpea mild mosaic virus detected in bean plants grown from been seeds harvested from one region in Tanzania, no other pathogenic seedborne viruses were detected. The low incidence of infections caused by pathogenic viruses transmitted via bean seeds may be attributable to new, virus-resistant CB varieties released by breeding programs in Nicaragua and Tanzania.

Viral small RNAs reveal the genomic variations of three grapevine vein clearing virus quasispecies populations

Viral small RNAs (vsRNAs) include viral small interfering RNAs (vsiRNAs) that are initiators and products of RNA silencing, and small RNAs that are derived from viral RNAs with function still unknown. Sequencing of vsRNAs allows assembling of viral genomes and revelation of viral population variations at genomic levels. Grapevine vein clearing virus (GVCV) is a new member of the family Caulimoviridae whose DNA genome is replicated by reverse transcription of pre-genomic RNA molecules. In this short report, three genomic sequences of GVCV were assembled from vsRNAs that were isolated and sequenced from three individual grapevines in commercial vineyards and compared to the GVCV-CHA reference genome. Profiles of single nucleotide polymorphism among three viral populations indicated a closer relatedness between two populations in different grape cultivars at the same location than those in the same grape cultivar at different locations, suggesting the spread of GVCV populations among vineyards of close proximity. Classic types of vsiRNAs (21-nt, 22-nt, and 24-nt) were found in the three GVCV vsiRNA populations, but these did not produce alignment hotspots on the GVCV-CHA reference genome. The number of 36-nt reads is the highest among vsRNAs, the role of these vsRNAs remains unclear. The analysis of vsRNAs provides a first holistic picture of genomic variations among GVCV viral quasispecies populations that help monitor epidemics and evolution of GVCV populations, an emerging virus that is becoming a threat to grape production in the Midwest region of the USA.

Deep sequencing for discovery and evolutionary analysis of plant viruses

The advent of next generation sequencing (NGS), or deep sequencing, has allowed great advances to be made in discovery, diagnostics, and evolutionary studies in plant viruses. Various methods have been used for enrichment for virus-specific nucleic acids, each of which have some drawbacks. Many novel viruses have been discovered in plants by NGS technologies, and there is a good deal of promise for more comprehensive studies in virus evolution. However, each aspect of using NGS has its caveats that need to be considered, and there is still a need for better tools of analysis, as well as method for validation of sequence variation.

Different Virus-Derived siRNAs Profiles between Leaves and Fruits in Cucumber Green Mottle Mosaic Virus-Infected Lagenaria siceraria Plants

RNA silencing is an evolutionarily conserved antiviral mechanism, through which virus-derived small interfering RNAs (vsiRNAs) playing roles in host antiviral defense are produced in virus-infected plant. Deep sequencing technology has revolutionized the study on the interaction between virus and plant host through the analysis of vsiRNAs profile. However, comparison of vsiRNA profiles in different tissues from a same host plant has been rarely reported. In this study, the profiles of vsiRNAs from leaves and fruits of Lagenaria siceraria plants infected with Cucumber green mottle mosaic virus (CGMMV) were comprehensively characterized and compared. Many more vsiRNAs were present in infected leaves than in fruits. vsiRNAs from both leaves and fruits were mostly 21- and 22-nt in size as previously described in other virus-infected plants. Interestingly, vsiRNAs were predominantly produced from the viral positive strand RNAs in infected leaves, whereas in infected fruits they were derived equally from the positive and negative strands. Many leaf-specific positive vsiRNAs with lengths of 21-nt (2058) or 22-nt (3996) were identified but only six (21-nt) and one (22-nt) positive vsiRNAs were found to be specific to fruits. vsiRNAs hotspots were only present in the 5'-terminal and 3'-terminal of viral positive strand in fruits, while multiple hotspots were identified in leaves. Differences in GC content and 5'-terminal nucleotide of vsiRNAs were also observed in the two organs. To our knowledge, this provides the first high-resolution comparison of vsiRNA profiles between different tissues of the same host plant.

Host-specific accumulation and temperature effects on the generation of dimeric viral RNA species derived from the S-RNA of members of the Tospovirus genus

Polygonum ringspot virus (PolRSV) is a recently characterized Tospovirus reported in Italy. Northern blot analyses of PolRSV infections in Nicotiana benthamiana and tomato plants showed that a viral RNA species with nearly twice the length of the Small genomic RNA (S-RNA) accumulated abundantly in the former host, but was not detected in the latter. Additional assays confirmed that biogenesis of this novel RNA species was common to all PolRSV isolates tested and also to an isolate of Tomato spotted wilt virus (TSWV). Given its size, we hypothesized that the novel RNA species was a dimer molecule and we confirmed this hypothesis by RNA sequencing (RNAseq) analysis and reverse transcription (RT)-PCR of putative predicted dimer junction sites in RNA extracts of N. benthamiana challenged with PolRSV isolates Plg6 and Plg13/2. We also confirmed that these molecules are derived from head-to-tail dimers and often contain deletions at their junction sites. We named these novel molecules imperfect dimer RNAs (IMPD-RNAs). PolRSV IMPD-RNAs systemic accumulation in a range of host plants was restricted to N. benthamiana and Nicotiana occidentalis. Notably, IMPD-RNAs accumulation was modulated by temperature and their generation was restricted to late stages of systemic infection (12 days post-inoculation) in N. benthamiana. Differently from all other PolRSV isolates used in this study, Plg13/2 generated more IMPD-RNAs coupled with low amounts of genomic S-RNA and maintained them even at 18 °C, besides having lost the ability to infect tomato plants. This is the first characterization of S-RNA dimers for Tospovirus, and of occurrence of dimers of genomic segments at the whole organism level for Bunyaviridae.

MISIS-2: A bioinformatics tool for in-depth analysis of small RNAs and representation of consensus master genome in viral quasispecies

In most eukaryotes, small RNA (sRNA) molecules such as miRNAs, siRNAs and piRNAs regulate gene expression and repress transposons and viruses. AGO/PIWI family proteins sort functional sRNAs based on size, 5'-nucleotide and other sequence features. In plants and some animals, viral sRNAs are extremely diverse and cover the entire viral genome sequences, which allows for de novo reconstruction of a complete viral genome by deep sequencing and bioinformatics analysis of viral sRNAs. Previously, we have developed a tool MISIS to view and analyze sRNA maps of viruses and cellular genome regions which spawn multiple sRNAs. Here we describe a new release of MISIS, MISIS-2, which enables to determine and visualize a consensus sequence and count sRNAs of any chosen sizes and 5'-terminal nucleotide identities. Furthermore we demonstrate the utility of MISIS-2 for identification of single nucleotide polymorphisms (SNPs) at each position of a reference sequence and reconstruction of a consensus master genome in evolving viral quasispecies. MISIS-2 is a Java standalone program. It is freely available along with the source code at the website http://www.fasteris.com/apps.

Using Small RNA Deep Sequencing Data to Detect Human Viruses

Small RNA sequencing (sRNA-seq) can be used to detect viruses in infected hosts without the necessity to have any prior knowledge or specialized sample preparation. The sRNA-seq method was initially used for viral detection and identification in plants and then in invertebrates and fungi. However, it is still controversial to use sRNA-seq in the detection of mammalian or human viruses. In this study, we used 931 sRNA-seq runs of data from the NCBI SRA database to detect and identify viruses in human cells or tissues, particularly from some clinical samples. Six viruses including HPV-18, HBV, HCV, HIV-1, SMRV, and EBV were detected from 36 runs of data. Four viruses were consistent with the annotations from the previous studies. HIV-1 was found in clinical samples without the HIV-positive reports, and SMRV was found in Diffuse Large B-Cell Lymphoma cells for the first time. In conclusion, these results suggest the sRNA-seq can be used to detect viruses in mammals and humans.

Biogenesis, Function, and Applications of Virus-Derived Small RNAs in Plants

RNA silencing, an evolutionarily conserved and sequence-specific gene-inactivation system, has a pivotal role in antiviral defense in most eukaryotic organisms. In plants, a class of exogenous small RNAs (sRNAs) originating from the infecting virus called virus-derived small interfering RNAs (vsiRNAs) are predominantly responsible for RNA silencing-mediated antiviral immunity. Nowadays, the process of vsiRNA formation and the role of vsiRNAs in plant viral defense have been revealed through deep sequencing of sRNAs and diverse genetic analysis. The biogenesis of vsiRNAs is analogous to that of endogenous sRNAs, which require diverse essential components including dicer-like (DCL), argonaute (AGO), and RNA-dependent RNA polymerase (RDR) proteins. vsiRNAs trigger antiviral defense through post-transcriptional gene silencing (PTGS) or transcriptional gene silencing (TGS) of viral RNA, and they hijack the host RNA silencing system to target complementary host transcripts. Additionally, several applications that take advantage of the current knowledge of vsiRNAs research are being used, such as breeding antiviral plants through genetic engineering technology, reconstructing of viral genomes, and surveying viral ecology and populations. Here, we will provide an overview of vsiRNA pathways, with a primary focus on the advances in vsiRNA biogenesis and function, and discuss their potential applications as well as the future challenges in vsiRNAs research.