A comprehensive genome-wide study on tissue-specific and abiotic stress-specific miRNAs in Triticum aestivum

Cell-Type Specific miRNA Regulatory Network Responses to ABA Stress Revealed by Time Series Transcriptional Atlases in Arabidopsis

In plants, microRNAs (miRNAs) participate in complex gene regulatory networks together with the transcription factors (TFs) in response to biotic and abiotic stresses. To date, analyses of miRNAs-induced transcriptome remodeling are at the whole plant or tissue levels. Here, Arabidopsis's ABA-induced single-cell RNA-seq (scRNA-seq) is performed at different stages of time points-early, middle, and late. Single-cell level primary miRNAs (pri-miRNAs) atlas supported the rapid, dynamic, and cell-type specific miRNA responses under ABA treatment. MiRNAs respond rapidly and prior to target gene expression dynamics, and these rapid response miRNAs are highly cell-type specific, especially in mesophyll and vascular cells. MiRNA-TF-mRNA regulation modules are identified by identifying miRNA-contained feed-forward loops (M-FFLs) in the regulatory network, and regulatory networks with M-FFLs have higher co-expression and clustering coefficient (CC) values than those without M-FFLs, suggesting the hub role of miRNAs in regulatory networks. The cell-type-specific M-FFLs are regulated by these hub miRNAs rather than TFs through sc-RNA-seq network analysis. MiR858a-FBH3-MYB module inhibited the expression of MYB63 and MYB20, which related to the formation of plant secondary wall and the production of lignin, through M-FFL specifically in vascular. These results can provide prominent insights into miRNAs' dynamic and cell-type-specific roles in plant development and stress responses.

Characterization of tae-miR156(s) and their response to abiotic stress in wheat (Triticum aestivum L.)

The microRNA156 (miR156) has been widely studied in plants, however, the characterization of the miR156 family of genes in wheat and their expression patterns under abiotic stress are not completely clear. In this study, a total of 20 miR156 family members, referred to as tae-miR156a to tae-miR156t, were identified in wheat with their loci mapped to various chromosomes. These members were divided into five subgroups: miR156a/b/c/d/e/f, miR156g/h/i, miR156j/k, miR156l/m/n/o/p/q, and miR156r/s/t. They were highly conserved during evolution. The prediction of cis-elements in the tae-MIR156(s) promoter region revealed that the tae-MIR156(s) had diverse cis-acting elements; of these, 15 tae-MIR156(s) and 6 tae-MIR156(s) were found to be drought-responsive elements and cold-responsive elements, respectively. And the prediction target genes of tae-miR156(s) are mainly SPL transcription factor genes. Expression analysis based on quantitative real-time polymerase chain reaction (qRT‒PCR) showed that miR156(s) have different expression levels in the various wheat tissues, and the subgroups' response to abiotic stress varied. Among them, miR156g/h/i were strongly induced in the root of cold and heat stress, and miR156a/b/c/d/e/f were significantly increased in roots after drought stress, whereas miR156r/s/t were highly inhibited in leaves and roots after salt stress. These findings imply that tae-miR156(s) are involved in stress response in wheat, and they provide new fundamental knowledge for further analysis of the function of miR156 and its regulatory mechanism in response to abiotic stress.

Expression patterns of candidate genes for the Lr46/Yr29 "slow rust" locus in common wheat (Triticum aestivum L.) and associated miRNAs inform of the gene conferring the Puccinia triticina resistance trait

Leaf rust caused by Puccinia triticina (Pt) is one of the most impactful diseases causing substantial losses in common wheat (Triticum aestivum L.) crops. In adult plants resistant to Pt, a horizontal adult plant resistance (APR) is observed: APR protects the plant against multiple pathogen races and is distinguished by durable persistence under production conditions. The Lr46/Yr29 locus was mapped to chromosome 1B of common wheat genome, but the identity of the underlying gene has not been demonstrated although several candidate genes have been proposed. This study aimed to analyze the expression of nine candidate genes located at the Lr46/Yr29 locus and their four complementary miRNAs (tae-miR5384-3p, tae-miR9780, tae-miR9775, and tae-miR164), in response to Pt infection. The plant materials tested included five reference cultivars in which the molecular marker csLV46G22 associated with the Lr46/Yr29-based Pt resistance was identified, as well as one susceptible control cultivar. Biotic stress was induced in adult plants by inoculation with fungal spores under controlled conditions. Plant material was sampled before and at 6, 12, 24, 48 hours post inoculation (hpi). Differences in expression of candidate genes at the Lr46/Yr29 locus were analyzed by qRT-PCR and showed that the expression of the genes varied at the analyzed time points. The highest expression of Lr46/Yr29 candidate genes (Lr46-Glu1, Lr46-Glu2, Lr46-Glu3, Lr46-RLK1, Lr46-RLK2, Lr46-RLK3, Lr46-RLK4, Lr46-Snex, and Lr46-WRKY) occurred at 12 and 24 hpi and such expression profiles were obtained only for one candidate gene among the nine genes analyzed (Lr46-Glu2), indicating that it may be a contributing factor in the resistance response to Pt infection.

Identification and characterization of sulphotransferase (SOT) genes for tolerance against drought and heat in wheat and six related species

BACKGROUND

Sulphotransferase (SOT) enzyme (encoded by a conserved family of SOT genes) is involved in sulphonation of a variety of compounds, through transfer of a sulphuryl moiety from 3'phosphoadenosine- 5'phosphosulphate (PAPS) to a variety of secondary metabolites. The PAPS itself is derived from 3'adenosine-5'phosphosulphate (APS) that is formed after uptake of sulphate ions from the soil. The process provides tolerance against abiotic stresses like drought and heat in plants. Therefore, a knowledge of SOT genes in any crop may help in designing molecular breeding methods for improvement of tolerance for drought and heat.

METHODS

Sequences of rice SOT genes and SOT domain (PF00685) of corresponding proteins were both used for identification of SOT genes in wheat and six related species (T. urartu, Ae. tauschii, T. turgidum, Z. mays, B. distachyon and Hordeum vulgare), although detailed analysis was conducted only in wheat. The wheat genes were mapped on individual chromosomes and also subjected to synteny and collinearity analysis. The proteins encoded by these genes were examined for the presence of a complete SOT domain using 'Conserved Domain Database' (CDD) search tool at NCBI.

RESULTS

In wheat, 107 TaSOT genes, ranging in length from 969 bp to 7636 bp, were identified and mapped onto individual chromosomes. SSRs (simple sequence repeats), microRNAs, long non-coding RNAs (lncRNAs) and their target sites were also identified in wheat SOT genes. SOT proteins were also studied in detail. An expression assay of TaSOT genes via wheat RNA-seq data suggested engagement of these genes in growth, development and responses to various hormones and biotic/abiotic stresses.

CONCLUSIONS

The results of the present study should help in further functional characterization of SOT genes in wheat and other related crops.

