Identification of four functionally important microRNA families with contrasting differential expression profiles between drought-tolerant and susceptible rice leaf at vegetative stage

Comparative miRNome and transcriptome analyses reveal the expression of novel miRNAs in the panicle of rice implicated in sustained agronomic performance under terminal drought stress

MAIN CONCLUSION

Differential expression of 128 known and 111 novel miRNAs in the panicle of Nagina 22 under terminal drought stress targeting transcription factors, stress-associated genes, etc., enhances drought tolerance and helps sustain agronomic performance under terminal drought stress. Drought tolerance is a complex multigenic trait, wherein the genes are fine-tuned by coding and non-coding components in mitigating deleterious effects. MicroRNA (miRNA) controls gene expression at post-transcriptional level either by cleaving mRNA (transcript) or by suppressing its translation. miRNAs are known to control developmental processes and abiotic stress tolerance in plants. To identify terminal drought-responsive novel miRNA in contrasting rice cultivars, we constructed small RNA (sRNA) libraries from immature panicles of drought-tolerant rice [Nagina 22 (N 22)] and drought-sensitive (IR 64) cultivars grown under control and terminal drought stress. Our analysis of sRNA-seq data resulted in the identification of 169 known and 148 novel miRNAs in the rice cultivars. Among the novel miRNAs, 68 were up-regulated while 43 were down-regulated in the panicle of N 22 under stress. Interestingly, 31 novel miRNAs up-regulated in N 22 were down-regulated in IR 64, whereas 4 miRNAs down-regulated in N 22 were up-regulated in IR 64 under stress. To detect the effects of miRNA on mRNA expression level, transcriptome analysis was performed, while differential expression of miRNAs and their target genes was validated by RT-qPCR. Targets of the differentially expressed miRNAs include transcription factors and stress-associated genes involved in cellular/metabolic/developmental processes, response to abiotic stress, programmed cell death, photosynthesis, panicle/seed development, and grain yield. Differential expression of the miRNAs could be validated in an independent set of the samples. The findings might be useful in genetic improvement of drought-tolerant rice.

Defense Responses of Different Rice Varieties Affect Growth Performance and Food Utilization of Cnaphalocrocis medinalis Larvae

Rice leaf folder, Cnaphalocrocis medinalis (Guenée), is one of the most serious pests on rice. At present, chemical control is the main method for controlling this pest. However, the indiscriminate use of chemical insecticides has non-target effects and may cause environmental pollution. Besides, leaf curling behavior by C. medinalis may indirectly reduce the efficacy of chemical spray. Therefore, it is crucial to cultivate efficient rice varieties resistant to this pest. Previous studies have found that three different rice varieties, Zhongzao39 (ZZ39), Xiushui134 (XS134), and Yongyou1540 (YY1540), had varying degrees of infestation by C. medinalis. However, it is currently unclear whether the reason for this difference is related to the difference in defense ability of the three rice varieties against the infestation of C. medinalis. To explore this issue, the current study investigated the effects of three rice varieties on the growth performance and food utilization capability of the 4th instar C. medinalis. Further, it elucidated the differences in defense responses among different rice varieties based on the differences in leaf physiological and biochemical indicators and their impact on population occurrence. The results showed that the larval survival rate was the lowest, and the development period was significantly prolonged after feeding on YY1540. This was not related to the differences in leaf wax, pigments, and nutritional components among the three rice varieties nor to the feeding preferences of the larvae. The rate of superoxide anion production, hydrogen peroxide content, and the activity of three protective enzymes were negatively correlated with larval survival rate, and they all showed the highest in YY1540 leaves. Compared to other tested varieties, although the larvae feeding on YY1540 had higher conversion efficiency of ingested food and lower relative consumption rate, their relative growth was faster, indicating stronger food utilization capability. However, they had a lower accumulation of protein. This suggests that different rice varieties had different levels of oxidative stress after infestation by C. medinalis. The defense response of YY1540 was more intense, which was not conducive to the development of the larvae population. These results will provide new insights into the interaction mechanism between different rice varieties and C. medinalis and provide a theoretical basis for cultivating rice varieties resistant to this pest.

Genome-wide identification, characterization, and expression analysis unveil the roles of pseudouridine synthase (PUS) family proteins in rice development and stress response

Pseudouridylation, the conversion of uridine (U) to pseudouridine (Ѱ), is one of the most prevalent and evolutionary conserved RNA modifications, which is catalyzed by pseudouridine synthase (PUS) enzymes. Ѱs play a crucial epitranscriptomic role by regulating attributes of cellular RNAs across diverse organisms. However, the precise biological functions of PUSs in plants remain largely elusive. In this study, we identified and characterized 21 members in the rice PUS family which were categorized into six distinct subfamilies, with RluA and TruA emerging as the most extensive. A comprehensive analysis of domain structures, motifs, and homology modeling revealed that OsPUSs possess all canonical features of true PUS proteins, essential for substrate recognition and catalysis. The exploration of OsPUS promoters revealed presence of cis-acting regulatory elements associated with hormone and abiotic stress responses. Expression analysis of OsPUS genes showed differential expression at developmental stages and under stress conditions. Notably, OsTruB3 displayed high expression in salt, heat, and drought stresses. Several OsRluA members showed induction in heat stress, while a significant decline in expression was observed for various OsTruA members in drought and salinity. Furthermore, miRNAs predicted to target OsPUSs were themselves responsive to variable stresses, adding an additional layer of regulatory complexity of OsPUSs. Study of protein-protein interaction networks provided substantial support for the potential regulatory role of OsPUSs in numerous cellular and stress response pathways. Conclusively, our study provides functional insights into the OsPUS family, contributing to a better understanding of their crucial roles in shaping the development and stress adaptation in rice.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01396-4.

Function identification of miR159a, a positive regulator during poplar resistance to drought stress

Drought seriously affects the growth and development of plants. MiR159 is a highly conserved and abundant microRNA family that plays a crucial role in plant growth and stress responses. However, studies of its function in woody plants are still lacking. Here, the expression of miR159a was significantly upregulated after drought treatment in poplar, and the overexpression of miR159a (OX159a) significantly reduced the open area of the stomata and improved water-use efficiency in poplar. After drought treatment, OX159a lines had better scavenging ability of reactive oxygen species and damage of the membrane system was less than that in wild-type lines. MYB was the target gene of miR159a, as verified by psRNATarget prediction, RT-qPCR, degradome sequencing, and 5' rapid amplification of cDNA ends (5' RACE). Additionally, miR159a-short tandem target mimic suppression (STTM) poplar lines showed increased sensitivity to drought stress. Transcriptomic analysis comparing OX159a lines with wild-type lines revealed upregulation of a series of genes related to response to water deprivation and metabolite synthesis. Moreover, drought-responsive miR172d and miR398 were significantly upregulated and downregulated respectively in OX159a lines. This investigation demonstrated that miR159a played a key role in the tolerance of poplar to drought by reducing stomata open area, increasing the number and total area of xylem vessels, and enhancing water-use efficiency, and provided new insights into the role of plant miR159a and crucial candidate genes for the molecular breeding of trees with tolerance to drought stress.

Integration of miRNA dynamics and drought tolerant QTLs in rice reveals the role of miR2919 in drought stress response

