Deep-sequencing of plant viral small RNAs reveals effective and widespread targeting of viral genomes

The plant siRNA landscape

Whereas micro (mi)RNAs are considered the clean, noble side of the small RNA world, small interfering (si)RNAs are often seen as a noisy set of molecules whose barbarian acronyms reflect a large diversity of often elusive origins and functions. Twenty-five years after their discovery in plants, however, new classes of siRNAs are still being identified, sometimes in discrete tissues or at particular developmental stages, making the plant siRNA world substantially more complex and subtle than originally anticipated. Focusing primarily on the model Arabidopsis, we review here the plant siRNA landscape, including transposable elements (TE)-derived siRNAs, a vast array of non-TE-derived endogenous siRNAs, as well as exogenous siRNAs produced in response to invading nucleic acids such as viruses or transgenes. We primarily emphasize the extraordinary sophistication and diversity of their biogenesis and, secondarily, the variety of their known or presumed functions, including via non-cell autonomous activities, in the sporophyte, gametophyte, and shortly after fertilization.

Enhanced Susceptibility to Tomato Chlorosis Virus (ToCV) in Hsp90- and Sgt1-Silenced Plants: Insights from Gene Expression Dynamics

The emerging whitefly-transmitted crinivirus tomato chlorosis virus (ToCV) causes substantial economic losses by inducing yellow leaf disorder in tomato crops. This study explores potential resistance mechanisms by examining early-stage molecular responses to ToCV. A time-course transcriptome analysis compared naïve, mock, and ToCV-infected plants at 2, 7, and 14 days post-infection (dpi). Gene expression changes were most notable at 2 and 14 dpi, likely corresponding to whitefly feeding and viral infection. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed key genes and pathways associated with ToCV infection, including those related to plant immunity, flavonoid and steroid biosynthesis, photosynthesis, and hormone signaling. Additionally, virus-derived small interfering RNAs (vsRNAs) originating from ToCV predominantly came from RNA2 and were 22 nucleotides in length. Furthermore, two genes involved in plant immunity, Hsp90 (heat shock protein 90) and its co-chaperone Sgt1 (suppressor of the G2 allele of Skp1) were targeted through viral-induced gene silencing (VIGS), showing a potential contribution to basal resistance against viral infections since their reduction correlated with increased ToCV accumulation. This study provides insights into tomato plant responses to ToCV, with potential implications for developing effective disease control strategies.

A Jasmonic Acid-Related Mechanism Affects ARGONAUTE5 Expression and Antiviral Defense Against Potato Virus X in Arabidopsis thaliana

During virus infection, Argonaute (AGO) proteins bind to Dicer-produced virus small interfering RNAs and target viral RNA based on sequence complementarity, thereby limiting virus proliferation. The Arabidopsis AGO2 protein is important for resistance to multiple viruses, including potato virus X (PVX). In addition, AGO5 is important in systemic defense against PVX. Normally AGO5 is expressed only in reproductive tissues, and its induction by virus infection is thought to be important for its participation in antiviral defense. However, it is unclear what mechanisms induce AGO5 expression in response to virus infection. Here, we show that dde2-2, a mutant compromised in jasmonic acid (JA) biosynthesis, displays constitutive upregulation of AGO5. This mutant also showed increased resistance to PVX and this resistance was dependent on a functional AGO5 gene. Furthermore, methyl jasmonate treatment ablated AGO5 expression in leaves during virus infection and resulted in increased susceptibility to virus. Our results further support a role for AGO5 in antiviral RNA silencing and a negative regulation by JA, a plant hormone associated with defense against plant-feeding arthropods, which are often the vectors of plant viruses. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Infection Defects of RNA and DNA Viruses Induced by Antiviral RNA Interference

Immune recognition of viral genome-derived double-stranded RNA (dsRNA) molecules and their subsequent processing into small interfering RNAs (siRNAs) in plants, invertebrates, and mammals trigger specific antiviral immunity known as antiviral RNA interference (RNAi). Immune sensing of viral dsRNA is sequence-independent, and most regions of viral RNAs are targeted by virus-derived siRNAs which extensively overlap in sequence. Thus, the high mutation rates of viruses do not drive immune escape from antiviral RNAi, in contrast to other mechanisms involving specific virus recognition by host immune proteins such as antibodies and resistance (R) proteins in mammals and plants, respectively. Instead, viruses actively suppress antiviral RNAi at various key steps with a group of proteins known as viral suppressors of RNAi (VSRs). Some VSRs are so effective in virus counter-defense that potent inhibition of virus infection by antiviral RNAi is undetectable unless the cognate VSR is rendered nonexpressing or nonfunctional. Since viral proteins are often multifunctional, resistance phenotypes of antiviral RNAi are accurately defined by those infection defects of VSR-deletion mutant viruses that are efficiently rescued by host deficiency in antiviral RNAi. Here, we review and discuss in vivo infection defects of VSR-deficient RNA and DNA viruses resulting from the actions of host antiviral RNAi in model systems.

Divergent RNA viruses infecting sea lice, major ectoparasites of fish

Sea lice, the major ectoparasites of fish, have significant economic impacts on wild and farmed finfish, and have been implicated in the decline of wild salmon populations. As blood-feeding arthropods, sea lice may also be reservoirs for viruses infecting fish. However, except for two groups of negative-strand RNA viruses within the order Mononegavirales, nothing is known about viruses of sea lice. Here, we used transcriptomic data from three key species of sea lice (Lepeophtheirus salmonis, Caligus clemensi, and Caligus rogercresseyi) to identify 32 previously unknown RNA viruses. The viruses encompassed all the existing phyla of RNA viruses, with many placed in deeply branching lineages that likely represent new families and genera. Importantly, the presence of canonical virus-derived small interfering RNAs (viRNAs) indicates that most of these viruses infect sea lice, even though in some cases their closest classified relatives are only known to infect plants or fungi. We also identified both viRNAs and PIWI-interacting RNAs (piRNAs) from sequences of a bunya-like and two qin-like viruses in C. rogercresseyi. Our analyses showed that most of the viruses found in C. rogercresseyi occurred in multiple life stages, spanning from planktonic to parasitic stages. Phylogenetic analysis revealed that many of the viruses infecting sea lice were closely related to those that infect a wide array of eukaryotes with which arthropods associate, including fungi and parasitic tapeworms, implying that over evolutionary time there has been cross-phylum and cross-kingdom switching of viruses between arthropods and other eukaryotes. Overall, this study greatly expands our view of virus diversity in crustaceans, identifies viruses that infect and replicate in sea lice, and provides evidence that over evolutionary time, viruses have switched between arthropods and eukaryotic hosts in other phyla and kingdoms.

Viromic approach reveals differences in the composition, diversity and relative abundance of pumpkin viruses across main growing regions of China

China is the leading country for pumpkin production in the world. As other cucurbits, diseases caused by viruses are among the serious threats to pumpkin production, but our knowledge on the virus species infecting pumpkin plants is fragmentary. To understand the occurrence of viral diseases on pumpkin, we determined the geographical distribution characteristics, relative abundance and evolutionary relationship of pumpkin infected viruses by meta-transcriptome sequencing (RNA-seq) and viromic analysis of 159 samples exhibited typical viral symptoms collected across China in this study. Totally, 11 known and 3 new viruses were identified. Interestingly, 3 new viruses identified in this study should be positive-sense single-stranded RNA virus whose hosts are prokaryotes. The viruses identified in different sampling locations exhibit significant variations in term of virus species and relative abundance. Here, the results provide valuable information for understanding the virus species and their diversity in cultivated pumpkin across major growing regions of China.

Looking beyond Virus Detection in RNA Sequencing Data: Lessons Learned from a Community-Based Effort to Detect Cellular Plant Pathogens and Pests

