Data mining for miRNAs and their targets in the Triticeae

The miR319/TaGAMYB3 module regulates plant architecture and improves grain yield in common wheat (Triticum aestivum)

Plant architecture is a key determinant of crop productivity and adaptation. The highly conserved microRNA319 (miR319) family functions in various biological processes, but little is known about how miR319 regulates plant architecture in wheat (Triticum aestivum). Here, we determined that the miR319/TaGAMYB3 module controls plant architecture and grain yield in common wheat. Repressing tae-miR319 using short tandem target mimics resulted in favorable plant architecture traits, including increased plant height, reduced tiller number, enlarged spikes and flag leaves, and thicker culms, as well as enhanced grain yield in field plot tests. Overexpressing tae-miR319 had the opposite effects on plant architecture and grain yield. Although both TaPCF8 and TaGAMYB3 were identified as miR319 target genes, genetic complementation assays demonstrated that only miR319-resistant TaGAMYB3 (rTaGAMYB3) abolished tae-miR319-mediated growth inhibition of flag leaves and spikes. TaGAMYB3 functions as a transcriptional activator of downstream genes, including TaPSKR1, TaXTH23, TaMADS5 and TaMADS51, by binding to their promoters. Furthermore, TaGAMYB3 physically interacts with TaBA1, an important regulator of spike development, to additively activate the transcription of downstream genes such as TaMADS5. Our findings provide insight into how the miR319/TaGAMYB3 module regulates plant architecture and improves grain yield in common wheat.

Combined drought and heat stresses trigger different sets of miRNAs in contrasting potato cultivars

MicroRNAs are small, non-coding RNAs that are responsible for regulation of gene expression during plant growth and development. Although there are many studies on miRNAs in other plants, little work has been done to understand the role of miRNAs in abiotic stress tolerance in potatoes. This study investigates changes in miRNA profiles of two different potato cultivars (tolerant, Unica and susceptible, Russet Burbank) in response to heat, drought and their combination. Transcriptomic studies revealed that miRNA profiles depend on the susceptibility and tolerance of the cultivar and also the stress conditions. Large number of miRNAs were expressed in Unica, whereas Russet Burbank indicated lesser number of changes in miRNA expression. Physiological and transcriptional results clearly supported that Unica cultivar is tolerant to combined drought and heat stress compared to Russet Burbank. Moreover, psRNATarget analysis predicted that major miRNAs identified were targeting genes playing important roles in response to drought and heat stress and their important roles in genetic and post-transcriptional regulation, root development, auxin responses and embryogenesis were also observed. This study focused on eight miRNAs (Novel_8, Novel_9, Novel_105, miR156d-3p, miR160a-5p, miR162a-3p, miR172b-3p and miR398a-5p) and their putative targets where results indicate that they may play a vital role at different post-transcriptional levels against drought and heat stresses. We suggest that miRNA overexpression in plants can lead to increased tolerance against abiotic stresses; furthermore, there should be more emphasis on the studies to investigate the role of miRNAs in combined abiotic stress in plants.

Current status of microRNA-mediated regulation of drought stress responses in cereals

Drought is one of the most important abiotic stress factors impeding crop productivity. With the uncovering of their role as potential regulators of gene expression, microRNAs (miRNAs) have been recognized as new targets for developing stress resistance. MicroRNAs are small noncoding RNAs whose abundance is significantly altered under stress conditions. Interestingly, plant miRNAs predominantly targets transcription factors (TFs), and some of which are also the most critical drought-responsive genes that in turn could regulate the expression of numerous loci with drought-adaptive potential. The phytohormone ABA plays important roles in regulating stomatal conductance and in initiating an adaptive response to drought stress. miRNAs are implicated in regulating ABA-(abscisic acid) and non-ABA-mediated drought resistance pathways. For instance, miR159-MYB module and miR169-NFYA module participates in an ABA-dependent pathway, whereas several other ABA-independent miRNA-target modules (miR156-SPL; miR393-TIR1; miR160-ARF10, ARF16, ARF17; miR167-ARF6 and ARF8; miR390/TAS3siRNA-ARF2, ARF3, ARF4) collectively regulate drought responses in plants. Overall, miRNA-mediated drought response manifests diverse molecular, biochemical and physiological processes. Because of their immense role in controlling gene expression, miRNA manipulation has significant potential to augment plant tolerance to drought stress. This review compiles the current understanding of drought-responsive miRNAs in major cereals. Also, potential miRNA manipulation strategies currently in use along with the challenges and future perspectives are discussed.

Novel approaches on identification of conserved miRNAs for broad-spectrum Potyvirus control measures