To combat drought stress in rice, a major threat to global food security, three major quantitative trait loci for 'yield under drought stress' (qDTYs) were successfully exploited in the last decade. However, their molecular basis still remains unknown. To understand the role of secondary regulation by miRNA in drought stress response and their relation, if any, with the three qDTYs, the miRNA dynamics under drought stress was studied at booting stage in two drought tolerant (Sahbaghi Dhan and Vandana) and one drought sensitive (IR 20) cultivars. In total, 53 known and 40 novel differentially expressed (DE) miRNAs were identified. The primary drought responsive miRNAs were Osa-MIR2919, Osa-MIR3979, Osa-MIR159f, Osa-MIR156k, Osa-MIR528, Osa-MIR530, Osa-MIR2091, Osa-MIR531a, Osa-MIR531b as well as three novel ones. Sixty-one target genes that corresponded to 11 known and 4 novel DE miRNAs were found to be co-localized with the three qDTYs, out of the 1746 target genes identified. We could validate miRNA-mRNA expression under drought for nine known and three novel miRNAs in eight different rice genotypes showing varying degree of tolerance. From our study, Osa-MIR2919, Osa-MIR3979, Osa-MIR528, Osa-MIR2091-5p and Chr01_11911S14Astr and their target genes LOC_Os01g72000, LOC_Os01g66890, LOC_Os01g57990, LOC_Os01g56780, LOC_Os01g72834, LOC_Os01g61880 and LOC_Os01g72780 were identified as the most promising candidates for drought tolerance at booting stage. Of these, Osa-MIR2919 with 19 target genes in the qDTYs is being reported for the first time. It acts as a negative regulator of drought stress tolerance by modulating the cytokinin and brassinosteroid signalling pathway.

Unveiling the biosynthesis, mechanisms, and impacts of miRNAs in drought stress resilience in plants

Drought stress is one of the most serious threats to sustainable agriculture and is predicted to be further intensified in the coming decades. Therefore, understanding the mechanism of drought stress tolerance and the development of drought-resilient crops are the major goals at present. In recent years, noncoding microRNAs (miRNAs) have emerged as key regulators of gene expressions under drought stress conditions and are turning out to be the potential candidates that can be targeted to develop drought-resilient crops in the future. miRNAs are known to target and decrease the expression of various genes to govern the drought stress response in plants. In addition, emerging evidence also suggests a regulatory role of long non-coding RNAs (lncRNAs) in the regulation of miRNAs and the expression of their target genes by a process referred as miRNA sponging. In this review, we present the regulatory roles of miRNAs in the modulation of drought-responsive genes along with discussing their biosynthesis and action mechanisms. Additionally, the interactive roles of miRNAs with phytohormone signaling components have also been highlighted to present the global view of miRNA functioning under drought-stress conditions.

miRNAs for crop improvement

Climate change significantly impacts crop production by inducing several abiotic and biotic stresses. The increasing world population, and their food and industrial demands require focused efforts to improve crop plants to ensure sustainable food production. Among various modern biotechnological tools, microRNAs (miRNAs) are one of the fascinating tools available for crop improvement. miRNAs belong to a class of small non-coding RNAs playing crucial roles in numerous biological processes. miRNAs regulate gene expression by post-transcriptional target mRNA degradation or by translation repression. Plant miRNAs have essential roles in plant development and various biotic and abiotic stress tolerance. In this review, we provide propelling evidence from previous studies conducted around miRNAs and provide a one-stop review of progress made for breeding stress-smart future crop plants. Specifically, we provide a summary of reported miRNAs and their target genes for improvement of plant growth and development, and abiotic and biotic stress tolerance. We also highlight miRNA-mediated engineering for crop improvement and sequence-based technologies available for the identification of miRNAs associated with stress tolerance and plant developmental events.

Comparative Genome-Wide Analysis of MicroRNAs and Their Target Genes in Roots of Contrasting Indica Rice Cultivars under Reproductive-Stage Drought

Recurrent occurrence of drought stress in varying intensity has become a common phenomenon in the present era of global climate change, which not only causes severe yield losses but also challenges the cultivation of rice. This raises serious concerns for sustainable food production and global food security. The root of a plant is primarily responsible to perceive drought stress and acquire sufficient water for the survival/optimal growth of the plant under extreme climatic conditions. Earlier studies reported the involvement/important roles of microRNAs (miRNAs) in plants' responses to environmental/abiotic stresses. A number (738) of miRNAs is known to be expressed in different tissues under varying environmental conditions in rice, but our understanding of the role, mode of action, and target genes of the miRNAs are still elusive. Using contrasting rice [IR-64 (reproductive-stage drought sensitive) and N-22 (drought-tolerant)] cultivars, imposed with terminal (reproductive-stage) drought stress, we demonstrate differential expression of 270 known and 91 novel miRNAs in roots of the contrasting rice cultivars in response to the stress. Among the known miRNAs, osamiR812, osamiR166, osamiR156, osamiR167, and osamiR396 were the most differentially expressed miRNAs between the rice cultivars. In the root of N-22, 18 known and 12 novel miRNAs were observed to be exclusively expressed, while only two known (zero novels) miRNAs were exclusively expressed in the roots of IR-64. The majority of the target gene(s) of the miRNAs were drought-responsive transcription factors playing important roles in flower, grain development, auxin signaling, root development, and phytohormone-crosstalk. The novel miRNAs identified in this study may serve as good candidates for the genetic improvement of rice for terminal drought stress towards developing climate-smart rice for sustainable food production.

Role of non-coding RNAs against salinity stress in Oryza species: Strategies and challenges in analyzing miRNAs, tRFs and circRNAs

Salinity is an imbalanced concentration of mineral salts in the soil or water that causes yield loss in salt-sensitive crops. Rice plant is vulnerable to soil salinity stress at seedling and reproductive stages. Different non-coding RNAs (ncRNAs) post-transcriptionally regulate different sets of genes during different developmental stages under varying salinity tolerance levels. While microRNAs (miRNAs) are well known small endogenous ncRNAs, tRNA-derived RNA fragments (tRFs) are an emerging class of small ncRNAs derived from tRNA genes with a demonstrated regulatory role, like miRNAs, in humans but unexplored in plants. Circular RNA (circRNA), another ncRNA produced by back-splicing events, acts as target mimics by preventing miRNAs from binding with their target mRNAs, thereby reducing the miRNA's action upon its target. Same may hold true between circRNAs and tRFs. Hence, the work done on these ncRNAs was reviewed and no reports were found for circRNAs and tRFs under salinity stress in rice, either at seedling or reproductive stages. Even the reports on miRNAs are restricted to seedling stage only, in spite of severe effects on rice crop production due to salt stress during reproductive stage. Moreover, this review sheds light on strategies to predict and analyze these ncRNAs in an effective manner.

Integrative analysis of microRNAs and mRNAs reveals the regulatory networks of triterpenoid saponin metabolism in Soapberry (Sapindus mukorossi Gaertn.)

Triterpenoid saponin are important secondary metabolites and bioactive constituents of soapberry (Sapindus mukorossi Gaertn.) and are widely used in medicine and toiletry products. However, little is known about the roles of miRNAs in the regulation of triterpenoid saponin biosynthesis in soapberry. In this study, a total of 3036 miRNAs were identified, of which 1372 miRNAs were differentially expressed at different stages of pericarp development. Important KEGG pathways, such as terpenoid backbone biosynthesis, sesquiterpenoid and triterpenoid biosynthesis, and basal transcription factors were highlighted, as well the roles of some key miRNAs, such as ath-miR5021, han-miR3630-3p, and ppe-miR858, which may play important roles in regulating triterpenoid saponin biosynthesis. In addition, 58 miRNAs might participate in saponin biosynthesis pathways by predicting the targets of those miRNAs to 53 saponin biosynthesis structural genes. And 75 miRNAs were identified to potentially play vital role in saponin accumulation by targeting transcript factor genes, bHLH, bZIP, ERF, MYB, and WRKY, respectively, which are candidate regulatory genes in the pathway of saponin biosynthesis. The results of weighted gene coexpression network analysis (WGCNA) suggested that two saponin-specific miRNA modules and 10 hub miRNAs may participate in saponin biosynthesis. Furthermore, multiple miRNA-mRNA regulatory networks potentially involved in saponin biosynthesis were generated, e.g., ath-miR5021-SmIDI2/SmGPS5/SmbAS1/SmCYP71D-3/SmUGT74G-2, han-miR3630-3p-SmCYP71A-14/SmbHLH54/SmMYB135/SmWRKY32, and ppe-miR858-SmMYB5/SmMYB32. qRT-PCR analysis validated the expression patterns of nine miRNAs and 12 corresponding target genes. This study represents the first comprehensive analysis of miRNAs in soapberry and lays the foundation for further understanding of miRNA-based regulation in triterpenoid saponin biosynthesis.