High-throughput sequencing (HTS), more specifically RNA sequencing of plant tissues, has become an indispensable tool for plant virologists to detect and identify plant viruses. During the data analysis step, plant virologists typically compare the obtained sequences to reference virus databases. In this way, they are neglecting sequences without homologies to viruses, which usually represent the majority of sequencing reads. We hypothesized that traces of other pathogens might be detected in this unused sequence data. In the present study, our goal was to investigate whether total RNA-seq data, as generated for plant virus detection, is also suitable for the detection of other plant pathogens and pests. As proof of concept, we first analyzed RNA-seq datasets of plant materials with confirmed infections by cellular pathogens in order to check whether these non-viral pathogens could be easily detected in the data. Next, we set up a community effort to re-analyze existing Illumina RNA-seq datasets used for virus detection to check for the potential presence of non-viral pathogens or pests. In total, 101 datasets from 15 participants derived from 51 different plant species were re-analyzed, of which 37 were selected for subsequent in-depth analyses. In 29 of the 37 selected samples (78%), we found convincing traces of non-viral plant pathogens or pests. The organisms most frequently detected in this way were fungi (15/37 datasets), followed by insects (13/37) and mites (9/37). The presence of some of the detected pathogens was confirmed by independent (q)PCRs analyses. After communicating the results, 6 out of the 15 participants indicated that they were unaware of the possible presence of these pathogens in their sample(s). All participants indicated that they would broaden the scope of their bioinformatic analyses in future studies and thus check for the presence of non-viral pathogens. In conclusion, we show that it is possible to detect non-viral pathogens or pests from total RNA-seq datasets, in this case primarily fungi, insects, and mites. With this study, we hope to raise awareness among plant virologists that their data might be useful for fellow plant pathologists in other disciplines (mycology, entomology, bacteriology) as well.

Transgene Silencing, RNA Interference, and the Antiviral Defense Mechanism Directed by Small Interfering RNAs

One important discovery in plant pathology over recent decades is the natural antiviral defense mechanism mediated by RNA interference (RNAi). In antiviral RNAi, virus infection triggers Dicer processing of virus-specific double-stranded RNA into small interfering RNAs (siRNAs). Frequently, further amplified by host enzyme and cofactors, these virus-derived siRNAs direct specific virus clearance in an Argonaute protein-containing effector complex. The siRNAs derived from viruses and viroids accumulate to very high levels during infection. Because they overlap extensively in nucleotide sequence, this allows for deep sequencing and bioinformatics assembly of total small RNAs for rapid discovery and identification of viruses and viroids. Antiviral RNAi acts as the primary defense mechanism against both RNA and DNA viruses in plants, yet viruses still successfully infect plants. They do so because all currently recognized plant viruses combat the RNAi response by encoding at least one protein as a viral suppressor of RNAi (VSR) required for infection, even though plant viruses have small genome sizes with a limited coding capacity. This review article will recapitulate the key findings that have revealed the genetic pathway for the biogenesis and antiviral activity of viral siRNAs and the specific role of VSRs in infection by antiviral RNAi suppression. Moreover, early pioneering studies on transgene silencing, RNAi, and virus-plant/virus-virus interactions paved the road to the discovery of antiviral RNAi.

Small RNA Profiling of Cucurbit Yellow Stunting Disorder Virus from Susceptible and Tolerant Squash (Cucurbita pepo) Lines

RNA silencing is a crucial mechanism of the antiviral immunity system in plants. Small RNAs guide Argonaut proteins to target viral RNA or DNA, preventing virus accumulation. Small RNA profiles in Cucurbita pepo line PI 420328 with tolerance to cucurbit yellow stunting disorder virus (CYSDV) were compared with those in Gold Star, a susceptible cultivar. The lower CYSDV symptom severity in PI 420328 correlated with lower virus titers and fewer sRNAs derived from CYSDV (vsRNA) compared to Gold Star. Elevated levels of 21- and 22-nucleotide (nt) size class vsRNAs were observed in PI 420328, indicating more robust and efficient RNA silencing in PI 420328. The distribution of vsRNA hotspots along the CYSDV genome was similar in both PI 420328 and Gold Star. However, the 3’ UTRs, CPm, and p26 were targeted at a higher frequency in PI 420328.

Expanding known viral diversity in plants: virome of 161 species alongside an ancient canal

BACKGROUND

Since viral metagenomic approach was applied to discover plant viruses for the first time in 2006, many plant viruses had been identified from cultivated and non-cultivated plants. These previous researches exposed that the viral communities (virome) of plants have still largely uncharacterized. Here, we investigated the virome in 161 species belonging to 38 plant orders found in a riverside ecosystem.

RESULTS

We identified 245 distinct plant-associated virus genomes (88 DNA and 157 RNA viruses) belonging to 27 known viral families, orders, or unclassified virus groups. Some viral genomes were sufficiently divergent to comprise new species, genera, families, or even orders. Some groups of viruses were detected that currently are only known to infect organisms other than plants. It indicates a wider host range for members of these clades than previously recognized theoretically. We cannot rule out that some viruses could be from plant contaminating organisms, although some methods were taken to get rid of them as much as possible. The same viral species could be found in different plants and co-infections were common.

CONCLUSIONS

Our data describe a complex viral community within a single plant ecosystem and expand our understanding of plant-associated viral diversity and their possible host ranges.

A comprehensive review and meta-analysis of recent advances in biotechnology for plant virus research and significant accomplishments in human health and the pharmaceutical industry

Secondary metabolites made by plants and used through their metabolic routes are today's most reliable and cost-effective way to make pharmaceuticals and improve health. The concept of genetic engineering is used for molecular pharming. As more people use plants as sources of nanotechnology systems, they are adding to this. These systems are made up of viruses-like particles (VLPs) and virus nanoparticles (VNPs). Due to their superior ability to be used as plant virus expression vectors, plant viruses are becoming more popular in pharmaceuticals. This has opened the door for them to be used in research, such as the production of medicinal peptides, antibodies, and other heterologous protein complexes. This is because biotechnological approaches have been linked with new bioinformatics tools. Because of the rise of high-throughput sequencing (HTS) and next-generation sequencing (NGS) techniques, it has become easier to use metagenomic studies to look for plant virus genomes that could be used in pharmaceutical research. A look at how bioinformatics can be used in pharmaceutical research is also covered in this article. It also talks about plant viruses and how new biotechnological tools and procedures have made progress in the field.

Exogenous Application of dsRNA for the Control of Viruses in Cucurbits

The recurrent emergence of viral diseases in intensive horticultural crops requires alternative control strategies. The topical application of double-stranded RNA (dsRNA) molecules homologous to pathogens has been proposed as a tool for virus control in plants. These dsRNAs induce the silencing mechanism, the RNA interference (RNAi), that degrades homologous dsRNAs. Cucumber green mottle mosaic virus (CGMMV) represents a serious threat to cucurbit crops. Since genetic resistance to the virus is not yet available in commercial varieties, we aimed to control this virus by RNAi. For this purpose, we obtained constructions both for expressing dsRNA in bacteria to treat cucumber plants by topical application and for agroinoculation in experiments done in the growth chamber. Besides, greenhouse tests were performed in spring and in summer when plants were challenged with the virus, and differences in several parameters were investigated, including the severity of symptoms, dry weight, total height, virus accumulation, and virus-derived small interfering RNAs (vsiRNAs). Spraying of plants with dsRNA reduced significatively CGMMV symptoms in the plants in growth chamber tests. Agroinfiltration experiments done under identical conditions were also effective in limiting the progress of CGMMV disease. In the greenhouse assay performed in spring, symptoms were significantly reduced in dsRNA-sprayed plants, and the development of the plants improved with respect to non-treated plants. Virus titers and vsiRNAs were clearly reduced in dsRNA-treated plants. The effect of protection of the dsRNA was less evident in the greenhouse assay carried out in the summer. Besides, we investigated the mobility of long (ds)RNA derived from spraying or agroinfiltrated dsRNA and found that it could be detected in local, close distal, and far distal points from the site of application. VsiRNAs were also detected in local and distal points and the differences in accumulation were compared. In parallel, we investigated the capacity of dsRNAs derived from genes of tomato leaf curl New Delhi virus (ToLCNDV), another economically important virus in cucurbits, to limit the disease in zucchini, both by agroinfiltration or by direct spraying, but found no protective effect. In view of the results, the topical application of dsRNAs is postulated as a promising strategy for CGMMV control in the cucumber.

