Identification of chilling-responsive microRNAs and their targets in vegetable soybean (Glycine max L.)

Progress and Prospects of the Molecular Basis of Soybean Cold Tolerance

Cold stress is a major factor influencing the geographical distribution of soybean growth and causes immense losses in productivity. Understanding the molecular mechanisms that the soybean has undergone to survive cold temperatures will have immense value in improving soybean cold tolerance. This review focuses on the molecular mechanisms involved in soybean response to cold. We summarized the recent studies on soybean cold-tolerant quantitative trait loci (QTLs), transcription factors, associated cold-regulated (COR) genes, and the regulatory pathways in response to cold stress. Cold-tolerant QTLs were found to be overlapped with the genomic region of maturity loci of E1, E3, E4, pubescence color locus of T, stem growth habit gene locus of Dt1, and leaf shape locus of Ln, indicating that pleiotropic loci may control multiple traits, including cold tolerance. The C-repeat responsive element binding factors (CBFs) are evolutionarily conserved across species. The expression of most GmDREB1s was upregulated by cold stress and overexpression of GmDREB1B;1 in soybean protoplast, and transgenic Arabidopsis plants can increase the expression of genes with the DRE core motif in their promoter regions under cold stress. Other soybean cold-responsive regulators, such as GmMYBJ1, GmNEK1, GmZF1, GmbZIP, GmTCF1a, SCOF-1 and so on, enhance cold tolerance by regulating the expression of COR genes in transgenic Arabidopsis. CBF-dependent and CBF-independent pathways are cross-talking and work together to activate cold stress gene expression. Even though it requires further dissection for precise understanding, the function of soybean cold-responsive transcription factors and associated COR genes studied in Arabidopsis shed light on the molecular mechanism of cold responses in soybeans and other crops. Furthermore, the findings may also provide practical applications for breeding cold-tolerant soybean varieties in high-latitude and high-altitude regions.

Research progress on plant noncoding RNAs in response to low-temperature stress

Low temperature (LT) is an important factor limiting plant growth and distribution. Plants have evolved sophisticated adaptive mechanisms to cope with hypothermia. RNA silencing is the orchestrator of these cellular responses. RNA silencing, which modifies gene expression through noncoding RNAs (ncRNAs), is a strategy used by plants to combat environmental stress. ncRNAs, which have very little protein-coding capacity, work by binding reverse complementary endogenous transcripts. In plants, ncRNAs include small non-coding RNAs (sncRNAs), medium-sized non-coding RNAs (mncRNAs), and long non-coding RNAs (lncRNAs). Apart from describing the biogenesis of different ncRNAs (miRNAs, siRNAs, and lncRNAs), we thoroughly discuss the functions of these ncRNAs during cold acclimation. Two major classes of sncRNAs, microRNAs and siRNAs, play essential regulatory roles in cold response processes through the posttranscriptional gene silencing (PTGS) pathway or transcriptional gene silencing (TGS) pathway. Microarray or transcriptome sequencing analysis can reveal a large number of cold-responsive miRNAs in plants. In this review, the cold-response patterns of miRNAs verified by Northern blotting or quantitative PCR in Arabidopsis thaliana, rice, and many other important crops are discussed. The detailed molecular mechanisms of several miRNAs in Arabidopsis (miR397, miR408, miR402, and miR394) and rice (Osa-miR156, Osa-miR319, and Osa-miR528) that regulate plant cold resistance are elucidated. In addition, the regulatory mechanism of the lncRNA SVALKA in the cold signaling pathway is explained in detail. Finally, we present the challenges for understanding the roles of small ncRNAs in cold signal transduction.

Promoter of Vegetable Soybean GmTIP1;6 Responds to Diverse Abiotic Stresses and Hormone Signals in Transgenic Arabidopsis

Tonoplast intrinsic proteins (TIPs), a sub-family of aquaporins (AQPs), are known to play important roles in plant abiotic stress responses. However, evidence for the promoters of TIPs involvement in abiotic stress processes remains scarce. In this study, the promoter of the vegetable soybean GmTIP1;6 gene, which had the highest similarity to TIP1-type AQPs from other plants, was cloned. Expression pattern analyses indicated that the GmTIP1;6 gene was dramatically induced by drought, salt, abscisic acid (ABA), and methyl jasmonate (MeJA) stimuli. Promoter analyses revealed that the GmTIP1;6 promoter contained drought, ABA, and MeJA cis-acting elements. Histochemical staining of the GmTIP1;6 promoter in transgenic Arabidopsis corroborated that it was strongly expressed in the vascular bundles of leaves, stems, and roots. Beta-glucuronidase (GUS) activity assays showed that the activities of the GmTIP1;6 promoter were enhanced by different concentrations of polyethylene glycol 6000 (PEG 6000), NaCl, ABA, and MEJA treatments. Integrating these results revealed that the GmTIP1;6 promoter could be applied for improving the tolerance to abiotic stresses of the transgenic plants by promoting the expression of vegetable soybean AQPs.

