Identification of Cold-Responsive miRNAs and Their Target Genes in Nitrogen-Fixing Nodules of Soybean

Identification and Expression Analysis of the Nucleotidyl Transferase Protein (NTP) Family in Soybean (Glycine max) under Various Abiotic Stresses

Nucleotidyl transferases (NTPs) are common transferases in eukaryotes and play a crucial role in nucleotide modifications at the 3' end of RNA. In plants, NTPs can regulate RNA stability by influencing 3' end modifications, which in turn affect plant growth, development, stress responses, and disease resistance. Although the functions of NTP family members have been extensively studied in Arabidopsis, rice, and maize, there is limited knowledge about NTP genes in soybeans. In this study, we identified 16 members of the NTP family in soybeans, including two subfamilies (G1 and G2) with distinct secondary structures, conserved motifs, and domain distributions at the protein level. Evolutionary analysis of genes in the NTP family across multiple species and gene collinearity analysis revealed a relatively conserved evolutionary pattern. Analysis of the tertiary structure of the proteins showed that NTPs have three conserved aspartic acids that bind together to form a possible active site. Tissue-specific expression analysis indicated that some NTP genes exhibit tissue-specific expression, likely due to their specific functions. Stress expression analysis showed significant differences in the expression levels of NTP genes under high salt, drought, and cold stress. Additionally, RNA-seq analysis of soybean plants subjected to salt and drought stress further confirmed the association of soybean NTP genes with abiotic stress responses. Subcellular localization experiments revealed that GmNTP2 and GmNTP14, which likely have similar functions to HESO1 and URT1, are located in the nucleus. These research findings provide a foundation for further investigations into the functions of NTP family genes in soybeans.

miRNAs for crop improvement

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

Genome-Wide Identification and Characterization of the Soybean Snf2 Gene Family and Expression Response to Rhizobia

Sucrose nonfermenting 2 (Snf2) family proteins are the core component of chromatin remodeling complexes that can alter chromatin structure and nucleosome position by utilizing the energy of ATP, playing a vital role in transcription regulation, DNA replication, and DNA damage repair. Snf2 family proteins have been characterized in various species including plants, and they have been found to regulate development and stress responses in Arabidopsis. Soybean (Glycine max) is an important food and economic crop worldwide, unlike other non-leguminous crops, soybeans can form a symbiotic relationship with rhizobia for biological nitrogen fixation. However, little is known about Snf2 family proteins in soybean. In this study, we identified 66 Snf2 family genes in soybean that could be classified into six groups like Arabidopsis, unevenly distributed on 20 soybean chromosomes. Phylogenetic analysis with Arabidopsis revealed that these 66 Snf2 family genes could be divided into 18 subfamilies. Collinear analysis showed that segmental duplication was the main mechanism for expansion of Snf2 genes rather than tandem repeats. Further evolutionary analysis indicated that the duplicated gene pairs had undergone purifying selection. All Snf2 proteins contained seven domains, and each Snf2 protein had at least one SNF2_N domain and one Helicase_C domain. Promoter analysis revealed that most Snf2 genes had cis-elements associated with jasmonic acid, abscisic acid, and nodule specificity in their promoter regions. Microarray data and real-time quantitative PCR (qPCR) analysis revealed that the expression profiles of most Snf2 family genes were detected in both root and nodule tissues, and some of them were found to be significantly downregulated after rhizobial infection. In this study, we conducted a comprehensive analysis of the soybean Snf2 family genes and demonstrated their responsiveness to Rhizobia infection. This provides insight into the potential roles of Snf2 family genes in soybean symbiotic nodulation.

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.

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.

Innovation and appropriation in mycorrhizal and rhizobial Symbioses

Most land plants benefit from endosymbiotic interactions with mycorrhizal fungi, including legumes and some nonlegumes that also interact with endosymbiotic nitrogen (N)-fixing bacteria to form nodules. In addition to these helpful interactions, plants are continuously exposed to would-be pathogenic microbes: discriminating between friends and foes is a major determinant of plant survival. Recent breakthroughs have revealed how some key signals from pathogens and symbionts are distinguished. Once this checkpoint has been passed and a compatible symbiont is recognized, the plant coordinates the sequential development of two types of specialized structures in the host. The first serves to mediate infection, and the second, which appears later, serves as sophisticated intracellular nutrient exchange interfaces. The overlap in both the signaling pathways and downstream infection components of these symbioses reflects their evolutionary relatedness and the common requirements of these two interactions. However, the different outputs of the symbioses, phosphate uptake versus N fixation, require fundamentally different components and physical environments and necessitated the recruitment of different master regulators, NODULE INCEPTION-LIKE PROTEINS, and PHOSPHATE STARVATION RESPONSES, for nodulation and mycorrhization, respectively.

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.

