Skip to main content

Construction of small RNA-mediated gene regulatory networks in the roots of rice (Oryza sativa)

Abstract

Background

The root systems play essential roles for plants to anchorage to the soil, and to exploit the mineral and water resources. The molecular mechanisms underlying root development have been extensively studied to improve root system architecture, especially for the crops. Several microRNA (miRNA) families have been demonstrated to be involved in plant root development. However, whether the other small RNA (sRNA) species, which occupy a dominant portion of the plant endogenous sRNA population, possess potential roles in root development remains unclear.

Results

In this study, by using sRNA high-throughput sequencing data, we made a comparison of the sRNA accumulation levels between the rice root tips and the whole roots. The sRNAs highly accumulated in the root tips and in the whole roots were extracted respectively. After Argonaute 1 (AGO1) enrichment analysis, the sRNAs with great potential of performing target cleavages were included for target prediction and degradome sequencing data-based validation. As a result, lists of the targets regulated by the AGO1-enriched sRNAs were obtained for both the root tips and the whole roots. Further evidences were identified from microarray data of the target genes to support some of the sRNA—target interactions. Specifically, the expression patterns of certain target genes in the root tips and the whole roots were contrary to those of the regulating sRNAs. Besides, several targets were indicated to play important roles in root development based on literature mining.

Conclusions

Taken together, the regulatory networks mediated by the sRNAs highly accumulated in the root tips or in the whole roots could advance our current understanding of the sRNA-involved molecular mechanisms underlying rice root development. And, the sRNA—target lists could serve as the basis for further functional investigations.

Background

As one of the important organs, the root system helps a plant to anchorage to the soil, and to absorb mineral nutrients and water for normal development and growth during its life cycle. Since the root system architecture greatly affects the efficiency of nutrient and water absorption, the molecular factors influencing root development have been extensively investigated in order to facilitate genetic modification of the agronomic traits of the plants, especially for the crops [1, 2].

MicroRNAs (miRNAs), a small RNA (sRNA) species, were identified to be involved in various biological processes of plants during the past ten years [3, 4]. Certain miRNA families are indispensable for the establishment of the root system architecture [5], and are implicated in the transduction of multiple signals, such as nutrition, hormone and stress [6]. For examples, miR160 is required for root cap formation by targeting Auxin Response Factor (ARF) 10 and 16 in Arabidopsis[7]. miR164-mediated regulation of NAM/ATAF/CUC 1 (NAC1) is involved in Arabidopsis lateral root development through downregulating auxin signals [8]. Both miR160 and miR167 play important roles in adventitious rooting in Arabidopsis[9].

Based on high-throughput sequencing (HTS), a much more huge sRNA population was uncovered in plants in addition to the miRNAs. Whether some of these sRNAs play potential roles in root development remains unclear, although a few experimental hints point to this possibility recently [10]. To partially address this issue, we took advantage of the public sRNA HTS data originated from the root tips and the whole roots of rice (Oryza sativa). We made a comparison of the sRNA accumulation levels between the two biological samples to identify the sRNAs highly accumulated in the root tips and in the whole roots respectively. Then, these sRNAs were subjected to Argonaute 1 (AGO1) enrichment analysis. And, those enriched in AGO1 were considered to have potential of performing target cleavages, and were included in target prediction and degradome sequencing data-based validation. Based on the target lists, regulatory networks mediated by the AGO1-enriched sRNAs were established for both the root tips and the whole roots. Notably, within the networks, the expression patterns of several target genes, based on the microarray data, were contrary to those of the regulating sRNAs, supporting the interactions between the sRNAs and the targets. Besides, literature mining indicated that certain target genes played potential roles in root development. These evidences indicate the reliability of the networks. Our results could advance the current understanding of the regulatory roles of the sRNAs, in addition to the miRNAs, in plant root growth and development.

Results and discussion

Identification of the sRNAs highly accumulated in the root tips or in the whole roots of rice

We took advantage of the public sRNA HTS data of rice. The data sets could be divided into two groups: the ones (GSM409313, GSM409314 and GSM409315) prepared from the rice root tips (~250 μm) and the ones (RCn1I and RCn2D) from the whole roots (Figure 1). According to the previous reports [11, 12], all of these data sets originated from ~14-day-old rice (Oryza sativa L. ssp. japonica) seedlings, thus making the two data groups more comparable. To extract the sRNAs highly accumulated in the root tips, the following rules were applied: (1) the sRNA should be present in at least one of the “root tip” data sets (GSM409313, GSM409314 and GSM409315); (2) its normalized accumulation levels should be 3 RPM (reads per million; the normalization step is similar to that for the degradome sequencing data; see Methods for detail) or higher, and should be three times or more than the levels of the same sequences in all the “whole root” data sets (RCn1I and RCn2D). Similar rules were applied to identify the sRNAs highly accumulated in the whole roots. As a result, 91,265 sRNAs including eight miRNA sequences (osa-miR169i-5p.2, osa-miR390-3p, osa-miR531b, osa-miR1430, osa-miR1847.2, osa-miR1865-5p, osa-miR1874-3p and osa-miR5082) showed high abundances in the root tips (Additional file 1), and 7,446 sRNAs including 24 miRNA sequences (osa-miR156l-5p, osa-miR164d, osa-miR164e, osa-miR166k-3p/osa-miR166l-3p, osa-miR169f.2, osa-miR171f-5p, osa-miR393a, osa-miR408-5p, osa-miR435, osa-miR528-5p, osa-miR529b, osa-miR812i/osa-miR812j/osa-miR812g/osa-miR812h, osa-miR812r, osa-miR812s, osa-miR1318-3p, osa-miR1320-5p, osa-miR1425-3p, osa-miR1427, osa-miR1432, osa-miR1861h/osa-miR1861j, osa-miR3979-3p, osa-miR5788, osa-miR5795 and osa-miR5802) were highly accumulated in the whole roots (Additional file 2). Since the “whole root” includes the “root tip” cells, we reasoned that the high levels of the 7,446 sRNAs in the whole roots were largely attributed to their high abundances within the elongation and the maturation zones of the rice roots including the main roots, the adventitious roots and the lateral roots.

