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

Transcriptome variations in hybrids of wild emmer wheat (Triticum turgidum ssp. dicoccoides)

BACKGROUND

Wild emmer wheat is a great candidate to revitalize domesticated wheat genetic diversity. Recent years have seen intensive investigation into the evolution and domestication of wild emmer wheat, including whole-genome DNA and transcriptome sequencing. However, the impact of intraspecific hybridization on the transcriptome of wild emmer wheat has been poorly studied. In this study, we assessed changes in methylation patterns and transcriptomic variations in two accessions of wild emmer wheat collected from two marginal populations, Mt. Hermon and Mt. Amasa, and in their stable F4 hybrid.

RESULTS

Methylation-Sensitive Amplified Polymorphism (MSAP) detected significant cytosine demethylation in F4 hybrids vs. parental lines, suggesting potential transcriptome variation. After a detailed analysis, we examined nine RNA-Seq samples, which included three biological replicates from the F4 hybrid and its parental lines. RNA-Seq databases contained approximately 200 million reads, with each library consisting of 15 to 25 million reads. There are a total of 62,490 well-annotated genes in these databases, with 6,602 genes showing differential expression between F4 hybrid and parental lines Mt. Hermon and Mt. Amasa. The differentially expressed genes were classified into four main categories based on their expression patterns. Gene ontology (GO) analysis revealed that differentially expressed genes are associated with DNA/RNA metabolism, photosynthesis, stress response, phosphorylation and developmental processes.

CONCLUSION

This study highlights the significant transcriptomic changes resulting from intraspecific hybridization within natural plant populations, which might aid the nascent hybrid in adapting to various environmental conditions.

Identification of drought-tolerant hub genes in Iranian KC-2226 genotype of Aegilops tauschii using transcriptomic analysis

Aegilops tauschii, as a donor of D genome to the bread wheat with a valuable source of resistance to different biotic and abiotic stresses, is used to improve the quality of wheat cultivars. Every genotype has a specific genetic content, the investigation of which can lead to the identification of useful genes such as stress tolerance genes, including drought. Therefore, 23 genotypes of Ae. tauschii were selected to evaluate their morphological and physiological traits under greenhouse conditions. Among them, a superior tolerant genotype (KC-2226) was chosen for transcriptomic analysis. Our result showed that 5007 and 3489 genes were deferentially up- and downregulated, respectively. Upregulated genes were involved in photosynthesis, glycolysis/gluconeogenesis, and amino acid biosynthesis whereas downregulated genes were often engaged in DNA synthesis, replication, repair and topological changes. The result of protein-protein interaction network analysis showed that AT1G76550 (1.46), AT1G20950 (1.42), IAR4 (1.19), and PYD2 (1.16) among upregulated genes and THY-1 (44), PCNA1 (41) and TOPII (22) among down-regulated genes had the highest interactions with other genes. In conclusion, Ae. tauschii employs elevated transcription of specific genes involved in photosynthesis, glycolysis and gluconeogenesis and amino acid biosynthesis pathways rather than genes active in DNA synthesis and repair to provide the energy needed for the plant to survive under stress conditions.

Aegilops tauschii genome assembly Aet v5.0 features greater sequence contiguity and improved annotation

Aegilops tauschii is the donor of the D subgenome of hexaploid wheat and an important genetic resource. The reference-quality genome sequence Aet v4.0 for Ae. tauschii acc. AL8/78 was therefore an important milestone for wheat biology and breeding. Further advances in sequencing acc. AL8/78 and release of the Aet v5.0 sequence assembly are reported here. Two new optical maps were constructed and used in the revision of pseudomolecules. Gaps were closed with Pacific Biosciences long-read contigs, decreasing the gap number by 38,899. Transposable elements and protein-coding genes were reannotated. The number of annotated high-confidence genes was reduced from 39,635 in Aet v4.0 to 32,885 in Aet v5.0. A total of 2245 biologically important genes, including those affecting plant phenology, grain quality, and tolerance of abiotic stresses in wheat, was manually annotated and disease-resistance genes were annotated by a dedicated pipeline. Disease-resistance genes encoding nucleotide-binding site domains, receptor-like protein kinases, and receptor-like proteins were preferentially located in distal chromosome regions, whereas those encoding transmembrane coiled-coil proteins were dispersed more evenly along the chromosomes. Discovery, annotation, and expression analyses of microRNA (miRNA) precursors, mature miRNAs, and phasiRNAs are reported, including miRNA target genes. Other small RNAs, such as hc-siRNAs and tRFs, were characterized. These advances enhance the utility of the Ae. tauschii genome sequence for wheat genetics, biotechnology, and breeding.

