Microarray-based analysis of stress-regulated microRNAs in Arabidopsis thaliana

MicroRNA-mRNA expression profiles and their potential role in cadmium stress response in Brassica napus

BACKGROUND

Oilseed rape is an excellent candidate for phytoremediation of cadmium (Cd) contaminated soils given its advantages of high biomass, fast growth, moderate metal accumulation, ease of harvesting, and metal tolerance, but the cadmium response pathways in this species (Brassica napus) have yet to be fully elucidated. A combined analysis of miRNA and mRNA expression to infer Cd-induced regulation has not been reported in B. napus.

RESULTS

We characterized concurrent changes in miRNA and mRNA profiles in the roots and shoots of B. napus seedlings after 10 days of 10 mg/L Cd²⁺ treatment. Cd treatment significantly affected the expression of 22 miRNAs belonging to 11 families in the root and 29 miRNAs belonging to 14 miRNA families in the shoot. Five miRNA families (MIR395, MIR397, MIR398, MIR408 and MIR858) and three novel miRNAs were differentially expressed in both tissues. A total of 399 differentially expressed genes (DEGs) in the root and 389 DEGs in the shoot were identified, with very little overlap between tissue types. Eight anti-regulation miRNA-mRNA interaction pairs in the root and eight in the shoot were identified in response to Cd and were involved in key plant stress response pathways: for example, four genes targeted by miR398 were involved in a pathway for detoxification of superoxide radicals. Cd stress significantly impacted the photosynthetic pathway. Transcription factor activation, antioxidant response pathways and secondary metabolic processes such as glutathione (GSH) and phenylpropanoid metabolism were identified as major components for Cd-induced response in both roots and shoots.

CONCLUSIONS

Combined miRNA and mRNA profiling revealed miRNAs, genes and pathways involved in Cd response which are potentially critical for adaptation to Cd stress in B. napus. Close crosstalk between several Cd-induced miRNAs and mRNAs was identified, shedding light on possible mechanisms for response to Cd stress in underground and aboveground tissues in B. napus. The pathways, genes, and miRNAs identified here will be valuable targets for future improvement of cadmium tolerance in B. napus.

Dynamic architecture and regulatory implications of the miRNA network underlying the response to stress in melon

miRNAs are small RNAs that regulate mRNAs at both transcriptional and posttranscriptional level. In plants, miRNAs are involved in the regulation of different processes including development and stress-response. Elucidating how stress-responsive miRNAs are regulated is key to understand the global response to stress but also to develop efficient biotechnological tools that could help to cope with stress. Here, we describe a computational approach based on sRNA sequencing, transcript quantification and degradome data to analyse the accumulation, function and structural organization of melon miRNAs reactivated under seven biotic and abiotic stress conditions at two and four days post-treatment. Our pipeline allowed us to identify fourteen stress-responsive miRNAs (including evolutionary conserved such as miR156, miR166, miR172, miR319, miR398, miR399, miR894 and miR408) at both analysed times. According to our analysis miRNAs were categorized in three groups showing a broad-, intermediate- or narrow- response range. miRNAs reactive to a broad range of environmental cues appear as central components in the stress-response network. The strictly coordinated response of miR398 and miR408 (broad response-range) to the seven stress treatments during the period analysed here reinforces this notion. Although both, the amplitude and diversity of the miRNA-related response to stress changes during the exposition time, the architecture of the miRNA-network is conserved. This organization of miRNA response to stress is also conserved in rice and soybean supporting the conservation of miRNA-network organization in other crops. Overall, our work sheds light into how miRNA networks in plants organize and function during stress.

Diurnal Regulation of Plant Epidermal Wax Synthesis through Antagonistic Roles of the Transcription Factors SPL9 and DEWAX

Plant surface waxes form an outer barrier that protects the plant from many forms of environmental stress. The deposition of cuticular waxes on the plant surface is regulated by external environmental changes, including light and dark cycles. However, the underlying molecular mechanisms controlling light regulation of wax production are still poorly understood, especially at the posttranscriptional level. In this paper, we report the regulation of cuticular wax production by the miR156-SPL9 (SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 9) module in Arabidopsis (Arabidopsis thaliana). When compared with wild-type plants, miR156 and SPL9 mutants showed significantly altered cuticular wax amounts in both stems and leaves. Furthermore, it was found that SPL9 positively regulates gene expression of the alkane-forming enzyme ECERIFERUM1 (CER1), as well as the primary (1-) alcohol-forming enzyme ECERIFERUM4 (CER4), to enhance alkane and 1-alcohol synthesis, respectively. Our results indicate that complex formation of SPL9 with a negative regulator of wax synthesis, DEWAX, will hamper SPL9 DNA binding ability, possibly by interfering with SPL9 homodimerization. Combined with their diurnal gene and protein expressions, this dynamic repression-activation transcriptional module defines a dynamic mechanism that may allow plants to optimize wax synthesis during daily cycles. These findings provide a regulatory framework for environmental signal integration in the regulation of wax synthesis.

The miR396-GRF Regulatory Module Controls the Embryogenic Response in Arabidopsis via an Auxin-Related Pathway

In plants, microRNAs have been indicated to control various developmental processes, including somatic embryogenesis (SE), which is triggered in the in vitro cultured somatic cells of plants. Although a transcriptomic analysis has indicated that numerous MIRNAs are differentially expressed in the SE of different plants, the role of specific miRNAs in the embryogenic reprogramming of the somatic cell transcriptome is still poorly understood. In this study, we focused on performing a functional analysis of miR396 in SE given that the transcripts of MIR396 genes and the mature molecules of miR396 were found to be increased during an SE culture of Arabidopsis. In terms of miR396 in embryogenic induction, we observed the SE-associated expression pattern of MIR396b in explants of the β-glucuronidase (GUS) reporter line. In order to gain insight into the miR396-controlled mechanism that is involved in SE induction, the embryogenic response of mir396 mutants and the 35S:MIR396b overexpressor line to media with different 2,4-Dichlorophenoxyacetic acid (2,4-D) concentrations was evaluated. The results suggested that miR396 might contribute to SE induction by controlling the sensitivity of tissues to auxin treatment. Within the targets of miR396 that are associated with SE induction, we identified genes encoding the GROWTH-REGULATING FACTOR (GRF) transcription factors, including GRF1, GRF4, GRF7, GRF8, and GRF9. Moreover, the study suggested a regulatory relationship between miR396, GRF, and the PLETHORA (PLT1 and PLT2) genes during SE induction. A complex regulatory relationship within the miR396–GRF1/4/8/9–PLT1/2 module that involves the negative and positive control of GRFs and PLT (respectively) by miR396 might be assumed.

Small RNA and degradome deep sequencing reveal respective roles of cold-related microRNAs across Chinese wild grapevine and cultivated grapevine

Background

Chinese wild grapevine (Vitis amurensis) has remarkable cold stress tolerance, exceeding that of the common cultivated grapevine (Vitis vinifera L.).

Result

