Characterization of small RNAs and their target genes in wheat seedlings using sequencing-based approaches

Genome-wide expression analysis of novel heat-responsive microRNAs and their targets in contrasting wheat genotypes at reproductive stage under terminal heat stress

Introduction

Heat stress at terminal stage of wheat is critical and leads to huge yield losses worldwide. microRNAs (miRNAs) play significant regulatory roles in gene expression associated with abiotic and biotic stress at the post-transcriptional level.

Methods

In the present study, we carried out a comparative analysis of miRNAs and their targets in flag leaves as well as developing seeds of heat tolerant (RAJ3765) and heat susceptible (HUW510) wheat genotypes under heat stress and normal conditions using small RNA and degradome sequencing.

Results and discussion

A total of 84 conserved miRNAs belonging to 35 miRNA families and 93 novel miRNAs were identified in the 8 libraries. Tae-miR9672a-3p, tae-miR9774, tae-miR9669-5p, and tae-miR5048-5p showed the highest expression under heat stress. Tae-miR9775, tae-miR9662b-3p, tae-miR1120a, tae-miR5084, tae-miR1122a, tae-miR5085, tae-miR1118, tae-miR1130a, tae-miR9678-3p, tae-miR7757-5p, tae-miR9668-5p, tae-miR5050, tae-miR9652-5p, and tae-miR9679-5p were expressed only in the tolerant genotype, indicating their role in heat tolerance. Comparison between heat-treated and control groups revealed that 146 known and 57 novel miRNAs were differentially expressed in the various tissues. Eight degradome libraries sequence identified 457 targets of the differentially expressed miRNAs. Functional analysis of the targets indicated their involvement in photosynthesis, spliceosome, biosynthesis of nucleotide sugars and protein processing in the endoplasmic reticulum, arginine and proline metabolism and endocytosis.

Conclusion

This study increases the number of identified and novel miRNAs along with their roles involved in heat stress response in contrasting genotypes at two developing stages of wheat.

Differential regulation of miRNAs involved in the susceptible and resistance responses of wheat cultivars to wheat streak mosaic virus and Triticum mosaic virus

BACKGROUND

Wheat streak mosaic virus (WSMV) and Triticum mosaic virus (TriMV) are components of the wheat streak mosaic virus disease complex in the Great Plains region of the U.S.A. and elsewhere. Co-infection of wheat with WSMV and TriMV causes synergistic interaction with more severe disease symptoms compared to single infections. Plants are equipped with multiple antiviral mechanisms, of which regulation of microRNAs (miRNAs) is a potentially effective constituent. In this investigation, we have analyzed the total and relative expression of miRNA transcriptome in two wheat cultivars, Arapahoe (susceptible) and Mace (temperature-sensitive-resistant), that were mock-inoculated or inoculated with WSMV, TriMV, or both at 18 °C and 27 °C.

RESULTS

Our results showed that the most abundant miRNA family among all the treatments was miRNA166, followed by 159a and 168a, although the order of the latter two changed depending on the infections. When comparing infected and control groups, twenty miRNAs showed significant upregulation, while eight miRNAs were significantly downregulated. Among them, miRNAs 9670-3p, 397-5p, and 5384-3p exhibited the most significant upregulation, whereas miRNAs 319, 9773, and 9774 were the most downregulated. The comparison of infection versus the control group for the cultivar Mace showed temperature-dependent regulation of these miRNAs. The principal component analysis confirmed that less abundant miRNAs among differentially expressed miRNAs were strongly correlated with the inoculated symptomatic wheat cultivars. Notably, miRNAs 397-5p, 398, and 9670-3p were upregulated in response to WSMV and TriMV infections, an observation not yet reported in this context. The significant upregulation of these three miRNAs was further confirmed with RT-qPCR analysis; in general, the RT-qPCR results were in agreement with our computational analysis. Target prediction analysis showed that the miRNAs standing out in our analysis targeted genes involved in defense response and regulation of transcription.

CONCLUSION

Investigation into the roles of these miRNAs and their corresponding targets holds promise for advancing our understanding of the mechanisms of virus infection and possible manipulation of these factors for developing durable virus resistance in crop plants.

MeRIP-Seq for Identifying Stress-Responsive Transcriptome-Wide m6A Profiles in Plants

Recent advancements in detection and mapping methods have enabled researchers to uncover the biological importance of RNA chemical modifications, which play a vital role in post-transcriptional gene regulation. Although numerous types of RNA modifications have been identified in higher eukaryotes, only a few have been extensively studied for their biological functions. Of these, N6-methyladenosine (m6A) is the most prevalent and important mRNA modification that influences various aspects of RNA metabolism, including mRNA stability, degradation, splicing, alternative polyadenylation, export, and localization, as well as translation. Thus, they have implications for a variety of biological processes, including growth, development, and stress responses. The m6A deposition or removal on transcripts is dynamic and is altered in response to internal and external cues. Because this mark can alter gene expression under stress conditions, it is essential to identify the transcripts that can acquire or lose this epitranscriptomic mark upon exposure to stress conditions. Here we describe a step-by-step protocol for identifying stress-responsive transcriptome-wide m6A changes using RNA immunoprecipitation followed by high-throughput sequencing (MeRIP-seq).

Molecular and genetic perspectives of cold tolerance in wheat

Environmental variation is the most crucial problem as it is causing food insecurity and negatively impacts food availability, utilization, assessment, and stability. Wheat is the largest and extensively cultivated staple food crop for fulfilling global food requirements. Abiotic stresses including salinity, heavy metal toxicity, drought, extreme temperatures, and oxidative stresses being the primary cause of productivity loss are a serious threat to agronomy. Cold stress is a foremost ecological constraint that is extremely influencing plant development, and yield. It is extremely hampering the propagative development of plant life. The structure and function of plant cells depend on the cell's immune system. The stresses due to cold, affect fluid in the plasma membrane and change it into crystals or a solid gel phase. Plants being sessile in nature have evolved progressive systems that permit them to acclimatize the cold stress at the physiological as well as molecular levels. The phenomenon of acclimatisation of plants to cold stress has been investigated for the last 10 years. Studying cold tolerance is critical for extending the adaptability zones of perennial grasses. In the present review, we have elaborated the current improvement of cold tolerance in plants from molecular and physiological viewpoints, such as hormones, the role of the posttranscriptional gene, micro RNAs, ICE-CBF-COR signaling route in cold acclimatization and how they are stimulating the expression of underlying genes encoding osmoregulatory elements and strategies to improve cold tolerance in wheat.

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.

Comprehensive, integrative genomic analysis of microRNA expression profiles in different tissues of two wheat cultivars with different traits

