Identification of virus-derived siRNAs and their targets in RBSDV-infected rice by deep sequencing

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.

Contribution of Small RNA Pathway to Interactions of Rice with Pathogens and Insect Pests

Small RNAs (sRNAs) are mainly classified into microRNAs (miRNAs) and small interfering RNAs (siRNAs) according to their origin. miRNAs originate from single-stranded RNA precursors, whereas siRNAs originate from double-stranded RNA precursors that are synthesized by RNA-dependent RNA polymerases. Both of single-stranded and double-stranded RNA precursors are processed into sRNAs by Dicer-like proteins. Then, the sRNAs are loaded into ARGONAUTE proteins, forming RNA-induced silencing complexes (RISCs). The RISCs repress the expression of target genes with sequences complementary to the sRNAs through the cleavage of transcripts, the inhibition of translation or DNA methylation. Here, we summarize the recent progress of sRNA pathway in the interactions of rice with various parasitic organisms, including fungi, viruses, bacteria, as well as insects. Besides, we also discuss the hormone signal in sRNA pathway, and the emerging roles of circular RNAs and long non-coding RNAs in rice immunity. Obviously, small RNA pathway may act as a part of rice innate immunity to coordinate with growth and development.

Barley yellow dwarf virus-GAV-derived vsiRNAs are involved in the production of wheat leaf yellowing symptoms by targeting chlorophyll synthase

BACKGROUND

Wheat yellow dwarf virus disease is infected by barley yellow dwarf virus (BYDV), which causes leaf yellowing and dwarfing symptoms in wheat, thereby posing a serious threat to China's food production. The infection of plant viruses can produce large numbers of vsiRNAs, which can target host transcripts and cause symptom development. However, few studies have been conducted to explore the role played by vsiRNAs in the interaction between BYDV-GAV and host wheat plants.

METHODS

In this study, small RNA sequencing was conducted to profile vsiRNAs in BYDV-GAV-infected wheat plants. The putative targets of vsiRNAs were predicted by the bioinformatics software psRNATarget. RT-qPCR and VIGS were employed to identify the function of selected target transcripts. To confirm the interaction between vsiRNA and the target, 5' RACE was performed to analyze the specific cleavage sites.

RESULTS

From the sequencing data, we obtained a total of 11,384 detected vsiRNAs. The length distribution of these vsiRNAs was mostly 21 and 22 nt, and an A/U bias was observed at the 5' terminus. We also observed that the production region of vsiRNAs had no strand polarity. The vsiRNAs were predicted to target 23,719 wheat transcripts. GO and KEGG enrichment analysis demonstrated that these targets were mostly involved in cell components, catalytic activity and plant-pathogen interactions. The results of RT-qPCR analysis showed that most chloroplast-related genes were downregulated in BYDV-GAV-infected wheat plants. Silencing of a chlorophyll synthase gene caused leaf yellowing that was similar to the symptoms exhibited by BYDV-GAV-inoculated wheat plants. A vsiRNA from an overlapping region of BYDV-GAV MP and CP was observed to target chlorophyll synthase for gene silencing. Next, 5' RACE validated that vsiRNA8856 could cleave the chlorophyll synthase transcript in a sequence-specific manner.

CONCLUSIONS

This report is the first to demonstrate that BYDV-GAV-derived vsiRNAs can target wheat transcripts for symptom development, and the results of this study help to elucidate the molecular mechanisms underlying leaf yellowing after viral infection.

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.

Nitric Oxide Enhances Rice Resistance to Rice Black-Streaked Dwarf Virus Infection

BACKGROUND

Rice black-streaked dwarf virus (RBSDV) causes one of the most important rice virus diseases of plants in East Asia. However, molecular mechanism(s)controlling rice resistance to infection is largely unknown.

RESULTS

In this paper, we showed that RBSDV infection in rice significantly induced nitric oxide (NO) production. This finding was further validated through a genetic approach using a RBSDV susceptible (Nipponbare) and a RBSDV resistant (15HPO187) cultivar. The production of endogenous NO was muchhigher in the 15HPO187 plants, leading to a much lower RBSDV disease incidence. Pharmacological studies showed that the applications of NO-releasingcompounds (i.e., sodium nitroprusside [SNP] and nitrosoglutathione [GSNO]) to rice plants reduced RBSDV disease incidence. After RBSDV infection, the levels of OsICS1, OsPR1b and OsWRKY 45 transcripts were significantly up-regulated by NO in Nipponbare. The increased salicylic acid contents were also observed. After the SNP treatment, protein S-nitrosylation in rice plants was also increased, suggesting that the NO-triggered resistance to RBSDV infection was partially mediated at the post-translational level. Although Osnia2 mutant rice produced less endogenous NO after RBSDV inoculation and showed a higher RBSDV disease incidence, its RBSDV susceptibility could be reduced by SNP treatment.