Figure 1
figure 1

Workflow for the identification of the rice small RNAs (sRNAs) highly accumulated in the root tips and the whole roots based on high-throughput sequencing data. The total numbers of the identified sRNAs including microRNAs are shown in the corresponding boxes. The numbers in the brackets indicate the numbers of the microRNA sequences identified in this study. The rules for the two-step filtering are also shown.

In plants, some of the sRNAs, such as the miRNAs, exert their regulatory activities through target cleavages [3, 4, 13]. The sRNAs incorporate into specific AGO-associated RNA-induced silencing complexes (RISCs), and guide the complexes to the target transcripts to perform cleavages. It has been reported that AGO1 is an RNA slicer selectively recruiting miRNAs and short interfering RNAs in plants [14, 15]. In this consideration, we set out to identify AGO1-enriched sRNAs from the sRNAs highly accumulated in the root tips and the whole roots obtained above. The following rules were applied: (1) the sRNAs should be detectable in at least one of the AGO1-associated sRNA HTS data sets (GSM455962, GSM455963 and GSM455964); (2) its normalized accumulation levels should be 3 RPM or higher, and should be three times or more than the levels of the same sequence in the control set (GSM455965). As a result, 68 AGO1-enriched sRNAs were extracted from the sRNAs highly accumulated in the root tips, and 102 AGO1-enriched sRNAs including osa-miR528-5p were identified from the sRNAs highly accumulated in the whole roots (Figure 1, Additional file 3 and Additional file 4).

Next, the sequence characteristics of the sRNAs before and after the AGO1 enrichment filtering were examined separately. For both the sRNAs highly accumulated in the root tips and those highly abundant in the whole roots, no significant bias of the 5′ terminal nucleotide compositions was observed before the AGO1 enrichment filtering (Figure 2A). However, for those enriched in AGO1 complexes, the proportion of 5′ U (uridine)-started sRNAs was increased. Specifically, for the sRNAs highly accumulated in the whole roots, the percentage of 5′ U-started sRNAs increased from 28.68% to 39.60%. For those highly accumulated in the root tips, the percentage increased from 23.36% to 55.88% (Figure 2B). Different from the miRNAs, 21 nt (nucleotide) is not the dominant length for the other sRNAs highly accumulated in the whole roots. More interestingly, 24 nt is the dominant length for both the sRNAs and the miRNAs highly accumulated in the root tips (Figure 2C). After the AGO1 enrichment filtering, the sRNAs ranging from 19 nt to 21 nt occupied a dominant portion (96.04%) of the sRNAs abundant in the whole roots. And, for the sRNAs highly accumulated in the root tips, the 21 nt ones increased significantly (from 2.49% to 54.41%) (Figure 2D). Considering the previous report on the sequence characteristics of AGO1-associated sRNAs [16], our AGO1 enrichment filtering tended to selectively extract 5′ U-started, 19 to 21 nt sRNAs from the sRNAs highly accumulated in the root tips and the whole roots, respectively.

Figure 2
figure 2

Sequence characteristics of the small RNAs (sRNAs) and the microRNAs (miRNAs) highly accumulated in the root tips and the whole roots of rice. (A) 5 terminal compositions (measured by percentages) of the sRNAs and the miRNAs highly accumulated in the root tips (RT_sRNA and RT_miRNA) and the whole roots (WR_sRNA and WR_miRNA). (B) v5 terminal compositions (measured by percentages) of the Argonaute 1 (AGO1)-enriched sRNAs and miRNAs highly accumulated in the root tips (RT_sRNA_AGO1 and RT_miRNA_AGO1) and the whole roots (WR_sRNA_AGO1 and WR_miRNA_AGO1). (C) Sequence length distributions (measured by percentages) of the sRNAs and the miRNAs highly accumulated in the root tips and the whole roots. The occupancies of the sRNAs and the miRNAs of 21 nt and 24 nt are highlighted by yellow and red boxes respectively. (D) Sequence length distributions (measured by percentages) of the AGO1-enriched sRNAs and miRNAs highly accumulated in the root tips and the whole roots.

Identification of the targets of the AGO1-enriched sRNAs and network construction

As mentioned above, the sRNAs associated with the AGO1 silencing complexes are more likely to modulate gene expression through target cleavages. In this analysis, only the AGO1-enriched sRNAs were subjected to target prediction. miRU algorithm [17, 18] with default parameters was employed for this prediction. As a result, 1,699 target binding sites were predicted for the 68 AGO1-enriched sRNAs highly expressed in the root tips, and 1,979 binding sites were predicted for the 102 AGO1-enriched sRNAs (including osa-miR528-5p) highly expressed in the whole roots.