Here, we surveyed the expression dynamics of microRNAs (miRNAs) across Chinese wild grapevine (cv. Beibinghong) and cultivated grapevine (cv. Cabernet Sauvignon) under cold stress using high-throughput sequencing. We identified 186 known miRNAs in cultivated grape and 427 known miRNAs in Beibinghong. Of the identified miRNAs, 59 are conserved miRNAs orthologous in Cabernet Sauvignon and Beibinghong. In addition, 105 and 129 novel miRNAs were identified in Cabernet Sauvignon and Beibinghong, respectively. The expression of some miRNAs was related to cold stress both in Cabernet Sauvignon and Beibinghong. Many cold-related miRNAs in Cabernet Sauvignon and Beibinghong were predicted to target stress response-related genes such as MYB, WRKY, bHLH transcription factor genes, and heat shock protein genes. However, the expression tendency under cold treatment of many of these miRNAs was different between Cabernet Sauvignon and Beibinghong. Different modes of expression of cultivated and Chinese wild grape miRNAs were indicated in key pathways under cold stress by degradome, target prediction, GO, and KEGG analyses.

Conclusion

Our study indicated three likely reasons that led to the different cold stress tolerance levels of Cabernet Sauvignon and Beibinghong. Specifically, there may be (1) differential expression of orthologous miRNAs between cultivated grapevine and Chinese wild grape; (2) species-specific miRNAs or target genes; or (3) different regulatory models of miRNAs in cultivated and Chinese wild grape in some key pathways.

Knockdown of ghAlba_4 and ghAlba_5 Proteins in Cotton Inhibits Root Growth and Increases Sensitivity to Drought and Salt Stresses

We found 33, 17, and 20 Alba genes in Gossypium hirsutum, Gossypium arboretum, and Gossypium raimondii, respectively. The Alba protein lengths ranged from 62 to 312 aa, the molecular weight (MW) from 7.003 to 34.55 kDa, grand average hydropathy values of -1.012 to 0.609 and isoelectric (pI) values of -3 to 11. Moreover, miRNAs such as gra-miR8770 targeted four genes, gra-miR8752 and gra-miR8666 targeted three genes, and each and gra-miR8657 a, b, c, d, e targeted 10 genes each, while the rests targeted 1 to 2 genes each. Similarly, various cis-regulatory elements were detected with significant roles in enhancing abiotic stress tolerance, such as CBFHV (RYCGAC) with a role in cold stress acclimation among others. Two genes, Gh_D01G0884 and Gh_D01G0922, were found to be highly induced under water deficit and salt stress conditions. Through virus-induced gene silencing (VIGS), the VIGS cotton plants were found to be highly susceptible to both water deficit and salt stresses; the VIGS plants exhibited a significant reduction in root growth, low cell membrane stability (CMS), saturated leaf weight (SLW), chlorophyll content levels, and higher excised leaf water loss (ELWL). Furthermore, the stress-responsive genes and ROS scavenging enzymes were significantly reduced in the VIGS plants compared to either the wild type (WT) and or the positively controlled plants. The VIGS plants registered higher concentration levels of hydrogen peroxide and malondialdehyde, with significantly lower levels of the various antioxidants evaluated an indication that the VIGS plants were highly affected by salt and drought stresses. This result provides a key foundation for future exploration of the Alba proteins in relation to abiotic stress.

The miRNA-mRNA Networks Involving Abnormal Energy and Hormone Metabolisms Restrict Tillering in a Wheat Mutant dmc

Tillers not only determine plant architecture but also influence crop yield. To explore the miRNA regulatory network restraining tiller development in a dwarf-monoculm wheat mutant (dmc) derived from Guomai 301 (wild type, WT), we employed miRNome and transcriptome integrative analysis, real-time qRT-PCR, histochemistry, and determinations of the key metabolites and photosynthesis parameters. A total of 91 differentially expressed miRNAs (DEMs) were identified between dmc and WT. Among them, 40 key DEMs targeted 45 differentially expressed genes (DEGs) including the key DEGs encode growth-regulating factors (GRF), auxin response factors (ARF), and other proteins involved in the metabolisms of hormones and carbohydrates, etc. Compared with WT, both the chlorophyll contents and the photosynthesis rate were lower in dmc. The contents of glucose, sucrose, fructose, and maltose were lower in dmc. The contents of auxin (IAA) and zeatin (ZA) were significantly lower, but gibberellin (GA) was significantly higher in the tiller tissues of dmc. This research demonstrated that the DEMs regulating hormone and carbohydrate metabolisms were important causes for dmc to not tiller. A primary miRNA-mRNA regulatory model for dmc tillering was established. The lower photosynthesis rate, insufficient energy, and abnormal hormone metabolisms restrict tillering in dmc.

Deficiency in the double-stranded RNA binding protein HYPONASTIC LEAVES1 increases sensitivity to the endoplasmic reticulum stress inducer tunicamycin in Arabidopsis

Objective

microRNA (miRNA) is a small non-coding RNA that regulates gene expression by sequence-dependent binding to protein-coding mRNA in eukaryotic cells. In plants, miRNA plays important roles in a plethora of physiological processes, including abiotic and biotic stress responses. The present study was conducted to investigate whether miRNA-mediated regulation is important for the endoplasmic reticulum (ER) stress response in Arabidopsis.

Results

We found that hyl1 mutant plants are more sensitive to tunicamycin, an inhibitor of N-linked glycosylation that causes ER stress than wild-type plants. Other miRNA-related mutants, se and ago1, exhibited similar sensitivity to the wild-type, indicating that the hypersensitive phenotype is attributable to the loss-of-function of HYL1, rather than deficiency in general miRNA biogenesis and function. However, the transcriptional response of select ER stress-responsive genes in hyl1 mutant plants was indistinguishable from that of wild-type plants, suggesting that the loss-of-function of HYL1 does not affect the ER stress signaling pathways.

Overexpression of soybean miR169c confers increased drought stress sensitivity in transgenic Arabidopsis thaliana

The miR169 family, a large-scale microRNA gene family conserved in plants, is involved in stress responses, although how soybean miR169 functions in response to drought stress remains unclear. We show that gma-miR169c exerts a negative regulatory role in the response to drought stress by inhibiting the expression of its target gene, nuclear factor Y-A (NF-YA). A real-time RT-PCR analysis indicated that gma-miR169c is widely expressed in soybean tissues and induced by polyethylene glycol (PEG), high salt, cold stress and abscisic acid (ABA). Histochemical ß-glucuronidase (GUS) staining showed that the gma-miR169c promoter drives GUS reporter gene expression in various transgenic Arabidopsis tissues, and the stress-induced pattern was confirmed in transgenic Arabidopsis and transgenic soybean hairy roots. Arabidopsis overexpressing gma-miR169c is more sensitive to drought stress, with reduced survival, accelerated leaf water loss, and shorter root length than wild-type plants. We identified a precise cleavage site for 10 gma-miR169c targets and found reduced transcript levels of the AtNFYA1 and AtNFYA5 transcription factors in gma-miR169c-overexpressing Arabidopsis and reduced expression of the stress response genes AtRD29A, AtRD22, AtGSTU25 and AtCOR15A. These results indicate that gma-miR169c plays a negative regulatory role in drought stress and is a candidate miRNA for improving plant drought adaptation.

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.

Biology and Function of miR159 in Plants