Recognizing the hidden half in wheat: root system attributes associated with drought tolerance

Improving drought tolerance in wheat is crucial for maintaining productivity and food security. Roots are responsible for the uptake of water from soil, and a number of root traits are associated with drought tolerance. Studies have revealed many quantitative trait loci and genes controlling root development in plants. However, the genetic dissection of root traits in response to drought in wheat is still unclear. Here, we review crop root traits associated with drought, key genes governing root development in plants, and quantitative trait loci and genes regulating root system architecture under water-limited conditions in wheat. Deep roots, optimal root length density and xylem diameter, and increased root surface area are traits contributing to drought tolerance. In view of the diverse environments in which wheat is grown, the balance among root and shoot traits, as well as individual and population performance, are discussed. The known functions of key genes provide information for the genetic dissection of root development of wheat in a wide range of conditions, and will be beneficial for molecular marker development, marker-assisted selection, and genetic improvement in breeding for drought tolerance.

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.

Genome evolution during bread wheat formation unveiled by the distribution dynamics of SSR sequences on chromosomes using FISH

BACKGROUND

During the bread wheat speciation by polyploidization, a series of genome rearrangement and sequence recombination occurred. Simple sequence repeat (SSR) sequences, predominately located in heterochromatic regions of chromosomes, are the effective marker for tracing the genomic DNA sequence variations. However, to date the distribution dynamics of SSRs on chromosomes of bread wheat and its donors, including diploid and tetraploid Triticum urartu, Aegilops speltoides, Aegilops tauschii, Triticum turgidum ssp. dicocoides, reflecting the genome evolution events during bread wheat formation had not been comprehensively investigated.

RESULTS

The genome evolution was studied by comprehensively comparing the distribution patterns of (AAC)_n, (AAG)_n, (AGC)_n and (AG)_n in bread wheat Triticum aestivum var. Chinese Spring and its progenitors T. urartu, A. speltoides, Ae. tauschii, wild tetroploid emmer wheat T. dicocoides, and cultivated emmer wheat T. dicoccum. Results indicated that there are specific distribution patterns in different chromosomes from different species for each SSRs. They provided efficient visible markers for identification of some individual chromosomes and SSR sequence evolution tracing from the diploid progenitors to hexaploid wheat. During wheat speciation, the SSR sequence expansion occurred predominately in the centromeric and pericentromeric regions of B genome chromosomes accompanied by little expansion and elimination on other chromosomes. This result indicated that the B genome might be more sensitive to the "genome shock" and more changeable during wheat polyplodization.

CONCLUSIONS

During the bread wheat evolution, SSRs including (AAC)_n, (AAG)_n, (AGC)_n and (AG)_n in B genome displayed the greatest changes (sequence expansion) especially in centromeric and pericentromeric regions during the polyploidization from Ae. speltoides S genome, the most likely donor of B genome. This work would enable a better understanding of the wheat genome formation and evolution and reinforce the viewpoint that B genome was originated from S genome.

RNAseq analysis reveals drought-responsive molecular pathways with candidate genes and putative molecular markers in root tissue of wheat

Drought is one of the major impediments in wheat productivity. Traditional breeding and marker assisted QTL introgression had limited success. Available wheat genomic and RNA-seq data can decipher novel drought tolerance mechanisms with putative candidate gene and marker discovery. Drought is first sensed by root tissue but limited information is available about how roots respond to drought stress. In this view, two contrasting genotypes, namely, NI5439 41 (drought tolerant) and WL711 (drought susceptible) were used to generate ~78.2 GB data for the responses of wheat roots to drought. A total of 45139 DEGs, 13820 TF, 288 miRNAs, 640 pathways and 435829 putative markers were obtained. Study reveals use of such data in QTL to QTN refinement by analysis on two model drought-responsive QTLs on chromosome 3B in wheat roots possessing 18 differentially regulated genes with 190 sequence variants (173 SNPs and 17 InDels). Gene regulatory networks showed 69 hub-genes integrating ABA dependent and independent pathways controlling sensing of drought, root growth, uptake regulation, purine metabolism, thiamine metabolism and antibiotics pathways, stomatal closure and senescence. Eleven SSR markers were validated in a panel of 18 diverse wheat varieties. For effective future use of findings, web genomic resources were developed. We report RNA-Seq approach on wheat roots describing the drought response mechanisms under field drought conditions along with genomic resources, warranted in endeavour of wheat productivity.