CONCLUSIONS

Collectively, our genetic and molecular evidence revealed that endogenous NO was a vital signal responsible for rice resistance to RBSDV infection.

Plant virome reconstruction and antiviral RNAi characterization by deep sequencing of small RNAs from dried leaves

In plants, RNA interference (RNAi) generates small interfering (si)RNAs from entire genomes of viruses, satellites and viroids. Therefore, deep small (s)RNA sequencing is a universal approach for virome reconstruction and RNAi characterization. We tested this approach on dried barley leaves from field surveys. Illumina sequencing of sRNAs from 2 plant samples identified in both plants Hordeum vulgare endornavirus (HvEV) and barley yellow mosaic bymovirus (BaYMV) and, additionally in one plant, a novel strain of Japanese soil-borne wheat mosaic furovirus (JSBWMV). De novo and reference-based sRNA assembly yielded complete or near-complete genomic RNAs of these viruses. While plant sRNAs showed broad size distribution, viral sRNAs were predominantly 21 and 22 nucleotides long with 5′-terminal uridine or adenine, and were derived from both genomic strands. These bona fide siRNAs are presumably processed from double-stranded RNA precursors by Dicer-like (DCL) 4 and DCL2, respectively, and associated with Argonaute 1 and 2 proteins. For BaYMV (but not HvEV, or JSBWMV), 24-nucleotide sRNAs represented the third most abundant class, suggesting DCL3 contribution to anti-bymovirus defence. Thus, viral siRNAs are well preserved in dried leaf tissues and not contaminated by non-RNAi degradation products, enabling both complete virome reconstruction and inference of RNAi components mediating antiviral defense.

siRNAs Derived from Cymbidium Mosaic Virus and Odontoglossum Ringspot Virus Down-modulated the Expression Levels of Endogenous Genes in Phalaenopsis equestris

Interplay between Cymbidium mosaic virus (CymMV)/Odontoglossum ringspot virus (ORSV) and its host plant Phalaenopsis equestris remain largely unknown, which led to deficiency of effective measures to control disease of P. equestris caused by infecting viruses. In this study, for the first time, we characterized viral small interfering RNAs (vsiRNAs) profiles in P. equestris co-infected with CymMV and ORSV through small RNA sequencing technology. CymMV and ORSV small interfering RNAs (siRNAs) demonstrated several general and specific/new characteristics. vsiRNAs, with A/U bias at the first nucleotide, were predominantly 21-nt long and they were derived predominantly (90%) from viral positive-strand RNA. 21-nt siRNA duplexes with 0-nt overhangs were the most abundant 21-nt duplexes, followed by 2-nt overhangs and then 1-nt overhangs 21-nt duplexes in infected P. equestris. Continuous but heterogeneous distribution and secondary structures prediction implied that vsiRNAs originate predominantly by direct Dicer-like enzymes cleavage of imperfect duplexes in the most folded regions of the positive strand of both viruses RNA molecular. Furthermore, we totally predicted 54 target genes by vsiRNAs with psRNATarget server, including disease/stress response-related genes, RNA interference core components, cytoskeleton-related genes, photosynthesis or energy supply related genes. Gene Ontology classification showed that a majority of the predicted targets were related to cellular components and cellular processes and performed a certain function. All target genes were down-regulated with different degree by vsiRNAs as shown by real-time reverse transcription polymerase chain reaction. Taken together, CymMV and ORSV siRNAs played important roles in interplay with P. equestris by down modulating the expression levels of endogenous genes in host plant.

Rice black-streaked dwarf virus P10 acts as either a synergistic or antagonistic determinant during superinfection with related or unrelated virus

Rice black-streaked dwarf virus (RBSDV), a member of the genus Fijivirus, is a devastating pathogen of crop plants. RBSDV S10 encodes a capsid protein (P10) that is an important component of the double-layered particle. However, little information is available on the roles of RBSDV P10 in viral infection or in interactions with other viruses. Here, we demonstrate that the expression of P10 in plants alleviates the symptoms of both RBSDV and the closely related Southern rice black-streaked dwarf virus (SRBSDV), and reduces the disease incidence, but renders the plants more susceptible to the unrelated Rice stripe virus (RSV). Further experiments suggest that P10-mediated resistance to RBSDV and SRBSDV operates at the protein level, rather than the RNA level, and is not a result of post-transcriptional gene silencing. Transcriptomic data reveal that the expression of P10 in plants significantly suppresses the expression of rice defence-related genes, which may play important roles in resistance to RSV infection. After infection with RBSDV, plants are more resistant to subsequent challenge by SRBSDV, but more susceptible to RSV. Overall, these results indicate that P10 acts as an important effector in virus interactions.

