New wheat microRNA using whole-genome sequence

Computational screening of miRNAs and their targets in saffron (Crocus sativus L.) by transcriptome mining

A robust workflow for the identification of miRNAs and their targets in saffron was developed. MicroRNA-mediated gene regulation in saffron is potentially involved in several biological processes, including the biosynthesis of highly valuable apocarotenoids. Saffron (Crocus sativus L.) is the most expensive spice in the world and a major source of apocarotenoids. Even though miRNAs (20-24 nt non-coding small RNAs) are important regulators of gene expression at transcriptional and post-transcriptional levels, their role in saffron has not been thoroughly investigated. As a result, a workflow for computational identification of miRNAs and their targets can be useful to uncover the regulatory networks underlying biological processes in this valuable plant. The efficiency of several assembly tools such as Trans-ABySS, Trinity, Bridger, rnaSPAdes, and EvidentialGene was evaluated based on both reference-based and reference-free metrics using transcriptome data. A reliable workflow for computational identification of miRNAs and their targets in saffron was described. The EvidentialGene was found to be the most efficient de novo transcriptome assembler for saffron as a complex triploid model, followed by the Trinity. In total, 66 miRNAs from 19 different families that target 2880 genes, including several transcription factors involved in the flowering transition, were identified. Three of the identified targets were involved in the terpenoids backbone biosynthesis. CsCCD and CsUGT genes involved in the apocarotenoids biosynthetic pathway were targeted by csa-miR156g and csa-miR156b-3p, revealing a unique post-transcriptional regulation dynamic in saffron. The identified miRNAs and their targets add to our understanding of the many biological roles of miRNAs in saffron and shed new light on the control of the apocarotenoid biosynthetic pathway in this valuable plant.

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.

Pan-Genome miRNomics in Brachypodium

Pan-genomes are efficient tools for the identification of conserved and varying genomic sequences within lineages of a species. Investigating genetic variations might lead to the discovery of genes present in a subset of lineages, which might contribute into beneficial agronomic traits such as stress resistance or yield. The content of varying genomic regions in the pan-genome could include protein-coding genes as well as microRNA(miRNAs), small non-coding RNAs playing key roles in the regulation of gene expression. In this study, we performed in silico miRNA identification from the genomic sequences of 54 lineages of Brachypodium distachyon, aiming to explore varying miRNA contents and their functional interactions. A total of 115 miRNA families were identified in 54 lineages, 56 of which were found to be present in all lineages. The miRNA families were classified based on their conservation among lineages and potential mRNA targets were identified. Obtaining information about regulatory mechanisms stemming from these miRNAs offers strong potential to provide a better insight into the complex traits that were potentially present in some lineages. Future work could lead us to introduce these traits to different lineages or other economically important plant species in order to promote their survival in different environmental conditions.

Comparative Small RNA Profiling and Functional Exploration on Wheat With High- and Low-Cadmium Accumulation

Cadmium is a toxic metal widely found in workplaces and plant soil because of extensive industrialization. Wheat is an important source of food generated from plant soil. The different responses of wheat against different omic levels of cadmium have been observed and widely studied worldwide. With the development of high-throughput sequencing, micro-level biological research has extended to the microRNA level. In this study, high-cadmium-accumulating wheat cultivars (Annong9267) and low-cadmium-accumulating wheat cultivars (Qian 102032) were used as experimental models. The two cultivars were treated by Cd for 2 h to explore the microRNA profiles in root and leaf tissues through small RNA sequencing. Important small RNAs, such as tae-miR9663-5p and tae-miR6201, and potential small RNA-mediated mechanisms associated with cadmium accumulation were identified by summarizing specific microRNA profiling patterns and their respective target genes. At the wheat roots and leaves, differentially expressed small RNAs related to cadmium accumulation in different plant tissues (roots or leaves) were identified, and functional enrichment analyses on target genes of differentially expressed miRNAs in low- and high-cadmium-accumulating wheat cultivars in different plant tissues (roots or leaves) obtained some known mature miRNAs and new miRNAs. The identified miRNA will be regarded as a potential screening biomarker for low-cadmium-accumulating wheat.