An advanced systems biology framework of feature engineering for cold tolerance genes discovery from integrated omics and non-omics data in soybean

Soybean is sensitive to low temperatures during the crop growing season. An urgent demand for breeding cold-tolerant cultivars to alleviate the production loss is apparent to cope with this scenario. Cold-tolerant trait is a complex and quantitative trait controlled by multiple genes, environmental factors, and their interaction. In this study, we proposed an advanced systems biology framework of feature engineering for the discovery of cold tolerance genes (CTgenes) from integrated omics and non-omics (OnO) data in soybean. An integrative pipeline was introduced for feature selection and feature extraction from different layers in the integrated OnO data using data ensemble methods and the non-parameter random forest prioritization to minimize uncertainties and false positives for accuracy improvement of results. In total, 44, 143, and 45 CTgenes were identified in short-, mid-, and long-term cold treatment, respectively, from the corresponding gene-pool. These CTgenes outperformed the remaining genes, the random genes, and the other candidate genes identified by other approaches in an independent RNA-seq database. Furthermore, we applied pathway enrichment and crosstalk network analyses to uncover relevant physiological pathways with the discovery of underlying cold tolerance in hormone- and defense-related modules. Our CTgenes were validated by using 55 SNP genotype data of 56 soybean samples in cold tolerance experiments. This suggests that the CTgenes identified from our proposed systematic framework can effectively distinguish cold-resistant and cold-sensitive lines. It is an important advancement in the soybean cold-stress response. The proposed pipelines provide an alternative solution to biomarker discovery, module discovery, and sample classification underlying a particular trait in plants in a robust and efficient way.

Genome-Wide Association Study Reveals a Genetic Mechanism of Salt Tolerance Germinability in Rice (Oryza sativa L.)

Salt stress is one of the factors that limits rice production, and an important task for researchers is to cultivate rice with strong salt tolerance. In this study, 211 rice accessions were used to determine salt tolerance germinability (STG) indices and conduct a genome-wide association study (GWAS) using 36,727 SNPs. The relative germination energy (RGE), relative germination index (RGI), relative vigor index (RVI), relative mean germination time (RMGT), relative shoot length (RSL), and relative root length (RRL) were used to determine the STG indices in rice. A total of 43 QTLs, including 15 for the RGE, 6 for the RGI, 7 for the RVI, 3 for the RMGT, 1 for the RSL, and 11 for the RRL, were identified on nine chromosome regions under 60 and 100 mM NaCl conditions. For these STG-related QTLs, 18 QTLs were co-localized with previous studies, and some characterized salt-tolerance genes, such as OsCOIN, OsHsp17.0, and OsDREB2A, are located in these QTL candidates. Among the 25 novel QTLs, qRGE60-1-2 co-localized with qRGI60-1-1 on chromosome 1, and qRGE60-3-1 and qRVI60-3-1 co-localized on chromosome 3. According to the RNA-seq database, 16 genes, including nine for qRGE60-1-2 (qRGI60-1-1) and seven for qRGE60-3-1 (qRVI60-3-1), were found to show significant differences in their expression levels between the control and salt treatments. Furthermore, the expression patterns of these differentially expressed genes were analyzed, and nine genes (five for qRGE60-1-2 and four for qRGE60-3-1) were highly expressed in embryos at the germination stage. Haplotype analysis of these nine genes showed that the rice varieties with elite haplotypes in the LOC_Os03g13560, LOC_Os03g13840, and LOC_Os03g14180 genes had high STG. GWAS validated the known genes underlying salt tolerance and identified novel loci that could enrich the current gene pool related to salt tolerance. The resources with high STG and significant loci identified in this study are potentially useful in breeding for salt tolerance.

Evolutionary and Regulatory Pattern Analysis of Soybean Ca2+ ATPases for Abiotic Stress Tolerance