In silico determination of transposon-derived miRNAs and targets in Aegilops species

Transposable elements (TEs) are found almost in all living organism, shaping organisms' genomes. miRNAs are noncoding RNA types which are especially important in gene expression regulations. Many previously determined plant miRNAs are identical/homologous to transposons (TE-MIR). The aim of this study was computational characterization of novel TE-related miRNAs and their targets in Aegilops genome by using stringent criteria. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed by BLAST2GO. Seventeen novel TE-related miRNAs in Aegilops genome were identified for the first time. GO analyses indicated that 40 targets played different roles in biological processes, cellular components and molecular functions. Moreover, these genes were involved in 10 metabolic pathways such as purine metabolism, nitrogen metabolism, oxidative phosphorylation, etc. as a result of KEGG analyses. Identification of miRNAs and their targets are significant to understand miRNA-TEs relationships and even how TEs affect plant growth and development. Obtaining results of this study are expected to provide possible new insight into Aegilops and its related species, wheat, with respect to miRNAs evolution and domestication.Communicated by Ramaswamy H. Sarma.

MicroRNA-guided regulation of heat stress response in wheat

Background

With rising global temperature, understanding plants’ adaptation to heat stress has implications in plant breeding. MicroRNAs (miRNAs) are small, non-coding, regulatory RNAs guiding gene expression at the post-transcriptional level. In this study, small RNAs and the degradome (parallel analysis of RNA ends) of leaf tissues collected from control and heat-stressed wheat plants immediately at the end of the stress period and 1 and 4 days later were analysed.

Results

Sequencing of 24 small RNA libraries produced 55.2 M reads while 404 M reads were obtained from the corresponding 24 PARE libraries. From these, 202 miRNAs were ascertained, of which mature miRNA evidence was obtained for 104 and 36 were found to be differentially expressed after heat stress. The PARE analysis identified 589 transcripts targeted by 84 of the ascertained miRNAs. PARE sequencing validated the targets of the conserved members of miRNA156, miR166 and miR393 families as squamosa promoter-binding-like, homeobox leucine-zipper and transport inhibitor responsive proteins, respectively. Heat stress responsive miRNA targeted superoxide dismutases and an array of homeobox leucine-zipper proteins, F-box proteins and protein kinases. Query of miRNA targets to interactome databases revealed a predominant association of stress responses such as signalling, antioxidant activity and ubiquitination to superoxide dismutases, F-box proteins, pentatricopeptide repeat-containing proteins and mitochondrial transcription termination factor-like proteins.

Conclusion

The interlaced data set generated in this study identified and validated heat stress regulated miRNAs and their target genes associated with thermotolerance. Such accurate identification and validation of miRNAs and their target genes are essential to develop novel regulatory gene-based breeding strategies.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5799-6) contains supplementary material, which is available to authorized users.

Identification and Characterization of Salt-Responsive MicroRNAs in Vicia faba by High-Throughput Sequencing

Salt stress has detrimental effects on plant growth and development. MicroRNAs (miRNAs) are a class of noncoding RNAs that are involved in post-transcriptional gene expression regulation. In this study, small RNA sequencing was employed to identify the salt stress-responsive miRNAs of the salt-sensitive Hassawi-3 and the salt-tolerant ILB4347 genotypes of faba bean, growing under salt stress. A total of 527 miRNAs in Hassawi-3 plants, and 693 miRNAs in ILB4347 plants, were found to be differentially expressed. Additionally, 284 upregulated and 243 downregulated miRNAs in Hassawi-3, and 298 upregulated and 395 downregulated miRNAs in ILB4347 plants growing in control and stress conditions were recorded. Target prediction and annotation revealed that these miRNAs regulate specific salt-responsive genes, which primarily included genes encoding transcription factors and laccases, superoxide dismutase, plantacyanin, and F-box proteins. The salt-responsive miRNAs and their targets were functionally enriched by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, which showed that the miRNAs were involved in salt stress-related biological pathways, including the ABC transporter pathway, MAPK signaling pathway, plant hormone signal transduction, and the phosphatidylinositol signaling system, among others, suggesting that the miRNAs play an important role in the salt stress tolerance of the ILB4347 genotype. These results offer a novel understanding of the regulatory role of miRNAs in the salt response of the salt-tolerant ILB4347 and the salt-sensitive Hassawi-3 faba bean genotypes.