Systematic Comparison of Nanopore and Illumina Sequencing for the Detection of Plant Viruses and Viroids Using Total RNA Sequencing Approach

High-throughput sequencing (HTS) has become an important tool for plant virus detection and discovery. Nanopore sequencing has been rapidly developing in the recent years and offers new possibilities for fast diagnostic applications of HTS. With this in mind, a study was completed, comparing the most established HTS platform (MiSeq benchtop sequencer-Illumina), with the MinION sequencer (Oxford Nanopore Technologies) for the detection of plant viruses and viroids. Method comparisons were performed on five selected samples, containing two viroids, which were sequenced using nanopore technology for the first time and 11 plant viruses with different genome organizations. For all samples, sequencing libraries for the MiSeq were prepared from ribosomal RNA-depleted total RNA (rRNA-depleted totRNA) and for MinION sequencing, direct RNA sequencing of totRNA was used. Moreover, for one of the samples, which contained five different plant viruses and a viroid, three additional variations of sample preparation for MinION sequencing were also used: direct RNA sequencing of rRNA-depleted totRNA, cDNA-PCR sequencing of totRNA, and cDNA-PCR sequencing of rRNA-depleted totRNA. Whilst direct RNA sequencing of total RNA was the quickest of the tested approaches, it was also the least sensitive: using this approach, we failed to detect only one virus that was present in a sample at an extremely low titer. All other MinION sequencing approaches showed improved performance with outcomes similar to Illumina sequencing, with cDNA-PCR sequencing of rRNA-depleted totRNA showing the best performance amongst tested nanopore MinION sequencing approaches. Moreover, when enough sequencing data were generated, high-quality consensus viral genome sequences could be reconstructed from MinION sequencing data, with high identity to the ones generated from Illumina data. The results of this study implicate that, when an appropriate sample and library preparation are selected, nanopore MinION sequencing could be used for the detection of plant viruses and viroids with similar performance as Illumina sequencing. Taken as a balance of practicality and performance, this suggests that MinION sequencing may be an ideal tool for fast and affordable virus diagnostics.

Hotspot siRNA Confers Plant Resistance against Viral Infection

Simple Summary

A hallmark of antiviral RNAi is the production of viral siRNA (vsiRNA). Profiling of vsiRNAs indicates that certain hotspot regions of viral genome or transcribed viral RNAs are more prone to RNAi-mediated cleavage. However, the biological relevance of hotspot vsiRNAs to the host innate defence remains to be elucidated. Here, we show that direct targeting a hotspot by synthetic vsiRNA confers plant resistance to virus infection. Hotspot and coldspot vsiRNAs, based on vsiRNA profile of the African cassava mosaic virus (ACMV), were synthesised. However, only the double-stranded hotspot vsiRNA protected plants from ACMV infection with undetectable levels of viral DNA replication and viral mRNA. We further demonstrated that the hotspot vsiRNA-mediated virus resistance had a threshold effect and required an active RDR6. These data show that hotspot vsiRNAs bear a functional significance on antiviral RNAi, suggesting that they may have the potential as exogenous protection agents for controlling destructive plant viral diseases.

Abstract

A hallmark of antiviral RNA interference (RNAi) is the production of viral small interfering RNA (vsiRNA). Profiling of vsiRNAs indicates that certain regions of viral RNA genome or transcribed viral RNA, dubbed vsiRNA hotspots, are more prone to RNAi-mediated cleavage for vsiRNA biogenesis. However, the biological relevance of hotspot vsiRNAs to the host innate defence against pathogens remains to be elucidated. Here, we show that direct targeting a hotspot by a synthetic vsiRNA confers host resistance to virus infection. Using Northern blotting and RNAseq, we obtained a profile of vsiRNAs of the African cassava mosaic virus (ACMV), a single-stranded DNA virus. Sense and anti-sense strands of small RNAs corresponding to a hotspot and a coldspot vsiRNA were synthesised. Co-inoculation of Nicotiana benthamiana with the double-stranded hotspot siRNA protected plants from ACMV infection, where viral DNA replication and accumulation of viral mRNA were undetectable. The sense or anti-sense strand of this hotspot vsiRNA, and the coldspot vsiRNA in both double-stranded and single-stranded formats possessed no activity in viral protection. We further demonstrated that the hotspot vsiRNA-mediated virus resistance had a threshold effect and required an active RDR6. These data show that hotspot vsiRNAs bear a functional significance on antiviral RNAi, suggesting that they may have the potential as an exogenous protection agent for controlling destructive viral diseases in plants.

Characterization of small interfering RNAs derived from pepper mild mottle virus in infected pepper plants by high-throughput sequencing

Pepper mild mottle virus (PMMoV) infects pepper plants and induces severe yield losses in China. However, the molecular interaction between PMMoV and pepper plants is largely unknown. RNA silencing is a eukaryotically conserved mechanism against viruses mediated by virus-derived small interfering RNAs (vsiRNAs) in plants. In this study, the profiles of vsiRNAs from PMMoV in infected pepper plants were obtained by high-throughput sequencing. The results showed that vsiRNAs were predominantly 21 and 22 nucleotides (nts) in length, and had a U bias at the 5'-terminal. The single-nucleotide resolution maps revealed that vsiRNAs were heterogeneously distributed throughout PMMoV genomic RNAs and hotspots of sense and antisense strands were mainly located in the RdRp and CP coding regions. The host transcripts targeted by vsiRNAs were predicted and they are mainly involved in physiological pathways related to stress response, cell regulation, and metabolism process. In addition, PMMoV infection induced significant up-regulation of CaAGO1a/1b/2, CaDCL2 and CaRDR1 gene transcripts in pepper plants, which are important components involved in antiviral RNA silencing pathway. Taken together, our results suggest the possible roles of vsiRNAs in PMMoV-pepper interactions.

In Silico Methods for the Identification of Viral-Derived Small Interfering RNAs (vsiRNAs) and Their Application in Plant Genomics

The current era of high-throughput sequencing (HTS) technology has expedited the detection and diagnosis of viruses and viroids in the living system including plants. HTS data has become vital to study the etiology of the infection caused by both known as well as novel viral elements in planta, and their impact on overall crop health and productivity. Viral-derived small interfering RNAs are generated as a result of defence response by the host via RNAi machinery. They are immensely exploited for performing exhaustive viral investigations in plants using bioinformatics as well as experimental approaches.This chapter briefly presents the basics of virus-derived small interfering RNAs (vsiRNAs ) biology in plants and their applications in plant genomics and highlights in silico strategies exploited for virus/viroid detection. It gives a systematic pipeline for vsiRNAs identification using currently available bioinformatics tools and databases. This will surely work as a quick beginner's recipe for the in silico revelation of plant vsiRNAs as well as virus/viroid diagnosis using high-throughput sequencing data.

Reprogramming of RNA silencing triggered by cucumber mosaic virus infection in Arabidopsis

BACKGROUND

RNA silencing has an important role mediating sequence-specific virus resistance in plants. The complex interaction of viruses with RNA silencing involves the loading of viral small interfering RNAs (vsiRNAs) into its host ARGONAUTE (AGO) proteins. As a side effect of their antiviral activity, vsiRNAs loading into AGO proteins can also mediate the silencing of endogenous genes. Here, we analyze at the genome-wide level both aspects of the interference of cucumber mosaic virus (CMV) with the RNA silencing machinery of Arabidopsis thaliana.