Degradome sequencing data is a valuable resource for large-scale validation of the sRNA—target interactions [19, 20]. Four degradome sequencing data sets of rice (see detail in Methods) were recruited to do a comprehensive validation of the predicted targets according to the analytical workflow proposed previously [21]. Considering the reports that some evident cleavage signals out of the canonical regions (10th to 11th nt of the regulating sRNAs) were observed [2225], the binding sites with prominent cleavage signals resided within 8th to 12th nt of the regulating sRNAs were retained. After target plot-based screening, 22 interactions involving 16 target genes and six sRNAs highly accumulated in the root tips, and 171 interactions involving 29 genes and 27 sRNAs (including osa-miR528-5p) highly abundant in the whole roots passed the validation (Additional file 5: Figure S1 and Additional file 6: Figure S2). Among the validated sRNA—target pairs, we observed that many binding sites were supported by evident degradome signatures resided within the canonical cleavage sites (10th to 11th nt of the regulating sRNAs) (Figure 3). In addition to osa-miR528-5p (Figure 3K), several sRNAs, such as RT_High_sRNA58165, RT_High_sRNA73918, RT_High_RNA75473, WR_High_sRNA6511 and WR_High_sRNA7168, cleaved the corresponding transcripts at the canonical sites. For WR_High_sRNA6511, we observed that the sRNA possessed two binding sites on the transcripts LOC_Os01g48060.1, LOC_Os01g54990.1 and LOC_Os05g48870.1, separately (Figure 3 F, 3G and 3H). Notably, all of these sRNAs are distinct from the current plant miRNA registries in miRBase (release 19) [26]. The result indicates that the sRNAs, which remain to be classified, could guide AGO1 complexes to perform miRNA-like target cleavages in the root systems of rice.

Figure 3
figure 3

Degradome sequencing data-based validation of the targets regulated by the small RNAs (sRNAs) highly accumulated in the root tips and the whole roots of rice. For all the figure panels (A to K), the IDs of the target transcripts and the regulating sRNAs are listed on the top. The y axes measure the intensity (in RPM, reads per million) of the degradome signals, and the x axes indicate the positions of the cleavage signals on the target transcripts. The binding sites of the sRNA regulators on their target transcripts were denoted by gray horizontal lines, and the dominant cleavage signals were marked by gray arrowheads. The figure keys for the degradome signatures from four libraries are shown on the bottom right. Please note, for (F) to (H), each target transcript possesses two binding sites of WR_High_sRNA6511.

Based on the target lists obtained above, the sRNA-mediated regulatory networks were constructed in the root tips and the whole roots by using Cytoscape [27] (Additional file 7: Figure S3 and Additional file 8: Figure S4). At first glance, RT_High_sRNA58165 regulates 16 transcripts encoded by ten different genes (Additional file 7: Figure S3), pointing to the possibility that this sRNA might play an essential role in root tip formation. More complex network was observed for those sRNAs highly accumulated in the whole roots. In some cases, one target is shared by two or more sRNAs. For example, LOC_Os07g29820.1 is regulated by WR_High_sRNA7168 and WR_High_sRNA6608. LOC_Os09g23730.1 is regulated by WR_High_sRNA317, WR_High_sRNA1124, WR_High_sRNA1878, WR_High_sRNA2359 and WR_High_sRNA2909 (Additional file 8: Figure S4).

Biological meanings of the networks supported by target expression- and literature mining-based evidences

Considering the repressive role of the AGO1-enriched sRNAs in gene expression, we deduced that for the sRNA—target pairs validated above, the target genes should display expression patterns contrary to those of the sRNA regulators. For example, the target gene of a sRNA highly accumulated in the root tips was supposed to have higher expression level in the whole roots than in the root tips. To support this notion, 24 microarray data sets contributed by a previous report [28] were utilized. These data sets enabled us to investigate the expression levels of all the target genes in different parts of the roots (eight different regions: root cap, division zone, elongation zone, maturation zone I, maturation zone II, maturation zone III, maturation zone IV and maturation zone V; each with three replicates). According to the report by Takehisa et al.[28], in this analysis, the microarray data sets of the root cap and the division zone were grouped into the “root tip” category (~200 μm), which was anatomically comparable to the sRNA HTS library group (GSM409313, GSM409314 and GSM409315) prepared from the root tips (~250 μm). All the 24 microarray data sets were treated as the “whole root” category in order to make it anatomically consistent with the sRNA HTS library group (RCn1I and RCn2D) prepared from the whole roots. All the microarray data sets were generated from the 10-day-old rice seedlings, which was largely consistent with the 14-day-old seedlings for sRNA library preparation. For the expression level analysis, the mean value of the three replicates was calculated for each target gene.

The microarray data-based analysis provided supportive evidences for some of the sRNA—target interactions. Specifically, the level of LOC_Os02g43370.1 regulated by RT_High_sRNA55881 in the whole roots is 4.53 times higher than its level in the root tips (0.125 vs. 0.028) (Figure 4B and Additional file 9: Table S1). For LOC_Os07g41280.1/LOC_Os07g41280.2/LOC_Os07g41280.3 cleaved by RT_High_sRNA58165, their expression value in the whole roots is 5.547 versus 0.670 in the root tips. Also regulated by RT_High_sRNA58165, the level of LOC_Os09g27500.1 in the whole roots is 0.317 versus 0.023 in the root tips (Figure 4A and Additional file 9: Table S1). For LOC_Os04g57610.1/LOC_Os04g57610.3 regulated by WR_High_sRNA6608, their expression value in the whole roots is 1.355, which is lower than that (3.909) in the root tips. For LOC_Os06g46410.1/LOC_Os06g46410.2 also regulated by WR_High_sRNA6608, their expression value in the whole roots is 0.643 versus 1.463 in the root tips. Also for WR_High_sRNA6608, the level of its target transcript LOC_Os12g41950.1 in the whole roots is 1.736, which is lower than that (4.185) in the root tips (Figure 5B and Additional file 9: Table S1). For LOC_Os09g23730.1 regulated by WR_High_sRNA317, WR_High_sRNA1124, WR_High_sRNA1878, WR_High_sRNA2359 and WR_High_sRNA2909, its expression levels in the whole roots and the root tips are 5.430 and 8.273 respectively (Figure 5C and Additional file 9: Table S1). Taken together, the results demonstrated that at least some of the target genes displayed expression patterns contrary to those of the regulating sRNAs. Failure to validate the regulatory relationships of the remaining sRNA—target pairs by microarray data analysis might attribute to several embedded factors. The most likely reason might be that the expression of the target genes was not only modulated by the identified regulatory sRNAs, but also regulated by the other factors such as transcription factors and epigenetic factors which might disturb the contrary accumulation patterns between the sRNAs and the targets in the root systems. From this point of view, we should recognize that the regulation of the expression of a gene is quite complex, and could not be determined by only one regulator in most cases. Another potential reason is that the target gene might be specifically expressed in the root part in which the regulating sRNAs was highly abundant. For example, the sRNA was highly abundant in the root tips, and the target gene was expressed in the root tips but not in the remaining parts of the rice roots. It will lead to the failure to observe the contrary accumulation pattern between this sRNA—target pair.