MicroR159 (miR159) is ancient, being present in the majority of land plants where it targets a class of regulatory genes called GAMYB or GAMYB-like via highly conserved miR159-binding sites. These GAMYB genes encode R2R3 MYB domain transcription factors that transduce the gibberellin (GA) signal in the seed aleurone and the anther tapetum. Here, GAMYB plays a conserved role in promoting the programmed cell death of these tissues, where miR159 function appears weak. By contrast, GAMYB is not involved in GA-signaling in vegetative tissues, but rather its expression is deleterious, leading to the inhibition of growth and development. Here, the major function of miR159 is to mediate strong silencing of GAMYB to enable normal growth. Highlighting this requirement of strong silencing are conserved RNA secondary structures associated with the miR159-binding site in GAMYB mRNA that promotes miR159-mediated repression. Although the miR159-GAMYB pathway in vegetative tissues has been implicated in a number of different functions, presently no conserved role for this pathway has emerged. We will review the current knowledge of the different proposed functions of miR159, and how this ancient pathway has been used as a model to help form our understanding of miRNA biology in plants.

The genome-wide impact of cadmium on microRNA and mRNA expression in contrasting Cd responsive wheat genotypes

BACKGROUND

Heavy metal ATPases (HMAs) are responsible for Cd translocation and play a primary role in Cd detoxification in various plant species. However, the characteristics of HMAs and the regulatory mechanisms between HMAs and microRNAs in wheat (Triticum aestivum L) remain unknown.

RESULTS

By comparative microRNA and transcriptome analysis, a total three known and 19 novel differentially expressed microRNAs (DEMs) and 1561 differentially expressed genes (DEGs) were found in L17 after Cd treatment. In H17, by contrast, 12 known and 57 novel DEMs, and only 297 Cd-induced DEGs were found. Functional enrichments of DEMs and DEGs indicate how genotype-specific biological processes responded to Cd stress. Processes found to be involved in microRNAs-associated Cd response include: ubiquitin mediated proteolysis, tyrosine metabolism, and carbon fixation pathways and thiamine metabolism. For the mRNA response, categories including terpenoid backbone biosynthesis and phenylalanine metabolism, and photosynthesis - antenna proteins and ABC transporters were enriched. Moreover, we identified 32 TaHMA genes in wheat. Phylogenetic trees, chromosomal locations, conserved motifs and expression levels in different tissues and roots under Cd stress are presented. Finally, we infer a microRNA-TaHMAs expression network, indicating that miRNAs can regulate TaHMAs.

CONCLUSION

Our findings suggest that microRNAs play important role in wheat under Cd stress through regulation of targets such as TaHMA2;1. Identification of these targets will be useful for screening and breeding low-Cd accumulation wheat lines.

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.

Expression of miR159 Is Altered in Tomato Plants Undergoing Drought Stress

In a scenario of global climate change, water scarcity is a major threat for agriculture, severely limiting crop yields. Therefore, alternatives are urgently needed for improving plant adaptation to drought stress. Among them, gene expression reprogramming by microRNAs (miRNAs) might offer a biotechnologically sound strategy. Drought-responsive miRNAs have been reported in many plant species, and some of them are known to participate in complex regulatory networks via their regulation of transcription factors involved in water stress signaling. We explored the role of miR159 in the response of Solanum lycopersicum Mill. plants to drought stress by analyzing the expression of sly-miR159 and its target SlMYB transcription factor genes in tomato plants of cv. Ailsa Craig grown in deprived water conditions or in response to mechanical damage caused by the Colorado potato beetle, a devastating insect pest of Solanaceae plants. Results showed that sly-miR159 regulatory function in the tomato plants response to distinct stresses might be mediated by differential stress-specific MYB transcription factor targeting. sly-miR159 targeting of SlMYB33 transcription factor transcript correlated with accumulation of the osmoprotective compounds proline and putrescine, which promote drought tolerance. This highlights the potential role of sly-miR159 in tomato plants’ adaptation to water deficit conditions.

Overexpression of SlGRAS7 Affects Multiple Behaviors Leading to Confer Abiotic Stresses Tolerance and Impacts Gibberellin and Auxin Signaling in Tomato

Abiotic stresses remain the key environmental issues that reduce plant development and therefore affect crop production. Transcription factors, such as the GRAS family, are involved in various functions of abiotic stresses and plant growth. The GRAS family of tomato (Solanum lycopersicum), SlGRAS7, is described in this study. We produced overexpressing SlGARS7 plants to learn more about the GRAS transcription factors. Plants overexpressing SlGARS7 (SlGRAS7-OE) showed multiple phenotypes related to many behaviors, including plant height, root and shoot length, and flowering time. We observed that many genes in the SlGRAS7-OE seedlings that are associated with auxin and gibberellin (GA) are downregulated and have altered sensitivity to GA₃/IAA. SlGRAS7 was upregulated during abiotic stresses following treatment with sodium chloride (NaCl) and D-mannitol in the wild-type (WT) tomato. Tomato plants overexpressing SlGRAS7 showed more resistance to drought and salt stress comparison with WT. Our study of SlGRAS7 in tomato demonstrates how GRAS showed an integrative role, improving resistance to abiotic stresses and enhancing gibberellin/auxin signaling through reproductive as well as vegetative processes.

Determination of the MiRNAs Related to Bean Pyralid Larvae Resistance in Soybean Using Small RNA and Transcriptome Sequencing

Soybean is one of the most important oil crops in the world. Bean pyralid is a major leaf-feeding insect of soybean. In order to screen out the functional genes and regulatory pathways related to the resistance for bean pyralid larvae, the small RNA and transcriptome sequencing were performed based on the highly resistant material (Gantai-2-2) and highly susceptible material (Wan 82-178) of soybean. The results showed that, when comparing 48 h feeding with 0 h feeding, 55 differentially expressed miRNAs were identified in Gantai-2-2 and 58 differentially expressed miRNAs were identified in Wan82-178. When comparing Gantai-2-2 with Wan82-178, 77 differentially expressed miRNAs were identified at 0 h feeding, and 70 differentially expressed miRNAs were identified at 48 h feeding. The pathway analysis of the predicted target genes revealed that the plant hormone signal transduction, RNA transport, protein processing in the endoplasmic reticulum, zeatin biosynthesis, ubiquinone and other terpenoid-quinone biosynthesis, and isoquinoline alkaloid biosynthesis may play important roles in soybean's defense against the stress caused by bean pyralid larvae. According to conjoint analysis of the miRNA/mRNA, a total of 20 differentially expressed miRNAs were negatively correlated with 26 differentially expressed target genes. The qRT-PCR analysis verified that the small RNA sequencing results were credible. According to the analyses of the differentially expressed miRNAs, we speculated that miRNAs are more likely to play key roles in the resistance to insects. Gma-miR156q, Gma-miR166u, Gma-miR166b, Gma-miR166j-3p, Gma-miR319d, Gma-miR394a-3p, Gma-miR396e, and so on-as well as their negatively regulated differentially expressed target genes-may be involved in the regulation of soybean resistance to bean pyralid larvae. These results laid a foundation for further in-depth research regarding the action mechanisms of insect resistance.

Characterization on the conservation and diversification of miRNA156 gene family from lower to higher plant species based on phylogenetic analysis at the whole genomic level