GmSPX8, a nodule-localized regulator confers nodule development and nitrogen fixation under phosphorus starvation in soybean

BACKGROUND

Biological nitrogen fixation (BNF) is an important nitrogen source for legume plants, and highly efficient nitrogen fixation requires sufficient phosphorus (P). However, the mechanism of maintaining nitrogen fixation of the legume nodules under low P concentration remains largely unknown.

RESULTS

A nodule-localized SPX protein, GmSPX8, was discovered by transcriptome and functional analysis of its role in N2 fixation was characterized in soybean nodules. GmSPX8 was preferentially expressed in nodules and its expression was gradually increased during nodule development. And also the expression pattern was investigated using reporter gene β-glucuronidase (GUS) driven by the promoter of GmSPX8. GmSPX8 was greatly induced and the GUS activity was increased by 12.2% under P deficiency. Overexpression of GmSPX8 in transgenic plants resulted in increased nodule number, nodule fresh weight and nitrogenase activity by 15.0%, 16.0%, 42.5%, subsequently leading to increased N and P content by 17.0% and 19.0%, while suppression of GmSPX8 showed significantly impaired nodule development and nitrogen fixation efficiency under low P stress. These data indicated that GmSPX8 conferred nodule development and nitrogen fixation under low P condition. By yeast two-hybrid screening, GmPTF1 was identified as a potential interacting protein of GmSPX8, which was further confirmed by BiFC, Y2H and pull down assay. Transcript accumulation of GmPTF1 and its downstream genes such as GmEXLB1 and EXPB2 were increased in GmSPX8 overexpressed transgenic nodules, and in the presence of GmSPX8, the transcriptional activity of GmPTF1 in yeast cells and tobacco leaves was greatly enhanced.

CONCLUSIONS

In summary, these findings contribute novel insights towards the role of GmSPX8 in nodule development and nitrogen fixation partly through interacting with GmPTF1 in soybean under low P condition.

Differential microRNA expression, microRNA arm switching, and microRNA:long noncoding RNA interaction in response to salinity stress in soybean

BACKGROUND

Soybean is a major legume crop with high nutritional and environmental values suitable for sustainable agriculture. Noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), are important regulators of gene functions in eukaryotes. However, the interactions between these two types of ncRNAs in the context of plant physiology, especially in response to salinity stress, are poorly understood.

RESULTS

Here, we challenged a cultivated soybean accession (C08) and a wild one (W05) with salt treatment and obtained their small RNA transcriptomes at six time points from both root and leaf tissues. In addition to thoroughly analyzing the differentially expressed miRNAs, we also documented the first case of miRNA arm-switching (miR166m), the swapping of dominant miRNA arm expression, in soybean in different tissues. Two arms of miR166m target different genes related to salinity stress (chloroplastic beta-amylase 1 targeted by miR166m-5p and calcium-dependent protein kinase 1 targeted by miR166m-3p), suggesting arm-switching of miR166m play roles in soybean in response to salinity stress. Furthermore, two pairs of miRNA:lncRNA interacting partners (miR166i-5p and lncRNA Gmax_MSTRG.35921.1; and miR394a-3p and lncRNA Gmax_MSTRG.18616.1) were also discovered in reaction to salinity stress.

CONCLUSIONS

This study demonstrates how ncRNA involves in salinity stress responses in soybean by miRNA arm switching and miRNA:lncRNA interactions. The behaviors of ncRNAs revealed in this study will shed new light on molecular regulatory mechanisms of stress responses in plants, and hence provide potential new strategies for crop improvement.

Comprehensive Identification and Profiling of miRNAs Involved in Terpenoid Synthesis of Gleditsia sinensis Lam