Wheat is one of the most important food sources on Earth. MicroRNAs (miRNA) play important roles in wheat productivity. To identify wheat miRNAs, we constructed and sequenced sRNA libraries from leaves and roots of two wheat cultivars (RAC875 and Kukri) with many different traits. Given that available miRNA wheat complement in the plant-specific database PmiREN ( https://pmiren.com ) does not include root tissues and root-associated miRNAs might thus be missing, we performed first the prediction of novel miRNAs using the sRNAbench tool. We found a total of 150 putatively novel miRNA genes with expression of both arms from 289 unique mature sequences and nearly 30% of all miRNA reads in roots corresponded to novel miRNAs. In contrast, this figure in leaves dropped to under 3%, confirming the undersampling of roots in the complement of known miRNAs. By using 120 publicly available wheat datasets, 598 Zea mays small RNA libraries, 64 plant species genomes, wheat degradome library, and functional enrichment analysis, a subset of novel miRNAs were confirmed as bona-fide miRNAs. Of the total 605 miRNAs identified in this study inclusive of 316 known miRNAs, 528 miRNAs were shared by both cultivars, 429 miRNAs were shared by both root tissues and 329 miRNAs were shared by both leaf tissues. In addition, 32 miRNAs were specific to Kukri while 45 miRNAs were specific to RAC875. These miRNAs had diverse functions, such as regulation of gene transcription, protein translation, energy metabolism, and cell cycle progression. Our data provide a genome-wide miRNA expression profile in these two wheat cultivars and help functional studies of wheat genomics.

Characterization of small interfering RNAs derived from pepper mild mottle virus in infected pepper plants by high-throughput sequencing

Pepper mild mottle virus (PMMoV) infects pepper plants and induces severe yield losses in China. However, the molecular interaction between PMMoV and pepper plants is largely unknown. RNA silencing is a eukaryotically conserved mechanism against viruses mediated by virus-derived small interfering RNAs (vsiRNAs) in plants. In this study, the profiles of vsiRNAs from PMMoV in infected pepper plants were obtained by high-throughput sequencing. The results showed that vsiRNAs were predominantly 21 and 22 nucleotides (nts) in length, and had a U bias at the 5'-terminal. The single-nucleotide resolution maps revealed that vsiRNAs were heterogeneously distributed throughout PMMoV genomic RNAs and hotspots of sense and antisense strands were mainly located in the RdRp and CP coding regions. The host transcripts targeted by vsiRNAs were predicted and they are mainly involved in physiological pathways related to stress response, cell regulation, and metabolism process. In addition, PMMoV infection induced significant up-regulation of CaAGO1a/1b/2, CaDCL2 and CaRDR1 gene transcripts in pepper plants, which are important components involved in antiviral RNA silencing pathway. Taken together, our results suggest the possible roles of vsiRNAs in PMMoV-pepper interactions.

Temporal expression study of miRNAs in the crown tissues of winter wheat grown under natural growth conditions

BACKGROUND

Winter wheat requires prolonged exposure to low temperature to initiate flowering (vernalization). Shoot apical meristem of the crown is the site of cold perception, which produces leaf primordia during vegetative growth before developing into floral primordia at the initiation of the reproductive phase. Although many essential genes for winter wheat cold acclimation and floral initiation have been revealed, the importance of microRNA (miRNA) meditated post-transcriptional regulation in crowns is not well understood. To understand the potential roles of miRNAs in crown tissues, we performed a temporal expression study of miRNAs in crown tissues at the three-leaf stage, winter dormancy stage, spring green-up stage, and jointing stage of winter wheat grown under natural growth conditions.

RESULTS

In total, 348 miRNAs belonging to 298 miRNA families, were identified in wheat crown tissues. Among them, 92 differentially expressed miRNAs (DEMs) were found to be significantly regulated from the three-leaf stage to the jointing stage. Most of these DEMs were highly expressed at the three-leaf stage and winter dormancy stage, and then declined in later stages. Six DEMs, including miR156a-5p were markedly induced during the winter dormancy stage. Eleven DEMs, including miR159a.1, miR390a-5p, miR393-5p, miR160a-5p, and miR1436, were highly expressed at the green-up stage. Twelve DEMs, such as miR172a-5p, miR394a, miR319b-3p, and miR9676-5p were highly induced at the jointing stage. Moreover, 14 novel target genes of nine wheat or Pooideae-specific miRNAs were verified using RLM-5' RACE assay. Notably, six mTERFs and two Rf1 genes, which are associated with mitochondrial gene expression, were confirmed as targets of three wheat-specific miRNAs.

CONCLUSIONS

The present study not only confirmed the known miRNAs associated with phase transition and floral development, but also identified a number of wheat or Pooideae-specific miRNAs critical for winter wheat cold acclimation and floral development. Most importantly, this study provided experimental evidence that miRNA could regulate mitochondrial gene expression by targeting mTERF and Rf1 genes. Our study provides valuable information for further exploration of the mechanism of miRNA mediated post-transcriptional regulation during winter wheat vernalization and inflorescent initiation.

PVsiRNAPred: Prediction of plant exclusive virus-derived small interfering RNAs by deep convolutional neural network

Plant exclusive virus-derived small interfering RNAs (vsiRNAs) regulate various biological processes, especially important in antiviral immunity. The identification of plant vsiRNAs is important for understanding the biogenesis and function mechanisms of vsiRNAs and further developing anti-viral plants. In this study, we extracted plant vsiRNA sequences from the PVsiRNAdb database. We then utilized deep convolutional neural network (CNN) to develop a deep learning algorithm for predicting plant vsiRNAs based on vsiRNA sequence composition, known as PVsiRNAPred. The key part of PVsiRNAPred is the CNN module, which automatically learns hierarchical representations of vsiRNA sequences related to vsiRNA profiles in plants. When evaluated using an independent testing dataset, the accuracy of the model was 65.70%, which was higher than those of five conventional machine learning method-based classifiers. In addition, PVsiRNAPred obtained a sensitivity of 67.11%, specificity of 64.26% and Matthews correlation coefficient (MCC) of 0.31, and the area under the receiver operating characteristic (ROC) curve (AUC) of PVsiRNAPred was 0.71 in the independent test. The permutation test with 1000 shuffles resulted in a p value of<0.001. The above results reveal that PVsiRNAPred has favorable generalization capabilities. We hope PVsiRNAPred, the first bioinformatics algorithm for predicting plant vsiRNAs, will allow efficient discovery of new vsiRNAs.

Analysis of miRNAs in Two Wheat Cultivars Infected With Puccinia striiformis f. sp. tritici

MicroRNAs are small RNAs that regulate gene expression in eukaryotes. In this study, we analyzed the small RNA profiles of two cultivars that exhibit different reactions to stripe rust infection: one susceptible, the other partially resistant. Using small RNA libraries prepared from the two wheat cultivars infected with stripe rust fungus (Puccinia striiformis f. sp. tritici), we identified 182 previously known miRNAs, 91 variants of known miRNAs, and 163 candidate novel wheat miRNAs. Known miRNA loci were usually copied in all three wheat sub-genomes, whereas novel miRNA loci were often specific to a single sub-genome. DESeq2 analysis of differentially expressed microRNAs revealed 23 miRNAs that exhibit cultivar-specific differences. TA078/miR399b showed cultivar-specific differential regulation in response to infection. Using different target prediction algorithms, 145 miRNAs were predicted to target wheat genes, while 69 miRNAs were predicted to target fungal genes. We also confirmed reciprocal expression of TA078/miR399b and tae-miR9664 and their target genes in different treatments, providing evidence for miRNA-mediated regulation during infection. Both known and novel miRNAs were predicted to target fungal genes, suggesting trans-kingdom regulation of gene expression. Overall, this study contributes to the current repository of wheat miRNAs and provides novel information on the yet-uncharacterized roles for miRNAs in the wheat-stripe rust pathosystem.

Identification and Temporal Expression Analysis of Conserved and Novel MicroRNAs in the Leaves of Winter Wheat Grown in the Field

Cold acclimation and vegetative/reproductive transition are two important evolutionary adaptive mechanisms for winter wheat surviving the freezing temperature in winter and successful seeds setting in the next year. MicroRNA (miRNA) is a class of regulatory small RNAs (sRNAs), which plays critical roles in the growth and development of plants. However, the regulation mechanism of miRNAs during cold acclimation and vegetative/reproductive transition of winter wheat is not much understood. In this study, four sRNA libraries from leaves of winter wheat grown in the field at the three-leaf stage, winter dormancy stage, spring green-up stage, and jointing stage were analyzed to identify known and novel miRNAs and to understand their potential roles in the growth and development of winter wheat. We examined miRNA expression using a high-throughput sequencing technique. A total of 373 known, 55 novel, and 27 putative novel miRNAs were identified. Ninety-one miRNAs were found to be differentially expressed at the four stages. Among them, the expression of six known and eight novel miRNAs was significantly suppressed at the winter dormancy stage, whereas the expression levels of seven known and eight novel miRNAs were induced at this stage; three known miRNAs and three novel miRNAs were significantly induced at the spring green-up stage; six known miRNAs were induced at the spring green-up stage and reached the highest expression level at the jointing stage; and 20 known miRNAs and 10 novel miRNAs were significantly induced at the jointing stage. Expression of a number of representative differentially expressed miRNAs was verified using quantitative real-time polymerase chain reaction (qRT-PCR). Potential target genes for known and novel miRNAs were predicted. Moreover, six novel target genes for four Pooideae species-specific miRNAs and two novel miRNAs were verified using the RNA ligase-mediated 5'-rapid amplification of cDNA ends (RLM-5'RACE) technique. These results indicate that miRNAs are key non-coding regulatory factors modulating the growth and development of wheat. Our study provides valuable information for in-depth understanding of the regulatory mechanism of miRNAs in cold acclimation and vegetative/reproductive transition of winter wheat grown in the field.

