Figures
Abstract
RNA silencing, mediated by small RNAs including microRNAs (miRNAs) and small interfering RNAs (siRNAs), is a potent antiviral or antibacterial mechanism, besides regulating normal cellular gene expression critical for development and physiology. To gain insights into host small RNA metabolism under infections by different viruses, we used Solexa/Illumina deep sequencing to characterize the small RNA profiles of rice plants infected by two distinct viruses, Rice dwarf virus (RDV, dsRNA virus) and Rice stripe virus (RSV, a negative sense and ambisense RNA virus), respectively, as compared with those from non-infected plants. Our analyses showed that RSV infection enhanced the accumulation of some rice miRNA*s, but not their corresponding miRNAs, as well as accumulation of phased siRNAs from a particular precursor. Furthermore, RSV infection also induced the expression of novel miRNAs in a phased pattern from several conserved miRNA precursors. In comparison, no such changes in host small RNA expression was observed in RDV-infected rice plants. Significantly RSV infection elevated the expression levels of selective OsDCLs and OsAGOs, whereas RDV infection only affected the expression of certain OsRDRs. Our results provide a comparative analysis, via deep sequencing, of changes in the small RNA profiles and in the genes of RNA silencing machinery induced by different viruses in a natural and economically important crop host plant. They uncover new mechanisms and complexity of virus-host interactions that may have important implications for further studies on the evolution of cellular small RNA biogenesis that impact pathogen infection, pathogenesis, as well as organismal development.
Author Summary
Small RNA-mediated RNA silencing is a widespread antiviral or antibacterial mechanism in different organisms. Although the host and pathogen factors involved in this mode of host defense and pathogen counter-defense have been extensively investigated, much less is known about how a pathogen alters the small RNA metabolism in a host. To help fill this knowledge gap, we first used deep sequencing to characterize the miRNA and siRNA profiles of rice plants infected by two distinct viruses, Rice dwarf virus (RDV) and Rice stripe virus (RSV), respectively. Our analyses showed that these two viruses had distinct impacts on rice small RNA metabolism. More significantly, RSV infection, but not RDV infection, enhanced the accumulation of some rice miRNA*s from conserved miRNA precursors and accumulation of phased siRNAs from a particular precursor. Furthermore, RSV infection also induced the expression of novel phased miRNAs from several conserved miRNA precursors. While RSV infection significantly elevated the expression of certain OsDCLs and OsAGOs, RDV infection only affected the expression of certain OsRDRs. These data uncover new mechanisms of virus-host interactions that affect host small RNA metabolism.
Citation: Du P, Wu J, Zhang J, Zhao S, Zheng H, Gao G, et al. (2011) Viral Infection Induces Expression of Novel Phased MicroRNAs from Conserved Cellular MicroRNA Precursors. PLoS Pathog 7(8): e1002176. https://doi.org/10.1371/journal.ppat.1002176
Editor: Raul Andino, University of California San Francisco, United States of America
Received: October 31, 2010; Accepted: June 7, 2011; Published: August 25, 2011
Copyright: © 2011 Du et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work is supported by grants from the Natural Science Foundation of China (31030005 and 30910103904) and Transgenic Program of China (2009ZX08009-059B, 2008ZX08010-001, 2008ZX08009-003). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
RNA-mediated gene silencing is a widespread mechanism of host defense against viral [1]–[6] and bacterial [7] infections. The 21–24 nt small RNAs, produced from DICER processing of double-stranded RNAs (dsRNAs) or RNA transcripts with stem-loop structures, are broadly defined as small interfering RNAs (siRNAs) and microRNAs (miRNAs), respectively [8]–[10]. They are incorporated into ARGONAUTES (AGOs) to form an RNA-INDUCED SILENCING COMPLEX (RISC). The RISC then recognizes its target RNA/DNA sequences through specific base pairing, to activate RNA cleavage or translation repression or DNA methylation [9], [11]–[14]. In plants, the miRNA precursors are processed into miRNA/miRNA* duplexes mostly by DICER-LIKE 1 (DCL1) with 2-nt 3′overhangs [9], [15], [16]. After methylation at the 3′ end, the miRNA sequences are preferentially incorporated into RISC to regulate gene expression, whereas the miRNA* sequences are usually degraded [17], [18].
Viruses encode dedicated proteins that function as viral suppressors of RNA silencing (VSRs), or other multi-functional proteins, to defend against host RNA silencing by interfering with distinct steps of the silencing pathways [19]. Previous studies, based on RNA gel blots, showed that transgenic expression of VSRs from many different plant viruses often caused reduced accumulation of conserved miRNAs [20]–[26]. The Tobacco mosaic virus movement and coat protein interactions also alter accumulation of tobacco miRNAs [27]. The biological function of this down-regulated miRNA accumulation for viral infection or plant defense remains to be understood.
The siRNAs are produced via processing of dsRNAs derived from distinct sources and are classified into three types: trans-acting siRNAs (ta-siRNAs), natural antisense transcript-derived siRNAs (nat-siRNAs) and repeat-associated siRNAs (ra-siRNAs). The ta-siRNAs are generated in a phased pattern through DCL4-processing of dsRNA substrates formed via the activity of RNA-DEPENDENT RNA POLYMERASE 6 (RDR6) [28]–[31]. The nat-sRNAs are produced from dsRNAs formed by natural antisense cis-transcript pairs [32], [33]. The ra-siRNAs are derived from transposons and repetitive elements [34], [35]. Transgenic expression of VSRs from some plant viruses also often led to reduced accumulation of siRNAs, likely as a means of dampening host silencing against viral infection [20], [23], [24].
There is evidence that some cellular miRNAs play anti-viral roles against animal viruses, although a particular miRNA is exploited to support viral infection [4], [5]. Many animal viruses cause generally down-regulated expression of host miRNAs as shown by microarrays [36]–[39] and quantitative real-time PCR [40]. Deep sequencing also revealed a similar pattern, and additionally identified a few new miRNAs induced only in virus-infected cells [41]–[43]. Some new miRNAs are induced in an organ-specific manner [41].
To gain further insights into viral interactions with the host RNA silencing pathways, we used deep sequencing to characterize the small RNA profiles of rice plants infected by two different RNA viruses together with microarray and quantitative RT-PCR analyses of the expression patterns of RNA silencing pathway genes. Rice is one of the most important crop plants and emerging models for RNA silencing studies [44]–[49]. Rice dwarf virus (RDV), which causes millions of dollars crop losses each year, is a member of Phytoreovirus whose genome consists of 12 dsRNAs that encode 12 proteins. The RDV non-structural protein Pns10 has been identified as a VSR, which has siRNA-duplex binding activities [50], [51]. Rice stripe virus (RSV), another devastating rice pathogen, is a member of Tenuivirus whose genome consists of four negative sense and ambisense single-stranded RNAs that encode seven proteins. The nonstructural protein NS3 functions as a VSR that also has siRNA-duplex binding activities [52]. Both viruses are transmitted via insect vector in a persistent manner and the eggs from viruliferous female adults also carry viruses and spread diseases. RSV and RDV are transmitted to rice plants by the small brown plant hopper (Laodelphax striatellus) and leafhopper (Nephotettix cincticeps or Resilia dorsalis), respectively.
Our analyses showed that RSV and RDV infections differentially affected rice small RNA profiles. RSV infection enhanced the accumulation of some rice miRNA*s, but not their corresponding miRNAs, as well as accumulation of phased siRNAs from a particular precursor. Significantly, RSV infection also induced the expression of novel miRNAs in a phased pattern from several conserved miRNA precursors. In comparison, no such changes in host small RNA expression was observed in RDV-infected rice plants. RSV infection significantly elevated the expression levels of selective OsDCLs and OsAGOs, whereas RDV infection only affected the expression of certain OsRDRs. Our results provide a comparative analysis, via deep sequencing, of changes in the small RNA profiles and in the genes of RNA silencing machinery induced by different viruses in a natural and economically important crop host plant. They uncover new mechanisms and complexity of virus-host interactions that may have important implications for further studies on the evolution of cellular small RNA biogenesis that impact pathogen infection, pathogenesis organismal development, and crop-protection technology development.
Results
Deep sequencing of small RNAs
We sequenced the small RNAs from rice plants infected with RDV and RSV, respectively, and from plants mock-inoculated controls by using the Solexa/Illumina deep sequencing method. The three-week-old rice seedlings were exposed, respectively, to viruliferous leafhopper, N. cincticeps (RDV), virus-free N. cincticeps (RDV, mock), viruliferous planthopper, L. Striatellus (RSV) and virus-free L. Striatellus (RSV, mock). After three weeks, when virus-induced symptoms appeared in the systemic leaves of virus-infected plants, the plants from each treatment (i.e., RDV-infected, RDV mock, RSV-infected and RSV mock) were pooled to prepare an RNA library for sequencing. We performed sequencing with three biological replicates for each treatment, with approximately 15 plants pooled in each replicate. (All the sequencing data can be available from the website: http://www.cbi.pku.edu.cn/download/rdsv/rdsv.html.) From each library, we obtained more than 50% of small RNA sequences that had at least one perfect match in the rice genome and no more than one mismatch in the virus genome (Table 1). The similar percentages of small RNA sequences matched to the rice and virus genomes in all replicates for a given treatment indicate similar quality of RNA preparation and sequencing. Because of the large volume of total sequence data and our primary focus in this study on learning about how the two viruses would affect small RNA profiles and the RNA silencing machinery in a common host, the virus-derived small RNA profiles and their biological implications will be analyzed and presented in a future report. For the rice small RNAs, we normalized the total sequence reads of each library to one million, and then used the average read value of unique sequences from all replicates in each treatment for further analysis.
Infection by RSV, but not by RDV, led to altered expression of selective rice miRNAs
We compared the total reads of rice miRNAs among the four libraries, using pooled data from the three independent biological replicates. There were 4060, 4852, 8449 and 5014 unique sequences, with 89036, 91409, 118975 and 90214 reads, that match perfectly to miRNA precursors from RDV-infected, mock (RDV)-inoculated, RSV-infected, and mock (RSV)-inoculated rice plants, respectively (Table 2). Notably, the number of unique sequences mapped to miRNA precursors from the RSV-infected rice plants was nearly twice of those from each of the other three types of plants.
Of all the sequence reads mapped to miRNA precursors, 94.7% (84344), 93.1% (85112), 68.9% (81975) and 91.0% (82086) were mature miRNA sequences, and 1.6% (1437), 2.0% (1813), 27.4% (32600) and 1.9% (1758) were miRNA* sequences (Table 2). Interestingly, the reads of miRNA* sequences in the RSV-infected rice small RNA library were 15 times higher than those from the other three libraries, and was nearly half of the miRNA reads from the same library. The sequence data from the three biological replicates are presented in Supplemental Table S1, which reproducibly show more than about 10-times higher accumulation of miRNA* sequences in RSV-infected rice plants than in the other three groups of plants.
For many miRNAs there were no obvious differences in expression levels between RSV-infected and mock-inoculated rice plants. However, a number of miRNAs showed significant changes in expression in RSV-infected plants. These were defined as those having reads of 100 or more and showing at least two-fold changes in reads between RSV-infected and mock-inoculated plants. As shown in Table 3, expressions of miR156b, miR159a1, miR164a, miR166, miR167a, miR1884b, miR393b, miR396e and miR528 were down regulated, whereas miR535, miR390 and miR171 were up regulated in the infected plants. Although data from the three independent biological replicates were reproducible (Supplemental Table S1 and S2), we nonetheless used RNA gel blots to verify the altered expressions of some miRNAs with higher reads. As shown in Figure 1A, miR156, miR166 and miR167 were down regulated, whereas miR172 showed no obvious changes in accumulation. The expression of miR168 showed a two-fold increase on RNA gel blots, in near agreement with the sequence data (Table 3). By contrast, fewer miRNAs showed changes in accumulation levels in RDV-infected plant as compared to those in the mock-inoculated plants. Specifically, miR167a, miR171 and miR1863 were down regulated and only miR393 was induced in RDV-infected rice plants (Table 3). We confirmed the down-regulation of miR167 in RDV-infected rice in northern blots (Figure 1A). Thus, different viruses can differentially affect the expression of some miRNAs in a common host.
(A) RNA gel blots showing expression of miR156, miR164, miR166, miR167, miR168 and miR172 in virus infected rice plants. Rice U6 was reprobed as a control. The expression levels from mock (RDV)-inoculated plants are set to a value of 1.00 and the other three plants are expressed relative to this reference value. (B) Expression analysis of miRNA targets from RDV-infected plants and (C) from RSV-infected plants by quantitative real-time RT-PCR analysis. The expression levels of the assayed genes were normalized to the expression level of OsEF-1α.Os12g41680, Os03g43930, Os04g57610 and Os03g60430 are the targets of miR164, miR166, miR167 and miR172, respectively [54].
Using quantitative real-time RT-PCR, we analyzed the expression of one target gene for each miRNA shown in Figure 1A. As shown in Figure 1B, the expression of AGO1b, a target of miR168, was induced during RSV infection (Figure 1C). This mimics the situation in Arabidopsis thaliana where the expressions of miR168 and AGO1 are transcriptionally co-regulated [53]. Os03g43930, a HD-Zip transcription factor and a target of miR166 [54], was up-regulated in agreement with the down-regulation of miR166 expression during RSV infection (Figure 1C). The expression of Os03g60430 (target of miR172) [55] showed no changes, whereas the expression of Os04g57610 (target of miR167) [54] increased in RSV infected rice (Figure 1C). These patterns correlated well with unaltered expression miR172 and down-regulated expression of miR167. Notably, none of these genes showed significant changes in expression levels, as did their cognitive miRNAs, in RDV-infected plants (Figure 1B). Using four different pairs of primers, we were unable to obtain conclusive data about the expression of potential miR156 targets (LOC_Os04g46580 and LOC_Os07g 32170).
RSV infection enhanced accumulation of specific rice miRNA*
In RSV-infected rice plants, many miRNA*s accumulated to high levels, whereas their corresponding miRNA sequences did not show any obvious changes compared to the mock control (Table 3). These include miRNA*s for some members of the miR160 family (miR160a–f), miR166 family (miR166a–e, g–l and n), miR167 family (miR167a, c–e, h and i), miR171 family (miR171c–f and miR171i), miR396 family (miR396a–c, e and f) and miRNA* of miR1318, miR1425, miR159a, miR168, miR172d, miR390, miR444b.2 and miR528.
These data were not attributed to sequencing bias. First, the vast majority of miRNA* in RSV-infected plants were present at low levels as in mock control plants. Second, the levels of both miRNA and miRNA* in RDV-infected plants were not different from those in the mock controls. Third, northern blots confirmed that miR1425*, miR160* and miR171* were significantly up regulated, whereas their corresponding miRNAs were not (Figure 2A). Thus, we concluded that RSV infection specifically enhanced the accumulation of certain miRNA* sequences, but not their corresponding miRNAs. We notice shorter sequences for miR160* and miR171* in RSV-infected rice (Figure 2A), which maybe 20-nt variants of miRNA* through sequencing data (Table 3; see also Table S2 for all sequencing data from three biological replicates).
(A) Expression analysis of miR1425/miR1425*, miR171/miR171* and miR160/miR160* in RDV infected, mock (RDV)-inoculated, RSV-infected as well as mock (RSV)-inoculated rice plants by RNA gel blots. Rice U6 was reprobed as a control. (B) RNA gel blot analysis of the expression of Os11g15060, a predicted target of miR1425*. The 28S rRNA stained with ethidium bromide was used as a loading control. (C) Mapping of the cleavage site in Os11g15060 by 5′-RACE. The numbers above the arrows indicate the frequencies of sequenced RACE clones corresponding to each inferred cleavage site.
We predicted the putative targets for some of these miRNAs*. Os11g15060, a SAM (S-adenosyl-L-Met)-dependent carboxyl methyltransferase, was predicted to be a target of miR1425*. Northern blots showed down-regulated expression of Os11g15060 in RSV-infected plants (Figure 2B). A cleavage product was detected (Figure 2B), and 5′-RACE (Rapid Amplification of 5′ Complementary DNA Ends) identified the cleavage site in the miR1425* binding region (Figure 2C), providing direct evidence that a plant miRNA* could direct the cleavage of its target mRNA. We also identified a second cleavage site outside the miR1425* binding region (Figure 2C). How this second site was derived is not clear, but many miRNA targets have more than one cleavage site as reported in other plants [56]–[58] and even in a green alga [59]. Whether this second cleavage results from the action of another yet-to-be identified small RNA induced by viral infection remains to be further investigated. Intriguingly, the expression of Os11g15060 was also down-regulated in RDV-infected plants (Figure 2B). However, absence of a cleavage product and failed 5′ RACE to identify a cleavage site (data not shown) suggests that this down-regulation in the RDV-infected plants was not caused by RNA silencing, or caused by partial silencing as well as another mechanism that remains to be identified. We also analyzed the expression of some potential targets of miR160* and miR171* including Os11g38140 and Os02g49240 (potential targets of miR160*), and Os03g38170 (potential target of miR171*). The expression of the three genes decreased in RSV-infected rice as compared with that in RDV-infected rice (Figure S1).
RSV infection induced expression of new phased miRNAs from conserved precursors
The small RNA libraries of RSV-infected rice contained many unique sequences, absent from the other three libraries, which are mapped to several conserved miRNA precursors (Table 2). Some of these belong to the miR159 family, whose precursors are approximately 200 nt in length with a stem structure of 80–90 nt (Figure 3). As shown in Figure 3A–D (see also Supplemental Table S3 for all sequence data), besides the reported miRNA/miRNA* pair for each precursor of the family (i.e., miR159a.1/miR159a.1*, miR159a.2/miR159a.2*, miR159b.1/miR159b.1*, miR159c.1/miR159c.1* and miR159f.1/miR159f.1*, labeled red and blue respectively for each pair), new pairs of miRNA/miRNA* were produced from each precursor. These included miR159a.3/miR159a.3*, miR159b.2/miR159b.2*, miR159b.3/miR159b.3*, miR159c.2/miR159c.2*, miR159c.3/miR159c.3*, miR159f.2/miR159f.2* and miR159f.3/miR159f.3*. These pairs were generated in a phased pattern and often detected at higher levels than the reported pairs. The phased miRNA-miRNA*s have the characteristics of 2 nt overhangs at the 3′ end. These observations, together with the observation that such phased production of new miRNA/miRNA* were absent from RDV-infected or any mock-infected plants, ruled out the possibility that they were degradation products or sequencing errors.
(A–E) Precursors of phased miRNAs. The sequences shown in red and blue are reported miRNAs and miRNA*s, respectively. The histograms indicate the reads of corresponding sequences on the precursors shown in black bars. (F) RNA gel blot showing the expression of phased miRNAs in RSV-infected rice. Rice U6 was reprobed as a control. (G) Expression analysis of the predicted targets of phased miRNAs by real time RT-PCR. Os03g02240 was previously validated as a target of miR159.2. [61]. Os12g03530 and Os02g4700 are predicted targets of miR159a.3 and miR394.2, respectively. The expression levels of the assayed genes were normalized to the expression level of OsEF-1α.
In addition to the above three-duplex phase forms of miRNA/miRNA*s, we also found two-duplex phase forms derived from some other precursors. The precursor of miR394 in the miRBase database (http://microrna.sanger.ac.uk/sequences, version 12.0) is about 100 nt in length. However, we found that the actual precursor is longer and contains a 27-nt extension from the 5′ and 3′ ends of the reported precursor, respectively (Figure 3E and Figure S2). In this longer precursor, a new miR394.2/394.2* duplex appeared at the distal end of the stem structure, in phase with the reported 394.1/394.1* (Figure 3E and Figure S2).
Northern blots confirmed the expression of miR159.2/miR159.2*, miR159.3 and miR394.2 (Figure 3F). The resolution of northern blots did not permit distinction between family members, so that each band could contain multiple members of a family of miRNAs/miRNA*s.
Using clustalW [60], we analyzed the conservation of precursor sequences of miR159, miR319, and miR394 in different plants (Figure S3). We found that, compared with the reported mature miRNA/miRNA* sequences, the newly identified phased miRNA/miRNA* sequences are less well conserved.
The target of miR159a.2 is Os03g02240, a GT2 transcription factor, which has been verified by 5′ RACE [61]. Quantitative real-time RT-PCR analysis showed that RSV infection down-regulated the expression of Os03g02240 (Figure 3G). Os12g03530 and Os02g47000, the predicted targets of miR159a.3 and miR394.2, respectively, were down-regulated in RSV-infected rice plants (Figure 3G). However, in RDV-infected rice plants, the expression of these genes was unchanged (data not shown).
RSV infection enhanced the accumulation of distinct rice phased siRNA from a pre-existing RNA precursor
The transcripts encoded by the AK120922 locus of rice genome can fold into long inverted repeat structures, producing 12 21-nt phased small RNA duplexes [45], [62], [63]. From the stem region proximal to the terminal loop to the distal end of the RNA secondary structure, the 12 small RNA duplexes produced are named P1-12_5′ on the 5′ arm and P1-12_3′ on the 3′ arm (Figure 4A). One of these duplexes was initially characterized as miR436/miR436* duplex [62], but further studies demonstrated all small RNAs, including the so-called miR436, are DCL4-dependent siRNAs [45]. Analysis of our sequencing data showed that the reads of some AK120922-derived siRNAs in the RSV-infected rice plants were much higher than those in the mock-inoculated plants. In particular, the reads of P5_3′ increased by at least 100-fold. The higher expressions of these siRNAs were confirmed by northern blots (Figure 4B). No such changes were observed in RDV-infected plants (data not shown). Using realtime PCR, we analyzed the expression of the potential targets of P5_3′ and P10_3′. Surprisingly, we found that both of them were up-regulated in RSV-infected plants (Figure S4).
(A) Accumulation of siRNAs in RSV-infected rice compared that in mock-inoculated rice. The upper structure illustrates the long RNA hairpin. The read numbers below the hairpin, which were each normalized to one million with total reads of the corresponding library, indicate the expression levels of siRNAs at corresponding positions in the 3′ arm of the hairpin. (B) RNA gel blots confirm accumulation of the highly induced P10_3′ and P5_3′, hybridized with the complementary 32P-labeled probes.
RDV and RSV infections differentially modified the expression of rice RNA silencing pathway genes
To gain additional insights into the effects of RDV and RSV infection on the RNA silencing pathways/machinery, we characterized the expression profiles of various genes involved in the biogenesis/function of miRNAs/siRNAs by using the same plant materials as used for small RNA deep sequencing. These genes include 8 homologs of OsDCLs, five OsRDRs and 19 AGOs that have been identified in rice [44]–[49].
We first analyzed the expression profiles of OsDCLs, OsRDRs and OsAGOs, according to the annotations of Kapoor et al [48]. Figure 5A presents microarray data showing that RDV and RSV infections affected the expression of different members of OsDCL, OsRDR and OsAGO families. The microarray data were further verified by quantitative real-time RT-PCR measurements (Figure 5B and C). Of the 8 DCL homologs in the rice genome [48], OsDCL3a (LOC_Os01g68120) and OsDCL3b (LOC_Os10g34430) were significantly down-regulated, and OsDCL2 (LOC_Os03g38740) significantly up-regulated, in RSV-infected plants (Figure 5C). In contrast, RDV infection had almost no influence on the expression of OsDCLs in rice plants (Figure 5B).
(A) Expression profiles from microarrays. The bar on the bottom represents the scale of relative expression levels of genes (Log2). (B) Expression profiles from RDV-infected plants and (C) Expression profiles from RSV-infected plant based on quantitative real-time RT-PCR analysis. The expression levels of the assayed genes were normalized to the expression level of OsEF-1α. The expression level in the corresponding mock rice was set to 1.0. The nomenclature of these genes is based on Kapoor et al. (2008).
In rice, there are four AGO1 homologs: OsAGO1a (LOC_Os02g45070), OsAGO1b (LOC_Os04g47870), OsAGO1c (LOC_Os02g58490) and OsAGO1d (LOC_Os06g51310) [32], [33]. The expression of OsAGO1a-c as well as OsAGO2 (LOC_Os04g52540) and OsAGO18 (LOC_Os07g28850) increased significantly in RSV-infected rice plants (Figure 5C). In RDV-infected plants, only the expression of OsAGO11 increased significantly (Figure 5B).
The expression of OsRDR1 (LOC_Os02g50330) and OsRDR4 (LOC_Os01g10140) was significantly enhanced in RDV-infected plants (Figure 5B), but not in RSV-infected plants (Figure 5C). OsRDR2 did not change expression levels in plants infected by either virus.
These data demonstrate that the two viruses have distinct effects on the expression of different genes of the RNA silencing pathways in the common host rice.
Discussion
Recent studies used deep sequencing to obtain profiles of viral siRNAs [64]–[68] and viroid small RNAs [69], [70] in infected plants. Such data laid a foundation for further investigations on the biogenesis mechanisms and functions of viral and viroid small RNAs. Our present work provides the first deep sequencing analysis of plant small RNA profiles under viral infection conditions. This analysis not only filled a critical knowledge gap in RNA silencing-based virus-host interactions, but provided novel insights into the impact of viral infection on host small RNA biogenesis. Our results showed down-regulated as well as up-regulated accumulation of certain rice miRNAs, with up-regulation being more extensive. Most significantly, RSV infection, but not RDV infection, induced the expression of novel phased miRNAs from several families of conserved cellular miRNA precursors. Real-time RT-PCR experiments showed reduced accumulations of the predicted target mRNAs for these phased miRNAs (i.e., miR159.2 families, miR159.3 and miR394.2), indicating that the induced phased miRNAs have regulatory functions. How the reduced accumulation of these target mRNAs contributes to the establishment of viral infection and/or host defense will be an important focus of future research.
As compared to the short animal miRNA precursors, which are usually 70–80 nt, plant miRNA precursors are generally much longer with most known precursors to be 200–300 nt [15]. These long plant precursors make it possible to produce multiple miRNAs. Indeed, multiple, phased miRNAs are produced from some miRNA precursors in the single-cell alga, Chlamydomonas reinhardtii [59], [71] and in higher plants [72]. Some of these miRNAs are differentially regulated by bacterial infection in A. thaliana [72]. Our finding that new phased miRNAs are induced during the infection of a plant virus significantly broadens the landscape of phased miRNA biogenesis during pathogen infection. Further analyses of the host small RNA profiles involving different pathogens and hosts may lead to additional examples, and an understanding of the broad biological significance, of pathogen infection-induced expression of phased and other novel miRNAs.
The enhanced production of miRNA*s during RSV infection may be attributed to inhibition of RISC activity and miRNA/miRNA* unwinding by the RSV VSR, via direct or indirect interaction with the miRNA/miRNA* duplexes, as have been shown for VSRs of other viruses [22], [24]. If this were the case, one would expect elevated accumulations of miRNAs and miRNA*s at the same time. However, our deep sequencing showed that only the accumulations of miRNA*s, not the corresponding miRNAs from many miRNA precursors were enhanced. Such special accumulation of miRNA*, but not miRNAs, may be explained if RSV enhanced the activities of certain RISCs specially associated with miRNA*s or interfered with loading of some miRNAs into RISCs. Since we do not have data yet to support this postulation, we need to leave other alternative possibilities open for exploration. Schnettler et al. (2010) reported stabilization of miR171c/miR171c* complex by several tospoviruses in infected N. benthamiana, with elevated accumulation of miR171c* in the infected plants. This stabilization appears to be due to the activity of the viral NSs proteins [73].
The specific accumulation of miRNA*s of several miRNA families and the conservation of the miRNA* sequences within the families suggest that there may be certain base preferences among the miRNA*s accumulated during RSV infection. Using WebLogo [74], we found that there was an ‘A’ bias in the 19th nucleotide from the 5′ terminus (Figure S5). Considering that the A. thaliana AGO2 preferentially produces miRNA* sequences [75] and the up-regulation of OsAGO2 and OsAGO18 during RSV infection of rice, we propose that the ‘A’ bias may direct certain miRNA* sequences into these OsAGOs. Still we cannot exclude the potential influence of OsAGO1, as the ‘A’ bias in the 19th nucleotide from the 5′ terminus of miRNA* corresponds to the 5′ terminal ‘U’ of miRNA. In Drosophila melanogaster, miRNA*s were reported to have regulatory functions [76]. Here, by northern blots and 5′ RACE, we demonstrated that Os11g15060, a predicted target of miR1425*, was specifically down regulated by cleavage during RSV infection (Figure 2C). This is the first direct experimental demonstration, to the best of our knowledge, that a plant miRNA* can regulate the stability of its predicted target mRNA.
Although many siRNAs were produced in rice [77], our current analysis showed that RSV infection enhanced mainly the accumulation of some phased siRNAs derived from the AK120922 transcripts. The specificity of this enhancement was supported by the observation that RDV infection did not have such an effect. Previous studies showed that accumulation of siRNAs decreased in transgenic plants expressing the VSRs of several plant viruses [20], [23], [24]. The inhibitory effect of VSRs on siRNA accumulation, via binding of VSRs with siRNAs, is one of the mechanisms of viral counter-defense [1], [2], [19]. The enhanced accumulation of phased siRNAs from AK120922 during RSV infection, in contrast to the generally observed decrease in siRNA accumulations, may be biologically significant. Whether the VSR or other proteins encoded by RSV play a role in this enhancement is an outstanding mechanistic question to be answered in the future.
The induced expression of new phased miRNAs and the enhanced production of selective miRNA*s and phased siRNAs in rice plants infected by RSV, but not by RDV, indicate that such changes were not a general response to viral infection. Rather, they were caused by distinct virus-host interactions. One of the primary consequences of such interactions was conceivably altered expression and/or function of the RNA silencing machinery, which further leads to altered small RNA profiles. Indeed, our microarray and quantitative real-time RT-PCR analyses demonstrated that the expression files of OsDCLs, OsAGOs and OsRDRs were differentially altered in rice infected with the two viruses. During RDV infection, with the exception of OsRDR1 being enhanced, there were no significant changes in the expression of RNA silencing pathway genes. In contrast, during RSV infection, the expression levels of OsDCL3a and OsDCL3b were down-regulated, whereas the expression of OsDCL2 was enhanced. Three of the four OsAGO1 homologs, OsAGO1a, OsAGO1b and OsAGO1c, were up-regulated in RSV-infected plants. Whether the altered expression patterns of the OsDCLs, OsRDRs and OsAGOs are responsible for the altered miRNA/siRNA biogenesis/accumulations in the RSV- and RDV-infected rice plants remains to be determined. Current data indicate that OsDCL1 participates in the production of 21-nt miRNAs and DCL4 mainly produces some siRNAs [28], [29]. DCL3 was recently reported to produce 24-nt miRNAs [14]. The four small RNA libraries we generated comprised mainly 24-nt sequences. However, in the RSV-infected rice small RNA libraries, the number of 24-nt small RNAs was reduced to about three quarters of the other three libraries (Figure S6). This correlated with a decreased level of DCL3 expression in RSV-infected rice. These results suggested that OsDCL3 may be primarily responsible for producing 24 nt small RNAs, just as the A. thaliana homolog does (Figure 5) [78], [79]. The production of 21 nt and 24 nt phased miRNAs suggests involvement of DCL1 and potentially also DCL3.
In summary, our data indicate that at least some viruses may have evolved mechanisms to induce expression of phased miRNAs from well-conserved cellular miRNA precursors. Whether such new miRNAs play a role in host defense or in viral infection will be an important question to be investigated in the future. While this manuscript was under review, Hu et al. (2011) reported that Oilseed rape mosaic tobamovirus infection of A. thaliana also led to elevated accumulation of host siRNAs and some miRNA-like small RNAs (ml-sRNAs). These ml-sRNAs are derived in phase with known miRNAs from miRNA-precursors [80]. They are different from the sequences we report here. Altogether, as for many previous discoveries, virus-induced production of phased miRNAs and ml-sRNAs may provide yet another useful model system to study the molecular mechanisms underlying the evolution and biogenesis of miRNAs. It remains to be seen whether virus-induced biogenesis of phased miRNAs is more widespread in plants and other organisms. Another important question is how different viruses affect the production of different small RNAs and how such differences impact specific mechanisms of viral infection and host responses.
Materials and Methods
Plant growth, virus infection and small RNA sequencing
RDV Fujian isolate and RSV Jiangsu isolate, China, were maintained in “Oryza sativa spp. japonica” rice plants grown in greenhouses at 25±3°C, 55±5% RH and under natural sunlight. Insects (Nephotettix cincticeps and Laodelphax striatellus) were maintained in cages that contained rice seedlings in greenhouses at 25±3°C. To obtain high viruliferous insects, nymphs were reared on virus-infected rice plants for 1 week, and insects were maintained up to the adult stage with occasional replacement of seedlings by healthy rice seedlings. Rice seedlings were grown in a greenhouse at 25–28°C and 70±5% relative humidity under natural sunlight. Three-week-old seedlings were placed individually in single tubes of 4-cm in diameter and 25-cm in height that each contained 15–20 ml of liquid nutrient medium at the bottom. The viruliferous insects of N. cincticeps (carrying RDV), L. Striatellus (carrying RSV) as well as virus-free N. cincticeps (mock for RDV) and L. Striatellus (mock for RSV) were added to each tube. Each tube was then sealed with a nylon net. After 2 days in growth chambers with a 14-h/10-h light/dark cycle, 70±5% relative humidity and a temperature regime of 28°C (day)/25°C (night), the insects were removed and the rice seedlings were returned to the greenhouse to grow under the above greenhouse conditions. After approximately three weeks of growth, when the newly developed leaves started to show viral symptoms, the whole seedlings were harvested. With each sample, at least 15 rice seedlings were pooled for RNA extractions. Total RNAs were extracted using Trizol reagent (Invitrogen, Carlsbad, CA, USA) for RT-PCR, small RNA sequencing, microarray analysis, and northern blotting. RT-PCR was used to test individual rice seedlings for infection with RDV or RSV. All PCR primers are listed in Table S4. The small RNA library construction and Illumina 1G sequencing were carried out at BGI-Shenzhen (Shenzhen, Guangdong, China) using standard Solexa/Illumina protocols. Briefly, the total RNA was separated through 15% TBE urea denaturing PAGE gels and small RNAs of 15–30 nt were recovered. Then, 5′ and 3′ RNA adaptors were ligated to these small RNAs followed by reverse transcription into cDNAs. These cDNAs were finally amplified by PCR and subjected to Solexa/Illumina sequencing.
Bioinformatics
After removing the adaptor sequences, small RNA sequences with 18–28 nt in length were used for further analysis. BOAT provided by CBI (http://boat.cbi.pku.edu.cn/) was used for mapping small RNAs to the O. sativa genome sequences (TIGR Rice Annotation Release 5.0, ftp://ftp.plantbiology.msu.edu/pub/data/Eukaryotic_Projects/o_sativa/annotation_dbs/) as well as RDV (ftp://ftp.ncbi.nih.gov/genomes/Viruses/Rice_dwarf_virus_uid14797/) and RSV (ftp://ftp.ncbi.nih.gov/genomes/Viruses/Rice_stripe_virus_uid14795/) Rice stripe virus) genome sequences. Small RNAs with perfect genomic matches were used for further analysis. The small RNAs were annotated with reference to the following databases: miRBase (http://microrna.sanger.ac.uk/sequences, version 12.0) for miRNA sequences, Rfam (http://www.sanger.ac.uk/Software/Rfam/) for noncoding RNA (rRNAs, tRNAs, snoRNAs, and snRNAs) sequences and Repbase (http://www.girinst.org) for transposons and repeats.
WebLogo [74] was used for analyzing of relative frequencies of nucleotides at each position of small the sRNAs, and the mfold program [81] was used for predicting the stem-loop structures. Target genes were predicted by using the miRU web server (http://bioinfo3.noble.org/miRNA/miRU.htm) and our own Perl script [59]. For submission to the web server, we chose the default parameters (score for each 20 nt≤3, G: U pairs≤6, indels≤1 and other mismatches ≤3) and the TIGR Rice Genome mRNA dataset (OSA1 release 5, 01/23/2007). For our own Perl script, based on the score standard mentioned before [59], targets with score ≤5 were chosen. Identical search results from both methods were considered potential targets of newly identified miRNAs or phased miRNAs. For miRNA*s induced during RSV infection, their targets were validated based on the microarray data. ClustalW was used for the alignment of miRNA precursors and miRNA* sequences.
Microarray hybridization and data analysis
Hybridization of GeneChip rice genome array (Affymetrix), scanning and analysis were performed by the Affymetrix custom service (CapitalBio, Beijing, China) following standard protocols (http://www.affymetrix.com/support/technical/manual/expression_manual.affx). Three biological replicates were conducted for RNA from each type of plant samples. Data analysis and comparison of the samples was finished using the standard Affymetrix protocol (CapitalBio). Cluster3.0 software was used for cluster analysis. The expression profile of some mentioned gene through microarray analysis was shown in Table S5.
Northern blot hybridization
Total RNA was used for mRNA and small RNA northern blot hybridization. For small RNA blots, 10–60 µg of total RNA (depending on the relative expression levels from deep sequencing) of each sample was separated on 15% polyacrylamide denaturing gels and then transferred to Hybond-N+ membranes (Amersham BioScience, Piscataway, NJ, USA). The membranes were cross-linked by UV transillumination and dried at 120°C for 30 min. DNA oligonucleotides complementary to small RNAs, which were labeled with γ-32P-ATP by T4 polynucleotide kinase (New England Biolabs, Beverly, MA, USA), were used as probes for hybridization. Membranes were prehybridized with buffer for 2 h followed by hybridization with the DNA probes overnight at 40°C in 5X SSC, 20 mM Na2HPO4 (PH 7.2), 7% SDS, 2X Denhardt's Solution, 100 µg/ml salmon sperm DNA. After washing twice at 40°C with 2X SSC and 0.1% SDS, the blots were imaged with a PhosphorImager (PerkinElmer, Shelton, CT, USA). The membranes were stripped with 0.1X SSC and 1% SDS at 100°C and reprobed. Probes complementary to U6 sequences were used as a loading control.
For mRNA blots, 15 µg of total RNA was separated by 1% formaldehyde agarose gel and transferred to Hybond-N+ membranes that were then cross-linked and dried as described above. Prehybridization and hybridization solution was 5X SSC, 1% SDS, 5X Denhardt's Solution, 100 µg/ml salmon sperm DNA and 50% formamide. The probes were labeled with α-32P-dCTP by DNA polymerase 1 large (Klenow) fragment (Promega, Madison, Wisconsin, USA) and were complementary to a 500-bp fragment corresponding to the 5′ partial sequence of Os11g15060 (target of miR1425*).
Quantitative real-time RT-PCR
Total RNA was treated with RNase-free DNase I (TAKARA Biotechnology, Dalian, China) at 37°C for 30 min. After phenol/chloroform extraction and isopropanol precipitation, the RNA was quantified with a UV/ visible spectrophotometer (Amersham). Two µ\g of total RNA was reverse-transcribed using poly (T) adapter with SuperScript Reverse Transcriptase (Invitrogen). qPCR was performed using SYBR Green Realtime PCR Master Mix (Toyobo) OsEF-1a gene was used as an internal control, with primers CX1597 (59-GCACGCTCTTCTTGCTTTCACTCT-39) and CX1598 (59-AAAGGTCACCACCATACCAGGCTT-39) [45]. Three independent biological replicates were conducted. These data were further normalized with the normalized expression values obtained from qRT-PCRs and bar charts plotted by using Microsoft Excel and SPSS (Statistical Product and Service Solutions) software (IBM; Version.10.0). All the other primers used are listed in Supplemental Table S4.
Target gene validation by 5′ RACE
Validation of target genes by mapping the cleavage sites was conducted with 5′ RACE by following the GeneRacer Kit (Invitrogen) protocols described previously [57]. Total RNA of RSV-infected rice was directly ligated to GeneRacer RNA Oligo adaptor without any modifications. RT-PCR was used to synthesize cDNAs using GeneRacer Oligo dT primer. GeneRacer 5′ Primer (5′CGACTGGAGCACGAGGACACTGA3′) and the target gene-specific outer primers (Table S4) were used for the first round of amplification. Then the GeneRacer 5′ Nested Primer (5′ GGACACTGACATGGACTGAAGGAGTA) and gene-specific inner primers (Table S4) were used for the second round of amplification. The PCR products were cloned and sequenced to identify the cleavage sites.
Supporting Information
Figure S1.
Realtime RT-PCR tested the expression of targets of miR160* and miR171* during RDV infection (A) and RSV infection (B). The expression levels of the assayed genes were normalized to the expression level of OsEF-1α. Os11g38140 and Os02g49240 were the potential targets of miR160*, and Os03g38170 was a potential target of miR171*.
https://doi.org/10.1371/journal.ppat.1002176.s001
(DOC)
Figure S2.
Structure of miR394 precursor and small RNAs produced from miR394 precursor (A), The secondary structure of newly-identified miR394 precursor (165nt), which was done using Mfold. (B), Small RNAs generated from the new-identified miR394 precursors. The reported miRNA and miRNA* were shown in red and blue, and the new phased miRNA and miRNA* were shown in purple and indigo, respectively. The numbers following the small RNA sequences were the reads of the corresponding sequences in RSV-infected and mock (RSV) rice small RNA libraries of the first biological repeat.
https://doi.org/10.1371/journal.ppat.1002176.s002
(PDF)
Figure S3.
Alignment of genomic sequences of phased miRNA precursors in different plants. The miRNA precursor sequences used for the alignment come from the miRBase database (http://microrna.sanger.ac.uk/sequences, version 12.0). Osa: Oryza sativa; sof: Saccharum officinarum; sbi: Sorghum bicolor; zma: Zea mays; ptc: Populus trichocarpa; ath: Arabidopsis thaliana; gma: Glycine max; vvi: Vitis vinifera; sly: Solanum lycopersicum; mtr: Medicago truncatula.
https://doi.org/10.1371/journal.ppat.1002176.s003
(PDF)
Figure S4.
Realtime RT-PCR tested the expression of targets of P5_3′ and P10_3′ during RDV infection (A) and RSV infection (B). Os03g52630 was the potential targets of P5_3′, and Os04g51070 was a potential target of P10_3′. The expression levels of the assayed genes were normalized to the expression level of OsEF-1α.
https://doi.org/10.1371/journal.ppat.1002176.s004
(DOC)
Figure S5.
The nucleotide bias of RSV induced miRNA* sequences calculated based on the unique sequences (A) and reads (B).
https://doi.org/10.1371/journal.ppat.1002176.s005
(DOC)
Figure S6.
Size distribution of rice small RNAs in virus-infected and mock-inoculated rice. Proportion of unique sequences of different sizes in the total rice unique sequences of the four libraries.
https://doi.org/10.1371/journal.ppat.1002176.s006
(DOC)
Table S1.
Summary of total small RNA reads mapped to known rice miRNA precursors in virus-infected and mock-inoculated rice plants. The footnotes of the table are as follows: a Reads were normalized to one million with the total reads of each library. b Perfect match to sense miRNA precursor sequences from the miRBase database (http://microrna.sanger.ac.uk/sequences, version 12.0). c Encompasses the defined miRNA/miRNA* sequence ±2 nt on each side. d Red indicates the up-regulated small RNAs in RSV infected rice plants, but not in the other plants.
https://doi.org/10.1371/journal.ppat.1002176.s007
(DOC)
Table S2.
Sequencing reads of known miRNAs and miRNA*s in four libraries from three biological replicates. The footnotes of the table are as follows: a Reads were normalized to one million with the total sequence reads of each library. b The sequences are potential miRNA*s, some of which were not detected during our sequencing. c The miR396f* is also reported as miR396f-3p (miRBase database: http://www.mirbase.org/). d Red and green indicate miRNAs whose reads (> = 50 for miRNA in at least one sample) increased (red) or decreased (green) by at least two-fold in comparison with those from mock-inoculated plants in at least two biological replicates. They also indicate miRNA* sequences whose reads (> = 10 in at least one sample) increased (red) or decreased (green) by at least two-fold in virus-infected plants in comparison with those from mock-infected plants in at least two biological replicates. e These sequences have 20-nt variants which have high reads in RSV-infected rice libraries.
https://doi.org/10.1371/journal.ppat.1002176.s008
(DOC)
Table S3.
Newly-identified miRNAs from the known miRNA precursors and their target genes. The footnotes of the table are as follows: a Reads was average values of three-repeat reads of each libraries, which normallized to 1 million with the total sequencing reads. b The precursor of miR394 is longer than the one in miRBase database, and miR394.2 have higher reads than miR394 only in RSV-infected rice library. c From the same position,other varient sequences with different length was produced. d The newly-identified miRNAs only have high reads in RSV-infected rice library, and form the phased-miRNA with the reported miRNA before . e The newly-identified miRNA are in tandem configuration with the reported miRNAs. f The sequences have much higher reads than the reported ones. g The miRNA sequences, reads number of which were outstanding shown in red and blue, were greatly accumulated and reduced only during RSV infection, respectively.
https://doi.org/10.1371/journal.ppat.1002176.s009
(DOC)
Table S4.
Primers and probes used in this paper. * indicates the corresponding miRNA star sequences.
https://doi.org/10.1371/journal.ppat.1002176.s010
(DOC)
Table S5.
Expressional analysis of RNA silencing pathway genes by microarray. The footnote of the table is as follows: a Fold change = (normalized signal intensity of virus-infected rice samples)/ (normalized signal intensity of mock-inoculated rice samples).
https://doi.org/10.1371/journal.ppat.1002176.s011
(DOC)
Acknowledgments
We thank Professors Gongyin Ye (Zhejiang University, China) and Yijun Zhou (Institute of Plant Protection, Jiangsu Academy of Agricultural Sciences, China) for providing Nephotettix cincticeps and Laodelphax Striatellus, respectively, and for helping with virus inoculations. We also thank Drs. Yijun Qi, Zicai Liang and Xiaoming Zhang for critical reading of the manuscript.
Author Contributions
Conceived and designed the experiments: PD YL. Performed the experiments: PD JW. Analyzed the data: PD JZ HZ GG LW YL. Contributed reagents/materials/analysis tools: SZ. Wrote the paper: PD YL.
References
- 1. Li F, Ding SW (2006) Virus counterdefense: Diverse strategies for evading the RNA-silencing immunity. Annu Rev Microbiol 60: 503–531.
- 2. Ding SW, Voinnet O (2007) Antiviral immunity directed by small RNAs. Cell 130: 413–426.
- 3. Mlotshwa S, Pruss GJ, Vance V (2008) Small RNAs in viral infection and host defense. Trends Plant Sci 13: 375–382.
- 4. Umbach JL, Cullen BR (2009) The role of RNAi and microRNAs in animal virus replication and antiviral immunity. Genes Dev 23: 1151–1164.
- 5. Cullen BR (2010) Five Questions about Viruses and MicroRNAs. PLoS Pathog 6: e1000787.
- 6. Ding SW (2010) RNA-based antiviral immunity. Nat Rev Immunol 10: 632–644.
- 7. Padmanabhan C, Zhang XM, Jin HL (2009) Host small RNAs are big contributors to plant innate immunity. Curr Opin Plant Biol12: 465–472.
- 8. Ramachandran V, Chen XM (2008) Degradation of microRNAs by a family of exoribonucleases in Arabidopsis. Science 321: 1490–1492.
- 9. Voinnet O (2009) Origin, Biogenesis, and Activity of Plant MicroRNAs. Cell 136: 669–687.
- 10. Malone CD, Hannon GJ (2009) Small RNAs as Guardians of the Genome. Cell 136: 656–668.
- 11. Baulcombe D (2004) RNA silencing in plants. Nature 431: 356–363.
- 12. Brodersen P, Sakvarelidze-Achard L, Bruun-Rasmussen M, Dunoyer P, Yamamoto YY, et al. (2008) Widespread translational inhibition by plant miRNAs and siRNAs. Science 320: 1185–1190.
- 13. Zamore PD, Haley B (2005) Ribo-gnome: The big world of small RNAs. Science 309: 1519–1524.
- 14. Wu L, Zhou HY, Zhang QQ, Zhang JG, Ni FR, et al. (2010) DNA Methylation Mediated by a MicroRNA Pathway. Mol Cell 38: 465–475.
- 15. Bartel DP (2004) MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 116: 281–297.
- 16. Kim VN, Nam JW (2006) Genomics of microRNA. Trends Genet 22: 165–173.
- 17. Baumberger N, Baulcombe DC (2005) Arabidopsis ARGONAUTE1 is an RNA Slicer that selectively recruits rnicroRNAs and short interfering RNAs. Proc Natl Acad Sci U S A 102: 11928–11933.
- 18. Yu B, Yang ZY, Li JJ, Minakhina S, Yang MC, et al. (2005) Methylation as a crucial step in plant microRNA biogenesis. Science 307: 932–935.
- 19. Diaz-Pendon JA, Ding SW (2008) Direct and indirect roles of viral suppressors of RNA silencing in pathogenesis. Annu Rev Phytopathol 46: 303–326.
- 20. Mallory AC, Reinhart BJ, Bartel D, Vance VB, Bowman LH (2002) A viral suppressor of RNA silencing differentially regulates the accumulation of short interfering RNAs and micro-RNAs in tobacco. Proc Natl Acad Sci U S A 99: 15228–15233.
- 21. Kasschau KD, Xie Z, Allen E, Llave C, Chapman EJ, et al. (2003) P1/HC-Pro, a viral suppressor of RNA silencing, interferes with Arabidopsis development and miRNA unction. Dev Cell 4: 205–217.
- 22. Chapman EJ, Prokhnevsky AI, Gopinath K, Dolja VV, Carrington JC (2004) Viral RNA silencing suppressors inhibit the microRNA pathway at an intermediate step (vol 18, pg 1179, 2004). Genes Dev 18: 1510–1510.
- 23. Chen J, Li WX, Xie D, Peng JR, Ding SW (2004) Viral virulence protein suppresses RNA silencing-mediated defense but upregulates the role of microrna in host gene expression. Plant Cell 16: 1302–1313.
- 24. Dunoyer P, Lecellier CH, Parizotto EA, Himber C, Voinnet O (2004) Probing the microRNA and small interfering RNA pathways with virus-encoded suppressors of RNA silencing. Plant Cell 16: 1235–1250.
- 25. Chellappan P, Vanitharani R, Fauquet CM (2005) MicroRNA-binding viral protein interferes with Arabidopsis development. Proc Natl Acad Sci U S A 102: 10381–10386.
- 26. Mlotshwa S, Schauer SE, Smith TH, Mallory AC, Herr JM Jr, et al. (2005) Ectopic DICER-LIKE1 expression in P1/HC-Pro Arabidopsis rescues phenotypic anomalies but not defects in microRNA and silencing pathways. Plant Cell 17: 2873–2885.
- 27. Bazzini AA, Hopp HE, Beachy RN, Asurmendi S (2007) Infection and coaccumulation of tobacco mosaic virus proteins alter microRNA levels, correlating with symptom and plant development. Proc Natl Acad Sci U S A 104: 12157–12162.
- 28. Peragine A, Yoshikawa M, Wu G, Albrecht HL, Poethig RS (2004) SGS3 and SGS2/SDE1/RDR6 are required for juvenile development and the production of trans-acting siRNAs in Arabidopsis. Genes Dev 18: 2368–2379.
- 29. Vazquez F, Vaucheret H, Rajagopalan R, Lepers C, Gasciolli V, et al. (2004) Endogenous trans-acting siRNAs regulate the accumulation of Arabidopsis mRNAs. Mol Cell 16: 69–79.
- 30. Yoshikawa M, Peragine A, Park MY, Poethig RS (2005) A pathway for the biogenesis of trans-acting siRNAs in Arabidopsis. Genes Dev 19: 2164–2175.
- 31. Allen E, Xie ZX, Gustafson AM, Carrington JC (2005) microRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell 121: 207–221.
- 32. Borsani O, Zhu JH, Verslues PE, Sunkar R, Zhu JK (2005) Endogenous siRNAs derived from a pair of natural cis-antisense transcripts regulate salt tolerance in Arabidopsis. Cell 123: 1279–1291.
- 33. Katiyar-Agarwal S, Gao S, Vivian-Smith A, Jin H (2007) A novel class of bacteria-induced small RNAs in Arabidopsis. Genes Dev 21: 3123–3134.
- 34. Vazquez F (2006) Arabidopsis endogenous small RNAs: highways and byways. Trends Plant Sci 11: 460–468.
- 35. Henderson IR, Jacobsen SE (2007) Epigenetic inheritance in plants. Nature 447: 418–424.
- 36. Godshalk SE, Bhaduri-McIntosh S, Slack FJ (2008) Epstein-Barr virus-mediated dysregulation of human microRNA expression. Cell Cycle 7: 3595–3600.
- 37. Wang FZ, Weber F, Croce C, Liu CG, Liao XD, et al. (2008) Human cytomegalovirus infection alters the expression of cellular microRNA species that affect its replication. J Virol 82: 9065–9074.
- 38. Li Y, Chan EY, Li J, Ni C, Peng X, et al. (2010) MicroRNA expression and virulence in pandemic influenza virus-infected mice. J Virol 84: 3023–3032.
- 39. Liu X, Wang T, Wakita T, Yang W (2010) Systematic identification of microRNA and messenger RNA profiles in hepatitis C virus-infected human hepatoma cells. Virology 398: 57–67.
- 40. Buck AH, Perot J, Chisholm MA, Kumar DS, Tuddenham L, et al. (2010) Post-transcriptional regulation of miR-27 in murine cytomegalovirus infection. RNA 16: 307–315.
- 41. Wang Y, Brahmakshatriya V, Zhu H, Lupiani B, Reddy SM, et al. (2009) Identification of differentially expressed miRNAs in chicken lung and trachea with avian influenza virus infection by a deep sequencing approach. BMC Genomics 10: 512.
- 42. Parameswaran P, Sklan E, Wilkins C, Burgon T, Samuel MA, et al. (2010) Six RNA Viruses and Forty-One Hosts: Viral Small RNAs and Modulation of Small RNA Repertoires in Vertebrate and Invertebrate Systems. PLoS Pathog 6: e1000764.
- 43. Cui L, Guo X, Qi Y, Qi X, Ge Y, et al. (2010) Identification of microRNAs involved in the host response to enterovirus 71 infection by a deep sequencing approach. J Biomed Biotechnol 425939.
- 44. Liu B, Li PC, Li X, Liu CY, Cao SY, et al. (2005) Loss of function of OsDCL1 affects microRNA accumulation and causes developmental defects in rice. Plant Physiol 139: 296–305.
- 45. Liu B, Chen ZY, Song XW, Liu CY, Cui X, et al. (2007) Oryza sativa dicer-like4 reveals a key role for small interfering RNA silencing in plant development. Plant Cell 19: 2705–2718.
- 46. Nagasaki H, Itoh JI, Hayashi K, Hibara KI, Satoh-Nagasawa N, et al. (2007) The small interfering RNA production pathway is required for shoot meristern initiation in rice. Proc Natl Acad Sci U S A 104: 14867–14871.
- 47. Nonomura KI, Morohoshi A, Nakano M, Eiguchi M, Miyao A, et al. (2007) A germ cell-specific gene of the ARGONAUTE family is essential for the progression of premeiotic mitosis and meiosis during sporogenesis in rice. Plant Cell 19: 2583–2594.
- 48. Kapoor M, Arora R, Lama T, Nijhawan A, Khurana JP, et al. (2008) Genome-wide identification, organization and phylogenetic analysis of Dicer-like, Argonaute and RNA-dependent RNA Polymerase gene families and their expression analysis during reproductive development and stress in rice. BMC Genomics 9: 1471–2164.
- 49. Wu L, Zhang QQ, Zhou HY, Ni FR, Wu XY, et al. (2009) Rice MicroRNA Effector Complexes and Targets. Plant Cell 21: 3421–3435.
- 50. Cao XS, Zhou P, Zhang XM, Zhu SF, Zhong XH, et al. (2005) Identification of an RNA silencing suppressor from a plant double-stranded RNA virus. J Virol 79: 13018–13027.
- 51. Ren B, Guo YY, Gao F, Zhou P, et al. (2010) Multiple Functions of Rice Dwarf Phytoreovirus Pns10 in Suppressing Systemic RNA Silencing. J Virol 84: 12914–12923.
- 52. Xiong RY, Wu JX, Zhou YJ, Zhou XP (2009) Characterization and subcellular localization of an RNA silencing suppressor encoded by Rice stripe tenuivirus. Virology 387: 29–40.
- 53. Vaucheret H, Mallory AC, Bartel DP (2006) AGO1 homeostasis entails coexpression of MIR168 and AGO1 and preferential stabilization of miR168 by AGO1. Mol Cell 22: 129–136.
- 54. Liu Q, Zhang YC, Wang CY, Luo YC, Huang QJ, et al. (2009) Expression analysis of phytohormone-regulated microRNAs in rice, implying their regulation roles in plant hormone signaling. FEBS Lett 583: 723–728.
- 55. Zhou M, Gu L, L P, Song X, Wei L, et al. (2010) Degradome sequencing reveals endogenous small RNA targets in rice (Oryza sativa L . ssp . indica). Front Biol 5: 67–90.
- 56. Jones-Rhoades MW, Bartel DP (2004) Computational identification of plant MicroRNAs and their targets, including a stress-induced miRNA. Mol Cell 14: 787–799.
- 57. Llave C, Xie ZX, Kasschau KD, Carrington JC (2002) Cleavage of Scarecrow-like mRNA targets directed by a class of Arabidopsis miRNA. Science 297: 2053–2056.
- 58. Allen E, Xie ZX, Gustafson AM, Sung GH, Spatafora JW, et al. (2004) Evolution of microRNA genes by inverted duplication of target gene sequences in Arabidopsis thaliana. Nat Genet 36: 1282–1290.
- 59. Zhao T, Li GL, Mi SJ, Li S, Hannon GJ, et al. (2007) A complex system of small RNAs in the unicellular green alga Chlamydomonas reinhardtii. Genes Dev 21: 1190–1203.
- 60. Thompson JD, Higgins DG, Gibson TJ (1994) Clustal-W - Improving the Sensitivity of Progressive Multiple Sequence Alignment through Sequence Weighting, Position-Specific Gap Penalties and Weight Matrix Choice. Nucleic Acids Res 22: 4673–4680.
- 61. Lacombe S, Nagasaki H, Santi C, Duval D, Piegu B, et al. (2008) Identification of precursor transcripts for 6 novel miRNAs expands the diversity on the genomic organisation and expression of miRNA genes in rice. BMC Plant Bio 8: 123.
- 62. Sunkar R, Girke T, Jain PK, Zhu JK (2005) Cloning and characterization of MicroRNAs from rice. Plant Cell 17: 1397–1411.
- 63. Sunkar R, Girke T, Zhu JK (2005) Identification and characterization of endogenous small interfering RNAs from rice. Nucleic Acids Res 33: 4443–4454.
- 64. Qi XP, Bao FS, Xie ZX (2009) Small RNA Deep Sequencing Reveals Role for Arabidopsis thaliana RNA-Dependent RNA Polymerases in Viral siRNA Biogenesis. PloS One 4: e4971.
- 65. Wang XB, Wu QF, Ito T, Cillo F, Li WX, et al. (2010) RNAi-mediated viral immunity requires amplification of virus-derived siRNAs in Arabidopsis thaliana. Proc Natl Acad Sci U S A 107: 484–489.
- 66. Szittya G, Moxon S, Pantaleo V, Toth G, Rusholme Pilcher RL, et al. (2010) Structural and functional analysis of viral siRNAs. PLoS Pathog 6: e1000838.
- 67. Garcia-Ruiz H, Takeda A, Chapman EJ, Sullivan CM, Fahlgren N, et al. (2010) Arabidopsis RNA-dependent RNA polymerases and dicer-like proteins in antiviral defense and small interfering RNA biogenesis during Turnip Mosaic Virus infection. Plant Cell 22: 481–496.
- 68. Lin KY, Cheng CP, Chang BC, Wang WC, Huang YW, et al. (2010) Global analyses of small interfering RNAs derived from Bamboo mosaic virus and its associated satellite RNAs in different plants. PloS One 5: e11928.
- 69. Navarro B, Pantaleo V, Gisel A, Moxon S, Dalmay T, et al. (2009) Deep sequencing of viroid-derived small RNAs from grapevine provides new insights on the role of RNA silencing in plant-viroid interaction. PloS One 4: e7686.
- 70. Bolduc F, Hoareau C, St-Pierre P, Perreault JP (2010) In-depth sequencing of the siRNAs associated with peach latent mosaic viroid infection. BMC Mol Biol 11: 16.
- 71. Molnar A, Schwach F, Studholme DJ, Thuenemann EC, Baulcombe DC (2007) miRNAs control gene expression in the single-cell alga Chlamydomonas reinhardtii. Nature 447: 1126–1129.
- 72. Zhang W, Gao S, Zhou X, Xia J, Chellappan P, et al. (2010) Multiple distinct small RNAs originate from the same microRNA precursors. Genome Biol 11: R81.
- 73. Schnettler E, Hemmes H, Huismann R, Goldbach R, Prins M, et al. (2010) Diverging affinity of tospovirus RNA silencing suppressor proteins, NSs, for various RNA duplex molecules. J Virol 84: 11542–54.
- 74. Crooks GE, Hon G, Chandonia JM, Brenner SE (2004) WebLogo: A sequence logo generator. Genome Res 14: 1188–1190.
- 75. Mi SJ, Cai T, Hu YG, Chen Y, Hodges E, et al. (2008) Sorting of small RNAs into Arabidopsis argonaute complexes is directed by the 5′ terminal nucleotide. Cell 133: 116–127.
- 76. Okamura K, Phillips MD, Tyler DM, Duan H, Chou YT, et al. (2008) The regulatory activity of microRNA star species has substantial influence on microRNA and 3′ UTR evolution. Nat Struct Mol Biol 15: 354–363.
- 77. Nobuta K, Venu RC, Lu C, Belo A, Vemaraju K, et al. (2007) An expression atlas of rice mRNAs and small RNAs. Nat Biotechnol 25: 473–477.
- 78. Blevins T, Rajeswaran R, Shivaprasad PV, Beknazariants D, Si-Ammour A, et al. (2006) Four plant Dicers mediate viral small RNA biogenesis and DNA virus induced silencing. Nucleic Acids Res 34: 6233–6246.
- 79. Deleris A, Gallego-Bartolome J, Bao JS, Kasschau KD, Carrington JC, et al. (2006) Hierarchical action and inhibition of plant Dicer-like proteins in antiviral defense. Science 313: 68–71.
- 80. Hu Q, Hollunder J, Niehl A, Kørner CJ, Gereige D, et al. (2006) Specific impact of tobamovirus infection on the Arabidopsis small RNA profile. PLoS One 6: e19549.
- 81. Zuker M (2003) Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res 31: 3406–3415.