Identification of novel MiRNAs and MiRNA expression profiling during grain development in indica rice

Identification of miRNAs associated with Aspergillus flavus infection and their targets in groundnut (Arachis hypogaea L.)

BACKGROUND

The quality of groundnut produce is adversely impacted due to aflatoxin contamination by the fungus Aspergillus flavus. Although the transcriptomic control is not fully understood, the interaction between long non-coding RNAs and microRNAs in regulating A. flavus and aflatoxin contamination remains unclear. This study was carried out to identify microRNAs (miRNAs) to enhance the understanding of in vitro seed colonization (IVSC) resistance mechanism in groundnut.

RESULT

In this study, resistant (J 11) and susceptible (JL 24) varieties of groundnut were treated with toxigenic A. flavus (strain AF-11-4), and total RNA was extracted at 1 day after inoculation (1 DAI), 2 DAI, 3 DAI and 7 DAI. Seeds of JL 24 showed higher mycelial growth than J 11 at successive days after inoculation. A total of 208 known miRNAs belonging to 36 miRNA families, with length varying from 20-24 nucleotides, were identified, along with 27 novel miRNAs, with length varying from 20-22 nucleotides. Using psRNATarget server, 952 targets were identified for all the miRNAs. The targeted genes function as disease resistant proteins encoding, auxin responsive proteins, squamosa promoter binding like proteins, transcription factors, pentatricopeptide repeat-containing proteins and growth regulating factors. Through differential expression analysis, seven miRNAs (aly-miR156d-3p, csi-miR1515a, gma-miR396e, mtr-miR2118, novo-miR-n27, ptc-miR482d-3p and ppe-miR396a) were found common among 1 DAI, 2 DAI, 3 DAI and 7 DAI in J 11, whereas ten miRNAs (csi-miR159a-5p, csi-miR164a-3p, novo-miR-n17, novo-miR-n2, osa-miR162b, mtr-miR2118, ptc-miR482d-3p, ptc-miR167f-3p, stu-miR319-3p and zma-miR396b-3p) were found common among 1 DAI, 2 DAI, 3 DAI and 7 DAI in JL 24. Two miRNAs, ptc-miR482d-3p and mtr-miR2118, showed contrasting expression at different time intervals between J 11 and JL 24. These two miRNAs were found to target those genes with NBS-LRR function, making them potential candidates for marker development in groundnut breeding programs aimed at enhancing resistance against A. flavus infection.

CONCLUSION

This study enhances our understanding of the involvement of two miRNAs namely, ptc-miR482d-3p and mtr-miR2118, along with their NBS-LRR targets, in conferring resistance against A. flavus-induced aflatoxin contamination in groundnut under in vitro conditions.

Integrated mRNA and miRNA analysis reveals the regulatory network of oxidative stress and inflammation in Coilia nasus brains during air exposure and salinity mitigation

BACKGROUND

Air exposure is an inevitable source of stress that leads to significant mortality in Coilia nasus. Our previous research demonstrated that adding 10‰ NaCl to aquatic water could enhance survival rates, albeit the molecular mechanisms involved in air exposure and salinity mitigation remained unclear. Conversely, salinity mitigation resulted in decreased plasma glucose levels and improved antioxidative activity. To shed light on this phenomenon, we characterized the transcriptomic changes in the C. nasus brain upon air exposure and salinity mitigation by integrated miRNA-mRNA analysis.

RESULTS

The plasma glucose level was elevated during air exposure, whereas it decreased during salinity mitigation. Antioxidant activity was suppressed during air exposure, but was enhanced during salinity mitigation. A total of 629 differentially expressed miRNAs (DEMs) and 791 differentially expressed genes (DEGs) were detected during air exposure, while 429 DEMs and 1016 DEGs were identified during salinity mitigation. GO analysis revealed that the target genes of DEMs and DEGs were enriched in biological process and cellular component during air exposure and salinity mitigation. KEGG analysis revealed that the target genes of DEMs and DEGs were enriched in metabolism. Integrated analysis showed that 24 and 36 predicted miRNA-mRNA regulatory pairs participating in regulating glucose metabolism, Ca2+ transport, inflammation, and oxidative stress. Interestingly, most of these miRNAs were novel miRNAs.

CONCLUSION

In this study, substantial miRNA-mRNA regulation pairs were predicted via integrated analysis of small RNA sequencing and RNA-Seq. Based on predicted miRNA-mRNA regulation and potential function of DEGs, miRNA-mRNA regulatory network involved in glucose metabolism and Ca2+ transport, inflammation, and oxidative stress in C. nasus brain during air exposure and salinity mitigation. They regulated the increased/decreased plasma glucose and inhibited/promoted antioxidant activity during air exposure and salinity mitigation. Our findings would propose novel insights to the mechanisms underlying fish responses to air exposure and salinity mitigation.

High Daytime Temperature Responsive MicroRNA Profiles in Developing Grains of Rice Varieties with Contrasting Chalkiness

High temperature impairs starch biosynthesis in developing rice grains and thereby increases chalkiness, affecting the grain quality. Genome encoded microRNAs (miRNAs) fine-tune target transcript abundances in a spatio-temporal specific manner, and this mode of gene regulation is critical for a myriad of developmental processes as well as stress responses. However, the role of miRNAs in maintaining rice grain quality/chalkiness during high daytime temperature (HDT) stress is relatively unknown. To uncover the role of miRNAs in this process, we used five contrasting rice genotypes (low chalky lines Cyp, Ben, and KB and high chalky lines LaGrue and NB) and compared the miRNA profiles in the R6 stage caryopsis samples from plants subjected to prolonged HDT (from the onset of fertilization through R6 stage of caryopsis development). Our small RNA analysis has identified approximately 744 miRNAs that can be grouped into 291 families. Of these, 186 miRNAs belonging to 103 families are differentially regulated under HDT. Only two miRNAs, Osa-miR444f and Osa-miR1866-5p, were upregulated in all genotypes, implying that the regulations greatly varied between the genotypes. Furthermore, not even a single miRNA was commonly up/down regulated specifically in the three tolerant genotypes. However, three miRNAs (Osa-miR1866-3p, Osa-miR5150-3p and canH-miR9774a,b-3p) were commonly upregulated and onemiRNA (Osa-miR393b-5p) was commonly downregulated specifically in the sensitive genotypes (LaGrue and NB). These observations suggest that few similarities exist within the low chalky or high chalky genotypes, possibly due to high genetic variation. Among the five genotypes used, Cypress and LaGrue are genetically closely related, but exhibit contrasting chalkiness under HDT, and thus, a comparison between them is most relevant. This comparison revealed a general tendency for Cypress to display miRNA regulations that could decrease chalkiness under HDT compared with LaGrue. This study suggests that miRNAs could play an important role in maintaining grain quality in HDT-stressed rice.

Cytological, transcriptome and miRNome temporal landscapes decode enhancement of rice grain size

BACKGROUND

Rice grain size (GS) is an essential agronomic trait. Though several genes and miRNA modules influencing GS are known and seed development transcriptomes analyzed, a comprehensive compendium connecting all possible players is lacking. This study utilizes two contrasting GS indica rice genotypes (small-grained SN and large-grained LGR). Rice seed development involves five stages (S1-S5). Comparative transcriptome and miRNome atlases, substantiated with morphological and cytological studies, from S1-S5 stages and flag leaf have been analyzed to identify GS proponents.

RESULTS