Genotype- and tissue-specific miRNA profiles and their targets in three alfalfa (Medicago sativa L) genotypes

BACKGROUND

Alfalfa (Medicago sativa L.) is a forage legume with significant agricultural value worldwide. MicroRNAs (miRNAs) are key components of post-transcriptional gene regulation and essentially regulate many aspects of plant growth and development. Although miRNAs were reported in alfalfa, their expression profiles in different tissues and the discovery of novel miRNAs as well as their targets have not been described in this plant species.

RESULTS

To identify tissue-specific miRNA profiles in whole plants, shoots and roots of three different alfalfa genotypes (Altet-4, NECS-141and NF08ALF06) were used. Small RNA libraries were generated and sequenced using a high-throughput sequencing platform. Analysis of these libraries enabled identification of100 miRNA families; 21 of them belong to the highly conserved families while the remaining 79 families are conserved at the minimum between M. sativa and the model legume and close relative, M. truncatula. The profiles of the six abundantly expressed miRNA families (miR156, miR159, miR166, miR319, miR396 and miR398) were relatively similar between the whole plants, roots and shoots of these three alfalfa genotypes. In contrast, robust differences between shoots and roots for miR160 and miR408 levels were evident, and their expression was more abundant in the shoots. Additionally, 17 novel miRNAs were identified and the relative abundance of some of these differed between tissue types. Further, the generation and analysis of degradome libraries from the three alfalfa genotypes enabled confirmation of 69 genes as targets for 31 miRNA families in alfalfa.

CONCLUSIONS

The miRNA profiles revealed both similarities and differences in the expression profiles between tissues within a genotype as well as between the genotypes. Among the highly conserved miRNA families, miR166 was the most abundantly expressed in almost all tissues from the three genotypes. The identification of conserved and novel miRNAs as well as their targets in different tissues of multiple genotypes increased our understanding of miRNA-mediated gene regulation in alfalfa and could provide valuable insights for practical research and plant improvement applications in alfalfa and related legume species.

Genome-wide identification of conserved and novel microRNAs in one bud and two tender leaves of tea plant (Camellia sinensis) by small RNA sequencing, microarray-based hybridization and genome survey scaffold sequences

BACKGROUND

MicroRNAs (miRNAs) are important for plant growth and responses to environmental stresses via post-transcriptional regulation of gene expression. Tea, which is primarily produced from one bud and two tender leaves of the tea plant (Camellia sinensis), is one of the most popular non-alcoholic beverages worldwide owing to its abundance of secondary metabolites. A large number of miRNAs have been identified in various plants, including non-model species. However, due to the lack of reference genome sequences and/or information of tea plant genome survey scaffold sequences, discovery of miRNAs has been limited in C. sinensis.

RESULTS

Using small RNA sequencing, combined with our recently obtained genome survey data, we have identified and analyzed 175 conserved and 83 novel miRNAs mainly in one bud and two tender leaves of the tea plant. Among these, 93 conserved and 18 novel miRNAs were validated using miRNA microarray hybridization. In addition, the expression pattern of 11 conserved and 8 novel miRNAs were validated by stem-loop-qRT-PCR. A total of 716 potential target genes of identified miRNAs were predicted. Further, Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that most of the target genes were primarily involved in stress response and enzymes related to phenylpropanoid biosynthesis. The predicted targets of 4 conserved miRNAs were further validated by 5'RLM-RACE. A negative correlation between expression profiles of 3 out of 4 conserved miRNAs (csn-miR160a-5p, csn-miR164a, csn-miR828 and csn-miR858a) and their targets (ARF17, NAC100, WER and MYB12 transcription factor) were observed.

