Construction of gene regulatory networks mediated by vegetative and reproductive stage-specific small RNAs in rice (Oryza sativa)

Degradome sequencing reveals an integrative miRNA-mediated gene interaction network regulating rice seed vigor

BACKGROUND

It is well known that seed vigor is essential for agricultural production and rice (Oryza sativa L.) is one of the most important crops in the world. Though we previously reported that miR164c regulates rice seed vigor, but whether and how other miRNAs cooperate with miR164c to regulate seed vigor is still unknown.

RESULTS

Based on degradome data of six RNA samples isolated from seeds of the wild-type (WT) indica rice cultivar 'Kasalath' as well as two modified lines in 'Kasalath' background (miR164c-silenced line [MIM164c] and miR164c overexpression line [OE164c]), which were subjected to either no aging treatment or an 8-day artificial aging treatment, 1247 different target transcripts potentially cleaved by 421 miRNAs were identified. The miRNA target genes were functionally annotated via GO and KEGG enrichment analyses. By STRING database assay, a miRNA-mediated gene interaction network regulating seed vigor in rice was revealed, which comprised at least four interconnected pathways: the miR5075-mediated oxidoreductase related pathway, the plant hormone related pathway, the miR164e related pathway, and the previously reported RPS27AA related pathway. Knockout and overexpression of the target gene Os02g0817500 of miR5075 decreased and enhanced seed vigor, respectively. By Y2H assay, the proteins encoded by five seed vigor-related genes, Os08g0295100, Os07g0633100, REFA1, OsPER1 and OsGAPC3, were identified to interact with Os02g0817500.

CONCLUSIONS

miRNAs cooperate to regulate seed vigor in rice via an integrative gene interaction network comprising miRNA target genes and other functional genes. The result provided a basis for fully understanding the molecular mechanisms of seed vigor regulation.

A rice transcription factor, OsMADS57, positively regulates high salinity tolerance in transgenic Arabidopsis thaliana and Oryza sativa plants

MADS-box transcription factors (TFs) play indispensable roles in various aspects of plant growth, development as well as in response to environmental stresses. Several MADS-box genes have been reported to be involved in the salt tolerance in different plant species. However, the role of the transcription factor OsMADS57 under salinity stress is still unknown. Here, the results of this study showed that OsMADS57 was mainly expressed in roots and leaves of rice plants (Oryza sativa). Gene expression pattern analysis revealed that OsMADS57 was induced by NaCl. Overexpression of OsMADS57 in both Arabidopsis thaliana (A. thaliana) and rice could improve their salt tolerance, which was demonstrated by higher germination rates, longer root length and better growth status of overexpression plants than wild type (WT) under salinity conditions. In contrast, RNA interference (RNAi) lines of rice showed more sensitivity towards salinity. Moreover, less reactive oxygen species (ROS) accumulated in OsMADS57 overexpressing lines when exposed to salt stress, as measured by 3, 3'-diaminobenzidine (DAB) or nitroblue tetrazolium (NBT) staining. Further experiments exhibited that overexpression of OsMADS57 in rice significantly increased the tolerance ability of plants to oxidative damage under salt stress, mainly by increasing the activities of antioxidative enzymes such as superoxide dismutase (SOD) and peroxidase (POD), reducing malonaldehyde (MDA) content and improving the expression of stress-related genes. Taken together, these results demonstrated that OsMADS57 plays a positive role in enhancing salt tolerance by activating the antioxidant system.

Reproductive phasiRNAs regulate reprogramming of gene expression and meiotic progression in rice

Plant spermatogenesis is a complex process that directly affects crop breeding. A rapid change in gene abundance occurs at early meiosis prophase, when gene regulation is selective. However, how these genes are regulated remains unknown. Here, we show that rice reproductive phasiRNAs are essential for the elimination of a specific set of RNAs during meiotic prophase I. These phasiRNAs cleave target mRNAs in a regulatory manner such that one phasiRNA can target more than one gene, and/or a single gene can be targeted by more than one phasiRNA to efficiently silence target genes. Our investigation of phasiRNA-knockdown and PHAS-edited transgenic plants demonstrates that phasiRNAs and their nucleotide variations are required for meiosis progression and fertility. This study highlights the importance of reproductive phasiRNAs for the reprogramming of gene expression during meiotic progression and establishes a basis for future studies on the roles of phasiRNAs with a goal of crop improvement.