miRNA156 family members (miR156s) participate in regulating the transition of plant vegetative and reproductive growth, flower development, and formation of berry skin color by negatively modulating their target gene SPLs. However, the evolution and functional diversification of miR156s in plants remain elusive. Phylogenetic analysis on 310 miR156s from 51 plant species on miRBase 21.0 showed that only miR156a could be conserved in the 51 plant species, but their sequences exhibited variation; another set of miR156s, such as miR156m/n/o/p/q/r/s/t/u/v/w/x/y/z, was identified only in certain special plant species (Glycine max and Malus); also, all base variations in the sequences of 310 miR156s occurred within one miR156 seed sequence, "TGACAGAAGAGAGTGAGCAC," and the changed base sites were mainly located at the 11th and 14th bases from the 5' end of the miR156 seed sequence, in which some base variations of miR156s resulted in a difference in miR156 targeting modes; by contrast, miR156 precursor sequences are highly divergent across diverse species. Similarly, cis-regulatory motifs on the promoter sequence of MIR156s in various plants also exhibited significant discrepancy. The intragenic MIR156 genes overlapped their target SBP genes, thereby suggesting that some microRNAs (miRNAs) originate from duplication of target genes. These traits might be the reasons of the conservation and diversification of miR156 gene family. This study identified the conserved seed sequence "TGACAGAAGAGAGTGAGCAC," and the sequence variation characterization, of miR156 family evolution, also investigated the varied traits of their promoters, precursors, and mature sequences in sequence evolutions and found some miRNAs might originate from duplication of target genes. Our findings will contribute to our understanding of the functional diversification of miRNAs and the interactions of miRNA/target pairs based on the evolutionary history of miRNA genes.

Genome-Wide Analysis of the miRNA-mRNAs Network Involved in Cold Tolerance in Populus simonii × P. nigra

Background: Cold tolerance is important for plants’ geographical distribution and survival in extreme seasonal variations of climate. However, Populus simonii × P. nigra shows wide adaptability and strong cold resistance. Transcriptional and post-transcriptional regulation of cold-responsive genes is crucial for cold tolerance in plants. To understand the roles of regulatory RNAs under cold induction in Populus simonii × P. nigra, we constructed cDNA and small RNA libraries from leaf buds treated or not with −4 °C for 8 h for analysis. Results: Through high-throughput sequencing and differential expression analysis, 61 miRNAs and 1229 DEGs were identified under cold induction condition in Populus simonii × P. nigra. The result showed that miR167a, miR1450, miR319a, miR395b, miR393a-5p, miR408-5p, and miR168a-5p were downregulated, whereas transcription level of miR172a increased under the cold treatment. Thirty-one phased-siRNA were also obtained (reads ≥ 4) and some of them proceeded from TAS3 loci. Analysis of the differentially expressed genes (DEGs) showed that transcription factor genes such as Cluster-15451.2 (putative MYB), Cluster-16493.29872 (putative bZIP), Cluster-16493.29175 (putative SBP), and Cluster-1378.1 (putative ARF) were differentially expressed in cold treated and untreated plantlets of Populus simonii × P. nigra. Integrated analysis of miRNAs and transcriptome showed miR319, miR159, miR167, miR395, miR390, and miR172 and their target genes, including MYB, SBP, bZIP, ARF, LHW, and ATL, were predicted to be involved in ARF pathway, SPL pathway, DnaJ related photosystem II, and LRR receptor kinase, and many of them are known to resist chilling injury. Particularly, a sophisticated regulatory model including miRNAs, phasiRNAs, and targets of them was set up. Conclusions: Integrated analysis of miRNAs and transcriptome uncovered the complicated regulation of the tolerance of cold in Populus simonii × P. nigra. MiRNAs, phasiRNAs, and gene-encoded transcription factors were characterized at a whole genome level and their expression patterns were proved to be complementary. This work lays a foundation for further research of the pathway of sRNAs and regulatory factors involved in cold tolerance.

Integration of mRNA and miRNA Analysis Reveals the Molecular Mechanism Underlying Salt and Alkali Stress Tolerance in Tobacco

Salinity is one of the most severe forms of abiotic stress and affects crop yields worldwide. Plants respond to salinity stress via a sophisticated mechanism at the physiological, transcriptional and metabolic levels. However, the molecular regulatory networks involved in salt and alkali tolerance have not yet been elucidated. We developed an RNA-seq technique to perform mRNA and small RNA (sRNA) sequencing of plants under salt (NaCl) and alkali (NaHCO3) stress in tobacco. Overall, 8064 differentially expressed genes (DEGs) and 33 differentially expressed microRNAs (DE miRNAs) were identified in response to salt and alkali stress. A total of 1578 overlapping DEGs, which exhibit the same expression patterns and are involved in ion channel, aquaporin (AQP) and antioxidant activities, were identified. Furthermore, genes involved in several biological processes, such as "photosynthesis" and "starch and sucrose metabolism," were specifically enriched under NaHCO3 treatment. We also identified 15 and 22 miRNAs that were differentially expressed in response to NaCl and NaHCO3, respectively. Analysis of inverse correlations between miRNAs and target mRNAs revealed 26 mRNA-miRNA interactions under NaCl treatment and 139 mRNA-miRNA interactions under NaHCO3 treatment. This study provides new insights into the molecular mechanisms underlying the response of tobacco to salinity stress.

Stress-responsive miRNAome of Glycine max (L.) Merrill: molecular insights and way forward

MAIN CONCLUSION

Analysis of stress-associated miRNAs of Glycine max (L.) Merrill reveals wider ramifications of small RNA-mediated (conserved and legume-specific miRNAs) gene regulatory foot prints in molecular adaptive responses. MicroRNAs (miRNAs) are indispensable components of gene regulatory mechanism of plants. Soybean is a crop of immense commercial potential grown worldwide for its edible oil and soy meal. Intensive research efforts, using the next generation sequencing and bioinformatics techniques, have led to the identification and characterization of numerous small RNAs, especially microRNAs (miRNAs), in soybean. Furthermore, studies have unequivocally demonstrated the significance of miRNAs during the developmental processes and various stresses in soybean. In this review, we summarize the current state of understanding of miRNA-based abiotic and biotic stress responses in soybean. In addition, the molecular insights gained from the stress-related soybean miRNAs have been compared to the miRNAs of other crops, especially legumes, and the core commonalities have been highlighted, though differences among them were not ignored. Nature of response of soybean-derived conserved miRNAs during various stresses was also analyzed to gain deeper insights regarding sRNAome-based defense responses. This review further provides way forward in legume small RNA transcriptomics based on the adaptive responses of soybean and other legume-derived miRNAs.

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

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

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

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

Meta-Analysis of Salt Stress Transcriptome Responses in Different Rice Genotypes at the Seedling Stage

Rice (Oryza sativa L.) is one of the most important staple food crops worldwide, while its growth and productivity are threatened by various abiotic stresses, especially salt stress. Unraveling how rice adapts to salt stress at the transcription level is vital. It can provide valuable information on enhancing the salt stress tolerance performance of rice via genetic engineering technologies. Here, we conducted a meta-analysis of different rice genotypes at the seedling stage based on 96 public microarray datasets, aiming to identify the key salt-responsive genes and understand the molecular response mechanism of rice under salt stress. In total, 5559 genes were identified to be differentially expressed genes (DEGs) under salt stress, and 3210 DEGs were identified during the recovery process. The Gene Ontology (GO) enrichment results revealed that the salt-response mechanisms of shoots and roots were different. A close-knit signaling network, consisting of the Ca²⁺ signal transduction pathway, the mitogen-activated protein kinase (MAPK) cascade, multiple hormone signals, transcription factors (TFs), transcriptional regulators (TRs), protein kinases (PKs), and other crucial functional proteins, plays an essential role in rice salt stress response. In this study, many unreported salt-responsive genes were found. Besides this, MapMan results suggested that TNG67 can shift to the fermentation pathway to produce energy under salt stress and may enhance the Calvin cycle to repair a damaged photosystem during the recovery stage. Taken together, these findings provide novel insights into the salt stress molecular response and introduce numerous candidate genes for rice salt stress tolerance breeding.