CONCLUSION

In summary, the present study is one of few such studies on miRNA detection and identification in the tea plant. The predicted target genes of majority of miRNAs encoded enzymes, transcription factors, and functional proteins. The miRNA-target transcription factor gene interactions may provide important clues about the regulatory mechanism of these miRNAs in the tea plant. The data reported in this study will make a huge contribution to knowledge on the potential miRNA regulators of the secondary metabolism pathway and other important biological processes in C. sinensis.

Isolation of Early-Responsive MicroRNA From Diuraphis noxia (Hemiptera: Aphididae)-Resistant Wheat

The Russian wheat aphid (Diuraphis noxia Kurdjumov) is an economically important pest of small grains in many countries. The past decades have seen the deployment of resistance-carrying wheat (Triticum aestivum L.) cultivars to control D. noxia. However, the emergence of resistance-breaking biotypes is negating this strategy. The role that noncoding RNA (ncRNA) molecules play in the wheat-D. noxia interaction has not been studied to date. This study aimed to isolate differentially regulated microRNA from a resistant and susceptible near-isogenic wheat line after aphid infestation. Twenty-seven identified miRNA were mostly related to stress-linked miRNA, and their predicted targets were linked with known D. noxia-feeding regulated proteins. These included transcription factors, signaling proteins, carbohydrate metabolism, and disease resistance pathways. Gene expression of three putative miRNAs and a predicted nucleotide-binding leucine-rich repeat gene with an identified miRNA target site in the NB-ARC domain displayed differential regulation between the resistant and susceptible plants. This study marks the initial investigation into understanding the role of ncRNA in a D. noxia-resistant wheat line after infestation and reports a correlation between a miRNA and its putative target for this interaction.

isomiR2Function: An Integrated Workflow for Identifying MicroRNA Variants in Plants

In plants, post transcriptional regulation by non-coding RNAs (ncRNAs), in particular miRNAs (19-24 nt) has been involved in modulating the transcriptional landscape in developmental, biotic and abiotic interactions. In past few years, considerable focus has been leveraged on delineating and deciphering the role of miRNAs and their canonical isomiRs in plants. However, proper classification and accurate prediction of plant isomiRs taking into account the relative features by which we define isomiRs, such as templated or non-templated is still lacking. In the present research, we present isomiR2Function, a standalone easily deployable tool that allows for the robust and high-throughput discovery of templated and non-templated isomiRs. Additionally, isomiR2Function allows for identification of differentially expressed isomiRs and in parallel target prediction based on both transcripts or PARE-Seq either using Targetfinder or Cleaveland. isomiR2Function allows for the functional enrichment of the detected targets using TopGO package. Benchmarking of isomiR2Function revealed highly accurate prediction and classification of isomiRs as compared to the previously developed isomiR prediction tools. Additionally, the downstream implementation of additional features allows isomiR2Function to be classified as a single standalone tool for isomiR profiling from discovery to functional roles. All in all, isomiR2Function allows the streamline processing of the miRNA-seq for the identification and characterization of isomiRs with minimal efforts. isomiR2Function can be accessed through: https://github.com/347033139/isomiR2Function.

Small RNA profiles in soybean primary root tips under water deficit

BACKGROUND

Soybean (Glycine max) production is significantly hampered by frequent droughts in many regions of the world including the United States. Identifying microRNA (miRNA)-controlled posttranscriptional gene regulation under drought will enhance our understanding of molecular basis of drought tolerance in this important cash crop. Indeed, miRNA profiles in soybean exposed to drought were studied but not from the primary root tips, which is not only a main zone of water uptake but also critical for water stress sensing and signaling.

METHODS

Here we report miRNA profiles specifically from well-watered and water-stressed primary root tips (0 to 8 mm from the root apex) of soybean. Small RNA sequencing confirmed the expression of vastly diverse miRNA (303 individual miRNAs) population, and, importantly several conserved miRNAs were abundantly expressed in primary root tips.

RESULTS

Notably, 12 highly conserved miRNA families were differentially regulated in response to water-deficit; six were upregulated while six others were downregulated at least by one fold (log2) change. Differentially regulated soybean miRNAs are targeting genes include auxin response factors, Cu/Zn Superoxide dismutases, laccases and plantacyanin and several others.

CONCLUSIONS

These results highlighted the importance of miRNAs in primary root tips both under control and water-deficit conditions; under control conditions, miRNAs could be important for cell division, cell elongation and maintenance of the root apical meristem activity including quiescent centre whereas under water stress differentially regulated miRNAs could decrease auxin signaling and oxidative stress as well as other metabolic processes that save energy and water.

RNA-seq studies using wheat PHYTOCHROME B and PHYTOCHROME C mutants reveal shared and specific functions in the regulation of flowering and shade-avoidance pathways

BACKGROUND

In cereal crops such as wheat, an optimal timing of developmental transitions is required to maximize grain yield. Many of these developmental changes are precisely regulated by changes in the duration, intensity or quality of light. Phytochromes are dimeric photoreceptors that absorb light maximally in the red and far-red wavelengths and induce large-scale transcriptional changes in response to variation in light quality. In wheat, PHYC is required for early flowering under long days. However, it is currently unknown whether this function requires the presence of PHYB. In this study, we characterized the role of PHYB in wheat development and used RNA-seq to analyze and compare the transcriptomes of phyB-null and phyC-null TILLING mutants.

RESULTS

Under long-day photoperiods, phyB-null plants exhibit a severe delay in flowering comparable to the delay observed in phyC-null plants. These results demonstrate that both genes are required for the induction of wheat flowering under long days. Using replicated RNA-seq studies we identified 82 genes that are significantly up or down regulated in both the phyB-null and phyC-null mutant relative to their respective wild-type controls. Among these genes are several well-characterized positive regulators of flowering, including PPD1, FT1 and VRN1. Eight-fold more genes were differentially regulated only in the phyB-null mutant (2202) than only in the phyC-null mutant (261). The PHYB-regulated genes were enriched in components of the auxin, gibberellin and brassinosteroid biosynthesis and signaling pathways, and in transcription factors with putative roles in regulating vegetative development and shade-avoidance responses. Several genes involved in abiotic stress tolerance pathways were also found to be regulated by PHYB.

CONCLUSIONS