Small RNA-Seq to Characterize Viruses Responsible of Lettuce Big Vein Disease in Spain

The emerging lettuce big-vein disease (LBVD) is causing losses in lettuce production ranging from 30 to 70% worldwide. Several studies have associated this disease with Mirafiori lettuce big-vein virus (MiLBVV) alone or in mixed infection with lettuce big-vein associated virus (LBVaV). We used Illumina small RNA sequencing (sRNA-seq) to identify viruses present in symptomatic lettuce plants from commercial fields in Southern Spain. Data analysis using the VirusDetect tool showed the consistent presence of MiLBVV and LBVaV in diseased plants. Populations of MiLBVV and LBVaV viral small RNAs (sRNAs) were characterized, showing features essentially similar to those of other viruses, with the peculiarity of an uneven asymmetric distribution of MiLBVV virus-derived small RNAs (vsRNAs) for the different polarities of genomic RNA4 vs. RNAs1 to 3. Sanger sequencing of coat protein genes was used to study MiLBVV and LBVaV phylogenetic relationships and population genetics. The Spanish MiLBVV population was composed of isolates from three well-differentiated lineages and reflected almost all of the diversity reported for the MiLBVV species, whereas the LBVaV population showed very little genetic differentiation at the regional scale but lineage differentiation at a global geographical scale. Universal primers were used to detect and quantify the accumulation of MiLBVV and LBVaV in field samples; both symptomatic and asymptomatic plants from affected fields carried equal viral loads, with LBVaV accumulating at higher levels than MiLBVV.

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.

Impact of Two Reoviruses and Their Coinfection on the Rice RNAi System and vsiRNA Production

Both Southern rice black-streaked dwarf virus (SRBSDV) and Rice ragged stunt virus (RRSV) belong to the family Reoviridae, and synergistic infection of these two viruses commonly occurs in the field. This study revealed that both SRBSDV and RRSV affect the RNA interference (RNAi) pathway and form different virus-derived interfering RNA (vsiRNA) profiles in rice. Co-infection of rice by SRBSDV and RRSV up-regulated the expression of rice DICER-like (DCL) proteins but down-regulated the expression of rice RNA-dependent RNA polymerases (RDRs), and the accumulation of vsiRNAs of either RBSDV or RRSV was decreased compared with that in singly infected plants. The majority of SRBSDV vsiRNAs were 21 nt or 22 nt in length, whether plants were singly infected with SRBSDV or co-infected with RRSV. On the other hand, the majority of RRSV vsiRNAs were 20 nt, 21 nt, or 22 nt in length, among which those 20 nt in length accounted for the largest proportion; co-infection with SRBSDV further increased the proportion of 20 nt vsiRNAs and decreased the proportion of 21 nt vsiRNAs. Co-infection had no effects on the strand favoritism and hot spots of the vsiRNAs, but changed the bias of the 5′ terminal nucleotide significantly. This study provides a reference for further study on the pathogenesis and synergistic mechanism of SRBSDV and RRSV.

Profile of siRNAs derived from green fluorescent protein (GFP)-tagged Papaya leaf distortion mosaic virus in infected papaya plants

We used green fluorescent protein (GFP)-tagged Papaya leaf distortion mosaic virus (PLDMV-GFP) to track PLDMV infection by fluorescence. The virus-derived small interfering RNAs (vsiRNAs) of PLDMV-GFP were characterized from papaya plants by next-generation sequencing. The foreign GFP gene inserted into the PLDMV genome was also processed as a viral gene into siRNAs by components involved in RNA silencing. The siRNAs derived from PLDMV-GFP accumulated preferentially as 21- and 22-nucleotide (nt) lengths, and most of the 5'-terminal ends were biased towards uridine (U) and adenosine (A). The single-nucleotide resolution map revealed that vsiRNAs were heterogeneously distributed throughout the PLDMV-GFP genome, and vsiRNAs derived from the sense strand were more abundant than those from the antisense strand. The hotspots were mainly distributed in the P1 and GFP coding region of the antisense strand. In addition, 979 papaya genes targeted by the most abundant 1000 PLDMV-GFP vsiRNAs were predicted and annotated using GO and KEGG classification. Results suggest that vsiRNAs play key roles in PLDMV-papaya interactions. These data on the characterization of PLDMV-GFP vsiRNAs will help to provide insight into the function of vsiRNAs and their host target regulation patterns.