A potential endogenous gibberellin-mediated signaling cascade regulated floral transition in Magnolia × soulangeana 'Changchun'

The floral transition is a critical developmental switch in plants, and has profound effects on the flower production and yield. Magnolia × soulangeana 'Changchun' is known as a woody ornamental plant, which can bloom in spring and summer, respectively. In this study, anatomical observation, physiological measurement, transcriptome, and small RNA sequencing were performed to investigate potential endogenous regulatory mechanisms underlying floral transition in 'Changchun'. Transition of the shoot apical meristem from vegetative to reproductive growth occurred between late April and early May. During this specific developmental process, a total of 161,645 unigenes were identified, of which 73,257 were significantly differentially expressed, while a number of these two categories of miRNAs were 299 and 148, respectively. Further analysis of differentially expressed genes (DEGs) revealed that gibberellin signaling could regulate floral transition in 'Changchun' in a DELLA-dependent manner. In addition, prediction and analysis of miRNA targeted genes suggested that another potential molecular regulatory module was mediated by the miR172 family and other several novel miRNAs (Ms-novel_miR139, Ms-novel_miR229, and Ms-novel_miR232), with the participation of up- or down-regulating genes, including MsSVP, MsAP2, MsTOE3, MsAP1, MsGATA6, MsE2FA, and MsMDS6. Through the integrated analysis of mRNA and miRNA, our research results will facilitate the understanding of the potential molecular mechanism underlying floral transition in 'Changchun', and also provide basic experimental data for the plant germplasm resources innovation in Magnolia.

Mining Magnaporthe oryzae sRNAs With Potential Transboundary Regulation of Rice Genes Associated With Growth and Defense Through Expression Profile Analysis of the Pathogen-Infected Rice

In recent years, studies have shown that phytopathogenic fungi possess the ability of cross-kingdom regulation of host plants through small RNAs (sRNAs). Magnaporthe oryzae, a causative agent of rice blast, introduces disease by penetrating the rice tissues through appressoria. However, little is known about the transboundary regulation of M. oryzae sRNAs during the interaction of the pathogen with its host rice. Therefore, investigation of the regulation of M. oryzae through sRNAs in the infected rice plants has important theoretical and practical significance for disease control and production improvement. Based on the high-throughput data of M. oryzae sRNAs and the mixed sRNAs during infection, the differential expressions of sRNAs in M. oryzae before and during infection were compared, it was found that expression levels of 366 M. oryzae sRNAs were upregulated significantly during infection. We trained a SVM model which can be used to predict differentially expressed sRNAs, which has reference significance for the prediction of differentially expressed sRNAs of M. oryzae homologous species, and can facilitate the research of M. oryzae in the future. Furthermore, fifty core targets were selected from the predicted target genes on rice for functional enrichment analysis, the analysis reveals that there are nine biological processes and one KEGG pathway associated with rice growth and disease defense. These functions correspond to thirteen rice genes. A total of fourteen M. oryzae sRNAs targeting the rice genes were identified by data analysis, and their authenticity was verified in the database of M. oryzae sRNAs. The 14 M. oryzae sRNAs may participate in the transboundary regulation process and act as sRNA effectors to manipulate the rice blast process.

Identification and expression analysis of microRNAs during ovule development in rice (Oryza sativa) by deep sequencing

MicroRNA (miRNA) expression profiles during rice ovule development revealed the possible miRNA-mediated regulation between ovule sporophytic tissue and female gametophyte and the involvement of miRNAs in programmed cell death. MiRNAs are 20-24-nucleotide small RNAs that play key roles in the regulation of many growth and developmental processes in plants. Rice ovule development comprises a series of biological events, which are regulated by complex molecular mechanisms. To gain insight into miRNA-mediated regulation of rice ovule development, Illumina sequencing was used to examine the expression of miRNAs from the megaspore mother cell meiosis stage to the fertilized ovule stage. Based on the sequencing data, 486 known and 204 novel miRNAs were identified during rice ovule development. Moreover, 56, 65 and 11 differentially expressed miRNAs between adjacent developmental stages were identified. By analyzing transcriptome and degradome data, we identified 41, 65 and 12 coherent target genes for the differentially expressed miRNAs in ovule development. We found that changes in the expression of plant hormone-related miRNAs may play important roles in embryo sac development, providing evidence for cross-talk communication between sporophytic tissue and the female gametophyte. Additionally, we revealed that miRNAs may be involved in programmed cell death after fertilization. Finally, we constructed miRNA-mediated regulatory networks that are active during rice ovule development.