PHYB and PHYC are both required for the photoperiodic induction of wheat flowering, whereas PHYB alone regulates a large number of genes involved in hormone biosynthesis and signaling, shade-avoidance response, and abiotic stress tolerance. Our analysis provides a comprehensive overview of the PHYB- and PHYC-mediated transcriptional changes during light signaling, and an initial step towards the dissection of this regulatory gene network in wheat. This further dissection will be required to explore the individual phytochrome-mediated developmental responses and to evaluate their potential to improve wheat adaptation to changing environments.

Characterization of microRNAs and their targets in wild barley (Hordeum vulgare subsp. spontaneum) using deep sequencing

MicroRNAs (miRNA) are a class of small, endogenous RNAs that play a negative regulatory role in various developmental and metabolic processes of plants. Wild barley (Hordeum vulgare subsp. spontaneum), as the progenitor of cultivated barley (Hordeum vulgare subsp. vulgare), has served as a valuable germplasm resource for barley genetic improvement. To survey miRNAs in wild barley, we sequenced the small RNA library prepared from wild barley using the Illumina deep sequencing technology. A total of 70 known miRNAs and 18 putative novel miRNAs were identified. Sequence analysis revealed that all of the miRNAs identified in wild barley contained the highly conserved hairpin sequences found in barley cultivars. MiRNA target predictions showed that 12 out of 52 miRNA families were predicted to target transcription factors, including 8 highly conserved miRNA families in plants and 4 wheat-barley conserved miRNA families. In addition to transcription factors, other predicted target genes were involved in diverse physiological and metabolic processes and stress defense. Our study for the first time reported the large-scale investigation of small RNAs in wild barley, which will provide essential information for understanding the regulatory role of miRNAs in wild barley and also shed light on future practical utilization of miRNAs for barley improvement.

Characterization of Small RNAs Derived from tRNAs, rRNAs and snoRNAs and Their Response to Heat Stress in Wheat Seedlings

Small RNAs (sRNAs) derived from non-coding RNAs (ncRNAs), such as tRNAs, rRNAs and snoRNAs, have been identified in various organisms. Several observations have indicated that cleavage of tRNAs and rRNAs is induced by various stresses. To clarify whether sRNAs in wheat derived from tRNAs (stRNAs), rRNAs (srRNAs) and snoRNAs (sdRNAs) are produced specifically in association with heat stress responses, we carried out a bioinformatic analysis of sRNA libraries from wheat seedlings and performed comparisons between control and high-temperature-treated samples to measure the differential abundance of stRNAs, srRNAs and sdRNAs. We found that the production of sRNAs from tRNAs, 5.8S rRNAs, and 28S rRNAs was more specific than that from 5S rRNAs and 18S rRNAs, and more than 95% of the stRNAs were processed asymmetrically from the 3’ or 5’ ends of mature tRNAs. We identified 333 stRNAs and 8,822 srRNAs that were responsive to heat stress. Moreover, the expression of stRNAs derived from tRNA-Val-CAC, tRNA-Thr-UGU, tRNA-Tyr-GUA and tRNA-Ser-UGA was not only up-regulated under heat stress but also induced by osmotic stress, suggesting that the increased cleavage of tRNAs might be a mechanism that developed in wheat seedlings to help them cope with adverse environmental conditions.

Dissecting miRNAs in Wheat D Genome Progenitor, Aegilops tauschii

As the post-transcriptional regulators of gene expression, microRNAs or miRNAs comprise an integral part of understanding how genomes function. Although miRNAs have been a major focus of recent efforts, miRNA research is still in its infancy in most plant species. Aegilops tauschii, the D genome progenitor of bread wheat, is a wild diploid grass exhibiting remarkable population diversity. Due to the direct ancestry and the diverse gene pool, A. tauschii is a promising source for bread wheat improvement. In this study, a total of 87 Aegilops miRNA families, including 51 previously unknown, were computationally identified both at the subgenomic level, using flow-sorted A. tauschii 5D chromosome, and at the whole genome level. Predictions at the genomic and subgenomic levels suggested A. tauschii 5D chromosome as rich in pre-miRNAs that are highly associated with Class II DNA transposons. In order to gain insights into miRNA evolution, putative 5D chromosome miRNAs were compared to its modern ortholog, Triticum aestivum 5D chromosome, revealing that 48 of the 58 A. tauschii 5D miRNAs were conserved in orthologous T. aestivum 5D chromosome. The expression profiles of selected miRNAs (miR167, miR5205, miR5175, miR5523) provided the first experimental evidence for miR5175, miR5205 and miR5523, and revealed differential expressional changes in response to drought in different genetic backgrounds for miR167 and miR5175. Interestingly, while miR5523 coding regions were present and expressed as pre-miR5523 in both T. aestivum and A. tauschii, the expression of mature miR5523 was observed only in A. tauschii under normal conditions, pointing out to an interference at the downstream processing of pre-miR5523 in T. aestivum. Overall, this study expands our knowledge on the miRNA catalog of A. tauschii, locating a subset specifically to the 5D chromosome, with ample functional and comparative insight which should contribute to and complement efforts to develop drought tolerant wheat varieties.

Identification and characterization of microRNAs in Eucheuma denticulatum by high-throughput sequencing and bioinformatics analysis

Eucheuma denticulatum, an economically and industrially important red alga, is a valuable marine resource. Although microRNAs (miRNAs) play an essential role in gene post-transcriptional regulation, no research has been conducted to identify and characterize miRNAs in E. denticulatum. In this study, we identified 134 miRNAs (133 conserved miRNAs and one novel miRNA) from 2,997,135 small-RNA reads by high-throughput sequencing combined with bioinformatics analysis. BLAST searching against miRBase uncovered 126 potential miRNA families. A conservation and diversity analysis of predicted miRNA families in different plant species was performed by comparative alignment and homology searching. A total of 4 and 13 randomly selected miRNAs were respectively validated by northern blotting and stem-loop reverse transcription PCR, thereby demonstrating the reliability of the miRNA sequencing data. Altogether, 871 potential target genes were predicted using psRobot and TargetFinder. Target genes classification and enrichment were conducted based on Gene Ontology analysis. The functions of target gene products and associated metabolic pathways were predicted by Kyoto Encyclopedia of Genes and Genomes pathway analysis. A Cytoscape network was constructed to explore the interrelationships of miRNAs, miRNA-target genes and target genes. A large number of miRNAs with diverse target genes will play important roles for further understanding some essential biological processes in E. denticulatum. The uncovered information can serve as an important reference for the protection and utilization of this unique red alga in the future.

Biogenesis, Function, and Applications of Virus-Derived Small RNAs in Plants