Salinity-associated microRNAs and their potential roles in mediating salt tolerance in rice colonized by the endophytic root fungus Piriformospora indica

Piriformospora indica (P. indica), an endophytic root fungus, supports the growth and enhanced tolerance of plants to biotic and abiotic stresses. Several recent studies showed the significant role of small RNA (sRNA) molecules including microRNAs (miRNAs) in plant adaption to environmental stress, but little is known concerning the symbiosis-mediated salt stress tolerance regulated at miRNAs level. The overarching goal of this research is to elucidate the impact of miRNAs in regulating the P. indica-mediated salt tolerance in rice. Applying sRNA-seq analysis led to identify a set of 547 differentially abundant miRNAs in response to P. indica inoculation and salt stress. These included 206 rice-specific and 341 previously known miRNAs from other plant species. In silico analysis of miRNAs predictions of the differentially abundant miRNAs led to identifying of 193 putatively target genes, most of which were encoded either genes or transcription factors involved in nutrient uptake, sodium ion transporters, growth regulators, and auxin- responsive proteins. The rice-specific miRNAs targeted the transcription factors involved in the import of potassium ions into the root cells, the export of sodium ions, and plant growth and development. Interestingly, P. indica affected the differential abundance of miRNAs regulated genes and transcription factors linked to salt stress tolerance. Our data helps to understand the molecular basis of salt stress tolerance mediated by symbionts in plant and the potential impact of miRNAs for genetic improvement of rice varieties for tolerance to salt stress.

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.

Identification of conserved miRNA molecules in einkorn wheat (Triticum monococcum subsp. monococcum) by using small RNA sequencing analysis

Triticum monococcum subsp. monococcum as a first cultivated diploid wheat species possesses desirable agronomic and quality characteristics. Drought and salinity are the most dramatic environmental stress factors that have serious impact on yield and quality of crops; however, plants can use alternative defense mechanisms against these stresses. The posttranscriptional alteration of gene expression by microRNAs (miRNAs) is one of the most conserved mechanisms. In plant species including wheat genomes, miRNAs have been implicated in the management of salt and drought stress; however, studies on einkorn wheat (Triticum monococcum subsp. monococcum) are not yet available. In this study, we aimed to identify conserved miRNAs in einkorn wheat using next generation sequencing technology and bioinformatics analysis. In order to include a larger set of miRNAs, small RNA molecules from pooled plant samples grown under normal, drought, and salinity conditions were used for the library preparation and sequence analysis. After bioinformatics analysis, we identified 167 putative mature miRNA sequences belonging to 140 distinct miRNA families. We also presented a comparative analysis to propose that miRNAs and their target genes were involved in salt and drought stress control in addition to a comprehensive analysis of the scanned target genes in the T. aestivum genome.

High-throughput sequencing of small RNAs revealed the diversified cold-responsive pathways during cold stress in the wild banana (Musa itinerans)

BACKGROUND

Cold stress is one of the most severe abiotic stresses affecting the banana production. Although some miRNAs have been identified, little is known about the role of miRNAs in response to cold stress in banana, and up to date, there is no report about the role of miRNAs in the response to cold stress in the plants of the cultivated or wild bananas.

RESULT

Here, a cold-resistant line wild banana (Musa itinerans) from China was used to profile the cold-responsive miRNAs by RNA-seq during cold stress. Totally, 265 known mature miRNAs and 41 novel miRNAs were obtained. Cluster analysis of differentially expressed (DE) miRNAs indicated that some miRNAs were specific for chilling or 0 °C treated responses, and most of them were reported to be cold-responsive; however, some were seldom reported to be cold-responsive in response to cold stress, e.g., miR395, miR408, miR172, suggesting that they maybe play key roles in response to cold stress. The GO and KEGG pathway enrichment analysis of DE miRNAs targets indicated that there existed diversified cold-responsive pathways, and miR172 was found likely to play a central coordinating role in response to cold stress, especially in the regulation of CK2 and the circadian rhythm. Finally, qPCR assays indicated the related targets were negatively regulated by the tested DE miRNAs during cold stress in the wild banana.

CONCLUSIONS