Diversity of Expression Patterns of Lr34, Lr67, and Candidate Genes towards Lr46 with Analysis of Associated miRNAs in Common Wheat Hybrids in Response to Puccinia triticina Fungus

Leaf rust caused by Puccinia triticina (Pt) is one of the most dangerous diseases causing significant losses in common wheat crops. In adult plants resistant to rust, a horizontal adult plant resistance (APR) type is observed, which protects the plant against multiple pathogen races and is distinguished by greater persistence under production conditions. Crucial pleiotropic slow-rust genes such as Lr34, Lr46, Lr67, and Lr68, in combination with other genes of lesser influence, continue to increase durable resistance to rust diseases. Based on our previous results, we selected four candidate genes for Lr46 out of ten candidates and analysed them for expression before and after inoculation by P. triticina. As part of our study, we also investigated the expression patterns of miRNA molecules complementary to Lr34 and the candidate genes. The aim of the study was to analyse the expression profiles of candidate genes for the Lr46 gene and the Lr34 and Lr67 genes responsible for the differential leaf-rust resistance of hybrid forms of the F1 generation resulting from crosses between the Glenlea cultivar and cultivars from Polish breeding companies. In addition, the expression of five miRNAs (tae-miR9653b, tae-miR5384-3p, tae-miR9780, tae-miR9775 and tae-miR164), complementary to Lr34, and selected candidate genes were analysed using stem-loop RT-PCR and ddPCR. Biotic stress was induced in adult plants by inoculation with Pt fungal spores, under controlled conditions. Plant material was collected before and 6, 12, 24, and 48 h after inoculation (hpi). Differences in expression patterns of Lr34, Lr67, and candidate genes (for Lr46) were analysed by qRT-PCR and showed that gene expression changed at the analysed time points. Identification of molecular markers coupled to the Lr genes studied was also carried out to confirm the presence of these genes in wheat hybrids. qRT-PCR was used to examine the expression levels of the resistance genes. The highest expression of Lr46/Yr29 genes (Lr46-Glu2, Lr46-RLK1, Lr46-RLK2, and Lr46-RLK3) occurred at 12 and 24 hpi, and such expression profiles were obtained for only one candidate gene among the four genes analysed (Lr46-Glu2), indicating that it may be involved in resistance mechanisms of response to Pt infection.

Comprehensive, integrative genomic analysis of microRNA expression profiles in different tissues of two wheat cultivars with different traits