Molecular Aspects of MicroRNAs and Phytohormonal Signaling in Response to Drought Stress: A Review

Phytohormones play an essential role in plant growth and development in response to environmental stresses. However, plant hormones require a complex signaling network combined with other signaling pathways to perform their proper functions. Thus, multiple phytohormonal signaling pathways are a prerequisite for understanding plant defense mechanism against stressful conditions. MicroRNAs (miRNAs) are master regulators of eukaryotic gene expression and are also influenced by a wide range of plant development events by suppressing their target genes. In recent decades, the mechanisms of phytohormone biosynthesis, signaling, pathways of miRNA biosynthesis and regulation were profoundly characterized. Recent findings have shown that miRNAs and plant hormones are integrated with the regulation of environmental stress. miRNAs target several components of phytohormone pathways, and plant hormones also regulate the expression of miRNAs or their target genes inversely. In this article, recent developments related to molecular linkages between miRNAs and phytohormones were reviewed, focusing on drought stress.

Comparative transcriptomic analysis and antioxidant defense mechanisms in clusterbean (Cyamopsis tetragonoloba (L.) Taub.) genotypes with contrasting drought tolerance

To understand drought tolerance mechanism(s) in clusterbean (Cyamopsis tetragonoloba), we conducted physiological, biochemical, and de novo comparative transcriptome analysis of drought-tolerant (RGC-1002) and drought-sensitive (RGC-1066) genotypes subjected to 30 days of drought stress. Relative water content (RWC) was maintained in tolerant genotype but was reduced in sensitive genotype. Leaf pigment concentrations were higher in tolerant genotype. Net photosynthesis was significantly decreased in sensitive genotype but insignificant reduction was found in tolerant genotype. Enzymatic antioxidant (GR, APX, DHAR) activities were enhanced in tolerant genotype, while there were insignificant changes in these enzymes in sensitive genotype. The ratios of antioxidant molecules (ASC/DHA and GSH/GSSG) were higher in tolerant genotype as compared to sensitive genotype. In sensitive genotype, 6625 differentially expressed genes (DEGs) were upregulated and 5365 genes were downregulated. In tolerant genotype, 5206 genes were upregulated and 2793 genes were downregulated. In tolerant genotype, transketolase family protein, phosphoenolpyruvate carboxylase 3, temperature-induced lipocalin, and cytochrome oxidase were highly upregulated. Moreover, according to Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the drought tolerance may be attributed to upregulated starch and sucrose metabolism-related genes in tolerant genotype. Finally, quantitative real-time PCR confirmed the reproducibility of the RNA-seq data.

Physiological and Transcriptomic Analyses Reveal Exogenous Trehalose Is Involved in the Responses of Wheat Roots to High Temperature Stress

High temperature stress seriously limits the yield and quality of wheat. Trehalose, a non-reducing disaccharide, has been shown involved in regulating plant responses to a variety of environmental stresses. This study aimed to explore the molecular regulatory network of exogenous trehalose to improve wheat heat tolerance through RNA-sequencing technology and physiological determination. The physiological data and RNA-seq showed that trehalose reduced malondialdehyde content and relative conductivity in wheat roots, and affecting the phenylpropane biosynthesis, starch and sucrose metabolism, glutathione metabolism, and other pathways. Our results showed that exogenous trehalose alleviates the oxidative damage caused by high temperature, coordinating the effect of wheat on heat stress by re-encoding the overall gene expression, but two wheat varieties showed different responses to high temperature stress after trehalose pretreatment. This study preliminarily revealed the effect of trehalose on gene expression regulation of wheat roots under high temperature stress, which provided a reference for the study of trehalose.

miRPreM and tiRPreM: Improved methodologies for the prediction of miRNAs and tRNA-induced small non-coding RNAs for model and non-model organisms

In recent years, microRNAs (miRNAs) and tRNA-derived RNA fragments (tRFs) have been reported extensively following different approaches of identification and analysis. Comprehensively analyzing the present approaches to overcome the existing variations, we developed a benchmarking methodology each for the identification of miRNAs and tRFs, termed as miRNA Prediction Methodology (miRPreM) and tRNA-induced small non-coding RNA Prediction Methodology (tiRPreM), respectively. We emphasized the use of respective genome of organism under study for mapping reads, sample data with at least two biological replicates, normalized read count support and novel miRNA prediction by two standard tools with multiple runs. The performance of these methodologies was evaluated by using Oryza coarctata, a wild rice species as a case study for model and non-model organisms. With organism-specific reference genome approach, 98 miRNAs and 60 tRFs were exclusively found. We observed high accuracy (13 out of 15) when tested these genome-specific miRNAs in support of analyzing the data with respective organism. Such a strong impact of miRPreM, we have predicted more than double number of miRNAs (186) as compared with the traditional approaches (79) and with tiRPreM, we have predicted all known classes of tRFs within the same small RNA data. Moreover, the methodologies presented here are in standard form in order to extend its applicability to different organisms rather than restricting to plants. Hence, miRPreM and tiRPreM can fulfill the need of a comprehensive methodology for miRNA prediction and tRF identification, respectively, for model and non-model organisms.

The sweetpotato β-amylase gene IbBAM1.1 enhances drought and salt stress resistance by regulating ROS homeostasis and osmotic balance

Abiotic stressors, such as drought and high salinity, seriously affect plant growth, productivity, and quality. Maintaining reactive oxygen species (ROS) homeostasis and osmotic balance plays a crucial role in abiotic stress tolerance. β-amylase (BAM) hydrolyzes α-1,4-glycosidic bonds by releasing maltose from starch in the regulation of soluble sugars. However, the function and mechanism of BAMs related to abiotic stress resistance remain unclear in sweetpotato (Ipomoea batatas (L.) Lam.). In this study, we isolated a novel β-amylase gene IbBAM1.1, which was strongly induced by PEG6000, NaCl, and maltose treatments in sweetpotato variety Yanshu25. Overexpression of IbBAM1.1 conferred enhanced tolerance to the drought and high salinity stressors in Arabidopsis thaliana. The activity of β-amylase and the degradation of starch were promoted under drought or salt stress. Accordingly, the contents of osmoprotectants, including maltose and proline were significantly higher in the transgenic lines than those in wild type (WT) plants. Less ROS, such as H₂O₂ and O₂^-, accumulated in the overexpressing lines than in WT plants. Superoxide dismutase activity was strongly enhanced and the level of malondialdehyde was lower under the drought or salt treatment in transgenic plants. Taken together, these results demonstrate that IbBAM1.1 acted as a positive regulator, at least in part, by regulating the level of osmoprotectants to balance the osmotic pressure and activate the scavenging system to maintain ROS homeostasis in the 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.

High-throughput miRNA deep sequencing in response to drought stress in sugarcane

Drought is a major factor which reduces cane growth and productivity. In the present study, we sequenced drought susceptible (V1) and drought tolerant (V2) sugarcane varieties using high-throughput miRNA deep sequencing method to study the regulation of gene expression by miRNAs during drought stress in sugarcane. A total of 1224 conserved miRNAs which belong to 89 miRNA families were identified and 38% of the differentially regulated miRNAs were common for both varieties. Additionally 435 novel miRNAs were also identified from four small RNA libraries. We identified 145 miRNAs that were differentially expressed in susceptible variety (V1-31) and 143 miRNAs differentially expressed in the tolerant variety (V2-31). Target prediction revealed that the genes mainly encoded transcription factors, proteins, phosphatase and kinases involved in signal transduction pathways, integral component of membrane and inorganic ion transport metabolism, enzymes involved in carbohydrate transport and metabolism and drought-stress-related proteins involved in defense mechanisms. Pathway analysis of targets revealed that "General function prediction only" was the most significant pathway observed in both tolerant and susceptible genotypes followed by "signal transduction mechanisms". Functional annotation of the transcripts revealed genes like calcium-dependent protein kinase, respiratory burst oxidase, caffeic acid 3-O-methyltransferase, peroxidase, calmodulin, glutathione S-transferase and transcription factors like MYB, WRKY that are involved in drought tolerant pathways. qRT-PCR was used to verify the expression levels of miRNAs and their potential targets obtained from RNA sequencing results.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-021-02857-x.