Gleditsia sinensis Lam. is a tree with worldwide distribution and important economic and medicinal values; its pods contain terpenoids including gleditsioside, thiamine, and brassinosteroids. However, thus far, there are few studies on the terpenoid regulation of G. sinensis at the molecular level. microRNA (miRNA) is a class of small RNAs with conserved and crucial roles in the regulation of diverse biological processes during plant growth and development. To identify the miRNAs of G. sinensis and evaluate their involvement in terpenoid synthesis, this investigation quantified the content changes in saponins in pods at three developmental stages: May (pod-setting stage), July (elongation stage), and September (browning stage), and then we performed genome-wide miRNA profiles during the three development stages of the G. sinensis pods. A total of 351 conserved miRNAs belonging to 216 families were identified, among which 36 conserved miRNAs exist specifically in legumes. Through target analysis, 708 unigenes were predicted to be candidate targets of 37 differentially expressed miRNAs. The targets of miR838-3p and miR2093-5p were involved in the derived branches of monoterpenes and gleditsioside, in brassinosteroid biosynthesis (BRB), and in indole alkaloid biosynthesis (IAB). Intriguingly, the targets of miR829-3p.1 were predicted to take part in thiamine biosynthesis, and the targets of miR4414b and miR5037a were involved in the main process of cytokinin synthesis. The corresponding targets participated in BRB, IAB, and terpenoid backbone biosynthesis, which were enriched significantly, suggesting that miR2093-5p, miR4414b, miR5037a, miR829-3p.1, and miR838-3p play indispensable roles in the regulation of triterpenoid saponin and monoterpenoid biosynthesis. To date, this is the first report of miRNA identification in G. sinensis and miRNA expression profiles at different developmental stages of G. sinensis pods, which provides a basis for further uncovering the molecular regulation of terpenoid synthesis in G. sinensis and new insights into the role of miRNAs in legumes.

Genome-Wide Characterization of Salt-Responsive miRNAs, circRNAs and Associated ceRNA Networks in Tomatoes

Soil salinization is a major environmental stress that causes crop yield reductions worldwide. Therefore, the cultivation of salt-tolerant crops is an effective way to sustain crop yield. Tomatoes are one of the vegetable crops that are moderately sensitive to salt stress. Global market demand for tomatoes is huge and growing. In recent years, the mechanisms of salt tolerance in tomatoes have been extensively investigated; however, the molecular mechanism through which non-coding RNAs (ncRNAs) respond to salt stress is not well understood. In this study, we utilized small RNA sequencing and whole transcriptome sequencing technology to identify salt-responsive microRNAs (miRNAs), messenger RNAs (mRNAs), and circular RNAs (circRNAs) in roots of M82 cultivated tomato and Solanum pennellii (S. pennellii) wild tomato under salt stress. Based on the theory of competitive endogenous RNA (ceRNA), we also established several salt-responsive ceRNA networks. The results showed that circRNAs could act as miRNA sponges in the regulation of target mRNAs of miRNAs, thus participating in the response to salt stress. This study provides insights into the mechanisms of salt tolerance in tomatoes and serves as an effective reference for improving the salt tolerance of salt-sensitive cultivars.

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.

Shoot-derived miR2111 controls legume root and nodule development

Legumes control their nodule numbers through the autoregulation of nodulation (AON). Rhizobia infection stimulates the production of root-derived CLE peptide hormones that are translocated to the shoot where they regulate a new signal. We used soybean to demonstrate that this shoot-derived signal is miR2111, which is transported via phloem to the root where it targets transcripts of Too Much Love (TML), a negative regulator of nodulation. Shoot perception of rhizobia-induced CLE peptides suppresses miR2111 expression, resulting in TML accumulation in roots and subsequent inhibition of nodule organogenesis. Feeding synthetic mature miR2111 via the petiole increased nodule numbers per plant. Likewise, elevating miR2111 availability by over-expression promoted nodulation, while target mimicry of TML induced the opposite effect on nodule development in wild-type plants and alleviated the supernodulating and stunted root growth phenotypes of AON-defective mutants. Additionally, in non-nodulating wild-type plants, ectopic expression of miR2111 significantly enhanced lateral root emergence with a decrease in lateral root length and average root diameter. In contrast, hairy roots constitutively expressing the target mimic construct exhibited reduced lateral root density. Overall, these findings demonstrate that miR2111 is both the critical shoot-to-root factor that positively regulates root nodule development and also acts to shape root system architecture.

Dynamics of miRNA mediated regulation of legume symbiosis