Wheat is one of the most important food sources on Earth. MicroRNAs (miRNA) play important roles in wheat productivity. To identify wheat miRNAs, we constructed and sequenced sRNA libraries from leaves and roots of two wheat cultivars (RAC875 and Kukri) with many different traits. Given that available miRNA wheat complement in the plant-specific database PmiREN ( https://pmiren.com ) does not include root tissues and root-associated miRNAs might thus be missing, we performed first the prediction of novel miRNAs using the sRNAbench tool. We found a total of 150 putatively novel miRNA genes with expression of both arms from 289 unique mature sequences and nearly 30% of all miRNA reads in roots corresponded to novel miRNAs. In contrast, this figure in leaves dropped to under 3%, confirming the undersampling of roots in the complement of known miRNAs. By using 120 publicly available wheat datasets, 598 Zea mays small RNA libraries, 64 plant species genomes, wheat degradome library, and functional enrichment analysis, a subset of novel miRNAs were confirmed as bona-fide miRNAs. Of the total 605 miRNAs identified in this study inclusive of 316 known miRNAs, 528 miRNAs were shared by both cultivars, 429 miRNAs were shared by both root tissues and 329 miRNAs were shared by both leaf tissues. In addition, 32 miRNAs were specific to Kukri while 45 miRNAs were specific to RAC875. These miRNAs had diverse functions, such as regulation of gene transcription, protein translation, energy metabolism, and cell cycle progression. Our data provide a genome-wide miRNA expression profile in these two wheat cultivars and help functional studies of wheat genomics.

Rapid sequence and functional diversification of a miRNA superfamily targeting calcium signaling components in seed plants

MicroRNA (miRNA)-directed posttranscriptional gene silencing (miR-PTGS) is an integral component of gene regulatory networks governing plant development and responses to the environment. The sequence homology between Sly-miR4376, a miRNA common to Solanaceae and reported to target autoinhibited Ca²⁺ -ATPase 10 (ACA10) messenger RNA (mRNA) in tomato, and Arabidopsis miR391 (Ath-miR391), previously annotated as a nonconserved member of the deeply conserved miR390 family, has prompted us to revisit the function of Ath-miR391, as well as its regulatory conservation. A combination of genetic, molecular, and bioinformatic analyses revealed a hidden conservation for miR-PTGS of ACA10 homologs in spermatophytes. We found that the Arabidopsis ACA10 mRNA undergoes miR391-directed cleavage in vivo. Furthermore, transgenic overexpression of miR391 recapitulated the compact inflorescence (cif) phenotypes characteristic of ACA10 loss-of-function mutants, due to miR391-directed PTGS of ACA10. Significantly, comprehensive data mining revealed robust evidence for widespread PTGS of ACA10 homologs directed by a superfamily of related miRNAs sharing a conserved sequence core. Intriguingly, the ACA-targeting miRNAs in Poaceae also direct PTGS for calmodulin-like proteins which are putative Ca²⁺ sensors. The PTGS of ACA10 homologs is therefore directed by a miRNA superfamily that is of ancient origin and has undergone rapid sequence diversification associated with functional innovation.

Small RNAs: The Essential Regulators in Plant Thermotolerance

Small RNAs (sRNAs) are a class of non-coding RNAs that consist of 21-24 nucleotides. They have been extensively investigated as critical regulators in a variety of biological processes in plants. sRNAs include two major classes: microRNAs (miRNAs) and small interfering RNAs (siRNAs), which differ in their biogenesis and functional pathways. Due to global warming, high-temperature stress has become one of the primary causes for crop loss worldwide. Recent studies have shown that sRNAs are involved in heat stress responses in plants and play essential roles in high-temperature acclimation. Genome-wide studies for heat-responsive sRNAs have been conducted in many plant species using high-throughput sequencing. The roles for these sRNAs in heat stress response were also unraveled subsequently in model plants and crops. Exploring how sRNAs regulate gene expression and their regulatory mechanisms will broaden our understanding of sRNAs in thermal stress responses of plant. Here, we highlight the roles of currently known miRNAs and siRNAs in heat stress responses and acclimation of plants. We also discuss the regulatory mechanisms of sRNAs and their targets that are responsive to heat stress, which will provide powerful molecular biological resources for engineering crops with improved thermotolerance.

fIdentification of B. napus small RNAs responsive to infection by a necrotrophic pathogen

BACKGROUND

Small RNAs are short non-coding RNAs that are key gene regulators controlling various biological processes in eukaryotes. Plants may regulate discrete sets of sRNAs in response to pathogen attack. Sclerotinia sclerotiorum is an economically important pathogen affecting hundreds of plant species, including the economically important oilseed B. napus. However, there are limited studies on how regulation of sRNAs occurs in the S. sclerotiorum and B. napus pathosystem.

RESULTS

We identified different classes of sRNAs from B. napus using high throughput sequencing of replicated mock and infected samples at 24 h post-inoculation (HPI). Overall, 3999 sRNA loci were highly expressed, of which 730 were significantly upregulated during infection. These 730 up-regulated sRNAs targeted 64 genes, including disease resistance proteins and transcriptional regulators. A total of 73 conserved miRNA families were identified in our dataset. Degradome sequencing identified 2124 cleaved mRNA products from these miRNAs from combined mock and infected samples. Among these, 50 genes were specific to infection. Altogether, 20 conserved miRNAs were differentially expressed and 8 transcripts were cleaved by the differentially expressed miRNAs miR159, miR5139, and miR390, suggesting they may have a role in the S. sclerotiorum response. A miR1885-triggered disease resistance gene-derived secondary sRNA locus was also identified and verified with degradome sequencing. We also found further evidence for silencing of a plant immunity related ethylene response factor gene by a novel sRNA using 5'-RACE and RT-qPCR.

CONCLUSIONS

The findings in this study expand the framework for understanding the molecular mechanisms of the S. sclerotiorum and B. napus pathosystem at the sRNA level.

Tuning Beforehand: A Foresight on RNA Interference (RNAi) and In Vitro-Derived dsRNAs to Enhance Crop Resilience to Biotic and Abiotic Stresses

Crop yield is severely affected by biotic and abiotic stresses. Plants adapt to these stresses mainly through gene expression reprogramming at the transcriptional and post-transcriptional levels. Recently, the exogenous application of double-stranded RNAs (dsRNAs) and RNA interference (RNAi) technology has emerged as a sustainable and publicly acceptable alternative to genetic transformation, hence, small RNAs (micro-RNAs and small interfering RNAs) have an important role in combating biotic and abiotic stresses in plants. RNAi limits the transcript level by either suppressing transcription (transcriptional gene silencing) or activating sequence-specific RNA degradation (post-transcriptional gene silencing). Using RNAi tools and their respective targets in abiotic stress responses in many crops is well documented. Many miRNAs families are reported in plant tolerance response or adaptation to drought, salinity, and temperature stresses. In biotic stress, the spray-induced gene silencing (SIGS) provides an intelligent method of using dsRNA as a trigger to silence target genes in pests and pathogens without producing side effects such as those caused by chemical pesticides. In this review, we focus on the potential of SIGS as the most recent application of RNAi in agriculture and point out the trends, challenges, and risks of production technologies. Additionally, we provide insights into the potential applications of exogenous RNAi against biotic stresses. We also review the current status of RNAi/miRNA tools and their respective targets on abiotic stress and the most common responsive miRNA families triggered by stress conditions in different crop species.

Identification and characterization of heat-responsive microRNAs at the booting stage in two rice varieties, 9311 and Nagina 22

MicroRNAs (miRNAs) are small, non-coding, regulatory RNAs that play important roles in abiotic stress responses in plants, but their regulatory roles in the adaptive response to heat stress at the booting stage in two rice varieties, 9311 and Nagina 22, remain largely unknown. In this study, 464 known miRNAs and 123 potential novel miRNAs were identified. Of these miRNAs, a total of 90 differentially expressed miRNAs were obtained with 9311 libraries as the control group, of which 54 were upregulated and 36 were downregulated. To gain insight into functional significance, 2773 potential target genes of these 90 differentially expressed miRNAs were predicted. GO enrichment analysis showed that the predicted target genes of differentially expressed miRNAs included NACs, LACs, CSD, and Hsp40. KEGG pathway analysis showed that the target genes of these differentially expressed miRNAs were significantly enriched in the plant hormone signal transduction pathway. The expression levels of 10 differentially expressed miRNAs and their target genes obtained by qRT-PCR were largely consistent with the sequencing results. This study lays a foundation for the elucidation of the miRNA-mediated regulatory mechanisms in rice at elevated temperatures.

Deciphering the role of microRNAs during Pi54 gene mediated Magnaporthe oryzae resistance response in rice

The broad-spectrum resistance gene Pi54 confers resistance to multiple isolates of Magnaporthe oryzae in rice. In order to decipher the molecular mechanism underlying the Pi54 mediated resistance in rice line Taipei309 ^Pi54 (carrying Pi54), miRNAome study was performed at 24 h post-inoculation (hpi) with M. oryzae. A total of 222 known miRNAs representing 101 miRNA families were found in this study. Of these, 29 and 24 miRNAs were respectively up- and down-regulated in the resistant Taipei309 ^Pi54 . Defence response (DR) genes, like, NBSGO35, and OsWAK129b, and genes related to transcription factors were up-regulated in Taipei309 ^Pi54 line. The vast array of miRNA candidates identified here are miR159c, miR167c, miR2100, miR2118o, miR2118l, miR319a, miR393, miR395l, miR397a, miR397b, miR398, miR439g, miR531b, miR812f, and miR815c, and they manifest their role in balancing the interplay between various DR genes during Pi54 mediated resistance. We also validated miRNA/target gene pairs involved in hormone signalling, and cross-talk among hormone pathways regulating the rice immunity. This study suggests that the Pi54 gene mediated blast resistance is influenced by several microRNAs through PTI and ETI components in the rice line Taipei309 ^Pi54 , leading to incompatible host-pathogen interaction.

Northern Blotting Technique for Detection and Expression Analysis of mRNAs and Small RNAs

Northern analysis is a conventional but gold standard method for detection and quantification of gene expression changes. It not only detects the presence of a transcript but also indicates size and relative comparison of transcript abundance on a single membrane. In recent years it has been aptly adapted to validate and study the size and expression of small noncoding RNAs. Here, we describe protocols employed in our laboratory for conventional northern analysis with total RNA/mRNA to study gene expression and validation of small noncoding RNAs using low molecular weight fraction of RNAs. A brief account on the recent advancements for improving the sensitivity and efficiency of northern blot detection is also included in this chapter.

Stress induced dynamic adjustment of conserved miR164:NAC module

Aims including the rationale

Salinity and drought are the two major stresses limiting the productivity of economically important crops such as Glycine max (soybean). The incidence of these stresses during the pod development stages affects the quality and quantity of seeds, which compromise the yield of soybean. The miR164:NAC module has been shown to play a critical role in regulating the response to salt and drought stress in several plant species. However, biological role of miR164:NAC module in salt stress in soybean is not fully understood.

Methods

In this study, we identified 215 salt responsive miRNAs, using miScript miRNA array with a sensitive and a tolerant soybean genotype, William82 and INCASoy36, respectively. The targets of these salt regulated miRNAs were searched in the degradome datasets.

Key results

It was found that four salt stress deregulated miRNAs targeted the NAC transcription factor and among these miR164k and miR408d showed antagonistic expression in the two soybean genotypes. The expression of miR164k was higher in salt tolerant INCASoy36 as compared to salt sensitive William82, under unstressed conditions. However under salt stress, miR164k was downregulated in INCASoy36 (-2.65 fold), whereas it was upregulated in William82 (4.68 fold). A transient co-expression assay validated that gma-miR164k directs the cleavage of GmNAC1 transcript. Bioinformatics analysis revealed that the regulation of NAC transcription factor family by members of miR164 family is conserved across many species. The dynamic expression profiles of miR164 and NAC-TFs were captured in different tissues of rice, tobacco, and two soybean genotypes under drought and salt stress conditions.

Main conclusion

Collectively, our results suggest that genetically determined dynamic modulation of the conserved miR164:NAC-TF module may play an important role in determining the adaptive response of plants to stress.

Research Tools for the Functional Genomics of Plant miRNAs During Zygotic and Somatic Embryogenesis

During early plant embryogenesis, some of the most fundamental decisions on fate and identity are taken making it a fascinating process to study. It is no surprise that higher plant embryogenesis was intensively analysed during the last century, while somatic embryogenesis is probably the most studied regeneration model. Encoded by the MIRNA, short, single-stranded, non-coding miRNAs, are commonly present in all Eukaryotic genomes and are involved in the regulation of the gene expression during the essential developmental processes such as plant morphogenesis, hormone signaling, and developmental phase transition. During the last few years dedicated to miRNAs, analytical methods and tools have been developed, which have afforded new opportunities in functional analyses of plant miRNAs, including (i) databases for in silico analysis; (ii) miRNAs detection and expression approaches; (iii) reporter and sensor lines for a spatio-temporal analysis of the miRNA-target interactions; (iv) in situ hybridisation protocols; (v) artificial miRNAs; (vi) MIM and STTM lines to inhibit miRNA activity, and (vii) the target genes resistant to miRNA. Here, we attempted to summarise the toolbox for functional analysis of miRNAs during plant embryogenesis. In addition to characterising the described tools/methods, examples of the applications have been presented.

Analysis of miRNAs in Two Wheat Cultivars Infected With Puccinia striiformis f. sp. tritici

MicroRNAs are small RNAs that regulate gene expression in eukaryotes. In this study, we analyzed the small RNA profiles of two cultivars that exhibit different reactions to stripe rust infection: one susceptible, the other partially resistant. Using small RNA libraries prepared from the two wheat cultivars infected with stripe rust fungus (Puccinia striiformis f. sp. tritici), we identified 182 previously known miRNAs, 91 variants of known miRNAs, and 163 candidate novel wheat miRNAs. Known miRNA loci were usually copied in all three wheat sub-genomes, whereas novel miRNA loci were often specific to a single sub-genome. DESeq2 analysis of differentially expressed microRNAs revealed 23 miRNAs that exhibit cultivar-specific differences. TA078/miR399b showed cultivar-specific differential regulation in response to infection. Using different target prediction algorithms, 145 miRNAs were predicted to target wheat genes, while 69 miRNAs were predicted to target fungal genes. We also confirmed reciprocal expression of TA078/miR399b and tae-miR9664 and their target genes in different treatments, providing evidence for miRNA-mediated regulation during infection. Both known and novel miRNAs were predicted to target fungal genes, suggesting trans-kingdom regulation of gene expression. Overall, this study contributes to the current repository of wheat miRNAs and provides novel information on the yet-uncharacterized roles for miRNAs in the wheat-stripe rust pathosystem.

Lineage-Specific Evolved MicroRNAs Regulating NB-LRR Defense Genes in Triticeae

Disease resistance genes encoding proteins with nucleotide binding sites and Leucine-Rich Repeat (NB-LRR) domains include many members involved in the effector-triggered immunity pathway in plants. The transcript levels of these defense genes are negatively regulated by diverse microRNAs (miRNAs) in angiosperms and gymnosperms. In wheat, using small RNA expression datasets and degradome datasets, we identified five miRNA families targeting NB-LRR defense genes in monocots, some of which arose in the Triticeae species era. These miRNAs regulate different types of NB-LRR genes, most of them with coil-coiled domains, and trigger the generation of secondary small interfering RNAs (siRNA) as a phased pattern in the target site regions. In addition to acting in response to biotic stresses, they are also responsive to abiotic stresses such as heat, drought, salt, and light stress. Their copy number and expression variation in Triticeae suggest a rapid birth and death frequency. Altogether, non-conserved miRNAs as conserved transcriptional regulators in gymnosperms and angiosperms regulating the disease resistance genes displayed quick plasticity including the variations of sequences, gene copy number, functions, and expression level, which accompanied with NB-LRR genes may be tune-regulated to plants in natural environments with various biotic and abiotic stresses.

Exploring miRNAs for developing climate-resilient crops: A perspective review

Climate changes and environmental stresses have significant implications on global crop production and necessitate developing crops that can withstand an array of climate changes and environmental perturbations such as irregular water-supplies leading to drought or water-logging, hyper soil-salinity, extreme and variable temperatures, ultraviolet radiations and metal stress. Plants have intricate molecular mechanisms to cope with these dynamic environmental changes, one of the most common and effective being the reprogramming of expression of stress-responsive genes. Plant microRNAs (miRNAs) have emerged as key post-transcriptional and translational regulators of gene-expression for modulation of stress implications. Recent reports are establishing their key roles in epigenetic regulations of stress/adaptive responses as well as in providing plants genome-stability. Several stress responsive miRNAs are being identified from different crop plants and miRNA-driven RNA-interference (RNAi) is turning into a technology of choice for improving crop traits and providing phenotypic plasticity in challenging environments. Here we presents a perspective review on exploration of miRNAs as potent targets for engineering crops that can withstand multi-stress environments via loss-/gain-of-function approaches. This review also shed a light on potential roles plant miRNAs play in genome-stability and their emergence as potent target for genome-editing. Current knowledge on plant miRNAs, their biogenesis, function, their targets, and latest developments in bioinformatics approaches for plant miRNAs are discussed. Though there are recent reviews discussing primarily the individual miRNAs responsive to single stress factors, however, considering practical limitation of this approach, special emphasis is given in this review on miRNAs involved in responses and adaptation of plants to multi-stress environments including at epigenetic and/or epigenomic levels.

Deep sequencing identified potential miRNAs involved in defence response, stress and plant growth characteristics of wild genotypes of cardamom

Cardamom has long been used as a food flavouring agent and in ayurvedic medicines for mouth ulcers, digestive problems and even depression. Extensive occurrence of pests and diseases adversely affect its cultivation and result in substantial reductions in total production and productivity. Numerous studies revealed the significant role of miRNAs in plant biotic stress responses. In the current study, miRNA profiling of cultivar and wild cardamom genotypes was performed using an Ion Proton sequencer. We identified 161 potential miRNAs representing 42 families, including monocot/tissue-specific and 14 novel miRNAs in both genotypes. Significant differences in miRNA family abundance between the libraries were observed in read frequencies. A total of 19 miRNAs (from known miRNAs) displayed a twofold difference in expression between wild and cultivar genotypes. We found 1168 unique potential targets for 40 known miRNA families in wild and 1025 potential targets for 42 known miRNA families in cultivar genotypes. The differential expression analysis revealed that most miRNAs identified were highly expressed in cultivars and, furthermore, lower expression of miR169 and higher expression of miR529 in wild cardamom proved evidence that wild genotypes have stronger drought stress tolerance and floral development than cultivars. Potential targets predicted for the newly identified miRNAs from the miRNA libraries of wild and cultivar cardamom genotypes involved in metabolic and developmental processes and in response to various stimuli. qRT-PCR confirmed miRNAs were differentially expressed between wild and cultivar genotypes. Furthermore, four target genes were validated experimentally to confirm miRNA-mRNA target pairing using RNA ligase-mediated 5' Rapid Amplification of cDNA Ends (5'RLM-RACE) PCR.