Homeostasis in the soybean miRNA396-GRF network is essential for productive soybean cyst nematode infections

Heterodera glycines, the soybean cyst nematode, penetrates soybean roots and migrates to the vascular cylinder where it forms a feeding site called the syncytium. MiRNA396 (miR396) targets growth-regulating factor (GRF) genes, and the miR396-GRF1/3 module is a master regulator of syncytium development in model cyst nematode H. schachtii infection of Arabidopsis. Here, we investigated whether this regulatory system operates similarly in soybean roots and is likewise important for H. glycines infection. We found that a network involving nine MIR396 and 23 GRF genes is important for normal development of soybean roots and that GRF function is specified in the root apical meristem by miR396. All MIR396 genes are down-regulated in the syncytium during its formation phase while, specifically, 11 different GRF genes are up-regulated. The switch to the syncytium maintenance phase coincides with up-regulation of MIR396 and down-regulation of the 11 GRF genes specifically via post-transcriptional regulation by miR396. Furthermore, interference with the miR396-GRF6/8-13/15-17/19 regulatory network, through either overexpression or knockdown experiments, does not affect the number of H. glycines juveniles that enter the vascular cylinder to initiate syncytia, but specifically inhibits efficient H. glycines development to adult females. Therefore, homeostasis in the miR396-GRF6/8-13/15-17/19 regulatory network is essential for productive H. glycines infections.

Profiling the Abiotic Stress Responsive microRNA Landscape of Arabidopsis thaliana

It is well established among interdisciplinary researchers that there is an urgent need to address the negative impacts that accompany climate change. One such negative impact is the increased prevalence of unfavorable environmental conditions that significantly contribute to reduced agricultural yield. Plant microRNAs (miRNAs) are key gene expression regulators that control development, defense against invading pathogens and adaptation to abiotic stress. Arabidopsis thaliana (Arabidopsis) can be readily molecularly manipulated, therefore offering an excellent experimental system to alter the profile of abiotic stress responsive miRNA/target gene expression modules to determine whether such modification enables Arabidopsis to express an altered abiotic stress response phenotype. Towards this goal, high throughput sequencing was used to profile the miRNA landscape of Arabidopsis whole seedlings exposed to heat, drought and salt stress, and identified 121, 123 and 118 miRNAs with a greater than 2-fold altered abundance, respectively. Quantitative reverse transcriptase polymerase chain reaction (RT-qPCR) was next employed to experimentally validate miRNA abundance fold changes, and to document reciprocal expression trends for the target genes of miRNAs determined abiotic stress responsive. RT-qPCR also demonstrated that each miRNA/target gene expression module determined to be abiotic stress responsive in Arabidopsis whole seedlings was reflective of altered miRNA/target gene abundance in Arabidopsis root and shoot tissues post salt stress exposure. Taken together, the data presented here offers an excellent starting platform to identify the miRNA/target gene expression modules for future molecular manipulation to generate plant lines that display an altered response phenotype to abiotic stress.

Genome-wide identification of CAMTA gene family members in Phaseolus vulgaris L. and their expression profiling during salt stress

The calmodulin-binding transcriptional activator (CAMTA) family was first observed in tobacco (NtER1) during a screening for the CaM-binding proteins, which are known to be one of the fast response stress proteins. Due to the increased importance of plant transcription factors in recent years; genome-wide identification of CAMTA genes has been performed in several plant species, except for Phaseolus vulgaris. Therefore, our aim was to identify and characterize CAMTA genes in P. vulgaris via in silico genome-wide analysis approach. Our results showed a total of eight CAMTA genes that were identified and observed on five out of 11 chromosomes of P. vulgaris. Four gene couples were found to be segmentally-duplicated and these segmental duplication events were shown to occur from 29.97 to 92.06 MYA. The phylogenetic tree of CAMTA homologs from P. vulgaris, A. thaliana, and G. max. revealed three groups based on their homology and the intron numbers of Pvul-CAMTA genes, ranged from 11 to 12. According to the syteny analysis; CAMTA genes of P. vulgaris and G. max revealed higher similarity, because they have highly similar genomes compared to A. thaliana. All Pvul-CAMTA genes were targeted by miRNAs, which play a role in response mechanism of salt stress. To detect expression levels in different plant tissues, mRNA analysis of Pvul-CAMTA genes were performed using publicly available expression data in Phytozome v12.1. In addition, responses of Pvul-CAMTA genes to salt stress, were also examined via both RNAseq and qRT-PCR analysis. To identify and to obtain insight into biological functions of CAMTA genes in the genome of P. vulgaris, several analyses were conducted using many online and offline bioinformatic tools, genome databases and qRT-PCR analyses. Due to this study being the first in the identification of CAMTA genes in P. vulgaris, this study could be considered as an useful source for future CAMTA genes studies in either P. vulgaris or comparative different plant species.

Modulation of miRNA expression in natural populations of A. thaliana along a wide altitudinal gradient of Indian Himalayas

Plant populations growing along an altitudinal gradient are exposed to different environmental conditions. They are excellent resources to study regulatory mechanisms adopted by plants to respond to different environmental stresses. Regulation by miRNA is one of such strategies. Here, we report how different miRNAs are preferentially expressed in the three natural populations of A. thaliana originating from a wide altitudinal range. The expression level of miRNAs was mostly governed by temperature and radiation. Majority of the identified miRNAs expressed commonly in the three populations. However, 30 miRNAs expressed significantly at different level between the low and the high altitude populations. Most of these miRNAs regulate the genes associated with different developmental processes, abiotic stresses including UV, cold, secondary metabolites, etc. Further, the expression of miR397 and miR858 involved in lignin biosynthesis and regulation of secondary metabolites respectively, may be regulated by light intensity. A few miRNAs expressed at increasing level with the increase in the altitude of the site indicating environment driven tight regulation of these miRNAs. Further, several novel miRNAs and isomiR diversity specific to the Himalayas are reported which might have an adaptive advantage. To the best of our knowledge, this is the first report on miRNA expression from natural plant populations.

Banana sRNAome and degradome identify microRNAs functioning in differential responses to temperature stress

Background

Temperature stress is a major environmental factor affecting not only plant growth and development, but also fruit postharvest life and quality. MicroRNAs (miRNAs) are a class of non-coding small RNAs that play important roles in various biological processes. Harvested banana fruit can exhibit distinct symptoms in response to different temperature stresses, but the underlying miRNA-mediated regulatory mechanisms remained unknown.

Results

Here, we profiled temperature-responsive miRNAs in banana, using deep sequencing and computational and molecular analyses. In total 113 known miRNAs and 26 novel banana-specific miRNAs were identified. Of these miRNAs, 42 miRNAs were expressed differentially under cold and heat stresses. Degradome sequencing identified 60 target genes regulated by known miRNAs and half of these targets were regulated by 15 temperature-responsive miRNAs. The correlative expression patterns between several miRNAs and their target genes were further validated via qRT-PCR. Our data showed that miR535 and miR156 families may derive from a common ancestor during evolution and jointly play a role in fine-tuning SPL gene expression in banana. We also identified the miRNA-triggered phased secondary siRNAs in banana and found miR393-TIR1/AFB phasiRNA production displaying cold stress-specific enrichment.