RESULTS

We observe CMV-derived vsiRNAs affect the levels of endogenous sRNA classes. Furthermore, we analyze the incorporation of vsiRNAs into AGO proteins with a described antiviral role and the viral suppressor of RNA silencing (VSR) 2b, by combining protein immunoprecipitation with sRNA high-throughput sequencing. Interestingly, vsiRNAs represent a substantial percentage of AGO-loaded sRNAs and displace other endogenous sRNAs. As a countermeasure, the VSR 2b loaded vsiRNAs and mRNA-derived siRNAs, which affect the expression of the genes they derive from. Additionally, we analyze how vsiRNAs incorporate into the endogenous RNA silencing pathways by exploring their target mRNAs using parallel analysis of RNA end (PARE) sequencing, which allow us to identify vsiRNA-targeted genes genome-wide.

CONCLUSIONS

This work exemplifies the complex relationship of RNA viruses with the endogenous RNA silencing machinery and the multiple aspects of virus resistance and virulence that this interaction induces.

Use of deep sequencing to profile small RNAs derived from tomato spotted wilt orthotospovirus and hippeastrum chlorotic ringspot orthotospovirus in infected Capsicum annuum

Virus-derived small RNAs are one of the key factors of RNA silencing in plant defence against viruses. We obtained virus-derived small interfering RNA profiles from Tomato spotted wilt orthotospovirus and Hippeastrum chlorotic ringspot orthotospovirus infected Capsicum annuum XX19 and XY11 by deep sequencing one day after inoculation. The vsiRNAs data were mapped to the TSWV and HCRV genomes, and the results showed that the vsiRNAs measured 19-24 nucleotides in length. Most of the vsiRNAs were mapped to the S segment of the viral genome. For XX19 and XY11 infected with HCRV, the distribution range of vsiRNAs in S RNA was 52.06-55.20%, while for XX19 and XY11 infected with TSWV, the distribution range of vsiRNAs in S RNA was 87.76-89.07%. The first base at the 5' end of the siRNA from TSWV and HCRV was primarily biased towards A, U, or C. Compared with mock-inoculated XX19 and XY11, the expression level of CaRDR1 was upregulated in TSWV- and HCRV-inoculated XX19 and XY11. CaAGO2 and CaAGO5 were upregulated in XY11 against HCRV infection, and CaRDR2 was downregulated in TSWV-infected XY11 and XX19. The profile of HCRV and TSWV vsiRNA verified in this study could be useful for selecting key vsiRNA such as those in disease-resistant varieties by artificially synthesizing amiRNA.

Small RNA and Transcriptome Sequencing of a Symptomatic Peony Plant Reveals Mixed Infections with Novel Viruses

To identify the viruses in tree peony plants associated with the symptoms of yellowing, leaf rolling, stunted growth, and decline, high-throughput sequencing of small RNA and mRNA was conducted from a single symptomatic plant. Bioinformatic analyses and reconstruction of viral genomes indicated mixed viral infections involving cycas necrotic stunt virus, apple stem grooving virus, lychnis mottle virus, grapevine line pattern virus, and three new viruses designated as peony yellowing-associated citrivirus (PYaCV, Citrivirus in Betaflexiviridae), peony betaflexivirus 1 (PeV1, unclassified in Betaflexiviridae), and peony leafroll-associated virus (PLRaV, Ampelovirus in Closteroviridae). PYaCV was 8,666 nucleotides (nt) in length, comprising three open reading frames (ORFs), and shared 63.8 to 75.9% nt sequence identity with citrus leaf blotch virus (CLBV) isolates. However, the ORF encoding the replication-associated protein (REP) shared 57 and 52% sequence identities at the nt and amino acid (aa) level, respectively, with other reported CLBV isolates, which were below the criterion for species classification within the family Betaflexiviridae. Recombination analysis identified putative recombination sites in PYaCV, which originated from CLBV. PeV1, only identified from the transcriptome data, was 8,124 nt in length, with five ORFs encoding the REP (ORF1), triple gene block (ORF2 to 4) and coat protein (CP, ORF5). Phylogenetic analysis and sequence comparison showed that PeV1 clustered with an unassigned member, the garlic yellow mosaic-associated virus within the Betaflexiviridae family, into a separate clade. Partial genome sequence analysis of PLRaV (12,545 nt) showed it contained seven ORFs encoding the partial polyprotein 1a, the RNA-dependent RNA polymerase (RdRp), two small hydrophobic proteins p11 and p6, HSP70h, p55, and a CP duplicate, which shared low aa sequence identity with Closteroviridae family members. Phylogenetic analysis based on the aa sequences of RdRp or HSP70h indicated that PLRaV clustered with grapevine leafroll-associated virus 1 (GLRaV-1) and GLRaV-13 in the Ampelovirus genus. Field investigation confirmed the wide distribution of these viruses, causing mixed infections of peony plants in Beijing.

Non-Coding RNAs in Response to Drought Stress

Drought stress causes changes in the morphological, physiological, biochemical and molecular characteristics of plants. The response to drought in different plants may vary from avoidance, tolerance and escape to recovery from stress. This response is genetically programmed and regulated in a very complex yet synchronized manner. The crucial genetic regulations mediated by non-coding RNAs (ncRNAs) have emerged as game-changers in modulating the plant responses to drought and other abiotic stresses. The ncRNAs interact with their targets to form potentially subtle regulatory networks that control multiple genes to determine the overall response of plants. Many long and small drought-responsive ncRNAs have been identified and characterized in different plant varieties. The miRNA-based research is better documented, while lncRNA and transposon-derived RNAs are relatively new, and their cellular role is beginning to be understood. In this review, we have compiled the information on the categorization of non-coding RNAs based on their biogenesis and function. We also discuss the available literature on the role of long and small non-coding RNAs in mitigating drought stress in plants.

Atypical molecular features of RNA silencing against the phloem-restricted polerovirus TuYV

In plants and some animal lineages, RNA silencing is an efficient and adaptable defense mechanism against viruses. To counter it, viruses encode suppressor proteins that interfere with RNA silencing. Phloem-restricted viruses are spreading at an alarming rate and cause substantial reduction of crop yield, but how they interact with their hosts at the molecular level is still insufficiently understood. Here, we investigate the antiviral response against phloem-restricted turnip yellows virus (TuYV) in the model plant Arabidopsis thaliana. Using a combination of genetics, deep sequencing, and mechanical vasculature enrichment, we show that the main axis of silencing active against TuYV involves 22-nt vsiRNA production by DCL2, and their preferential loading into AGO1. Moreover, we identify vascular secondary siRNA produced from plant transcripts and initiated by DCL2-processed AGO1-loaded vsiRNA. Unexpectedly, and despite the viral encoded VSR P0 previously shown to mediate degradation of AGO proteins, vascular AGO1 undergoes specific post-translational stabilization during TuYV infection. Collectively, our work uncovers the complexity of antiviral RNA silencing against phloem-restricted TuYV and prompts a re-assessment of the role of its suppressor of silencing P0 during genuine infection.

A virus-derived siRNA activates plant immunity by interfering with ROS scavenging

Virus-derived small interference RNAs (vsiRNAs) not only suppress virus infection in plants via induction of RNA silencing but also enhance virus infection by regulating host defensive gene expression. However, the underlying mechanisms that control vsiRNA-mediated host immunity or susceptibility remain largely unknown. In this study, we generated several transgenic wheat lines using four artificial microRNA expression vectors carrying vsiRNAs from Wheat yellow mosaic virus (WYMV) RNA1. Laboratory and field tests showed that two transgenic wheat lines expressing amiRNA1 were highly resistant to WYMV infection. Further analyses showed that vsiRNA1 could modulate the expression of a wheat thioredoxin-like gene (TaAAED1), which encodes a negative regulator of reactive oxygen species (ROS) production in the chloroplast. The function of TaAAED1 in ROS scavenging could be suppressed by vsiRNA1 in a dose-dependent manner. Furthermore, transgenic expression of amiRNA1 in wheat resulted in broad-spectrum disease resistance to Chinese wheat mosaic virus, Barley stripe mosaic virus, and Puccinia striiformis f. sp. tritici infection, suggesting that vsiRNA1 is involved in wheat immunity via ROS signaling. Collectively, these findings reveal a previously unidentified mechanism underlying the arms race between viruses and plants.

