Stress responsive miRNAs and isomiRs in cereals

Deep-sequence profiling of miRNAs and their target prediction in Monotropa hypopitys

Myco-heterotroph Monotropa hypopitys is a widely spread perennial herb used to study symbiotic interactions and physiological mechanisms underlying the development of non-photosynthetic plant. Here, we performed, for the first time, transcriptome-wide characterization of M. hypopitys miRNA profile using high throughput Illumina sequencing. As a result of small RNA library sequencing and bioinformatic analysis, we identified 55 members belonging to 40 families of known miRNAs and 17 putative novel miRNAs unique for M. hypopitys. Computational screening revealed 206 potential mRNA targets for known miRNAs and 31 potential mRNA targets for novel miRNAs. The predicted target genes were described in Gene Ontology terms and were found to be involved in a broad range of metabolic and regulatory pathways. The identification of novel M. hypopitys-specific miRNAs, some with few target genes and low abundances, suggests their recent evolutionary origin and participation in highly specialized regulatory mechanisms fundamental for non-photosynthetic biology of M. hypopitys. This global analysis of miRNAs and their potential targets in M. hypopitys provides a framework for further investigation of miRNA role in the evolution and establishment of non-photosynthetic myco-heterotrophs.

High-throughput sequencing reveals differential regulation of miRNAs in fenoxaprop-P-ethyl-resistant Beckmannia syzigachne

Non-target site resistance (NTSR) to herbicides is an increasing concern for weed control. The majority of previous studies have focused on metabolic resistance mechanisms of NTSR, but no research exists on gene regulation mechanisms behind herbicide resistance, such as microRNA (miRNA). Here, we identified 3 American sloughgrass (Beckmannia syzigachne Steud.) populations containing fenoxaprop-P-ethyl-resistant plants. We then constructed small RNA libraries and subjected them to deep sequencing and bioinformatics analyses. Forty known and 36 potentially novel, predicted miRNAs were successfully identified. Of these, we identified 3 conserved, predicted candidate NTSR-determinant miRNAs and their potential corresponding target genes, as well as 4 novel potential miRNAs with high count. Target gene prediction and annotation indicated that these 7 differentially expressed miRNAs potentially play a role in regulating specific stress-responsive genes, very likely related to herbicide resistance. Expression profiles were determined with quantitative real-time PCR. The present study is a novel, large-scale characterization of weed miRNAs. The results should further our understanding of miRNA expression profiles associated with herbicide resistance, allowing for the development of more effective weed management strategies.

In silico search and biological validation of microRNAs related to drought response in peach and almond

Plant responses to drought stress are regulated at the transcriptional and post-transcriptional levels through noncoding endogenous microRNAs. These microRNAs play key roles in gene expression, mainly by down-regulating target mRNAs. In this work, an in silico search and validation for microRNAs related to drought response in peach ('G.H. Hill'), almond ('Sefied') and an interspecific peach-almond hybrid ('GN 15') has been performed. We used qPCR to analyse the gene expression of several miRNAs described as being related to drought response in peach, including miR156, miR159, miR160, miR167, miR171, miR172, miR398, miR403, miR408, miR842 and miR2275 under mild and severe water deficit. These miRNAs were in silico selected on the basis of previous works, their conservation in plants and their drought response. qPCR analysis confirmed the implication of these miRNAs in the dehydration stress response in the three assayed genotypes. Comparison of miRNA expression patterns in the three evaluated genotypes indicated that the hybrid GN 15 showed higher expression levels of specific miRNAs which should be related to the observed drought tolerance. mRNA target transcripts of the miRNAs studied were predicted using the Rose database, which includes transcription factors that regulate plant growth and development. In addition, results showed that the promoter region contains responsive elements to hormone-mediated regulatory elements. Network analysis not only unravelled the interaction between miRNAs and their predicted gene targets but also highlighted the roles of miRNAs in response to drought stress.

Characterization of wheat miRNAs and their target genes responsive to cadmium stress

A increasing number of microRNAs have been shown to play important regulatory roles in plant responses to various metal stresses. However, little information about miRNAs especially miRNAs responsive to cadmium (Cd) stress is available in wheat. To investigate the role of miRNAs in responses to Cd stress, wheat seedlings were subjected to 250 μM Cd solution for 6, 12, 24 and 48 h, and analyses of morphological and physiological changes as well as the expression of five miRNAs and their corresponding targets were carried out. Our results demonstrated that miRNAs and their targets were differentially expressed in leaves and roots of wheat seedlings exposed to Cd stress. Furthermore, miR398 may involve in oxidative stress tolerance by regulating its target CSD to participate in Cd stress. Among ten miRNA-target pairs studied, nine pairs showed complex regulation relationship in leaves and roots of wheat seedlings exposed to Cd stress. These findings suggested that miRNAs are involved in the mediation of Cd stress signaling responses in wheat. The characterization of the miRNAs and the associated targets in responses to Cd exposure provides a framework for understanding the molecular mechanism of heavy metal tolerance in plants.