RNA silencing, an evolutionarily conserved and sequence-specific gene-inactivation system, has a pivotal role in antiviral defense in most eukaryotic organisms. In plants, a class of exogenous small RNAs (sRNAs) originating from the infecting virus called virus-derived small interfering RNAs (vsiRNAs) are predominantly responsible for RNA silencing-mediated antiviral immunity. Nowadays, the process of vsiRNA formation and the role of vsiRNAs in plant viral defense have been revealed through deep sequencing of sRNAs and diverse genetic analysis. The biogenesis of vsiRNAs is analogous to that of endogenous sRNAs, which require diverse essential components including dicer-like (DCL), argonaute (AGO), and RNA-dependent RNA polymerase (RDR) proteins. vsiRNAs trigger antiviral defense through post-transcriptional gene silencing (PTGS) or transcriptional gene silencing (TGS) of viral RNA, and they hijack the host RNA silencing system to target complementary host transcripts. Additionally, several applications that take advantage of the current knowledge of vsiRNAs research are being used, such as breeding antiviral plants through genetic engineering technology, reconstructing of viral genomes, and surveying viral ecology and populations. Here, we will provide an overview of vsiRNA pathways, with a primary focus on the advances in vsiRNA biogenesis and function, and discuss their potential applications as well as the future challenges in vsiRNAs research.

Small RNAs from the wheat stripe rust fungus (Puccinia striiformis f.sp. tritici)

BACKGROUND

Wheat stripe rust, caused by Puccinia striiformis f. sp. tritici, is a costly global disease that burdens farmers with yield loss and high fungicide expenses. This sophisticated biotrophic parasite infiltrates wheat leaves and develops infection structures inside host cells, appropriating nutrients while suppressing the plant defense response. Development in most eukaryotes is regulated by small RNA molecules, and the success of host-induced gene silencing technology in Puccinia spp. implies the existence of a functional RNAi system. However, some fungi lack this capability, and small RNAs have not yet been reported in rust fungi. The objective of this study was to determine whether P. striiformis carries an endogenous small RNA repertoire.

RESULTS

We extracted small RNA from rust-infected wheat flag leaves and performed high-throughput sequencing. Two wheat cultivars were analyzed: one is susceptible; the other displays partial high-temperature adult plant resistance. Fungal-specific reads were identified by mapping to the P. striiformis draft genome and removing reads present in uninfected control libraries. Sequencing and bioinformatics results were verified by RT-PCR. Like other RNAi-equipped fungi, P. striiformis produces large numbers of 20-22 nt sequences with a preference for uracil at the 5' position. Precise post-transcriptional processing and high accumulation of specific sRNA sequences were observed. Some predicted sRNA precursors possess a microRNA-like stem-loop secondary structure; others originate from much longer inverted repeats containing gene sequences. Finally, sRNA-target prediction algorithms were used to obtain a list of putative gene targets in both organisms. Predicted fungal target genes were enriched for kinases and small secreted proteins, while the list of wheat targets included homologs of known plant resistance genes.

CONCLUSIONS

This work provides an inventory of small RNAs endogenous to an important plant pathogen, enabling further exploration of gene regulation on both sides of the host/parasite interaction. We conclude that small RNAs are likely to play a role in regulating the complex developmental processes involved in stripe rust pathogenicity.

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.

Computational analysis of miRNA-target community network reveals cross talk among different metabolisms

To date, only a few conserved miRNAs have been predicted in hexaploid (AABBDD) bread wheat and till now community behavior among miRNA is still in dark. Analysis of publically available 1287279 ESTs from NCBI resulted 262 putative pre-miRNAs and 39 novel mature miRNAs. A total 22,468 targets were identified on 21 chromosomes. MiRNA target community was identified for genomes with different levels of cross talks. Gene ontology of these community targets suggests their differential involvement in different metabolisms along with common and stringent involvement in nitrogen metabolism.

An integrative approach to identify hexaploid wheat miRNAome associated with development and tolerance to abiotic stress

Background

Wheat is a major staple crop with broad adaptability to a wide range of environmental conditions. This adaptability involves several stress and developmentally responsive genes, in which microRNAs (miRNAs) have emerged as important regulatory factors. However, the currently used approaches to identify miRNAs in this polyploid complex system focus on conserved and highly expressed miRNAs avoiding regularly those that are often lineage-specific, condition-specific, or appeared recently in evolution. In addition, many environmental and biological factors affecting miRNA expression were not yet considered, resulting still in an incomplete repertoire of wheat miRNAs.

Results

We developed a conservation-independent technique based on an integrative approach that combines machine learning, bioinformatic tools, biological insights of known miRNA expression profiles and universal criteria of plant miRNAs to identify miRNAs with more confidence. The developed pipeline can potentially identify novel wheat miRNAs that share features common to several species or that are species specific or clade specific. It allowed the discovery of 199 miRNA candidates associated with different abiotic stresses and development stages. We also highlight from the raw data 267 miRNAs conserved with 43 miRBase families. The predicted miRNAs are highly associated with abiotic stress responses, tolerance and development. GO enrichment analysis showed that they may play biological and physiological roles associated with cold, salt and aluminum (Al) through auxin signaling pathways, regulation of gene expression, ubiquitination, transport, carbohydrates, gibberellins, lipid, glutathione and secondary metabolism, photosynthesis, as well as floral transition and flowering.

Conclusion

This approach provides a broad repertoire of hexaploid wheat miRNAs associated with abiotic stress responses, tolerance and development. These valuable resources of expressed wheat miRNAs will help in elucidating the regulatory mechanisms involved in freezing and Al responses and tolerance mechanisms as well as for development and flowering. In the long term, it may help in breeding stress tolerant plants.

Electronic supplementary material

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

Transcriptome profiling of root microRNAs reveals novel insights into taproot thickening in radish (Raphanus sativus L.)

BACKGROUND

Radish (Raphanus sativus L.) is an economically important root vegetable crop, and the taproot-thickening process is the most critical period for the final productivity and quality formation. MicroRNAs (miRNAs) are a family of non-coding small RNAs that play an important regulatory function in plant growth and development. However, the characterization of miRNAs and their roles in regulating radish taproot growth and thickening remain largely unexplored. A Solexa high-throughput sequencing technology was used to identify key miRNAs involved in taproot thickening in radish.

RESULTS