Quality Assessment and Validation of High-Throughput Sequencing for Grapevine Virus Diagnostics

Development of High-Throughput Sequencing (HTS), also known as next generation sequencing, revolutionized diagnostic research of plant viruses. HTS outperforms bioassays and molecular diagnostic assays that are used to screen domestic and quarantine grapevine materials in data throughput, cost, scalability, and detection of novel and highly variant virus species. However, before HTS-based assays can be routinely used for plant virus diagnostics, performance specifications need to be developed and assessed. In this study, we selected 18 virus-infected grapevines as a test panel for measuring performance characteristics of an HTS-based diagnostic assay. Total nucleic acid (TNA) was extracted from petioles and dormant canes of individual samples and constructed libraries were run on Illumina NextSeq 500 instrument using a 75-bp single-end read platform. Sensitivity was 98% measured over 264 distinct virus and viroid infections with a false discovery rate (FDR) of approximately 1 in 5 positives. The results also showed that combining a spring petiole test with a fall cane test increased sensitivity to 100% for this TNA HTS assay. To evaluate extraction methodology, these results were compared to parallel dsRNA extractions. In addition, in a more detailed dilution study, the TNA HTS assay described here consistently performed well down to a dilution of 5%. In that range, sensitivity was 98% with a corresponding FDR of approximately 1 in 5. Repeatability and reproducibility were assessed at 99% and 93%, respectively. The protocol, criteria, and performance levels described here may help to standardize HTS for quality assurance and accreditation purposes in plant quarantine or certification programs.

Identification of Plant Viruses Infecting Pear Using RNA Sequencing

Asian pear (Pyrus pyrifolia) is a widely cultivated and commercially important fruit crop, which is occasionally subject to severe economic losses due to latent viral infections. Thus, the aim of the present study was to examine and provide a comprehensive overview of virus populations infecting a major pear cultivar ('Singo') in Korea. From June 2017 to October 2019, leaf samples (n = 110) of pear trees from 35 orchards in five major pear-producing regions were collected and subjected to RNA sequencing. Most virus-associated contigs matched the sequences of known viruses, including apple stem grooving virus (ASGV) and apple stem pitting virus (ASPV). However, some contigs matched the sequences of apple green crinkle-associated virus and cucumber mosaic virus. In addition, three complete or nearly complete genomes were constructed based on transcriptome data and subjected to phylogenetic analyses. Based on the number of virus-associated reads, ASGV and ASPV were identified as the dominant viruses of 'Singo.' The present study describes the virome of a major pear cultivar in Korea, and looks into the diversity of viral communities in this cultivar. This study can provide valuable information on the complexity of genetic variability of viruses infecting pear trees.

Global Advances in Tomato Virome Research: Current Status and the Impact of High-Throughput Sequencing

Viruses cause a big fraction of economically important diseases in major crops, including tomato. In the past decade (2011–2020), many emerging or re-emerging tomato-infecting viruses were reported worldwide. In this period, 45 novel viral species were identified in tomato, 14 of which were discovered using high-throughput sequencing (HTS). In this review, we first discuss the role of HTS in these discoveries and its general impact on tomato virome research. We observed that the rate of tomato virus discovery is accelerating in the past few years due to the use of HTS. However, the extent of the post-discovery characterization of viruses is lagging behind and is greater for economically devastating viruses, such as the recently emerged tomato brown rugose fruit virus. Moreover, many known viruses still cause significant economic damages to tomato production. The review of databases and literature revealed at least 312 virus, satellite virus, or viroid species (in 22 families and 39 genera) associated with tomato, which is likely the highest number recorded for any plant. Among those, here, we summarize the current knowledge on the biology, global distribution, and epidemiology of the most important species. Increasing knowledge on tomato virome and employment of HTS to also study viromes of surrounding wild plants and environmental samples are bringing new insights into the understanding of epidemiology and ecology of tomato-infecting viruses and can, in the future, facilitate virus disease forecasting and prevention of virus disease outbreaks in tomato.

dsRNA Molecules From the Tobacco Mosaic Virus p126 Gene Counteract TMV-Induced Proteome Changes at an Early Stage of Infection

Exogenous application of double-stranded RNA (dsRNA) in the tobacco-Tobacco mosaic virus (TMV) pathosystem was shown previously to induce resistance against TMV providing an alternative approach to transgenesis. In the present study, we employed proteomics technology to elucidate the effect of TMV on tobacco as well as the effect of exogenous application of TMV p126 dsRNA molecules (dsRNAp126) at an early stage of the tobacco-TMV interaction. The proteome of tobacco leaf at 15 min post inoculation (mpi) in the presence or absence of dsRNAp126 molecules was studied. Thirty-six tobacco proteins were differentially accumulated in TMV-infected vs. healthy tobacco leaf tissue. The identified main differential TMV-responsive proteins were found to be involved in photosynthesis, energy metabolism, stress, and defense responses. Most of the virus-induced changes in the tobacco leaf proteome were not observed in the leaves treated with dsRNAp126 + TMV. The results indicated that the protein changes induced by TMV infection were counteracted by the exogenous application of dsRNAp126 molecules. Moreover, using small RNA sequencing, we showed that the exogenously applied dsRNAp126 was efficiently processed in tobacco as early as 15 min post application (mpa) to produce small interfering RNAs (siRNAs); the dicing pattern was not affected by the presence of TMV. The presence of dsRNAp126 reduced TMV p126 RNA abundance suggesting virus titer reduction via a sequence-specific mechanism, since a non-homologous dsRNA did not protect from TMV infection nor affect TMV accumulation.

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.

Non-perfectly matching small RNAs can induce stable and heritable epigenetic modifications and can be used as molecular markers to trace the origin and fate of silencing RNAs

Short non-coding RNA molecules (sRNAs) play a fundamental role in gene regulation and development in higher organisms. They act as molecular postcodes and guide AGO proteins to target nucleic acids. In plants, sRNA-targeted mRNAs are degraded, reducing gene expression. In contrast, sRNA-targeted DNA sequences undergo cytosine methylation referred to as RNA-directed DNA methylation (RdDM). Cytosine methylation can suppress transcription, thus sRNAs are potent regulators of gene expression. sRNA-mediated RdDM is involved in genome stability through transposon silencing, mobile signalling for epigenetic gene control and hybrid vigour. Since cytosine methylation can be passed on to subsequent generations, RdDM contributes to transgenerational inheritance of the epigenome. Using a novel approach, which can differentiate between primary (inducer) and secondary (amplified) sRNAs, we show that initiation of heritable RdDM does not require complete sequence complementarity between the sRNAs and their nuclear target sequences. sRNAs with up to four regularly interspaced mismatches are potent inducers of RdDM, however, the number and disruptive nature of nucleotide polymorphisms negatively correlate with their efficacy. Our findings contribute to understanding how sRNA can directly shape the epigenome and may be used in designing the next generation of RNA silencing constructs.

Characterization of Hibiscus Latent Fort Pierce Virus-Derived siRNAs in Infected Hibiscus rosa-sinensis in China

