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

A Survey of the Transcriptomic Resources in Durum Wheat: Stress Responses, Data Integration and Exploitation

Durum wheat (Triticum turgidum subsp. durum (Desf.) Husn.) is an allotetraploid cereal crop of worldwide importance, given its use for making pasta, couscous, and bulgur. Under climate change scenarios, abiotic (e.g., high and low temperatures, salinity, drought) and biotic (mainly exemplified by fungal pathogens) stresses represent a significant limit for durum cultivation because they can severely affect yield and grain quality. The advent of next-generation sequencing technologies has brought a huge development in transcriptomic resources with many relevant datasets now available for durum wheat, at various anatomical levels, also focusing on phenological phases and environmental conditions. In this review, we cover all the transcriptomic resources generated on durum wheat to date and focus on the corresponding scientific insights gained into abiotic and biotic stress responses. We describe relevant databases, tools and approaches, including connections with other "omics" that could assist data integration for candidate gene discovery for bio-agronomical traits. The biological knowledge summarized here will ultimately help in accelerating durum wheat breeding.

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

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

Temporal expression study of miRNAs in the crown tissues of winter wheat grown under natural growth conditions

BACKGROUND

Winter wheat requires prolonged exposure to low temperature to initiate flowering (vernalization). Shoot apical meristem of the crown is the site of cold perception, which produces leaf primordia during vegetative growth before developing into floral primordia at the initiation of the reproductive phase. Although many essential genes for winter wheat cold acclimation and floral initiation have been revealed, the importance of microRNA (miRNA) meditated post-transcriptional regulation in crowns is not well understood. To understand the potential roles of miRNAs in crown tissues, we performed a temporal expression study of miRNAs in crown tissues at the three-leaf stage, winter dormancy stage, spring green-up stage, and jointing stage of winter wheat grown under natural growth conditions.

RESULTS

In total, 348 miRNAs belonging to 298 miRNA families, were identified in wheat crown tissues. Among them, 92 differentially expressed miRNAs (DEMs) were found to be significantly regulated from the three-leaf stage to the jointing stage. Most of these DEMs were highly expressed at the three-leaf stage and winter dormancy stage, and then declined in later stages. Six DEMs, including miR156a-5p were markedly induced during the winter dormancy stage. Eleven DEMs, including miR159a.1, miR390a-5p, miR393-5p, miR160a-5p, and miR1436, were highly expressed at the green-up stage. Twelve DEMs, such as miR172a-5p, miR394a, miR319b-3p, and miR9676-5p were highly induced at the jointing stage. Moreover, 14 novel target genes of nine wheat or Pooideae-specific miRNAs were verified using RLM-5' RACE assay. Notably, six mTERFs and two Rf1 genes, which are associated with mitochondrial gene expression, were confirmed as targets of three wheat-specific miRNAs.

CONCLUSIONS

The present study not only confirmed the known miRNAs associated with phase transition and floral development, but also identified a number of wheat or Pooideae-specific miRNAs critical for winter wheat cold acclimation and floral development. Most importantly, this study provided experimental evidence that miRNA could regulate mitochondrial gene expression by targeting mTERF and Rf1 genes. Our study provides valuable information for further exploration of the mechanism of miRNA mediated post-transcriptional regulation during winter wheat vernalization and inflorescent initiation.

microRNA 166: an evolutionarily conserved stress biomarker in land plants targeting HD-ZIP family

MicroRNAs (miRNAs) are significant class of noncoding RNAs having analytical investigating and modulatory roles in various signaling mechanisms in plants related to growth, development and environmental stress. Conserved miRNAs are an affirmation of land plants evolution and adaptation. They are a proof of indispensable roles of endogenous gene modulators that mediate plant survival on land. Out of such conserved miRNA families, is one core miRNA known as miR166 that is highly conserved among land plants. This particular miRNA is known to primarily target HD ZIP-III transcription factors. miR166 has roles in various developmental processes, as well as regulatory roles against biotic and abiotic stresses in major crop plants. Major developmental roles indirectly modulated by miR166 include shoot apical meristem and vascular differentiation, leaf and root development. In terms of abiotic stress, it has decisive regulatory roles under drought, salinity, and temperature along with biotic stress management. miR166 and its target genes are also known for their beneficial synergy with microorganisms in leguminous crops in relation to lateral roots and nodule development. Hence it is important to study the roles of miR166 in different crop plants to understand its defensive roles against environmental stresses and improve plant productivity by reprogramming several gene functions at molecular levels. This review is hence a summary of different regulatory roles of miR166 with its target HD-ZIP III and its modulatory and fine tuning against different environmental stresses in various plants.

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

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