Conclusions

Our results provide a foundation for understanding the miRNA-dependent temperature stress response in banana. The characterized correlations between miRNAs and their response to temperature stress could serve as markers in the breeding programs or tools for improving temperature tolerance of banana.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5395-1) contains supplementary material, which is available to authorized users.

Functional Characterization of microRNA171 Family in Tomato

Deeply conserved plant microRNAs (miRNAs) function as pivotal regulators of development. Nevertheless, in the model crop Solanum lycopersicum (tomato) several conserved miRNAs are still poorly annotated and knowledge about their functions is lacking. Here, the tomato miR171 family was functionally analyzed. We found that the tomato genome contains at least 11 SlMIR171 genes that are differentially expressed along tomato development. Downregulation of sly-miR171 in tomato was successfully achieved by transgenic expression of a short tandem target mimic construct (STTM171). Consequently, sly-miR171-targeted mRNAs were upregulated in the silenced plants. Target upregulation was associated with irregular compound leaf development and an increase in the number of axillary branches. A prominent phenotype of STTM171 expressing plants was their male sterility due to a production of a low number of malformed and nonviable pollen. We showed that sly-miR171 was expressed in anthers along microsporogenesis and significantly silenced upon STTM171 expression. Sly-miR171-silenced anthers showed delayed tapetum ontogenesis and reduced callose deposition around the tetrads, both of which together or separately can impair pollen development. Collectively, our results show that sly-miR171 is involved in the regulation of anther development as well as shoot branching and compound leaf morphogenesis.

Cre-miR914-regulated RPL18 is involved with UV-B adaptation in Chlamydomonas reinhardtii

UV radiation is a serious threat to life, and algae have developed highly efficient adaptations to UV radiation through the course of evolution. To date, studies investigating the mechanisms of UV adaptation in algae have focused on physiological regulation and associated protein coding genes, with only a few reports on associated protein non-coding genes. In a previous study, we found that Cre-miR914 was significantly down-regulated in Chlamydomonas reinhardtii in response to heat shock. In the present study, we aimed to determine whether Cre-miR914 plays a role in response to UV-B radiation. Our bioinformatics analysis indicated that the potential target gene of Cre-miR914 is ribosomal protein L18 (RPL18). We also measured the expression of Cre-miR914 and RPL18 in response to UV-B radiation through qPCR analysis. Then, we constructed cell lines overexpressing Cre-miR914 or RPL18, and performed survival experiments under UV-B stress. The results showed that Cre-miR914 overexpression decreased resistance while RPL18 overexpression enhanced tolerance to UV-B radiation. These results indicate that Cre-miR914 and its potential target gene RPL18 are involved in the adaptation to UV-B in C. reinhardtii.

miR393s regulate salt stress response pathway in Arabidopsis thaliana through scaffold protein RACK1A mediated ABA signaling pathways

Scaffold protein Receptor for Activated C Kinase 1 (RACK1) is a negative regulator of plant stress hormone - abscisic acid (ABA) mediated pathways. RACK1 has been reported to regulate global miRNA biogenesis pathway in C. elegans, humans, and in Arabidopsis. RACK1 regulates different steps of miRNA biogenesis and stability in response to different stimuli in plants. miR393s is implicated in salt stress response pathway through an antagonistic response between the stress hormone ABA-mediated salt stress and growth hormone auxin. Specifically, the known auxin receptor clade transcripts TIR1/AFB2 are the target for the miR393s. By down-regulating the auxin signaling pathways, the miR393s inhibit the regulation of salt tolerance by auxin. Here we show that genetic loss of RACK1A- the predominant member of the three genes family of RACK1 in Arabidopsis, results in the inhibition of miR393 level causing the same salt sensitivities as the individual mir393a or mir393b or the double mutant mir393ab phenotypes. We propose that down-regulation of auxin signaling through RACK1A induced miR393 biogenesis potentially regulates the Arabidopsis acclimation to salinity. Our findings fill up a molecular gap in our understanding of the role of miR393 mediated ABA and auxin-regulated salt stress responses.

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.

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

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

Regulation through MicroRNAs in Response to Low-Energy N+ Ion Irradiation in Oryza sativa

MicroRNAs (miRNAs) are a non-coding regulatory RNAs that play significant roles in plant growth and development, especially in the stress response. Low-energy ion radiation, a type of environmental stress, can cause multiple biological effects. To understand the roles of miRNAs in response to low-energy N⁺ ion radiation in Oryza sativa, high-throughput sequencing of small RNAs was carried out to detect the expression of miRNAs in the shoots of the rice after 2 × 10¹⁷ N+/cm² irradiation. The differentially expressed 28 known miRNAs were identified, 17 of these identified miRNAs were validated by real-time quantitative fluorescent PCR (q-PCR), including 9 up-regulated miRNAs (miR1320-3p, miR1320-5p, miR156b-3p, miR156c-5p, miR156c-3p/g-3p, miR1561-5p, miR398b and miR6250) and 8 down-regulated miRNAs (miR156a/e/i, miR156k, miR160f-5p, miR166j-5p, miR1846e and miR399d). In addition, 45 novel radiationresponsive miRNAs were predicted, and 8 of them were verified by q-PCR. The target genes of radiation-responsive miRNAs were predicted and gene function enrichment analysis was performed with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). The expression of 9 targets of 4 known miRNA families (miR156, miR399, miR1320 and miR398) and 2 targets of 2 novel miRNAs were quantified by q-PCR, and a strong negative regulation relation between miRNAs and their targets were observed. Those targets including SQUAMOSA promoterbinding-like protein (SPL) genes, copper/zinc superoxide dismutase (Cu/Zn-SOD), copper chaperone for SOD (CCS1) and electron transporter/ heat-shock protein binding protein (HSP), which are involved in growth and defense against various stresses, especially associated with reactive oxygen species (ROS) scavenging. This work provides important information for understanding the ROS generation and elimination mechanisms closely related to miRNAs in rice seedlings after low-energy N⁺ radiation exposure.

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

BACKGROUND

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

RESULT

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

CONCLUSIONS

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

Physiological studies and genome-wide microRNA profiling of cold-stressed Brassica napus

Temperature extremes, including cold, adversely impact plant growth and development. Plant responses to cold stress (CS) are regulated at both transcriptional and post-transcriptional levels. MicroRNAs (miRNAs), small non-coding RNAs, are known to be involved in post-transcriptional regulation of various developmental processes and metal stress in Brassica napus L. (canola), however, their role in response to CS is largely unknown. In this study, changes in various physiological parameters and endogenous abundance of miRNAs were characterized in spring canola seedlings (DH12075) exposed to 4 °C for 0-48 h. Cold stress induced electrolyte leakage, increased the levels of malondialdheyde and antioxidant enzymes and reduced photosynthetic efficiency. Using small RNA sequencing, 70 known and 126 novel miRNAs were identified in CS leaf tissues and among these, 25 known and 104 novel miRNAs were differentially expressed. Quantitative real-time (qRT) PCR analysis of eight selected miRNAs confirmed their CS responsiveness. Furthermore, the expression of six out of eight miRNAs exhibited an opposite trend in a winter variety of canola, 'Mendel', when compared to 'DH12075'. This first study on the B. napus miRNAome provides a framework for further functional analysis of these miRNAs and their targets in response to CS which may contribute towards the future development of cold resilient crops.