SA-Mediated Regulation and Control of Abiotic Stress Tolerance in Rice

Environmental or abiotic stresses are a common threat that remains a constant and common challenge to all plants. These threats whether singular or in combination can have devastating effects on plants. As a semiaquatic plant, rice succumbs to the same threats. Here we systematically look into the involvement of salicylic acid (SA) in the regulation of abiotic stress in rice. Studies have shown that the level of endogenous salicylic acid (SA) is high in rice compared to any other plant species. The reason behind this elevated level and the contribution of this molecule towards abiotic stress management and other underlying mechanisms remains poorly understood in rice. In this review we will address various abiotic stresses that affect the biochemistry and physiology of rice and the role played by SA in its regulation. Further, this review will elucidate the potential mechanisms that control SA-mediated stress tolerance in rice, leading to future prospects and direction for investigation.

Pan-Genome miRNomics in Brachypodium

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

Plant Non-Coding RNAs: Origin, Biogenesis, Mode of Action and Their Roles in Abiotic Stress

As sessile species, plants have to deal with the rapidly changing environment. In response to these environmental conditions, plants employ a plethora of response mechanisms that provide broad phenotypic plasticity to allow the fine-tuning of the external cues related reactions. Molecular biology has been transformed by the major breakthroughs in high-throughput transcriptome sequencing and expression analysis using next-generation sequencing (NGS) technologies. These innovations have provided substantial progress in the identification of genomic regions as well as underlying basis influencing transcriptional and post-transcriptional regulation of abiotic stress response. Non-coding RNAs (ncRNAs), particularly microRNAs (miRNAs), short interfering RNAs (siRNAs), and long non-coding RNAs (lncRNAs), have emerged as essential regulators of plants abiotic stress response. However, shared traits in the biogenesis of ncRNAs and the coordinated cross-talk among ncRNAs mechanisms contribute to the complexity of these molecules and might play an essential part in regulating stress responses. Herein, we highlight the current knowledge of plant microRNAs, siRNAs, and lncRNAs, focusing on their origin, biogenesis, modes of action, and fundamental roles in plant response to abiotic stresses.

Barley microRNAs as metabolic sensors for soil nitrogen availability

Barley (Hordeum vulgare) is one of the most important crops in the world, ranking 4th in the worldwide production. Crop breeders are facing increasing environmental obstacles in the field, such as drought, salinity but also toxic over fertilization which not only impacts quality of the grain but also an yield. One of the most prevalent mechanisms of gene expression regulation in plants is microRNA-mediated silencing of target genes. We identified 13 barley microRNAs and 2 microRNAs* that are nitrogen excess responsive. Four microRNAs respond only in root, eight microRNAs only in shoot and one displays broad response in roots and shoots. We demonstrate that 2 microRNAs* are induced in barley shoot by nitrogen excess. For all microRNAs we identified putative target genes and confirmed microRNA-guided cleavage sites for ten out of thirteen mRNAs. None of the identified microRNAs or their target genes is known as nitrogen excess responsive. Analysis of expression pattern of thirteen target mRNAs and their cognate microRNAs showed expected correlations of their levels. The plant microRNAs analyzed are also known to respond to nitrogen deprivation and exhibit the opposite expression pattern when nitrogen excess/deficiency conditions are compared. Thus, they can be regarded as metabolic sensors of the regulation of nitrogen homeostasis in plants.

Characterization of microRNA genes from Pigeonpea (Cajanus cajan L.) and understanding their involvement in drought stress

MicroRNAs (miRNA) are non-coding 20-24 nucleotide long RNAs regulating gene expression. In this study, we have characterized and analysed expression of miRNAs in Pigeonpea by using bioinformatics and experimental tools. We identified 116 miRNAs belonging to 32 phylogenetic families. Further, transcription start sites of miRNA genes revealed abundance and unique arrangement of adenine at +1 and thymine at -1 position. Promoter analysis exhibited presence of 19 most prevalent motifs which comprises majorly of TATA box and MYC domains. In total, 252 miRNA-targets were identified and found to be involved in various developmental processes and stress responses. Moreover, genome-wide localization studies demonstrated clustering of cca-miRNA 395 and 169 genes. The tandem and segmental duplication events were observed suggesting miRNA genes have been originated parallelly with protein coding genes. The expression analysis revealed induction of cca-miR169a, 398a and 408 miRNAs under drought stress highlighting their involvement. Conversely, down-regulation of their putative targets (NFYA, SOD, and UCLA, respectively) confirmed regulatory role of miRNAs in their expression emphasising the negative relationship between these miRNAs and targets in Pigeonpea. This study reports vast repertoire of miRNA genes which further can be experimentally characterized to elucidate their functions in various biological processes and can be recommended for Pigeonpea improvement programs.

ROS Homeostasis in Abiotic Stress Tolerance in Plants

Climate change-induced abiotic stress results in crop yield and production losses. These stresses result in changes at the physiological and molecular level that affect the development and growth of the plant. Reactive oxygen species (ROS) is formed at high levels due to abiotic stress within different organelles, leading to cellular damage. Plants have evolved mechanisms to control the production and scavenging of ROS through enzymatic and non-enzymatic antioxidative processes. However, ROS has a dual function in abiotic stresses where, at high levels, they are toxic to cells while the same molecule can function as a signal transducer that activates a local and systemic plant defense response against stress. The effects, perception, signaling, and activation of ROS and their antioxidative responses are elaborated in this review. This review aims to provide a purview of processes involved in ROS homeostasis in plants and to identify genes that are triggered in response to abiotic-induced oxidative stress. This review articulates the importance of these genes and pathways in understanding the mechanism of resistance in plants and the importance of this information in breeding and genetically developing crops for resistance against abiotic stress in plants.

Temporal responses of conserved miRNAs to drought and their associations with drought tolerance and productivity in rice

BACKGROUND

Plant miRNAs play crucial roles in responses to drought and developmental processes. It is essential to understand the association of miRNAs with drought-tolerance (DT), as well as their impacts on growth, development, and reproduction (GDP). This will facilitate our utilization of rice miRNAs in breeding.

RESULTS

In this study, we investigated the time course of miRNA responses to a long-term drought among six rice genotypes by high-throughput sequencing. In total, 354 conserved miRNAs were drought responsive, representing obvious genotype- and stage-dependent patterns. The drought-responsive miRNAs (DRMs) formed complex regulatory network via their coexpression and direct/indirect impacts on the rice transcriptome. Based on correlation analyses, 211 DRMs were predicted to be associated with DT and/or GDP. Noticeably, 14.2% DRMs were inversely correlated with DT and GDP. In addition, 9 pairs of mature miRNAs, each derived from the same pre-miRNAs, were predicted to have opposite roles in regulating DT and GDP. This suggests a potential yield penalty if an inappropriate miRNA/pre-miRNA is utilized. miRNAs have profound impacts on the rice transcriptome reflected by great number of correlated drought-responsive genes. By regulating these genes, a miRNA could activate diverse biological processes and metabolic pathways to adapt to drought and have an influence on its GDP.

CONCLUSION

Based on the temporal pattern of miRNAs in response to drought, we have described the complex network between DRMs. Potential associations of DRMs with DT and/or GDP were disclosed. This knowledge provides valuable information for a better understanding in the roles of miRNAs play in rice DT and/or GDP, which can facilitate our utilization of miRNA in breeding.

Transcriptome analysis reveals the roles of stem nodes in cadmium transport to rice grain

BACKGROUND

Node is the central organ of transferring nutrients and ions in plants. Cadmium (Cd) induced crop pollution threatens the food safety. Breeding of low Cd accumulation cultivar is a chance to resolve this universal problem. This study was performed to identify tissue specific genes involved in Cd accumulation in different rice stem nodes. Panicle node and the first node under panicle (node I) were sampled in two rice cultivars: Xiangwanxian No. 12 (low Cd accumulation cultivar) and Yuzhenxiang (high Cd accumulation cultivar). RNA-seq analysis was performed to identify differentially expressed genes (DEGs) and microRNAs.