Comparative profiling of roots small RNA expression and corresponding gene ontology and pathway analyses for low- and high-cadmium-accumulating genotypes of wheat in response to cadmium stress

MicroRNAs (miRNAs) participate in multiple biological processes in plant. Cd accumulation ability differs among varieties in wheat, but little is known about miRNAs and their function in Cd accumulation of wheat under Cd stress. Therefore, the present study detected small RNAs responsible for differential Cd accumulation between two contrasting wheat genotypes (low-Cd accumulation one L17 and high-Cd accumulation one H17) to identify novel targets to further study Cd stress in wheat. Under normal conditions, 139 miRNAs were differentially expressed between L17 and H17, while this value reached 142 after Cd exposure. For Cd-induced DEMs, total 25 miRNAs were differentially expressed in L17 after Cd treatment, while, 70 Cd-induced DEMs were found in H17. Moreover, GO analysis revealed that target genes of DEMs related to lipid biosynthetic process and chlorophyll binding are uniquely enriched in L17, target genes of DEMs related to ribosome biogenesis and sucrose alpha-glucosidase activity are uniquely enriched in H17. By pathway analysis, target genes of DEMs related to PI3K-Akt signaling pathway was specifically enriched in L17, target genes of DEMs related to carbohydrate digestion and absorption pathway was uniquely enriched in H17. In addition, miRNA-gene co-expression showed that tae-miR9774 was uniquely expressed between L17Cd and L17CK, while tae-miR398 was specially expressed between H17Cd and H17CK. Our results suggested that Cd-accumulating ability of L17 and H17 varied from the expression of induced unique miRNA, such as expression of tae-miR-9774 and tae-miR-398. Our study not only provide the foundation for further exploring the miRNAs-induced molecular mechanisms of Cd accumulation in wheat but also supply novel strategies for improving phytoremediation ability of food plants through genetic engineering.

Discovery and profiling of small RNAs from Puccinia triticina by deep sequencing and identification of their potential targets in wheat

Cross-kingdom RNAi is a well-documented phenomenon where sRNAs generated by host and pathogens may govern resistance or susceptible phenotypes during host-pathogen interaction. With the first example of the direct involvement of fungal generated sRNAs in virulence of plant pathogenic fungi Botrytis cinerea and recently from Puccinia striiformis f. sp. tritici, we attempted to identify sRNAs in Puccinia triticina (P. triticina). Four sRNA libraries were prepared and sequenced using Illumina sequencing technology and a total of ~ 1-1.28 million potential sRNAs and two microRNA-like small RNA (mil-RNAs) candidates were identified. Computational prediction of targets using a common set of sRNAs and P. triticina mil-RNAs (pt-mil-RNAs) within P. triticina and wheat revealed the majority of the targets as repetitive elements in P. triticina whereas in wheat, the target genes were identified to be involved in many biological processes including defense-related pathways. We found 9 receptor-like kinases (RLKs) and 14 target genes of each related to reactive oxygen species (ROS) pathway and transcription factors respectively, including significant numbers of target genes from various other categories. Expression analysis of twenty selected sRNAs, targeting host genes pertaining to ROS related, disease resistance, metabolic processes, transporter, apoptotic inhibitor, and transcription factors along with two pt-mil-RNAs by qRT-PCR showed distinct patterns of expression of the sRNAs in urediniospore-specific libraries. In this study, for the first time, we report identification of novel sRNAs identified in P. triticina including two pt-mil-RNAs that may play an important role in biotrophic growth and pathogenicity.

Chromosome-based survey sequencing reveals the genome organization of wild wheat progenitor Triticum dicoccoides