Wheat miRNA ancestors: evident by transcriptome analysis of A, B, and D genome donors

MicroRNAs are critical players of post-transcriptional gene regulation with profound effects on the fundamental processes of cellular life. Their identification and characterization, together with their targets, hold great significance in exploring and exploiting their roles on a functional context, providing valuable clues into the regulation of important biological processes, such as stress tolerance or environmental adaptation. Wheat is a hardy crop, extensively harvested in temperate regions, and is a major component of the human diet. With the advent of the next generation sequencing technologies considerably decreasing sequencing costs per base-pair, genomic, and transcriptomic data from several wheat species, including the progenitors and wild relatives have become available. In this study, we performed in silico identification and comparative analysis of microRNA repertoires of bread wheat (Triticum aestivum L.) and its diploid progenitors and relatives, Aegilops sharonensis, Aegilops speltoides, Aegilops tauschii, Triticum monococcum, and Triticum urartu through the utilization of publicly available transcriptomic data. Over 200 miRNA families were identified, majority of which have not previously been reported. Ancestral relationships expanded our understanding of wheat miRNA evolution, while T. monococcum miRNAs delivered important clues on the effects of domestication on miRNA expression. Comparative analyses on wild Ae. sharonensis accessions highlighted candidate miRNAs that can be linked to stress tolerance. The miRNA repertoires of bread wheat and its diploid progenitors and relatives provide important insight into the diversification and distribution of miRNA genes, which should contribute to the elucidation of miRNA evolution of Poaceae family. A thorough understanding of the convergent and divergent expression profiles of miRNAs in different genetic backgrounds can provide unique opportunities to modulation of gene regulation for better crop performance.

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.

Topological characteristics of target genes regulated by abiotic-stress-responsible miRNAs in a rice interactome network

A great number of microRNAs (miRNAs) have been identified in responding and acting in gene regulatory networks associated with plant tolerance to abiotic stress conditions, such as drought, salinity, and high temperature. The topological exploration of target genes regulated by abiotic-stress-responsible miRNAs (ASRmiRs) in a network facilitates to discover the molecular basis of plant abiotic stress response. This study was based on the staple food rice (Oryza sativa) in which ASRmiRs were manually curated. After having compared the topological properties of target genes (stress-miR-targets) with those (non-stress-miR-targets) not regulated by ASRmiRs in a rice interactome network, we found that stress-miR-targets exhibited distinguishable topological properties. The interaction probability analysis and k-core decomposition showed that stress-miR-targets preferentially interacted with non-stress-miR-targets and located at the peripheral positions in the network. Our results indicated an obvious topological distinction between the two types of genes, reflecting the specific mechanisms of action of stress-miR-targets in rice abiotic stress response. Also, the results may provide valuable clues to elucidate molecular mechanisms of crop response to abiotic stress.

OMICS Technologies and Applications in Sugar Beet

Sugar beet is a species of the Chenopodiaceae family. It is an important sugar crop that supplies approximately 35% of the sugar in the world. Sugar beet M14 line is a unique germplasm that contains genetic materials from Beta vulgaris L. and Beta corolliflora Zoss. And exhibits tolerance to salt stress. In this review, we have summarized OMICS technologies and applications in sugar beet including M14 for identification of novel genes, proteins related to biotic and abiotic stresses, apomixes and metabolites related to energy and food. An OMICS overview for the discovery of novel genes, proteins and metabolites in sugar beet has helped us understand the complex mechanisms underlying many processes such as apomixes, tolerance to biotic and abiotic stresses. The knowledge gained is valuable for improving the tolerance of sugar beet and other crops to biotic and abiotic stresses as well as for enhancing the yield of sugar beet for energy and food production.

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.

A Comprehensive Prescription for Plant miRNA Identification

microRNAs (miRNAs) are tiny ribo-regulatory molecules involved in various essential pathways for persistence of cellular life, such as development, environmental adaptation, and stress response. In recent years, miRNAs have become a major focus in molecular biology because of their functional and diagnostic importance. This interest in miRNA research has resulted in the development of many specific software and pipelines for the identification of miRNAs and their specific targets, which is the key for the elucidation of miRNA-modulated gene expression. While the well-recognized importance of miRNAs in clinical research pushed the emergence of many useful computational identification approaches in animals, available software and pipelines are fewer for plants. Additionally, existing approaches suffers from mis-identification and annotation of plant miRNAs since the miRNA mining process for plants is highly prone to false-positives, particularly in cereals which have a highly repetitive genome. Our group developed a homology-based in silico miRNA identification approach for plants, which utilizes two Perl scripts "SUmirFind" and "SUmirFold" and since then, this method helped identify many miRNAs particularly from crop species such as Triticum or Aegliops. Herein, we describe a comprehensive updated guideline by the implementation of two new scripts, "SUmirPredictor" and "SUmirLocator," and refinements to our previous method in order to identify genuine miRNAs with increased sensitivity in consideration of miRNA identification problems in plants. Recent updates enable our method to provide more reliable and precise results in an automated fashion in addition to solutions for elimination of most false-positive predictions, miRNA naming and miRNA mis-annotation. It also provides a comprehensive view to genome/transcriptome-wide location of miRNA precursors as well as their association with transposable elements. The "SUmirPredictor" and "SUmirLocator" scripts are freely available together with a reference high-confidence plant miRNA list.