miRNomes involved in imparting thermotolerance to crop plants

Thermal stress is one of the challenges to crop plants that negatively impacts crop yield. To overcome this ever-growing problem, utilization of regulatory mechanisms, especially microRNAs (miRNAs), that provide efficient and precise regulation in a targeted manner have been found to play determining roles. Besides their roles in plant growth and development, many recent studies have shown differential regulation of several miRNAs during abiotic stresses including heat stress (HS). Thus, understanding the underlying mechanism of miRNA-mediated gene expression during HS will enable researchers to exploit this regulatory mechanism to address HS responses. This review focuses on the miRNAs and regulatory networks that were involved in physiological, metabolic and morphological adaptations during HS in plant, specifically in crops. Illustrated examples including, the miR156-SPL, miR169-NF-YA5, miR395-APS/AST, miR396-WRKY, etc., have been discussed in specific relation to the crop plants. Further, we have also discussed the available plant miRNA databases and bioinformatics tools useful for miRNA identification and study of their regulatory role in response to HS. Finally, we have briefly discussed the future prospects about the miRNA-related mechanisms of HS for improving thermotolerance in crop plants.

Ascophyllum nodosum extract mitigates salinity stress in Arabidopsis thaliana by modulating the expression of miRNA involved in stress tolerance and nutrient acquisition