In this study, the profiling of miRNAs by RNA-seq in response to cold stress in the plants of the wild banana (Musa itinerans) was reported for the first time. The results showed that there existed diversified cold-responsive pathways, which provided insight into the roles of miRNAs during cold stress, and would be helpful for alleviating cold stress and cold-resistant breeding in bananas.

Hotspots in the genomic architecture of field drought responses in wheat as breeding targets

Wheat can adapt to most agricultural conditions across temperate regions. This success is the result of phenotypic plasticity conferred by a large and complex genome composed of three homoeologous genomes (A, B, and D). Although drought is a major cause of yield and quality loss in wheat, the adaptive mechanisms and gene networks underlying drought responses in the field remain largely unknown. Here, we addressed this by utilizing an interdisciplinary approach involving field water status phenotyping, sampling, and gene expression analyses. Overall, changes at the transcriptional level were reflected in plant spectral traits amenable to field-level physiological measurements, although changes in photosynthesis-related pathways were found likely to be under more complex post-transcriptional control. Examining homoeologous genes with a 1:1:1 relationship across the A, B, and D genomes (triads), we revealed a complex genomic architecture for drought responses under field conditions, involving gene homoeolog specialization, multiple gene clusters, gene families, miRNAs, and transcription factors coordinating these responses. Our results provide a new focus for genomics-assisted breeding of drought-tolerant wheat cultivars.

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.

Transcriptional regulation of osmotic stress tolerance in wheat (Triticum aestivum L.)

The current review provides an updated, new insights into the regulation of transcription mediated underlying mechanisms of wheat plants to osmotic stress perturbations. Osmotic stress tolerance mechanisms being complex are governed by multiple factors at physiological, biochemical and at the molecular level, hence approaches like "OMICS" that can underpin mechanisms behind osmotic tolerance in wheat is of paramount importance. The transcription factors (TFs) are a class of molecular proteins, which are involved in regulation, modulation and orchestrating the responses of plants to a variety of environmental stresses. Recent reports have provided novel insights on the role of TFs in osmotic stress tolerance via direct molecular links. However, our knowledge on the regulatory role TFs during osmotic stress tolerance in wheat remains limited. The present review in its first part sheds light on the importance of studying the role of osmotic stress tolerance in wheat plants and second aims to decipher molecular mechanisms of TFs belonging to several classes, including DREB, NAC, MYB, WRKY and bHLH, which have been reported to engage in osmotic stress mediated gene expression in wheat and third part covers the systems biology approaches to understand the transcriptional regulation of osmotic stress and the role of long non-coding RNAs in response to osmotic stress with special emphasis on wheat. The current concept may lead to an understanding in molecular regulation and signalling interaction of TFs under osmotic stress to clarify challenges and problems for devising potential strategies to improve complex regulatory events involved in plant tolerance to osmotic stress adaptive pathways in wheat.

Transcriptomic analysis of Aegilops tauschii during long-term salinity stress

Aegilops tauschii is the diploid progenitor of the bread wheat D-genome. It originated from Iran and is a source of abiotic stress tolerance genes. However, little is known about the molecular events of salinity tolerance in Ae. tauschii. This study investigates the leaf transcriptional changes associated with long-term salt stress. Total RNA extracted from leaf tissues of control and salt-treated samples was sequenced using the Illumina technology, and more than 98 million high-quality reads were assembled into 255,446 unigenes with an average length of 1398 bp and an N50 of 2269 bp. Functional annotation of the unigenes showed that 93,742 (36.69%) had at least a significant BLAST hit in the SwissProt database, while 174,079 (68.14%) showed significant similarity to proteins in the NCBI nr database. Differential expression analysis identified 4506 salt stress-responsive unigenes. Bioinformatic analysis of the differentially expressed unigenes (DEUs) revealed a number of biological processes and pathways involved in the establishment of ion homeostasis, signaling processes, carbohydrate metabolism, and post-translational modifications. Fine regulation of starch and sucrose content may be important features involved in salt tolerance in Ae. tauschii. Moreover, 82% of DEUs mapped to the D-subgenome, including known QTL for salt tolerance, and these DEUs showed similar salt stress responses in other accessions of Ae. tauschii. These results could provide fundamental insight into the regulatory process underlying salt tolerance in Ae. tauschii and wheat and facilitate identification of genes involved in their salt tolerance mechanisms.

Exploration of miRNAs and target genes of cytoplasmic male sterility line in cotton during flower bud development