Although limited progress have been made about pathogen system of Hibiscus rosa-sinensis and Hibiscus latent Fort Pierce virus (HLFPV), interaction between plant host and pathogen remain largely unknown, which led to deficiency of effective measures to control disease of hibiscus plants caused by HLFPV. In this study, infection of HLFPV in Hibiscus rosa-sinensis was firstly confirmed for the first time by traditional electron microscopy, modern reverse transcription polymerase chain reaction and RNA-seq methods in China (HLFPV-Ch). Sequence properties analyzing suggested that the full-length sequences (6,465 nt) of HLFPV-Ch had a high sequence identity and a similar genomic structure with other tobamoviruses. It includes a 5'-terminal untranslated region (UTR), followed by four open reading frames encoding for a 128.5-kDa replicase, a 186.5-kDa polymerase, a 31-kDa movement protein, 17.6-kDa coat protein, and the last a 3'-terminal UTR. Furthermore, HLFPV-Ch-derived virus-derived siRNAs (vsiRNAs) ant its putative target genes, reported also for the first time, were identified and characterized from disease Hibiscus rosa-sinensis through sRNA-seq and Patmatch server to investigate the interaction in this pathogen systems. HLFPV-Ch-derived vsiRNAs demonstrated several general and specific characteristics. Gene Ontology classification revealed predicted target genes by vsiRNAs are involved in abroad range of cellular component, molecular function and biological processes. Taken together, for first time, our results certified the HLFPV infection in China and provide an insight into interaction between HLFPV and Hibiscus rosa-sinensis.

Biochemical characterization of the dicing activity of Dicer-like 2 in the model filamentous fungus Neurospora crassa

Dicing of double-stranded RNA (dsRNA) into small RNA is an essential process to trigger transcriptional and post-transcriptional gene silencing. Using cell-free extracts of the model filamentous fungus Neurospora crassa, we successfully detected the dicing activity of one of two N. crassa Dicers NcDCL2. The predominant 23-nucleotide (nt) cleavage product was always detected from 30-nt to 130-nt dsRNA substrates, and additional products of approximately 18 to 28 nt were occasionally produced. The enzymatic properties of NcDCL2 are different from those of insect and plant small interfering RNA (siRNA)-producing Dicers, Drosophila melanogaster Dicer-2 and Arabidopsis thaliana DCL3 and DCL4 (AtDCL3 and AtDCL4). Whereas AtDCL3 and AtDCL4 preferentially cleave short and long dsRNAs, respectively, NcDCL2 cleaved both short and long dsRNAs. These results suggest that N. crassa has a single siRNA-producing Dicer NcDCL2, which is a prototype of plant siRNA-producing Dicers with distinct functions in diverse RNA silencing pathways. The dicing assay reported here is convenient to detect and biochemically characterize the dicing activities of both plant and fungal Dicers, and is likely applicable to other organisms.

In vitro production of synthetic viral RNAs and their delivery into mammalian cells and the application of viral RNAs in the study of innate interferon responses

Mammalian cells express different types of RNA molecules that can be classified as protein coding RNAs (mRNA) and non-coding RNAs (ncRNAs) the latter of which have housekeeping and regulatory functions in cells. Cellular RNAs are not recognized by cellular pattern recognition receptors (PRRs) and innate immunity is not activated. RNA viruses encode and express RNA molecules that usually differ from cell-specific RNAs and they include for instance 5'capped and 5'mono- and triphosphorylated RNAs, small viral RNAs and viral RNA-protein complexes called vRNPs. These molecules are recognized by certain members of Toll-like receptor (TLR) and RIG-I-like receptor (RLR) families leading to activation of innate immune responses and the production of antiviral cytokines, such as type I and type III interferons (IFNs). Virus-specific ssRNA and dsRNA molecules that mimic the viral genomic RNAs or their replication intermediates can efficiently be produced by bacteriophage T7 DNA-dependent RNA polymerase and bacteriophage phi6 RNA-dependent RNA polymerase, respectively. These molecules can then be delivered into mammalian cells and the mechanisms of activation of innate immune responses can be studied. In addition, synthetic viral dsRNAs can be processed to small interfering RNAs (siRNAs) by a Dicer enzyme to produce a swarm of antiviral siRNAs. Here we describe the biology of RNAs, their in vitro production and delivery into mammalian cells as well as how these molecules can be used to inhibit virus replication and to study the mechanisms of activation of the innate immune system.

What does it take to be antiviral? An Argonaute-centered perspective on plant antiviral defense

RNA silencing is a major mechanism of constitutive antiviral defense in plants, mediated by a number of proteins, including the Dicer-like (DCL) and Argonaute (AGO) endoribonucleases. Both DCL and AGO protein families comprise multiple members. In particular, the AGO protein family has expanded considerably in different plant lineages, with different family members having specialized functions. Although the general mode of action of AGO proteins is well established, the properties that make different AGO proteins more or less efficient at targeting viruses are less well understood. In this report, we review methodologies used to study AGO antiviral activity and current knowledge about which AGO family members are involved in antiviral defense. In addition, we discuss what is known about the different properties of AGO proteins thought to be associated with this function.

Small RNA profiling analysis of two recombinant strains of potato virus Y in infected tobacco plants

Plant viral infections lead to accumulation of virus-derived small interfering RNAs (vsiRNAs) as a result of host defense mechanisms. High-throughput sequencing technology enables vsiRNA profiling analyses from virus infected plants, which provide important insights into virus-host interactions. Potato virus Y (PVY) is a detrimental plant pathogen that can infect a variety of solanaceous crops, e.g., potato, tobacco, tomato, and pepper. We analyzed and characterized vsiRNAs derived from Nicotiana tabacum cv. Samsun infected with two recombinant PVY strains, N-Wi and NTN. We observed that the average percentage of vsiRNAs derived from plants infected with N-Wi was higher than from plants infected with NTN, indicating that N-Wi invokes a stronger host response than NTN in tobacco. The size distribution pattern and polarity of vsiRNAs were similar between both virus strains with the 21 and 22 nucleotide (nt) vsiRNA classes as most predominant and the sense/antisense vsiRNAs ratio nearly equal in the 20-24 nt class. However, the percentage of sense vsiRNAs was significantly higher in the 25-26 nt long vsiRNAs. Distinct vsiRNA hotspots, identifying highly abundant reads of different unique vsiRNA sequences, were observed in both viral genomes. Previous studies found an A or U bias at the 5' terminal nucleotide position of 21 nt vsiRNAs; in contrast, our analysis revealed a C and U nucleotide bias. This study provides insights that will help further elucidate differential processing of vsiRNAs in plant antiviral defense.

High-Throughput Sequencing for Deciphering the Virome of Alfalfa (Medicago sativa L.)

Alfalfa (Medicago sativa L.), also known as lucerne, is a major forage crop worldwide. In the United States, it has recently become the third most valuable field crop, with an estimated value of over $9.3 billion. Alfalfa is naturally infected by many different pathogens, including viruses, obligate parasites that reproduce only inside living host cells. Traditionally, viral infections of alfalfa have been considered by breeders, growers, producers and researchers to be diseases of limited importance, although they are widespread in all major cultivation areas. However, over the past few years, due to the rapid development of high-throughput sequencing (HTS), viral metagenomics, bioinformatics tools for interpreting massive amounts of HTS data and the increasing accessibility of public data repositories for transcriptomic discoveries, several emerging viruses of alfalfa with the potential to cause serious yield losses have been described. They include alfalfa leaf curl virus (family Geminiviridae), alfalfa dwarf virus (family Rhabdoviridae), alfalfa enamovirus 1 (family Luteoviridae), alfalfa virus S (family Alphaflexiviridae) and others. These discoveries have called into question the assumed low economic impact of viral diseases in alfalfa and further suggested their possible contribution to the severity of complex infections involving multiple pathogens. In this review, we will focus on viruses of alfalfa recently described in different laboratories on the basis of the above research methodologies.

High-Throughput Sequencing Facilitates Discovery of New Plant Viruses in Poland