Comparative Small RNA Analysis of Pollen Development in Autotetraploid and Diploid Rice

MicroRNAs (miRNAs) play key roles in plant reproduction. However, knowledge on microRNAome analysis in autotetraploid rice is rather limited. Here, high-throughput sequencing technology was employed to analyze miRNAomes during pollen development in diploid and polyploid rice. A total of 172 differentially expressed miRNAs (DEM) were detected in autotetraploid rice compared to its diploid counterpart, and 57 miRNAs were specifically expressed in autotetraploid rice. Of the 172 DEM, 115 and 61 miRNAs exhibited up- and down-regulation, respectively. Gene Ontology analysis on the targets of up-regulated DEM showed that they were enriched in transport and membrane in pre-meiotic interphase, reproduction in meiosis, and nucleotide binding in single microspore stage. osa-miR5788 and osa-miR1432-5p_R+1 were up-regulated in meiosis and their targets revealed interaction with the meiosis-related genes, suggesting that they may involve in the genes regulation associated with the chromosome behavior. Abundant 24 nt siRNAs associated with transposable elements were found in autotetraploid rice during pollen development; however, they significantly declined in diploid rice, suggesting that 24 nt siRNAs may play a role in pollen development. These findings provide a foundation for understanding the effect of polyploidy on small RNA expression patterns during pollen development that cause pollen sterility in autotetraploid rice.

A network analysis of miRNA mediated gene regulation of rice: crosstalk among biological processes

To understand the network architecture of miRNA mediated regulations at the genomic and functional levels of rice, we have made an unambiguous annotation of the experimentally verified miRNAs, predicted their targets and the possible biological functions they can affect. Some functions, namely translational and protein modifications and photosynthesis are targeted by higher percentage of miRNA. Using transformation procedures, we constructed a genome scale miRNA-miRNA functional synergistic network (MFSN). The analysis of MFSN modules help to identify miRNAs co-regulating target genes having several interrelated biological processes. Some of these target genes are also co-expressed under particular conditions. For example, the genes co-expressed under drought conditions as well as those targeted by miRNAs present in a MFSN module have interdependent biological processes namely, photosynthesis, cell-wall biogenesis, root development and xylan synthesis. The stress-induced miRNAs and their distributions, and the presence of transcription factors in the target set of MFSN modules were also analyzed.

Combined Small RNA and Degradome Sequencing Reveals Novel MiRNAs and Their Targets in the High-Yield Mutant Wheat Strain Yunong 3114

Wheat is one of the main food sources worldwide; large amount studies have been conducted to improve wheat production. MicroRNAs (miRNAs) with about 20-30 nucleotide are a class of regulatory small RNAs (sRNAs), which could regulate gene expression through sequence-specific base pairing with target mRNAs, playing important roles in plant growth. An ideal plant architecture (IPA) is crucial to enhance yield in bread wheat. In this study, the high-yield wheat strain Yunong 3114 was EMS-mutagenesis from the wild-type strain Yunong 201, exhibiting a preferable plant structure compared with the wild-type strain. We constructed small RNA and degradome libraries from Yunong 201 and Yunong 3114, and performed small RNA sequencing of these libraries in order identify miRNAs and their targets related to IPA in wheat. Totally, we identified 488 known and 837 novel miRNAs from Yunong 3114 and 391 known and 533 novel miRNAs from Yunong 201. The number of miRNAs in the mutant increased. A total of 37 known and 432 putative novel miRNAs were specifically expressed in the mutant strain; furthermore, 23 known and 159 putative novel miRNAs were specifically expressed in the wild-type strain. A total of 150 known and 100 novel miRNAs were differentially expressed between mutant and wild-type strains. Among these differentially expressed novel miRNAs, 4 and 8 predict novel miRNAs were evidenced by degradome sequencing and showed up-regulated and down-regulated expressions in the mutant strain Yunong 3114, respectively. Targeted gene annotation and previous results indicated that this set of miRNAs is related to plant structure. Our results further suggested that miRNAs may be necessary to obtain an optimal wheat structure.