Potyviridae comprises more than 200 ssRNA viruses, many of which have a broad host range and geographical distributions. Potyvirids (members of Potyviridae) infect several economically important plants such as saffron, cardamom, cucumber, pepper, potato, tomato, yam, etc. Cumulatively, potyvirids cause a substantial economic loss. The major bottleneck in developing an efficient antiviral strategy is that viruses quickly evade host immunity owing to their higher mutation and recombination rates. Due to this reason, the emergence of newer and improved broad-spectrum approaches to combat viral infections is essential. The use of microRNA's (miRNA) to circumvent viral infection against animal viruses has been successfully employed. Fewer studies reported the development of efficient miRNA-based antivirus resistant strategies against plant viruses and none focused on multiple virus resistance. We focused on potyviruses since studies are limited and identification of conserved miRNAs among various host plants would be an initiative to design broad-spectrum antivirus strategies. In this study, we predicted evolutionarily conserved miRNAs by BLAST searching of reported miRNAs from 15 plants against the GSS and EST sequences of banana. A total of nine miRNAs were predicted and screened in nine diverse potyvirids' hosts (Banana, Tomato, Green gram, Jasmine, Chilli, Coriander, Onion, Rose and Colocasia) belonging to eight different orders (Zingiberales, Solanales, Fabales, Lamiales, Apiales, Asperagales, Rosales and Alismatales). Results suggested that miR168 and miR162 are conserved among all the selected plants. This comprehensive study laid the foundations to design broad-spectrum antivirus resistance using miRNAs. To conclude miR168 and miR162 are conserved among many plants and play a crucial role in evading virus infection which could be used as a potential candidate for developing antiviral strategies against potyvirid infections.

Identification of a redox-dependent regulatory network of miRNAs and their targets in wheat

Reactive oxygen species and antioxidants have an important role in the regulation of plant growth and development under both optimal and stress conditions. In this study, we investigate a possible redox control of miRNAs in wheat (Triticum aestivum ssp. aestivum). Treatment of seedlings with 10 mM H2O2 via the roots for 24 h resulted in decreased glutathione content, increased half-cell reduction potential of the glutathione disulphide/glutathione redox pair, and greater ascorbate peroxidase activity compared to the control plants. These changes were accompanied by alterations in the miRNA transcript profile, with 70 miRNAs being identified with at least 1.5-fold difference in their expression between control and treated (0, 3, 6 h) seedlings. Degradome sequencing identified 86 target genes of these miRNAs, and 6722 possible additional target genes were identified using bioinformatics tools. The H2O2-responsiveness of 1647 target genes over 24 h of treatment was also confirmed by transcriptome analysis, and they were mainly found to be related to the control of redox processes, transcription, and protein phosphorylation and degradation. In a time-course experiment (0-24 h of treatment) a correlation was found between the levels of glutathione, other antioxidants, and the transcript levels of the H2O2-responsive miRNAs and their target mRNAs. This relationship together with bioinformatics modelling of the regulatory network indicated glutathione-related redox control of miRNAs and their targets, which allows the adjustment of the metabolism to changing environmental conditions.

Global gene expression profiling related to temperature-sensitive growth abnormalities in interspecific crosses between tetraploid wheat and Aegilops tauschii

Triploid wheat hybrids between tetraploid wheat and Aegilops tauschii sometimes show abnormal growth phenotypes, and the growth abnormalities inhibit generation of wheat synthetic hexaploids. In type II necrosis, one of the growth abnormalities, necrotic cell death accompanied by marked growth repression occurs only under low temperature conditions. At normal temperature, the type II necrosis lines show grass-clump dwarfism with no necrotic symptoms, excess tillers, severe dwarfism and delayed flowering. Here, we report comparative expression analyses to elucidate the molecular mechanisms of the temperature-dependent phenotypic plasticity in the triploid wheat hybrids. We compared gene and small RNA expression profiles in crown tissues to characterize the temperature-dependent phenotypic plasticity. No up-regulation of defense-related genes was observed under the normal temperature, and down-regulation of wheat APETALA1-like MADS-box genes, considered to act as flowering promoters, was found in the grass-clump dwarf lines. Some microRNAs, including miR156, were up-regulated, whereas the levels of transcripts of the miR156 target genes SPLs, known to inhibit tiller and branch number, were reduced in crown tissues of the grass-clump dwarf lines at the normal temperature. Unusual expression of the miR156/SPLs module could explain the grass-clump dwarf phenotype. Dramatic alteration of gene expression profiles, including miRNA levels, in crown tissues is associated with the temperature-dependent phenotypic plasticity in type II necrosis/grass-clump dwarf wheat hybrids.