Symbiotic nitrogen fixation in legume nodules is important in soils with low nitrogen availability. The initiation and sustainability of symbiosis require cellular reprogramming that involves the miRNA-mediated inhibition or activation of specific nodulation genes. The high-throughput sequencing of small RNA libraries has identified miRNAs and their targets, which are the major players in the post-transcriptional gene regulation (PTGS) of the different stages of legume-rhizobia symbiosis ranging from bacterial colonization and organogenesis to symbiotic nitrogen fixation. Here, we present an overview of information obtained from the miRNA libraries from nodulating tissues that have been sequenced to date. The functional analysis of miRNAs has revealed roles in phytohormone homeostasis and spatio-temporal regulation, as well as the mobility of miRNAs and their functions in shoot to root signalling that affects diverse functions, including bacterial entry, meristem division and differentiation, nitrogen fixation and senescence. Furthermore, small RNA fragments of rhizobial origin repress complementary plant mRNAs. We also consider the roles of miRNAs in determinate or indeterminate nodules. Taken together, this overview confirms that miRNAs are master regulators of the legume-rhizobia symbiosis.

GmPAP12 Is Required for Nodule Development and Nitrogen Fixation Under Phosphorus Starvation in Soybean

Nodulation process in legume plants is essential for biological nitrogen fixation during which process a large amount of phosphorus (P) is required. Under P deficiency, nodule formation is greatly affected, and induction of purple acid phosphatases (PAPs) is an adaptive strategy for nodules to acquire more P. However, regulation roles of PAPs in nodules remain largely understood. In this study, by transcriptome sequencing technology, five PAP genes were found to be differentially expressed, which led to the greatly increased acid phosphatase (APase) and phytase activities in soybean mature nodules under P starvation conditions; and among the five PAP genes, GmPAP12 had the highest transcript level, and RT-PCR indicated expression of GmPAP12 was gradually increasing during nodule development. GUS activity driven by GmPAP12 promoter was also significantly induced in low phosphorus conditions. Further functional analysis showed that under low phosphorus stress, overexpression of GmPAP12 resulted in higher nodule number, fresh weight, and nitrogenase activity as well as the APase activity than those of control plant nodules, whereas the growth performance and APase activity of nodules on hairy roots were greatly lower when GmPAP12 was suppressed, indicating that GmPAP12 may promote P utilization in soybean nodules under low P stress, which thus played an important role in nodulation and biological nitrogen fixation. Moreover, P1BS elements were found in the promoter of GmPAP12, and yeast one-hybrid experiment further proved the binding of P1BS by transcription factor GmPHR1 in the promoter of GmPAP12. At last, overexpression and suppression of GmPHR1 in nodules indeed caused highly increased and decreased expression of GmPAP12, respectively, indicating that GmPAP12 is regulated by GmPHR1 in soybean nodules. Taken together, these data suggested that GmPAP12 was a novel soybean PAP involved in the P utilization and metabolism in soybean root nodules and played an important role in the growth and development of root nodules and biological nitrogen fixation.

Small RNA and Transcriptome Sequencing Reveals miRNA Regulation of Floral Thermogenesis in Nelumbo nucifera

The sacred lotus (Nelumbo nucifera Gaertn.) can produce heat autonomously and maintain a relatively stable floral chamber temperature for several days when blooming. Floral thermogenesis is critical for flower organ development and reproductive success. However, the regulatory role of microRNA (miRNA) underlying floral thermogenesis in N. nucifera remains unclear. To comprehensively understand the miRNA regulatory mechanism of thermogenesis, we performed small RNA sequencing and transcriptome sequencing on receptacles from five different developmental stages. In the present study, a total of 172 known miRNAs belonging to 39 miRNA families and 126 novel miRNAs were identified. Twenty-nine thermogenesis-related miRNAs and 3024 thermogenesis-related mRNAs were screened based on their expression patterns. Of those, seventeen differentially expressed miRNAs (DEMs) and 1765 differentially expressed genes (DEGs) had higher expression during thermogenic stages. The upregulated genes in the thermogenic stages were mainly associated with mitochondrial function, oxidoreductase activity, and the energy metabolism process. Further analysis showed that miR156_2, miR395a_5, miR481d, and miR319p may play an important role in heat-producing activity by regulating cellular respiration-related genes. This study provides comprehensive miRNA and mRNA expression profile of receptacle during thermogenesis in N. nucifera, which advances our understanding on the regulation of floral thermogenesis mediated by miRNA.

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.

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.