Histology shows prolonged endosperm development and cell enlargement in LGR. Stand-alone and comparative RNAseq analyses manifest S3 (5-10 days after pollination) stage as crucial for GS enhancement, coherently with cell cycle, endoreduplication, and programmed cell death participating genes. Seed storage protein and carbohydrate accumulation, cytologically and by RNAseq, is shown to be delayed in LGR. Fourteen transcription factor families influence GS. Pathway genes for four phytohormones display opposite patterns of higher expression. A total of 186 genes generated from the transcriptome analyses are located within GS trait-related QTLs deciphered by a cross between SN and LGR. Fourteen miRNA families express specifically in SN or LGR seeds. Eight miRNA-target modules display contrasting expressions amongst SN and LGR, while 26 (SN) and 43 (LGR) modules are differentially expressed in all stages.

CONCLUSIONS

Integration of all analyses concludes in a "Domino effect" model for GS regulation highlighting chronology and fruition of each event. This study delineates the essence of GS regulation, providing scope for future exploits. The rice grain development database (RGDD) ( www.nipgr.ac.in/RGDD/index.php ; https://doi.org/10.5281/zenodo.7762870 ) has been developed for easy access of data generated in this paper.

Dynamic Regulation of Grapevine's microRNAs in Response to Mycorrhizal Symbiosis and High Temperature

MicroRNAs (miRNAs) are non-coding small RNAs that play crucial roles in plant development and stress responses and can regulate plant interactions with beneficial soil microorganisms such as arbuscular mycorrhizal fungi (AMF). To determine if root inoculation with distinct AMF species affected miRNA expression in grapevines subjected to high temperatures, RNA-seq was conducted in leaves of grapevines inoculated with either Rhizoglomus irregulare or Funneliformis mosseae and exposed to a high-temperature treatment (HTT) of 40 °C for 4 h per day for one week. Our results showed that mycorrhizal inoculation resulted in a better plant physiological response to HTT. Amongst the 195 identified miRNAs, 83 were considered isomiRs, suggesting that isomiRs can be biologically functional in plants. The number of differentially expressed miRNAs between temperatures was higher in mycorrhizal (28) than in non-inoculated plants (17). Several miR396 family members, which target homeobox-leucine zipper proteins, were only upregulated by HTT in mycorrhizal plants. Predicted targets of HTT-induced miRNAs in mycorrhizal plants queried to STRING DB formed networks for Cox complex, and growth and stress-related transcription factors such as SQUAMOSA promoter-binding-like-proteins, homeobox-leucine zipper proteins and auxin receptors. A further cluster related to DNA polymerase was found in R. irregulare inoculated plants. The results presented herein provide new insights into miRNA regulation in mycorrhizal grapevines under heat stress and can be the basis for functional studies of plant-AMF-stress interactions.

MicroRNAs modulate ethylene induced retrograde signal for rice endosperm starch biosynthesis by default expression of transcriptome

Control of stage specific spike in ethylene production at anthesis has been a vauable route to potentially enhance genetic ceiling for grain filling of rice spikelet. A number of genes controlling ethylene homeostasis and starch synthesis have been identified so long, but lack of credible information on master modulation of gene expression by miRNAs and their target genes associated with hormonal dynamics obfuscate mechanisms controlling genotype difference in quantum of grain filling. The confusion accounts for consequent shrinkage of options for yield manipulation. In a two by two factorial design, miRNA regulation of spikelet specific grain development in low against high sterile recombinant inbred lines of rice Oryza sativa L. namely CR 3856-62-11-3-1-1-1-1-1-1 (SR 157) and CR 3856-63-1-1-1-1-1-1-1 (SR 159) respectively, and inferior verses superior spikelets were compared during first 10 days after anthesis. Grain filling was poorer in SR159 than SR157 and inferior spikelets in the former were most vulnerable. Between the cultivars, overall expression of unique miRNAs with targets on ethylene pathway genes was higher in SR159 than SR157 and the situation was opposite for auxin pathway genes. Precision analysis in psTarget server database identified up-regulation of MIR2877 and MIR530-5p having Os11t0141000-02 and Os07t0239400-01 (PP2A regulatory subunit-like protein and ethylene-responsive small GTP-binding proteins) and MIR396h having Os01t0643300-02 (an auxin efflux carrier protein) and Os01t0643300-01 (a PIN1-like auxin transport protein), as targets with highest probability at anthesis and 5 days after anthesis respectively, in the inferior spikelet and the fold change values of DGE matched with pattern of gene expression (relative transcript level) in the qRT-PCR studies conducted for relevant miRNAs and protein factors for ethylene and auxin signalling. In conclusion, epigenetic regulation of both auxin and ethylene homeostasis control grain filling of rice spikelet was established, but evidences were more robust for the latter.

Study of expressions of miRNAs in the spikelets based on their spatial location on panicle in rice cultivars provided insight into their influence on grain development

Development of rice cultivars bearing numerous spikelets by breeding approach to increase the yearly production of rice to approximately 800 million metric tons to feed the ever increasing population of the world accompanies poor grain filling in the inferior spikelets preventing achievement of the yield potential. As the initial stages of caryopses development are of much importance for grain filling, spatio-temporal expressions of the miRNAs were studied during these periods in the spikelets of a compact-panicle rice cultivar, Oryza sativa cv. Mahalaxmi, bearing numerous spikelets per panicle to understand the reason of poor grain filling at the level of the initial biochemical events. Differential expression of several known miRNAs between the superior and inferior spikelets suggested great difference in metabolism related to grain filling in the spikelets based on their spatial location on compact panicle. Expressions of five known and four novel miRNAs were validated by Northern. Their targets included the enzymes directly involved in starch biosynthesis like sucrose synthase, starch synthase and pullulanase, besides others. Spatio-temporal expression studies of these miRNAs in the spikelets of Mahalaxmi revealed a pattern of mostly a greater expression in the inferior spikelets compared with the superior ones concomitant with an inverse expression of the target genes, which was not observed in the lax-panicle cultivar Upahar. The study thus revealed that the grain filling in rice is greatly regulated by miRNAs, and these miRNAs or their target genes could be considered for biotechnological interventions for improving grain filling in the rice cultivars of interest.

MiR396 regulatory network and its expression during grain development in wheat

Wheat contains the largest number of miR396 family with 17 miR396 in Poaceae. MiR396 regulatory network underlying wheat grain development has not comprehensively been explored. Our results showed that precursor miR396 family in Poaceae exhibited not only conservativeness but also diversification especially in wheat. Five haplotypes were detected in Poaceae species, while 4 haplotypes in wheat with Hap-4 (miR396a) and Hap-5 (miR396n) unique to wheat. GO enrichment analysis of target genes showed that the first 20 enrichment functions of miR396a and miR396n are completely different from each other, and also completely different from miR396(b-g), miR396(h-m), and miR396(o-q). Functional annotation on the 18 target genes shared by miR396(b-g), miR396(h-m), and miR396(o-q) found that 11 of the 18 target genes are growth-regulating factor (GRF) genes. Our results indicated that, during the grain filling stage of wheat, miR396 is involved in the development of grains by regulating the expression of GRF genes (GRF1, GRF6, and GRF9). Although the enrichment function of miR396(b-g), miR396(h-m), and miR396(o-q) is the same, the gene functional networks they formed differ greatly. Our results indicated that polyploidization enriches not only the diversity of miR396 family and its target genes but also gene functional networks in wheat. These results laid foundation for further elucidating function of miR396 gene family underlying wheat grain development.

Multi-Omics Analysis of Small RNA, Transcriptome, and Degradome in T. turgidum-Regulatory Networks of Grain Development and Abiotic Stress Response

Crop reproduction is highly sensitive to water deficit and heat stress. The molecular networks of stress adaptation and grain development in tetraploid wheat (Triticum turgidum durum) are not well understood. Small RNAs (sRNAs) are important epigenetic regulators connecting the transcriptional and post-transcriptional regulatory networks. This study presents the first multi-omics analysis of the sRNAome, transcriptome, and degradome in T. turgidum developing grains, under single and combined water deficit and heat stress. We identified 690 microRNAs (miRNAs), with 84 being novel, from 118 sRNA libraries. Complete profiles of differentially expressed miRNAs (DEMs) specific to genotypes, stress types, and different reproductive time-points are provided. The first degradome sequencing report for developing durum grains discovered a significant number of new target genes regulated by miRNAs post-transcriptionally. Transcriptome sequencing profiled 53,146 T. turgidum genes, swith differentially expressed genes (DEGs) enriched in functional categories such as nutrient metabolism, cellular differentiation, transport, reproductive development, and hormone transduction pathways. miRNA-mRNA networks that affect grain characteristics such as starch synthesis and protein metabolism were constructed on the basis of integrated analysis of the three omics. This study provides a substantial amount of novel information on the post-transcriptional networks in T. turgidum grains, which will facilitate innovations for breeding programs aiming to improve crop resilience and grain quality.

miR164c and miR168a regulate seed vigor in rice

MicroRNAs (miRNAs) are key regulators of gene expression in many important biological processes of plants. However, few miRNAs have been shown to regulate seed vigor. Here, we conducted microarray assays to analyze miRNA expression levels in seeds of the rice (Oryza sativa L.) cultivar ZR02. Results showed significant differences in the expression of 11 miRNAs between artificially aged and untreated control seeds. Among these, osa-miR164c was transcriptionally upregulated, while osa-miR168a was downregulated in artificially aged seeds; this was verified by quantitative real-time PCR analysis. Under the same aging condition, osa-miR164c overexpression in OE164c transgenic seeds and osa-miR168a silencing in MIM168a transgenic seeds of the rice cultivar Kasalath led to lower germination rates, whereas osa-miR164c silencing in MIM164c and osa-miR168a overexpression in OE168a resulted in higher seed germination rates compared with wild-type seeds. Meanwhile, changes in cytomembrane permeability of seeds and in the expression level of osa-miR164c target genes (OsPM27 and OsPSK5) and osa-miR168a target genes (OsAGO1 and OsPTR2) under aging conditions coincided with changes in seed vigor induced by osa-miR164c and osa-miR168a. Thus, genetic manipulation of miRNAs has important implications in the development of crop cultivars with high vigor and extended life span of seeds.

Comparative Analysis between Wild and Cultivated Cucumbers Reveals Transcriptional Changes during Domestication Process

The cultivated cucumber (Cucumis sativus L.) was reported to have been developed from a wild cucumber (Cucumis hystrix Chakrav.), nevertheless, these two organisms exhibit noteworthy differences. For example, the wild cucumber is known for its high resistance to different biotic and abiotic stresses. Moreover, the leaves and fruits of the wild cucumber have a bitter taste compared to the cultivated cucumber. These differences could be attributed mainly to the differences in gene expression levels. In the present investigation, we analyzed the RNA-sequencing data to show the differentially expressed genes (DEGs) between the wild and cultivated cucumbers. The identified DEGs were further utilized for Gene Ontology (GO) and pathway enrichment analysis and for identification of transcription factors and regulators. In the results, several enriched GO terms in the biological process, cellular component, and molecular functions categories were identified and various enriched pathways, especially the biosynthesis pathways of secondary products were recognized. Plant-specific transcription factor families were differentially expressed between the wild and cultivated cucumbers. The results obtained provide preliminary evidence for the transcriptional differences between the wild and cultivated cucumbers which developed during the domestication process as a result of natural and/or artificial selection, and they formulate the basis for future genetic research and improvement of the cultivated cucumber.

A missense mutation in Large Grain Size 1 increases grain size and enhances cold tolerance in rice

Grain shape is controlled by quantitative trait loci (QTLs) in rice (Oryza sativa L.). A rice mutant (JF178) with long and large grains has been used in a breeding program for over a decade, but its genetic basis has been unclear. Here, a semi-dominant QTL, designated Large Grain Size 1 (LGS1), was cloned and the potential molecular mechanism of LGS1 function was studied. Near-isogenic lines (NILs) and a map-based approach were employed to clone the LGS1 locus. LGS1 encodes the OsGRF4 transcription factor and contains a 2 bp missense mutation in the coding region that coincides with the putative pairing site of miRNA396. The LGS1 transcript levels in the mutant line were found to be higher than the lgs1 transcript levels in the control plants, suggesting that the mutation might disrupt the pairing of the LGS1 mRNA with miR396. In addition to producing larger grains, LGS1 also enhanced cold tolerance at the seedling stage and increased the survival rate of seedlings after cold stress treatment. These findings indicate that the mutation in LGS1 appears to disturb the GRF4-miR396 stress response network and results in the development of enlarged grains and enhancement of cold tolerance in rice.

Molecular insights into the regulation of rice kernel elongation

A large number of rice agronomic traits are complex, multi factorial and polygenic. As the mechanisms and genes determining grain size and yield are largely unknown, the identification of regulatory genes related to grain development remains a preeminent approach in rice genetic studies and breeding programs. Genes regulating cell proliferation and expansion in spikelet hulls and participating in endosperm development are the main controllers of rice kernel elongation and grain size. We review here and discuss recent findings on genes controlling rice grain size and the mechanisms, epialleles, epigenomic variation, and assessment of controlling genes using genome-editing tools relating to kernel elongation.

A Resource for Inactivation of MicroRNAs Using Short Tandem Target Mimic Technology in Model and Crop Plants

microRNAs (miRNAs) are endogenous small non-coding RNAs that bind to mRNAs and target them for cleavage and/or translational repression, leading to gene silencing. We previously developed short tandem target mimic (STTM) technology to deactivate endogenous miRNAs in Arabidopsis. Here, we created hundreds of STTMs that target both conserved and species-specific miRNAs in Arabidopsis, tomato, rice, and maize, providing a resource for the functional interrogation of miRNAs. We not only revealed the functions of several miRNAs in plant development, but also demonstrated that tissue-specific inactivation of a few miRNAs in rice leads to an increase in grain size without adversely affecting overall plant growth and development. RNA-seq and small RNA-seq analyses of STTM156/157 and STTM165/166 transgenic plants revealed the roles of these miRNAs in plant hormone biosynthesis and activation, secondary metabolism, and ion-channel activity-associated electrophysiology, demonstrating that STTM technology is an effective approach for studying miRNA functions. To facilitate the study and application of STTM transgenic plants and to provide a useful platform for storing and sharing of information about miRNA-regulated gene networks, we have established an online Genome Browser (https://blossom.ffr.mtu.edu/designindex2.php) to display the transcriptomic and miRNAomic changes in STTM-induced miRNA knockdown plants.

Small RNA and degradome profiling involved in seed development and oil synthesis of Brassica napus

MicroRNAs (miRNAs) play a prominent role in post-transcriptional gene expression regulation and have been involved in various biological and metabolic processes to regulate gene expression. For Brassica napus, improving seed-weight and oil-content is the main breeding goal. In order to better understand the regulation mechanism of miRNAs during seed-weight formation and oil-content accumulation in B. napus, in this study, a high-throughput sequencing technology was used to profile miRNAs expression of Brassica napus immature seeds from one to six weeks after flowering. A total of 1,276 miRNAs, including 1,248 novel and 28 known miRNAs, were obtained from both the high-seed-weight with low-oil-content RNA pool (S03) and the low-seed-weight with high-oil-content RNA pool (S04). Analysis of their expression profiles disclosed that 300 novel and two known miRNAs were differentially expressed between S03 and S04. For degradome analysis, 57 genes with 64 degradation sites were predicted to be targeted for degradation by these miRNAs. Further bioinformatics analysis indicated that these differentially expressed miRNAs might participate in regulation of myriad cellular and molecular processes, during seed development and oil synthesis. Finally, 6 target genes with potential roles in regulation of seed development and 9 other targets in seed oil synthesis, were further confirmed as candidate genes from small RNA and degradome sequencing.

V. AM, Solanke AU, Singh NK, Sharma TR, Gaikwad K. Genome-wide discovery of tissue-specific miRNAs in clusterbean (Cyamopsis tetragonoloba) indicates their association with galactomannan biosynthesis

Owing to the presence of 80% soluble dietary fibre, high protein content and high value gum, clusterbean (Cyamopsis tetragonoloba) has recently emerged as an economically important legume. The developing clusterbean seeds accumulate 90% galactomannans in the endosperm and, therefore, can be used as a model crop to understand galactomannan biosynthesis and its regulation. miRNAs are tiny master regulators of their corresponding target genes, resulting in variations in the amounts of their metabolic end products. To understand the role of these regulators in galactomannan biosynthesis regulation, small RNA libraries were prepared and sequenced from five tissues of clusterbean genotype RGC-936, and miRanalyzer and DSAP programs were used to identify conserved miRNAs and novel small RNAs. A total of 187 known and 171 novel miRNAs were found to be differentially expressed, of which 10 miRNAs were validated. A complicated network topology and 35% sharing of the target mRNAs between known and novel miRNAs suggest random evolution of novel miRNAs. The gene ontology (GO) annotation of potential target genes revealed the genes coding for signalling and carbohydrate metabolism (50.10%), kinases and other enzymes (20.75%), transcription factors (10.20%), transporters (8.35%) and other targets (10.6%). Two novel unigenes were annotated as ManS (mannosyltransferase/mannan synthase) and UGE (UDP- D-glucose 4-epimerase) and validated as targets for three novel miRNAs, that is Ct-miR3130, Ct-miR3135 and Ct-miR3157. Our findings reveal that these novel miRNAs could play an important role in the regulation of the galactomannan pathway in C. tetragonoloba and possibly other galactomannan-producing species.

Unravelling miRNA regulation in yield of rice (Oryza sativa) based on differential network model

Rice (Oryza sativa L.) is one of the essential staple food crops and tillering, panicle branching and grain filling are three important traits determining the grain yield. Although miRNAs have been reported being regulating yield, no study has systematically investigated how miRNAs differentially function in high and low yield rice, in particular at a network level. This abundance of data from high-throughput sequencing provides an effective solution for systematic identification of regulatory miRNAs using developed algorithms in plants. We here present a novel algorithm, Gene Co-expression Network differential edge-like transformation (GRN-DET), which can identify key regulatory miRNAs in plant development. Based on the small RNA and RNA-seq data, miRNA-gene-TF co-regulation networks were constructed for yield of rice. Using GRN-DET, the key regulatory miRNAs for rice yield were characterized by the differential expression variances of miRNAs and co-variances of miRNA-mRNA, including osa-miR171 and osa-miR1432. Phytohormone cross-talks (auxin and brassinosteroid) were also revealed by these co-expression networks for the yield of rice.

Comparative analysis of miRNA expression profiles in transgenic and non-transgenic rice using miRNA-Seq

Safety assessment for genetically modified organisms (GMOs) is required before their release. To date, miRNAs that play important roles in eukaryotic gene regulation have not been considered in the current assessment system. In this study, we identified 6 independent Bt and EPSPS GM rice lines using PCR and immune strip. We analyzed the expression levels of Cry1Ac and EPSPS using quantitative real-time PCR and western blot. Further, miRNAs from the developing seeds of the 6 GM rice lines and the wild-type line were investigated using deep sequencing and bioinformatic approaches. Although these GM lines have different types of integration sites, copy numbers, and levels of gene expression, 21 differentially expressed miRNAs have been found compared to wild type. There is no correlation between transgenic protein expression level and the quantity of differentially expressed miRNAs. This study provides useful data about the miRNA composition of GM plants, and it might be helpful for future risk assessments of miRNA-based GM plants.

Comprehensive analysis of NAC transcription factors and their expression during fruit spine development in cucumber (Cucumis sativus L.)

The cucumber (Cucumis sativus L.) is an important vegetable crop worldwide, and fruit trichomes or spines are an important trait for external fruit quality. The mechanisms underlying spine formation are not well understood, but the plant-specific NAC family of transcription factors may play important roles in fruit spine initiation and development. In this study, we conducted a genome-wide survey and identified 91 NAC gene homologs in the cucumber genome. Clustering analysis classified these genes into six subfamilies; each contained a varying number of NAC family members with a similar intron-exon structure and conserved motifs. Quantitative real-time PCR analysis revealed tissue-specific expression patterns of these genes, including 10 and 12 that exhibited preferential expression in the stem and fruit, respectively. Thirteen of the 91 NAC genes showed higher expression in the wild-type plant than in its near-isogenic trichome mutant, suggesting their important roles in fruit spine development. Exogenous application of four plant hormones promoted spine formation and increased spine density on the cucumber fruits; several NAC genes showed differential expression over time in response to phytohormone treatments on cucumber fruit, implying their essential roles in fruit-trichome development. Among the NAC genes identified, 12 were found to be targets of 13 known cucumber micro-RNAs. Collectively, these findings provide a useful resource for further analysis of the interactions between NAC genes and genes underlying trichome organogenesis and development during fruit spine development in cucumber.

microRNAs participate in gene expression regulation and phytohormone cross-talk in barley embryo during seed development and germination

BACKGROUND

Small RNA and degradome sequencing have identified a large number of miRNA-target pairs in plant seeds. However, detailed spatial and temporal studies of miRNA-mediated regulation, which can reflect links between seed development and germination are still lacking.

RESULTS

In this study, we extended our investigation on miRNAs-involved gene regulation by a combined analysis of seed maturation and germination in barley. Through bioinformatics analysis of small RNA sequencing data, a total of 1324 known miRNA families and 448 novel miRNA candidates were identified. Of those, 16 known miRNAs with 40 target genes, and three novel miRNAs with four target genes were confirmed based on degradome sequencing data. Conserved miRNA families such as miR156, miR168, miR166, miR167, and miR894 were highly expressed in embryos of developing and germinating seeds. A barley-specific miRNA, miR5071, which was predicted to target an OsMLA10-like gene, accumulated at a high level, suggesting its involvement in defence response during these two developmental stages. Based on target prediction and Kyoto Encyclopedia of Genes and Genomes analysis of putative targets, nine highly expressed miRNAs were found to be related to phytohormone signalling and hormone cross-talk. Northern blot and qRT-PCR analysis showed that these miRNAs displayed differential expression patterns during seed development and germination, indicating their different roles in hormone signalling pathways. In addition, we showed that miR393 affected seed development through targeting two genes encoding the auxin receptors TIR1/AFBs in barley, as over-expression of miR393 led to an increased length-width ratio of seeds, whereas target mimic (MIM393)-mediated inhibition of its activity decreased the 1000-grain weight of seeds. Furthermore, the expression of auxin-responsive genes, abscisic acid- and gibberellic acid-related genes was altered in miR393 misexpression lines during germination and early seedling growth.

CONCLUSIONS

Our work indicates that miRNA-target pairs participate in gene expression regulation and hormone interaction in barley embryo and provides evidence that miR393-mediated auxin response regulation affects grain development and influences gibberellic acid and abscisic acid homeostasis during germination.

De novo assembly and annotation of the Zhe-Maidong (Ophiopogon japonicus (L.f.) Ker-Gawl) transcriptome in different growth stages

Zhe-Maidong (Ophiopogon japonicus (L.f.) Ker-Gawl) is a traditional medicinal herb in the family Liliaceae that has significant pharmacological effects on immunity and cardiovascular disease. In this study, three different growth stages of Zhe-Maidong were investigated using RNA-seq, and a total of 16.4 Gb of raw data was obtained. After filtering and assembling, 96,738 unigenes with an average length of 605.3 bp were ultimately generated. A total of 77,300 unigenes were annotated using information from five databases, including the NT, NR, SwissProt, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) databases. Additionally, the mechanisms of flavonoid, saponin and polysaccharide biosynthesis and of accumulation at different stages of tuber development were also characterized. From the first to third years, the contents of flavonoids, saponins and polysaccharides all increased, whereas the expression levels of related genes decreased in the flavonoid and saponin pathways and first increased and then decreased in the polysaccharide pathway. The results of this study provide the most comprehensive expressed sequence resource for Zhe-Maidong and will expand the available O. japonicus gene library and facilitate further genome-wide research and analyses of this species.

The Sequences of 1504 Mutants in the Model Rice Variety Kitaake Facilitate Rapid Functional Genomic Studies

The availability of a whole-genome sequenced mutant population and the cataloging of mutations of each line at a single-nucleotide resolution facilitate functional genomic analysis. To this end, we generated and sequenced a fast-neutron-induced mutant population in the model rice cultivar Kitaake (Oryza sativa ssp japonica), which completes its life cycle in 9 weeks. We sequenced 1504 mutant lines at 45-fold coverage and identified 91,513 mutations affecting 32,307 genes, i.e., 58% of all rice genes. We detected an average of 61 mutations per line. Mutation types include single-base substitutions, deletions, insertions, inversions, translocations, and tandem duplications. We observed a high proportion of loss-of-function mutations. We identified an inversion affecting a single gene as the causative mutation for the short-grain phenotype in one mutant line. This result reveals the usefulness of the resource for efficient, cost-effective identification of genes conferring specific phenotypes. To facilitate public access to this genetic resource, we established an open access database called KitBase that provides access to sequence data and seed stocks. This population complements other available mutant collections and gene-editing technologies. This work demonstrates how inexpensive next-generation sequencing can be applied to generate a high-density catalog of mutations.

The pivotal role of small non-coding RNAs in the regulation of seed development

Seeds represent a crucial stage of the seed plants life cycle. It is during seed development that the foundations of the future plant body, and the ability to give rise to a new plant capable of growing under sometimes adverse environmental conditions, are established. Small non-coding RNAs are major regulators of gene expression both at the post-transcriptional and transcriptional levels and, not surprisingly, these elements play major roles in seed development and germination. We review here the current knowledge about small RNA expression and functions in seed development, going from the morphogenesis phase comprehending embryo development and patterning, to the several steps of the maturation phase, ending in the transition to the germination. A special focus is given to the small RNAs for which functional studies have been conducted and their participation in regulatory networks operating in seeds. Many challenges remain ahead for dissecting the complex small RNA landscape in seeds, but this is a highly relevant issue in plant biology and advances in this area will most certainly impact plant breeding.

Detection of dietetically absorbed maize-derived microRNAs in pigs

MicroRNAs are a class of small RNAs that are important in post-transcriptional gene regulation in animals and plants. These single-stranded molecules are widely distributed in organisms and influence fundamental biological processes. Interestingly, recent studies have reported that diet-derived plant miRNAs could regulate mammalian gene expression, and these studies have broadened our view of cross-kingdom communication. In the present study, we evaluated miRNA levels in cooked maize-containing chow diets, and found that plant miRNAs were resistant to the harsh cooking conditions to a certain extent. After feeding fresh maize to pigs (7 days), maize-derived miRNAs could be detected in porcine tissues and serum, and the authenticity of these plant miRNAs was confirmed by using oxidization reactions. Furthermore, in vivo and in vitro experiments demonstrated that dietary maize miRNAs could cross the gastrointestinal tract and enter the porcine bloodstream. In the porcine cells, we found that plant miRNAs are very likely to specifically target their endogenous porcine mRNAs and influence gene expression in a fashion similar to that of mammalian miRNAs. Our results indicate that maize-derived miRNAs can cross the gastrointestinal tract and present in pigs, and these exogenous miRNAs have the potential to regulate mammalian gene expression.

OsGRF4 controls grain shape, panicle length and seed shattering in rice

Traits such as grain shape, panicle length and seed shattering, play important roles in grain yield and harvest. In this study, the cloning and functional analysis of PANICLE TRAITS 2 (PT2), a novel gene from the Indica rice Chuandali (CDL), is reported. PT2 is synonymous with Growth-Regulating Factor 4 (OsGRF4), which encodes a growth-regulating factor that positively regulates grain shape and panicle length and negatively regulates seed shattering. Higher expression of OsGRF4 is correlated with larger grain, longer panicle and lower seed shattering. A unique OsGRF4 mutation, which occurs at the OsmiRNA396 target site of OsGRF4, seems to be associated with high levels of OsGRF4 expression, and results in phenotypic difference. Further research showed that OsGRF4 regulated two cytokinin dehydrogenase precursor genes (CKX5 and CKX1) resulting in increased cytokinin levels, which might affect the panicle traits. High storage capacity and moderate seed shattering of OsGRF4 may be useful in high-yield breeding and mechanized harvesting of rice. Our findings provide additional insight into the molecular basis of panicle growth.

MicroRNA Signatures of Drought Signaling in Rice Root

Background

Drought stress is one of the most important abiotic stresses and the main constraint to rice agriculture. MicroRNA-mediated post-transcriptional gene regulation is one of the ways to establish drought stress tolerance in plants. MiRNAs are 20–24-nt regulatory RNAs that play an important role in regulating plant gene expression upon exposure to biotic and abiotic stresses.

Methodology/Principal Findings

In this study, we applied a partial root drying system as well as a complete root drying system to identify miRNAs involved in conditions of drought stress, drought signaling and wet signaling using high-throughput sequencing. To this end, we produced four small RNA libraries: (1) fully-watered (WW), (2) fully-droughted (WD), and split-root systems where (3) one-half was well watered (SpWW) and (4) the other half was water-deprived (SpWD). Our analysis revealed 10,671 and 783 unique known and novel miRNA reads in all libraries, respectively. We identified, 65 (52 known + 13 novel), 72 (61 known + 11 novel) and 51 (38 known + 13 novel) miRNAs that showed differential expression under conditions of drought stress, drought signaling and wet signaling, respectively. The results of quantitative real-time PCR showed expression patterns similar to the high-throughput sequencing results. Furthermore, our target prediction led to the identification of 244, 341 and 239 unique target genes for drought-stress-, drought-signaling- and wet-signaling-responsive miRNAs, respectively.

Conclusions/Significance

Our results suggest that miRNAs that are responsive under different conditions could play different roles in the regulation of abscisic acid signaling, calcium signaling, detoxification and lateral root formation.

Quantitative phosphoproteomic analysis of early seed development in rice (Oryza sativa L.)

Rice (Oryza sativa L.) seed serves as a major food source for over half of the global population. Though it has been long recognized that phosphorylation plays an essential role in rice seed development, the phosphorylation events and dynamics in this process remain largely unknown so far. Here, we report the first large scale identification of rice seed phosphoproteins and phosphosites by using a quantitative phosphoproteomic approach. Thorough proteomic studies in pistils and seeds at 3, 7 days after pollination resulted in the successful identification of 3885, 4313 and 4135 phosphopeptides respectively. A total of 2487 proteins were differentially phosphorylated among the three stages, including Kip related protein 1, Rice basic leucine zipper factor 1, Rice prolamin box binding factor and numerous other master regulators of rice seed development. Moreover, differentially phosphorylated proteins may be extensively involved in the biosynthesis and signaling pathways of phytohormones such as auxin, gibberellin, abscisic acid and brassinosteroid. Our results strongly indicated that protein phosphorylation is a key mechanism regulating cell proliferation and enlargement, phytohormone biosynthesis and signaling, grain filling and grain quality during rice seed development. Overall, the current study enhanced our understanding of the rice phosphoproteome and shed novel insight into the regulatory mechanism of rice seed development.

Small RNA and degradome sequencing reveals microRNAs and their targets involved in tomato pedicel abscission

We constructed small RNA and degradome sequencing libraries to identify miRNAs and targets involved in tomato pedicel abscission, and confirmed their roles via quantitative real-time PCR. MicroRNAs (miRNAs) are endogenous small RNAs which play crucial negatively regulatory roles at both the transcriptional and post-transcriptional levels in plants; however, limited knowledge is available on the expression profiles of miRNAs and their target genes during tomato pedicel abscission. Taking advantage of small RNA (sRNA) and degradome sequencing technology, a total of 56 known and 11 novel candidate miRNAs targeting 223 mRNA genes were confirmed during pedicel abscission. Gene ontology annotation and KEGG pathway analysis showed that these target genes were significantly enriched in intracellular, membrane-bounded organelle-related biological processes as well as in metabolic, plant-pathogen interaction and hormone signaling pathways. We screened 17 miRNA/target pairs for further analysis and performed quantitative real-time PCR to identify the roles. Cluster analysis of selected miRNAs revealed that the expression profiles of miRNAs varied in different stages of abscission and could be impacted by ethylene treatment. In the present study, the correlations between miRNAs and targets suggested a complex regulatory network of miRNA-mediated target interaction during pedicel abscission. Additionally, the expression profiles of miRNAs and their targets changed by ethylene might be a considerable reason why ethylene promotes pedicel abscission. Our study provides new insights into the expression and regulatory profiles of miRNAs during tomato pedicel abscission.

Temporal small RNA transcriptome profiling unraveled partitioned miRNA expression in developing maize endosperms between reciprocal crosses

In angiosperms, the endosperm nurtures the embryo and provides nutrients for seed germination. To identify the expression pattern of small interfering RNA in the developing maize endosperm, we have performed high-throughput small RNA transcriptome sequencing of kernels at 0, 3, and 5 days after pollination (DAP) and endosperms at 7, 10, and 15 DAP using B73 and Mo17 reciprocal crosses in previous study. Here, we further explored these small RNA-seq data to investigate the potential roles of miRNAs in regulating the gene expression process. In total, 57 conserved miRNAs and 18 novel miRNAs were observed highly expressed in maize endosperm. Temporal expression profiling indicated that these miRNAs exhibited dynamic and partitioned expression patterns at different developmental stages between maize reciprocal crosses, and quantitative RT-PCR results further confirmed our observation. In addition, we found a subset of distinct tandem miRNAs are generated from a single stem-loop structure in maize that might be conserved in monocots. Furthermore, a SNP variation of Zma-miR408-5p at 11th base position was characterized between B73 and Mo17 which might lead to completely different functions in repressing targets. More interestingly, Zma-miR408-5p exhibited B73-biased expression pattern in the B73 and Mo17 reciprocal hybrid endosperms at 7, 10, and 15 DAP according to the reads abundance with SNPs and CAPS experiment. Together, this study suggests that miRNA plays a crucial role in regulating endosperm development, and exhibited distinct expression patterns in developing endosperm between maize reciprocal crosses.

Identification of four functionally important microRNA families with contrasting differential expression profiles between drought-tolerant and susceptible rice leaf at vegetative stage

Background

Developing drought-tolerant rice varieties with higher yield under water stressed conditions provides a viable solution to serious yield-reduction impact of drought. Understanding the molecular regulation of this polygenic trait is crucial for the eventual success of rice molecular breeding programmes. microRNAs have received tremendous attention recently due to its importance in negative regulation. In plants, apart from regulating developmental and physiological processes, microRNAs have also been associated with different biotic and abiotic stresses. Hence here we chose to analyze the differential expression profiles of microRNAs in three drought treated rice varieties: Vandana (drought-tolerant), Aday Sel (drought-tolerant) and IR64 (drought-susceptible) in greenhouse conditions via high-throughput sequencing.

Results

Twenty-six novel microRNA candidates involved in the regulation of diverse biological processes were identified based on the detection of miRNA*. Out of their 110 predicted targets, we confirmed 16 targets from 5 novel microRNA candidates. In the differential expression analysis, mature microRNA members from 49 families of known Oryza sativa microRNA were differentially expressed in leaf and stem respectively with over 28 families having at least a similar mature microRNA member commonly found to be differentially expressed between both tissues. Via the sequence profiling data of leaf samples, we identified osa-miR397a/b, osa-miR398b, osa-miR408-5p and osa-miR528-5p as being down-regulated in two drought-tolerant rice varieties and up-regulated in the drought-susceptible variety. These microRNAs are known to be involved in regulating starch metabolism, antioxidant defence, respiration and photosynthesis. A wide range of biological processes were found to be regulated by the target genes of all the identified differentially expressed microRNAs between both tissues, namely root development (5.3–5.7 %), cell transport (13.2–18.4 %), response to stress (10.5–11.3 %), lignin catabolic process (3.8–5.3 %), metabolic processes (32.1–39.5 %), oxidation-reduction process (9.4–13.2 %) and DNA replication (5.7–7.9 %). The predicted target genes of osa-miR166e-3p, osa-miR166h-5p*, osa-miR169r-3p* and osa-miR397a/b were found to be annotated to several of the aforementioned biological processes.

Conclusions

The experimental design of this study, which features rice varieties with different drought tolerance and tissue specificity (leaf and stem), has provided new microRNA profiling information. The potentially regulatory importance of the microRNA genes mentioned above and their target genes would require further functional analyses.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1851-3) contains supplementary material, which is available to authorized users.