Ascophyllum nodosum extract (ANE) contains bioactive compounds that improve the growth of Arabidopsis in experimentally-induced saline conditions; however, the molecular mechanisms through which ANE elicits tolerance to salinity remain largely unexplored. Micro RNAs (miRNAs) are key regulators of gene expression, playing crucial roles in plant growth, development, and stress tolerance. Next generation sequencing of miRNAs from leaves of control Arabidopsis and from plants subjected to three treatments (ANE, NaCl and ANE+NaCl) was used to identify ANE-responsive miRNA in the absence and presence of saline conditions. Differential gene expression analysis revealed that ANE had a strong effect on miRNAs expression in both conditions. In the presence of salinity, ANE tended to reduce the up-regulation or the down-regulation trend induced caused by NaCl in miRNAs such as ath-miR396a-5p, ath-miR399, ath-miR2111b and ath-miR827. To further uncover the effects of ANE, the expression of several target genes of a number of ANE-responsive miRNAs was analyzed by qPCR. NaCl, but not ANE, down-regulated miR396a-5p, which negatively regulated the expression of AtGRF7 leading to a higher expression of AtDREB2a and AtRD29 in the presence of ANE+NaCl, as compared to ANE alone. ANE+NaCl initially reduced and then enhanced the expression of ath-miR169g-5p, while the expression of the target genes AtNFYA1 and ATNFYA2, known to be involved in the salinity tolerance mechanism, was increased as compared to ANE or to NaCl treatments. ANE and ANE+NaCl modified the expression of ath-miR399, ath-miR827, ath-miR2111b, and their target genes AtUBC24, AtWAK2, AtSYG1 and At3g27150, suggesting a role of ANE in phosphate homeostasis. In vivo and in vitro experiments confirmed the improved growth of Arabidopsis in presence of ANE, in saline conditions and in phosphate-deprived medium, further substantiating the influence of ANE on a variety of essential physiological processes in Arabidopsis including salinity tolerance and phosphate uptake.

Plant reference genes for development and stress response studies

Many reference genes are used by different laboratories for gene expression analyses to indicate the relative amount of input RNA/DNA in the experiment. These reference genes are supposed to show least variation among the treatments and with the control sets in a given experiment. However, expression of reference genes varies significantly from one set of experiment to the other. Thus, selection of reference genes depends on the experimental conditions. Sometimes the average expression of two or three reference genes is taken as standard. This review consolidated the details of about 120 genes attempted for normalization during comparative expression analysis in 16 different plants. Plant species included in this review are Arabidopsis thaliana, cotton (Gossypium hirsutum), tobacco (Nicotiana benthamiana and N. tabacum), soybean (Glycine max), rice (Oryza sativa), blueberry (Vaccinium corymbosum), tomato (Solanum lycopersicum), wheat (Triticum aestivum), potato (Solanum tuberosum), sugar cane (Saccharum sp.), carrot (Daucus carota), coffee (Coffea arabica), cucumber (Cucumis sativus), kiwi (Actinidia deliciosa) and grape (Vitis vinifera). The list includes model and cultivated crop plants from both monocot and dicot classes. We have categorized plant-wise the reference genes that have been used for expression analyses in any or all of the four different conditions such as biotic stress, abiotic stress, developmental stages and various organs and tissues, reported till date. This review serves as a guide during the reference gene hunt for gene expression analysis studies.

Hybrid incompatibilities in interspecific crosses between tetraploid wheat and its wild diploid relative Aegilops umbellulata

Hybrid abnormalities, severe growth abortion and grass-clump dwarfism, were found in the tetraploid wheat/Aegilops umbellulata hybrids, and the gene expression changes were conserved in the hybrids with those in other wheat synthetic hexaploids. Aegilops umbellulata Zhuk., a diploid goatgrass species with a UU genome, has been utilized as a genetic resource for wheat breeding. Here, we examine the reproductive barriers between tetraploid wheat cultivar Langdon (Ldn) and various Ae. umbellulata accessions by conducting interspecific crossings. Through systematic cross experiments, three types of hybrid incompatibilities were found: seed production failure in crosses, hybrid growth abnormalities and sterility in the ABU hybrids. Hybrid incompatibilities were widely distributed over the entire range of the natural species, and in about 50% of the cross combinations between tetraploid Ldn and Ae. umbellulata accessions, ABU F₁ hybrids showed one of two abnormal growth phenotypes: severe growth abortion (SGA) or grass-clump dwarfism. Expression of the shoot meristem maintenance-related and cell cycle-related genes was markedly repressed in crown tissues of hybrids showing SGA, suggesting dysfunction of mitotic cell division in the shoot apices. The grass-clump dwarf phenotype may be explained by down-regulation of wheat APETALA1-like MADS box genes, which act as flowering promoters, and altered expression in crown tissues of the miR156/SPLs module, which controls tiller number and branching. These gene expression changes in growth abnormalities were well conserved between the Ldn/Ae. umbellulata plants and interspecific hybrids from crosses of Ldn and wheat D-genome progenitor Ae. tauschii.