Nitrogen Starvation-Responsive MicroRNAs Are Affected by Transgenerational Stress in Durum Wheat Seedlings

Stress events have transgenerational effects on plant growth and development. In Mediterranean regions, water-deficit and heat (WH) stress is a frequent issue that negatively affects crop yield and quality. Nitrogen (N) is an essential plant macronutrient and often a yield-limiting factor for crops. Here, the response of durum wheat seedlings to N starvation under the transgenerational effects of WH stress was investigated in two genotypes. Both genotypes showed a significant reduction in seedling height, leaf number, shoot and root weight (fresh and dry), primary root length, and chlorophyll content under N starvation stress. However, in the WH stress-tolerant genotype, the percentage reduction of most traits was lower in progeny from the stressed parents than progeny from the control parents. Small RNA sequencing identified 1534 microRNAs in different treatment groups. Differentially expressed microRNAs (DEMs) were characterized subject to N starvation, parental stress and genotype factors, with their target genes identified in silico. GO and KEGG enrichment analyses revealed the biological functions, associated with DEM-target modules in stress adaptation processes, that could contribute to the phenotypic differences observed between the two genotypes. The study provides the first evidence of the transgenerational effects of WH stress on the N starvation response in durum wheat.

Blocking miR530 Improves Rice Resistance, Yield, and Maturity

MicroRNAs fine-tune plant growth and resistance against multiple biotic and abiotic stresses. The trade-off between biomass and resistance can penalize crop yield. In this study, we have shown that rice miR530 regulates blast disease resistance, yield, and growth period. While the overexpression of miR530 results in compromised blast disease resistance, reduced grain yield, and late maturity, blocking miR530 using a target mimic (MIM530) leads to enhanced resistance, increased grain yield, and early maturity. Further study revealed that the accumulation of miR530 was decreased in both leaves and panicles along with the increase of age. Such expression patterns were accordant with the enhanced resistance from seedlings to adult plants, and the grain development from panicle formation to fully-filled seeds. Divergence analysis of miR530 precursor with upstream 1,000-bp promoter sequence in 11 rice species revealed that miR530 was diverse in Oryza sativa japonica and O. sativa indica group, which was consistent with the different accumulation of miR530 in japonica accessions and indica accessions. Altogether, our results indicate that miR530 coordinates rice resistance, yield, and maturity, thus providing a potential regulatory module for breeding programs aiming to improve yield and disease resistance.

Identification, evolution, expression, and docking studies of fatty acid desaturase genes in wheat (Triticum aestivum L.)

Backgrounds

Fatty acid desaturases (FADs) introduce a double bond into the fatty acids acyl chain resulting in unsaturated fatty acids that have essential roles in plant development and response to biotic and abiotic stresses. Wheat germ oil, one of the important by-products of wheat, can be a good alternative for edible oils with clinical advantages due to the high amount of unsaturated fatty acids. Therefore, we performed a genome-wide analysis of the wheat FAD gene family (TaFADs).

Results

68 FAD genes were identified from the wheat genome. Based on the phylogenetic analysis, wheat FADs clustered into five subfamilies, including FAB2, FAD2/FAD6, FAD4, DES/SLD, and FAD3/FAD7/FAD8. The TaFADs were distributed on chromosomes 2A-7B with 0 to 10 introns. The Ka/Ks ratio was less than one for most of the duplicated pair genes revealed that the function of the genes had been maintained during the evolution. Several cis-acting elements related to hormones and stresses in the TaFADs promoters indicated the role of these genes in plant development and responses to environmental stresses. Likewise, 72 SSRs and 91 miRNAs in 36 and 47 TaFADs have been identified. According to RNA-seq data analysis, the highest expression in all developmental stages and tissues was related to TaFAB2.5, TaFAB2.12, TaFAB2.15, TaFAB2.17, TaFAB2.20, TaFAD2.1, TaFAD2.6, and TaFAD2.8 genes while the highest expression in response to temperature stress was related to TaFAD2.6, TaFAD2.8, TaFAB2.15, TaFAB2.17, and TaFAB2.20. Furthermore, docking simulations revealed several residues in the active site of TaFAD2.6 and TaFAD2.8 in close contact with the docked oleic acid that could be useful in future site-directed mutagenesis studies to increase the catalytic efficiency of them and subsequently improve agronomic quality and tolerance of wheat against environmental stresses.