Identification and characterization of microRNAs at different flowering developmental stages in moso bamboo (Phyllostachys edulis) by high-throughput sequencing

Researching moso bamboo flowering has been difficult because of its unknown flowering interval and the rarity of florescent samples. To identify microRNAs (miRNAs) and study their expression patterns during the flower developmental process of moso bamboo, small RNAs from non-flowering leaves and four flower developmental periods were sequenced using Illumina technology. In total, 409 known miRNAs and 492 differentially expressed novel miRNAs were identified in moso bamboo. Of the known miRNAs that were differentially expressed between non-flowering and flowering samples, 64 were predicted to have a total of 308 targets. Among the miRNAs, seven known and five novel miRNAs were selected, as were four of their target genes, and their expression profiles were validated using qRT-PCR. The results indicated that the miRNA expression levels were negatively correlated with those of their targets. The research comprehensively revealed that the differentially expressed miRNAs and their targets participated in diverse biological pathways and played significant regulatory roles in moso bamboo flowering. The data provide a significant resource for understanding the molecular mechanisms in moso bamboo flowering and senescence, and serve as the primary foundation for further studies on metabolic regulatory networks that involve miRNAs.

Stress responsive miRNAs and isomiRs in cereals

Abiotic and biotic stress conditions are vital determinants in the production of cereals, the major caloric source in human nutrition. Small RNAs, miRNAs and isomiRs are central to post-transcriptional regulation of gene expression in a variety of cellular processes including development and stress responses. Several miRNAs have been identified using new technologies and have roles in stress responses in plants, including cereals. The overall knowledge about the cereal miRNA repertoire, as well as an understanding of complex miRNA mediated mechanisms of target regulation in response to stress conditions, is far from complete. Ongoing efforts that add to our understanding of complex miRNA machinery have implications in plant response to stress conditions. Additionally, sequence variants of miRNAs (isomiRNAs or isomiRs), regulation of their expression through dissection of upstream regulatory elements, the role of Processing-bodies (P-bodies) in miRNA exerted gene regulation and yet unveiled organellar plant miRNAs are newly emerging topics, which will contribute to the elucidation of the miRNA machinery and its role in cereal tolerance against abiotic and biotic stresses.