Identification of drought-responsive miRNAs and physiological characterization of tea plant (Camellia sinensis L.) under drought stress

BACKGROUND

Drought stress is one of the major natural challenges in the main tea-producing regions of China. The tea plant (Camellia sinensis) is a traditional beverage plant whose growth status directly affects tea quality. Recent studies have revealed that microRNAs (miRNAs) play key functions in plant growth and development. Although some miRNAs have been identified in C. sinensis, little is known about their roles in the drought stress response of tea plants.

RESULTS

Physiological characterization of Camellia sinensis 'Tieguanyin' under drought stress showed that the malondialdehyde concentration and electrical conductivity of leaves of drought-stressed plants increased when the chlorophyll concentration decreased under severe drought stress. We sequenced four small-RNA (sRNA) libraries constructed from leaves of plants subjected to four different treatments, normal water supply (CK); mild drought stress (T1); moderate drought stress (T2) and severe drought stress (T3). A total of 299 known mature miRNA sequences and 46 novel miRNAs were identified. Gene Ontology enrichment analysis revealed that most of the differentially expressed-miRNA target genes were related to regulation of transcription. Kyoto Encyclopedia of Genes and Genomes analysis revealed that the most highly enriched pathways under drought stress were D-alanine metabolism, sulfur metabolism, and mineral absorption pathways. Real-time quantitative PCR (qPCR) was used to validate the expression patterns of 21 miRNAs (2 up-regulated and 19 down-regulated under drought stress). The observed co-regulation of the miR166 family and their targets ATHB-14-like and ATHB-15-like indicate the presence of negative feedback regulation in miRNA pathways.

CONCLUSIONS

Analyses of drought-responsive miRNAs in tea plants showed that most of differentially expressed-miRNA target genes were related to regulation of transcription. The results of study revealed that the expressions of phase-specific miRNAs vary with morphological, physiological, and biochemical changes. These findings will be useful for research on drought resistance and provide insights into the mechanisms of drought adaptation and resistance in C. sinensis.

Novel Insights into miRNA Regulation of Storage Protein Biosynthesis during Wheat Caryopsis Development under Drought Stress

Drought stress is a significant abiotic stress factor that affects wheat yield and quality. MicroRNA (miRNA) plays an important role in regulating caryopsis development in response to drought stress. However, little is known about the expression characteristics of miRNAs and how they regulate protein accumulation in wheat caryopsis under drought stress. To address this, two small RNA libraries of wheat caryopsis under control and drought stress conditions were constructed and sequenced. A total of 125 miRNAs were identified in the two samples, of which 110 were known and 15 were novel. A total of 1,981 miRNA target genes were predicted and functional annotations were obtained from various databases for 1,641 of them. Four miRNAs were identified as differential expression under drought stress, and the expression patterns of three of them were consistent with results obtained by reverse transcription polymerase chain reaction (RT-PCR) and reverse transcription quantitative polymerase chain reaction (RT-qPCR). Moreover, three miRNA-target pairs showed negative regulation tendency, as revealed by RT-qPCR. Functional enrichment and pathway analysis revealed that four pathways might be involved in storage protein biosynthesis. Furthermore, drought stress significantly increased the accumulation of protein bodies and protein content in wheat endosperm. In summary, our findings suggest that drought stress may enhance storage protein by regulating the expression of miRNAs and their target genes.

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.

Global gene expression profiling related to temperature-sensitive growth abnormalities in interspecific crosses between tetraploid wheat and Aegilops tauschii

Triploid wheat hybrids between tetraploid wheat and Aegilops tauschii sometimes show abnormal growth phenotypes, and the growth abnormalities inhibit generation of wheat synthetic hexaploids. In type II necrosis, one of the growth abnormalities, necrotic cell death accompanied by marked growth repression occurs only under low temperature conditions. At normal temperature, the type II necrosis lines show grass-clump dwarfism with no necrotic symptoms, excess tillers, severe dwarfism and delayed flowering. Here, we report comparative expression analyses to elucidate the molecular mechanisms of the temperature-dependent phenotypic plasticity in the triploid wheat hybrids. We compared gene and small RNA expression profiles in crown tissues to characterize the temperature-dependent phenotypic plasticity. No up-regulation of defense-related genes was observed under the normal temperature, and down-regulation of wheat APETALA1-like MADS-box genes, considered to act as flowering promoters, was found in the grass-clump dwarf lines. Some microRNAs, including miR156, were up-regulated, whereas the levels of transcripts of the miR156 target genes SPLs, known to inhibit tiller and branch number, were reduced in crown tissues of the grass-clump dwarf lines at the normal temperature. Unusual expression of the miR156/SPLs module could explain the grass-clump dwarf phenotype. Dramatic alteration of gene expression profiles, including miRNA levels, in crown tissues is associated with the temperature-dependent phenotypic plasticity in type II necrosis/grass-clump dwarf wheat hybrids.

Long noncoding miRNA gene represses wheat β-diketone waxes

The cuticle of terrestrial plants functions as a protective barrier against many biotic and abiotic stresses. In wheat and other Triticeae, β-diketone waxes are major components of the epicuticular layer leading to the bluish-white glaucous trait in reproductive-age plants. Glaucousness in durum wheat is controlled by a metabolic gene cluster at the WAX1 (W1) locus and a dominant suppressor INHIBITOR of WAX1 (Iw1) on chromosome 2B. The wheat D subgenome from progenitor Aegilops tauschii contains W2 and Iw2 paralogs on chromosome 2D. Here we identify the Iw1 gene from durum wheat and demonstrate the unique regulatory mechanism by which Iw1 acts to suppress a carboxylesterase-like protein gene, W1-COE, within the W1 multigene locus. Iw1 is a long noncoding RNA (lncRNA) containing an inverted repeat (IR) with >80% identity to W1-COE The Iw1 transcript forms a miRNA precursor-like long hairpin producing a 21-nt predominant miRNA, miRW1, and smaller numbers of related sRNAs associated with the nonglaucous phenotype. When Iw1 was introduced into glaucous bread wheat, miRW1 accumulated, W1-COE and its paralog W2-COE were down-regulated, and the phenotype was nonglaucous and β-diketone-depleted. The IR region of Iw1 has >94% identity to an IR region on chromosome 2 in Ae. tauschii that also produces miRW1 and lies within the marker-based location of Iw2 We propose the Iw loci arose from an inverted duplication of W1-COE and/or W2-COE in ancestral wheat to form evolutionarily young miRNA genes that act to repress the glaucous trait.