Three small RNA libraries from 'NAU-YH' taproot collected at pre-cortex splitting stage, cortex splitting stage and expanding stage were constructed. In all, 175 known and 107 potential novel miRNAs were discovered, from which 85 known and 13 novel miRNAs were found to be significantly differentially expressed during taproot thickening. Furthermore, totally 191 target genes were identified for the differentially expressed miRNAs. These target genes were annotated as transcription factors and other functional proteins, which were involved in various biological functions including plant growth and development, metabolism, cell organization and biogenesis, signal sensing and transduction, and plant defense response. RT-qPCR analysis validated miRNA expression patterns for five miRNAs and their corresponding target genes.

CONCLUSIONS

The small RNA populations of radish taproot at different thickening stages were firstly identified by Solexa sequencing. Totally 98 differentially expressed miRNAs identified from three taproot libraries might play important regulatory roles in taproot thickening. Their targets encoding transcription factors and other functional proteins including NF-YA2, ILR1, bHLH74, XTH16, CEL41 and EXPA9 were involved in radish taproot thickening. These results could provide new insights into the regulatory roles of miRNAs during the taproot thickening and facilitate genetic improvement of taproot in radish.

Identification and comparative analysis of differentially expressed miRNAs in leaves of two wheat (Triticum aestivum L.) genotypes during dehydration stress

BACKGROUND

MicroRNAs (miRNAs) play critical roles in the processes of plant growth and development, but little is known of their functions during dehydration stress in wheat. Moreover, the mechanisms by which miRNAs confer different levels of dehydration stress tolerance in different wheat genotypes are unclear.

RESULTS

We examined miRNA expressions in two different wheat genotypes, Hanxuan10, which is drought-tolerant, and Zhengyin1, which is drought-susceptible. Using a deep-sequencing method, we identified 367 differentially expressed miRNAs (including 46 conserved miRNAs and 321 novel miRNAs) and compared their expression levels in the two genotypes. Among them, 233 miRNAs were upregulated and 10 were downregulated in both wheat genotypes after dehydration stress. Interestingly, 13 miRNAs exhibited opposite patterns of expression in the two wheat genotypes, downregulation in the drought-tolerant cultivar and upregulation in the drought-susceptible cultivar. We also identified 111 miRNAs that were expressed predominantly in only one or the other genotype after dehydration stress. We verified the expression patterns of a number of representative miRNAs using qPCR analysis and northern blot, which produced results consistent with those of the deep-sequencing method. Moreover, monitoring the expression levels of 10 target genes by qPCR analysis revealed negative correlations with the levels of their corresponding miRNAs.

CONCLUSIONS

These results indicate that differentially expressed patterns of miRNAs between these two genotypes may play important roles in dehydration stress tolerance in wheat and may be a key factor in determining the levels of stress tolerance in different wheat genotypes.

From Genetics to Functional Genomics: Improvement in Drought Signaling and Tolerance in Wheat

Drought being a yield limiting factor has become a major threat to international food security. It is a complex trait and drought tolerance response is carried out by various genes, transcription factors (TFs), microRNAs (miRNAs), hormones, proteins, co-factors, ions, and metabolites. This complexity has limited the development of wheat cultivars for drought tolerance by classical breeding. However, attempts have been made to fill the lost genetic diversity by crossing wheat with wild wheat relatives. In recent years, several molecular markers including single nucleotide polymorphisms (SNPs) and quantitative trait loci (QTLs) associated with genes for drought signaling pathways have been reported. Screening of large wheat collections by marker assisted selection (MAS) and transformation of wheat with different genes/TFs has improved drought signaling pathways and tolerance. Several miRNAs also provide drought tolerance to wheat by regulating various TFs/genes. Emergence of OMICS techniques including transcriptomics, proteomics, metabolomics, and ionomics has helped to identify and characterize the genes, proteins, metabolites, and ions involved in drought signaling pathways. Together, all these efforts helped in understanding the complex drought tolerance mechanism. Here, we have reviewed the advances in wide hybridization, MAS, QTL mapping, miRNAs, transgenic technique, genome editing system, and above mentioned functional genomics tools for identification and utility of signaling molecules for improvement in wheat drought tolerance.

Discovery and comparative profiling of microRNAs in representative monopodial bamboo (Phyllostachys edulis) and sympodial bamboo (Dendrocalamus latiflorus)

Background

According to the growth pattern of bamboo, sympodial bamboo and monopodial bamboo are considered as two mainly kinds of bamboo. They have different phenotypes and different characteristics in developmental stage. Much attention had been paid on the study of bamboo cultivation, processing, physiology, biochemistry and molecular biology, which had made great progresses in the last decade, especially for the highlighted achievement of the bamboo genomics. However, there is no information available on concerning comparative profiling of miRNAs between sympodial bamboo and monopodial bamboo, which might play important roles in the regulation of bamboo development.

Methodology/Principal Findings

We identified the profiles of small RNAs using leaf tissues from one sympodial bamboo i.e. moso bamboo (Phyllostachys edulis) and another monopodial bamboo i.e. ma bamboo (Dendrocalamus latiflorus). The result showed that there were 19,295,759 and 11,513,888 raw sequence reads, in which 92 and 69 conserved miRNAs, as well as 95 and 62 novel miRNAs were identified in moso bamboo and ma bamboo, respectively. The ratio of high conserved miRNA families in ma bamboo is more than that in moso bamboo. In addition, a total of 49 and 106 potential targets were predicted in moso bamboo and ma bamboo, respectively, in which several targets for novel miRNAs are transcription factors that play important roles in plant development. More importantly, annotation of differentially expressed target genes was performed based on the analysis of pathway and gene ontology terms enrichment.

Conclusions/Significance

This study provides the first large-scale sight of discovery and comparative characterization of miRNAomes between two representative bamboos belonged to sympodial bamboo and monopodial bamboo, respectively. Although it will be necessary to validate the function of miRNAs through more experimental research in further, these results lay a foundation for unraveling the miRNA-mediated molecular processes in different kinds of bamboo.

A comprehensive genome-wide study on tissue-specific and abiotic stress-specific miRNAs in Triticum aestivum

Productivity of wheat crop is largely dependent on its growth and development that, in turn, is mainly regulated by environmental conditions, including abiotic stress factors. miRNAs are key regulators of gene expression networks involved in diverse aspects of development and stress responses in plants. Using high-throughput sequencing of eight small RNA libraries prepared from diverse abiotic stresses and tissues, we identified 47 known miRNAs belonging to 20 families, 49 true novel and 1030 candidate novel miRNAs. Digital gene expression analysis revealed that 257 miRNAs exhibited tissue-specific expression and 74 were associated with abiotic stresses. Putative target genes were predicted for miRNAs identified in this study and their grouping into functional categories indicated that the putative targets were involved in diverse biological processes. RLM-RACE of predicted targets of three known miRNAs (miR156, miR160 and miR164) confirmed their mRNA cleavage, thus indicating their regulation at post-transcriptional level by the corresponding miRNAs. Mapping of the sequenced data onto the wheat progenitors and closely related monocots revealed a large number of evolutionary conserved miRNAs. Additional expression profiling of some of these miRNAs in other abiotic stresses underline their involvement in multiple stresses. Our findings provide valuable resource for an improved understanding of the role of miRNAs in stress tolerance as well as plant development.