Identification and Characterization of microRNAs during Maize Grain Filling

The grain filling rate is closely associated with final grain yield of maize during the period of maize grain filling. To identify the key microRNAs (miRNAs) and miRNA-dependent gene regulation networks of grain filling in maize, a deep-sequencing technique was used to research the dynamic expression patternsof miRNAs at four distinct developmental grain filling stages in Zhengdan 958, which is an elite hybrid and cultivated widely in China. The sequencing result showed that the expression amount of almost all miRNAs was changing with the development of the grain filling and formed in seven groups. After normalization, 77 conserved miRNAs and 74 novel miRNAs were co-detected in these four samples. Eighty-one out of 162 targets of the conserved miRNAs belonged to transcriptional regulation (81, 50%), followed by oxidoreductase activity (18, 11%), signal transduction (16, 10%) and development (15, 9%). The result showed that miRNA 156, 393, 396 and 397, with their respective targets, might play key roles in the grain filling rate by regulating maize growth, development and environment stress response. The result also offered novel insights into the dynamic change of miRNAs during the developing process of maize kernels and assistedin the understanding of how miRNAs are functioning about the grain filling rate.

MicroRNA-based biotechnology for plant improvement

MicroRNAs (miRNAs) are an extensive class of newly discovered endogenous small RNAs, which negatively regulate gene expression at the post-transcription levels. As the application of next-generation deep sequencing and advanced bioinformatics, the miRNA-related study has been expended to non-model plant species and the number of identified miRNAs has dramatically increased in the past years. miRNAs play a critical role in almost all biological and metabolic processes, and provide a unique strategy for plant improvement. Here, we first briefly review the discovery, history, and biogenesis of miRNAs, then focus more on the application of miRNAs on plant breeding and the future directions. Increased plant biomass through controlling plant development and phase change has been one achievement for miRNA-based biotechnology; plant tolerance to abiotic and biotic stress was also significantly enhanced by regulating the expression of an individual miRNA. Both endogenous and artificial miRNAs may serve as important tools for plant improvement.