Cytoplasmic male sterility (CMS) lines provide crucial material to harness heterosis for crop plants, which serves as an important strategy for hybrid seed production. However, the molecular mechanism remains obscure. Although microRNAs (miRNAs) play important roles in vegetative growth and reproductive growth, there are few reports on miRNAs regulating the development of male sterility in Upland cotton. In present study, 12 small RNA libraries were constructed and sequenced for two development stages of flower buds from a CMS line and its maintainer line. Based on the results, 256 novel miRNAs were allocated to 141 new miRNA families, and 77 known miRNAs belonging to 54 conserved miRNA families were identified as well. Comparative analysis revealed that 61 novel and 10 conserved miRNAs were differentially expressed. Further transcriptome analysis identified 232 target genes for these miRNAs, which participated in cellular developmental process, cell death, pollen germination, and sexual reproduction. In addition, expression patterns of typical miRNA and the negatively regulated target genes, such as PPR, ARF, AP2, and AFB, were verified by qRT-PCR in cotton flower buds. These targets were previously reported to be related to reproduction development and male sterility, suggesting that miRNAs might act as regulators of CMS occurrence. Some miRNAs displayed specific expression profiles in special developmental stages of CMS line and its fertile hybrid (F₁). Present study offers new information on miRNAs and their related target genes in exploiting CMS mechanism, and revealing the miRNA regulatory networks in Upland cotton.

Deep sequencing and in silico analysis of small RNA library reveals novel miRNA from leaf Persicaria minor transcriptome

In current era, majority of microRNA (miRNA) are being discovered through computational approaches which are more confined towards model plants. Here, for the first time, we have described the identification and characterization of novel miRNA in a non-model plant, Persicaria minor (P. minor) using computational approach. Unannotated sequences from deep sequencing were analyzed based on previous well-established parameters. Around 24 putative novel miRNAs were identified from 6,417,780 reads of the unannotated sequence which represented 11 unique putative miRNA sequences. PsRobot target prediction tool was deployed to identify the target transcripts of putative novel miRNAs. Most of the predicted target transcripts (mRNAs) were known to be involved in plant development and stress responses. Gene ontology showed that majority of the putative novel miRNA targets involved in cellular component (69.07%), followed by molecular function (30.08%) and biological process (0.85%). Out of 11 unique putative miRNAs, 7 miRNAs were validated through semi-quantitative PCR. These novel miRNAs discoveries in P. minor may develop and update the current public miRNA database.

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.

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.

Genome-Wide Identification and Characterization of Salinity Stress-Responsive miRNAs in Wild Emmer Wheat (Triticum turgidum ssp. dicoccoides)

MicroRNAs (miRNAs) are a class of endogenous small noncoding RNAs which regulate diverse molecular and biochemical processes at a post-transcriptional level in plants. As the ancestor of domesticated wheat, wild emmer wheat (Triticum turgidum ssp. dicoccoides) has great genetic potential for wheat improvement. However, little is known about miRNAs and their functions on salinity stress in wild emmer. To obtain more information on miRNAs in wild emmer, we systematically investigated and characterized the salinity-responsive miRNAs using deep sequencing technology. A total of 88 conserved and 124 novel miRNAs were identified, of which 50 were proven to be salinity-responsive miRNAs, with 32 significantly up-regulated and 18 down-regulated. miR172b and miR1120a, as well as mi393a, were the most significantly differently expressed. Targets of these miRNAs were computationally predicted, then Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis found that the targets of salinity-responsive miRNAs were enriched in transcription factors and stress-related proteins. Finally, we investigated the expression profiles of seven miRNAs ranging between salt-tolerant and sensitive genotypes, and found that they played critical roles in salinity tolerance in wild emmer. Our results systematically identified the salinity-responsive miRNAs in wild emmer, not only enriching the miRNA resource but also laying the foundation for further study on the biological functions and evolution of miRNAs in wild wheat and beyond.

MicroRNAs in model and complex organisms

Non-coding RNAs such as microRNAs (miRNAs) are very tiny ribonucleotides having an essential role in gene regulation at both post-transcriptional and translational levels. They are very conserved and expressed in worms, flies, plants, and mammals in a sequence-specific manner. Furthermore, it is now possible to clone miRNAs using the new genome editing tool CRISPR/cas9, which shows benefit in control of untargeted effect. In this special issue, we tried to cover researches associated with functional roles of miRNAs accross model and complex organisms.

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.