Conclusions

This study provides comprehensive information that can lead to the detection of candidate genes for wheat genetic modification.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07199-1.

Integrated Analysis of Small RNA, Transcriptome, and Degradome Sequencing Reveals the Water-Deficit and Heat Stress Response Network in Durum Wheat

Water-deficit and heat stress negatively impact crop production. Mechanisms underlying the response of durum wheat to such stresses are not well understood. With the new durum wheat genome assembly, we conducted the first multi-omics analysis with next-generation sequencing, providing a comprehensive description of the durum wheat small RNAome (sRNAome), mRNA transcriptome, and degradome. Single and combined water-deficit and heat stress were applied to stress-tolerant and -sensitive Australian genotypes to study their response at multiple time-points during reproduction. Analysis of 120 sRNA libraries identified 523 microRNAs (miRNAs), of which 55 were novel. Differentially expressed miRNAs (DEMs) were identified that had significantly altered expression subject to stress type, genotype, and time-point. Transcriptome sequencing identified 49,436 genes, with differentially expressed genes (DEGs) linked to processes associated with hormone homeostasis, photosynthesis, and signaling. With the first durum wheat degradome report, over 100,000 transcript target sites were characterized, and new miRNA-mRNA regulatory pairs were discovered. Integrated omics analysis identified key miRNA-mRNA modules (particularly, novel pairs of miRNAs and transcription factors) with antagonistic regulatory patterns subject to different stresses. GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis revealed significant roles in plant growth and stress adaptation. Our research provides novel and fundamental knowledge, at the whole-genome level, for transcriptional and post-transcriptional stress regulation in durum wheat.

Identification and Temporal Expression Analysis of Conserved and Novel MicroRNAs in the Leaves of Winter Wheat Grown in the Field

Cold acclimation and vegetative/reproductive transition are two important evolutionary adaptive mechanisms for winter wheat surviving the freezing temperature in winter and successful seeds setting in the next year. MicroRNA (miRNA) is a class of regulatory small RNAs (sRNAs), which plays critical roles in the growth and development of plants. However, the regulation mechanism of miRNAs during cold acclimation and vegetative/reproductive transition of winter wheat is not much understood. In this study, four sRNA libraries from leaves of winter wheat grown in the field at the three-leaf stage, winter dormancy stage, spring green-up stage, and jointing stage were analyzed to identify known and novel miRNAs and to understand their potential roles in the growth and development of winter wheat. We examined miRNA expression using a high-throughput sequencing technique. A total of 373 known, 55 novel, and 27 putative novel miRNAs were identified. Ninety-one miRNAs were found to be differentially expressed at the four stages. Among them, the expression of six known and eight novel miRNAs was significantly suppressed at the winter dormancy stage, whereas the expression levels of seven known and eight novel miRNAs were induced at this stage; three known miRNAs and three novel miRNAs were significantly induced at the spring green-up stage; six known miRNAs were induced at the spring green-up stage and reached the highest expression level at the jointing stage; and 20 known miRNAs and 10 novel miRNAs were significantly induced at the jointing stage. Expression of a number of representative differentially expressed miRNAs was verified using quantitative real-time polymerase chain reaction (qRT-PCR). Potential target genes for known and novel miRNAs were predicted. Moreover, six novel target genes for four Pooideae species-specific miRNAs and two novel miRNAs were verified using the RNA ligase-mediated 5'-rapid amplification of cDNA ends (RLM-5'RACE) technique. These results indicate that miRNAs are key non-coding regulatory factors modulating the growth and development of wheat. Our study provides valuable information for in-depth understanding of the regulatory mechanism of miRNAs in cold acclimation and vegetative/reproductive transition of winter wheat grown in the field.

Identification of microRNA-target modules from rice variety Pusa Basmati-1 under high temperature and salt stress