P₂-type Ca²⁺ ATPases are responsible for cellular Ca²⁺ transport, which plays an important role in plant development and tolerance to biotic and abiotic stresses. However, the role of P₂-type Ca²⁺ ATPases in stress response and stomatal regulation is still elusive in soybean. In this study, a total of 12 P₂-type Ca²⁺ ATPases genes (GmACAs and GmECAs) were identified from the genome of Glycine max. We analyzed the evolutionary relationship, conserved motif, functional domain, gene structure and location, and promoter elements of the family. Chlorophyll fluorescence imaging analysis showed that vegetable soybean leaves are damaged to different extents under salt, drought, cold, and shade stresses. Real-time quantitative PCR (RT-qPCR) analysis demonstrated that most of the GmACAs and GmECAs are up-regulated after drought, cold, and NaCl treatment, but are down-regulated after shading stress. Microscopic observation showed that different stresses caused significant stomatal closure. Spatial location and temporal expression analysis suggested that GmACA8, GmACA9, GmACA10, GmACA12, GmACA13, and GmACA11 might promote stomatal closure under drought, cold, and salt stress. GmECA1 might regulate stomatal closure in shading stress. GmACA1 and GmECA3 might have a negative function on cold stress. The results laid an important foundation for further study on the function of P₂-type Ca²⁺ ATPase genes GmACAs and GmECAs for breeding abiotic stress-tolerant vegetable soybean.

A Multi-Level Iterative Bi-Clustering Method for Discovering miRNA Co-regulation Network of Abiotic Stress Tolerance in Soybeans

Although growing evidence shows that microRNA (miRNA) regulates plant growth and development, miRNA regulatory networks in plants are not well understood. Current experimental studies cannot characterize miRNA regulatory networks on a large scale. This information gap provides an excellent opportunity to employ computational methods for global analysis and generate valuable models and hypotheses. To address this opportunity, we collected miRNA-target interactions (MTIs) and used MTIs from Arabidopsis thaliana and Medicago truncatula to predict homologous MTIs in soybeans, resulting in 80,235 soybean MTIs in total. A multi-level iterative bi-clustering method was developed to identify 483 soybean miRNA-target regulatory modules (MTRMs). Furthermore, we collected soybean miRNA expression data and corresponding gene expression data in response to abiotic stresses. By clustering these data, 37 MTRMs related to abiotic stresses were identified, including stress-specific MTRMs and shared MTRMs. These MTRMs have gene ontology (GO) enrichment in resistance response, iron transport, positive growth regulation, etc. Our study predicts soybean MTRMs and miRNA-GO networks under different stresses, and provides miRNA targeting hypotheses for experimental analyses. The method can be applied to other biological processes and other plants to elucidate miRNA co-regulation mechanisms.

Using the Knowledge of Post-transcriptional Regulations to Guide Gene Selections for Molecular Breeding in Soybean

The omics approaches allow the scientific community to successfully identify genomic regions associated with traits of interest for marker-assisted breeding. Agronomic traits such as seed color, yield, growth habit, and stress tolerance have been the targets for soybean molecular breeding. Genes governing these traits often undergo post-transcriptional modifications, which should be taken into consideration when choosing elite genes for molecular breeding. Post-transcriptional regulations of genes include transcript regulations, protein modifications, and even the regulation of the translational machinery. Transcript regulations involve elements such as microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) for the maintenance of transcript stability or regulation of translation efficiency. Protein modifications involve molecular modifications of target proteins and the alterations of their interacting partners. Regulations of the translational machinery include those on translation factors and the ribosomal protein complex. Post-transcriptional regulations usually involve a set of genes instead of a single gene. Such a property may facilitate molecular breeding. In this review, we will discuss the post-transcriptional modifications of genes related to favorable agronomic traits such as stress tolerance, growth, and nutrient uptake, using examples from soybean as well as other crops. The examples from other crops may guide the selection of genes for marker-assisted breeding in soybean.

Transcriptomic analysis of Vigna radiata in response to chilling stress and uniconazole application

Background

Chilling injury of mung bean (Vigna radiata (L.)) during the blooming and podding stages is a major agricultural threat in Northeast China. Uniconazole (UNZ) can alleviate water deficit stress in soybean and waterlogging stress in mung bean. However, there has been no report on the effect of UNZ application on the growth and transcriptomic profile of mung bean under chilling stress.

Results