Expanding frontiers in plant transcriptomics in aid of functional genomics and molecular breeding

The transcript pool of a plant part, under any given condition, is a collection of mRNAs that will pave the way for a biochemical reaction of the plant to stimuli. Over the past decades, transcriptome study has advanced from Northern blotting to RNA sequencing (RNA-seq), through other techniques, of which real-time quantitative polymerase chain reaction (PCR) and microarray are the most significant ones. The questions being addressed by such studies have also matured from a solitary process to expression atlas and marker-assisted genetic enhancement. Not only genes and their networks involved in various developmental processes of plant parts have been elucidated, but also stress tolerant genes have been highlighted. The transcriptome of a plant with altered expression of a target gene has given information about the downstream genes. Marker information has been used for breeding improved varieties. Fortunately, the data generated by transcriptome analysis has been made freely available for ample utilization and comparison. The review discusses this wide variety of transcriptome data being generated in plants, which includes developmental stages, abiotic and biotic stress, effect of altered gene expression, as well as comparative transcriptomics, with a special emphasis on microarray and RNA-seq. Such data can be used to determine the regulatory gene networks, which can subsequently be utilized for generating improved plant varieties.

Differentially expressed microRNA cohorts in seed development may contribute to poor grain filling of inferior spikelets in rice

BACKGROUND