Transcriptome and miRNAs analyses enhance our understanding of the evolutionary advantages of polyploidy

Polyploid organisms have more than two sets of chromosomes, including autopolyploid via intraspecific genome doubling, and allopolyploid via merging genomes of distinct species by hybridization. Polyploid organisms are widespread in plants, indicating that polyploidy has some evolutionary advantages over its diploid ancestor. Actually, polyploidy is always tightly associated with hybrid vigor and adaptation to adverse environmental conditions. However, why polyploidy can develop such advantages is poorly known. MicroRNAs (miRNAs) are endogenous ∼21 nt small RNAs which can play important regulatory roles in animals and plants by targeting mRNAs for cleavage or translational repression. MicroRNAs are essential for cell development, differentiation, signal transduction, and show an adaptive response to biotic and abiotic stresses. Environmental stresses cause plants to over- or under-express certain miRNAs or synthesize new miRNAs to cope with stress. We have here reviewed our current knowledge on the molecular mechanisms, which can account for the evolutionary advantages of polyploidy over its diploid ancestor from genome-wide gene expression and microRNAs expression perspectives.

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.

Elevated carbon dioxide and drought modulate physiology and storage-root development in sweet potato by regulating microRNAs

Elevated CO₂ along with drought is a serious global threat to crop productivity. Therefore, understanding the molecular mechanisms plants use to protect these stresses is the key for plant growth and development. In this study, we mimicked natural stress conditions under a controlled Soil-Plant-Atmosphere-Research (SPAR) system and provided the evidence for how miRNAs regulate target genes under elevated CO₂ and drought conditions. Significant physiological and biomass data supported the effective utilization of source-sink (leaf to root) under elevated CO₂. Additionally, elevated CO₂ partially rescued the effect of drought on total biomass. We identified both known and novel miRNAs differentially expressed during drought, CO₂, and combined stress, along with putative targets. A total of 32 conserved miRNAs belonged to 23 miRNA families, and 25 novel miRNAs were identified by deep sequencing. Using the existing sweet potato genome database and stringent analyses, a total of 42 and 22 potential target genes were predicted for the conserved and novel miRNAs, respectively. These target genes are involved in drought response, hormone signaling, photosynthesis, carbon fixation, sucrose and starch metabolism, etc. Gene ontology and KEGG ontology functional enrichment revealed that these miRNAs might target transcription factors (MYB, TCP, NAC), hormone signaling regulators (ARF, AP2/ERF), cold and drought factors (corA), carbon metabolism (ATP synthase, fructose-1,6-bisphosphate), and photosynthesis (photosystem I and II complex units). Our study is the first report identifying targets of miRNAs under elevated CO₂ levels and could support the molecular mechanisms under elevated CO₂ in sweet potato and other crops in the future.

Microsynteny analysis to understand evolution and impact of polyploidization on MIR319 family within Brassicaceae

The availability of a large number of whole-genome sequences allows comparative genomic analysis to reveal and understand evolution of regulatory regions and elements. The role played by events such as whole-genome and segmental duplications followed by genome fractionation in shaping genomic landscape and in expansion of gene families is crucial toward developing insights into evolutionary trends and consequences such as sequence and functional diversification. Members of Brassicaceae are known to have experienced several rounds of whole-genome duplication (WGD) that have been termed as paleopolyploidy, mesopolyploidy, and neopolyploidy. Such repeated events led to the creation and expansion of a large number of gene families. MIR319 is reported to be one of the most ancient and conserved plant MIRNA families and plays a role in growth and development including leaf development, seedling development, and embryo patterning. We have previously reported functional diversification of members of miR319 in Brassica oleracea affecting leaf architecture; however, the evolutionary history of the MIR319 gene family across Brassicaceae remains unknown and requires investigation. We therefore identified homologous and homeologous segments of ca. 100 kb, with or without MIR319, performed comparative synteny analysis and genome fractionation studies. We detected variable rates of gene retention across members of Brassicaceae when genomic blocks of MIR319a, MIR319b, and MIR319c were compared either between themselves or against Arabidopsis thaliana genome which was taken as the base genome. The highest levels of shared genes were found between A. thaliana and Capsella rubella in both MIR319b- and MIR319c-containing genomic segments, and with the closest species of A. thaliana, A. lyrata, only in MIR319a-containing segment. Synteny analysis across 12 genomes (with 30 sub-genomes) revealed MIR319c to be the most conserved MIRNA loci (present in 27 genomes/sub-genomes) followed by MIR319a (present in 23 genomes/sub-genomes); MIR319b was found to be frequently lost (present in 20 genomes/sub-genomes) and thus is under least selection pressure for retention. Genome fractionation revealed extensive and differential loss of MIRNA homeologous loci and flanking genes from various sub-genomes of Brassica species that is in accordance with their older history of polyploidy when compared to Camelina sativa, a recent neopolyploid, where the effect of genome fractionation was least. Finally, estimation of phylogenetic relationship using precursor sequences of MIR319 reveals MIR319a and MIR319b form sister clades, with MIR319c forming a separate clade. An intra-species synteny analysis between MIR319a-, MIR319b-, and MIR319c-containing genomic segments suggests segmental duplications at the base of Brassicaceae to be responsible for the origin of MIR319a and MIR319b.

Uncovering key small RNAs associated with gametocidal action in wheat

Gametocidal (Gc) chromosomes can kill gametes that lack them by causing chromosomal breakage to ensure their preferential transmission, and they have been exploited in genetic breeding. The present study investigated the possible roles of small RNAs (sRNAs) in Gc action. By sequencing two small RNA libraries from the anthers of Triticum aestivum cv. Chinese Spring (CS) and the Chinese Spring-Gc 3C chromosome monosomic addition line (CS-3C), we identified 239 conserved and 72 putative novel miRNAs, including 135 differentially expressed miRNAs. These miRNAs were predicted to target multiple genes with various molecular functions relevant to the features of Gc action, including sterility and genome instability. The transgenic overexpression of miRNA, which was up-regulated in CS-3C, reduced rice fertility. The CS-3C line exhibited a genome-wide reduction in 24 nt siRNAs compared with that of the CS line, particularly in transposable element (TE) and repetitive DNA sequences. Corresponding to this reduction, the bisulfite sequencing analysis of four retro-TE sequences showed a decrease in CHH methylation, typical of RNA-directed DNA methylation (RdDM). These results demonstrate that both miRNA-directed regulation of gene expression and siRNA-directed DNA methylation of target TE loci could play a role in Gc action.

miRNAs play important roles in aroma weakening during the shelf life of 'Nanguo' pear after cold storage