UNZ application before chilling stress at the R1 stage alleviated the decline in mung bean yield. UNZ delayed the decrease in leaf chlorophyll content under chilling stress at the R1 stage and accelerated the increase in leaf chlorophyll content during the recovery period. Eighteen separate RNA-Seq libraries were generated from RNA samples collected from leaves exposed to six different treatment schemes. The numbers of DEGs specific for UNZ treatment between D1 + S vs. D1 and D4 + S vs. D4 were 708 and 810, respectively. GO annotations showed that photosynthesis genes were obviously enriched among the genes affected by chilling stress and UNZ application. KEGG pathway enrichment analysis indicated that 4 pathways (cutin, suberin and wax biosynthesis; photosynthesis; porphyrin and chlorophyll metabolism; and ribosome) were downregulated, while plant–pathogen interaction was upregulated, by chilling stress. UNZ application effectively prevented the further downregulation of the gene expression of members of these 4 KEGG pathways under chilling stress.

Conclusions

UNZ application effectively delayed the decrease in photosynthetic pigment content under chilling stress and accelerated the increase in photosynthetic pigment content during the recovery period, thus effectively limiting the decline in mung bean yield. UNZ application effectively prevented the further downregulation of the gene expression of members of 4 KEGG pathways under chilling stress and increased mung bean tolerance to chilling stress.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08443-6.

Genome sequencing and population resequencing provide insights into the genetic basis of domestication and diversity of vegetable soybean

Vegetable soybean is one of the most important vegetables in China, and the demand for this vegetable has markedly increased worldwide over the past two decades. Here, we present a high-quality de novo genome assembly of the vegetable soybean cultivar Zhenong 6 (ZN6), which is one of the most popular cultivars in China. The 20 pseudochromosomes cover 94.57% of the total 1.01 Gb assembly size, with contig N50 of 3.84 Mb and scaffold N50 of 48.41 Mb. A total of 55 517 protein-coding genes were annotated. Approximately 54.85% of the assembled genome was annotated as repetitive sequences, with the most abundant long terminal repeat transposable elements. Comparative genomic and phylogenetic analyses with grain soybean Williams 82, six other Fabaceae species and Arabidopsis thaliana genomes highlight the difference of ZN6 with other species. Furthermore, we resequenced 60 vegetable soybean accessions. Alongside 103 previously resequenced wild soybean and 155 previously resequenced grain soybean accessions, we performed analyses of population structure and selective sweep of vegetable, grain, and wild soybean. They were clearly divided into three clades. We found 1112 and 1047 genes under selection in the vegetable soybean and grain soybean populations compared with the wild soybean population, respectively. Among them, we identified 134 selected genes shared between vegetable soybean and grain soybean populations. Additionally, we report four sucrose synthase genes, one sucrose-phosphate synthase gene, and four sugar transport genes as candidate genes related to important traits such as seed sweetness and seed size in vegetable soybean. This study provides essential genomic resources to promote evolutionary and functional genomics studies and genomically informed breeding for vegetable soybean.

Identification of the Complex Interplay Between Nematode-Related lncRNAs and Their Target Genes in Glycine max L

Soybean (Glycine max) is a major plant protein source and oilseed crop. However, plant-parasitic nematodes (PPNs) affect its annual yield. In the current study, in order to better understand the regulation of defense mechanism against PPNs in soybean, we investigated the role of long non-coding RNAs (lncRNAs) in response to two nematode species, Heterodera glycines (SCN: soybean cyst nematode) and Rotylenchulus reniformis (reniform). To this end, two publicly available RNA-seq data sets (SCN data set and RAD: reniform-associated data set) were employed to discover the lncRNAome profile of soybean under SCN and reniform infection, respectively. Upon identification of unannotated transcripts in these data sets, a seven-step pipeline was utilized to sieve these transcripts, which ended up in 384 and 283 potential lncRNAs in SCN data set and RAD, respectively. These transcripts were then used to predict cis and trans nematode-related targets in soybean genome. Computational prediction of target genes function, some of which were also among differentially expressed genes, revealed the involvement of putative nematode-responsive genes as well as enrichment of multiple stress responses in both data sets. Finally, 15 and six lncRNAs were proposed to be involved in microRNA-mediated regulation of gene expression in soybean in response to SNC and reniform infection, respectively. Collectively, this study provides a novel insight into the signaling and regulatory network of soybean-pathogen interactions and opens a new window for further research.

Enhancement of plant cold tolerance by soybean RCC1 family gene GmTCF1a

BACKGROUND