Wild emmer wheat (Triticum turgidum ssp. dicoccoides) is the progenitor of wheat. We performed chromosome-based survey sequencing of the 14 chromosomes, examining repetitive sequences, protein-coding genes, miRNA/target pairs and tRNA genes, as well as syntenic relationships with related grasses. We found considerable differences in the content and distribution of repetitive sequences between the A and B subgenomes. The gene contents of individual chromosomes varied widely, not necessarily correlating with chromosome size. We catalogued candidate agronomically important loci, along with new alleles and flanking sequences that can be used to design exome sequencing. Syntenic relationships and virtual gene orders revealed several small-scale evolutionary rearrangements, in addition to providing evidence for the 4AL-5AL-7BS translocation in wild emmer wheat. Chromosome-based sequence assemblies contained five novel miRNA families, among 59 families putatively encoded in the entire genome which provide insight into the domestication of wheat and an overview of the genome content and organization.

TaMIR1139: a wheat miRNA responsive to Pi-starvation, acts a critical mediator in modulating plant tolerance to Pi deprivation

Wheat miRNA member TaMIR1139 targets genes functional in various families and plays crucial roles in regulating plant Pi starvation tolerance. Through regulating target genes at posttranscriptional or translational level, plant miRNAs are involved in mediating diverse biological processes associated with growth, development, and responses to adverse stresses. In this study, we characterized the expression pattern and function of TaMIR1139, a miRNA member of wheat (T. aestivum) under Pi deprivation. TaMIR1139 precursor is also present in N. tabucum, suggesting the conserved nature of miR1139 across monocots and eudicots. TaMIR1139 targets seven genes within different families. The transcripts abundance of TaMIR1139 was induced upon Pi deprivation and the upregulated expression under Pi starvation was downregulated by the Pi recovery treatment, In contrast, the genes targeted by TaMIR1139 exhibited reduced transcripts upon Pi starvation and their downregulated expression was recovered by Pi-recovery condition, suggesting the regulation of them under TaMIR1139 through a cleavage mechanism. TaMIR1139 overexpression conferred the Pi-deprived plants improved phenotype, biomass, photosynthesis, and Pi acquisition. Transcriptome analysis identified numerous genes involving biological process, cellular components, and molecular function were differentially expressed in the TaMIR1139 overexpression lines, which suggests the TaMIR1139-mediated plant Pi starvation tolerance to be associated with the role of miRNA in extensively modulating the transcript profiling. A phosphate transporter (PT) gene NtPT showed significantly upregulated expression in TaMIR1139 overexpression lines; overexpression of it conferred plants improved Pi acquisition upon Pi starvation, suggesting its contribution to the TaMIR1139-mediated plant low-Pi stress resistance. Our investigation indicates that TaMIR1139 is critical in plant Pi starvation tolerance through transcriptionally regulating the target genes and modulating the Pi stress-defensiveness processes.

A large-scale multiomics analysis of wheat stem solidness and the wheat stem sawfly feeding response, and syntenic associations in barley, Brachypodium, and rice

The wheat stem sawfly (WSS), Cephus cinctus Norton (Hymenoptera: Cephidae), is an important pest of wheat and other cereals, threatening the quality and quantity of grain production. WSS larvae feed and develop inside the stem where they are protected from the external environment; therefore, pest management strategies primarily rely on host plant resistance. A major locus on the long arm of wheat chromosome 3B underlies most of the variation in stem solidness; however, the impact of stem solidness on WSS feeding has not been completely characterized. Here, we used a multiomics approach to examine the response to WSS in both solid- and semi-solid-stemmed wheat varieties. The combined transcriptomic, proteomic, and metabolomic data revealed that two important molecular pathways, phenylpropanoid and phosphate pentose, are involved in plant defense against WSS. We also detected a general downregulation of several key defense transcripts, including those encoding secondary metabolites such as DIMBOA, tricetin, and lignin, which suggested that the WSS larva might interfere with plant defense. We comparatively analyzed the stem solidness genomic region known to be associated with WSS tolerance in wild emmer, durum, and bread wheats, and described syntenic regions in the close relatives barley, Brachypodium, and rice. Additionally, microRNAs identified from the same genomic region revealed potential regulatory pathways associated with the WSS response. We propose a model outlining the molecular responses of the WSS–wheat interactions. These findings provide insight into the link between stem solidness and WSS feeding at the molecular level.

Redox and Ionic Homeostasis Regulations against Oxidative, Salinity and Drought Stress in Wheat (A Systems Biology Approach)