RESULTS

Xiangwanxian No. 12 had lower Cd concentration in panicle node, node I and grain compared with Yuzhenxiang, and node I had the highest Cd concentration in the two cultivars. RNA seq analysis identified 4535 DEGs and 70 miRNAs between the two cultivars. Most genesrelated to the "transporter activity", such as OsIRT1, OsNramp5, OsVIT2, OsNRT1.5A, and OsABCC1, play roles in blocking the upward transport of Cd. Among the genes related to "response to stimulus", we identified OsHSP70 and OsHSFA2d/B2c in Xiangwanxian No. 12, but not in Yuzhenxiang, were all down-regulated by Cd stimulus. The up-regulation of miRNAs (osa-miR528 and osa-miR408) in Xiangwanxian No. 12 played a potent role in lowering Cd accumulation via down regulating the expression of candidate genes, such as bZIP, ERF, MYB, SnRK1 and HSPs.

CONCLUSIONS

Both panicle node and node I of Xiangwanxian No. 12 played a key role in blocking the upward transportation of Cd, while node I played a critical role in Yuzhenxiang. Distinct expression patterns of various transporter genes such as OsNRT1.5A, OsNramp5, OsIRT1, OsVIT2 and OsABCC1 resulted in differential Cd accumulation in different nodes. Likewise, distinct expression patterns of these transporter genes are likely responsible for the low Cd accumulation in Xiangwanxian No. 12 cultivar. MiRNAs drove multiple transcription factors, such as OsbZIPs, OsERFs, OsMYBs, to play a role in Cd stress response.

Identification of drought-responsive miRNAs in Hippophae tibetana using high-throughput sequencing

MicroRNAs (miRNAs) play an important role in abiotic stress response in plants. However, the total miRNA profiles (miRNome) and drought-responsive miRNAs in H. tibetana have not been identified. In this study, we present the first report on the miRNome profiles of H. tibetana by high-throughput sequencing technology. 116 known and 4 predicted novel miRNAs were all identified in six H. tibetana samples. Moreover, to reveal the drought-responsive miRNAs in H. tibetana, we compared the miRNA profiles of H. tibetana grown under water sufficiency and drought stress. The results showed that 39 known miRNAs were up-regulated, while 34 miRNAs were downregulated under drought stress. Moreover, the expression of two novel miRNAs (novel_mir_24 and novel_mir_87) showed notable changes in response to drought stress. The target genes of these differentially expressed miRNAs were mainly enriched in cellular process, metabolic process, cell part, and response to stimulus. The identified drought-responsive miRNAs might be used for improving drought tolerance in H. tibetana and other plateau plants.

Differential expression of microRNAs in tomato leaves treated with different light qualities

Background

Light is the main source of energy and, as such, is one of the most important environmental factors for plant growth, morphogenesis, and other physiological responses. MicroRNAs (miRNAs) are endogenous non-coding RNAs that contain 21–24 nucleotides (nt) and play important roles in plant growth and development as well as stress responses. However, the role of miRNAs in the light response is less studied. We used tomato seedlings that were cultured in red light then transferred to blue light for 2 min to identify miRNAs related to light response by high-throughput sequencing.

Results

A total of 108 known miRNAs and 141 predicted novel miRNAs were identified in leaf samples from tomato leaves treated with the different light qualities. Among them, 15 known and 5 predicted novel miRNAs were differentially expressed after blue light treatment compared with the control (red light treatment). KEGG enrichment analysis showed that significantly enriched pathways included zeatin biosynthesis (ko00908), homologous recombination (ko03440), and plant hormone signal transduction (ko04075). Zeatin biosynthesis and plant hormone signal transduction are related to plant hormones, indicating that plant hormones play important roles in the light response.

Conclusion

Our results provide a theoretical basis for further understanding the role of miRNAs in the light response of plants.

Scion control of miRNA abundance and tree maturity in grafted avocado

BACKGROUND

Grafting is the common propagation method for avocado and primarily benefits orchard production by reducing the time to tree productivity. It also allows use of scions and rootstocks specifically selected for improved productivity and commercial acceptance. Rootstocks in avocado may be propagated from mature tree cuttings ('mature'), or from seed ('juvenile'). While the use of mature scion material hastens early bearing/maturity and economic return, the molecular factors involved in the role of the scion and/or rootstock in early bearing/reduced juvenility of the grafted tree are still unknown.

RESULTS

Here, we utilized juvenility and flowering associated miRNAs; miR156 and miR172 and their putative target genes to screen pre-graft and post-graft material in different combinations from avocado. The abundance of mature miR156, miR172 and the miR156 target gene SPL4, showed a strong correlation to the maturity of the scion and rootstock material in avocado. Graft transmissibility of miR156 and miR172 has been explored in annual plants. Here, we show that the scion may be responsible for grafted tree maturity involving these factors, while the rootstock maturity does not significantly influence miRNA abundance in the scion. We also demonstrate that the presence of leaves on cutting rootstocks supports graft success and contributes towards intergraft signalling involving the carbohydrate-marker TPS1.

CONCLUSION

Here, we suggest that the scion largely controls the molecular 'maturity' of grafted avocado trees, however, leaves on the rootstock not only promote graft success, but can influence miRNA and mRNA abundance in the scion. This constitutes the first study on scion and rootstock contribution towards grafted tree maturity using the miR156-SPL4-miR172 regulatory module as a marker for juvenility and reproductive competence.

Scion control of miRNA abundance and tree maturity in grafted avocado

BACKGROUND

Grafting is the common propagation method for avocado and primarily benefits orchard production by reducing the time to tree productivity. It also allows use of scions and rootstocks specifically selected for improved productivity and commercial acceptance. Rootstocks in avocado may be propagated from mature tree cuttings ('mature'), or from seed ('juvenile'). While the use of mature scion material hastens early bearing/maturity and economic return, the molecular factors involved in the role of the scion and/or rootstock in early bearing/reduced juvenility of the grafted tree are still unknown.

RESULTS

Here, we utilized juvenility and flowering associated miRNAs; miR156 and miR172 and their putative target genes to screen pre-graft and post-graft material in different combinations from avocado. The abundance of mature miR156, miR172 and the miR156 target gene SPL4, showed a strong correlation to the maturity of the scion and rootstock material in avocado. Graft transmissibility of miR156 and miR172 has been explored in annual plants. Here, we show that the scion may be responsible for grafted tree maturity involving these factors, while the rootstock maturity does not significantly influence miRNA abundance in the scion. We also demonstrate that the presence of leaves on cutting rootstocks supports graft success and contributes towards intergraft signalling involving the carbohydrate-marker TPS1.

CONCLUSION

Here, we suggest that the scion largely controls the molecular 'maturity' of grafted avocado trees, however, leaves on the rootstock not only promote graft success, but can influence miRNA and mRNA abundance in the scion. This constitutes the first study on scion and rootstock contribution towards grafted tree maturity using the miR156-SPL4-miR172 regulatory module as a marker for juvenility and reproductive competence.

Drought Response in Rice: The miRNA Story

As a semi-aquatic plant, rice requires water for proper growth, development, and orientation of physiological processes. Stress is induced at the cellular and molecular level when rice is exposed to drought or periods of low water availability. Plants have existing defense mechanisms in planta that respond to stress. In this review we examine the role played by miRNAs in the regulation and control of drought stress in rice through a summary of molecular studies conducted on miRNAs with emphasis on their contribution to drought regulatory networks in comparison to other plant systems. The interaction between miRNAs, target genes, transcription factors and their respective roles in drought-induced stresses is elaborated. The cross talk involved in controlling drought stress responses through the up and down regulation of targets encoding regulatory and functional proteins is highlighted. The information contained herein can further be explored to identify targets for crop improvement in the future.

Identification of microRNAS differentially regulated by water deficit in relation to mycorrhizal treatment in wheat