Characterization of miR061 and its target genes in grapevine responding to exogenous gibberellic acid

MicroRNAs (miRNAs), as an important growth regulator, are also involved in gibberellic acid (GA) signaling, revealing much relationship between miRNAs and GA in various plant responses. Grape is highly sensitive to GA₃, which plays a significant regulatory role in regulation of flower development, berry expansion, berry set, berry ripening, and seedlessness induction; further, it was found that grapevine miR061 (VvmiR061) is a GA₃ responsive miRNA. In this study, grapevine REV (VvREV) and HOX32 (VvHOX32), two target genes of VvmiR061, were predicted, verified, and cloned; homologous conservation was analyzed in various plants. The expression profiles of both VvmiR061 and its target genes (VvREV and VvHOX32) under GA₃ treatment were detected by qRT-PCR during grapevine flower and berry development. Results revealed that GA₃ treatment has upregulated the transcription of VvREV and VvHOX32, while it downregulated the expression of VvmiR061. The function of VvmiR061 in cleaving target genes VvREV and VvHOX32 was diminished by GA₃ treatment during flower developmental process. The results of this study exhibited the importance of VvmiR061 in regulating flower development and GA₃ signaling pathway and also contributed some to the knowledge of small RNA-mediated regulation in grape.

Response of microRNAs to cold treatment in the young spikes of common wheat

Background

MicroRNAs (miRNAs) are a class of small non-coding RNAs that play important roles in biotic and abiotic stresses by regulating their target genes. For common wheat, spring frost damage frequently occurs, especially when low temperature coincides with plants at early floral organ differentiation, which may result in significant yield loss. Up to date, the role of miRNAs in wheat response to frost stress is not well understood.

Results

We report here the sequencing of small RNA transcriptomes from the young spikes that were treated with cold stress and the comparative analysis with those of the control. A total of 192 conserved miRNAs from 105 families and nine novel miRNAs were identified. Among them, 34 conserved and five novel miRNAs were differentially expressed between the cold-stressed samples and the controls. The expression patterns of 18 miRNAs were further validated by quantitative real time polymerase chain reaction (qRT-PCR). Moreover, nearly half of the miRNAs were cross inducible by biotic and abiotic stresses when compared with previously published work. Target genes were predicted and validated by degradome sequencing. Gene Ontology (GO) enrichment analysis showed that the target genes of differentially expressed miRNAs were enriched for response to the stimulus, regulation of transcription, and ion transport functions. Since many targets of differentially expressed miRNAs were transcription factors that are associated with floral development such as ARF, SPB (Squamosa Promoter Binding like protein), MADS-box (MCM1, AG, DEFA and SRF), MYB, SPX (SYG1, Pho81 and XPR1), TCP (TEOSINTE BRANCHED, Cycloidea and PCF), and PPR (PentatricoPeptide Repeat) genes, cold-altered miRNA expression may cause abnormal reproductive organ development.

Conclusion

Analysis of small RNA transcriptomes and their target genes provide new insight into miRNA regulation in developing wheat inflorescences under cold stress. MiRNAs provide another layer of gene regulation in cold stress response that can be genetically manipulated to reduce yield loss in wheat.

Electronic supplementary material

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

Identification and profiling of conserved and novel microRNAs involved in oil and oleic acid production during embryogenesis in Carya cathayensis Sarg

MicroRNAs (miRNAs) are important regulators of plant development and fruit formation. Mature embryos of hickory (Carya cathayensis Sarg.) nuts contain more than 70% oil (comprising 90% unsaturated fatty acids), along with a substantial amount of oleic acid. To understand the roles of miRNAs involved in oil and oleic acid production during hickory embryogenesis, three small RNA libraries from different stages of embryogenesis were constructed. Deep sequencing of these three libraries identified 95 conserved miRNAs with 19 miRNA*s, 7 novel miRNAs (as well as their corresponding miRNA*s), and 26 potentially novel miRNAs. The analysis identified 15 miRNAs involved in oil and oleic acid production that are differentially expressed during embryogenesis in hickory. Among them, nine miRNA sequences, including eight conserved and one novel, were confirmed by qRT-PCR. In addition, 145 target genes of the novel miRNAs were predicted using a bioinformatic approach. Our results provide a framework for better understanding the roles of miRNAs during embryogenesis in hickory.

Genome-wide analysis of miRNAs and Tasi-RNAs in Zea mays in response to phosphate deficiency