Low temperature severely limits the growth, yield, and geographic distributions of soybean. Soybean plants respond to cold stress by reprogramming the expression of a series of cold-responsive genes. However, the intrinsic mechanism underlying cold-stress tolerance in soybean remains unclear. A. thaliana tolerant to chilling and freezing 1 (AtTCF1) is a regulator of chromosome condensation 1 (RCC1) family protein and regulates freezing tolerance through an independent C-repeat binding transcription factor (CBF) signaling pathway.

RESULTS

In this study, we identified a homologous gene of AtTCF1 in soybean (named GmTCF1a), which mediates plant tolerance to low temperature. Like AtTCF1, GmTCF1a contains five RCC1 domains and is located in the nucleus. GmTCF1a is strongly and specifically induced by cold stress. Interestingly, ectopic overexpression of GmTCF1a in Arabidopsis greatly increased plant survival rate and decreased electrolyte leakage under freezing stress. A cold-responsive gene, COR15a, was highly induced in the GmTCF1a-overexpressing transgenic lines.

CONCLUSIONS

GmTCF1a responded specifically to cold stress, and ectopic expression of GmTCF1a enhanced cold tolerance and upregulated COR15a levels. These results indicate that GmTCF1a positively regulates cold tolerance in soybean and may provide novel insights into genetic improvement of cold tolerance in crops.

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

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

Title: Hypermethylation of miRNA Genes During Nodule Development

DNA methylation has recently emerged as a powerful regulatory mechanism controlling the expression of key regulators of various developmental processes, including nodulation. However, the functional role of DNA methylation in regulating the expression of microRNA (miRNA) genes during the formation and development of nitrogen-fixing nodules remains largely unknown. In this study, we profiled DNA methylation patterns of miRNA genes during nodule formation, development, and early senescence stages in soybean (Glycine max) through the analysis of methylC-seq data. Absolute DNA methylation levels in the CG, CHH, and CHH sequence contexts over the promoter and primary transcript regions of miRNA genes were significantly higher in the nodules compared with the corresponding root tissues at these three distinct nodule developmental stages. We identified a total of 82 differentially methylated miRNAs in the nodules compared with roots. Differential DNA methylation of these 82 miRNAs was detected only in the promoter (69), primary transcript region (3), and both in the promoter and primary transcript regions (10). The large majority of these differentially methylated miRNAs were hypermethylated in nodules compared with the corresponding root tissues and were found mainly in the CHH context and showed stage-specific methylation patterns. Differentially methylated regions in the promoters of 25 miRNAs overlapped with transposable elements, a finding that may explain the vulnerability of miRNAs to DNA methylation changes during nodule development. Gene expression analysis of a set of promoter-differentially methylated miRNAs pointed to a negative association between DNA methylation and miRNA expression. Gene Ontology and pathways analyses indicate that changes in DNA methylation of miRNA genes are reprogrammed and contribute to nodule development through indirect regulation of genes involved in cellular processes and pathways with well-established roles in nodulation.

MicroRNA-Mediated Responses to Chromium Stress Provide Insight Into Tolerance Characteristics of Miscanthus sinensis

Chromium (Cr) is a heavy metal in nature, which poses a potential risk to toxicity to both animals and plants when releasing into the environment. However, the regulation of microRNA (miRNA)-mediated response to heavy metal Cr has not been studied in Miscanthus sinensis. In this study, based on high-throughput miRNA sequencing, a total of 104 conserved miRNAs and 158 nonconserved miRNAs were identified. Among them, there were 45 differentially expressed miRNAs in roots and 13 differentially expressed miRNAs in leaves. The hierarchical clustering analysis showed that these miRNAs were preferentially expressed in a certain tissue. There were 833 differentially expressed target genes of 45 miRNAs in roots and 280 differentially expressed target genes of 13 miRNA in leaves. After expression trend analysis, five significantly enriched modules were obtained in roots, and three significantly enriched trend blocks in leaves. Based on the candidate gene annotation and gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) function analysis, miR167a, novel_miR15, and novel_miR22 and their targets were potentially involved in Cr transportation and chelation. Besides, miR156a, miR164, miR396d, and novel_miR155 were identified as participating in the physiological and biochemical metabolisms and the detoxification of Cr of plants. The results demonstrated the critical role of miRNA-mediated responses to Cr treatment in M. sinensis, which involves ion uptake, transport, accumulation, and tolerance characteristics.

Identification of Differentially Methylated miRNA Genes During Compatible and Incompatible Interactions Between Soybean and Soybean Cyst Nematode