Uptake and impact of natural diet-derived small RNA in invertebrates: Implications for ecology and agriculture

The putative transfer and gene regulatory activities of diet-derived small RNAs (sRNAs) in ingesting animals are still debated. The existence of natural uptake of diet-derived sRNA by invertebrate species could have significant implication for our understanding of ecological relationships and could synergize with efforts to use RNA interference (RNAi) technology in agriculture. Here, we synthesize information gathered from studies in invertebrates using natural or artificial dietary delivery of sRNA and from studies of sRNA in vertebrate animals and plants to review our current understanding of uptake and impact of natural diet-derived sRNA on invertebrates. Our understanding has been influenced and sometimes confounded by the diversity of invertebrates and ingested plants studied, our limited insights into how gene expression may be modulated by dietary sRNAs at the mechanistic level, and the paucity of studies focusing directly on natural uptake of sRNA. As such, we suggest 2 strategies to investigate this phenomenon more comprehensively and thus facilitate the realization of its potentially broad impact on ecology and agriculture in the future.

Genome-wide changes in microRNA expression during short and prolonged heat stress and recovery in contrasting rice cultivars

MicroRNAs (miRNAs) are known to regulate expression of genes under stress. We report here the deep sequencing of small RNAs expressed during control, short and prolonged heat stress and recovery. Genome-wide identification of miRNAs in tolerant (Nagina 22) and susceptible (Vandana) rice cultivars was performed in 16 samples representing root and shoot of 13-day-old seedlings. The expression profile of miRNAs was analysed in 36 pairwise combinations to identify the genotype-, treatment- and tissue-dependent expression of miRNAs. Small-RNA sequencing of 16 libraries yielded ~271 million high-quality raw sequences; 162 miRNA families were identified. The highly expressed miRNAs in rice tissues were miR166, miR168, miR1425, miR529, mR162, miR1876, and miR1862. Expression of osa-miR1436, osa-miR5076, osa-miR5161, and osa-miR6253 was observed only in stressed tissue of both genotypes indicating their general role in heat stress response. Expression of osa-miR1439, osa-miR1848, osa-miR2096, osa-miR2106, osa-miR2875, osa-miR3981, osa-miR5079, osa-miR5151, osa-miR5484, osa-miR5792, and osa-miR5812 was observed only in Nagina 22 during high temperature, suggesting a specific role of these miRNAs in heat stress tolerance. This study provides details of the repertoire of miRNAs expressed in root and shoot of heat susceptible and tolerant rice genotypes under heat stress and recovery.

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.

SMARTER De-Stressed Cereal Breeding

In cereal breeding programs, improved yield potential and stability are ultimate goals when developing new varieties. To facilitate achieving these goals, reproductive success under stressful growing conditions is of the highest priority. In recent times, small RNA (sRNA)-mediated pathways have been associated with the regulation of genes involved in stress adaptation and reproduction in both model plants and several cereals. Reproductive and physiological traits such as flowering time, reproductive branching, and root architecture can be manipulated by sRNA regulatory modules. We review sRNA-mediated pathways that could be exploited to expand crop diversity with adaptive traits and, in particular, the development of high-yielding stress-tolerant cereals: SMARTER cereal breeding through 'Small RNA-Mediated Adaptation of Reproductive Targets in Epigenetic Regulation'.

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.

PmiRExAt: plant miRNA expression atlas database and web applications

High-throughput small RNA (sRNA) sequencing technology enables an entirely new perspective for plant microRNA (miRNA) research and has immense potential to unravel regulatory networks. Novel insights gained through data mining in publically available rich resource of sRNA data will help in designing biotechnology-based approaches for crop improvement to enhance plant yield and nutritional value. Bioinformatics resources enabling meta-analysis of miRNA expression across multiple plant species are still evolving. Here, we report PmiRExAt, a new online database resource that caters plant miRNA expression atlas. The web-based repository comprises of miRNA expression profile and query tool for 1859 wheat, 2330 rice and 283 maize miRNA. The database interface offers open and easy access to miRNA expression profile and helps in identifying tissue preferential, differential and constitutively expressing miRNAs. A feature enabling expression study of conserved miRNA across multiple species is also implemented. Custom expression analysis feature enables expression analysis of novel miRNA in total 117 datasets. New sRNA dataset can also be uploaded for analysing miRNA expression profiles for 73 plant species. PmiRExAt application program interface, a simple object access protocol web service allows other programmers to remotely invoke the methods written for doing programmatic search operations on PmiRExAt database.Database URL:http://pmirexat.nabi.res.in.

Identification of Arbuscular Mycorrhiza (AM)-Responsive microRNAs in Tomato

A majority of land plants can form symbiosis with arbuscular mycorrhizal (AM) fungi. MicroRNAs (miRNAs) have been implicated to regulate this process in legumes, but their involvement in non-legume species is largely unknown. In this study, by performing deep sequencing of sRNA libraries in tomato roots and comparing with tomato genome, a total of 700 potential miRNAs were predicted, among them, 187 are known plant miRNAs that have been previously deposited in miRBase. Unlike the profiles in other plants such as rice and Arabidopsis, a large proportion of predicted tomato miRNAs was 24 nt in length. A similar pattern was observed in the potato genome but not in tobacco, indicating a Solanum genus-specific expansion of 24-nt miRNAs. About 40% identified tomato miRNAs showed significantly altered expressions upon Rhizophagus irregularis inoculation, suggesting the potential roles of these novel miRNAs in AM symbiosis. The differential expression of five known and six novel miRNAs were further validated using qPCR analysis. Interestingly, three up-regulated known tomato miRNAs belong to a known miR171 family, a member of which has been reported in Medicago truncatula to regulate AM symbiosis. Thus, the miR171 family likely regulates AM symbiosis conservatively across different plant lineages. More than 1000 genes targeted by potential AM-responsive miRNAs were provided and their roles in AM symbiosis are worth further exploring.

Cross-talk between freezing response and signaling for regulatory transcriptions of MIR475b and its targets by miR475b promoter in Populus suaveolens

MicroRNAs (miRNAs) are small, non-coding RNAs that play important roles in post-transcriptional regulation of their target genes, yet the transcriptional regulation of plant miRNAs by promoter is poorly understood. Here, we firstly clone pri-miR475b cDNA and its native promoter from P. suaveolens, and characterize Psu-MIR475b as class-II gene transcribed by RNA polymerase II. By 5′ deletion analysis of Psu-miR475b promoter in a series of promoter-GUS chimeric vectors, we functionally identify three positive regulatory regions and multiple cis-acting elements responsible for Psu-miR475b promoter activity in response to freezing stress and exogenous hormone treatment. Moreover, the Psu-miR475b promoter activity displays a tissue-specific manner, negatively regulated by freezing stress and positively by MeJA, SA or GA treatment. Importantly, we comparatively analyze the time-course transcriptional profiles of Psu-miR475b and its targets in Psu-miR475b over-expression transgenic plants controlled by Psu-miR475b-specific promoter or CaMV 35S constitutive promoter, and explore the regulatory mechanism of Psu-miR475b promoter controlling transcriptional expressions of Psu-MIR475b and its targets in response to freezing stress and exogenous hormone treatment. Our results reveal that Psu-miR475b promoter-mediated transcriptions of Psu-MIR475b and its targets in response to freezing stress may be involved in a cross-talk between freezing response and stress signaling process.