Systems biology and omics has provided a comprehensive understanding about the dynamics of the genome, metabolome, transcriptome, and proteome under stress. In wheat, abiotic stresses trigger specific networks of pathways involved in redox and ionic homeostasis as well as osmotic balance. These networks are considerably more complicated than those in model plants, and therefore, counter models are proposed by unifying the approaches of omics and stress systems biology. Furthermore, crosstalk among these pathways is monitored by the regulation and streaming of transcripts and genes. In this review, we discuss systems biology and omics as a promising tool to study responses to oxidative, salinity, and drought stress in wheat.

Assembly and Annotation of Transcriptome Provided Evidence of miRNA Mobility between Wheat and Wheat Stem Sawfly

Wheat Stem Sawfly (WSS), Cephus Cinctus Norton (Hymenoptera: Cephidae), is one of the most important pests, causing yield and economic losses in wheat and barley. The lack of information about molecular mechanisms of WSS for defeating plant's resistance prevents application of effective pest control strategies therefore, it is essential to identify the genes and their regulators behind WSS infestations. Long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) are recognized with their regulatory functions on gene expression, tuning protein production by controlling transcriptional and post-transcriptional activities. A transcriptome-guided approach was followed in order to identify miRNAs, lncRNAs, and mRNA of WSS, and their interaction networks. A total of 1,893 were presented here as differentially expressed between larva and adult WSS insects. There were 11 miRNA families detected in WSS transcriptome. Together with the annotation of 1,251 novel mRNAs, the amount of genetic information available for WSS was expanded. The network between WSS miRNAs, lncRNAs, and mRNAs suggested miRNA-mediated regulatory roles of lncRNAs as competing endogenous RNAs. In the light of the previous evidence that small RNA molecules of a pathogen could suppress the immune response of host plant, we analyzed the putative interactions between larvae and wheat at the miRNA level. Overall, this study provides a profile of larva and adult WSS life stages in terms of coding and non-coding elements. These findings also emphasize the potential roles of wheat and larval miRNAs in wheat resistance to infestation and in the suppression of resistance which is critical for the development of effective pest control strategies.

RNA Sequencing and Co-expressed Long Non-coding RNA in Modern and Wild Wheats

There is an urgent need for the improvement of drought-tolerant bread and durum wheat. The huge and complex genome of bread wheat (BBAADD genome) stands as a vital obstruction for understanding the molecular mechanism underlying drought tolerance. However, tetraploid wheat (Triticum turgidum ssp., BBAA genome) is an ancestor of modern bread wheat and offers an important model for studying the drought response due to its less complex genome. Additionally, several wild relatives of tetraploid wheat have already shown a significant drought tolerance. We sequenced root transcriptome of three tetraploid wheat varieties with varying stress tolerance profiles, and built differential expression library of their transcripts under control and drought conditions. More than 5,000 differentially expressed transcripts were identified from each genotype. Functional characterization of transcripts specific to drought-tolerant genotype, revealed their association with osmolytes production and secondary metabolite pathways. Comparative analysis of differentially expressed genes and their non-coding RNA partners, long noncoding RNAs and microRNAs, provided valuable insight to gene expression regulation in response to drought stress. LncRNAs as well as coding transcripts share similar structural features in different tetraploid species; yet, lncRNAs slightly differ from coding transcripts. Several miRNA-lncRNA target pairs were detected as differentially expressed in drought stress. Overall, this study suggested an important pool of transcripts where their manipulations confer a better performance of wheat varieties under drought stress.

Abiotic stress miRNomes in the Triticeae