DNA methylation is a widespread epigenetic mark that affects gene expression and transposon mobility during plant development and stress responses. However, the role of DNA methylation in regulating the expression of microRNA (miRNA) genes remains largely unexplored. Here, we analyzed DNA methylation changes of miRNA genes using a pair of soybean (Glycine max) near-isogenic lines (NILs) differing in their response to soybean cyst nematode (SCN; Heterodera glycines). Differences in global DNA methylation levels over miRNA genes in response to SCN infection were observed between the isogenic lines. miRNA genes with significant changes in DNA methylation levels in the promoter and primary transcript-coding regions were detected in both lines. In the susceptible isogenic line (NIL-S), 82 differentially methylated miRNAs were identified in response to SCN infection whereas, in the resistant isogenic line (NIL-R), only 16 differentially methylated miRNAs were identified. Interestingly, gma-miR5032, gma-miR5043, gma-miR1520b, and gma-2107-ch16 showed opposite methylation patterns in the isogenic lines. In addition, the miRNA paralogs gma-miR5770a and gma-miR5770b showed hypermethylation and hypomethylation in NIL-S and NIL-R, respectively. Gene expression quantification of gma-miR5032, gma-miR5043, gma-miR1520b, and gma-miR5770a/b and their confirmed targets indicated a role of DNA methylation in regulating miRNA expression and, thus, their targets upon SCN infection. Furthermore, overexpression of these four miRNAs in NIL-S using transgenic hairy root system enhanced plant resistance to SCN to various degrees with a key role observed for miR5032. Together, our results provide new insights into the role of epigenetic mechanisms in controlling miRNA regulatory function during SCN-soybean interactions.[Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Stress induced dynamic adjustment of conserved miR164:NAC module

Aims including the rationale

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

Methods

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

Key results

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

Main conclusion

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

Integration of sRNA, degradome, transcriptome analysis and functional investigation reveals gma-miR398c negatively regulates drought tolerance via GmCSDs and GmCCS in transgenic Arabidopsis and soybean

BACKGROUND

Drought conditions adversely affect soybean growth, resulting in severe yield losses worldwide. Increasing experimental evidence indicates miRNAs are important post-transcriptional regulators of gene expression. However, the drought-responsive molecular mechanism underlying miRNA-mRNA interactions remains largely uncharacterized in soybean. Meanwhile, the miRNA-regulated drought response pathways based on multi-omics approaches remain elusive.

RESULTS

We combined sRNA, transcriptome and degradome sequencing to elucidate the complex regulatory mechanism mediating soybean drought resistance. One-thousand transcripts from 384 target genes of 365 miRNAs, which were enriched in the peroxisome, were validated by degradome-seq. An integrated analysis showed 42 miRNA-target pairs exhibited inversely related expression profiles. Among these pairs, a strong induction of gma-miR398c as a major gene negatively regulates multiple peroxisome-related genes (GmCSD1a/b, GmCSD2a/b/c and GmCCS). Meanwhile, we detected that alternative splicing of GmCSD1a/b might affect soybean drought tolerance by bypassing gma-miR398c regulation. Overexpressing gma-miR398c in Arabidopsis thaliana L. resulted in decreased percentage germination, increased leaf water loss, and reduced survival under water deficiency, which displayed sensitivity to drought during seed germination and seedling growth. Furthermore, overexpressing gma-miR398c in soybean decreased GmCSD1a/b, GmCSD2a/b/c and GmCCS expression, which weakened the ability to scavenge O2.-, resulting in increased relative electrolyte leakage and stomatal opening compared with knockout miR398c and wild-type soybean under drought conditions.

CONCLUSION

The study indicates that gma-miR398c negatively regulates soybean drought tolerance, and provides novel insights useful for breeding programs to improve drought resistance by CRISPR technology.

Elucidating the regulatory roles of microRNAs in maize (Zea mays L.) leaf growth response to chilling stress

MAIN CONCLUSION: miRNAs control leaf size of maize crop during chilling stress tolerance by regulating developmentally important transcriptional factors and sustaining redox homeostasis of cells. Chilling temperature (0-15 °C) is a major constraint for the cultivation of maize (Zea mays) which inhibits the early growth of maize leading to reduction in leaf size. Growth and development take place in meristem, elongation, and mature zones that are linearly located along the leaf base to tip. To prevent shortening of leaf caused by chilling, this study aims to elucidate the regulatory roles of microRNA (miRNA) genes in the controlling process switching between growth and developmental stages. In this respect, hybrid maize ADA313 seedlings were treated to the chilling temperature which caused 26% and 29% reduction in the final leaf length and a decline in cell production of the fourth leaf. The flow cytometry data integrated with the expression analysis of cell cycle genes indicated that the reason for the decline was a failure proceeding from G2/M rather than G1/S. Through an miRNome analysis of 321 known maize miRNAs, 24, 6, and 20 miRNAs were assigned to putative meristem, elongation, and mature zones, respectively according to their chilling response. To gain deeper insight into decreased cell production, in silico, target prediction analysis was performed for meristem specific miRNAs. Among the miRNAs, miR160, miR319, miR395, miR396, miR408, miR528, and miR1432 were selected for confirming the potential of negative regulation with their predicted targets by qRT-PCR. These findings indicated evidence for improvement of growth and yield under chilling stress of the maize.

Construction and analysis of degradome-dependent microRNA regulatory networks in soybean

BACKGROUND

Usually the microRNA (miRNA)-mediated gene regulatory network (GRN) is constructed from the investigation of miRNA expression profiling and target predictions. However, the higher/lower expression level of miRNAs does not always indicate the higher/lower level of cleavages and such analysis, thus, sometimes ignores the crucial cleavage events. In the present work, the degradome sequencing data were employed to construct the complete miRNA-mediated gene regulatory network in soybean, unlike the traditional approach starting with small RNA sequencing data.

RESULTS

We constructed the root-, cotyledon-, leaf- and seed-specific miRNA regulatory networks with the degradome sequencing data and the forthcoming verification of miRNA profiling analysis. As a result, we identified 205 conserved miRNA-target interactions (MTIs) involved with 6 conserved gma-miRNA families and 365 tissue-specific MTIs containing 24 root-specific, 45 leaf-specific, 63 cotyledon-specific and 225 seed-specific MTIs. We found a total of 156 miRNAs in tissue-specific MTIs including 18 tissue-specific miRNAs, however, only 3 miRNAs have consistent tissue-specific expression. Our study showed the degradome-dependent miRNA regulatory networks (DDNs) in four soybean tissues and explored their conservations and specificities.

CONCLUSIONS

The construction of DDNs may provide the complete miRNA-Target interactions in certain plant tissues, leading to the identification of the conserved and tissue-specific MTIs and sub-networks. Our work provides a basis for further investigation of the roles and mechanisms of miRNA-mediated regulation of tissue-specific growth and development in soybean.

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.

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.

Regulation of low temperature stress in plants by microRNAs

Low temperature is one of the most common environmental stresses that seriously affect the growth and development of plants. However, plants have the plasticity in their defence mechanisms enabling them to tolerate and, sometimes, even survive adverse environmental conditions. MicroRNAs (miRNAs) are small non-coding RNAs, approximately 18-24 nucleotides in length, and are being increasingly recognized as regulators of gene expression at the post-transcriptional level and have the ability to influence a broad range of biological processes. There is growing evidence in the literature that reprogramming of gene expression mediated through miRNAs is a major defence mechanism in plants enabling them to respond to stresses. To date, numerous studies have established the importance of miRNA-based regulation of gene expression under low temperature stress. Individual miRNAs can modulate the expression of multiple mRNA targets, and, therefore, the manipulation of a single miRNA has the potential to affect multiple biological processes. Numerous functional studies have attempted to identify the miRNA-target interactions and have elaborated the role of several miRNAs in cold-stress regulation. This review summarizes the current understanding of miRNA-mediated modulation of the expression of key genes as well as genetic and regulatory pathways, involved in low temperature stress responses in plants.

High throughput deep sequencing reveals the important roles of microRNAs during sweetpotato storage at chilling temperature

Sweetpotato (Impomoea batatas L.) is a globally important economic food crop with a potential of becoming a bioenergy and pharmaceutical crop. Thus, studying the molecular mechanism of tuberous root development and storage is very important. However, not too much progress has been made in this field. In this study, we employed the next generation high-throughput deep sequencing technology to sequence all small RNAs and degradome of sweetpotato for systematically investigating sweetpotato response to chilling stress during storage. A total of 190 known microRNAs (miRNAs) and 191 novel miRNAs were identified, and 428 transcripts were targeted by 184 identified miRNAs. More importantly, we identified 26 miRNAs differentially expressed between chilling stress and control conditions. The expression of these miRNAs and their targets was also confirmed by qRT-PCR. Integrated analysis of small RNAs and degradome sequencing reveals that miRNA-mediated SA signaling, ABA-dependent, and ROS response pathways are involved in sweetpotato root response to chilling stress during storage.

Genome-wide identification of soybean microRNA responsive to soybean cyst nematodes infection by deep sequencing

BACKGROUND

The soybean cyst nematode (SCN), Heterodera glycines, is one of the most devastating diseases limiting soybean production worldwide. It is known that small RNAs, including microRNAs (miRNAs) and small interfering RNAs (siRNAs), play important roles in regulating plant growth and development, defense against pathogens, and responses to environmental changes.

RESULTS

In order to understand the role of soybean miRNAs during SCN infection, we analyzed 24 small RNA libraries including three biological replicates from two soybean cultivars (SCN susceptible KS4607, and SCN HG Type 7 resistant KS4313N) that were grown under SCN-infested and -noninfested soil at two different time points (SCN feeding establishment and egg production). In total, 537 known and 70 putative novel miRNAs in soybean were identified from a total of 0.3 billion reads (average about 13.5 million reads for each sample) with the programs of Bowtie and miRDeep2 mapper. Differential expression analyses were carried out using edgeR to identify miRNAs involved in the soybean-SCN interaction. Comparative analysis of miRNA profiling indicated a total of 60 miRNAs belonging to 25 families that might be specifically related to cultivar responses to SCN. Quantitative RT-PCR validated similar miRNA interaction patterns as sequencing results.

CONCLUSION

These findings suggest that miRNAs are likely to play key roles in soybean response to SCN. The present work could provide a framework for miRNA functional identification and the development of novel approaches for improving soybean SCN resistance in future studies.

Conservation and diversification of the miR166 family in soybean and potential roles of newly identified miR166s

BACKGROUND

microRNA166 (miR166) is a highly conserved family of miRNAs implicated in a wide range of cellular and physiological processes in plants. miR166 family generally comprises multiple miR166 members in plants, which might exhibit functional redundancy and specificity. The soybean miR166 family consists of 21 members according to the miRBase database. However, the evolutionary conservation and functional diversification of miR166 family members in soybean remain poorly understood.

RESULTS

We identified five novel miR166s in soybean by data mining approach, thus enlarging the size of miR166 family from 21 to 26 members. Phylogenetic analyses of the 26 miR166s and their precursors indicated that soybean miR166 family exhibited both evolutionary conservation and diversification, and ten pairs of miR166 precursors with high sequence identity were individually grouped into a discrete clade in the phylogenetic tree. The analysis of genomic organization and evolution of MIR166 gene family revealed that eight segmental duplications and four tandem duplications might occur during evolution of the miR166 family in soybean. The cis-elements in promoters of MIR166 family genes and their putative targets pointed to their possible contributions to the functional conservation and diversification. The targets of soybean miR166s were predicted, and the cleavage of ATHB14-LIKE transcript was experimentally validated by RACE PCR. Further, the expression patterns of the five newly identified MIR166s and 12 target genes were examined during seed development and in response to abiotic stresses, which provided important clues for dissecting their functions and isoform specificity.

CONCLUSION

This study enlarged the size of soybean miR166 family from 21 to 26 members, and the 26 soybean miR166s exhibited evolutionary conservation and diversification. These findings have laid a foundation for elucidating functional conservation and diversification of miR166 family members, especially during seed development or under abiotic stresses.

Small RNA-mediated responses to low- and high-temperature stresses in cotton

MicroRNAs (miRNAs) are one class of endogenous non-coding RNAs modulating the expression of target genes involved in plant development and stress tolerance, by degrading mRNA or repressing translation. In this study, small RNA and mRNA degradome sequencing were used to identify low- and high-temperature stress-responsive miRNAs and their targets in cotton (Gossypium hirsutum). Cotton seedlings were treated under different temperature conditions (4, 12, 25, 35, and 42 °C) and then the effects were investigated. In total, 319 known miRNAs and 800 novel miRNAs were identified, and 168 miRNAs were differentially expressed between different treatments. The targets of these miRNAs were further analysed by degradome sequencing. Based on studies from Gene Ontology and Kyoto Encyclopedia of Genes and Genomes, the majority of the miRNAs are from genes that are likely involved in response to hormone stimulus, oxidation-reduction reaction, photosynthesis, plant-pathogen interaction and plant hormone signal transduction pathways. This study provides new insight into the molecular mechanisms of plant response to extreme temperature stresses, and especially the roles of miRNAs under extreme temperatures.