Transcriptome-wide identification of miRNA targets and a TAS3-homologous gene in Populus by degradome sequencing

BACKGROUND

Degradome sequencing has been applied to identify miRNA-directed mRNA cleavage and understand the biological function of miRNAs and their target genes in plants defense to stress. miRNAs involved in the response to cold stress have been identified in Populus, however, there are few reports about the validated targets of miRNAs in Populus under cold stress.

OBJECTIVES

The primary objective of this investigation was to globally identify and validate the targets of the miRNAs and regulatory components in Populus under cold stress.

METHODS

Populus plantlets grown in vitro were treated with cold (4 °C for 8 h) and total RNA was extracted using Trizol reagent. Approximately 200 µg total RNA was used for the construction of the degradome library, and degradome sequencing was conducted on an Illumina HiSeq 2000. The sequences were mapped to Populus genome using SOAP 2.0 and then were collected for degradome analysis. Additionally, trans-acting siRNA sequences from transacting siRNA gene 3 sequences and mature miRNAs cleaved from precursor miRNAs of Populus were analyzed. 5' RNA ligase-mediated-RACE (5'-RACE) were further conducted.

RESULTS

80 genes were experimentally determined to be the target of 51 unique miRNAs, including three down-regulated miRNAs (pto-miR156k, pto-miR169i-m, and pto-miR394a-5p/b-5p) and two up-regulated miRNAs (pto-miR167a-d and pto-miR167f/g). The specificity and diversity of cleavage sites of miRNA targets were validated through 5'-RACE experiment and the results were similar with that of degradome sequencing, further supporting the empirical cleavage of miRNAs on targets in vivo in Populus. Interestingly, the TAS-homologous gene pto-TAS3 (EF146176.1) was identified and 11 potential ta-siRNAs [D1(+)-D11(+)] and their possible biogenesis sites within the pto-TAS3 transcript sequence were predicted in Populus. In addition, the biosynthesis of miRNA from precursor miRNA (pre-miRNA) was also validated through the detection of a total of 17 pre-miRNAs.

CONCLUSION

Our investigation expands the application of degradome sequencing for evaluating miRNA regulatory elements and evidence of the miRNA synthesis process, and provides empirical evidence of bona fide cleavage of target genes by miRNAs in Populus, which might be used for the research of miRNA-mediated regulation mechanism and molecular improvement of resistance to cold stress.

Legume nodulation: The host controls the party

Global demand to increase food production and simultaneously reduce synthetic nitrogen fertilizer inputs in agriculture are underpinning the need to intensify the use of legume crops. The symbiotic relationship that legume plants establish with nitrogen-fixing rhizobia bacteria is central to their advantage. This plant-microbe interaction results in newly developed root organs, called nodules, where the rhizobia convert atmospheric nitrogen gas into forms of nitrogen the plant can use. However, the process of developing and maintaining nodules is resource intensive; hence, the plant tightly controls the number of nodules forming. A variety of molecular mechanisms are used to regulate nodule numbers under both favourable and stressful growing conditions, enabling the plant to conserve resources and optimize development in response to a range of circumstances. Using genetic and genomic approaches, many components acting in the regulation of nodulation have now been identified. Discovering and functionally characterizing these components can provide genetic targets and polymorphic markers that aid in the selection of superior legume cultivars and rhizobia strains that benefit agricultural sustainability and food security. This review addresses recent findings in nodulation control, presents detailed models of the molecular mechanisms driving these processes, and identifies gaps in these processes that are not yet fully explained.

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.

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.

sRNAome and transcriptome analysis provide insight into chilling response of cowpea pods

Cowpea is an important horticultural crop in tropical and subtropical areas of Asia, Africa, and Latin America, as well as parts of southern Europe and Central and South America. Chilling injury is a common physiological hazard of cowpea in cold chain logistics which reduce the cowpea pods nutritional quality and product value. However, the molecular mechanism involved in chilling injury remains unclear in cowpea pods. RNA-Seq and sRNA-Seq technologies were employed to decipher the miRNAs and mRNAs expression profiles and their regulatory networks in cowpea pods involved in chilling stress. Differentially expressed miRNAs and mRNA profiles were obtained based on cluster analysis, miRNAs and target genes were found to show coherent relationships in the regulatory networks of chilling injury. Furthermore, we found that numerous miRNAs and nat-siRNAs' targets were predicted to be key enzymes involved in the redox reactions such as POD, CAT, AO and LOX, energy metabolism such as ATPase, FAD and NAD related enzymes and different transcription factors such as WRKY, bHLH, MYB, ERF and NAC which play important roles in chilling injury.

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.