The inferior spikelets are defined to be those at portions where the grains receive less photosynthetic products during the seed development. The typical inferior spikelets are physically located on the proximal secondary branches in a rice panicle and traditionally characterized by a later flowering time and a slower grain-filling rate, compared to those so-called superior spikelets. Grains produced on the inferior spikelets are consequently under-developed and lighter in weight than those formed on the superior spikelets. MicroRNAs (miRNAs) are recognized as key players in regulating plant development through post-transcriptional gene regulations. We previously presented the evidence that miRNAs may influence grain-filling rate and played a role in determining the grain weight and yield in rice.

RESULTS

In this study, further analyses of the expressed small RNAs in superior and inferior spikelets were conducted at five distinct developmental stages of grain development. Totally, 457 known miRNAs and 13 novel miRNAs were analyzed, showing a differential expression of 141 known miRNAs between superior and inferior spikelets with higher expression levels of most miRNAs associated with the superior than the inferior spikelets during the early stage of grain filling. Genes targeted by those differentially expressed miRNAs (i.e. miR156, miR164, miR167, miR397, miR1861, and miR1867) were recognized to play roles in multiple developmental and signaling pathways related to plant hormone homeostasis and starch accumulation.

CONCLUSIONS

Our data established a complicated link between miRNA dynamics and the traditional role of hormones in grain filling and development, providing new insights into the widely accepted concepts of the so-called superior and inferior spikelets in rice production.