The continued growth in world population necessitates increases in both the quantity and quality of agricultural production. Triticeae members, particularly wheat and barley, make an important contribution to world food reserves by providing rich sources of carbohydrate and protein. These crops are grown over diverse production environments that are characterized by a range of environmental or abiotic stresses. Abiotic stresses such as drought, heat, salinity, or nutrient deficiencies and toxicities cause large yield losses resulting in economic and environmental damage. The negative effects of abiotic stresses have increased at an alarming rate in recent years and are predicted to further deteriorate due to climate change, land degradation, and declining water supply. New technologies have provided an important tool with great potential for improving crop tolerance to the abiotic stresses: microRNAs (miRNAs). miRNAs are small regulators of gene expression that act on many different molecular and biochemical processes such as development, environmental adaptation, and stress tolerance. miRNAs can act at both the transcriptional and post-transcriptional levels, although post-transcriptional regulation is the most common in plants where miRNAs can inhibit the translation of their mRNA targets via complementary binding and cleavage. To date, expression of several miRNA families such as miR156, miR159, and miR398 has been detected as responsive to environmental conditions to regulate stress-associated molecular mechanisms individually and/or together with their various miRNA partners. Manipulation of these miRNAs and their targets may pave the way to improve crop performance under several abiotic stresses. Here, we summarize the current status of our knowledge on abiotic stress-associated miRNAs in members of the Triticeae tribe, specifically in wheat and barley, and the miRNA-based regulatory mechanisms triggered by stress conditions. Exploration of further miRNA families together with their functions under stress will improve our knowledge and provide opportunities to enhance plant performance to help us meet global food demand.

Identification and characterization of durum wheat microRNAs in leaf and root tissues

MicroRNAs are a class of post-transcriptional regulators of plant developmental and physiological processes and responses to environmental stresses. Here, we present the study regarding the annotation and characterization of MIR genes conducted in durum wheat. We characterized the miRNAome of leaf and root tissues at tillering stage under two environmental conditions: irrigated with 100% (control) and 55% of evapotranspiration (early water stress). In total, 90 microRNAs were identified, of which 32 were classified as putative novel and species-specific miRNAs. In addition, seven microRNA homeologous groups were identified in each of the two genomes of the tetraploid durum wheat. Differential expression analysis highlighted a total of 45 microRNAs significantly differentially regulated in the pairwise comparisons leaf versus root. The miRNA families, miR530, miR395, miR393, miR5168, miR396 and miR166, miR171, miR319, and miR167, were the most expressed in leaves in comparison to roots. Putative microRNA targets were predicted for both five and three prime sequences derived from the stem-loop of the MIR gene. Gene ontology analysis showed significant overrepresented gene categories in microRNA targets belonging to transcription factors, phenylpropanoids, oxydases, and lipid binding-protein. This work represents one of the first genome wide characterization of MIR genes in durum wheat, identifying leaf and root tissue-specific microRNAs. This genomic identification of microRNAs together with the analysis of their expression profiles is a well-accepted starting point leading to a better comprehension of the role of MIR genes in the genus Triticum.

Uncovering leaf rust responsive miRNAs in wheat (Triticum aestivum L.) using high-throughput sequencing and prediction of their targets through degradome analysis

Deep sequencing identified 497 conserved and 559 novel miRNAs in wheat, while degradome analysis revealed 701 targets genes. QRT-PCR demonstrated differential expression of miRNAs during stages of leaf rust progression. Bread wheat (Triticum aestivum L.) is an important cereal food crop feeding 30 % of the world population. Major threat to wheat production is the rust epidemics. This study was targeted towards identification and functional characterizations of micro(mi)RNAs and their target genes in wheat in response to leaf rust ingression. High-throughput sequencing was used for transcriptome-wide identification of miRNAs and their expression profiling in retort to leaf rust using mock and pathogen-inoculated resistant and susceptible near-isogenic wheat plants. A total of 1056 mature miRNAs were identified, of which 497 miRNAs were conserved and 559 miRNAs were novel. The pathogen-inoculated resistant plants manifested more miRNAs compared with the pathogen infected susceptible plants. The miRNA counts increased in susceptible isoline due to leaf rust, conversely, the counts decreased in the resistant isoline in response to pathogenesis illustrating precise spatial tuning of miRNAs during compatible and incompatible interaction. Stem-loop quantitative real-time PCR was used to profile 10 highly differentially expressed miRNAs obtained from high-throughput sequencing data. The spatio-temporal profiling validated the differential expression of miRNAs between the isolines as well as in retort to pathogen infection. Degradome analysis provided 701 predicted target genes associated with defense response, signal transduction, development, metabolism, and transcriptional regulation. The obtained results indicate that wheat isolines employ diverse arrays of miRNAs that modulate their target genes during compatible and incompatible interaction. Our findings contribute to increase knowledge on roles of microRNA in wheat-leaf rust interactions and could help in rust resistance breeding programs.