Globally important cereal crop maize provides important nutritions and starch in dietary foods. Low phosphate (LPi) availability in the soil frequently limits the maize quality and yield across the world. Small non-coding RNAs (Snc-RNAs) play crucial roles in growth and adaptation of plants to the environment. Snc-RNAs like microRNAs (miRs) and trans-acting small interfering RNAs (Tasi-Rs) play important functions in posttranscriptional regulation of gene expression, which controls plant development, reproduction, and biotic/abiotic stress responses. In order to identify the miR and Tasi-R alterations in leaf and root of maize in response to sufficient phosphate and LPi at 3LS and 4LS, the snc-RNA population libraries for 0th, 1st, 2nd, 4th, and 8th day were constructed. These libraries were used for genome-wide alignment and RNA-fold analysis for possible prediction of potential miRs and Tasi-Rs. This study reported 174 known and conserved differentially expressed miRs of 27 miR families of maize plant. In addition, leaf and root specific potential novel miRs representing 155 new families were also discovered. Differentially expressed conserved as well as novel miR functions in root and leaf during early stage of Pi starvation were extensively discussed. Leaf and root specific miRs as well as common miRs with their target genes, participating in different biological, cellular, and metabolic processes were explored. Further, four miR390-directed Tasi-Rs which belong to TAS3 gene family along with other orthologs of Tasi-Rs were also identified. Finally, the study provides an insight into the composite regulatory mechanism of miRs in maize in response to Pi deficiency.

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.

A systemic identification approach for primary transcription start site of Arabidopsis miRNAs from multidimensional omics data

The 22-nucleotide non-coding microRNAs (miRNAs) are mostly transcribed by RNA polymerase II and are similar to protein-coding genes. Unlike the clear process from stem-loop precursors to mature miRNAs, the primary transcriptional regulation of miRNA, especially in plants, still needs to be further clarified, including the original transcription start site, functional cis-elements and primary transcript structures. Due to several well-characterized transcription signals in the promoter region, we proposed a systemic approach integrating multidimensional "omics" (including genomics, transcriptomics, and epigenomics) data to improve the genome-wide identification of primary miRNA transcripts. Here, we used the model plant Arabidopsis thaliana to improve the ability to identify candidate promoter locations in intergenic miRNAs and to determine rules for identifying primary transcription start sites of miRNAs by integrating high-throughput omics data, such as the DNase I hypersensitive sites, chromatin immunoprecipitation-sequencing of polymerase II and H3K4me3, as well as high throughput transcriptomic data. As a result, 93% of refined primary transcripts could be confirmed by the primer pairs from a previous study. Cis-element and secondary structure analyses also supported the feasibility of our results. This work will contribute to the primary transcriptional regulatory analysis of miRNAs, and the conserved regulatory pattern may be a suitable miRNA characteristic in other plant species.

Water-deficit stress-responsive microRNAs and their targets in four durum wheat genotypes

MicroRNAs (miRNAs) guide regulation at the post-transcriptional level by inducing messenger RNA (mRNA) degradation or translational inhibition of their target protein-coding genes. Durum wheat miRNAs may contribute to the genotypic water-deficit stress response in different durum varieties. Further investigation of the interactive miRNA-target regulatory modules and experimental validation of their response to water stress will contribute to our understanding of the small RNA-mediated molecular networks underlying stress adaptation in durum wheat. In this study, a comprehensive genome-wide in silico analysis using the updated Triticum transcriptome assembly identified 2055 putative targets for 113 conserved durum miRNAs and 131 targets for four novel durum miRNAs that putatively contribute to genotypic stress tolerance. Predicted mRNA targets encode various transcription factors, binding proteins and functional enzymes, which play vital roles in multiple biological pathways such as hormone signalling and metabolic processes. Quantitative PCR profiling further characterised 43 targets and 5 miRNAs with stress-responsive and/or genotype-dependent differential expression in two stress-tolerant and two stress-sensitive durum genotypes subjected to pre-anthesis water-deficit stress. Furthermore, a 5' RLM-RACE approach validated nine mRNA targets cleaved by water-deficit stress-responsive miRNAs, which, to our knowledge, has not been previously reported in durum wheat. The present study provided experimental evidence of durum miRNAs and target genes in response to water-deficit stress in contrasting durum varieties, providing new insights into the regulatory roles of the miRNA-guided RNAi mechanism underlying stress adaptation in durum wheat.

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.