Figure 4
figure 4

Certain subnetworks mediated by the small RNAs (sRNAs) highly accumulated in the root tips of rice. For both (A) and (B), the sRNAs are represented by blue nodes, and the targets by gray nodes. Notably, the target nodes in red color are the ones with expression patterns contrary to those of the regulating sRNAs, which are supported by microarray data. The target nodes with green outlines are the ones with potential biological functions in root development based on the literatures. The subnetworks were constructed by using Cytoscape [27].

Figure 5
figure 5

Certain subnetworks mediated by the small RNAs (sRNAs) highly accumulated in the whole roots of rice. For all the figure panels (A) to (D), the sRNAs are represented by blue nodes, and the targets by gray nodes. Notably, the target nodes in red color are the ones with expression patterns contrary to those of the regulating sRNAs, which are supported by microarray data. The target nodes with green outlines are the ones with potential biological functions in root development based on the literatures. The subnetworks were constructed by using Cytoscape [27].

We also intended to find some root-related evidences for the target genes based on functional annotations and literature mining. Fortunately, some positive hints were obtained. For examples, based on the annotation from TIGR rice (the rice genome annotation project established by the institute for genome research) [29], LOC_Os02g34950.1 regulated by RT_High_sRNA55881 (Figure 4B) is involved in post-embryonic development. And, its homologous gene in Arabidopsis, AT5G22370, participates in the organization of microtubules during cell division according to the TAIR (The Arabidopsis Information Resource) annotation [30]. Since the root tip includes a meristematic region in an active state of cell division, whether LOC_Os02g34950 plays an essential role in the formation and/or the maintenance of the root meristem remains to be investigated. Also in the root tips, LOC_Os11g41860 encoding an F-box domain-containing protein is regulated by RT_High_sRNA58165 based on the degradome sequencing data (Figure 4A). According to a recent report, an F-box gene LOC_Os02g44990 is involved in promoting root growth in rice [31], pointing to the possibility that LOC_Os11g41860 is functionally linked to root development. More supportive evidences were found for the target genes within the network mediated by the sRNAs abundant in the whole roots. LOC_Os03g01890, LOC_Os03g43930, LOC_Os10g33960 and LOC_Os12g41860, all of which encode START domain-containing proteins (homologous to Class III HD-ZIP proteins of Arabidopsis), were regulated by the same sRNA regulators, i.e. WR_High_sRNA5480, WR_High_sRNA6227, WR_High_sRNA6327, WR_High_sRNA6550, WR_High_sRNA6793, WR_High_sRNA6897, WR_High_sRNA6969, WR_High_sRNA7149 and WR_High_sRNA7172 (Figure 5A). In Arabidopsis and Medicago truncatula, miR166-mediated regulation of Class III HD-ZIP genes is important for root development [3234]. More interestingly, although the above nine sRNAs are not miRNAs, WR_High_sRNA6793 shares high sequence similarity with osa-miR166g-3p and osa-miR166h-3p of rice, and the other eight sRNAs share high sequence similarity with miR166 of the other plant species. In this regard, the four target genes regulated by the sRNAs homologous to miR166 might have a conserved role in rice root development. Another class of target genes including LOC_Os04g57610, LOC_Os06g46410 and LOC_Os12g41950 encode the rice ARF proteins (ARF12, ARF17 and ARF25 respectively) (Figure 5B). The three ARF genes share high sequence identity with each other, and LOC_Os04g57610 was reported to regulate root elongation in rice [35]. Moreover, the homologous ARF genes in Arabidopsis are regulated by miR167, and are essential for root system development [9, 36, 37]. Notably, the three ARF genes of rice were targeted by WR_High_sRNA6608 sharing high sequence similarity with osa-miR167. LOC_Os06g06050 targeted by osa-miR528-5p encodes OsFBL27, an F-box domain- and LRR-containing protein (Figure 5D). Its homologous gene AT2G42620 (MAX2, More Axillary Growth 2) is implicated in root development in Arabidopsis[3840]. And, miR528 was previously indicated to play a role in auxin signal-mediated rice root development [41]. LOC_Os09g23730 encoding an HMG-Y-related protein was targeted by five sRNAs highly accumulated in the whole roots (Figure 5C). Based on an early report, ectopic expression of the maize HMGB1 protein causes defects in root development of tobacco seedlings [42]. More direct evidence was found for LOC_Os07g33790 targeted by WR_High_sRNA7168 (Figure 5B) sharing sequence similarity with miR167 of the other plant species. This target gene encodes GLR (glutamate receptor) 3.4 based on the TIGR annotation. It was reported that GLR3.1 of rice was essential for the maintenance of the root apical meristem [43], and GLR3.2 and GLR3.4 of Arabidopsis were involved in lateral root initiation [44]. All these literature-based evidences support the potential roles of the target genes in rice root development.