Viruses cause epidemics on all major crops of agronomic importance, and a timely and accurate identification is essential for control. High throughput sequencing (HTS) is a technology that allows the identification of all viruses without prior knowledge on the targeted pathogens. In this paper, we used HTS technique for the detection and identification of different viral species occurring in single and mixed infections in plants in Poland. We analysed various host plants representing different families. Within the 20 tested samples, we identified a total of 13 different virus species, including those whose presence has not been reported in Poland before: clover yellow mosaic virus (ClYMV) and melandrium yellow fleck virus (MYFV). Due to this new finding, the obtained sequences were compared with others retrieved from GenBank. In addition, cucurbit aphid-borne yellows virus (CABYV) was also detected, and due to the recent occurrence of this virus in Poland, a phylogenetic analysis of these new isolates was performed. The analysis revealed that CABYV population is highly diverse and the Polish isolates of CABYV belong to two different phylogenetic groups. Our results showed that HTS-based technology is a valuable diagnostic tool for the identification of different virus species originating from variable hosts, and can provide rapid information about the spectrum of plant viruses previously not detected in a region.

sRIS: A Small RNA Illustration System for Plant Next-Generation Sequencing Data Analysis

Small RNA (sRNA), such as microRNA (miRNA) and short interfering RNA, are well-known to control gene expression based on degradation of target mRNA in plants. A considerable amount of research has applied next-generation sequencing (NGS) to reveal the regulatory pathways of plant sRNAs. Consequently, numerous bioinformatics tools have been developed for the purpose of analyzing sRNA NGS data. However, most methods focus on the study of sRNA expression profiles or novel miRNAs predictions. The analysis of sRNA target genes is usually not integrated into their pipelines. As a result, there is still no means available for identifying the interaction mechanisms between host and virus or the synergistic effects between two viruses. For the present study, a comprehensive system, called the Small RNA Illustration System (sRIS), has been developed. This system contains two main components. The first is for sRNA overview analysis and can be used not only to identify miRNA but also to investigate virus-derived small interfering RNA. The second component is for sRNA target prediction, and it employs both bioinformatics calculations and degradome sequencing data to enhance the accuracy of target prediction. In addition, this system has been designed so that figures and tables for the outputs of each analysis can be easily retrieved and accessed, making it easier for users to quickly identify and quantify their results. sRIS is available at http://sris.itps.ncku.edu.tw/.

The Tomato spotted wilt virus (TSWV) Genome is Differentially Targeted in TSWV-Infected Tomato (Solanum lycopersicum) with or without Sw-5 Gene

Tospoviruses cause significant losses to a wide range of agronomic and horticultural crops worldwide. The type member, Tomato spotted wilt tospovirus (TSWV), causes systemic infection in susceptible tomato cultivars, whereas its infection is localized in cultivars carrying the Sw-5 resistance gene. The response to TSWV infection in tomato cultivars with or without Sw-5 was determined at the virus small RNA level in the locally infected leaf. Predicted reads were aligned to TSWV reference sequences. The TSWV genome was found to be differentially processed among each of the three-viral genomic RNAs—Large (L), Medium (M) and Small (S)—in the Sw-5(+) compared to Sw-5(−) genotypes. In the Sw-5(+) cultivar, the L RNA had the highest number of viral small-interfering RNAs (vsiRNAs), whereas in the Sw-5(−) cultivar the number was higher in the S RNA. Among the three-viral genomic RNAs, the distribution of hotspots showed a higher number of reads per million reads of vsiRNAs of 21 and 22 nt class at the 5′ and 3′ ends of the L and the S RNAs, with less coverage in the M RNA. In the Sw-5(−) cultivar, the nature of the 5′ nucleotide-end in the siRNAs varied significantly; reads with 5′-adenine-end were most abundant in the mock control, whereas cytosine and uracil were more abundant in the infected plants. No such differences were seen in case of the resistant genotype. Findings provided insights into the response of tomato cultivars to TSWV infection.

Rice Stripe Mosaic Virus-Encoded P4 Is a Weak Suppressor of Viral RNA Silencing and Is Required for Disease Symptom Development

Viral suppressors of RNA silencing (VSRs) are a cluster of viral proteins that have evolved to counteract eukaryotic antiviral RNA silencing pathways, thereby contributing to viral pathogenicity. In this study, we revealed that the matrix protein P4 encoded by rice stripe mosaic virus (RSMV), which is an emerging cytoplasmic rhabdovirus, is a weak RNA silencing suppressor. By conducting yeast two-hybrid, bimolecular fluorescence complementation, and subcellular colocalization assays, we proved that P4 interacts with the rice endogenous suppressor of gene silencing 3 (OsSGS3). We also determined that P4 overexpression has no effect on OsSGS3 transcription. However, P4 can promote the degradation of OsSGS3 via ubiquitination and autophagy. Additionally, a potato virus X-based expression system was used to confirm that P4 enhances the development of mosaic symptoms on Nicotiana benthamiana leaves by promoting hydrogen peroxide accumulation but not cell death. To verify whether P4 is a pathogenicity factor in host plants, we generated transgenic P4-overexpressing rice plants that exhibited disease-related developmental defects including decreased plant height and excessive tillering. Our data suggest that RSMV-encoded P4 serves as a weak VSR that inhibits antiviral RNA silencing by targeting OsSGS3.

Genome-wide analysis of small RNAs from Odontoglossum ringspot virus and Cymbidium mosaic virus synergistically infecting Phalaenopsis

Cymbidium mosaic virus (CymMV) and Odontoglossum ringspot virus (ORSV) are the two most prevalent viruses infecting orchids and causing economic losses worldwide. Mixed infection of CymMV and ORSV could induce intensified symptoms as early at 10 days post-inoculation in inoculated Phalaenopsis amabilis, where CymMV pathogenesis was unilaterally enhanced by ORSV. To reveal the antiviral RNA silencing activity in orchids, we characterized the viral small-interfering RNAs (vsiRNAs) from CymMV and ORSV singly or synergistically infecting P. amabilis. We also temporally classified the inoculated leaf-tip tissues and noninoculated adjacent tissues as late and early stages of infection, respectively. Regardless of early or late stage with single or double infection, CymMV and ORSV vsiRNAs were predominant in 21- and 22-nt sizes, with excess positive polarity and under-represented 5'-guanine. While CymMV vsiRNAs mainly derived from RNA-dependent RNA polymerase-coding regions, ORSV vsiRNAs encompassed the coat protein gene and 3'-untranslated region, with a specific hotspot residing in the 3'-terminal pseudoknot. With double infection, CymMV vsiRNAs increased more than 5-fold in number with increasing virus titres. Most vsiRNA features remained unchanged with double inoculation, but additional ORSV vsiRNA hotspot peaks were prominent. The potential vsiRNA-mediated regulation of the novel targets in double-infected tissues thereby provides a different view of CymMV and ORSV synergism. Hence, temporally profiled vsiRNAs from taxonomically distinct CymMV and ORSV illustrate active antiviral RNA silencing in their natural host, Phalaenopsis, during both early and late stages of infection. Our findings provide insights into offence-defence interactions among CymMV, ORSV and orchids.

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.

Highly efficacious antiviral protection of plants by small interfering RNAs identified in vitro

In response to a viral infection, the plant’s RNA silencing machinery processes viral RNAs into a huge number of small interfering RNAs (siRNAs). However, a very few of these siRNAs actually interfere with viral replication. A reliable approach to identify these immunologically effective siRNAs (esiRNAs) and to define the characteristics underlying their activity has not been available so far. Here, we develop a novel screening approach that enables a rapid functional identification of antiviral esiRNAs. Tests on the efficacy of such identified esiRNAs of a model virus achieved a virtual full protection of plants against a massive subsequent infection in transient applications. We find that the functionality of esiRNAs depends crucially on two properties: the binding affinity to Argonaute proteins and the ability to access the target RNA. The ability to rapidly identify functional esiRNAs could be of great benefit for all RNA silencing-based plant protection measures against viruses and other pathogens.

