Subgenomic analysis of microRNAs in polyploid wheat

Comparative Analysis of Coding and Non-Coding Features within Insect Tolerance Loci in Wheat with Their Homologs in Cereal Genomes

Food insecurity and malnutrition have reached critical levels with increased human population, climate fluctuations, water shortage; therefore, higher-yielding crops are in the spotlight of numerous studies. Abiotic factors affect the yield of staple food crops; among all, wheat stem sawfly (Cephus cinctus Norton) and orange wheat blossom midge (Sitodiplosis mosellana) are two of the most economically and agronomically harmful insect pests which cause yield loss in cereals, especially in wheat in North America. There is no effective strategy for suppressing this pest damage yet, and only the plants with intrinsic tolerance mechanisms such as solid stem phenotypes for WSS and antixenosis and/or antibiosis mechanisms for OWBM can limit damage. A major QTL and a causal gene for WSS resistance were previously identified in wheat, and 3 major QTLs and a causal gene for OWBM resistance. Here, we present a comparative analysis of coding and non-coding features of these loci of wheat across important cereal crops, barley, rye, oat, and rice. This research paves the way for our cloning and editing of additional WSS and OWBM tolerance gene(s), proteins, and metabolites.

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.

Chromosome genomics uncovers plant genome organization and function

The identification of causal genomic loci and their interactions underlying various traits in plants has been greatly aided by progress in understanding the organization of the nuclear genome. This provides clues to the responses of plants to environmental stimuli at the molecular level. Apart from other uses, these insights are needed to fully explore the potential of new breeding techniques that rely on genome editing. However, genome analysis and sequencing is not straightforward in the many agricultural crops and their wild relatives that possess large and complex genomes. Chromosome genomics streamlines this task by dissecting the genome to single chromosomes whose DNA is then used instead of nuclear DNA. This results in a massive and lossless reduction in DNA sample complexity, reduces the time and cost of the experiment, and simplifies data interpretation. Flow cytometric sorting of condensed mitotic chromosomes makes it possible to purify single chromosomes in large quantities, and as the DNA remains intact this process can be coupled successfully with many techniques in molecular biology and genomics. Since the first experiments with flow cytometric sorting in the late 1980s, numerous applications have been developed, and chromosome genomics has been having a significant impact in many areas of research, including the sequencing of complex genomes of important crops and gene cloning. This review discusses these applications, describes their contribution to advancements in plant genome analysis and gene cloning, and outlines future directions.

Evolution of small RNA expression following hybridization and allopolyploidization: insights from Spartina species (Poaceae, Chloridoideae)

Differential expression of mi-RNAs targeting developmental processes and progressive downregulation of repeat-associated siRNAs following genome merger and genome duplication in the context of allopolyploid speciation in Spartina. The role of small RNAs on gene expression regulation and genome stability is arousing increased interest and is being explored in various plant systems. In spite of prominence of reticulate evolution and polyploidy that affects the evolutionary history of all plant lineages, very few studies analysed RNAi mechanisms with this respect. Here, we explored small RNAs diversity and expression in the context of recent allopolyploid speciation, using the Spartina system, which offers a unique opportunity to explore the immediate changes following hybridization and genome duplication. Small RNA-Seq analyses were conducted on hexaploid parental species (S. alterniflora and S. maritima), their F1 hybrid S. x townsendii, and the neoallododecaploid S. anglica. We identified 594 miRNAs, 2197 miRNA-target genes, and 3730 repeat-associated siRNAs (mostly targeting Class I/Copia-Ivana- Copia-SIRE and LINEs elements). For both mi- and ra-siRNAs, we detected differential expression patterns following genome merger and genome duplication. These misregulations include non-additive expression of miRNAs in the F1 hybrid and additional changes in the allopolyploid targeting developmental processes. Expression of repeat-associated siRNAs indicates a strengthen of transposable element repression during the allopolyploidization process. Altogether, these results confirm the central role small RNAs play in shaping regulatory changes in naturally formed recent allopolyploids.

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.

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.

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.

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.

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.

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.

Synteny and comparative analysis of miRNA retention, conservation, and structure across Brassicaceae reveals lineage- and sub-genome-specific changes