Identification of 23 novel conserved microRNAs in three rice cultivars

Plant microRNAs (miRNAs) play important roles as modulators of gene expression at the post-transcriptional level. Previous studies have shown that high-throughput sequencing is a powerful tool for the identification of miRNAs, and it is believed that many more miRNAs remain to be discovered. Here, we found 23 novel conserved miRNAs from three rice cultivars by high-throughput sequencing and further identified these through subsequent cloning and quantitative real-time polymerase chain reaction (qPCR). Eight of these novel miRNAs were detected with significant signals in the three rice cultivars by northern blotting assays. The quantitative analysis of their expression profiles showed that most of these miRNAs were perfectly or imperfectly negatively correlated with their target genes, which suggests that these miRNAs may play important roles during rice development. This is the first genome-wide investigation of miRNAs from different rice cultivars, and the data obtained expand the known rice miRNA inventory and provide further information about the regulatory roles played by miRNAs in rice development.

MicroRNAs responding to southern rice black-streaked dwarf virus infection and their target genes associated with symptom development in rice

Southern rice black-streaked dwarf virus (SRBSDV) is a recently emerged rice virus that has spread across Asia. This devastating virus causes rice plants to produce a variety of symptoms during different growth stages. MicroRNAs (miRNAs) comprise a large group of 21-24-nt RNA molecules that are important regulators of plant development processes and stress responses. In this study, we used microarray profiling to investigate rice miRNAs responding to SRBSDV infection at 3, 9, 15, and 20 days post-inoculation (dpi). Expression levels of 56 miRNAs were altered in SRBSDV-infected rice plants, with these changes classified into eight different regulation patterns according to their temporal expression dynamics. Fourteen miRNAs belonging to six families (miR164, R396, R530, R1846, R1858, and R2097) were significantly regulated at 20 dpi. We used RT-qPCR to search for expression level correlations between members of these families and their putative targets at 3, 9, and 15 dpi. Some members of the miR164, R396, R530, and R1846 families were found to be positively or negatively correlated with their respective targets during 3-15 days after SRBSDV infection, whereas in more cases the rice miRNAs were not in correlation with their targets along the post-inoculation period, suggesting that some additional factors may be involved in rice miRNA-target interactions. The reported functions of rice genes targeted by the miR164, R396, R530, R1846, and R1858 families indicated that these genes are associated with symptom development. These results provide insights into miRNA-mediated SRBSDV-rice interactions.