Virus-Based MicroRNA Silencing and Overexpressing in Common Wheat (Triticum aestivum L.)

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that arise from large RNA precursors with a stem-loop structure and play important roles in plant development and responses to environmental stress. Although a hundred and nineteen wheat miRNAs have been identified and registered in the miRBase (Release 21.0, June, 2014; http://www.mirbase.org), the functional characterization of these miRNAs in wheat growth and development is lagging due to lack of effective techniques to investigate endogenous miRNA functions in wheat. Here we report barley stripe mosaic virus(BSMV)-based miRNA overexpression and silence systems that can be applied to study miRNA functions in wheat. By utilizing the BSMV system, we successfully knocked down endogenous miR156 and miR166 levels and over-expressed endogenous miR156 and artificial miRNA against phytoene desaturase gene PDS (amiR-PDS) in wheat. amiR-PDS expression caused a great reduction in endogenous mRNA abundance of PDS gene in wheat plant, leading to leaf obviously photobleaching. miR156 silencing led to a great increase in the mRNA level of its target gene SPL2, resulting in a leaf-curl phenotype in wheat seedlings. In contrast, overexpression of miR156 led to a significant reduction in the mRNA level of SPL2 in wheat seedlings, resulting in serious delay of the vegetative phase transitions as well as booting and flowering in wheat. These confirmed that miR156 regulates wheat development and booting time through SPL genes. In summary, the BSMV-based miRNA overexpression and silence systems have extraordinary potential not only for functional study of protein-encoding genes but also for miRNA genes in wheat.

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.

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.

Global Identification of MicroRNAs and Their Targets in Barley under Salinity Stress

Salinity is a major limiting factor for agricultural production worldwide. A better understanding of the mechanisms of salinity stress response will aid efforts to improve plant salt tolerance. In this study, a combination of small RNA and mRNA degradome sequencing was used to identify salinity responsive-miRNAs and their targets in barley. A total of 152 miRNAs belonging to 126 families were identified, of which 44 were found to be salinity responsive with 30 up-regulated and 25 down-regulated respectively. The majority of the salinity-responsive miRNAs were up-regulated at the 8h time point, while down-regulated at the 3h and 27h time points. The targets of these miRNAs were further detected by degradome sequencing coupled with bioinformatics prediction. Finally, qRT-PCR was used to validate the identified miRNA and their targets. Our study systematically investigated the expression profile of miRNA and their targets in barley during salinity stress phase, which can contribute to understanding how miRNAs respond to salinity stress in barley and other cereal crops.

Root precursors of microRNAs in wild emmer and modern wheats show major differences in response to drought stress

MicroRNAs, small regulatory molecules with significant impacts on the transcriptional network of all living organisms, have been the focus of several studies conducted mostly on modern wheat cultivars. In this study, we investigated miRNA repertoires of modern durum wheat and its wild relatives, with differing degrees of drought tolerance, to identify miRNA candidates and their targets involved in drought stress response. Root transcriptomes of Triticum turgidum ssp. durum variety Kızıltan and two Triticum turgidum ssp. dicoccoides genotypes TR39477 and TTD-22 under control and drought conditions were assembled from individual RNA-Seq reads and used for in silico identification of miRNAs. A total of 66 miRNAs were identified from all species, across all conditions, of which 46 and 38 of the miRNAs identified from modern durum wheat and wild genotypes, respectively, had not been previously reported. Genotype- and/or stress-specific miRNAs provide insights into our understanding of the complex drought response. Particularly, miR1435, miR5024, and miR7714, identified only from drought-stress roots of drought-tolerant genotype TR39477, can be candidates for future studies to explore and exploit the drought response to develop tolerant varieties.

Harnessing NGS and Big Data Optimally: Comparison of miRNA Prediction from Assembled versus Non-assembled Sequencing Data--The Case of the Grass Aegilops tauschii Complex Genome

MicroRNAs (miRNAs) are small, endogenous, non-coding RNA molecules that regulate gene expression at the post-transcriptional level. As high-throughput next generation sequencing (NGS) and Big Data rapidly accumulate for various species, efforts for in silico identification of miRNAs intensify. Surprisingly, the effect of the input genomics sequence on the robustness of miRNA prediction was not evaluated in detail to date. In the present study, we performed a homology-based miRNA and isomiRNA prediction of the 5D chromosome of bread wheat progenitor, Aegilops tauschii, using two distinct sequence data sets as input: (1) raw sequence reads obtained from 454-GS FLX Titanium sequencing platform and (2) an assembly constructed from these reads. We also compared this method with a number of available plant sequence datasets. We report here the identification of 62 and 22 miRNAs from raw reads and the assembly, respectively, of which 16 were predicted with high confidence from both datasets. While raw reads promoted sensitivity with the high number of miRNAs predicted, 55% (12 out of 22) of the assembly-based predictions were supported by previous observations, bringing specificity forward compared to the read-based predictions, of which only 37% were supported. Importantly, raw reads could identify several repeat-related miRNAs that could not be detected with the assembly. However, raw reads could not capture 6 miRNAs, for which the stem-loops could only be covered by the relatively longer sequences from the assembly. In summary, the comparison of miRNA datasets obtained by these two strategies revealed that utilization of raw reads, as well as assemblies for in silico prediction, have distinct advantages and disadvantages. Consideration of these important nuances can benefit future miRNA identification efforts in the current age of NGS and Big Data driven life sciences innovation.

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.

Bioinformatic identification and expression analysis of banana microRNAs and their targets

MicroRNAs (miRNAs) represent a class of endogenous non-coding small RNAs that play important roles in multiple biological processes by degrading targeted mRNAs or repressing mRNA translation. Thousands of miRNAs have been identified in many plant species, whereas only a limited number of miRNAs have been predicted in M. acuminata (A genome) and M. balbisiana (B genome). Here, previously known plant miRNAs were BLASTed against the Expressed Sequence Tag (EST) and Genomic Survey Sequence (GSS), a database of banana genes. A total of 32 potential miRNAs belonging to 13 miRNAs families were detected using a range of filtering criteria. 244 miRNA:target pairs were subsequently predicted, most of which encode transcription factors or enzymes that participate in the regulation of development, growth, metabolism, and other physiological processes. In order to validate the predicted miRNAs and the mutual relationship between miRNAs and their target genes, qRT-PCR was applied to detect the tissue-specific expression levels of 12 putative miRNAs and 6 target genes in roots, leaves, flowers, and fruits. This study provides some important information about banana pre-miRNAs, mature miRNAs, and miRNA target genes and these findings can be applied to future research of miRNA functions.

Role of microRNAs in plant drought tolerance

Drought is a normal and recurring climate feature in most parts of the world and plays a major role in limiting crop productivity. However, plants have their own defence systems to cope with adverse climatic conditions. One of these defence mechanisms is the reprogramming of gene expression by microRNAs (miRNAs). miRNAs are small noncoding RNAs of approximately 22 nucleotides length, which have emerged as important regulators of genes at post-transcriptional levels in a range of organisms. Some miRNAs are functionally conserved across plant species and are regulated by drought stress. These properties suggest that miRNA-based genetic modifications have the potential to enhance drought tolerance in cereal crops. This review summarizes the current understanding of the regulatory mechanisms of plant miRNAs, involvement of plant miRNAs in drought stress responses in barley (Hordeum vulgare L.), wheat (Triticum spp.) and other plant species, and the involvement of miRNAs in plant-adaptive mechanisms under drought stress. Potential strategies and directions for future miRNA research and the utilization of miRNAs in the improvement of cereal crops for drought tolerance are also discussed.

Identification of novel miRNAs and miRNA expression profiling in wheat hybrid necrosis

MicroRNAs (miRNAs) play essential roles in a vast array of biological processes, including growth and development, defense against viral infection, and responses to environmental changes in plant. Wheat hybrid necrosis is an interesting genetic phenomenon observed frequency and it is lethal or semi lethal, resulting in gradual death or loss of productivity. However, the molecular basis and mechanisms associated with hybrid necrosis in wheat are still not well understood. Here, we report the population and expression profiles of miRNAs in wheat hybrid necrosis. We identified a total of 57 conserved miRNA families as well as 182 putative novel miRNAs. Expression profiling revealed that expression of 49 known miRNAs and 165 novel miRNAs was changed in hybrid necrosis. And the expression levels of some miRNAs and their predicated targets have been confirmed by qRT-PCR. These results indicate that these miRNAs, especially miR159, miR166, miR167 and miR5072 could be involved in the extensive regulation of gene expression in response to hybrid necrosis.

Novel and conserved heat-responsive microRNAs in wheat (Triticum aestivum L.)

MicroRNAs (miRNAs) are small endogenous RNAs of ~22 nucleotides that have been shown to play regulatory role by negatively affecting the expression of genes at the post-transcriptional level. Information of miRNAs on some important crops like soybean, Arabidopsis, and rice, etc. are available, but no study on heat-responsive novel miRNAs has yet been reported in wheat (Triticum aestivum L.). In the present investigation, a popular wheat cultivar HD2985 was used in small RNA library construction and Illumina HiSeq 2000 was used to perform high-throughput sequencing of the library after cluster generation; 110,896,604 and 87,743,861 reads were generated in the control (22 °C) and heat-treated (42 °C for 2 h) samples, respectively. Forty-four precursor and mature miRNAs were found in T. aestivum from miRBase v 19. The frequencies of the miRNA families varied from 2 (tae-miR1117) to 60,672 (tae-miR159b). We identify 1052 and 902 mature miRNA sequences in HD2985 control and HS-treated samples by mapping on reference draft genome of T. aestivum. Maximum identified miRNAs were located on IWGSC_CSS_3B_scaff (chromosome 3B). We could identify 53 and 46 mature miRNA in the control and HS samples and more than 516 target genes by mapping on the reference genome of Oryza sativa, Zea mays, and Sorghum bicolor. Using different pipelines and plant-specific criteria, 37 novel miRNAs were identified in the control and treated samples. Six novel miRNA were validated using qRT-PCR to be heat-responsive. A negative correlation was, however, observed between the expression of novel miRNAs and their targets. Target prediction and pathway analysis revealed their involvement in the heat stress tolerance. These novel miRNAs are new additions to miRNA database of wheat, and the regulatory network will be made use of in deciphering the mechanism of thermotolerance in wheat.

A chromosome-based draft sequence of the hexaploid bread wheat (Triticum aestivum) genome

An ordered draft sequence of the 17-gigabase hexaploid bread wheat (Triticum aestivum) genome has been produced by sequencing isolated chromosome arms. We have annotated 124,201 gene loci distributed nearly evenly across the homeologous chromosomes and subgenomes. Comparative gene analysis of wheat subgenomes and extant diploid and tetraploid wheat relatives showed that high sequence similarity and structural conservation are retained, with limited gene loss, after polyploidization. However, across the genomes there was evidence of dynamic gene gain, loss, and duplication since the divergence of the wheat lineages. A high degree of transcriptional autonomy and no global dominance was found for the subgenomes. These insights into the genome biology of a polyploid crop provide a springboard for faster gene isolation, rapid genetic marker development, and precise breeding to meet the needs of increasing food demand worldwide.

History and current status of wheat miRNAs using next-generation sequencing and their roles in development and stress

As small molecules that aid in posttranscriptional silencing, microRNA (miRNA) discovery and characterization have vastly benefited from the recent development and widespread application of next-generation sequencing (NGS) technologies. Several miRNAs were identified through sequencing of constructed small RNA libraries, whereas others were predicted by in silico methods using the recently accumulating sequence data. NGS was a major breakthrough in efforts to sequence and dissect the genomes of plants, including bread wheat and its progenitors, which have large, repetitive and complex genomes. Availability of survey sequences of wheat whole genome and its individual chromosomes enabled researchers to predict and assess wheat miRNAs both in the subgenomic and whole genome levels. Moreover, small RNA construction and sequencing-based studies identified several putative development- and stress-related wheat miRNAs, revealing their differential expression patterns in specific developmental stages and/or in response to stress conditions. With the vast amount of wheat miRNAs identified in recent years, we are approaching to an overall knowledge on the wheat miRNA repertoire. In the following years, more comprehensive research in relation to miRNA conservation or divergence across wheat and its close relatives or progenitors should be performed. Results may serve valuable in understanding both the significant roles of species-specific miRNAs and also provide us information in relation to the dynamics between miRNAs and evolution in wheat. Furthermore, putative development- or stress-related miRNAs identified should be subjected to further functional analysis, which may be valuable in efforts to develop wheat with better resistance and/or yield.

Bioinformatic identification and experimental validation of miRNAs from foxtail millet (Setaria italica)

MiRNAs are a novel group of non-coding small RNAs that negatively regulate gene expression. Many miRNAs have been identified and investigated extensively in plant species with sequenced genomes. However, few miRNAs have been identified in foxtail millet (Setaria italica), which is an ancient cereal crop of great importance for dry land agriculture. In this study, 271 foxtail millet miRNAs belonging to 44 families were identified using a bioinformatics approach. Twenty-three pairs of sense/antisense miRNAs belonging to 13 families, and 18 miRNA clusters containing members of 8 families were discovered in foxtail millet. We identified 432 potential targets for 38 miRNA families, most of which were predicted to be involved in plant development, signal transduction, metabolic pathways, disease resistance, and environmental stress responses. Gene ontology (GO) analysis revealed that 101, 56, and 23 target genes were involved in molecular functions, biological processes, and cellular components, respectively. We investigated the expression patterns of 43 selected miRNAs using qRT-PCR analysis. All of the miRNAs were expressed ubiquitously with many exhibiting different expression levels in different tissues. We validated five predicted targets of four miRNAs using the RNA ligase mediated rapid amplification of cDNA end (5'-RLM-RACE) method.

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.

Identification and characterization of a subset of microRNAs in wheat (Triticum aestivum L.)

MicroRNAs (miRNAs) represent a class of endogenous regulator for post-transcriptionally modulating gene expression. Elucidating complete miRNA repertoires for individual species is a long-desired goal in miRNA research. So far only 42 have been annotated for common wheat (Triticum aestivum) due to its large genome. Here, we employed miRDeep-P, a program developed previously for retrieving miRNAs from deep sequencing data in plants, to parse 14 sequenced small RNA libraries of wheat using expression sequence tags (ESTs) as the reference in lieu of a complete genome sequence. This effort identified 145 miRNAs along with the corresponding stem-looped precursors with many differentially expressed in development and associated with powdery mildew infection. Examination of the phylogenetic distribution of these miRNAs and their target genes revealed that many are conserved in monocots. The set of miRNAs identified in this study is thus useful to further miRNA research in wheat and other important cereal crop species.

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.

Unique and conserved microRNAs in wheat chromosome 5D revealed by next-generation sequencing

MicroRNAs are a class of short, non-coding, single-stranded RNAs that act as post-transcriptional regulators in gene expression. miRNA analysis of Triticum aestivum chromosome 5D was performed on 454 GS FLX Titanium sequences of flow-sorted chromosome 5D with a total of 3,208,630 good quality reads representing 1.34x and 1.61x coverage of the short (5DS) and long (5DL) arms of the chromosome respectively. In silico and structural analyses revealed a total of 55 miRNAs; 48 and 42 miRNAs were found to be present on 5DL and 5DS respectively, of which 35 were common to both chromosome arms, while 13 miRNAs were specific to 5DL and 7 miRNAs were specific to 5DS. In total, 14 of the predicted miRNAs were identified in wheat for the first time. Representation (the copy number of each miRNA) was also found to be higher in 5DL (1,949) compared to 5DS (1,191). Targets were predicted for each miRNA, while expression analysis gave evidence of expression for 6 out of 55 miRNAs. Occurrences of the same miRNAs were also found in Brachypodium distachyon and Oryza sativa genome sequences to identify syntenic miRNA coding sequences. Based on this analysis, two other miRNAs: miR1133 and miR167 were detected in B. distachyon syntenic region of wheat 5DS. Five of the predicted miRNA coding regions (miR6220, miR5070, miR169, miR5085, miR2118) were experimentally verified to be located to the 5D chromosome and three of them : miR2118, miR169 and miR5085, were shown to be 5D specific. Furthermore miR2118 was shown to be expressed in Chinese Spring adult leaves. miRNA genes identified in this study will expand our understanding of gene regulation in bread wheat.

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

Wheat is the most highly cultivated plant species for its grain production throughout the world. Because small RNA-dependent gene regulation is critical for successful completion of plant life cycle including its productivity, identification of not only miRNAs but also confirming their targets in wheat is important. To identify small RNAs including novel miRNAs as well as miRNA targets in wheat, we constructed small RNA and degradome libraries from wheat seedlings. Small RNA analysis resulted in identification of most conserved miRNAs including novel miRNAs that can be grouped into 32 miRNA families. The sequence analysis also led to the characterization of two abundantly expressed rRNA-derived small RNAs. To identify miRNA targets, degradome library was sequenced and the bioinformatic analysis confirmed 53 genes as targets for miRNAs and Tas3-siRNAs. Degradome analysis also confirmed a conserved fine-tuning mechanism of Tas3-siRNA abundance by siRNA-mediated silencing of TAS3 transcripts in diverse plant species. These findings added additional information to the small RNA knowledge-base in wheat.

Boron stress responsive microRNAs and their targets in barley

Boron stress is an environmental factor affecting plant development and production. Recently, microRNAs (miRNAs) have been found to be involved in several plant processes such as growth regulation and stress responses. In this study, miRNAs associated with boron stress were identified and characterized in barley. miRNA profiles were also comparatively analyzed between root and leave samples. A total of 31 known and 3 new miRNAs were identified in barley; 25 of them were found to respond to boron treatment. Several miRNAs were expressed in a tissue specific manner; for example, miR156d, miR171a, miR397, and miR444a were only detected in leaves. Additionally, a total of 934 barley transcripts were found to be specifically targeted and degraded by miRNAs. In silico analysis of miRNA target genes demonstrated that many miRNA targets are conserved transcription factors such as Squamosa promoter-binding protein, Auxin response factor (ARF), and the MYB transcription factor family. A majority of these targets were responsible for plant growth and response to environmental changes. We also propose that some of the miRNAs in barley such as miRNA408 might play critical roles against boron exposure. In conclusion, barley may use several pathways and cellular processes targeted by miRNAs to cope with boron stress.

Developmentally regulated expression and complex processing of barley pri-microRNAs

BACKGROUND

MicroRNAs (miRNAs) regulate gene expression via mRNA cleavage or translation inhibition. In spite of barley being a cereal of great economic importance, very little data is available concerning its miRNA biogenesis. There are 69 barley miRNA and 67 pre-miRNA sequences available in the miRBase (release 19). However, no barley pri-miRNA and MIR gene structures have been shown experimentally. In the present paper, we examine the biogenesis of selected barley miRNAs and the developmental regulation of their pri-miRNA processing to learn more about miRNA maturation in barely.

RESULTS

To investigate the organization of barley microRNA genes, nine microRNAs - 156g, 159b, 166n, 168a-5p/168a-3p, 171e, 397b-3p, 1120, and 1126 - were selected. Two of the studied miRNAs originate from one MIR168a-5p/168a-3p gene. The presence of all miRNAs was confirmed using a Northern blot approach. The miRNAs are encoded by genes with diverse organizations, representing mostly independent transcription units with or without introns. The intron-containing miRNA transcripts undergo complex splicing events to generate various spliced isoforms. We identified miRNAs that were encoded within introns of the noncoding genes MIR156g and MIR1126. Interestingly, the intron that encodes miR156g is spliced less efficiently than the intron encoding miR1126 from their specific precursors. miR397b-3p was detected in barley as a most probable functional miRNA, in contrast to rice where it has been identified as a complementary partner miRNA*. In the case of miR168a-5p/168a-3p, we found the generation of stable, mature molecules from both pre-miRNA arms, confirming evolutionary conservation of the stability of both species, as shown in rice and maize. We suggest that miR1120, located within the 3' UTR of a protein-coding gene and described as a functional miRNA in wheat, may represent a siRNA generated from a mariner-like transposable element.

CONCLUSIONS

Seven of the eight barley miRNA genes characterized in this study contain introns with their respective transcripts undergoing developmentally specific processing events prior to the dicing out of pre-miRNA species from their pri-miRNA precursors. The observed tendency to maintain the intron encoding miR156g within the transcript, and preferences in splicing the miR1126-harboring intron, may suggest the existence of specific regulation of the levels of intron-derived miRNAs in barley.

Sorting the wheat from the chaff: identifying miRNAs in genomic survey sequences of Triticum aestivum chromosome 1AL

Individual chromosome-based studies of bread wheat are beginning to provide valuable structural and functional information about one of the world's most important crops. As new genome sequences become available, identifying miRNA coding sequences is arguably as important a task as annotating protein coding sequences, but one that is not as well developed. We compared conservation-based identification of conserved miRNAs in 1.5× coverage survey sequences of wheat chromosome 1AL with a predictive method based on pre-miRNA hairpin structure alone. In total, 42 sequences expected to encode conserved miRNAs were identified on chromosome 1AL, including members of several miRNA families that have not previously been reported to be expressed in T. aestivum. In addition, we demonstrate that a number of sequences previously annotated as novel wheat miRNAs are closely related to transposable elements, particularly Miniature Inverted Terminal repeat Elements (MITEs). Some of these TE-miRNAs may well have a functional role, but separating true miRNA coding sequences from TEs in genomic sequences is far from straightforward. We propose a strategy for annotation to minimize the risk of mis-identifying TE sequences as miRNAs.

Identification and characterization of the microRNAs and their targets in Salmo salar

MicroRNAs (miRNAs) are small, non-coding and regulatory RNAs about 18 to 26 nucleotides long. Their conserved nature among the various organisms makes them a good source of new miRNAs discovery by comparative genomics approach. The study resulted in novel 75 precursor miRNAs containing 102 mature sequences belonging to 46 families in an important aquatic environmental monitoring fish (Salmo salar). All the miRNA families (let-7, mir-1, 7, 9, 21, 22, 92, 96, 122, 126, 128, 129, 132, 133, 142, 144, 147, 148, 196, 202, 212, 223, 375, 429, 430, 449, 451, 457, 466, 682, 700, 1388, 1594, 1600, 1607, 1616, 1642, 1681, 1701, 1720, 1772, 1782, 1787, 1814, 2189 and 3540) are found for the first time in S. salar. All 75 miRNA precursors form stable minimum free energy stem loop and the mature miRNAs reside in the stem portion of the stem loop structure. Their target proteins are involved in transcription factors (28%), metabolism (23%), signaling (18%), transportation (9%), immunity (8%), stress related activity (5%), cancer and tumor related activity (5%), growth and development (3%), and cell division (1%).

Computational identification of miRNAs in medicinal plant Senecio vulgaris (Groundsel)

RNAs Interference plays a very important role in gene silencing. In vitro identification of miRNAs is a slow process as it is difficult to isolate them. Nucleotide sequences of miRNAs are highly conserved among the plants and, this form the key feature behind the identification of miRNAs in plant species by homology alignment. In silico identification of miRNAs from EST database is emerging as a novel, faster and reliable approach. Here EST sequences of Senecio vulgaris (Groundsel) were searched against known miRNA sequences by using BLASTN tool. A total of 10 miRNAs were identified from 1956 EST sequences and 115 GSS sequences. The most stable miRNA identified is svu-mir-1. This approach will accelerate advance research in regulation of gene expression in Groundsel by interfering RNAs.

miRNA expression patterns of Triticum dicoccoides in response to shock drought stress

Drought is a major environmental stress factor that affects plant growth and development worldwide. Wild emmer wheat (Triticum turgidum ssp. dicoccoides), the ancestor of domesticated durum wheat (Triticum turgidum ssp. durum), has great potential for improving the understanding of the wheat drought response. MicroRNAs (miRNAs) are a recently discovered class of gene expression regulators that have also been linked to several plant stress responses; however, this relationship is just beginning to be understood. miRNA expression patterns of drought-resistant wild emmer wheat in response to drought stress were investigated using a plant miRNA microarray platform. Expression was detected to be 205 miRNAs in control and 438 miRNAs in drought-stressed leaf and root tissues. Of these miRNAs, the following 13 were differentially regulated in response to drought: miR1867, miR896, miR398, miR528, miR474, miR1450, miR396, miR1881, miR894, miR156, miR1432, miR166 and miR171. Regulation of miRNAs upon 4 and 8 h drought stress applications observed by qRT-PCR. Target transcripts of differentially regulated miRNAs were computationally predicted. In addition to miRNA microarray study, five new conserved T. turgidum miRNAs were identified through a homology-based approach, and their secondary structures and putative targets were predicted. These findings both computationally and experimentally highlight the presence of miRNAs in T. dicoccoides and further extend the role of miRNAs under shock drought stress conditions.

Discovery of barley miRNAs through deep sequencing of short reads

Background

MicroRNAs are important components of the regulatory network of biological systems and thousands have been discovered in both animals and plants. Systematic investigations performed in species with sequenced genomes such as Arabidopsis, rice, poplar and Brachypodium have provided insights into the evolutionary relationships of this class of small RNAs among plants. However, miRNAs from barley, one of the most important cereal crops, remain unknown.

Results

We performed a large scale study of barley miRNAs through deep sequencing of small RNAs extracted from leaves of two barley cultivars. By using the presence of miRNA precursor sequences in related genomes as one of a number of supporting criteria, we identified up to 100 miRNAs in barley. Of these only 56 have orthologs in wheat, rice or Brachypodium that are known to be expressed, while up to 44 appear to be specifically expressed in barley.

Conclusions

Our study, the first large scale investigation of small RNAs in barley, has identified up to 100 miRNAs. We demonstrate that reliable identification of miRNAs via deep sequencing in a species whose genome has not been sequenced requires a more careful analysis of sequencing errors than is commonly performed. We devised a read filtering procedure for dealing with errors. In addition, we found that the use of a large dataset of almost 35 million reads permits the use of read abundance distributions along putative precursor sequences as a practical tool for isolating miRNAs in a large background of reads originating from other non-coding and coding RNAs. This study therefore provides a generic approach for discovering novel miRNAs where no genome sequence is available.

Identification of plant microRNAs using expressed sequence tag analysis

microRNAs (miRNAs) are a new class of small endogenous noncoding regulatory RNAs, which play an important function in plant growth, development, phase change, and response to environmental stress. Identifying miRNAs is the first step for investigating miRNA-mediated gene regulation and miRNA function. In this chapter, we describe a comprehensive comparative genomics-based expressed sequence tag (EST) analysis for identifying miRNAs from a wide range of plant species. EST analysis is based on the conservation of miRNA sequences and the stem-loop hairpin secondary structures of miRNAs. In this method, potential miRNAs will first be identified by EST analysis followed by confirmation using TaqMan(®) MicroRNA qRT-PCR. This method is simple and reliable with high efficiency. This method has also been widely adopted by many scientists around the world and several hundreds of miRNAs have been identified in many plant species using this method.

A computational-based update on microRNAs and their targets in barley (Hordeum vulgare L.)

Background

Many plant species have been investigated in the last years for the identification and characterization of the corresponding miRNAs, nevertheless extensive studies are not yet available on barley (at the time of this writing). To extend and to update information on miRNAs and their targets in barley and to identify candidate polymorphisms at miRNA target sites, the features of previously known plant miRNAs have been used to systematically search for barley miRNA homologues and targets in the publicly available ESTs database. Matching sequences have then been related to Unigene clusters on which most of this study was based.

Results

One hundred-fifty-six microRNA mature sequences belonging to 50 miRNA families have been found to significantly match at least one EST sequence in barley. As expected on the basis of phylogenetic relations, miRNAs putatively orthologous to those of Triticum are significantly over-represented inside the set of identified barley microRNA mature sequences. Many previously known and several putatively new miRNA/target pairs have been identified. When the predicted microRNA targets were grouped into functional categories, biological processes previously known to be regulated by miRNAs, such as development and response to biotic and abiotic stress, have been highlighted and most of the target molecular functions were related to transcription regulation. Candidate microRNA coding genes have been reported and genetic variation (SNPs/indels) both in functional regions of putative miRNAs (mature sequence) and at miRNA target sites has been found.

Conclusions

This study has provided an update of the information on barley miRNAs and their targets representing a foundation for future studies. Many of previously known plant microRNAs have homologues in barley with expected important roles during development, nutrient deprivation, biotic and abiotic stress response and other important physiological processes. Putative polymorphisms at miRNA target sites have been identified and they can represent an interesting source for the identification of functional genetic variability.

Computational identification of miRNAs and their target genes from expressed sequence tags of tea (Camellia sinensis)

MicroRNAs (miRNAs) are a newly identified class of small non-protein-coding post-transcriptional regulatory RNA in both plants and animals. The use of computational homology based search for expressed sequence tags (ESTs) with the Ambros empirical formula and other structural feature criteria filter is a suitable combination towards the discovery and isolation of conserved miRNAs from tea and other plant species whose genomes are not yet sequenced. In the present study, we blasted the database of tea (Camellia sinensis) ESTs to search for potential miRNAs, using previously known plant miRNAs. For the first time, four candidate miRNAs from four families were identified in tea. Using the newly identified miRNA sequences, a total of 30 potential target genes were identified for 11 miRNA families; 6 of these predicted target genes encode transcription factors (20%), 16 target genes appear to play roles in diverse physiological processes (53%) and 8 target genes have hypothetical or unknown functions (27%). These findings considerably broaden the scope of understanding the functions of miRNA in tea.

Regulation of barley miRNAs upon dehydration stress correlated with target gene expression

We aim to identify conserved and dehydration responsive microRNAs (miRNAs) in Hordeum vulgare (barley). A total of 28 new barley miRNAs belonging to 18 distinct miRNA families were identified. Detailed nucleotide analyses revealed that barley pre-miRNAs are in the range of 46-114 nucleotides with average of 77.14. Using 28 newly detected miRNAs as queries, 445 potential target mRNAs were predicted. The predicted miRNAs were differentially expressed and some of them behaved similarly in leaf and root tissues upon stress treatment. Hvu-MIR156, Hvu-MIR166, Hvu-MIR171, and Hvu-MIR408 were detected as dehydration stress-responsive barley miRNAs. To discover target transcripts of barley miRNAs a modified 5' RLM-RACE was performed and seven cleaved miRNA transcripts were retrieved from drought stressed leaf samples. In silico analysis indicated 15 potential EST targets. Measurement of expression levels showed a positive correlation between levels of miRNA expression and suppression of their target mRNA transcripts in dehydration-stress-treated barley.

Identification and characterization of conserved microRNAs and their target genes in wheat (Triticum aestivum)

MicroRNAs (miRNAs) are non-coding small RNAs that regulate gene expression by translational repression or transcript degradation. A large number of miRNAs have been identified from model plant species; however, the character of conserved miRNAs is poorly understood. We studied 42 miRNA families that are conserved within the plant kingdom, using the miRBase database. Some conserved miRNA families were found to be preferentially expressed in dicots relative to monocots, especially miR403, miR472 and miR479. Using an improved homology search-based approach and the conserved miRNAs as the query set, 34 conserved miRNAs and the miR482 family were identified in wheat. Forty-six wheat mRNAs were predicted as their putative target genes. Most conserved wheat miRNAs were found to retain homologous target interactions and have analogous molecular functions. The miR172 displayed a wheat-specific function and was found to have an additional target interaction with succinyl-CoA ligase. We concluded that although miRNAs are conserved, the expression and function of some have drifted during long periods of plant evolution.

Computational identification of microRNAs and their targets in wheat (Triticum aestivum L.)

microRNAs (miRNAs) are a class of endogenous, non-coding, short (approximately 21 nt) RNAs directly involved in regulating gene expression at the post-transcriptional level. Previous reports have noted that plant miRNAs in the plant kingdom are highly conserved, which provides the foundation for identification of conserved miRNAs in other plant species through homology alignment. Conserved miRNAs in wheat are identified using EST (Expressed Sequence Tags) and GSS analysis. All previously known miRNAs in other plant species were blasted against wheat EST and GSS sequences to select novel miRNAs in wheat by a series of filtering criteria. From a total of 37 conserved miRNAs belonging to 18 miRNA families 10 conserved miRNAs comprising 4 families were reported in wheat. MiR395 is found to be a special family, because three members belonging to the same miR395 family are clustered together, similar to animal miRNAs. MiRNA targets are transcription factors involved in wheat growth and development, metabolism,and stress responses.