The recent availability of genome sequences together with syntenic block information for Brassicaceae offers an opportunity to study microRNA (miRNA) evolution across this family. We employed a synteny-based comparative genomics strategy to unambiguously identify miRNA homologs from the genome sequence of members of Brassicaceae. Such an analysis of miRNA across Brassicaceae allowed us to classify miRNAs as conserved, lineage-, karyotype- and sub-genome-specific. The differential loss of miRNA from sub-genomes in polyploid genomes of Brassica rapa and Brassica oleracea shows that miRNA also follows the rules of gene fractionation as observed in the case of protein-coding genes. The study of mature and miR* region of precursors revealed instances of in-dels and SNPs which reflect the evolutionary history of the genomes. High level of conservation in miR* regions in some cases points to their functional relevance which needs to be further investigated. We further show that sequence and length variability in precursor sequences can affect the free energy and foldback structure of miRNA which may ultimately affect their biogenesis and expression in the biological system.

Identification and Comparative Analysis of microRNA in Wheat (Triticum aestivum L.) Callus Derived from Mature and Immature Embryos during In vitro Culture

Feasible and efficient tissue culture plays an important role in plant genetic engineering. Wheat (Triticum aestivum L.) immature embryos (IMEs) are preferred for tissue culture to mature embryos (MEs) because IMEs easily generate embryogenic callus, producing large number of plants. The molecular mechanisms of regulation and the biological pathways involved in embryogenic callus formation in wheat remain unclear. Here, microRNAs (miRNAs) potentially involved in embryogenic callus formation and somatic embryogenesis were identified through deep sequencing of small RNAs (sRNAs) and analyzed with bioinformatics tools. Six sRNA libraries derived from calli of IMEs and MEs after 3, 6, or 15 d of culture (DC) were constructed and sequenced. A total of 85 known miRNAs were identified, of which 30, 33, and 18 were differentially expressed (P < 0.05) between the IME and ME libraries at 3, 6, and 15 DC, respectively. Additionally, 171 novel and 41 candidate miRNAs were also identified, of the novel miRNA, 69, 67, and 37 were differentially expressed (P < 0.05) between the two types of libraries at 3, 6, and 15 DC, respectively. The expression patterns of eight known and eight novel miRNAs were validated using quantitative real-time polymerase chain reaction. Gene ontology annotation of differentially expressed miRNA targets provided information regarding the underlying molecular functions, biological processes, and cellular components involved in embryogenic callus development. Functional miRNAs, such as miR156, miR164, miR1432, miR398, and miR397, differentially expressed in IMEs and MEs might be related to embryogenic callus formation and somatic embryogenesis. This study suggests that miRNA plays an important role in embryogenic callus formation and somatic embryogenesis in wheat, and our data provide a useful resource for further research.

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.

Flow Sorting Plant Chromosomes

Nuclear genomes of many important plant species are tremendously complicated to map and sequence. The ability to isolate single chromosomes, which represent small units of nuclear genome, is priceless in many areas of plant research including cytogenetics, genomics, and proteomics. Flow cytometry is the only technique which can provide large quantities of pure chromosome fractions suitable for downstream applications including physical mapping, preparation of chromosome-specific BAC libraries, sequencing, and optical mapping. Here, we describe step-by-step procedure of preparation of liquid suspensions of intact mitotic metaphase chromosomes and their flow cytometric analysis and sorting.

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.

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.

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.

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

Background

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

Results

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

Conclusion

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

Electronic supplementary material

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

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

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

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

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

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

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

New wheat microRNA using whole-genome sequence

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

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

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

Drought tolerance in modern and wild wheat

The genus Triticum includes bread (Triticum aestivum) and durum wheat (Triticum durum) and constitutes a major source for human food consumption. Drought is currently the leading threat on world's food supply, limiting crop yield, and is complicated since drought tolerance is a quantitative trait with a complex phenotype affected by the plant's developmental stage. Drought tolerance is crucial to stabilize and increase food production since domestication has limited the genetic diversity of crops including wild wheat, leading to cultivated species, adapted to artificial environments, and lost tolerance to drought stress. Improvement for drought tolerance can be achieved by the introduction of drought-grelated genes and QTLs to modern wheat cultivars. Therefore, identification of candidate molecules or loci involved in drought tolerance is necessary, which is undertaken by "omics" studies and QTL mapping. In this sense, wild counterparts of modern varieties, specifically wild emmer wheat (T. dicoccoides), which are highly tolerant to drought, hold a great potential. Prior to their introgression to modern wheat cultivars, drought related candidate genes are first characterized at the molecular level, and their function is confirmed via transgenic studies. After integration of the tolerance loci, specific environment targeted field trials are performed coupled with extensive analysis of morphological and physiological characteristics of developed cultivars, to assess their performance under drought conditions and their possible contributions to yield in certain regions. This paper focuses on recent advances on drought related gene/QTL identification, studies on drought related molecular pathways, and current efforts on improvement of wheat cultivars for drought tolerance.

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

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