Deep sequencing of wheat sRNA transcriptome reveals distinct temporal expression pattern of miRNAs in response to heat, light and UV

Understanding of plant adaptation to abiotic stresses has implications in plant breeding, especially in the context of climate change. MicroRNAs (miRNAs) and short interfering RNAs play a crucial role in gene regulation. Here, wheat plants were exposed to one of the following stresses: continuous light, heat or ultraviolet radiations over five consecutive days and leaf tissues from three biological replicates were harvested at 0, 1, 2, 3, 7 and 10 days after treatment (DAT). A total of 72 small RNA libraries were sequenced on the Illumina platform generating ~524 million reads corresponding to ~129 million distinct tags from which 232 conserved miRNAs were identified. The expression levels of 1, 2 and 79 miRNAs were affected by ultraviolet radiation, continuous light and heat, respectively. Approximately 55% of the differentially expressed miRNAs were downregulated at 0 and 1 DAT including miR398, miR528 and miR156 that control mRNAs involved in activation of signal transduction pathways and flowering. Other putative targets included histone variants and methyltransferases. These results suggest a temporal miRNA-guided post-transcriptional regulation that enables wheat to respond to abiotic stresses, particularly heat. Designing novel wheat breeding strategies such as regulatory gene-based marker assisted selection depends on accurate identification of stress induced miRNAs.

PeTMbase: A Database of Plant Endogenous Target Mimics (eTMs)

MicroRNAs (miRNA) are small endogenous RNA molecules, which regulate target gene expression at post-transcriptional level. Besides, miRNA activity can be controlled by a newly discovered regulatory mechanism called endogenous target mimicry (eTM). In target mimicry, eTMs bind to the corresponding miRNAs to block the binding of specific transcript leading to increase mRNA expression. Thus, miRNA-eTM-target-mRNA regulation modules involving a wide range of biological processes; an increasing need for a comprehensive eTM database arose. Except miRSponge with limited number of Arabidopsis eTM data no available database and/or repository was developed and released for plant eTMs yet. Here, we present an online plant eTM database, called PeTMbase (http://petmbase.org), with a highly efficient search tool. To establish the repository a number of identified eTMs was obtained utilizing from high-throughput RNA-sequencing data of 11 plant species. Each transcriptome libraries is first mapped to corresponding plant genome, then long non-coding RNA (lncRNA) transcripts are characterized. Furthermore, additional lncRNAs retrieved from GREENC and PNRD were incorporated into the lncRNA catalog. Then, utilizing the lncRNA and miRNA sources a total of 2,728 eTMs were successfully predicted. Our regularly updated database, PeTMbase, provides high quality information regarding miRNA:eTM modules and will aid functional genomics studies particularly, on miRNA regulatory networks.

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.

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.

Involvement of miR528 in the Regulation of Arsenite Tolerance in Rice (Oryza sativa L.)

Tens of miRNAs were previously established as being arsenic (As) stress responsive in rice. However, their functional role in As tolerance remains unclear. This study demonstrates that transgenic plants overexpressing miR528 (Ubi::MIR528) were more sensitive to arsenite [As(III)] compared with wild-type (WT) rice. Under normal and stress conditions, miR528-5p and -3p were highly up-regulated in both the roots and leaves of transgenic plants, which exhibited a negative correlation with the expression of seven target genes. Compared with WT plants, Ubi::MIR528 plants showed excessive oxidative stress generation and remarkable amino acid content changes in the roots and leaves upon As(III) exposure. Notably, the expression profiles of diverse functional genes were clearly different between WT and transgenic plants. Thus, the observed As(III) sensitivity of Ubi::MIR528 plants was likely due to the strong alteration of antioxidant enzyme activity and amino acid profiles and the impairment of the As(III) uptake, translocation, and tolerance systems of rice.

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.

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.

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.

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.