Arbuscular mycorrhizal fungi (AMF) are soil microrganisms that establish symbiosis with plants positively influencing their resistance to abiotic stresses. The aim of this work was to identify wheat miRNAs differentially regulated by water deficit conditions in presence or absence of AMF treatment. Small RNA libraries were constructed for both leaf and root tissues considering four conditions: control (irrigated) or water deficit in presence/absence of mycorrhizal (AMF) treatment. A total of 12 miRNAs were significantly regulated by water deficit in leaves: five in absence and seven in presence of AMF treatment. In roots, three miRNAs were water deficit-modulated in absence of mycorrhizal treatment while six were regulated in presence of it. The most represented miRNA family was miR167 that was regulated by water deficit in both leaf and root tissues. Interestingly, miR827-5p was differentially regulated in leaves in the absence of mycorrhizal treatment while it was water deficit-modulated in roots irrespective of AMF treatment. In roots, water deficit repressed miR827-5p, miR394, miR6187, miR167e-3p, and miR9666b-3p affecting transcription, RNA synthesis, protein synthesis, and protein modifications. In leaves, mycorrhizae modulated miR5384-3p and miR156e-3p affecting trafficking and cell redox homeostasis. DNA replication and transcription regulation should be targeted by the repression of miR1432-5p and miR166h-3p. This work provided interesting insights into the post-transcriptional mechanisms of wheat responses to water deficit in relation to mycorrhizal symbiosis.

Genome-wide identification and comparative analysis of drought-related microRNAs in two maize inbred lines with contrasting drought tolerance by deep sequencing

Drought has become one of the most serious abiotic stresses influencing crop production worldwide. Understanding the molecular regulatory networks underlying drought adaption and tolerance in crops is of great importance for future breeding. microRNAs (miRNAs), as important components of post-transcriptional regulation, play crucial roles in drought response and adaptation in plants. Here, we report a miRNome analysis of two maize inbred lines with contrasting levels of drought tolerance under soil drought in the field. Differential expression analysis showed 11 and 34 miRNAs were uniquely responded to drought in H082183 (drought tolerant) and Lv28 (drought sensitive), respectively, in leaves. In roots, 19 and 23 miRNAs uniquely responded to drought in H082183 and Lv28, respectively. Expression analysis of these drought-responsive miRNA-mRNA modules revealed miR164-MYB, miR164-NAC, miR159-MYB, miR156-SPL and miR160-ARF showed a negative regulatory relationship. Further analysis showed that the miR164-MYB and miR164-NAC modules in the tolerant line modulated the stress response in an ABA (abscisic acid)-dependent manner, while the miR156-SPL and miR160-ARF modules in the sensitive line participated in the inhibition of metabolism in drought-exposed leaves. Together, our results provide new insight into not only drought-tolerance-related miRNA regulation networks in maize but also key miRNAs for further characterization and improvement of maize drought tolerance.

OsmiR528 Enhances Cold Stress Tolerance by Repressing Expression of Stress Response-related Transcription Factor Genes in Plant Cells

Background:

MicroRNAs participate in many molecular mechanisms and signaling trans-duction pathways that are associated with plant stress tolerance by repressing expression of their target genes. However, how microRNAs enhance tolerance to low temperature stress in plant cells remains elusive.

Objective:

In this investigation, we demonstrated that overexpression of the rice microRNA528 (Os-miR528) increases cell viability, growth rate, antioxidants content, ascorbate peroxidase (APOX) activi-ty, and superoxide dismutase (SOD) activity and decreases ion leakage rate and thiobarbituric acid reac-tive substances (TBARS) under low temperature stress in Arabidopsis (Arabidopsis thaliana), pine (Pi-nus elliottii), and rice (Oryza sativa).

Methods:

To investigate the potential mechanism of OsmiR528 in increasing cold stress tolerance, we examined expression of stress-associated MYB transcription factors OsGAMYB-like1, OsMYBS3, OsMYB4, OsMYB3R-2, OsMYB5, OsMYB59, OsMYB30, OsMYB1R, and OsMYB20 in rice cells by qRT-PCR.

Results:

Our experiments demonstrated that OsmiR528 decreases expression of transcription factor OsMYB30 by targeting a F-box domain containing protein gene (Os06g06050), which is a positive regulator of OsMYB30. In OsmiR528 transgenic rice, reduced OsMYB30 expression results in in-creased expression of BMY genes OsBMY2, OsBMY6, and OsBMY10. The transcript levels of the OsBMY2, OsBMY6, and OsBMY10 were elevated by OsMYB30 knockdown, but decreased by Os-MYB30 overexpression in OsmiR528 transgenic cell lines, suggesting that OsmiR528 increases low temperature tolerance by modulating expression of stress response-related transcription factor.

Conclusion:

Our experiments provide novel information in increasing our understanding in molecular mechanisms of microRNAs-associated low temperature tolerance and are valuable in plant molecular breeding from monocotyledonous, dicotyledonous, and gymnosperm plants.

Differential Expression of MicroRNAs During Root Formation in Taxus Chinensis Var. mairei Cultivars

MicroRNAs (miRNAs) have been shown to play key roles in the regulation of plant growth and development by modifying the expression of their target genes. However, the influence of miRNAs on root formation and development in woody plants, such as Taxus chinensis, remains largely unknown. In the current study, we explored the phytohormone-response and nutrition-response miRNA expression profiles during T. chinensis rooting by quantitative real-time PCR (qPCR). We identified six phytohormone-response miRNAs, namely, miR164a, miR165, miR167a, miR171b, miR319, and miR391, and eight nutrition-response miRNAs, namely, miR169b, miR395a, miR399c, miR408, miR826, miR827, miR857, and miR2111a, that were differentially expressed at different rooting phases of T. chinensis. Using northern blot analysis of the putative target genes of these miRNAs, we detected the relative gene expression changes of the target genes. Taken together, our results suggest that miRNAs are involved in root formation of T. chinensis and that miRNAs may play important regulatory roles in primary root, crown root, and root hair formation by targeting phytohormone and/or nutrition response genes in T. chinensis. For the first time, these results expand our understanding of the molecular mechanisms of plant root formation and development in a conifer species.

Juvenility and Vegetative Phase Transition in Tropical/Subtropical Tree Crops

In plants, juvenile to adult phase transition is regulated by the sequential activity of two microRNAs: miR156 and miR172. A decline in miR156 and increase in miR172 abundance is associated with phase transition. There is very limited information on phase transition in economically important horticultural tree crops, which have a significantly long vegetative phase affecting fruit bearing. Here, we profiled various molecular cues known to be involved in phase transition and flowering, including the microRNAs miR156 and miR172, in three horticultural tree crops: avocado (Persea americana), mango (Mangifera indica), and macadamia (Macadamia integrifolia). We observed that miR156 expression decreases as these trees age and can potentially be used as a juvenility marker. Consistent with findings in annual plants, we also observed conserved regulation of the miR156-SPL3/4/5 regulatory module in these genetically distant tree crops, suggesting that this pathway may play a highly conserved role in vegetative identity. Meanwhile, the abundance of miR172 and its target AP2-like genes as well as the accumulation level of SPL9 transcripts were not related with plant age in these crops except in avocado where miR172 expression increased steadily. Finally, we demonstrate that various floral genes, including AP1 and SOC1 were upregulated in the reproductive phase and can be used as potential markers for the reproductive phase transition. Overall, this study provides an insight into the molecular associations of juvenility and phase transition in horticultural trees where crop breeding and improvement are encumbered by long juvenile phases.

Genome-wide identification of drought-responsive microRNAs in two sets of Malus from interspecific hybrid progenies