MicroRNAs As Potential Targets for Abiotic Stress Tolerance in Plants

The microRNAs (miRNAs) are small (20-24 nt) sized, non-coding, single stranded riboregulator RNAs abundant in higher organisms. Recent findings have established that plants assign miRNAs as critical post-transcriptional regulators of gene expression in sequence-specific manner to respond to numerous abiotic stresses they face during their growth cycle. These small RNAs regulate gene expression via translational inhibition. Usually, stress induced miRNAs downregulate their target mRNAs, whereas, their downregulation leads to accumulation and function of positive regulators. In the past decade, investigations were mainly aimed to identify plant miRNAs, responsive to individual or multiple environmental factors, profiling their expression patterns and recognizing their roles in stress responses and tolerance. Altered expressions of miRNAs implicated in plant growth and development have been reported in several plant species subjected to abiotic stress conditions such as drought, salinity, extreme temperatures, nutrient deprivation, and heavy metals. These findings indicate that miRNAs may hold the key as potential targets for genetic manipulations to engineer abiotic stress tolerance in crop plants. This review is aimed to provide recent updates on plant miRNAs, their biogenesis and functions, target prediction and identification, computational tools and databases available for plant miRNAs, and their roles in abiotic stress-responses and adaptive mechanisms in major crop plants. Besides, the recent case studies for overexpressing the selected miRNAs for miRNA-mediated enhanced abiotic stress tolerance of transgenic plants have been discussed.

PmiRExAt: plant miRNA expression atlas database and web applications

High-throughput small RNA (sRNA) sequencing technology enables an entirely new perspective for plant microRNA (miRNA) research and has immense potential to unravel regulatory networks. Novel insights gained through data mining in publically available rich resource of sRNA data will help in designing biotechnology-based approaches for crop improvement to enhance plant yield and nutritional value. Bioinformatics resources enabling meta-analysis of miRNA expression across multiple plant species are still evolving. Here, we report PmiRExAt, a new online database resource that caters plant miRNA expression atlas. The web-based repository comprises of miRNA expression profile and query tool for 1859 wheat, 2330 rice and 283 maize miRNA. The database interface offers open and easy access to miRNA expression profile and helps in identifying tissue preferential, differential and constitutively expressing miRNAs. A feature enabling expression study of conserved miRNA across multiple species is also implemented. Custom expression analysis feature enables expression analysis of novel miRNA in total 117 datasets. New sRNA dataset can also be uploaded for analysing miRNA expression profiles for 73 plant species. PmiRExAt application program interface, a simple object access protocol web service allows other programmers to remotely invoke the methods written for doing programmatic search operations on PmiRExAt database.Database URL:http://pmirexat.nabi.res.in.

Genome-wide identification of conserved microRNA and their response to drought stress in Dongxiang wild rice (Oryza rufipogon Griff.)

OBJECTIVES

To identify drought stress-responsive conserved microRNA (miRNA) from Dongxiang wild rice (Oryza rufipogon Griff., DXWR) on a genome-wide scale, high-throughput sequencing technology was used to sequence libraries of DXWR samples, treated with and without drought stress.

RESULTS

505 conserved miRNAs corresponding to 215 families were identified. 17 were significantly down-regulated and 16 were up-regulated under drought stress. Stem-loop qRT-PCR revealed the same expression patterns as high-throughput sequencing, suggesting the accuracy of the sequencing result was high. Potential target genes of the drought-responsive miRNA were predicted to be involved in diverse biological processes. Furthermore, 16 miRNA families were first identified to be involved in drought stress response from plants.

CONCLUSION

These results present a comprehensive view of the conserved miRNA and their expression patterns under drought stress for DXWR, which will provide valuable information and sequence resources for future basis studies.

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

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

Constitutive Expression of Rice MicroRNA528 Alters Plant Development and Enhances Tolerance to Salinity Stress and Nitrogen Starvation in Creeping Bentgrass

MicroRNA528 (miR528) is a conserved monocot-specific small RNA that has the potential of mediating multiple stress responses. So far, however, experimental functional studies of miR528 are lacking. Here, we report that overexpression of a rice (Oryza sativa) miR528 (Osa-miR528) in transgenic creeping bentgrass (Agrostis stolonifera) alters plant development and improves plant salt stress and nitrogen (N) deficiency tolerance. Morphologically, miR528-overexpressing transgenic plants display shortened internodes, increased tiller number, and upright growth. Improved salt stress resistance is associated with increased water retention, cell membrane integrity, chlorophyll content, capacity for maintaining potassium homeostasis, CATALASE activity, and reduced ASCORBIC ACID OXIDASE (AAO) activity; while enhanced tolerance to N deficiency is associated with increased biomass, total N accumulation and chlorophyll synthesis, nitrite reductase activity, and reduced AAO activity. In addition, AsAAO and COPPER ION BINDING PROTEIN1 are identified as two putative targets of miR528 in creeping bentgrass. Both of them respond to salinity and N starvation and are significantly down-regulated in miR528-overexpressing transgenics. Our data establish a key role that miR528 plays in modulating plant growth and development and in the plant response to salinity and N deficiency and indicate the potential of manipulating miR528 in improving plant abiotic stress resistance.

miRNA-based drought regulation in wheat

MicroRNAs (miRNAs) are a class of small non-coding regulatory RNAs that regulate gene expression by guiding target mRNA cleavage or translational inhibition. Drought is a common environmental stress influencing crop growth and development. To date, it has been reported that a number of plant miRNA are involved in drought stress response. In this study, we comparatively investigated drought stress-responsive miRNAs in the root and leaf of bread wheat (Triticum aestivum cv. Sivas 111/33) by miRNA microarray screening. miRNA microarray analysis showed that 285 miRNAs (207 upregulated and 78 downregulated) and 244 miRNAs (115 upregulated and 129 downregulated) were differentially expressed in leaf and root tissues, respectively. Among the differentially expressed miRNAs, 23 miRNAs were only expressed in the leaf and 26 miRNAs were only expressed in the root of wheat growth under drought stress. Upon drought treatment, expression of miR159, miR160, miR166, miR169, miR172, miR395, miR396, miR408, miR472, miR477, miR482, miR1858, miR2118, and miR5049 were found to be significantly differentiated in bread wheat. The regulatory network analysis showed that miR395 has connections with a number of target transcripts, and miR159 and miR319 share a number of target genes. Drought-tolerant and drought-sensitive wheat cultivars showed altered expression pattern upon drought stress in terms of investigated miRNA and their target transcript expression level.