MADS-box transcription factor OsMADS25 regulates root development through affection of nitrate accumulation in rice

MADS-box transcription factors are vital regulators participating in plant growth and development process and the functions of most of them are still unknown. ANR1 was reported to play a key role in controlling lateral root development through nitrate signal in Arabidopsis. OsMADS25 is one of five ANR1-like genes in Oryza Sativa and belongs to the ANR1 clade. Here we have investigated the role of OsMADS25 in the plant’s responses to external nitrate in Oryza Sativa. Our results showed that OsMADS25 protein was found in the nucleus as well as in the cytoplasm. Over-expression of OsMADS25 significantly promoted lateral and primary root growth as well as shoot growth in a nitrate-dependent manner in Arabidopsis. OsMADS25 overexpression in transgenic rice resulted in significantly increased primary root length, lateral root number, lateral root length and shoot fresh weight in the presence of nitrate. Down-regulation of OsMADS25 in transgenic rice exhibited significantly reduced shoot and root growth in the presence of nitrate. Furthermore, over-expression of OsMADS25 in transgenic rice promoted nitrate accumulation and significantly increased the expressions of nitrate transporter genes at high rates of nitrate supply while down-regulation of OsMADS25 produced the opposite effect. Taken together, our findings suggest that OsMADS25 is a positive regulator control lateral and primary root development in rice.

The effects of fluctuations in the nutrient supply on the expression of five members of the AGL17 clade of MADS-box genes in rice

The ANR1 MADS-box gene in Arabidopsis is a key gene involved in regulating lateral root development in response to the external nitrate supply. There are five ANR1-like genes in Oryza sativa, OsMADS23, OsMADS25, OsMADS27, OsMADS57 and OsMADS61, all of which belong to the AGL17 clade. Here we have investigated the responsiveness of these genes to fluctuations in nitrogen (N), phosphorus (P) and sulfur (S) mineral nutrient supply. The MADS-box genes have been shown to have a range of responses to the nutrient supply. The expression of OsMADS61 was transiently induced by N deprivation but was not affected by re-supply with various N sources. The expression of OsMADS25 and OsMADS27 was induced by re-supplying with NO3(-) and NH4NO3, but downregulated by NH4(+). The expression of OsMADS57 was significantly downregulated by N starvation and upregulated by 3 h NO3(-) re-supply. OsMADS23 was the only gene that showed no response to either N starvation nor NO3(-) re-supply. OsMADS57 was the only gene not regulated by P fluctuation whereas the expression of OsMADS23, OsMADS25 and OsMADS27 was downregulated by P starvation and P re-supply. In contrast, all five ANR1-related genes were significantly upregulated by S starvation. Our results also indicated that there were interactions among nitrate, sulphate and phosphate transporters in rice.

Trans-acting small interfering RNA4: key to nutraceutical synthesis in grape development?

The facility and versatility of microRNAs (miRNAs) to evolve and change likely underlies how they have become dominant constituents of eukaryotic genomes. In this opinion article I propose that trans-acting small interfering RNA gene 4 (TAS4) evolution may be important for biosynthesis of polyphenolics, arbuscular symbiosis, and bacterial pathogen etiologies. Expression-based and phylogenetic evidence shows that TAS4 targets two novel grape (Vitis vinifera L.) MYB transcription factors (VvMYBA6, VvMYBA7) that spawn phased small interfering RNAs (siRNAs) which probably function in nutraceutical bioflavonoid biosynthesis and fruit development. Characterization of the molecular mechanisms of TAS4 control of plant development and integration into biotic and abiotic stress- and nutrient-signaling regulatory networks has applicability to molecular breeding and the development of strategies for engineering healthier foods.