Drought stress can negatively impact apple fruit quality and yield. Apple microRNAs (miRNAs) participate in apple tree and fruit development, as well as in biotic stress tolerance; however, it is largely unknown whether these molecules are involved in the drought response. To identify drought-responsive miRNAs in Malus, we first examined the drought stress tolerance of ten F₁ progenies of R3 (M. × domestica) × M. sieversii. We performed Illumina sequencing on pooled total RNA from both drought-tolerant and drought-sensitive plants. The sequencing results identified a total of 206 known miRNAs and 253 candidate novel miRNAs from drought-tolerant plants and drought-sensitive plants under control or drought conditions. We identified 67 miRNAs that were differentially expressed in drought-tolerant plants compared with drought-sensitive plants under drought conditions. Under drought stress, 61 and 35 miRNAs were differentially expressed in drought-tolerant and drought-sensitive plants, respectively. We determined the expression levels of seven out of eight miRNAs by stem-loop qPCR analysis. We also predicted the target genes of all differentially expressed miRNAs and identified the expression of some genes. Gene Ontology analyses indicated that the target genes were mainly involved in stimulus response and cellular and metabolic processes. Finally, we confirmed roles of two miRNAs in apple response to mannitol. Our results reveal candidate miRNAs and their associated mRNAs that could be targeted for improving drought tolerance in Malus species, thus providing a foundation for understanding the molecular networks involved in the response of apple trees to drought stress.

Natural Variation in OsLG3 Increases Drought Tolerance in Rice by Inducing ROS Scavenging

Improving the performance of rice (Oryza sativa) under drought stress has the potential to significantly affect rice productivity. Here, we report that the ERF family transcription factor OsLG3 positively regulates drought tolerance in rice. In our previous work, we found that OsLG3 has a positive effect on rice grain length without affecting grain quality. In this study, we found that OsLG3 was more strongly expressed in upland rice than in lowland rice under drought stress conditions. By performing candidate gene association analysis, we found that natural variation in the promoter of OsLG3 is associated with tolerance to osmotic stress in germinating rice seeds. Overexpression of OsLG3 significantly improved the tolerance of rice plants to simulated drought, whereas suppression of OsLG3 resulted in greater susceptibility. Phylogenetic analysis indicated that the tolerant allele of OsLG3 may improve drought tolerance in cultivated japonica rice. Introgression lines and complementation transgenic lines containing the elite allele of OsLG3^IRAT109 showed increased drought tolerance, demonstrating that natural variation in OsLG3 contributes to drought tolerance in rice. Further investigation suggested that OsLG3 plays a positive role in drought stress tolerance in rice by inducing reactive oxygen species scavenging. Collectively, our findings reveal that natural variation in OsLG3 contributes to rice drought tolerance and that the elite allele of OsLG3 is a promising genetic resource for the development of drought-tolerant rice varieties.

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.

Identification of Blueberry miRNAs and Their Targets Based on High-Throughput Sequencing and Degradome Analyses

miRNAs are important regulators of plant gene expression. To better characterize their functions, we applied high-throughput sequencing and degradome analyses to investigate three blueberry (Vaccinium ashei) tissues. A total of 127 known and 101 novel miRNAs were identified. Moreover, 141 targets for 42 known and 19 novel miRNAs were experimentally validated by degradome sequencing. A functional analysis of these miRNA targets revealed they were associated with diverse biological activities and several pathways, e.g., anthocyanin biosynthesis and cytokinin signal transduction. The data presented herein expand our understanding of the regulation of blueberry miRNAs during floral and fruit development stages. They may also provide new insights into the roles of miRNAs during anthocyanin biosynthesis in blueberry fruits.

Formation of Protein Disulfide Bonds Catalyzed by OsPDIL1;1 is Mediated by MicroRNA5144-3p in Rice

Correct folding of proteins in the endoplasmic reticulum is important for their stability and function under stress. The protein disulfide isomerase (PDI) OsPDIL1;1 is a key protein-folding catalyst in rice (Oryza sativa L.). Here, microRNA5144 (osa-miR5144-3p) is reported to mediate the formation of protein disulfide bonds via targeting OsPDIL1;1 mRNA in rice seeds and seedlings during development and under conditions of abiotic stress, respectively. Expression analysis of transgenic rice and identification of cleavage sites showed that OsPDIL1;1 mRNA is a target of osa-miR5144-3p. Expression of osa-miR5144-3p and OsPDIL1;1 was shown to be inversely regulated in developing organs and under abiotic stress. The down-regulation of osa-miR5144-3p or overexpression of OsPDIL1;1 in transgenic rice showed increased total protein-disulfide bond content, compared with the wild type. This indicates that protein-disulfide bond formation is enhanced by down-regulation of osa-miR5144-3p or overexpression of OsPDIL1;1. These transgenic rice plants also displayed strong resistance to salinity and mercury stress, in comparison with the wild type. In contrast, the transgenic rice plants overexpressing osa-miR5144-3p or down-regulating OsPDIL1;1 had a lower protein-disulfide bond content; they were susceptible to abiotic stress and produced abnormal grains with small and loosely packed starch granules. These results indicate that protein-disulfide bond formation catalyzed by OsPDIL1;1 is modulated by osa-miR5144-3p in rice during development and is involved in resistance to abiotic stress.

Identification of miRNA-mediated drought responsive multi-tiered regulatory network in drought tolerant rice, Nagina 22

Comparative characterization of microRNA-mediated stress regulatory networks in contrasting rice cultivars is critical to decipher plant stress response. Consequently, a multi-level comparative analysis, using sRNA sequencing, degradome analysis, enzymatic and metabolite assays and metal ion analysis, in drought tolerant and sensitive rice cultivars was conducted. The study identified a group of miRNAs "Cultivar-specific drought responsive" (CSDR)-miRNAs (osa-miR159f, osa-miR1871, osa-miR398b, osa-miR408-3p, osa-miR2878-5p, osa-miR528-5p and osa-miR397a) that were up-regulated in the flag-leaves of tolerant cultivar, Nagina 22 (N22) and Vandana, but down-regulated in the sensitive cultivar, Pusa Basmati 1 (PB1) and IR64, during drought. Interestingly, CSDR-miRNAs target several copper-protein coding transcripts like plantacyanins, laccases and Copper/Zinc superoxide dismutases (Cu/Zn SODs) and are themselves found to be similarly induced under simulated copper-starvation in both N22 and PB1. Transcription factor OsSPL9, implicated in Cu-homeostasis also interacted with osa-miR408-3p and osa-miR528-5p promoters. Further, N22 flag leaves showed lower SOD activity, accumulated ROS and had a higher stomata closure. Interestingly, compared to PB1, internal Cu levels significantly decreased in the N22 flag-leaves, during drought. Thus, the study identifies the unique drought mediated dynamism and interplay of Cu and ROS homeostasis, in the flag leaves of drought tolerant rice, wherein CSDR-miRNAs play a pivotal role.

Identification of functionally important microRNAs from rice inflorescence at heading stage of a qDTY4.1-QTL bearing Near Isogenic Line under drought conditions

A cross between IR64 (high-yielding but drought-susceptible) and Aday Sel (drought-tolerant) rice cultivars yielded a stable line with enhanced grain yield under drought screening field trials at International Rice Research Institute. The major effect qDTY4.1 drought tolerance and yield QTL was detected in the IR77298-14-1-2-10 Backcrossed Inbred Line (BIL) and its IR87705-7-15-B Near Isogenic Line (NIL) with 93.9% genetic similarity to IR64. Although rice yield is extremely susceptible to water stress at reproductive stage, currently, there is only one report on the detection of drought-responsive microRNAs in inflorescence tissue of a Japonica rice line. In this study, more drought-responsive microRNAs were identified in the inflorescence tissues of IR64, IR77298-14-1-2-10 and IR87705-7-15-B via next-generation sequencing. Among the 32 families of inflorescence-specific non-conserved microRNAs that were identified, 22 families were up-regulated in IR87705-7-15-B. Overall 9 conserved and 34 non-conserved microRNA families were found as drought-responsive in rice inflorescence with 5 conserved and 30 non-conserved families induced in the IR87705-7-15-B. The observation of more drought-responsive non-conserved microRNAs may imply their prominence over conserved microRNAs in drought response mechanisms of rice inflorescence. Gene Ontology annotation analysis on the target genes of drought-responsive microRNAs identified in IR87705-7-15-B revealed over-representation of biological processes including development, signalling and response to stimulus. Particularly, four inflorescence-specific microRNAs viz. osa-miR5485, osa-miR5487, osa-miR5492 and osa-miR5517, and two non-inflorescence specific microRNAs viz. osa-miR169d and osa-miR169f.2 target genes that are involved in flower or embryonic development. Among them, osa-miR169d, osa-miR5492 and osa-miR5517 are related to flowering time control. It is also worth mentioning that osa-miR2118 and osa-miR2275, which are implicated in the biosynthesis of rice inflorescence-specific small interfering RNAs, were induced in IR87705-7-15-B but repressed in IR77298-14-1-2-10. Further, gene search within qDTY4.1 QTL region had identified multiple copies of NBS-LRR resistance genes (potential target of osa-miR2118), subtilisins and genes implicated in stomatal movement, ABA metabolism and cuticular wax biosynthesis.

