- Open Access
Characteristics of siRNAs derived from Southern rice black-streaked dwarf virus in infected rice and their potential role in host gene regulation
Virology Journalvolume 14, Article number: 27 (2017)
Virus-derived siRNAs (vsiRNAs)-mediated RNA silencing plays important roles in interaction between plant viruses and their hosts. Southern rice black-streaked dwarf virus (SRBSDV) is a newly emerged devastating rice reovirus with ten dsRNA genomic segments. The characteristics of SRBSDV-derived siRNAs and their biological implications in SRBSDV-rice interaction remain unexplored.
VsiRNAs profiling from SRBSDV-infected rice samples was done via small RNA deep sequencing. The putative rice targets of abundantly expressed vsiRNAs were bioinformatically predicted and subjected to functional annotation. Differential expression analysis of rice targets and RNA silencing components between infected and healthy samples was done using RT-qPCR.
The vsiRNA was barely detectable at 14 days post infection (dpi) but abundantly present along with elevated expression level of the viral genome at 28 dpi. From the 28-dpi sample, 70,878 reads of 18 ~ 30-nt vsiRNAs were recognized (which mostly were 21-nt and 22-nt), covering 75 ~ 91% of the length of the ten genomic segments respectively. 86% of the vsiRNAs had a <50% GC content and 79% of them were 5’-uridylated or adenylated. The production of vsiRNAs had no strand polarity but varied among segment origins. Each segment had a few hotspot regions where vsiRNAs of high abundance were produced. 151 most abundant vsiRNAs were predicted to target 844 rice genes, including several types of host resistance or pathogenesis related genes encoding F-box/LRR proteins, receptor-like protein kinases, universal stress proteins, tobamovirus multiplication proteins, and RNA silencing components OsDCL2a and OsAGO17 respectively, some of which showed down regulation in infected plants in RT-qPCR. GO and KEGG classification showed that a majority of the predicted targets were related to cell parts and cellular processes and involved in carbohydrate metabolism, translation, and signal transduction. The silencing component genes OsDCL2a, OsDCL2b, OsDCL4, and OsAGO18 were down regulated, while OsAGO1d, OsAGO2, OsRDR1 and OsRDR6 were up regulated, significantly, upon SRBSDV infection.
SRBSDV can regulate the expression of rice RNA silencing pathway components and the virus might compromise host defense and influence host pathogenesis via siRNA pathways.
RNA silencing or RNA interfering (RNAi) is an evolutionarily conserved gene inactivation mechanism that universally exists in eukaryotic organisms [1, 2] and plays critical roles in developmental regulation, response to stresses, and host defense against transposons and viruses [3–5]. In plants, RNA silencing is triggered by double-stranded (ds) or highly structured single-stranded (ss) RNAs, which are processed by Dicer-like (DCL) ribonucleases into two classes of small RNA molecules - microRNAs (miRNAs) sized about 21 nucleotides (nt) and small interfering RNAs (siRNAs) of 21 ~ 24 nt . Invasion of plant viruses may initiate DCL-mediated biogenesis of primary virus-derived siRNAs (vsiRNAs) from double-stranded viral replicative intermediates (RI) or secondary-structured ssRNA genomic regions [7–9], and then the vsiRNAs are incorporated into Argonaute (AGO) proteins to form the RNA-induced silencing complex (RISC) that cleaves the viral RNAs or host mRNAs in a sequence-specific manner [9–13]. The RNA silencing mechanism can be considered as a type of pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) . RNA silencing effects can be amplified by host RNA-dependent RNA polymerases (RDR), which are acquired to convert aberrant RNAs (including the RISC-cleaved products) into dsRNAs, the materials for biogenesis of secondary vsiRNAs directed by DCLs [14, 15]. For survival, plant viruses have evolved various viral suppressors of RNA silencing (VSRs) to overcome RNA silencing-mediated defense [5, 16]. For the counter-counter defense, plants constantly evolve disease resistance (R) genes which encode R proteins to interact with VSR, triggering effector-triggered immunity (ETI) and leading to hypersensitive response (HR)-based resistance [5, 17].
Many DCLs, AGOs and RDRs contribute to RNA silencing-based plant antiviral defense. DCL2 and DCL4 process dsRNA molecules into 22- and 21-nt vsiRNAs in a hierarchical and redundant manner respectively, which play essential roles in resistance to positive-strand RNA viruses, while DCL3 mainly functions to produce 24-nt vsiRNAs involved in plant resistance to DNA viruses [5, 14, 18, 19]. It has been demonstrated that specific AGOs preferentially bind different small RNAs dependent on their 5’-terminal nucleotides [5, 11, 20–22]. RDR1, RDR2 and RDR6 are the major effectors involved in defense of Arabidopsis thaliana specifically against various viruses [14, 15, 17, 19]. In rice, DCL4 confers the defense against Rice yellow mottle virus (RYMV, genus Sobemovirus, unassigned family), which can be inhibited by the VSR, P1 protein of the virus . AGO1 and AGO18 participate in the defense against Rice stripe virus (RSV, genus Tenuivirus, unassigned family) and Rice dwarf phytoreovirus (RDV, genus Phytoreovirus, family Reoviridae), wherein AGO1 is the direct antiviral effector, while AGO18, induced upon virus infection, combines with miR168 to sequester the down regulation of AGO1 by miR168 and thus confers a broad-spectrum antivirus resistance . Also, OsRDR6 has been proven to play role in host resistance to RSV and RDV [25, 26]. Although vsiRNAs-endowed host defense against viruses is universal, they may benefit from the RNA silencing machinery in some scenarios. Recent research has discovered a role of vsiRNAs in targeting certain host genes, suggesting that vsiRNAs may modulate host response and symptomology through the RNA silencing machinery to facilitate virus survival and spread [12, 27–30].
Southern rice black-streaked dwarf virus (SRBSDV) is a newly emerged rice reovirus (genus Fijivirus, family Reoviridae) as well as an insect virus efficiently transmitted by the long-distance migratory pest, the white-backed planthopper (WBPH, Sogatella furcifera) in a persistent, circulative propagative manner [31, 32]. The virus has overspread to vast rice-growing areas in China and some eastern Asian countries in recent years and caused severe loss to rice production there [31, 32]. The virus has a genome consisting of ten dsRNA segments named S1 to S10 according to their molecular weights from large to small. SRBSDV S1 ~ S4, S6, S8 and S10 each encodes one protein, whereas S5, S7 and S9 are bicistronic. SRBSDV P1 ~ P4, P8 and P10 are structural proteins and the other seven proteins are nonstructural . Sequence analysis suggested that its P1 ~ P4 are the viral RNA-dependent RNA polymerase (RdRP), major core structural protein, outer-shell B-spike protein and capping enzyme, respectively [31, 33, 34]. SRBSDV P6 is a VSR, and it has complicated interactions with P5-1 and P9-1 to form viroplasm matrices in WBPH cells [35–38]. The P7-1 was found to form tubule structures in WBPH cells, which may serve for intercellular virus transportation and spread in the insect vector [39, 40]. The P9-1 is the major component of viral viroplasm that is important to viral replication in insect vectors . The P8 and P10 are putative minor core protein and major outer capsid protein, respectively [31, 33, 34]. Previous studies revealed that SRBSDV infection may regulate the expression of a bunch of rice miRNAs  and some of the RNAi-related genes in WBPH under different temporal conditions, suggesting interactions of the virus and host silencing pathways . It has been demonstrated that the core component Dicer-2 of siRNA pathway has decisive influences on SRBSDV accumulation and dissemination in insect and vector competence . Most recently, SRBSDV-derived siRNAs were reported to originate equally from both strand orientations of the viral genome by the Dicer enzyme in viruliferous WBPH, however, no further information about the vsiRNAs was presented by the authors . To better understand SRBSDV-derived siRNAs and gain insight into their biological implications, in this study we characterized the vsiRNAs present in infected rice by deep sequencing and explored their potential role in host gene regulation, and the influence of SRBSDV infection on rice RNA silencing pathway was analyzed as well.
Plant cultivation and virus inoculation
Seedlings of rice (Sativa japonica L. cv. Nipponbare) used in this study were obtained using the water planting method described by Yoshida et al. . Briefly, rice seeds after soaked in water for 24 h were germinated on wet towels at 30 °C and then sown in a 2-L beaker half-filled with culture solution to grow into three- to four-leaf stage seedlings in a growth chamber under the conditions of 25 ± 0.5 °C, RH 75 ± 5%, and 12 h light/12 h dark. Twenty seedlings with uniform growth were selected and each of them were inoculated with five WBPH nymphs in third- or fourth-instar, which previously were hatched and fed on the SRBSDV-infected rice plants maintained in our lab. In parallel, a group of seedlings were mock-inoculated with WBPH nymphs born and fed on healthy plants. The insects were allowed to feed on the seedlings for 48 h before being removed manually. Infection of each inoculated plant was confirmed by RT-PCR using the previous method  at 10 days post inoculation/infection (dpi).
RNA extraction, small RNA sequencing and analysis of siRNAs
Leaf samples of ten SRBSDV-infected plants were harvested at 14 and 28 dpi, respectively, and pooled for total RNA extraction using a TRIzol Reagent (TaKaRa, Dalian, China). The obtained RNA products were treated by DNase I (TaKaRa) and their concentration, integrity and purity were verified by an Agilent Bioanalyzer 2000 system (Agilent, CA, USA). Deep sequencing of small RNAs was performed by Beijing Genome Institute (BGI, Shenzhen, China) using their Illumina HiSeq platform. In brief, sRNAs of 18 ~ 30 nt were recovered from a 15% denaturing polyacrylamide gel electrophoresis of total RNA, ligated with RNA adaptors, and reversely transcribed into cDNAs. These cDNAs were amplified by PCR and subjected to Solexa/Illumina sequencing. From the generated data, adapter sequences were trimmed and clean reads of small RNAs between 18 ~ 30 nt in length were extracted. These small RNA sequences were mapped to the SRBSDV genomic sequences (GenBank accession numbers: NC_014708 ~ NC_014717) and only those identical to viral sequences (of sense or antisense strands) within 2 nt mismatches were recognized as vsiRNAs. According to their abundance levels, single-nucleotide resolution maps of the vsiRNAs were drawn using the OriginPro 9.0 software to display the distribution hotspots along the viral genome.
Prediction and annotation of target genes
Target genes prediction for the vsiRNAs of over 40 reads was perform by the online analysis tool psRNATarget (http://plantgrn.noble.org/psRNATarget/)  using default parameters and a rice transcript database, MSU Rice Genome Annotation (version 7). The predicted target genes were subject to BLAST alignment and Gene Ontology (GO) annotation using BLAST2GO software . KEGG classification of the predicted target genes was done using KEGG Mapper - Annotate Sequence by BlastKOALA available on the Kyoto Encyclopedia of Genes and Genomes website (http://www.kegg.jp/kegg/tool/annotate_sequence.html).
Reverse transcription-quantitative PCR
The total RNAs were extracted from the SRBSDV-infected (14 and 28 dpi) and mock-inoculated samples, treated by the DNase I, and used as templates for RT-qPCR, which was done using a One Step SYBR® PrimeScript™ RT-PCR Kit II (Perfect Real Time) (TaKaRa) in a Thermal Cycler® Dice Real Time System TP800 (TaKaRa) following the following the manufacturers’ instructions. The one-step amplification cycling conditions were 42 °C for 5 min, 95 °C for 10 s, and 40 cycles of 95 °C for 5 s and 60 °C for 30 s. After completion of qPCR cycling, the melting curves were generated at 95 °C to verify the specificity of amplification. For quantification of SRBSDV genomic segments (S1 and S10) and rice transcripts, a U6 small nuclear RNA gene and an actin gene of rice were used as internal control, respectively. All reactions were performed in triplicate. Data normalization based on Ct values was done by the 2−△△Ct method  and t-test analysis was conducted using SPSS19.0 software. The RT-qPCR primers used in this study were listed in Additional file 1: Table S1.
Profiling of SRBSDV-derived siRNAs in infected rice
A total of 10,479,093 and 12,131,156 clean reads were obtained from the small RNA libraries of SRBSDV-infected samples at 14 and 28 dpi, respectively. Among these, only 1,238 reads (0.01%, containing 1,016 unique sequences of 18 ~ 27 nt long) from the 14-dpi sample were mapped to the SRBSDV genome (36 sequences with 2 nt mismatch and the rest with perfect match), and the two most abundant sequences had only 23 and 11 reads, respectively (Additional file 2: Table S2). The GC content of these vsiRNAs varied from 16.7% to 66.7%, and 91% (927/1016) of the unique sequences had a GC content lower than 50%. In the 28-dpi sample, 70,878 reads (0.42%, containing 19,860 unique sequences of 18 ~ 30 nt long) were recognized as vsiRNAs with the most abundant one having 336 reads. Their GC contents were 20.8 ~ 73.7% and 86% (17141/19860) of the unique sequences had a GC content lower than 50% (Additional file 3: Table S3). There were 366 vsiRNAs present in both samples (Additional file 4: Table S4). More than one half of total vsiRNA reads (55.6%) and unique reads (56.3%) originated from sense strands of the viral genome at 14 dpi, whereas the amount (<51%) of sense strands-derived vsiRNAs was nearly equal to that of antisense strands-derived vsiRNAs (Fig. 1a). The molecules of 21- and 22-nt represented the majority of the recognized vsiRNAs in both samples - in total, 77.7% for 14 dpi and 93.3% for 28 dpi respectively (Fig. 1b). This suggests that OsDCL4 and OsDCL2 are the predominant enzymes for biogenesis of SRBSDV-derived siRNAs. The vsiRNA accumulation level varied among SRBSDV genomic segments from which they were generated (Fig. 1c). At 14 dpi, the greatest proportion (one sixth) of the vsiRNAs reads came from genomic S1, and one eighth from both S2 and S4, while those produced from S10 were the least abundant (5.3%). At 28 dpi, S4-derived vsiRNAs were the most abundant (representing 16% of the total reads), followed by those produced from S1 (11.9%), S5 (11.2%), S6 (10.9%) and S2 (10.2%), while S3 contributed the smallest part (7.4%).
The 5’ terminal nucleotides of the obtained vsiRNAs were examined and the result showed that two thirds of the vsiRNA reads generated at 14 dpi had an adenosine (A, 35.54%) or uridine (U, 31.83%) at their 5’ ends, and those having 5’ guanine (G) or cytosine (C) were less abundant; while in the 28 dpi sample, near half of the vsiRNAs had 5’ U, near one third were 5’-adenylated, and 5’ G was the least abundant (Fig. 1d). These results were similar to that from many RNA plant viruses  and suggested that SRBSDV-derived siRNAs mainly are recruited by the OsAGO1 and OsAGO2 to form the RISC, according to the AGO preferences for small RNAs revealed by the previous study .
Comparison of the abundance of SRBSDV S1 and S10 between the infected and uninfected samples by RT-qPCR analysis showed that the expression levels of the two segments at 28 dpi were 25 ~ 30 times (4.96 fold for S1 and 4.69 fold for S10) higher than that at 14 dpi respectively (Fig. 1e). This implied that the trace amount of vsiRNA accumulation at 14 dpi was due to limited titer of the virus in the early stage of infection. With regard to this, our following analyses are based on the vsiRNA data obtained from the 28-dpi sample.
Distribution of vsiRNAs in the SRBSDV genome
To gain further insight of the origin of the vsiRNAs, we calculated the frequency of the viral genome’s each single nt position for being covered by the vsiRNA reads (Additional file 5: Table S5), and drew single nucleotide resolution maps for each SRBSDV genomic segment according to the “nt position-frequency of being covered” data (Fig. 2). The results indicated that the vsiRNAs were not always continuously distributed through the genomic segments. Instead, there were 9.3% (S3) to 25.4% (S1) of nt positions not being covered by the vsiRNA sequences from both sense and antisense orientations (Additional file 5: Table S5). As shown in Fig. 2, distribution of vsiRNAs along the viral genome was highly heterogeneous. For each of S1 ~ S10, we defined the region as a distribution hotspot when it contained at least 18 consecutive nucleotides and all their single-nucleotide read counts reached two thirds of the highest value in the segment. It is showed that the S1(−), S7 (−), S8 (−) and S10 (−), S4 (+) and S9 (+) strands each contained one hotspot, both S2 (+) and S5 (+) strands contained two hotspots, and both sense and antisense strands of S3 and S6 had a hotspot respectively (Fig. 2). Except the one from the S6 (−) strand that was located in the 3’-UTR region, all the hotspots resided in coding regions of the viral genes. However, of these 14 hotspots, four peaks (from the S2 (+), S4 (+) and S3 (−) strands respectively) were not corresponding to the most abundant vsiRNA sequences identified from their respective segments; instead, some abundant vsiRNAs from several segments formed lower peaks in the single-nucleotide resolution maps (Fig. 2). These four exceptions came from the fact that the maps were created based on the frequencies of covering by vsiRNAs for single nt positions, not for complete siRNA-generating spots (i.e. 21 ~ 24 nt long regions). Four hotspots were found to generate a 21-nt and a 22-nt vsiRNA simultaneously: vsiRS1_1 and vsiRS1_2, vsiRS5_1 and vsiRS5_2, vsiRS5_5 and vsiRS5_6, and vsiRS6_1 and vsiRS6_2 from the same site of the S1 (−), S5 (+) and S6 (+) strands, respectively. The abundance levels of the 22-nt vsiRNAs were two times (S1) or slightly higher (S5 and S6) than that of the 21-nt ones (Table 1). This result indicated that the same sites on the viral genome may be recognized and cleaved by DCL4 and DCL2, and possibly, by DCL2 in preference. The similar phenomenon can also be found in other plant viruses [29, 30].
Prediction of putative rice target genes of vsiRNAs
In this study, we used the small RNA target analysis server psRNATarget  to explore the rice genes putatively targeted by SRBSDV-derived siRNAs. Because of the large number of vsiRNA species identified, only the 168 vsiRNAs with over 40 reads were selected for target prediction, including seven sequences of 40 ~ 58 reads from S2, the only segment not producing vsiRNAs of over 100 reads (Table 1). In total, 975 putative pairs of vsiRNA-target were predicted, including 151 unique vsiRNAs and 844 individual rice genes (Additional file 6: Table S6). Twenty-five out of the 28 most abundant siRNAs from nine segments (excluding S2) were predicted to each target 1 ~ 32 host genes, respectively, while three of them (vsiRS5_5, vsiRS5_6 and vsiRS9_1) had no rice targets (Table 1). Notably, a vsiRNA coming from 3’-UTR region of S8 putatively had 32 targets, 26 of which were retrotransposons (Fig. 2; marked by an asterisk). Besides, the seven most abundant vsiRNAs from S2 were predicted to have 1 ~ 6 rice targets, respectively. Annotation of the 844 genes predicted was done referring to the rice transcript database MSU Rice Genome Annotation (version 7). The result showed that a number of annotated targets were related to disease/stress response, including 33 R genes (13 of which were annotated as NBS-LRR or NBS protein genes), 17 F-box/LRR or F-box domain-containing protein genes, 15 receptor-like protein kinase (RLK) genes, five universal stress protein (USP) genes, and six tobamovirus multiplication protein genes (TOM). A majority of plant R genes encode NBS-LRR proteins that play roles in pathogen sensing/detection and host defense against a variety of pathogens including bacteria, fungi, nematodes and viruses [51–53], and some R genes may interact with VSRs to induce HR-mediated resistance [5, 54, 55]. Many plant F-box proteins as the subunits of the SCF complex, a type of ubiquitin E3 ligase, are involved in multiple biological and developmental processes and stress responses [56, 57] and can be regulated by the RNA silencing pathways in plants . Also, it has been reported that F-box proteins can reduce the expression of AGO1 at protein level expression in Arabidopsis, and some plant viruses encode an F-box domain-containing VSR to degrade this antiviral component through the autophagy pathway [59, 60]. RLKs are a large family of proteins that function in various signal transduction pathways and participate in hormonal response pathways, cell differentiation, plant growth and development, self-incompatibility, and pathogen recognition and resistance [61–63]. USPs can increase drought tolerance of plants [64, 65]. TOM proteins are a type of plant transmembrane proteins that are required by tobamoviral genome replication and may suppress the host’s RNA silencing process when they are overexpressed . Our analysis suggested that the several kinds of stress resistance genes above-mentioned were likely down regulated by vsiRNAs, making the SRBSDV-infected host become more susceptible to other biotic or abiotic stimuli. Besides, two rice genes LOC_Os03g38740.1 and LOC_Os02g07310.1, putatively targeted by the vsiRS8_16 and vsiRS2_1, were annotated as OsDCL2a and OsAGO17, respectively (Additional file 6: Table S6). OsAGO17 is male gametophyte-specific and likely involved in rice pollen development regulated by miRNAs . Moreover, seven chloroplast envelope membrane protein or chloroplast precursor-related genes were annotated, implying that photosynthesis pathways might be influenced by vsiRNAs in infected rice.
In BLAST2GO analysis, a total of 409 predicted target genes of vsiRNAs were annotated and classified into three GO categories (biological process, cellular component, and molecular function with 366, 378, 294 individual genes respectively) and 34 subcategories (Fig. 3). The three most highly represented GO terms were “cell”, “cell part” and “organelle” under the cellular component category, comprising 75 ~ 92% of the annotated genes; five terms under the biological process section, “cellular process”, “metabolic process”, “single-organism process”, “response to stimuli” and “developmental process”, represented 40 ~ 60% of the classified targets; and the two most abundant GO terms under the molecular function category, “catalytic activity” and “binding”, included 48% and 37% of the annotated genes respectively (Fig. 3).
Furthermore, we conducted KEGG classification for the predicted target genes of vsiRNAs. In total, 129 genes were classified into six classes – metabolism, genetic information processing, environmental information processing, cellular processes, organismal systems, and human diseases – and then further assigned to 27 KEGG pathways (Fig. 4). It showed that the five most involved pathways were “carbohydrate metabolism”, “translation”, “signal transduction”, “folding, sorting and degradation” and “infectious diseases”, which represented 74% of the annotated genes.
To check whether the predicted target genes were down regulated by vsiRNAs, we carried out RT-qPCR confirmation for 18 selected rice genes, which putatively were targeted by the 12 vsiRNAs of high abundance (>100 reads, except vsiRS2_1) coming from every genomic segments (Table 2). The result showed that the all the 18 genes were down regulated (fold change = −0.96 ~ −11.88) and 11 of them had significantly decreased expression levels, in the SRBSDV-infected sample (28 dpi) compared to the uninfected control (Fig. 5). However, all the 18 genes had no differential expression between the infected and virus-free samples at 14 dpi. These results suggested that these highly expressed vsiRNAs can most likely target some certain host genes.
Differential expression analysis of several RNA silencing components
To investigate the influence of SRBSDV infection on the RNA silencing pathways in rice, we compared the expression levels of some representative silencing components between the SRBSDV-infected and mock-inoculated samples. The results revealed that OsDCL2a, OsDCL2b, OsDCL4, and OsAGO18 were significantly down regulated (among which OsDCL2a was predicted as a target of vsiRS8_16 in this study, as aforementioned), while OsAGO1d, OsAGO2, OsRDR1 and OsRDR6 significantly up regulated, in the infected sample compared with the no infection control (Fig. 6). Notably, the abundance of OsAGO2 became 5.35 fold higher in the infected sample, suggesting its key role implicated in the host’s response to and/or defense against the virus infection.
With the advantage of next-generation sequencing technology, vsiRNA populations have been characterized from a variety of origins, including animal and fungal viruses, and many plant viruses of ssRNA or DNA genomes [8, 9, 68–71], while the characteristics of dsRNA plant virus-derived siRNAs remain relatively less studied. Li et al. profiled the vsiRNAs in Laodelphax striatellus infected with RBSDV, the species most closely related to SRBSDV, and reported that they were predominantly 21 and 22 nt long, equally derived from both strand orientations, and frequently produced from hotspots located in 5’- or 3’-terminal regions of viral genomic segments . A most recent study focusing on WBPH miRNA profile also reported the biogenesis of SRBSDV siRNAs without strand favoritism in the viruliferous insect, but the authors did not provide further information about the vsiRNAs . In this study, we expanded the knowledge of dsRNA plant virus-derived siRNAs by presenting the first detailed profile of SRBSDV siRNAs in rice at early and later periods of infection (i.e. 14 and 28 dpi). Interestingly, we found that the regions covering up to 20 ~ 25% of the length of SRBSDV S1 ~ S3 did not produce vsiRNAs. Since a same given virus in different hosts may generate diverse vsiRNA profiles in terms of polarity in strand origins, population constitution, distribution hotspot and pattern in the viral genome, etc. , further study will reveal the diversity of SRBSDV siRNA profiles from different hosts like rice, maize and WBPH. Similar to RBSDV- or RSV-derived siRNAs, SRBSDV siRNA population underwent temporal alterations in constitution and abundance (Fig. 1E), and it was likely that increased expression of viral RNA provided more substrates for vsiRNA production [72, 73]. This indicated that the virus triggered the host’s RNA silencing machinery and overcame it to establish successful infection.
vsiRNAs are known to be the core component of RNA silencing that direct host AGOs to cleave viral mRNAs at the complementary sites . In general, vsiRNAs are abundantly produced and heterogeneously distributed throughout the viral genome with a number of biogenesis hotspots [8, 9, 15, 29, 30, 70, 73]. Researchers have found that the efficiencies of vsiRNAs in targeting viral RNA are not correlated with whether they come from hotspot regions , suggesting that hotspot vsiRNAs have some biological functions other than targeting viral genes. Evidences have demonstrated that vsiRNAs can regulate host genes by exploiting the RNA silencing mechanism, and therefore enhance their fitness and influence host pathogenesis. The siRNAs derived from the satellite RNA of cucumber mosaic virus-Y (genus Cucumovirus, family Bromoviridae) was found to down regulate the key gene (CHLI) involved in chlorophyll biosynthesis, accounting for the yellowing symptom in Nicotiana benthamiana . Peach latent mosaic viroid produces two siRNAs that target the chloroplastic heat shock protein 90 gene in peach to induce host albinism and potentially create an improved host environment for infection . A siRNA derived from potato spindle tuber viroid may silence two callose synthase genes, greatly affecting the disease severity and viroid accumulation . Also, the vsiRNAs present in the maize infected with sugarcane mosaic virus (SCMV, genus Potyvirus, family Potyviridae) and/or maize chlorotic mottle virus (MCMV, genus Machlomovirus, family Tombusviridae) were predicted to regulated multiple host genes [29, 30]. In our study, we identified a bunch of SRBSDV siRNAs potentially targeting plentiful rice genes involved in various processes and pathways (Figs. 3 and 4), especially the genes related to disease resistance/stress response, pathogenesis, and developmental processes of the host. We then confirmed the down regulation of a part of predicted rice target genes in infected samples via RT-qPCR (Table 2, Fig. 5). Our results suggested that SRBSDV infection might generally reduce the host’s tolerance to multiple kinds of other pathogens or abiotic stresses and thus aggravate the severity and epidemic of diseases in rice growing areas. In natural conditions, coinfection of SRBSDV and another rice reovirus, rice ragged stunt virus (genus Oryzavirus, family Reoviridae), has been found in southern China, and their synergism enhances symptom development, virus accumulation, and vector acquisition and transmission of both the viruses [74, 75]. The synergistic co-infection of SCMV and MCMV could increase accumulation of both viruses-derived siRNAs in maize . Whether similar phenomenon occurs in SRBSDV and RRSV co-infected rice and the role of vsiRNAs implicated in their synergism deserve to be uncovered. We also found that vsiRNAs might down regulate several chloroplast-related genes, the pollen development-related gene AGO17, and many carbohydrate metabolism-related genes (Table 2, Figs. 4 and 5), suggesting the potential roles of vsiRNAs in development of the leaf symptoms, barren grains and dwarfism caused by SRBSDV [31, 32].
To address the interactions of SRBSDV and rice RNA silencing pathway, we examined the differential expression of several representative silencing component genes crucial to plant antiviral defense. Down regulation of OsDCL2 and OsDCL4 in SRBSDV-infected rice plants may influence the production of 21- and 22-nt vsiRNAs and thus weaken host defense, allowing viruses to establish successful infection. OsDCL4 is known as the SHOOT ORGANIZATION 1 gene (SHO1), and studies have discovered that mutations of OsDCL4 and OsRDR6 can induce abnormal expression of Auxin Response Factors (ARFs) and result in many growth abnormities including defects in spikes and lateral organs [76–78]. The down regulation of OsDCL4 in SRBSDV-infected rice might partially account for the symptoms of small spikes and poorly-developed roots [31, 32]. Many VSRs can block biogenesis of vsiRNAs by interacting with plant DCLs ; for example, OsDCL4 is suppressed by RYMV P1 . It is yet to validate whether SRBSDV P6 interacts with the two OsDCLs. Plant AGO1 and AGO2 participate in host defense against many species of RNA viruses and can be suppressed by VSRs [79, 80]. In rice resistance against RDV and RSV, OsAGO18 protects OsAGO1 from being silenced by miR168, conferring a broad-spectrum antivirus resistance . In our study, SRBSDV infection induced up regulation of OsAGO1 and OsAGO2 but down regulation of OsAGO18, and in particular, OsAGO2 was drastically up regulated (Fig. 6), suggesting facilitated host antiviral defense. The dramatically activated AGO2 might not only act in vsiRNA binding, but also have additional biological functions, such as participating in miRNA-mediated host gene regulation responding to SRBSDV infection  or counterbalancing the suppression of AGO2 by AGO1 via miRNA pathways . Whether the down regulated OsAGO17 in SRBSDV infected plant affected RNA silencing-directed antiviral defense is unclear. OsRDR6 plays role in defense against RSV and RDV [25, 26], and it is crucial to biogenesis of 21-nt and 24-nt small RNAs and required for normal development of spikelets [77, 78]. In our study, OsRDR1 and OsRDR6 were up regulated in infected rice, implying their roles in host response to SRBSDV infection. In general, our results suggested the complicated interactions between the virus and rice RNA silencing pathway, which likely are modulated by the intertwined effects of multiple silencing components.
In this study, we profiled the vsiRNAs in SRBSDV-infected rice plants, and explored their potential roles in targeting a series of rice genes related to host defense, pathogenesis and symptomology. We also found that SRBSDV infection could affect the expression of several rice RNA silencing pathway components. The results obtained in this study provide an understanding of SRBSDV-derived siRNAs and an insight into the implications of RNA silencing pathway in reovirus-rice interaction.
Days post inoculation
RNA-dependent RNA polymerase
RNA-induced silencing complex
Reverse transcription-quantitative PCR
Southern rice black-streaked dwarf virus
Viral suppressor of RNA silencing
Shabalina SA, Koonin EV. Origins and evolution of eukaryotic RNA interference. Trends Ecol Evol. 2008;23:578–87. doi:10.1016/j.tree.2008.06.005.
Weiberg A, Jin H. Small RNAs – the secret agents in the plant-pathogen interactions. Curr Opin Plant Biol. 2015;26:87–94. doi:10.1016/j.pbi.2015.05.033.
Castel SE, Martienssen RA. RNA interference (RNAi) in the Nucleus: roles for small RNAs in transcription, epigenetics and beyond. Nat Rev Genet. 2013;14:100–12. doi:10.1038/nrg3355.
Martinez de Alba AE, Elvira-Matelot E, Vaucheret H. Gene silencing in plants: A diversity of pathways. Biochim Biophys Acta. 2013;1829:1300–8. doi:10.1016/j.bbagrm.2013.10.005.
Pumplin N, Voinnet O. RNA silencing suppression by plant pathogens: defence, counter-defence and counter-counter-defence. Nat Rev Microbiol. 2013;11:745–60. doi:10.1038/nrmicro3120.
Brodersen P, Voinnet O. The diversity of RNA silencing pathways in plants. Trends Genet. 2006;22:268–80. doi:10.1016/j.tig.2006.03.003.
Molnar A, Csorba T, Lakatos L, Varallyay E, Lacomme C, Burgyán J. Plant Virus-Derived Small Interfering RNAs Originate Predominantly from Highly Structured Single-Stranded Viral RNAs. J Virol. 2005;79:7812–8. doi:10.1128/JVI.79.12.7812-7818.2005.
Donaire L, Wang Y, Gonzalez-Ibeas D, Mayer KF, Aranda MA, Llave C. Deep-sequencing of plant viral small RNAs reveals effective and widespread targeting of viral genomes. Virology. 2009;392:203–14. doi:10.1016/j.virol.2009.07.005.
Szittya G, Moxon S, Pantaleo V, Toth G, Rusholme Pilcher RL, Moulton V, et al. Structural and Functional Analysis of Viral siRNAs. PLoS Pathog. 2010;6:e1000838. doi:10.1371/journal.ppat.1000838.
Baumberger N, Baulcombe DC. Arabidopsis ARGONAUTE1 is an RNA Slicer that selectively recruits microRNAs and short interfering RNAs. Proc Natl Acad Sci U S A. 2005;102:11928–33. doi:10.1073/pnas.0505461102.
Mi S, Cai T, Hu Y, Chen Y, Hodges E, Ni F, et al. Sorting of Small RNAs into Arabidopsis Argonaute Complexes Is Directed by the 5’ Terminal Nucleotide. Cell. 2008;133:116–27. doi:10.1016/j.cell.2008.02.034.
Shimura H, Pantaleo V, Ishihara T, Myojo N, Inaba J, Sueda K, et al. A Viral Satellite RNA Induces Yellow Symptoms on Tobacco by Targeting a Gene Involved in Chlorophyll Biosynthesis using the RNA Silencing Machinery. PLoS Pathog. 2011;7:e1002021. doi:10.1371/journal.ppat.1002021.
Smith NA, Eamens AL, Wang MB. Viral Small Interfering RNAs Target Host Genes to Mediate Disease Symptoms in Plants. PLoS Pathog. 2011;7:e1002022. doi:10.1371/journal.ppat.1002022.
Garcia-Ruiz H, Takeda A, Chapman EJ, Sullivan CM, Fahlgren N, Brempelis KJ, et al. Arabidopsis RNA-Dependent RNA Polymerases and Dicer-Like Proteins in Antiviral Defense and Small Interfering RNA Biogenesis during Turnip mosaic virus Infection. Plant Cell. 2010;22:481–96. doi:10.1105/tpc.15.00204.
Wang XB, Wu Q, Ito T, Cillo F, Li WX, Chen X, et al. RNAi-mediated viral immunity requires amplification of virus-derived siRNAs in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2010;107:484–9. doi:10.1073/pnas.0904086107.
Incarbone M, Dunoyer P. RNA silencing and its suppression: novel insights from in planta analyses. Trends Plant Sci. 2013;18:382–92. doi:10.1016/j.tplants.2013.04.001.
Garcia-Ruiz H, Carbonell A, Hoyer JS, Fahlgren N, Gilbert KB, Takeda A, et al. Roles and Programming of Arabidopsis ARGONAUTE Proteins during Turnip mosaic virus Infection. PLoS Pathog. 2015;11:e1004755. doi:10.1371/journal.ppat.1004755.
Ding SW. RNA-based antiviral immunity. Nat Rev Immunol. 2010;10:632–44. doi:10.1038/nri2824.
Qu F, Ye X, Morris TJ. Arabidopsis DRB4, AGO1, AGO7, and RDR6 participate in a DCL4-initiated antiviral RNA silencing pathway negatively regulated by DCL1. Proc Natl Acad Sci U S A. 2008;105:14732–7. doi:10.1073/pnas.0805760105.
Takeda A, Iwasaki S, Watanabe T, Utsumi M, Watanabe Y. The Mechanism Selecting the Guide Strand from Small RNA Duplexes is Different Among Argonaute proteins. Plant Cell Physiol. 2008;49:493–500. doi:10.1093/pcp/pcn043.
Wang XB, Jovel J, Udomporn P, Wang Y, Wu Q, Li WX, et al. The 21-Nucleotide, but Not 22-Nucleotide, Viral Secondary Small Interfering RNAs Direct Potent Antiviral Defense by Two Cooperative Argonautes in Arabidopsis thaliana. Plant Cell. 2011;23:1625–38. doi:10.1105/tpc.110.082305.
Zhang X, Niu D, Carbonell A, Wang A, Lee A, Tun V, et al. ARGONAUTE PIWI domain and microRNA duplex structure regulate small RNA sorting in Arabidopsis. Nat Commun. 2014;5:5468.
Lacombe S, Bangratz M, Vignols F, Brugidou C. The Rice yellow mottle virus P1 protein exhibits dual functions to suppress and activate gene silencing. Plant J Cell Mol Biol. 2010;61:371–82. doi:10.1111/j.1365-313X.2009.04062.x.
Wu J, Yang Z, Wang Y, Zheng L, Ye R, Ji Y, et al. Viral-inducible Argonaute18 confers broadspectrum virus resistance in rice by sequestering a host microRNA. eLife. 2015;4. doi:10.7554/eLife.05733.
Jiang L, Qian D, Zheng H, Meng LY, Chen J, Le WJ, et al. RNA-dependent RNA polymerase 6 of rice (Oryza sativa) plays role in host defense against negative-strand RNA virus, Rice stripe virus. Virus Res. 2012;163:512–9. doi:10.1016/j.virusres.2011.11.016.
Hong W, Qian D, Sun R, Jiang L, Wang Y, Wei C, et al. OsRDR6 plays role in host defense against double-stranded RNA virus, Rice Dwarf Phytoreovirus. Sci Rep. 2015;5:11324. doi:10.1038/srep11324.
Navarro B, Gisel A, Rodio ME, Delgado S, Flores R, Di Serio F. Small RNAs containing the pathogenic determinant of a chloroplast-replicating viroid guide the degradation of a host mRNA as predicted by RNA silencing. Plant J Cell Mol Biol. 2012;70:991–1003. doi:10.1111/j.1365-313X.2012.04940.x.
Adkar-Purushothama CR, Brosseau C, Giguere T, Sano T, Moffett P, Perreault JP. Small RNA Derived from the Virulence Modulating Region of the Potato spindle tuber viroid Silences callose synthase Genes of Tomato Plants. Plant Cell. 2015;27:2178–94. doi:10.1105/tpc.15.00523.
Xia Z, Peng J, Li Y, Chen L, Li S, Zhou T, et al. Characterization of Small Interfering RNAs Derived from Sugarcane Mosaic Virus in Infected Maize Plants by Deep Sequencing. PLoS One. 2014;9:e97013. doi:10.1371/journal.pone.0097013.
Xia Z, Zhao Z, Chen L, Li M, Zhou T, Deng C, et al. Synergistic infection of two viruses MCMV and SCMV increases the accumulations of both MCMV and MCMV-derived siRNAs in maize. Sci Rep. 2016;6:20520. doi:10.1038/srep20520.
Zhou GH, Wen JJ, Cai DJ, Li P, Xu DL, Zhang SG. Southern rice black-streaked dwarf virus: a new proposed Fijivirus species in the family Reoviridae. Chin Sci Bull. 2008;53:3677–85.
Zhou G, Xu D, Zhang M. Southern rice black-streaked dwarf virus: a white-backed planthopper-transmitted fijivirus threatening rice production in Asia. Front Microbiol. 2013;4:270. doi:10.3389/fmicb.2013.00270.
Wang Q, Yang J, Zhou GH, Zhang HM, Chen JP, Adams MJ. The Complete Genome Sequence of Two Isolates of Southern Rice black-streaked dwarf virus, a New Member of the Genus Fijivirus. J Phytopathol. 2010;158:733–7. doi:10.1111/j.1439-0434.2010.01679.x.
Zhang HM, Chen JP, Adams MJ. Molecular Characterisation of Segments 1 to 6 of Rice Black-Streaked Dwarf Virus from China Provides the Complete Genome. Arch Virol. 2001;146:2331–9.
Lu YH, Zhang JF, Xiong RY, Xu QF, Zhou YJ. Identification of an RNA Silencing Suppressor Encoded by Southern rice black-streaked dwarf virus S6. Sci Agric Sin. 2011;44:2909–17.
Li J, Xue J, Zhang HM, Yang J, Lv MF, Xie L, et al. Interactions between the P6 and P5-1 proteins of southern rice black-streaked dwarf fijivirus in yeast and plant cells. Arch Virol. 2013;158:1649–59. doi:10.1007/s00705-013-1660-4.
Wang Q, Tao T, Zhang Y, Wu W, Li D, Yu J, et al. Rice black-streaked dwarf virus P6 self-interacts to form punctate, viroplasm-like structures in the cytoplasm and recruits viroplasm-associated protein p 9–1. Virol J. 2011;8:24. doi:10.1186/1743-422X-8-24.
Mao Q, Zheng S, Han Q, Chen H, Ma Y, Jia D, et al. New Model for the Genesis and Maturation of Viroplasms Induced by Fijiviruses in Insect Vector Cells. J Virol. 2013;87:6819–28. doi:10.1128/JVI.00409-13.
Liu Y, Jia D, Chen H, Chen Q, Xie L, Wu Z, et al. The P7-1 protein of southern rice black-streaked dwarf virus, a fijivirus, induces the formation of tubular structures in insect cells. Arch Virol. 2011;156:1729–36. doi:10.1007/s00705-011-1041-9.
Jia D, Mao Q, Chen H, Wang A, Liu Y, Wang H, et al. Virus-Induced Tubule: a Vehicle for Rapid Spread of Virions through Basal Lamina from Midgut Epithelium in the Insect Vector. J Virol. 2014;88:10488–500. doi:10.1128/JVI.01261-14.
Jia D, Chen H, Zheng A, Chen Q, Liu Q, Xie L, et al. Development of an Insect Vector Cell Culture and RNA Interference System to Investigate the Functional Role of Fijivirus Replication Protein. J Virol. 2012;86:5800–7. doi:10.1128/JVI.07121-11.
Xu D, Mou G, Wang K, Zhou G. MicroRNAs responding to southern rice black-streaked dwarf virus infection and their target genes associated with symptom development in rice. Virus Res. 2014;190:60–8. doi:10.1016/j.virusres.2014.07.007.
Xu D, Zhong T, Feng W, Zhou G. Tolerance and responsive gene expression of Sogatella furcifera under extreme temperature stresses are altered by its vectored plant virus. Sci Rep. 2016;6:31521. doi:10.1038/srep31521.
Lan H, Chen H, Liu Y, Jiang C, Mao Q, Jia D, et al. Small Interfering RNA Pathway Modulates Initial Viral Infection in Midgut Epithelium of Insect after Ingestion of Virus. J Virol. 2016;90:917–29. doi:10.1128/JVI.01835-15.
Chang ZX, Tang N, Wang L, Zhang LQ, Akinyemi IA, Wu QF. Identification and characterization of microRNAs in the white-backed planthopper, Sogatella furcifera. Insect Sci. 2016;23:452–68. doi:10.1111/1744-7917.12343.
Yoshida S, Forno DA, Gomez KA. Laboratory Manual for Physiological Studies of Rice. 3rd ed. Manila: International Rice Research Institute; 1976. p. 61–4.
Tu Z, Ling B, Xu D, Zhang M, Zhou G. Effects of southern rice black-streaked dwarf virus on the development and fecundity of its vector, Sogatella furcifera. Virol J. 2013;10:145. doi:10.1186/1743-422X-10-145.
Dai X, Zhao PX. psRNATarget: a plant small RNA target analysis server. Nucleic Acids Res. 2011;39(Web Server issue):W155–9. doi:10.1093/nar/gkr319.
Conesa A, Götz S. Blast2GO: A Comprehensive Suite for Functional Analysis in Plant Genomics. Int J Plant Genomics. 2008;2008:1–13. doi:10.1155/2008/619832.
Livak KJ, Schmittgen TD. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCt Method. Methods. 2001;25:402–8. doi:10.1006/meth.2001.1262.
Moffett P, Farnham G, Peart J, Baulcombe DC. Interaction between domains of a plant NBS-LRR protein in disease resistance-related cell death. EMBO J. 2002;21:4511–9.
DeYoung BJ, Innes RW. Plant NBS-LRR proteins in pathogen sensing and host defense. Nat Immunol. 2006;7:1243–9. doi:10.1038/ni1410.
McHale L, Tan X, Koehl P, Michelmore RW. Plant NBS-LRR proteins: adaptable guards. Genome Biol. 2006;7:212. doi:10.1186/gb-2006-7-4-212.
Angel CA, Hsieh YC, Schoelz JE. Comparative Analysis of the Capacity of Tombusvirus P22 and P19 Proteins to Function as Avirulence Determinants in Nicotiana species. Mol Plant Microbe Interact. 2011;24:91–9. doi:10.1094/MPMI-04-10-0089.
de Ronde D, Butterbach P, Lohuis D, Hedil M, van Lent JW, Kormelink R. Tsw gene-based resistance is triggered by a functional RNA silencing suppressor protein of the Tomato spotted wilt virus. Mol Plant Pathol. 2013;14:405–15. doi:10.1111/mpp.12016.
Kim HS, Delaney TP. Arabidopsis SON1 Is an F-Box Protein That Regulates a Novel Induced Defense Response Independent of Both Salicylic Acid and Systemic Acquired Resistance. Plant Cell. 2002;14:1469–82.
Schumann N, Navarro-Quezada A, Ullrich K, Kuhl C, Quint M. Molecular Evolution and Selection Patterns of Plant F-Box Proteins with C-Terminal Kelch Repeats. Plant Physiol. 2011;155:835–50. doi:10.1104/pp.110.166579.
Jones-Rhoades MW, Bartel DP, Bartel B. MicroRNAs and Their Regulatory Roles in Plants. Annu Rev Plant Biol. 2006;57:19–53. doi:10.1146/annurev.arplant.57.032905.105218.
Earley K, Smith M, Weber R, Gregory B, Poethig R. An endogenous F-box protein regulates ARGONAUTE1 in Arabidopsis thaliana. Silence. 2010;1:15. doi:10.1186/1758-907X-1-15.
Derrien B, Baumberger N, Schepetilnikov M, Viotti C, De Cillia J, Ziegler-Graff V, et al. Degradation of the antiviral component ARGONAUTE1 by the autophagy pathway. Proc Natl Acad Sci U S A. 2012;109:15942–6. doi:10.1073/pnas.1209487109.
Shiu SH, Bleecker AB. Plant Receptor-Like Kinase Gene Family: Diversity, Function, and Signaling. Sci STKE. 2001;2001:re22.
Morris ER, Walker JC. Receptor-like protein kinases: the keys to response. Curr Opin Plant Biol. 2003;6:339–42.
Burdiak P, Rusaczonek A, Witoń D, Głów D, Karpiński S. Cysteine-rich receptor-like kinase CRK5 as a regulator of growth, development, and ultraviolet radiation responses in Arabidopsis thaliana. J Exp Bot. 2015;66:3325–37. doi:10.1093/jxb/erv143.
Isokpehi RD, Simmons SS, Cohly HH, Ekunwe SI, Begonia GB, Ayensu WK. Identification of Drought-Responsive Universal Stress Proteins in Viridiplantae. Bioinform Biol Insights. 2011;5:41–58. doi:10.4137/BBI.S6061.
Sinha P, Pazhamala LT, Singh VK, Saxena RK, Krishnamurthy L, Azam S, et al. Identification and Validation of Selected Universal Stress Protein Domain Containing Drought-Responsive Genes in Pigeonpea (Cajanus cajan L.). Front Plant Sci. 2016;6:1065. doi:10.3389/fpls.2015.01065.
Ishibashi K, Nishikiori M, Ishikawa M. Interactions Between Tobamovirus Replication Proteins and Cellular Factors: Their Impacts on Virus Multiplication. Mol Plant Microbe Interact. 2010;23:1413–9. doi:10.1094/MPMI-04-10-0102.
Peng H, Chun J, Ai TB, Tong YA, Zhang R, Zhao MM, et al. MicroRNA profiles and their control of male gametophyte development in rice. Plant Mol Biol. 2012;80:85–102. doi:10.1007/s11103-012-9898-x.
Llave C. Virus-derived small interfering RNAs at the core of plant-virus interactions. Trends Plant Sci. 2010;15:701–7. doi:10.1016/j.tplants.2010.09.001.
Ding S-W, Lu R. Virus-derived siRNAs and piRNAs in immunity and pathogenesis. Curr Opin Virol. 2012;1:533–44. doi:10.1016/j.coviro.2011.10.028.
Blevins T, Rajeswaran R, Aregger M, Borah BK, Schepetilnikov M, Baerlocher L, et al. Massive production of small RNAs from a non-coding region of Cauliflower mosaic virus in plant defense and viral counter-defense. Nucleic Acids Res. 2011;39:5003–14. doi:10.1093/nar/gkr119.
Wang J, Tang Y, Yang Y, Ma N, Ling X, Kan J, et al. Cotton Leaf Curl Multan Virus-Derived Viral Small RNAs Can Target Cotton Genes to Promote Viral Infection. Front Plant Sci. 2016;7:1162. doi:10.3389/fpls.2016.01162.
Li J, Andika IB, Shen J, Lv Y, Ji Y, Sun L, et al. Characterization of Rice Black-Streaked Dwarf Virus- and Rice Stripe Virus-Derived siRNAs in Singly and Doubly Infected Insect Vector Laodelphax striatellus. PLoS ONE. 2013;8:e66007. doi:10.1371/journal.pone.0066007.
Xu Y, Huang L, Fu S, Wu J, Zhou X. Population Diversity of Rice Stripe Virus-Derived siRNAs in Three Different Hosts and RNAi-Based Antiviral Immunity in Laodelphgax striatellus. PLoS One. 2012;7:e46238. doi:10.1371/journal.pone.0046238.
Li S, Wang H, Zhou G. Synergism Between Southern rice black-streaked dwarf virus and Rice ragged stunt virus Enhances Their Insect Vector Acquisition. Phytopathol. 2014;104:794–9. doi:10.1094/PHYTO-11-13-0319-R.
Wang H, Xu D, Pu L, Zhou G. Southern rice black-streaked dwarf virus Alters Insect Vectors' Host Orientation Preferences to Enhance Spread and Increase Rice ragged stunt virus Co-Infection. Phytopathol. 2014;104:196–201. doi:10.1094/PHYTO-08-13-0227-R.
Liu B, Chen Z, Song X, Liu C, Cui X, Zhao X, et al. Oryza sativa Dicer-like4 Reveals a Key Role for Small Interfering RNA Silencing in Plant Development. Plant Cell. 2007;19:2705–18.
Nagasaki H, Itoh J, Hayashi K, Hibara K, Satoh-Nagasawa N, Nosaka M, et al. The small interfering RNA production pathway is required for shoot meristem initiation in rice. Proc Natl Acad Sci U S A. 2007;104:14867–71. doi:10.1073/pnas.0704339104.
Song X, Wang D, Ma L, Chen Z, Li P, Cui X, et al. Rice RNA-dependent RNA polymerase 6 acts in small RNA biogenesis and spikelet development. Plant J. 2012;71:378–89. doi:10.1111/j.1365-313X.2012.05001.x.
Csorba T, Burgyán J. Antiviral Silencing and Suppression of Gene Silencing in Plants. In: Wang A, Zhou X, editors. Current Research Topics in Plant Virology. Switzerland: Springer International Publishing AG; 2016. p. 1–33.
Alvarado VY, Scholthof HB. AGO2: a new Argonaute compromising plant virus accumulation. Front Plant Sci. 2012;2:112. doi:10.3389/fpls.2011.00112.
We thank Jingxian Xue for her excellent work in growth chamber and lab maintenance.
This research was finacilly supported by the National Natural Science Foundation of China (31272012) and the Special Fund for Agroscientific Research in the Public Interest (201303021).
Availability of data and materials
GZ conceived and designed the experiments and revised the manuscript. DX performed the experiments, analyzed the data, and wrote the paper. Both authors read and approved the final manuscript.
The authors have no competing interests to report.
Consent for publication
Ethics approval and consent to participate
List of the RT-qPCR primers used in this study. (XLSX 10 kb)
List of the SRBSDV-derived siRNAs identified by deep sequencing from infected rice samples at 14 dpi. (XLSX 70 kb)
List of the SRBSDV-derived siRNAs identified by deep sequencing from infected rice samples at 28 dpi. (XLSX 1015 kb)
List of the SRBSDV-derived siRNAs identified at both 14 dpi and 28 dpi. (XLSX 28 kb)
List of counts of the nucleotides in SRBSDV-derived siRNAs against nt position on the viral genome. (XLSX 433 kb)
List of putative target genes of SRBSDV-derived siRNAs of >40 reads predicted by psRNATarget. (XLSX 78 kb)