Regulation of Isoflavone Biosynthesis by miRNAs in Two Contrasting Soybean Genotypes at Different Seed Developmental Stages

Owing to the presence of nutritionally important, health-promoting bioactive compounds, especially isoflavones, soybean has acquired the status of a functional food. miRNAs are tiny riboregulator of gene expression by either decreasing and/or increasing the expression of their corresponding target genes. Despite several works on identification and functional characterization of plant miRNAs, the role of miRNAs in the regulation of isoflavones metabolism is still a virgin field. In the present study, we identified a total of 31 new miRNAs along with their 245 putative target genes from soybean seed-specific ESTs using computational approach. The Kyoto Encyclopedia of Genes and Genomes pathway analyses indicated that miRNA putatively regulates metabolism and genetic information processing. Out of that, a total of 5 miRNAs (Gma-miRNA12, Gma-miRNA24, Gma-miRNA26, Gma-miRNA28, and Gma-miRNA29) were predicted and validated for their probable role during isoflavone biosynthesis. We also validated their five target genes using RA-PCR, which is as good as 5'RLM-RACE. Temporal regulation [35 days after flowering, 45, 55, and 65 DAF] of miRNAs and their targets showed differential expression schema. Differential expression of Gma-miR26 and Gma-miRNA28 along with their corresponding target genes (Glyma.10G197900 and Glyma.09G127200) showed a direct relationship with the total isoflavone content. Therefore, understanding the miRNA-based genetic regulation of isoflavone pathway would assist in selection and manipulation to get high-performing soybean genotypes with better isoflavone yield.

Using genomic information to improve soybean adaptability to climate change

Climate change has brought severe challenges to agriculture. It is anticipated that there will be a drop in crop yield - including that of soybean - due to climatic stress factors that include drastic fluctuations in temperature, drought, flooding and high salinity. Genomic information on soybean has been accumulating rapidly since initial publication of its reference genome, providing a valuable tool for the improvement of cultivated soybean. Not only are many molecular markers that are associated with important quantitative trait loci now identified, but we also have a more detailed picture of the genomic variations among soybean germplasms, enabling us to utilize these as tools to assist crop breeding. In this review, we will summarize and discuss the currently available soybean genomic approaches, including whole-genome sequencing, sequencing-based genotyping, functional genomics, proteomics, and epigenomics. The information uncovered through these techniques will help further pinpoint important gene candidates and genomic loci associated with adaptive traits, as well as achieving a better understanding of how soybeans cope with the changing climate.

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.

Breeding approaches and genomics technologies to increase crop yield under low-temperature stress

Improved knowledge about plant cold stress tolerance offered by modern omics technologies will greatly inform future crop improvement strategies that aim to breed cultivars yielding substantially high under low-temperature conditions. Alarmingly rising temperature extremities present a substantial impediment to the projected target of 70% more food production by 2050. Low-temperature (LT) stress severely constrains crop production worldwide, thereby demanding an urgent yet sustainable solution. Considerable research progress has been achieved on this front. Here, we review the crucial cellular and metabolic alterations in plants that follow LT stress along with the signal transduction and the regulatory network describing the plant cold tolerance. The significance of plant genetic resources to expand the genetic base of breeding programmes with regard to cold tolerance is highlighted. Also, the genetic architecture of cold tolerance trait as elucidated by conventional QTL mapping and genome-wide association mapping is described. Further, global expression profiling techniques including RNA-Seq along with diverse omics platforms are briefly discussed to better understand the underlying mechanism and prioritize the candidate gene (s) for downstream applications. These latest additions to breeders' toolbox hold immense potential to support plant breeding schemes that seek development of LT-tolerant cultivars. High-yielding cultivars endowed with greater cold tolerance are urgently required to sustain the crop yield under conditions severely challenged by low-temperature.

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