Identification of drought-responsive microRNAs in tomato using high-throughput sequencing

Drought is a major abiotic stress affecting crop productivity and quality. As a class of noncoding RNA, microRNA (miRNA) plays important roles in plant growth, development, and stress response. However, their response and roles in tomato drought stress is largely unknown. Here, by using high-throughput sequencing, we compared the miRNA profiles before and after drought treatment in two tomato genotypes: M82, a drought-sensitive cultivated tomato (Solanum lycopersicum), and IL2-5, a drought-tolerant introgression line derived from M82 and the tomato wild species S. pennellii (LA0716). A total of 108 conserved and 208 novel miRNAs were identified, among them, 32 and 68 were significantly changed in expression after stress. Further, 1936 putative target genes were predicted for those differentially-expressed miRNAs. Gene ontology and pathway analysis showed that many of the target genes were involved in stress resistance, such as genes in GO terms including response to stress, defense response, response to stimulus, phosphorylation, and signal transduction. Our results suggested that miRNAs play an essential role in the drought response of tomato. This work will help to further characterize specific miRNAs functioning in drought tolerance.

Comprehensive survey and evolutionary analysis of genome-wide miRNA genes from ten diploid Oryza species

BACKGROUND

MicroRNAs (miRNAs) are non-coding RNAs that play versatile roles in post-transcriptional gene regulation. Although much is known about their biogenesis, and gene regulation very little is known about their evolutionary relation among the closely related species.

RESULT

All the orthologous miRNA genes of Oryza sativa (japonica) from 10 different Oryza species were identified, and the evolutionary changes among these genes were analysed. Significant differences in the expansion of miRNA gene families were observed across the Oryza species. Analysis of the nucleotide substitution rates indicated that the mature sequences show the least substitution rates among the different regions of miRNA genes, and also show a very much less substitution rates as compared to that of all protein-coding genes across the Oryza species. Evolution of miRNA genes was also found to be contributed by transposons. A non-neutral selection was observed at 80 different miRNA loci across Oryza species which were estimated to have lost ~87% of the sequence diversity during the domestication. The phylogenetic analysis revealed that O. longistaminata diverged first among the AA-genomes, whereas O. brachyantha and O. punctata appeared as the eminent out-groups. The miR1861 family organised into nine distinct compact clusters in the studied Oryza species except O. brachyantha. Further, the expression analysis showed that 11 salt-responsive miRNAs were differentially regulated between O. coarctata and O. glaberrima.

CONCLUSION

Our study provides the evolutionary dynamics in the miRNA genes of 10 different Oryza species which will support more investigations about the structural and functional organization of miRNA genes of Oryza species.

Profiling of drought-responsive microRNA and mRNA in tomato using high-throughput sequencing

BACKGROUND

Abiotic stresses cause severe loss of crop production. Among them, drought is one of the most frequent environmental stresses, which limits crop growth, development and productivity. Plant drought tolerance is fine-tuned by a complex gene regulatory network. Understanding the molecular regulation of this polygenic trait is crucial for the eventual success to improve plant yield and quality. Recent studies have demonstrated that microRNAs play critical roles in plant drought tolerance. However, little is known about the microRNA in drought response of the model plant tomato. Here, we described the profiling of drought-responsive microRNA and mRNA in tomato using high-throughput next-generation sequencing.

RESULTS

Drought stress was applied on the seedlings of M82, a drought-sensitive cultivated tomato genotype, and IL9-1, a drought-tolerant introgression line derived from the stress-resistant wild species Solanum pennellii LA0716 and M82. Under drought, IL9-1 performed superior than M82 regarding survival rate, H2O2 elimination and leaf turgor maintenance. A total of four small RNA and eight mRNA libraries were constructed and sequenced using Illumina sequencing technology. 105 conserved and 179 novel microRNAs were identified, among them, 54 and 98 were differentially expressed upon drought stress, respectively. The majority of the differentially-expressed conserved microRNAs was up-regulated in IL9-1 whereas down-regulated in M82. Under drought stress, 2714 and 1161 genes were found to be differentially expressed in M82 and IL9-1, respectively, and many of their homologues are involved in plant stress, such as genes encoding transcription factor and protein kinase. Various pathways involved in abiotic stress were revealed by Gene Ontology and pathway analysis. The mRNA sequencing results indicated that most of the target genes were regulated by their corresponding microRNAs, which suggested that microRNAs may play essential roles in the drought tolerance of tomato.

CONCLUSION

In this study, numerous microRNAs and mRNAs involved in the drought response of tomato were identified using high-throughput sequencing, which will provide new insights into the complex regulatory network of plant adaption to drought stress. This work will also help to exploit new players functioning in plant drought-stress tolerance.

Identification of Drought-Responsive MicroRNAs from Roots and Leaves of Alfalfa by High-Throughput Sequencing

Alfalfa, an important forage legume, is an ideal crop for sustainable agriculture and a potential crop for bioenergy resources. Drought, one of the most common environmental stresses, substantially affects plant growth, development, and productivity. MicroRNAs (miRNAs) are newly discovered gene expression regulators that have been linked to several plant stress responses. To elucidate the role of miRNAs in drought stress regulation of alfalfa, a high-throughput sequencing approach was used to analyze 12 small RNA libraries comprising of four samples, each with three biological replicates. From the 12 libraries, we identified 348 known miRNAs belonging to 80 miRNA families, and 281 novel miRNAs, using Mireap software. Eighteen known miRNAs in roots and 12 known miRNAs in leaves were screened as drought-responsive miRNAs. With the exception of miR319d and miR157a which were upregulated under drought stress, the expression pattern of drought-responsive miRNAs was different between roots and leaves in alfalfa. This is the first study that has identified miR3512, miR3630, miR5213, miR5294, miR5368 and miR6173 as drought-responsive miRNAs. Target transcripts of drought-responsive miRNAs were computationally predicted. All 447 target genes for the known miRNAs were predicted using an online tool. This study provides a significant insight on understanding drought-responsive mechanisms of alfalfa.

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.

Identification of Submergence-Responsive MicroRNAs and Their Targets Reveals Complex MiRNA-Mediated Regulatory Networks in Lotus (Nelumbo nucifera Gaertn)

MicroRNAs (miRNAs) are endogenous non-coding RNAs with important regulatory functions in plant development and stress responses. However, their population abundance in lotus (Nelumbo nucifera Gaertn) has so far been poorly described, particularly in response to stresses. In this work, submergence-related miRNAs and their target genes were systematically identified, compared, and validated at the transcriptome-wide level using high-throughput sequencing data of small RNA, Mrna, and the degradome. A total of 128 known and 20 novel miRNAs were differentially expressed upon submergence. We identified 629 target transcripts for these submergence-responsive miRNAs. Based on the miRNA expression profiles and GO and KEGG annotation of miRNA target genes, we suggest possible molecular responses and physiological changes of lotus in response to submergence. Several metabolic, physiological and morphological adaptations-related miRNAs, i.e., NNU_far-miR159, NNU_gma-miR393h, and NNU_aly-miR319c-3p, were found to play important regulatory roles in lotus response to submergence. This work will contribute to a better understanding of miRNA-regulated adaption responses of lotus to submergence stress.