Conclusions

Summarily, this study on the sRNAs in the rice root systems could advance our knowledge of the regulatory activities and the biological functions of the plant sRNAs other than the miRNAs. The networks established in this study were supported by the degradome sequencing data, the microarray data and the literature-based evidences, indicating that the networks could provide a basis for in-depth analysis of sRNA-mediated gene regulation in the rice roots.

Methods

Data sources

High-throughput sequencing data sets of rice sRNAs were retrieved from GEO (Gene Expression Omnibus; http://www.ncbi.nlm.nih.gov/geo/) [45] and the Rice Small RNA Database 2 belonging to the Next-Gen Sequence Databases (http://mpss.udel.edu/rice_sRNA2/index.php?menu=../common/web/library_info.php?SITE=rice_sRNA2&showAll=true) [46]. Specifically, GSM409313, GSM409314 and GSM409315 downloaded from GEO were prepared from the root tips (~250 μm) of 13 to 14-day old rice seedlings [12]. RCn1I and RCn2D downloaded from the Next-Gen Sequence Databases were prepared from the whole roots of 14-day-old rice seedlings [11].

The AGO1-associated sRNA HTS data sets were downloaded from GEO. GSM455962, GSM455963 and GSM455964 were prepared from the rice AGO1a, AGO1b and AGO1c proteins, respectively. GSM455965 was prepared from the total protein extract of rice [47].

The rice degradome sequencing data sets (GSM434596, GSM455938, GSM455939 and GSM476257) were obtained from GEO.

The microarray data is a gift from a recent report by Takehisa et al.[28]. All the data sets could be obtained from GEO under the accession ID GSE30136. Specifically, the GSM ID for all the used microarray data sets are: GSM746016 (Root cap_rep1), GSM746017 (Root cap_rep2), GSM746018 (Root cap_rep3), GSM746019 (Division zone_rep1), GSM746020 (Division zone_rep2), GSM746021 (Division zone_rep3), GSM746022 (Elongation zone_rep1), GSM746023 (Elongation zone_rep2), GSM746024 (Elongation zone_rep3), GSM746025 (Maturation zone_I_rep1), GSM746026 (Maturation zone_I_rep2), GSM746027 (Maturation zone_I_rep3), GSM746028 (Maturation zone_II_rep1), GSM746029 (Maturation zone_II_rep2), GSM746030 (Maturation zone_II_rep3), GSM746031 (Maturation zone_III_rep1), GSM746032 (Maturation zone_III_rep2), GSM746033 (Maturation zone_III_rep3), GSM746034 (Maturation zone_IV_rep1), GSM746035 (Maturation zone_IV_rep2), GSM746036 (Maturation zone_IV_rep3), GSM746037 (Maturation zone_V_rep1), GSM746038 (Maturation zone_V_rep2) and GSM746039 (Maturation zone_V_rep3). GSM746016, GSM746017, GSM746018, GSM746019, GSM746020 and GSM746021 were used as the target gene expression profiling data of the root tips [the six libraries were prepared from the root tips (200 μm) of rice, which correlates well with the sRNA libraries of the root tips], and all the 24 GSM data sets were used as a whole to investigate the expression levels of the target genes within the whole roots.

The miRNA sequences were obtained from miRBase (release 19; http://www.mirbase.org/) [26]. The transcripts of rice genes and the gene annotations were retrieved from the rice genome annotation project established by the institute for genome research (TIGR rice, release 6.1; http://rice.plantbiology.msu.edu/index.shtml) [29].

Prediction and validation of the sRNA targets

Target prediction was performed by using miRU algorithm [17, 18] with default parameters. The degradome sequencing data were utilized to validate the predicted sRNA—target pairs. First, in order to allow cross-library comparison, the normalized read count (in RPM, reads per million) of a short sequence from a specific degradome library was calculated by dividing the raw count of this sequence by the total counts of the library, and then multiplied by 106. Then, two-step filtering was performed to extract the most likely sRNA—target pairs. During the first step, the predicted sRNA binding sites along with the 50-nt surrounding sequences at both ends were collected in order to reduce the BLAST time. For the BLAST, all the collected degradome data sets were utilized at the same time to do a comprehensive search. It was based on the scenario that a sRNA—target pair was considered to be the candidate once the cleavage signal(s) existed in any data set(s). The predicted targets met the following criteria were retained: (1) there must be perfectly matched degradome signatures with their 5′ ends resided within 8th to 12th nt region from the 5′ ends of the regulating sRNAs; and (2) for a specific position within the 8th to 12th nt region, which could be regarded as the potential cleavage site, there must be two or more distinct degradome signatures with accordant 5′ ends to support this position. These retained transcripts were subjected to a second BLAST, and the degradome signals along each transcript were obtained to provide a global view of the signal noise when compared to the signal intensity within a specific target binding site. Referring to our previous study [21], both global and local t-plots were drawn. Finally, exhaustive manual filtering was performed, and only the transcripts with cleavage signals easy to be recognized were extracted as the potential sRNA—target pairs.

Availability of supporting data

The data sets supporting the results of this article are included within the article and its additional files.

Abbreviations

miRNA:

MicroRNA

sRNA:

Small RNA

AGO1:

Argonaute 1

ARF:

Auxin response factor

HTS:

High-throughput sequencing

RPM:

Reads per million

RISC:

RNA-induced silencing complex

U:

Uridine

nt:

Nucleotide

TIGR:

The Institute for Genome Research

TAIR:

The Arabidopsis information resource

MAX2:

More axillary growth 2

GLR:

Glutamate receptor

GEO:

Gene expression omnibus.

References

  1. Petricka JJ, Winter CM, Benfey PN: Control of Arabidopsis root development. Annu Rev Plant Biol. 2012, 63: 563-590. 10.1146/annurev-arplant-042811-105501.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Coudert Y, Perin C, Courtois B, Khong NG, Gantet P: Genetic control of root development in rice, the model cereal. Trends Plant Sci. 2010, 15: 219-226. 10.1016/j.tplants.2010.01.008.

    Article  CAS  PubMed  Google Scholar 

  3. Chen X: Small RNAs and their roles in plant development. Annu Rev Cell Dev Biol. 2009, 25: 21-44. 10.1146/annurev.cellbio.042308.113417.

    Article  PubMed  Google Scholar 

  4. Jones-Rhoades MW, Bartel DP, Bartel B: MicroRNAs and their regulatory roles in plants. Annu Rev Plant Biol. 2006, 57: 19-53. 10.1146/annurev.arplant.57.032905.105218.

    Article  CAS  PubMed  Google Scholar 

  5. Khan GA, Declerck M, Sorin C, Hartmann C, Crespi M, Lelandais-Briere C: MicroRNAs as regulators of root development and architecture. Plant Mol Biol. 2011, 77: 47-58. 10.1007/s11103-011-9793-x.

    Article  CAS  PubMed  Google Scholar 

  6. Meng Y, Ma X, Chen D, Wu P, Chen M: MicroRNA-mediated signaling involved in plant root development. Biochem Biophys Res Commun. 2010, 393: 345-349. 10.1016/j.bbrc.2010.01.129.

    Article  CAS  PubMed  Google Scholar 

  7. Wang JW, Wang LJ, Mao YB, Cai WJ, Xue HW, Chen XY: Control of root cap formation by MicroRNA-targeted auxin response factors in Arabidopsis. Plant Cell. 2005, 17: 2204-2216. 10.1105/tpc.105.033076.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Guo HS, Xie Q, Fei JF, Chua NH: MicroRNA directs mRNA cleavage of the transcription factor NAC1 to downregulate auxin signals for Arabidopsis lateral root development. Plant Cell. 2005, 17: 1376-1386. 10.1105/tpc.105.030841.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Gutierrez L, Bussell JD, Pacurar DI, Schwambach J, Pacurar M, Bellini C: Phenotypic plasticity of adventitious rooting in Arabidopsis is controlled by complex regulation of AUXIN RESPONSE FACTOR transcripts and microRNA abundance. Plant Cell. 2009, 21: 3119-3132. 10.1105/tpc.108.064758.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Furuta K, Lichtenberger R, Helariutta Y: The role of mobile small RNA species during root growth and development. Curr Opin Cell Biol. 2012, 24: 211-216. 10.1016/j.ceb.2011.12.005.

    Article  CAS  PubMed  Google Scholar 

  11. Jeong DH, Park S, Zhai J, Gurazada SG, De Paoli E, Meyers BC, Green PJ: Massive analysis of rice small RNAs: mechanistic implications of regulated microRNAs and variants for differential target RNA cleavage. Plant Cell. 2011, 23: 4185-4207. 10.1105/tpc.111.089045.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Johnson C, Kasprzewska A, Tennessen K, Fernandes J, Nan GL, Walbot V, Sundaresan V, Vance V, Bowman LH: Clusters and superclusters of phased small RNAs in the developing inflorescence of rice. Genome Res. 2009, 19: 1429-1440. 10.1101/gr.089854.108.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Voinnet O: Origin, biogenesis, and activity of plant microRNAs. Cell. 2009, 136: 669-687. 10.1016/j.cell.2009.01.046.

    Article  CAS  PubMed  Google Scholar 

  14. Baumberger N, Baulcombe DC: Arabidopsis ARGONAUTE1 is an RNA Slicer that selectively recruits microRNAs and short interfering RNAs. Proc Natl Acad Sci USA. 2005, 102: 11928-11933. 10.1073/pnas.0505461102.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Vaucheret H: Plant ARGONAUTES. Trends Plant Sci. 2008, 13: 350-358. 10.1016/j.tplants.2008.04.007.

    Article  CAS  PubMed  Google Scholar 

  16. Mi S, Cai T, Hu Y, Chen Y, Hodges E, Ni F, Wu L, Li S, Zhou H, Long C, Chen S, Hannon GJ, Qi Y: Sorting of small RNAs into Arabidopsis argonaute complexes is directed by the 5′ terminal nucleotide. Cell. 2008, 133: 116-127. 10.1016/j.cell.2008.02.034.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Dai X, Zhao PX: psRNATarget: a plant small RNA target analysis server. Nucleic Acids Res. 2011, 39: W155-W159. 10.1093/nar/gkr319.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Zhang Y: miRU: an automated plant miRNA target prediction server. Nucleic Acids Res. 2005, 33: W701-W704. 10.1093/nar/gki383.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. German MA, Luo S, Schroth G, Meyers BC, Green PJ: Construction of Parallel Analysis of RNA Ends (PARE) libraries for the study of cleaved miRNA targets and the RNA degradome. Nat Protoc. 2009, 4: 356-362. 10.1038/nprot.2009.8.

    Article  CAS  PubMed  Google Scholar 

  20. German MA, Pillay M, Jeong DH, Hetawal A, Luo S, Janardhanan P, Kannan V, Rymarquis LA, Nobuta K, German R, De Paoli E, Lu C, Schroth G, Meyers BC, Green PJ: Global identification of microRNA-target RNA pairs by parallel analysis of RNA ends. Nat Biotechnol. 2008, 26: 941-946. 10.1038/nbt1417.

    Article  CAS  PubMed  Google Scholar 

  21. Meng Y, Shao C, Chen M: Toward microRNA-mediated gene regulatory networks in plants. Brief Bioinform. 2011, 12: 645-659. 10.1093/bib/bbq091.

    Article  CAS  PubMed  Google Scholar 

  22. Llave C, Xie Z, Kasschau KD, Carrington JC: Cleavage of Scarecrow-like mRNA targets directed by a class of Arabidopsis miRNA. Science. 2002, 297: 2053-2056. 10.1126/science.1076311.

    Article  CAS  PubMed  Google Scholar 

  23. Addo-Quaye C, Eshoo TW, Bartel DP, Axtell MJ: Endogenous siRNA and miRNA targets identified by sequencing of the Arabidopsis degradome. Curr Biol. 2008, 18: 758-762. 10.1016/j.cub.2008.04.042.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Addo-Quaye C, Snyder JA, Park YB, Li YF, Sunkar R, Axtell MJ: Sliced microRNA targets and precise loop-first processing of MIR319 hairpins revealed by analysis of the Physcomitrella patens degradome. RNA. 2009, 15: 2112-2121. 10.1261/rna.1774909.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Li YF, Zheng Y, Addo-Quaye C, Zhang L, Saini A, Jagadeeswaran G, Axtell MJ, Zhang W, Sunkar R: Transcriptome-wide identification of microRNA targets in rice. Plant J. 2010, 62: 742-759. 10.1111/j.1365-313X.2010.04187.x.

    Article  CAS  PubMed  Google Scholar 

  26. Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008, 36: D154-D158. 10.1093/nar/gkn221.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13: 2498-2504. 10.1101/gr.1239303.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Takehisa H, Sato Y, Igarashi M, Abiko T, Antonio BA, Kamatsuki K, Minami H, Namiki N, Inukai Y, Nakazono M, Nagamura Y: Genome-wide transcriptome dissection of the rice root system: implications for developmental and physiological functions. Plant J. 2012, 69: 126-140. 10.1111/j.1365-313X.2011.04777.x.

    Article  CAS  PubMed  Google Scholar 

  29. Yuan Q, Ouyang S, Liu J, Suh B, Cheung F, Sultana R, Lee D, Quackenbush J, Buell CR: The TIGR rice genome annotation resource: annotating the rice genome and creating resources for plant biologists. Nucleic Acids Res. 2003, 31: 229-233. 10.1093/nar/gkg059.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Huala E, Dickerman AW, Garcia-Hernandez M, Weems D, Reiser L, LaFond F, Hanley D, Kiphart D, Zhuang M, Huang W, Mueller LA, Bhattacharyya D, Bhaya D, Sobral BW, Beavis W, Meinke DW, Town CD, Somerville C, Rhee SY: The Arabidopsis Information Resource (TAIR): a comprehensive database and web-based information retrieval, analysis, and visualization system for a model plant. Nucleic Acids Res. 2001, 29: 102-105. 10.1093/nar/29.1.102.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Yan YS, Chen XY, Yang K, Sun ZX, Fu YP, Zhang YM, Fang RX: Overexpression of an F-box protein gene reduces abiotic stress tolerance and promotes root growth in rice. Mol Plant. 2011, 4: 190-197. 10.1093/mp/ssq066.

    Article  CAS  PubMed  Google Scholar 

  32. Boualem A, Laporte P, Jovanovic M, Laffont C, Plet J, Combier JP, Niebel A, Crespi M, Frugier F: MicroRNA166 controls root and nodule development in Medicago truncatula. Plant J. 2008, 54: 876-887. 10.1111/j.1365-313X.2008.03448.x.

    Article  CAS  PubMed  Google Scholar 

  33. Carlsbecker A, Lee JY, Roberts CJ, Dettmer J, Lehesranta S, Zhou J, Lindgren O, Moreno-Risueno MA, Vaten A, Thitamadee S, Campilho A, Sebastian J, Bowman JL, Helariutta Y, Benfey PN: Cell signalling by microRNA165/6 directs gene dose-dependent root cell fate. Nature. 2010, 465: 316-321. 10.1038/nature08977.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Grigg SP, Galinha C, Kornet N, Canales C, Scheres B, Tsiantis M: Repression of apical homeobox genes is required for embryonic root development in Arabidopsis. Curr Biol. 2009, 19: 1485-1490. 10.1016/j.cub.2009.06.070.

    Article  CAS  PubMed  Google Scholar 

  35. Qi Y, Wang S, Shen C, Zhang S, Chen Y, Xu Y, Liu Y, Wu Y, Jiang D: OsARF12, a transcription activator on auxin response gene, regulates root elongation and affects iron accumulation in rice (Oryza sativa). New Phytol. 2012, 193: 109-120. 10.1111/j.1469-8137.2011.03910.x.

    Article  CAS  PubMed  Google Scholar 

  36. Gifford ML, Dean A, Gutierrez RA, Coruzzi GM, Birnbaum KD: Cell-specific nitrogen responses mediate developmental plasticity. Proc Natl Acad Sci USA. 2008, 105: 803-808. 10.1073/pnas.0709559105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Tian CE, Muto H, Higuchi K, Matamura T, Tatematsu K, Koshiba T, Yamamoto KT: Disruption and overexpression of auxin response factor 8 gene of Arabidopsis affect hypocotyl elongation and root growth habit, indicating its possible involvement in auxin homeostasis in light condition. Plant J. 2004, 40: 333-343. 10.1111/j.1365-313X.2004.02220.x.

    Article  CAS  PubMed  Google Scholar 

  38. Kapulnik Y, Delaux PM, Resnick N, Mayzlish-Gati E, Wininger S, Bhattacharya C, Sejalon-Delmas N, Combier JP, Becard G, Belausov E, Beeckman T, Dor E, Hershenhorn J, Koltai H: Strigolactones affect lateral root formation and root-hair elongation in Arabidopsis. Planta. 2011, 233: 209-216. 10.1007/s00425-010-1310-y.

    Article  CAS  PubMed  Google Scholar 

  39. Mayzlish-Gati E, De-Cuyper C, Goormachtig S, Beeckman T, Vuylsteke M, Brewer PB, Beveridge CA, Yermiyahu U, Kaplan Y, Enzer Y, Wininger S, Resnick N, Cohen M, Kapulnik Y, Koltai H: Strigolactones are involved in root response to low phosphate conditions in Arabidopsis. Plant Physiol. 2012, 160: 1329-1341. 10.1104/pp.112.202358.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Rasmussen A, Mason MG, De Cuyper C, Brewer PB, Herold S, Agusti J, Geelen D, Greb T, Goormachtig S, Beeckman T, Beveridge CA: Strigolactones suppress adventitious rooting in Arabidopsis and pea. Plant Physiol. 2012, 158: 1976-1987. 10.1104/pp.111.187104.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Meng Y, Huang F, Shi Q, Cao J, Chen D, Zhang J, Ni J, Wu P, Chen M: Genome-wide survey of rice microRNAs and microRNA-target pairs in the root of a novel auxin-resistant mutant. Planta. 2009, 230: 883-898. 10.1007/s00425-009-0994-3.

    Article  CAS  PubMed  Google Scholar 

  42. Lichota J, Ritt C, Grasser KD: Ectopic expression of the maize chromosomal HMGB1 protein causes defects in root development of tobacco seedlings. Biochem Biophys Res Commun. 2004, 318: 317-322. 10.1016/j.bbrc.2004.03.193.

    Article  CAS  PubMed  Google Scholar 

  43. Li J, Zhu S, Song X, Shen Y, Chen H, Yu J, Yi K, Liu Y, Karplus VJ, Wu P, Deng XW: A rice glutamate receptor-like gene is critical for the division and survival of individual cells in the root apical meristem. Plant Cell. 2006, 18: 340-349. 10.1105/tpc.105.037713.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Vincill ED, Clarin AE, Molenda JN, Spalding EP: Interacting glutamate receptor-like proteins in phloem regulate lateral root initiation in Arabidopsis. Plant Cell. 2013, 10.1105/tpc.113.110668.

    Google Scholar 

  45. Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Muertter RN, Edgar R: NCBI GEO: archive for high-throughput functional genomic data. Nucleic Acids Res. 2009, 37: D885-890. 10.1093/nar/gkn764.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. Nakano M, Nobuta K, Vemaraju K, Tej SS, Skogen JW, Meyers BC: Plant MPSS databases: signature-based transcriptional resources for analyses of mRNA and small RNA. Nucleic Acids Res. 2006, 34: D731-735. 10.1093/nar/gkj077.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  47. Wu L, Zhang Q, Zhou H, Ni F, Wu X, Qi Y: Rice MicroRNA effector complexes and targets. Plant Cell. 2009, 21: 3421-3435. 10.1105/tpc.109.070938.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank all the publicly available datasets and the scientists behind them. This work was funded by the National Natural Science Foundation of China [31100937], [31271380] and [31125011], the Starting Grant funded by Hangzhou Normal University to Yijun Meng [2011QDL60].

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Chaogang Shao or Yijun Meng.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Conceived and designed the experiments: YM, CS, HW, YJ. Performed the experiments: XM, CS. Analyzed the data: XM, CS, YM. Contributed reagents/materials/analysis tools: XM, CS, YM. Wrote the paper: XM, YM, HW, YJ. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1:Data S1. List of sRNAs highly expressed in the root tips of rice. (TXT 6 MB)

Additional file 2:Data S2. List of sRNAs highly expressed in the whole roots of rice. (TXT 480 KB)

Additional file 3:Data S3. AGO1-enriched sRNAs highly expressed in the root tips of rice. (TXT 6 KB)

Additional file 4:Data S4. AGO1-enriched sRNAs highly expressed in the whole roots of rice. (TXT 9 KB)

12864_2013_5227_MOESM5_ESM.pdf

Additional file 5: Figure S1: Degradome sequencing data-based validation of the targets of the AGO1-enriched sRNAs highly expressed in the root tips of rice. (PDF 432 KB)

12864_2013_5227_MOESM6_ESM.pdf

Additional file 6: Figure S2: Degradome sequencing data-based validation of the targets of the AGO1-enriched sRNAs highly expressed in the whole roots of rice. (PDF 1 MB)

12864_2013_5227_MOESM7_ESM.pdf

Additional file 7: Figure S3: Network mediated by the AGO1-enriched sRNAs highly expressed in the root tips of rice. (PDF 5 KB)

12864_2013_5227_MOESM8_ESM.pdf

Additional file 8: Figure S4: Network mediated by the AGO1-enriched sRNAs highly expressed in the whole roots of rice. (PDF 20 KB)

12864_2013_5227_MOESM9_ESM.xls

Additional file 9: Table S1: Microarray data-based investigation of the expression patterns of the target genes. (XLS 110 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Ma, X., Shao, C., Wang, H. et al. Construction of small RNA-mediated gene regulatory networks in the roots of rice (Oryza sativa). BMC Genomics 14, 510 (2013). https://doi.org/10.1186/1471-2164-14-510

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-14-510

Keywords