Chilling stress is a major factor limiting the yield and quality of vegetable soybean (Glycine max L.) on a global scale. In the present study, systematic identification and functional analysis of miRNAs under chilling stress were carried out to clarify the molecular mechanism of chilling resistance. Two independent small RNA libraries from leaves of soybean were constructed and sequenced with the high-throughput Illumina Solexa system. A total of 434 known miRNAs and 3 novel miRNAs were identified. Thirty-five miRNAs were verified by qRT-PCR analysis. Furthermore, their gene targets were identified via high-throughput degradome sequencing. A total of 898 transcripts were targeted by 54 miRNA families attributed to five categories. More importantly, we identified 51 miRNAs differentially expressed between chilling stress and control conditions. The targets of these miRNAs were enriched in oxidation-reduction, signal transduction, and metabolic process functional categories. Our qRT-PCR analysis confirmed a negative relationship among the miRNAs and their targets under chilling stress. Our work thus provides comprehensive molecular evidence supporting the involvement of miRNAs in chilling-stress responses in vegetable soybean.

Genome-Wide Investigation of MicroRNAs and Their Targets in Response to Freezing Stress in Medicago sativa L., Based on High-Throughput Sequencing

Winter damage, especially in northern climates, is a major limitation of the utilization of perennial forages such as alfalfa. Therefore, improving freezing tolerance is imperative in alfalfa genetic breeding. However, freezing tolerance is a complex trait that is determined by many genes. To understand the complex regulation mechanisms of freezing tolerance in alfalfa, we performed small RNA sequencing analysis under cold (4°) and freezing (-8°) stress. The sequencing results revealed that 173 known, and 24 novel miRNAs were expressed, and that the expression of 35 miRNAs was affected by cold and/or freezing stress. Meanwhile, 105 target genes cleaved by these miRNAs were characterized by degradome sequencing. These targets were associated with biological regulation, cellular processes, metabolic processes, and response to stress. Interestingly, most of them were characterized as transcription factors (TFs), including auxin response factors, SBP, NAC, AP2/ERF, and GRF, which play important roles in plant abiotic responses. In addition, important miRNAs and mRNAs involved in nodulation were also identified, for example, the relationship between miR169 and the TF CCAAT (also named as NF-YA/HAP2), which suggested that nodulation has an important function in freezing tolerance in alfalfa. Our results provide valuable information to help determine the molecular mechanisms of freezing tolerance in alfalfa, which will aid the application of these miRNAs and their targets in the improvement of freezing tolerance in alfalfa and related plants.

Genome-wide identification of the Phaseolus vulgaris sRNAome using small RNA and degradome sequencing

BACKGROUND

MiRNAs and phasiRNAs are negative regulators of gene expression. These small RNAs have been extensively studied in plant model species but only 10 mature microRNAs are present in miRBase version 21, the most used miRNA database, and no phasiRNAs have been identified for the model legume Phaseolus vulgaris. Thanks to the recent availability of the first version of the common bean genome, degradome data and small RNA libraries, we are able to present here a catalog of the microRNAs and phasiRNAs for this organism and, particularly, we suggest new protagonists in the symbiotic nodulation events.

RESULTS

We identified a set of 185 mature miRNAs, including 121 previously unpublished sequences, encoded by 307 precursors and distributed in 98 families. Degradome data allowed us to identify a total of 181 targets for these miRNAs. We reveal two regulatory networks involving conserved miRNAs: those known to play crucial roles in the establishment of nodules, and novel miRNAs present only in common bean, suggesting a specific role for these sequences. In addition, we identified 125 loci that potentially produce phased small RNAs, with 47 of them having all the characteristics of being triggered by a total of 31 miRNAs, including 14 new miRNAs identified in this study.

CONCLUSIONS

We provide here a set of new small RNAs that contribute to the broader knowledge of the sRNAome of Phaseolus vulgaris. Thanks to the identification of the miRNA targets from degradome analysis and the construction of regulatory networks between the mature microRNAs, we present here the probable functional regulation associated with the sRNAome and, particularly, in N2-fixing symbiotic nodules.

Symbiotic nitrogen fixation in legume nodules: metabolism and regulatory mechanisms

The special issue “Symbiotic Nitrogen Fixation in Legume Nodules: Metabolism and Regulatory Mechanisms” aims to investigate the physiological and biochemical advances in the symbiotic process with an emphasis on nodule establishment, development and functioning. The original research articles included in this issue provide important information regarding novel aspects of nodule metabolism and various regulatory pathways, which could have important future implications. This issue also included one review article that highlights the importance of using legume trees in the production of renewable biofuels.