Whole-genome discovery of miRNAs and their targets in wheat (Triticum aestivum L.)

BACKGROUND

MicroRNAs (miRNAs) are small, non-coding RNAs playing essential roles in plant growth, development, and stress responses. Sequencing of small RNAs is a starting point for understanding their number, diversity, expression and possible roles in plants.

RESULTS

In this study, we conducted a genome-wide survey of wheat miRNAs from 11 tissues, characterizing a total of 323 novel miRNAs belonging to 276 families in wheat. A miRNA conservation analysis identified 191 wheat-specific miRNAs, 2 monocot-specific miRNAs, and 30 wheat-specific variants from 9 highly conserved miRNA families. To understand possible roles of wheat miRNAs, we determined 524 potential targets for 124 miRNA families through degradome sequencing, and cleavage of a subset of them was validated via 5' RACE. Based on the genome-wide identification and characterization of miRNAs and their associated target genes, we further identified 64 miRNAs preferentially expressing in developing or germinating grains, which could play important roles in grain development.

CONCLUSION

We discovered 323 wheat novel miRNAs and 524 target genes for 124 miRNA families in a genome-wide level, and our data will serve as a foundation for future research into the functional roles of miRNAs in wheat.

Identification of abiotic stress miRNA transcription factor binding motifs (TFBMs) in rice

Plant growth and yield are affected by many abiotic stresses like salinity, drought, cold and heavy metal; these stresses trigger up and down-regulate several genes through various transcription factors (TFs). Transcription factor binding motifs (TFBMs), located in the upstream region of the genes, associate with TFs to regulate the gene expression. Many factors, including the activation of miRNAs, which are encoded by genes having independent transcription units, regulate the gene expression. TFBMs in the regulatory region of miRNA sequences influence the miRNA expression, which in turn influences the expression of other genes in the cell. However, the current level of information available on TFBMs of miRNA involved in abiotic stress related defense pathway(s) is limited and in-depth studies in this direction may lead to a better understanding of their role in expression and regulation of defense responses in plants. In this study, various aspects related to genomic positions of pre-miRNA, prediction of TSS and TATA box positions and identification of known, unique motifs at regulatory regions of all the reported miRNAs of rice associated with different abiotic stresses are discussed. Sixteen motifs were identified in this study, of which nine are known cis-regulatory elements associated with various stresses, two strong motifs, (CGCCGCCG, CGGCGGCG) and five unique motifs which might play a vital role in the regulation of abiotic stresses related miRNA genes. Common motifs shared by miRNAs that are involved in more than one abiotic stresses were also identified. The motifs identified in this study will be a resource for further functional validation.

Evaluation of three RT-qPCR-based miRNA detection methods using seven rice miRNAs

Three frequently-used reverse transcription-quantitative polymerase chain reaction (RT-qPCR)-based miRNA detection methods, stem-loop RT-qPCR, poly(A)-tailing RT-qPCR, and miQPCR, were evaluated using seven selected rice miRNAs. The results revealed that miRNA abundance and sequence characteristics can affect capability of detection. The stem-loop amplification technique detected highly and moderately abundant miRNAs. The poly(A)-tailing method detected both highly abundant and sparsely present miRNAs, but failed to detect miRNAs with a hairpin structure. Only a few miRNAs were detectable by the miQPCR method. We suggest that a combination of methods should be used for reliable quantitative investigation of miRNAs.

Functional genomics based understanding of rice endosperm development

Seed development, especially the relevant regulatory mechanism and genetic network are of fundamental scientific interest. Seed development consists of the development of embryo and endosperm; and endosperm development of rice (model species of monocots) is closely related to grain yield and quality. Recent genetic studies, together with other approaches, including transcriptome and proteomics analysis, high-throughput sequencing (RNA-seq, ChIP-seq), revealed the crucial roles of genetic and epigenetic controls in rice endosperm development. Here we summarize and update the genetic networks involved in the regulation of endosperm initiation, cell cycle regulation, aleurone layer specification, starch synthesis, storage protein accumulation and endosperm size, and the interactions between embryo and endosperm.

Identification and expression analysis of microRNAs at the grain filling stage in rice(Oryza sativa L.)via deep sequencing

MicroRNAs (miRNAs) have been shown to play crucial roles in the regulation of plant development. In this study, high-throughput RNA-sequencing technology was used to identify novel miRNAs, and to reveal miRNAs expression patterns at different developmental stages during rice (Oryza sativa L.) grain filling. A total of 434 known miRNAs (380, 402, 390 and 392 at 5, 7, 12 and 17 days after fertilization, respectively.) were obtained from rice grain. The expression profiles of these identified miRNAs were analyzed and the results showed that 161 known miRNAs were differentially expressed during grain development, a high proportion of which were up-regulated from 5 to 7 days after fertilization. In addition, sixty novel miRNAs were identified, and five of these were further validated experimentally. Additional analysis showed that the predicted targets of the differentially expressed miRNAs may participate in signal transduction, carbohydrate and nitrogen metabolism, the response to stimuli and epigenetic regulation. In this study, differences were revealed in the composition and expression profiles of miRNAs among individual developmental stages during the rice grain filling process, and miRNA editing events were also observed, analyzed and validated during this process. The results provide novel insight into the dynamic profiles of miRNAs in developing rice grain and contribute to the understanding of the regulatory roles of miRNAs in grain filling.

Characterization and expression patterns of microRNAs involved in rice grain filling

MicroRNAs (miRNAs) are upstream gene regulators of plant development and hormone homeostasis through their directed cleavage or translational repression of the target mRNAs, which may play crucial roles in rice grain filling and determining the final grain weight and yield. In this study, high-throughput sequencing was performed to survey the dynamic expressions of miRNAs and their corresponding target genes at five distinct developmental stages of grain filling. In total, 445 known miRNAs and 45 novel miRNAs were detected with most of them expressed in a developmental stage dependent manner, and the majority of known miRNAs, which increased gradually with rice grain filling, showed negatively related to the grain filling rate. Detailed expressional comparisons revealed a clear negative correlation between most miRNAs and their target genes. It was found that specific miRNA cohorts are expressed in a developmental stage dependent manner during grain filling and the known functions of these miRNAs are involved in plant hormone homeostasis and starch accumulation, indicating that the expression dynamics of these miRNAs might play key roles in regulating rice grain filling.