Cold storage is commonly employed to delay senescence in 'Nanguo' pear after harvest. However, this technique also causes fruit aroma weakening. MicroRNAs are regulators of gene expression at the post-transcriptional level that play important roles in plant development and in eliciting responses to abiotic environmental stressors. In this study, the miRNA transcript profiles of the fruit on the first day (C0, LT0) move in and out of cold storage and the optimum tasting period (COTP, LTOTP) during shelf life at room temperature and after cold storage were analyzed, respectively. 314 known miRNAs were identified in 'Nanguo' pear; 176 and 135 miRNAs were significantly differentially expressed on the C0 vs. LT0 and on the COTP vs. LTOTP, respectively. After prediction the target genes of these differentially expressed miRNAs, 9 s-lipoxygenase (LOX2S), 13 s-lipoxygenase (LOX1_5), hydroperoxide lyase (HPL), and alcohol dehydrogenase (ADH1) were found differentially expressed, which were the key genes during aroma formation. The expression pattern of these target genes and the related miRNAs were identified by RT-PCR. mdm-miR172a-h, mdm-miR159a/b, mdm-miR160a-e, mdm-miR395a-i, mdm-miR399a/b/c, mdm/ppe-miR535a/b, and mdm-miR7120a/b may negatively regulate the target genes expression. These results indicate that miRNAs may play key roles in aroma weakening in cold storage 'Nanguo' pear and provide valuable information for studying the molecular mechanisms of miRNAs in the aroma weakening of fruit due to low temperature.

Genome-Wide Investigation of the Role of MicroRNAs in Desiccation Tolerance in the Resurrection Grass Tripogon loliiformis

Drought causes approximately two-thirds of crop and yield loss worldwide. To sustain future generations, there is a need to develop robust crops with enhanced water use efficiency. Resurrection plants are naturally resilient and tolerate up to 95% water loss with the ability to revive upon watering. Stress is genetically encoded and resilient species may garner tolerance by tightly regulating the expression of stress-related genes. MicroRNAs (miRNAs) post-transcriptionally regulate development and other stress response processes in eukaryotes. However, their role in resurrection plant desiccation tolerance is poorly understood. In this study, small RNA sequencing and miRNA expression profiling was conducted using Tripogon loliiformis plants subjected to extreme water deficit conditions. Differentially expressed miRNA profiles, target mRNAs, and their regulatory processes were elucidated. Gene ontology enrichment analysis revealed that development, stress response, and regulation of programmed cell death biological processes; Oxidoreductase and hydrolyase molecular activities; and SPL, MYB, and WRKY transcription factors were targeted by miRNAs during dehydration stress, indicating the indispensable regulatory role of miRNAs in desiccation tolerance. This study provides insights into the molecular mechanisms of desiccation tolerance in the resurrection plant T. loliiformis. This information will be useful in devising strategies for crop improvement on enhanced drought tolerance and water use efficiency.

Repetitive somatic embryogenesis induced cytological and proteomic changes in embryogenic lines of Pseudotsuga menziesii [Mirb.]

BACKGROUND

To explore poorly understood differences between primary and subsequent somatic embryogenic lines of plants, we induced secondary (2ry) and tertiary (3ry) lines from cotyledonary somatic embryos (SEs) of two Douglas-fir genotypes: SD4 and TD17. The 2ry lines exhibited significantly higher embryogenic potential (SE yields) than the 1ry lines initiated from zygotic embryos (SD4, 2155 vs 477; TD17, 240 vs 29 g- 1 f.w.). Moreover, we observed similar differences in yield between 2ry and 3ry lines of SD4 (2400 vs 3921 g- 1 f.w.). To elucidate reasons for differences in embryogenic potential induced by repetitive somatic embryogenesis we then compared 2ry vs 1ry and 2ry vs 3ry lines at histo-cytological (using LC-MS/MS) and proteomic levels.

RESULTS

Repetitive somatic embryogenesis dramatically improved the proliferating lines' cellular organization (genotype SD4's most strongly). Frequencies of singulated, bipolar SEs and compact polyembryogenic centers with elongated suspensors and apparently cleavable embryonal heads increased in 2ry and (even more) 3ry lines. Among 2300-2500 identified proteins, 162 and 228 were classified significantly differentially expressed between 2ry vs 1ry and 3ry vs 2ry lines, respectively, with special emphasis on "Proteolysis" and "Catabolic process" Gene Ontology categories. Strikingly, most of the significant proteins (> 70%) were down-regulated in 2ry relative to 1ry lines, but up-regulated in 3ry relative to 2ry lines, revealing a down-up pattern of expression. GO category enrichment analyses highlighted the opposite adjustments of global protein patterns, particularly for processes involved in chitin catabolism, lignin and L-phenylalanine metabolism, phenylpropanoid biosynthesis, oxidation-reduction, and response to karrikin. Sub-Network Enrichment Analyses highlighted interactions between significant proteins and both plant growth regulators and secondary metabolites after first (especially jasmonic acid, flavonoids) and second (especially salicylic acid, abscisic acid, lignin) embryogenesis cycles. Protein networks established after each induction affected the same "Plant development" and "Defense response" biological processes, but most strongly after the third cycle, which could explain the top embryogenic performance of 3ry lines.

CONCLUSIONS

This first report of cellular and molecular changes after repetitive somatic embryogenesis in conifers shows that each cycle enhanced the structure and singularization of EMs through modulation of growth regulator pathways, thereby improving the lines' embryogenic status.

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

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

Evolutionary and expression analysis of CAMTA gene family in Nicotiana tabacum yielded insights into their origin, expansion and stress responses

Calmodulin-binding transcription activators (CAMTAs) represent the novel gene family of transcriptional regulators, which play important biological functions. Though, the first ever plant CAMTA gene was evidenced in Nicotiana tabacum in 2002. But, the systematic identification, origin and function of this gene family has not been performed due to the lack of reference genome information until now. Here, we identified 29 CAMTA genes in four Nicotiana species, including thirteen NtabCAMTAs, six NsylCAMTAs, and five NtomCAMTAs and NbenCAMTAs. These CAMTA families were classified into five phylogenetic groups (I-V), among which, the group-IV CAMTAs probably emerged the earliest. The NtabCAMTA family genes have diverse structures, and are randomly localized on five chromosomes and scaffolds. N. tabacum acquired 11 copies of homolog CAMATA genes from the parental genomes of N. tomentosiformis and N. sylvestris, followed by expansion through polyploidization and duplication. The NtabCAMTA genes were differentially expressed in different plant parts, and showed sensitivity towards different abiotic and biotic stresses. Co-expression network analysis revealed that some NtabCAMTA subunits interact with each other, and co-expressed. The current study is the first report presenting a comprehensive overview of Nicotiana CAMTA families, and opens a new avenue for the improvement of the cultivated tobacco.

Identification and characterisation of a previously unknown drought tolerance-associated microRNA in barley

Drought is the most serious abiotic stress, and causes crop losses on a worldwide scale. The present study identified a previously unknown microRNA (designated as hvu-miRX) of 21 nucleotides (nt) in length in barley. Its precursor (designated pre-miRX) and primary transcript (designated pri-miRX) were also identified, with lengths of 73 and 559 nt, respectively. The identified upstream sequence of pri-miRX contained both the TATA box and the CAAT box, which are both required for initiation of transcription. Transient promoter activation assays showed that the core promoter region of pri-miRX ranged 500 nt from the transcription start site. In transgenic barley overexpression of the wheat DREB3 transcription factor (TaDREB3) caused hvu-miRX to be highly expressed as compared with the same miRNA in non-transgenic barley. However, the high expression was not directly associated with TaDREB3. Genomic analysis revealed that the hvu-miRX gene was a single copy located on the short arm of chromosome 2 and appeared to be only conserved in Triticeae, but not in other plant species. Notably, transgenic barley that overexpressed hvu-miRX showed drought tolerance. Degradome library analysis and other tests showed that hvu-miRX targeted various genes including transcription factors via the cleavage mode. Our data provides an excellent opportunity to develop drought stress tolerant cereals using hvu-miRX.