Rice stripe virus-derived siRNAs play different regulatory roles in rice and in the insect vector Laodelphax striatellus
BMC Plant Biology volume 18, Article number: 219 (2018)
Most plant viruses depend on vector insects for transmission. Upon viral infection, virus-derived small interfering RNAs (vsiRNAs) can target both viral and host transcripts. Rice stripe virus (RSV) is a persistent-propagative virus transmitted by the small brown planthopper (Laodelphax striatellus, Fallen) and can cause a severe disease on rice.
To investigate how vsiRNAs regulate gene expressions in the host plant and the insect vector, we analyzed the expression profiles of small RNAs (sRNAs) and mRNAs in RSV-infected rice and RSV-infected planthopper. We obtained 88,247 vsiRNAs in rice that were predominantly derived from the terminal regions of the RSV RNA segments, and 351,655 vsiRNAs in planthopper that displayed relatively even distributions on RSV RNA segments. 38,112 and 80,698 unique vsiRNAs were found only in rice and planthopper, respectively, while 14,006 unique vsiRNAs were found in both of them. Compared to mock-inoculated rice, 273 genes were significantly down-regulated genes (DRGs) in RSV-infected rice, among which 192 (70.3%) were potential targets of vsiRNAs based on sequence complementarity. Gene ontology (GO) analysis revealed that these 192 DRGs were enriched in genes involved in kinase activity, carbohydrate binding and protein binding. Similarly, 265 DRGs were identified in RSV-infected planthoppers, among which 126 (47.5%) were potential targets of vsiRNAs. These planthopper target genes were enriched in genes that are involved in structural constituent of cuticle, serine-type endopeptidase activity, and oxidoreductase activity.
Taken together, our results reveal that infection by the same virus can generate distinct vsiRNAs in different hosts to potentially regulate different biological processes, thus reflecting distinct virus-host interactions.
Many plant viruses, especially persistent-propagative plant viruses, infect a wide range of economical crops, including fruits, vegetables and food crops, representing a serious threat to agriculture production . Approximately 80% of plant viruses depend on insect vectors for transmission [2,3,4]. This insect dependent transmission can be divided into non-persistent, semi-persistent and persistent types, based primarily on the length of successful transmission of virus to the host plant . Persistent-propagative plant viruses can induce serious disease symptoms in agricultural crops , while they generally cause asymptomatic phenotypes in vector insects . Persistent-propagative plant viruses encounter different defense strategies when they replicate in vector insects and host plants . One of the antiviral immunity strategies is the RNA silencing machinery triggered by virus infection. This machinery targets both viral transcripts and host transcripts [9,10,11,12]. Virus-derived small interfering RNAs (vsiRNAs) that are involved in the RNAi-based antiviral response can also potentially alter the host transcriptome [13,14,15,16,17]. However, different defense mechanisms mediated by vsiRNAs in host plants and vector insects are poorly documented.
RNA silencing is a widespread antiviral mechanism [5, 18, 19]. vsiRNAs are processed by Dicer from double-stranded RNAs (dsRNAs) or structured single-stranded RNAs (ssRNAs) [20, 21]. In plants, the Dicer like (DCL) protein 4 (DCL4) primarily processes viral RNAs into 21-nt vsiRNAs. DCL2 is responsible for creation of 22-nt vsiRNAs when DCL4 is inactivated or suppressed [22,23,24]. In contrast, in insects, Dicer-2 and Argonaute 2 (AGO2) have substantial effects on viral replication. Viral dsRNAs accumulate during virus infection and are processed into vsiRNA duplexes by Dicer-2. These vsiRNAs are loaded onto insect AGO2 to target viral transcripts [25,26,27].
Rice stripe virus (RSV) causes rice stripe disease, one of the most destructive rice diseases in Eastern Asia [28, 29]. Since first discovered in 1963 in China, this disease has circulated widely and outbroken many times, causing up to 50% grain loss . In 2004, the infected area in Jiangsu province was 1,571,000 ha, accounting for 80% of the rice field . RSV belongs to the genus Tenuivirus and consists of four single-stranded RNA segments that encode seven proteins [20, 32, 33]. RSV is a typical persistent-propagative plant virus and is efficiently transmitted via one of the most economically important insects, the small brown planthopper (Laodelphax striatellus Fallen). vsiRNAs were previously identified in RSV-infected O. sativa, N. benthamiana and L. striatellus [21, 34]. There is large number of RSV-derived sRNAs (vsiRNAs) in these libraries, and of these vsiRNAs, O. sativa is the main source followed by N. benthamiana and then L. striatellus.
However, a comparative functional analysis of RSV vsiRNAs in RSV-infected host plants and vector insects was not reported. To gain further insights into the regulatory functions of vsiRNAs in host plants and an efficient insect vector, we studied the profiles of small RNAs (sRNAs) and mRNAs in RSV-infected rice and small brown planthopper. We show that RSV not only generates vsiRNAs that accumulate in both hosts, but also generates host-specific vsiRNAs that potentially control different defensive responses in the two hosts.
sRNA profiles in RSV-infected rice and small brown planthopper
sRNA libraries were constructed and sequenced using sRNAs isolated from rice plants and small brown planthoppers after being infected with RSV for 20 d and 5 d, respectively, with three biological replicates for each group. As controls, mock-inoculated rice and insect samples were prepared and sequenced in parallel. One RSV-infected rice sample and one mock-inoculated rice sample were discarded due to poor quality and only two replicates for each group were retained for analyses. After removal of adapter sequences, and low quality reads, we obtained, on average, 15,400,564, 13,876,277, 15,309,080, and 14,117,003 clean sRNA reads from RSV-infected rice, mock-inoculated rice, RSV-infected planthopper, and mock-inoculated planthopper samples, respectively (Additional file 1 and Additional file 2). After discarding reads mapped to structural RNAs (rRNAs, tRNAs, snoRNAs, snRNAs), the remaining reads were mapped perfectly to the host genome, and 9,208,240, 9,052,920, 6,727,312, and 5,842,694 host-derived sRNAs were obtained on average. These sRNAs were clustered into, on average, 2,468,202, 2,235,158, 1,330,659 and 1,141,043 unique sRNA reads.
The size distributions were quite different for rice- and insect-derived sRNAs. For rice sRNAs, the dominant sizes were 21-nt and 24-nt (Fig. 1a and Additional file 3), whereas the dominant sizes of sRNAs in planthopper were 22-nt and 26-nt to 28-nt (Fig. 1b and Additional file 3), reflecting different sRNA processing properties in these two hosts. The GC contents of sRNAs of different sizes in rice and planthopper samples were calculated. Rice-derived sRNAs had a higher average GC content than did insect-derived sRNAs (Additional file 4).
Because the 5′ terminal nucleotides of sRNAs are important for AGO sorting and function, we analyzed the 5′ nucleotide composition of sRNAs in each library. In rice, uridine (U) or adenosine (A) was the dominant 5′ terminal nucleotide for the 21-nt and 24-nt sRNAs, respectively (Fig. 1a and Additional file 3). In planthopper, U was the dominant nucleotide at the 5′ terminus for the 22-nt and 26-nt to 28-nt sRNAs (Fig. 1b and Additional file 5). In contrast, U and A were the most abundant 5′ terminal nucleotides for rice total sRNAs, while U was still the dominant 5′ terminal nucleotide for planthopper total sRNAs (Fig. 1c, d and Additional file 6).
We compared the abundance of host sRNAs in RSV-infected samples and mock-inoculated samples to identify host sRNAs that were regulated upon RSV infection. Compared to the mock-inoculated samples, 450 and 1,558 sRNAs were up- and down-regulated in RSV-infected rice, and 1,776 and 783 sRNAs were up- and down-regulated in RSV-infected planthopper (Fig. 1e and f).
RSV-derived siRNAs in RSV-infected rice and small brown planthopper
We mapped clean sRNA reads to the RSV genome and obtained RSV-derived siRNAs. We identified 38,122 and 80,698 unique vsiRNAs in RSV-infected rice and viruliferous planthopper, respectively, and 14,006 unique vsiRNAs were found in both hosts (Fig. 2a). In RSV-infected rice, the size of vsiRNAs ranged from 18-nt to 30-nt, with a peak at 21-nt (Fig. 2b, Additional file 7 and Additional file 8). In RSV-infected planthoppers, the size of vsiRNAs ranged predominantly from 18-nt to 24-nt, with a peak at 22-nt (Fig. 2c, Additional file 7 and Additional file 8). In addition, the preferential occurrence of A/U at the 5′ terminus was identified for the vsiRNAs in rice and planthoppers (Fig. 2d and Additional file 9).
We also looked in detail at the strand specificity and locations of origin of vsiRNAs. As shown in Fig. 2e, f and Additional file 10, the proportions of vsiRNAs that were generated from sense and antisense genome sequences were comparable in both hosts. Regarding the genomic regions from which vsiRNAs were generated, far fewer vsiRNAs were produced from RNA1 than from the other three RNA segments in rice (Fig. 2e). In contrast, RNA3 generated more vsiRNAs than the other three segments in planthopper (Fig. 2f). Moreover, accumulation peaks in the termini of viral RNA segments were observed in rice, but the vsiRNAs were more evenly distributed in planthopper (Fig. 2e, f).
The 8 most abundant vsiRNAs in rice and planthopper were selected for RT-PCR validation, respectively. The results showed that these vsiRNAs accumulated highly in RSV-infected rice/planthopper, demonstrating that RSV infection indeed generated highly abundant vsiRNAs in both host plants and vector insects (Fig. 3).
Transcriptome alterations of rice and planthopper upon RSV infection
To explore possible regulation mechanisms mediated by vsiRNAs in rice and small brown planthopper, we analyzed their transcriptomic responses to RSV infection. The same samples used for sRNA sequencing were subjected to mRNA sequencing and differentially expressed genes (DEGs) between RSV-infected and mock-inoculated samples were identified.
After infection with RSV, a total of 989 DEGs were identified in rice, among which 273 genes were down-regulated (Additional file 11 and Additional file 12). GO term enrichment analysis showed that 14 GO terms were enriched for the down-regulated genes in rice, including: response to stress, cell death, metabolic process, cell wall, extracellular region and cell as the most significantly enriched (p-value < 0.01) (Fig. 4a). Furthermore, 48 down-regulated genes were potential targets of vsiRNAs according to sequence complementarity. These target genes were enriched with genes involved in: kinase activity (p-value < 0.01), carbohydrate binding, protein binding, extracellular region, multi-organism reproductive process, pollen-pistil interaction, and single-organism cellular process (p-value < 0.05) (Fig. 4b). In the 716 up-regulated genes of rice, the most significantly enriched GO terms were response to biotic stimulus or stress, metabolic process, oxygen binding, catalytic activity and carbohydrate binding (p-value < 0.01) (Additional file 13).
In RSV-infected planthopper, 551 DEGs were identified, among which 265 genes were down-regulated (Additional file 11 and Additional file 14). 20 GO terms were enriched for the down-regulated genes: extracellular matrix and region, structural constituent of chitin-based larval cuticle, chitin binding, odorant binding, chitin-based cuticle development, and chitin metabolic process (Fig. 5a).
Seventy down-regulated genes were predicted to be the targets of vsiRNAs. These target genes were enriched in genes involved in: cuticle-related molecular functions or processes, signal recognition particle, extracellular matrix and region, oxidoreductase activity acting on paired donors, developmental pigmentation, SRP-dependent cotranslational protein targeting to membrane, and organic hydroxy compound biosynthetic process (Fig. 5b). In the 286 up-regulated genes, genes that were significantly enriched were involved in: response to gap junction, iron ion binding, cellular calcium ion homeostasis, and hexose catabolic process (Additional file 15).
Host-specific vsiRNAs targeted kinase activity genes in rice and pathogen/disease resistance genes in planthoppers
We then investigated the targets of host specific vsiRNAs in rice and planthoppers. Among the 48 potential targets of vsiRNAs in rice, 17 genes were found to be targets of rice specific vsiRNAs (Fig. 4c and Additional file 16). Thirteen out of the 17 genes had a kinase activity, such as protein kinases and receptor protein kinases. The other four potential targets took part in carbohydrate binding and pollen-pistil interaction, such as D-mannose binding lectin family protein and serine/threonine-protein kinase receptor precursor. These 17 genes were predicted to be targeted by from one to three rice-specific vsiRNAs (Fig. 4c). Six of these 17 genes were selected to validate differential expression patterns between RSV-infected and mock-inoculated rice using qRT-PCR. All of the 6 genes showed significantly lower transcript levels in RSV-infected rice than in mock-inoculated rice (Fig. 4d). These results are consistent with the transcriptome data.
In the 70 potential targets of vsiRNAs in planthoppers, 6 genes were found to be targets of planthopper specific vsiRNAs (Fig. 5c and Additional file 17). Some of these encoded pathogen/disease resistance proteins, such as trypsins, serpin-5, and peritrophic matrix protein (Fig. 5c). These 6 genes were regulated by from one to three planthopper specific vsiRNAs. Among the 70 genes, we investigated 6 potential targets, among which 4 had lower transcript levels in RSV-infected planthoppers and 2 had no significant differences between RSV-infected and mock-inoculated groups (Fig. 5d).
Common vsiRNAs regulated different targets in rice and planthopper
To determine whether same vsiRNAs can regulate similar pathways in host rice and vector planthopper upon RSV infection, we selected the 50 most abundant vsiRNAs from the 14,006 vsiRNAs that were expressed in rice and planthopper to predict their potential targets. A total of 39 and 58 down-regulated genes were predicted to have potential binding sites for one or more of the 50 vsiRNAs in rice and planthopper, respectively. GO analysis of the 39 predicted rice target genes showed that many of them encoded immune response and defense proteins that are involved in cell death and response to stress and biotic stimulus (Fig. 6a). In contrast, the 58 predicted planthopper target genes were most enriched in GO terms of extracellular matrix, structural constituent of cuticle, structural constituent of chitin−based larval cuticle, and chitin−based cuticle development (Fig. 6b). Therefore, it appears that RSV vsiRNAs may manipulate distinct biological processes in host plant and vector insect. To uncover whether vsiRNAs can target evolutionarily conserved genes in these two hosts, we used the 39 rice and 58 planthopper potential target genes as queries to search for homologous genes in the other host. We found that 11 rice target genes showed > 30% identity at the amino acid level with 12 planthopper target genes. However, these homolog-like target genes did not have similar functional annotations (Fig. 6c), indicating that they may perform different functions.
In this paper, we have examined the expression profiles of vsiRNAs in RSV-infected rice and planthopper, and showed that RSV generated either same or distinct vsiRNAs, potentially regulating different biological processes in the two hosts. A few studies have shown that vsiRNAs, as pathogenicity determinants, negatively regulate host mRNAs and effectively silence host genes [15, 35,36,37,38,39,40]. For example, chloroplast-related genes can be targeted by vsiRNA to regulate the symptoms of twisted-leaves and stunting . Our study suggests that a given virus can generate distinct vsiRNAs in different hosts to regulate different biological processes, reflecting distinct virus-host interactions.
The mechanisms of vsiRNA production may be different in plant and insect. We found that the 21-nt vsiRNA was the main size-class of vsiRNAs in rice while the 22-nt vsiRNA turned to be the most abundant class of vsiRNAs in planthopper, which was consistent with the findings from another study . Different sizes of vsiRNAs may be processed by different Dicers/DCLs. DCL4 is known to be mainly responsible for production of 21-nt vsiRNAs and DCL2 is responsible for the production of 22-nt vsiRNAs for rescuing silence against RNA viruses when DCL4 was inactivated or suppressed [22,23,24, 41]. Further investigation is needed to elucidate the role of rice DCL4 and DCL2 in producing RSV vsiRNAs.
Our study revealed a preferential occurrence of uridine or adenosine residues, as compared to cytosine and guanidine, at the 5′ terminus of vsiRNAs in rice and planthopper. Previous studies have indicated that the 5′-terminal nucleotides of sRNAs control the sorting of sRNAs to specific AGO complexes in plant [42, 43]. In Drosophila, sRNAs were sorted for AGO1 or AGO2 according to their duplex structures and 5′-terminal nucleotides [44, 45]. The 5′-terminal nucleotides of vsiRNAs from rice and planthopper were similar, suggesting that these nucleotides of 5′-terminal may have role(s) in loading of vsiRNAs into specific AGO complexes.
The hot spots of RSV-derived vsiRNA in rice are predominantly at the terminal regions of the RSV RNA segments. This phenomenon is similar to that of Cucumber mosaic virus (CMV) vsiRNAs, which are also enriched in the terminal regions of CMV genomes. Moreover, this enrichment is lost in RNA-dependent RNA polymerases (rdrs) mutant plants infected with CMV. Thus RDRs are supposed to be involved in amplification of vsiRNAs from CMV genomic RNAs . The enrichment of RSV-derived vsiRNAs in the terminal regions in rice may also originate from the host RDR-mediated vsiRNAs amplification.
We find that rice-derived sRNAs tend have a higher average GC content than do insect-derived sRNAs. The sRNA structural features of different GC values may affect accessibility, stability, or activity of the components in sRNAs biogenesis pathway. Higher GC content around the sRNA generating site may result in forming a relatively more stable secondary structure, which would be a preferable substrate for DCL processing . The difference of GC content of insect-derived sRNAs and rice-derived sRNAs point to the different sRNAs biogenesis pathway in insect and rice.
vsiRNA acts as an effective cross-talk molecule to defend the stress signaling components of hosts after RSV infection. The potential cross reaction of vsiRNAs with host genes is a systematic counter-offensive strategy deployed by virus to shift the virus/host ‘tug of war’ in the pathogen’s favor. Some pathogen/disease resistance genes such as genes encoding trypsins and peritrophic matrix protein are down-regulated by vsiRNAs in planthopper [48, 49]. The down-regulation of these resistance genes would favor RSV infection and replication in the planthopper. In rice, some crucial genes, such as those encoding protein kinases, mannose binding lectin family protein, act as hubs in the regulatory network in plant response to abiotic and biotic stress. Protein kinases are involved in defense of different plant virus infections in tomato and tobacco [50,51,52,53]. Mannose binding lectin genes are critical in developing resistance to various pathogenic organisms and insects [54,55,56,57]. Down-regulation of genes encoding protein kinases and mannose binding lectin family proteins by vsiRNAs may facilitate RSV pathogenicity in rice.
Some vsiRNA-derived sRNA common in host rice and vector planthopper were found to potentially target one or several host factors involved in various physiological pathways of different host. Similarly, vsiRNA derived from grapevine fleck virus (GFkV) and Grapevine rupestris stem pitting-associated virus (GRSPaV) were shown to target Vitis vinifera derived genes involved in ribosome biogenesis, biotic and abiotic stresses besides their role in plant’s defense mechanism . Accordingly, the targets of the corresponding common vsiRNAs from both host plant and vector insect are diverse not only in sequence complementarity but on the positions of the vsiRNA targets in mRNA, secondary structural features of the target mRNA , which could also play a significant role in small RNA-targeted gene silencing.
Virus, vector and host plant preparation
The viruliferous and non-viruliferous small brown planthopper strains used in this work were derived from a rice field population in Hai’an, Jiangsu Province, China, in 2009, and maintained in the laboratory until now. The lengthy period of maintenance and screening make sure that no other pathogens exist in the planthopper strains. The planthoppers were reared on 2 cm to 3 cm seedlings of rice (Oryza sativa Huangjinqing) in glass incubators, which were sealed with a nylon mesh, at 25 ± 1 °C, 80% relative humidity, and 16 h of light per day. The planthoppers were transferred to fresh rice seedlings every 8 days to ensure sufficient nutrition. The RSV-carrying rate of the viruliferous planthopper strains were maintained at no less than 90% through a purification selection every three months using a dot-ELISA assay with the monoclonal anti-Cp antibody .
Viruliferous fourth-instar planthoppers were transferred to fresh rice seedlings for RSV inoculation. After 24 h feeding, the planthoppers were removed and the inoculated seedlings were cultured for 20 d until disease symptom appeared. Each rice seedling was fed by three viruliferous planthoppers, and three seedlings were used in a replicate experiment. Healthy rice seedlings without inoculation were prepared and maintained at the same condition as the negative control. Leaves with typical disease symptoms from RSV-infected rice plants and leaves from healthy rice plants were collected at the same time and three biological replicates for each group were prepared for sequencing analysis.
Non-viruliferous fourth-instar planthopper nymphs were fed on an artificial diet containing the RSV crude preparations from the infected rice leaves as previously described . After feeding on RSV for 8 h, the nymphs were transferred to healthy rice seedlings and then collected after 5 d. Three biological replicates and 20 insects per replicate were prepared for sequencing analysis.
RNA isolation and sequencing
Total RNA was isolated from rice plants or planthoppers using the TRIzol Reagent (Invitrogen, Carlsbad, CA, USA) according to manufacturer’s instructions. Approximately 100 mg of rice leaves or 10 small brown planthoppers were ground in liquid nitrogen and mixed with 1 ml of TRIzol Reagent for RNA extraction. The concentration and quality of RNA were measured using NanoDrop spectrophotometer (Thermo Scientific, Waltham, MA, USA) and gel electrophoresis. Total RNA was used for mRNA library and sRNA library constructions. mRNA sequencing was performed on an Illumina Hiseq 2500 sequencer to generate paired-end 125 bp data and the sRNA sequencing was conducted on an Illumina Hiseq 4000 sequencer to generate single-end 50 bp sequences in Beijing Genome Institute (BGI, Shenzhen, China).
The rice (Oryza sativa ssp. japonica cv. Nipponbare) genome sequence (v. 7.0) and its gene annotations, including Gene Ontology (GO) annotations, were downloaded from the Rice Genome Annotation Project (http://rice.plantbiology.msu.edu). The genome sequence of the small brown planthopper and GO annotations have been recently published by our laboratory .
Analysis of sRNA sequencing data
The adapter sequence was removed using cutadapt (v1.16) . Reads that were 18 nt–30 nt and contain no ambiguous nucleotide (‘N’) were retained as clean reads, and were clustered into unique reads with an in-house python script. For rice sRNA datasets, clean reads were mapped to the ncRNA reference (ftp://ftp.ensemblgenomes.org/pub/plants/release-31/fasta/oryza_sativa/ncrna/Oryza_sativa.IRGSP-1.0.31.ncrna.fa.gz) downloaded from ensemble database to remove reads mapped to known rice rRNAs, tRNAs, snRNAs using bowtie software  with perfect matching. For planthopper sRNA datasets, clean reads mapped to the rRNAs, tRNAs, snoRNAs, snRNAs resolved by our lab  with perfect matching were discarded. The remaining sRNA reads were then mapped to the host genome (rice or planthopper genome) with perfect matching or to RSV genome, allowing up to one mismatch. Only sRNAs with at least one mapping location were regarded as genuine sRNAs and retained for down-stream analyses. EdgeR [63, 64] was used to identify differentially expressed sRNAs between RSV-infected rice or planthopper samples and the corresponding mock-inoculated samples using 2-fold change in expression and FDR < 0.05 as significance cutoffs. sRNAs with a match to the RSV genome (and without a match to the host genome) were regarded as vsiRNAs. We used psRobot  to identify potential target genes of vsiRNAs in rice, and used Miranda  to identify potential target genes of vsiRNAs in the planthopper.
Analysis of mRNA sequencing data
The quality of raw RNA sequencing reads was evaluated with FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc). To restrict to clean reads, low quality regions and adapter sequences were removed using SolexaQA (v2.2)  and cutadapt (v1.16), respectively. Reads shorter than 25 bases were discarded. Clean reads were mapped to the host (rice or planthopper, depending on the sample source) genome sequence using HISAT2 . Read counts of the annotated genes were summarized by HTSeq-count . EdgeR  was used to identify differentially expressed genes. The genes with at least a 2-fold change in expression level and an FDR < 0.05 were considered to be differentially expressed. GO term enrichment analysis was performed using the R package TopGO . Adrian Alexi’s improved weighted scoring algorithm and Fisher’s test were used to determine the significance of GO term enrichment. The significantly enriched GO terms were inferred with a cutoff p-value less than 0.05.
Reverse transcription PCR for vsiRNAs
Total RNAs was isolated from the viruliferous and nonviruliferous rice/planthoppers using TRIzol Reagent (Invitrogen, Carlsbad, CA, USA). The poly (A) was added to the 3′ end of sRNAs by E. coli Poly(A) Polymerase (NEB, Ipswich, MA, USA). The first-strand cDNAs of sRNAs were synthesized using Moloney murine leukemia virus (M-MLV) reverse transcriptase (NEB). The vsiRNAs were subjected to reverse transcription PCR (RT-PCR) using the SYBR Green sRNA expression assays. RT-PCR was amplified for 25–28 cycles. U6 snRNA was amplified as endogenous control. The RT-PCR primers are listed in Additional file 18.
In summary, we report here on virus/host interactions through vsiRNA-mRNA potential regulation in RSV-infected rice and planthoppers. We have demonstrated the possibility of vsiRNA-mediated downregulation of host genes using qPCR. Thus, it may be possible to use vsiRNA as effector molecules of silencing in artificial miRNA-mediated antiviral resistance. Further studies are needed to verify this interpretation.
Cucumber mosaic virus
differentially expressed genes
Moloney murine leukemia virus
National Center for Biotechnology Information
RNA-dependent RNA polymerases
Rice stripe virus
reverse transcription PCR
virus-derived small interfering RNAs
Soosaar JL, Burch-Smith TM, Dinesh-Kumar SP. Mechanisms of plant resistance to viruses. Nat Rev Microbiol. 2005;3(10):789–98.
Ng JC, Perry KL. Transmission of plant viruses by aphid vectors. Mol Plant Pathol. 2004;5(5):505–11.
Power AG. Insect transmission of plant viruses: a constraint on virus variability. Curr Opin Plant Biol. 2000;3(4):336–40.
Hohn T. Plant virus transmission from the insect point of view. Proc Natl Acad Sci U S A. 2007;104(46):17905–6.
Ding SW. RNA-based antiviral immunity. Nat Rev Immunol. 2010;10(9):632–44.
Gilbertson RL, Batuman O, Webster CG, Adkins S. Role of the insect supervectors Bemisia tabaci and Frankliniella occidentalis in the emergence and global spread of plant viruses. Annu Rev Virol. 2015;2(1):67–93.
Ramirez BC, Haenni AL. Molecular biology of tenuiviruses, a remarkable group of plant viruses. J Gen Virol. 1994;75(Pt 3):467–75.
de Haro LA, Dumon AD, Mattio MF, Arguello Caro EB, Llauger G, Zavallo D, Blanc H, Mongelli VC, Truol G, Saleh MC, et al. Mal de Rio Cuarto virus infection triggers the production of distinctive viral-derived siRNA profiles in wheat and its Planthopper vector. Front Plant Sci. 2017;8:766.
Ghoshal B, Sanfacon H. Symptom recovery in virus-infected plants: revisiting the role of RNA silencing mechanisms. Virology. 2015;479-480:167–79.
Ramesh SV, Williams S, Kappagantu M, Mitter N, Pappu HR. Transcriptome-wide identification of host genes targeted by tomato spotted wilt virus-derived small interfering RNAs. Virus Res. 2017;238:13–23.
Ratcliff F, Harrison BD, Baulcombe DC. A similarity between viral defense and gene silencing in plants. Science. 1997;276(5318):1558–60.
Huang J, Yang ML, Zhang XM. The function of small RNAs in plant biotic stress response. J Integr Plant Biol. 2016;58(4):312–27.
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(2):203–14.
Dunoyer P, Voinnet O. The complex interplay between plant viruses and host RNA-silencing pathways. Curr Opin Plant Biol. 2005;8(4):415–23.
Xia Z, Zhao Z, Chen L, Li M, Zhou T, Deng C, Zhou Q, Fan Z. 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.
Li L, Guo C, Wang B, Zhou T, Lei Y, Dai YH, He W, Liang C, Wang XF. RNAi-mediated transgenic rice resistance to Rice stripe virus. J Integr Agr. 2016;15(11):2539–49.
Huang J, Yang ML, Lu L, Zhang XM. Diverse functions of small RNAs in different plant-pathogen communications. Front Microbiol. 2016;7:1552.
Ding SW, Voinnet O. Antiviral immunity directed by small RNAs. Cell. 2007;130(3):413–26.
Mlotshwa S, Pruss GJ, Vance V. Small RNAs in viral infection and host defense. Trends Plant Sci. 2008;13(7):375–82.
Du P, Wu JG, Zhang JY, Zhao SQ, Zheng H, Gao G, Wei LP, Li Y. Viral infection induces expression of novel phased MicroRNAs from conserved cellular MicroRNA precursors. PLoS Pathog. 2011;7(8):e1002176.
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(9):e46238.
Bouche N, Lauressergues D, Gasciolli V, Vaucheret H. An antagonistic function for Arabidopsis DCL2 in development and a new function for DCL4 in generating viral siRNAs. EMBO J. 2006;25(14):3347–56.
Deleris A, Gallego-Bartolome J, Bao J, Kasschau KD, Carrington JC, Voinnet O. Hierarchical action and inhibition of plant dicer-like proteins in antiviral defense. Science. 2006;313(5783):68–71.
Fusaro AF, Matthew L, Smith NA, Curtin SJ, Dedic-Hagan J, Ellacott GA, Watson JM, Wang MB, Brosnan C, Carroll BJ, et al. RNA interference-inducing hairpin RNAs in plants act through the viral defence pathway. EMBO Rep. 2006;7(11):1168–75.
Galiana-Arnoux D, Dostert C, Schneemann A, Hoffmann JA, Imler JL. Essential function in vivo for Dicer-2 in host defense against RNA viruses in drosophila. Nat Immunol. 2006;7(6):590–7.
van Rij RP, Saleh MC, Berry B, Foo C, Houk A, Antoniewski C, Andino R. The RNA silencing endonuclease Argonaute 2 mediates specific antiviral immunity in Drosophila melanogaster. Genes Dev. 2006;20(21):2985–95.
Wang XH, Aliyari R, Li WX, Li HW, Kim K, Carthew R, Atkinson P, Ding SW. RNA interference directs innate immunity against viruses in adult Drosophila. Science. 2006;312(5772):452–4.
Falk BW, Tsai JH. Biology and molecular biology of viruses in the genus Tenuivirus. Annu Rev Phytopathol. 1998;36:139–63.
Xiong RY, Wu JX, Zhou YJ, Zhou XP. Identification of a movement protein of the Tenuivirus Rice stripe virus. J Virol. 2008;82(24):12304–11.
Li J, Xiang CY, Yang J, Chen JP, Zhang HM. Interaction of HSP20 with a viral RdRp changes its sub-cellular localization and distribution pattern in plants. Sci Rep. 2015;5:14016.
Wei TY, Yang JG, Liao FR, Gao FL, Lu LM, Zhang XT, Li F, Wu ZJ, Lin QY, Xie LH, et al. Genetic diversity and population structure of rice stripe virus in China. J Gen Virol. 2009;90(Pt 4):1025–34.
Shi BB, Lin L, Wang SH, Guo Q, Zhou H, Rong LL, Li JM, Peng JJ, Lu YW, Zheng HY, et al. Identification and regulation of host genes related to Rice stripe virus symptom production. New Phytol. 2016;209(3):1106–19.
Wang W, Zhao W, Li J, Luo L, Kang L, Cui F. The c-Jun N-terminal kinase pathway of a vector insect is activated by virus capsid protein and promotes viral replication. Elife. 2017;6:e26591.
Yan F, Zhang H, Adams MJ, Yang J, Peng J, Antoniw JF, Zhou Y, Chen J. Characterization of siRNAs derived from rice stripe virus in infected rice plants by deep sequencing. Arch Virol. 2010;155(6):935–40.
Angell SM, Baulcombe DC. Consistent gene silencing in transgenic plants expressing a replicating potato virus X RNA. EMBO J. 1997;16(12):3675–84.
Miozzi L, Gambino G, Burgyan J, Pantaleo V. Genome-wide identification of viral and host transcripts targeted by viral siRNAs in Vitis vinifera. Mol Plant Pathol. 2013;14(1):30–43.
Qi X, Bao FS, Xie Z. Small RNA deep sequencing reveals role for Arabidopsis thaliana RNA-dependent RNA polymerases in viral siRNA biogenesis. PLoS One. 2009;4(3):e4971.
Shimura H, Pantaleo V, Ishihara T, Myojo N, Inaba J, Sueda K, Burgyan J, Masuta C. 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(5):e1002021.
Smith NA, Eamens AL, Wang MB. Viral small interfering RNAs target host genes to mediate disease symptoms in plants. PLoS Pathog. 2011;7(5):e1002022.
Wang J, Tang Y, Yang Y, Ma N, Ling X, Kan J, He Z, Zhang B. Cotton leaf curl Multan virus-derived viral small RNAs can target cotton genes to promote viral infection. Front Plant Sci. 2016;7:1162.
Diaz-Pendon JA, Li F, Li WX, Ding SW. Suppression of antiviral silencing by cucumber mosaic virus 2b protein in Arabidopsis is associated with drastically reduced accumulation of three classes of viral small interfering RNAs. Plant Cell. 2007;19(6):2053–63.
Mi SJ, Cai T, Hu YG, Chen Y, Hodges E, Ni FR, Wu L, Li S, Zhou H, Long CZ, et al. Sorting of small RNAs into Arabidopsis argonaute complexes is directed by the 5 ' terminal nucleotide. Cell. 2008;133(1):116–27.
Montgomery TA, Howell MD, Cuperus JT, Li DW, Hansen JE, Alexander AL, Chapman EJ, Fahlgren N, Allen E, Carrington JC. Specificity of ARGONAUTE7-miR390 interaction and dual functionality in TAS3 trans-acting siRNA formation. Cell. 2008;133(1):128–41.
Khvorova A, Reynolds A, Jayasena SD. Functional siRNAs and miRNAs exhibit strand bias. Cell. 2003;115(2):209–16.
Ghildiyal M, Xu J, Seitz H, Weng ZP, Zamore PD. Sorting of Drosophila small silencing RNAs partitions microRNA* strands into the RNA interference pathway. Rna. 2010;16(1):43–56.
Wang XB, Wu QF, Ito T, Cillo F, Li WX, Chen XM, Yu JL, Ding SW. RNAi-mediated viral immunity requires amplification of virus-derived siRNAs in Arabidopsis thaliana. P Natl Acad Sci USA. 2010;107(1):484–9.
Ho T, Wang H, Pallett D, Dalmay T. Evidence for targeting common siRNA hotspots and GC preference by plant dicer-like proteins. FEBS Lett. 2007;581(17):3267–72.
Agrawal S, Kelkenberg M, Begum K, Steinfeld L, Williams CE, Kramer KJ, Beeman RW, Park Y, Muthukrishnan S, Merzendorfer H. Two essential peritrophic matrix proteins mediate matrix barrier functions in the insect midgut. Insect Biochem Mol Biol. 2014;49:24–34.
Kato N, Mueller CR, Fuchs JF, Mcelroy K, Wessely V, Higgs S, Christensen BM. Evaluation of the function of a type I Peritrophic matrix as a physical barrier for midgut epithelium invasion by mosquito-borne pathogens in Aedes aegypti. Vector-Borne Zoonot. 2008;8(5):701–12.
Gullner G, Juhasz C, Nemeth A, Barna B. Reactions of tobacco genotypes with different antioxidant capacities to powdery mildew and tobacco mosaic virus infections. Plant Physiol Biochem. 2017;119:232–9.
van Loon LC, Rep M, Pieterse CMJ. Significance of inducible defense-related proteins in infected plants. Annu Rev Phytopathol. 2006;44:135–62.
Huang Y, Li T, Xu ZS, Wang F, Xiong AS. Six NAC transcription factors involved in response to TYLCV infection in resistant and susceptible tomato cultivars. Plant Physiol Biochem. 2017;120:61–74.
Huang Y, Li MY, Wu P, Xu ZS, Que F, Wang F, Xiong AS. Members of WRKY group III transcription factors are important in TYLCV defense signaling pathway in tomato (Solanum lycopersicum). BMC Genomics. 2016;17(1):788.
Mondal HA, Chakraborti D, Majumder P, Roy P, Roy A, Bhattacharya SG, Das S. Allergenicity assessment of Allium sativum leaf agglutinin, a potential candidate protein for developing sap sucking insect resistant food crops. PLoS One. 2011;6(11):e27716.
Prasad V, Mishra SK, Srivastava S, Srivastava A. A virus inhibitory protein isolated from Cyamopsis tetragonoloba (L.) Taub. Upon induction of systemic antiviral resistance shares partial amino acid sequence homology with a lectin. Plant Cell Rep. 2014;33(9):1467–78.
Wang Z, Zhang K, Sun X, Tang K, Zhang J. Enhancement of resistance to aphids by introducing the snowdrop lectin gene gna into maize plants. J Biosci. 2005;30(5):627–38.
Yarasi B, Sadumpati V, Immanni CP, Vudem DR, Khareedu VR. Transgenic rice expressing Allium sativum leaf agglutinin (ASAL) exhibits high-level resistance against major sap-sucking pests. BMC Plant Biol. 2008;8:102.
Zhao W, Lu LX, Yang PC, Cui N, Kang L, Cui F. Organ-specific transcriptome response of the small brown planthopper toward rice stripe virus. Insect Biochem Molec. 2016;70:60–72.
Zhao W, Yang PC, Kang L, Cui F. Different pathogenicities of Rice stripe virus from the insect vector and from viruliferous plants. New Phytol. 2016;210(1):196–207.
Zhu JJ, Jiang F, Wang XH, Yang PC, Bao YY, Zhao W, Wang W, Lu H, Wang QS, Cui N, et al. Genome sequence of the small brown planthopper, Laodelphax striatellus. Gigascience. 2017;6(12):1–12.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet Journal. 2011;17:10–2.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
McCarthy DJ, Chen YS, Smyth GK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012;40(10):4288–97.
Wu HJ, Ma YK, Chen T, Wang M, Wang XJ. PsRobot: a web-based plant small RNA meta-analysis toolbox. Nucleic Acids Res. 2012;40(W1):W22–8.
John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS. Human MicroRNA targets. PLoS Biol. 2004;2(11):1862–79.
Cox MP, Peterson DA, Biggs PJ. SolexaQA: at-a-glance quality assessment of Illumina second-generation sequencing data. Bmc Bioinformatics. 2010;11:485.
Kim D, Landmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–U121.
Anders S, Pyl PT, Huber W. HTSeq--a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9.
Alexa AaR, J. topGO: Enrichment Analysis for Gene Ontology. R package version 2301 2016.
This work was funded by the National Key Plan for Scientific Research and Development of China (2017YFD0200406, 2016YFC1200603, 2017YFD0200904), Strategic Priority Research Program of the Chinese Academy of Sciences (CAS) (XDB11040000, XDB11050000), Ministry of Agriculture transgenic major projects (2016ZX08010001), Natural Science Foundation of China (NSFC 91540116, 31772162, 31472051, 31471782), and the Open Research Fund Program of State Key Laboratory of Integrated Management of Pest Insects and Rodents (Chinese IPM 1801, Chinese IPM 1605).
Availability of data and materials
The mRNA and sRNA sequencing data for both mock-inoculated and RSV-infected rice and small brown planthopper samples have been deposited to the Gene Expression Omnibus (GEO) database at the National Center for Biotechnology Information (NCBI) under accession number GSE113555. The genome of the RSV isolate kept in our lab was re-sequenced and and the assembled genome was deposited to the NCBI GenBank with accession numbers MF287953-MF287956.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Summary of sRNA sequencing data in mock-inoculated and RSV-infected rice samples. (XLSX 10 kb)
Table S2. Summary of sRNA sequencing data in mock-inoculated and RSV-infected planthopper samples. (XLSX 13 kb)
Figure S1. Size distribution of sRNAs in the mock-inoculated rice samples (A) and the second replicate of RSV-infected rice sample (B). Numbers of sRNAs with the same size but different 5′ terminal nucleotides were drawn in different colors. (XLS 28 kb)
Table S3. GC content of sRNAs in planthopper and rice samples. (PDF 220 kb)
Figure S2. Size distribution of sRNAs in the mock-inoculated planthopper samples (A) and other two replicates of RSV-infected planthopper samples (B). Numbers of sRNAs with the same size but different 5′ terminal nucleotides were drawn in different colors. (PDF 298 kb)
Figure S3. 5′ terminal nucleotide frequency of host-derived sRNAs in the mock-inoculated rice samples (A), the RSV-infected rice samples (B), the mock-inoculated planthopper samples (C), and the RSV-infected planthopper samples (D). (PDF 790 kb)
Figure S4. Size distribution of vsiRNAs in the second replicate of RSV-infected rice sample (A) and other two replicates of RSV-infected planthopper samples (B). Numbers of vsiRNAs with the same size but different 5′ terminal nucleotides were drawn in different colors. (PDF 306 kb)
Figure S5. Comparison of size distributions of vsiRNAs in RSV-infected rice and planthopper samples. (PDF 107 kb)
Figure S6. 5′ terminal nucleotide frequency of vsiRNAs in two RSV-infected rice samples (A) and three RSV-infected planthopper samples (B). (PDF 378 kb)
Figure S7. Percentage of vsiRNAs that were derived from the RSV genomic (+) and antigenomic (−) sequences in RSV-infected rice and small brown planthopper samples. (PDF 78 kb)
Figure S8. Number of up- and down-regulated genes in RSV-infected rice and planthopper samples compared with the corresponding mock-inoculated samples. (PDF 102 kb)
Figure S9. Volcano plot to show the fold change and error rates of differentially expressed genes (DEGs) in RSV-infected rice. Green dots represent DEGs. (PDF 118 kb)
Figure S10. Gene ontology (GO) enrichment analysis of the up-regulated genes in RSV-infected rice. (PDF 207 kb)
Figure S11. Volcano plot to show the fold change and error rates of differentially expressed genes (DEGs) in small brown planthopper after infected with RSV for 5 days. Green dots represent DEGs. (PDF 111 kb)
Figure S12. Gene ontology (GO) enrichment analysis for the up-regulated genes in small brown planthopper after infected with RSV for 5 days. (PDF 306 kb)
Table S4. List of 17 down-regulated rice genes that are potentially targeted by rice-specific vsiRNAs. (XLS 34 kb)
Table S5. List of 6 down-regulated planthopper genes that are potentially targeted by planthopper-specific vsiRNAs. (XLSX 54 kb)
Table S6. Primers used for RT-PCR analysis of vsiRNAs in RSV-infected rice and RSV-infected planthopper. (XLS 34 kb)
About this article
Cite this article
Yang, M., Xu, Z., Zhao, W. et al. Rice stripe virus-derived siRNAs play different regulatory roles in rice and in the insect vector Laodelphax striatellus. BMC Plant Biol 18, 219 (2018). https://doi.org/10.1186/s12870-018-1438-7