Genome-Wide Identification of MicroRNAs in Leaves and the Developing Head of Four Durum Genotypes during Water Deficit Stress

MicroRNAs (miRNAs) are small non-coding RNAs that play critical roles in plant development and abiotic stress responses. The miRNA transcriptome (miRNAome) under water deficit stress has been investigated in many plant species, but is poorly characterised in durum wheat (Triticum turgidum L. ssp. durum). Water stress during early reproductive stages can result in significant yield loss in durum wheat and this study describes genotypic differences in the miRNAome between water deficit tolerant and sensitive durum genotypes. Small RNA libraries (96 in total) were constructed from flag leaf and developing head tissues of four durum genotypes, with or without water stress to identify differentially abundant miRNAs. Illumina sequencing detected 110 conserved miRNAs and 159 novel candidate miRNA hairpins with 66 conserved miRNAs and five novel miRNA hairpins differentially abundant under water deficit stress. Ten miRNAs (seven conserved, three novel) were validated through qPCR. Several conserved and novel miRNAs showed unambiguous inverted regulatory profiles between the durum genotypes. Several miRNAs also showed differential abundance between two tissue types regardless of treatment. Predicted mRNA targets (130) of four novel durum miRNAs were characterised using Gene Ontology (GO) which revealed functions common to stress responses and plant development. Negative correlation was observed between several target genes and the corresponding miRNA under water stress. For the first time, we present a comprehensive study of the durum miRNAome under water deficit stress. The identification of differentially abundant miRNAs provides molecular evidence that miRNAs are potential determinants of water stress tolerance in durum wheat. GO analysis of predicted targets contributes to the understanding of genotypic physiological responses leading to stress tolerance capacity. Further functional analysis of specific stress responsive miRNAs and their interaction with targets is ongoing and will assist in developing future durum wheat varieties with enhanced water deficit stress tolerance.

Genome-wide identification of abiotic stress-regulated and novel microRNAs in mulberry leaf

As the most important food plant for sericultural industry, mulberry trees have to suffer from a wide range of abiotic and biotic stresses, such as drought and high salinity. MicroRNAs (miRNAs) have been proved to play important roles in abiotic stresses regulation in many plants. However, there are seldom reports on the miRNAs expression profiles upon abiotic challenges in mulberry. In this study, three small RNA libraries from mulberry leaf tissue with or without drought or salt treatment were constructed and deep sequenced. Total of 48 conserved miRNAs (including miRNA*) and 162 novel miRNAs were identified (processing precision value>0.1). A total of 270 and 1963 target genes were predicted for conserved miRNAs and novel miRNAs, respectively. 13 differentially expressed miRNAs were detected under drought or salt stresses by deep sequencing and qRT-PCR. 5' RLM-RACE validated Morus 013341 to be the target gene of miR-395a. Our results provided initial clue to further study molecular mechanism on abiotic stresses regulation in mulberry.

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.

MicroRNA nomenclature and the need for a revised naming prescription

A central environment and interface for microRNA (miRNA) registry and repository and a general standardized framework for their systematic annotation was established over a decade ago. However, the numbers of experimentally and computationally identified miRNAs are swiftly accumulating, and new aspects of miRNA-mediated gene regulation are being revealed. Currently, it is of great significance that the annotation framework should be redefined to include newly discovered miRNA species such as the variants of mature miRNAs (isomiRNAs), and organellar miRNAs: cipomiRNAs and mitomiRNAs. It is also of great importance that key terminology referring to the novelty, evolutionary history or biogenesis of miRNAs, as well as the confidence of their identification are standardized in the literature and disseminated in a central miRNA registry. Here, we review the status of miRNA nomenclature, curation and critical points of need for a revision of miRNA nomenclature and terminology.

Plant miRNAs: biogenesis, organization and origins

MicroRNAs, or miRNAs, are posttranscriptional regulators of gene expression. A wealth of observations and findings suggest highly complex, multicomponent, and intermingled pathways governing miRNA biogenesis and miRNA-mediated gene silencing. Plant miRNA genes are usually found as individual entities scattered around the intergenic and-to a much lesser extent-intragenic space, while miRNA gene clusters, formed by tandem or segmental duplications, also exist in plant genomes. Genome duplications are proposed to contribute to miRNA family expansions, as well. Evolutionarily young miRNAs retaining extensive homology to their loci of origin deliver important clues into miRNA origins and evolution. Additionally, imprecisely processed miRNAs evidence noncanonical routes of biogenesis, which may affect miRNA expression levels or targeting capabilities. Majority of the knowledge regarding miRNAs comes from model plant species. As ongoing research progressively expands into nonmodel systems, our understanding of miRNAs and miRNA-related pathways changes which opens up new perspectives and frontiers in miRNA research.

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.

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.