High temperature and salinity stress are major factors limiting the growth and productivity of rice crop on a global scale. It is therefore an essential prerequisite to understand the molecular genetic regulation of plant responses to dual stresses. MicroRNAs (miRs) are recognized as key controllers of gene expression which act mainly at the post-transcriptional level to regulate various aspects of plant development. The present study attempts to investigate the miR circuits that are modulated in response to high temperature and salinity stress in rice. To gain insights into the pathway, preliminary miR profiles were generated using the next-generation sequencing (NGS) datasets. The identified molecules were filtered on the basis of fold differential regulation under high temperature, and time kinetics of their expression under the two individual stresses was followed to capture the regulatory windows. The analysis revealed the involvement of common miR regulatory nodes in response to two different abiotic stresses, thereby broadening our perspective about the stress-mediated regulatory mechanisms operative in rice.

New insights into tomato microRNAs

Cultivated tomato, Solanum lycopersicum, is one of the most common fruits in the global food industry. Together with the wild tomato Solanum pennellii, it is widely used for developing better cultivars. MicroRNAs affect mRNA regulation, inhibiting its translation and/or promoting its degradation. Important proteins involved in these processes are ARGONAUTE and DICER. This study aimed to identify and characterize the genes involved in the miRNA processing pathway, miRNA molecules and target genes in both species. We validated the presence of pathway genes and miRNA in different NGS libraries and 6 miRNA families using quantitative RT-PCR. We identified 71 putative proteins in S. lycopersicum and 108 in S. pennellii likely involved in small RNAs processing. Of these, 29 and 32 participate in miRNA processing pathways, respectively. We identified 343 mature miRNAs, 226 pre-miRNAs in 87 families, including 192 miRNAs, which were not previously identified, belonging to 38 new families in S. lycopersicum. In S. pennellii, we found 388 mature miRNAs and 234 pre-miRNAs contained in 85 families. All miRNAs found in S. pennellii were unpublished, being identified for the first time in our study. Furthermore, we identified 2471 and 3462 different miRNA target in S. lycopersicum and S. pennellii, respectively.

Identification of miRNAs linked with the drought response of tef [Eragrostis tef (Zucc.) Trotter]

Tef [Eragrostis tef (Zucc.) Trotter], a staple food crop in the Horn of Africa and particularly in Ethiopia, has several beneficial agronomical and nutritional properties, including waterlogging and drought tolerance. In this study, we performed microRNA profiling of tef using the Illumina HiSeq 2500 platform, analyzing both shoots and roots of two tef genotypes, one drought-tolerant (Tsedey) and one drought-susceptible (Alba). We obtained more than 10 million filtered reads for each of the 24 sequenced small cDNA libraries. Reads mapping to known miRNAs were more abundant in the root than shoot tissues. Thirteen and 35 miRNAs were significantly modulated in response to drought, in Alba and Tsedey roots, respectively. One miRNA was upregulated under drought conditions in both genotypes. In shoots, nine miRNAs were modulated in common between the two genotypes and all showed similar trends of expression. One-hundred and forty-seven new miRNA mature sequences were identified in silico, 22 of these were detected in all relevant samples and seven were differentially regulated when comparing drought with normal watering. Putative targets of the miRNA regulated under drought in root and shoot tissues were predicted. Among the targets were transcription factors such as CCAAT-HAP2, MADS and NAC. Verification with qRT-PCR revealed that five of six potential targets showed a pattern of expression that was consistent with the correspondent miRNA amount measured by RNA-Seq. In general, candidate miRNAs involved in the post-transcriptional regulation of the tef response to drought could be included in next-generation breeding programs.

The contrasting microRNA content of a drought tolerant and a drought susceptible wheat cultivar

Drought stress represents one of the most common stresses affecting the productivity of crop plants. A rather recently discovered component of the plant response to drought is the cellular population of microRNAs. Here, the microRNA content was revealed of two bread wheat cultivars contrasting strongly with respect to the ability to withstand drought stress. A total of 1813 miRNAs was identified, grouped into 106 families. Some 104 of these miRNAs were predicted to match 212 novel miRNA precursors. In the drought tolerant cultivar (SM), 105 (33 known and 72 novel) miRNAs were altered in abundance by the imposition of drought stress, while the equivalent number in the more sensitive cultivar (SW) was 51 (20 and 31). An in silico analysis predicted that these miRNAs target at least 1959 genes in SM and 1111 in SW, suggesting their broad contribution to the drought stress response. Among the target genes were several known stress-related genes, encoding, for example, superoxide dismutase, various MYB transcription factors, various ABA signaling proteins and various MADS-box transcription factors. In many cases, the more susceptible cultivar SW behaved in a contrasting manner. The suggestion is that miRNAs represent an important aspect of the drought stress response, post-transcriptionally regulating a range of stress-related genes.