Extending the sRNAome of apple by next-generation sequencing

The global importance of apple as a fruit crop necessitates investigations into molecular aspects of the processes that influence fruit quality and yield, including plant development, fruit ripening and disease resistance. In order to study and understand biological processes it is essential to recognise the range of molecules, which influence these processes. Small non-coding RNAs are regulatory agents involved in diverse plant activities, ranging from development to stress response. The occurrence of these molecules in apple leaves was studied by means of next-generation sequencing. 85 novel microRNA (miRNA) gene loci were predicted and characterized along with known miRNA loci. Both cis- and trans-natural antisense transcript pairs were identified. Although the trans-overlapping regions were enriched in small RNA (sRNA) production, cis-overlaps did not seem to agree. More than 150 phased regions were also identified, and for a small subset of these, potential miRNAs that could initiate phasing, were revealed. Repeat-associated siRNAs, which are generated from repetitive genomic regions such as transposons, were also analysed. For this group almost all available repeat sequences, associated with the apple genome and present in Repbase, were found to produce siRNAs. Results from this study extend our current knowledge on apple sRNAs and their precursors significantly. A rich molecular resource has been created and is available to the research community to serve as a baseline for future studies.

Identification and characterization of microRNAs in the flag leaf and developing seed of wheat (Triticum aestivum L.)

BACKGROUND

MicroRNAs (miRNAs) regulate various biological processes in plants. Considerable data are available on miRNAs involved in the development of rice, maize and barley. In contrast, little is known about miRNAs and their functions in the development of wheat. In this study, five small RNA (sRNA) libraries from wheat seedlings, flag leaves, and developing seeds were developed and sequenced to identify miRNAs and understand their functions in wheat development.

RESULTS

Twenty-four known miRNAs belonging to 15 miRNA families were identified from 18 MIRNA loci in wheat in the present study, including 15 miRNAs (9 MIRNA loci) first identified in wheat, 13 miRNA families (16 MIRNA loci) being highly conserved and 2 (2 MIRNA loci) moderately conserved. In addition, fifty-five novel miRNAs were also identified. The potential target genes for 15 known miRNAs and 37 novel miRNAs were predicted using strict criteria, and these target genes are involved in a wide range of biological functions. Four of the 15 known miRNA families and 22 of the 55 novel miRNAs were preferentially expressed in the developing seeds with logarithm (log2) of the fold change of 1.0 ~ 7.6, and half of them were seed-specific, suggesting that they participate in regulating wheat seed development and metabolism. From 5 days post-anthesis to 20 days post-anthesis, miR164 and miR160 increased in abundance in the developing seeds, whereas miR169 decreased, suggesting their coordinating functions in the different developmental stages of wheat seed. Moreover, 8 known miRNA families and 28 novel miRNAs exhibited tissue-biased expression in wheat flag leaves, with the logarithm of the fold changes of 0.1 ~ 5.2. The putative targets of these tissue-preferential miRNAs were involved in various metabolism and biological processes, suggesting complexity of the regulatory networks in different tissues. Our data also suggested that wheat flag leaves have more complicated regulatory networks of miRNAs than developing seeds.

CONCLUSIONS

Our work identified and characterised wheat miRNAs, their targets and expression patterns. This study is the first to elucidate the regulatory networks of miRNAs involved in wheat flag leaves and developing seeds, and provided a foundation for future studies on specific functions of these miRNAs.

New wheat microRNA using whole-genome sequence

MicroRNAs are post-transcriptional regulators of gene expression, taking roles in a variety of fundamental biological processes. Hence, their identification, annotation and characterization are of great significance, especially in bread wheat, one of the main food sources for humans. The recent availability of 5× coverage Triticum aestivum L. whole-genome sequence provided us with the opportunity to perform a systematic prediction of a complete catalogue of wheat microRNAs. Using an in silico homology-based approach, stem-loop coding regions were derived from two assemblies, constructed from wheat 454 reads. To avoid the presence of pseudo-microRNAs in the final data set, transposable element related stem-loops were eliminated by repeat analysis. Overall, 52 putative wheat microRNAs were predicted, including seven, which have not been previously published. Moreover, with distinct analysis of the two different assemblies, both variety and representation of putative microRNA-coding stem-loops were found to be predominant in the intergenic regions. By searching available expressed sequences and small RNA library databases, expression evidence for 39 (out of 52) putative wheat microRNAs was provided. Expression of three of the predicted microRNAs (miR166, miR396 and miR528) was also comparatively quantified with real-time quantitative reverse transcription PCR. This is the first report on in silico prediction of a whole repertoire of bread wheat microRNAs, supported by the wet-lab validation.

Identification and characterization of microRNAs in the leaf of ma bamboo (Dendrocalamus latiflorus) by deep sequencing

MicroRNAs (miRNAs), a class of non-coding small endogenous RNAs of approximately 22 nucleotides, regulate gene expression at the post-transcriptional levels by targeting mRNAs for degradation or by inhibiting protein translation. Thousands of miRNAs have been identified in many species. However, there is no information available concerning miRNAs in ma bamboo (Dendrocalamus latiflorus), one of the most important non-timber forest products, which has essential ecological roles in forests. To identify miRNAs in D. latiflorus, a small RNA library was constructed from leaf tissues. Using next generation high-throughput sequencing technology and bioinformatics analysis, we obtained 11,513,607 raw sequence reads and identified 84 conserved miRNAs (54 mature miRNAs and 30 star miRNAs) belonging to 17 families, and 81 novel miRNAs (76 mature miRNAs and five star miRNAs) in D. latiflorus. One hundred and sixty-two potential targets were identified for the 81 novel bamboo miRNAs. Several targets for the novel miRNAs are transcription factors that play important roles in plant development. Among the novel miRNAs, 30 were selected and their expression profiles in response to different light conditions were validated by qRT-PCR. This study provides the first large-scale cloning and characterization of miRNAs in D. latiflorus. Eighty-four conserved and 81 novel miRNAs were identified in D. latiflorus. Our results present a broad survey of bamboo miRNAs based on experimental and bioinformatics analysis. Although it will be necessary to validate the functions of miRNAs by further experimental research, these results represent a starting point for future research on D. latiflorus and related species.