Genome-wide identification of the Dicer-like family in cotton and analysis of the DCL expression modulation in response to biotic stress in two contrasting commercial cultivars

BACKGROUND

Dicer-like proteins (DCLs) are essential players in RNA-silencing mechanisms, acting in gene regulation via miRNAs and in antiviral protection in plants and have also been associated to other biotic and abiotic stresses. To the best of our knowledge, despite being identified in some crops, cotton DCLs haven't been characterized until now. In this work, we characterized the DCLs of three cotton species and analyzed their expression profiles during biotic stress.

RESULTS

As main results, 11 DCLs in the allotetraploid cotton Gossypium hirsutum, 7 and 6 in the diploid G. arboreum and G. raimondii, were identified, respectively. Among some DCLs duplications observed in these genomes, the presence of an extra DCL3 in the three cotton species were detected, which haven't been found in others eudicots. All the DCL types identified by in silico analysis in the allotetraploid cotton genome were able to generate transcripts, as observed by gene expression analysis in distinct tissues. Based on the importance of DCLs for plant defense against virus, responses of cotton DCLs to virus infection and/or herbivore attack using two commercial cotton cultivars (cv.), one susceptible (FM966) and another resistant (DO) to polerovirus CLRDV infection, were analyzed. Both cvs. Responded differently to virus infection. At the inoculation site, the resistant cv. showed strong induction of DCL2a and b, while the susceptible cv. showed a down-regulation of these genes, wherever DCL4 expression was highly induced. A time course of DCL expression in aerial parts far from inoculation site along infection showed that DCL2b and DCL4 were repressed 24 h after infection in the susceptible cotton. As CLRDV is aphid-transmitted, herbivore attack was also checked. Opposite expression pattern of DCL2a and b and DCL4 was observed for R and S cottons, showing that aphid feeding alone may induce DCL modulation.

CONCLUSIONS

Almost all the DCLs of the allotetraploide G. hirsutum cotton were found in their relative diploids. Duplications of DCL2 and DCL3 were found in the three species. All four classes of DCL responded to aphid attack and virus infection in G. hirsutum. DCLs initial responses against the virus itself and/or herbivore attack may be contributing towards virus resistance.

Host factors against plant viruses

Plant virus genome replication and movement is dependent on host resources and factors. However, plants respond to virus infection through several mechanisms, such as autophagy, ubiquitination, mRNA decay and gene silencing, that target viral components. Viral factors work in synchrony with pro-viral host factors during the infection cycle and are targeted by antiviral responses. Accordingly, establishment of virus infection is genetically determined by the availability of the pro-viral factors necessary for genome replication and movement, and by the balance between plant defence and viral suppression of defence responses. Sequential requirement of pro-viral factors and the antagonistic activity of antiviral factors suggest a two-step model to explain plant-virus interactions. At each step of the infection process, host factors with antiviral activity have been identified. Here we review our current understanding of host factors with antiviral activity against plant viruses.

Effect of asymptomatic infection with southern tomato virus on tomato plants

Southern tomato virus (STV) is often found infecting healthy tomato plants (Solanum lycopersicum). In this study, we compared STV-free and STV-infected plants of cultivar M82 to determine the effect of STV infection on the host plant. STV-free plants exhibited a short and bushy phenotype, whereas STV-infected plants were taller. STV-infected plants produced more fruit than STV-free plants, and the germination rate of seeds from STV-infected plants was higher than that of seeds from STV-free plants. This phenotypic difference was also observed in progeny plants (siblings) derived from a single STV-infected plant in which the transmission rate of STV to progeny plants via the seeds was approximately 86%. These results suggest that the interaction between STV and host plants is mutualistic. Transcriptome analysis revealed that STV infection affects gene expression in the host plant and results in downregulation of genes involved in ethylene biosynthesis and signaling. STV-infected tomato plants might thus be artificially selected due to their superior traits as a crop.

Computational Workflow for Small RNA Profiling in Virus-Infected Plants

In this chapter we describe a series of computational pipelines for the in silico analysis of small RNAs (sRNA) produced in response to viral infections in plants. Our workflow is primarily focused on the analysis of sRNA populations derived from known or previously undescribed viruses infecting host plants. Furthermore, we provide an additional pipeline to examine host-specific endogenous sRNAs activated or specifically expressed during viral infections in plants. We present some key points for a successful and cost-efficient processing of next generation sequencing sRNA libraries, from purification of high quality RNA to guidance for library preparation and sequencing strategies. We report a series of free available tools and programs as well as in-house Perl scripts to perform customized sRNA-seq data mining. Previous bioinformatic background is not required, but experience with basic Unix commands is desirable.

Long DCL4-substrate dsRNAs efficiently induce RNA interference in plant cells

RNA interference (RNAi) is induced by the direct transfer of double-stranded RNAs (dsRNAs) into protoplasts prepared from Arabidopsis thaliana seedlings. In this protoplast RNAi system, we compared the efficacies of various-sized dsRNAs (between 21 and 139 nucleotides [nt]) for inducing RNAi and assessed the dsRNA-cleaving activities of Dicer-like 3 (DCL3) and 4 (DCL4). After the direct transfer of dsRNAs into protoplasts, cleaved RNA products of 21 nt were detected from long 130- or 500-nt dsRNAs by DCL4 but not from 37-nt dsRNAs. These results indicate that DCL4 preferentially cleaves long dsRNAs in protoplasts, consistent with our previous biochemical data regarding the substrate specificity of DCL4. Direct transfer of long dsRNAs of approximately 130 nt into protoplasts induces RNAi much more effectively (by approximately 60- to 400-fold) than direct transfer of short 37-nt dsRNAs. Although transfer of 21-nt dsRNAs into protoplasts induced RNAi without DCL4 activity, the induction of RNAi was less effective (by approximately 0.01-fold) compared with long dsRNAs. These results indicate that cleavage of long dsRNAs exceeding 100 nt by DCL4 into 21-nt dsRNAs is essential for efficient induction of RNAi in plant cells.

Molecular characterization of virus-derived small RNAs in Nicotiana benthamiana plants infected with tobacco curly shoot virus and its β satellite

Tobacco curly shoot virus (TbCSV) is a monopartite DNA virus of the genus Begomovirus, which causes leaf curl symptoms in tobacco and tomato. The β satellite of TbCSV (TbCSB induces more severe symptoms and enhanced virus accumulation when co-infects the host plants with TbCSV. Small interfering RNAs derived from virus(vsiRNAs) induce disease symptoms and promote virus invasion by target and guide the degradation of host transcripts The vsiRNAs derived from TbCSV and TbCSV + TbCSB remained to be explored to elucidate the molecular mechanism of symptoms development in plants. In the present work, two libraries of small RNA from TbCSV-infected and TbCSV + TbCSB-infected N. benthamiana plants were constructed and the vsiRNAs in both samples shared the same characteristics. The size of the vsiRNAs ranged from 18 to 30 nucleotides (nt), with most of them being 21 or 22 nt, which accounted for 29.11% and 23.22% in TbCSV plants and 29.39% and 21.82% in TbCSV + TbCSV plants, respectively. The vsiRNAs with A/U bias at the first site were abundant in both the TbCSV-treated and TbCSV + TbCSB-treated plants. It is discovered that the vsiRNAs continuously, but heterogeneously, distributed through bothe the TbCSV and TbCSB sequences. And the distribution profiles were similar in both the treatments such as mainly in the overlapping region of the AC2/AC3 coding sequences. The host transcripts targeted by vsiRNAs were predicted, and the targeted genes were found to be involved in varied biological processes. It is indicated that the presence of TbCSB does not significantly affect the production of vsiRNAs from TbCSV in plants, the distribution hotsopt of TbCSV vsiRNAs could be useful in designing effective targets for TbCSV resistance exploiting RNA interference.

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.