Identification of microRNAs and their mRNA targets during soybean nodule development: functional analysis of the role of miR393j-3p in soybean nodulation

GmBES1-1 dampens the activity of GmNSP1/2 to mediate brassinosteroid inhibition of nodulation in soybean

Soybean (Glycine max) forms root nodules to house rhizobial bacteria for biological nitrogen fixation. The development of root nodules is intricately regulated by endogenous and exogenous cues. The phytohormones brassinosteroids (BRs) have been shown to negatively regulate nodulation in soybean, but the underlying genetic and molecular mechanisms remain largely unknown. Here, we performed transcriptomic analyses and revealed that BR signaling negatively regulates nodulation factor (NF) signaling. We found that BR signaling inhibits nodulation through its signaling component GmBES1-1 by dampening NF signaling and nodule formation. In addition, GmBES1-1 can directly interact with both GmNSP1 and GmNSP2 to inhibit their interaction and the DNA-binding activity of GmNSP1. Furthermore, BR-induced nuclear accumulation of GmBES1-1 is essential for inhibiting nodulation. Taken together, our results demonstrate that regulation of GmBES1-1 subcellular localization by BRs plays a key role in the legume-rhizobium symbiosis and plant development, indicating a crosstalk mechanism between phytohormone and symbiosis signaling pathways.

MicroRNA 4407 modulates nodulation in soybean by repressing a root-specific ISOPENTENYLTRANSFERASE (GmIPT3)

MicroRNAs (miRNAs) are important regulators of plant biological processes, including soybean nodulation. One miRNA, miR4407, was identified in soybean roots and nodules. However, the function of miR4407 in soybean is still unknown. MiR4407, unique to soybean, positively regulates lateral root emergence and root structures and represses a root-specific ISOPENTENYLTRANSFERASE (GmIPT3). By altering the expression of miR4407 and GmIPT3, we investigated the role of miR4407 in lateral root and nodule development. Both miR4407 and GmIPT3 are expressed in the inner root cortex and nodule primordia. Upon rhizobial inoculation, miR4407 was downregulated while GmIPT3 was upregulated. Overexpressing miR4407 reduced the number of nodules in transgenic soybean hairy roots while overexpressing the wild-type GmIPT3 or a miR4407-resistant GmIPT3 mutant (mGmIPT3) significantly increased the nodule number. The mechanism of miR4407 and GmIPT3 functions was also linked to autoregulation of nodulation (AON), where miR4407 overexpression repressed miR172c and activated its target, GmNNC1, turning on AON. Exogenous CK mimicked the effects of GmIPT3 overexpression on miR172c, supporting the notion that GmIPT3 regulates nodulation by enhancing root-derived CK. Overall, our data revealed a new miRNA-mediated regulatory mechanism of nodulation in soybean. MiR4407 showed a dual role in lateral root and nodule development.

MsSPL9 Modulates Nodulation under Nitrate Sufficiency Condition in Medicago sativa

Nodulation in Leguminous spp. is induced by common environmental cues, such as low nitrogen availability conditions, in the presence of the specific Rhizobium spp. in the rhizosphere. Medicago sativa (alfalfa) is an important nitrogen-fixing forage crop that is widely cultivated around the world and relied upon as a staple source of forage in livestock feed. Although alfalfa's relationship with these bacteria is one of the most efficient between rhizobia and legume plants, breeding for nitrogen-related traits in this crop has received little attention. In this report, we investigate the role of Squamosa-Promoter Binding Protein-Like 9 (SPL9), a target of miR156, in nodulation in alfalfa. Transgenic alfalfa plants with SPL9-silenced (SPL9-RNAi) and overexpressed (35S::SPL9) were compared to wild-type (WT) alfalfa for phenotypic changes in nodulation in the presence and absence of nitrogen. Phenotypic analyses showed that silencing of MsSPL9 in alfalfa caused an increase in the number of nodules. Moreover, the characterization of phenotypic and molecular parameters revealed that MsSPL9 regulates nodulation under a high concentration of nitrate (10 mM KNO₃) by regulating the transcription levels of the nitrate-responsive genes Nitrate Reductase1 (NR1), NR2, Nitrate transporter 2.5 (NRT2.5), and a shoot-controlled autoregulation of nodulation (AON) gene, Super numeric nodules (SUNN). While MsSPL9-overexpressing transgenic plants have dramatically increased transcript levels of SUNN, NR1, NR2, and NRT2.5, reducing MsSPL9 caused downregulation of these genes and displayed a nitrogen-starved phenotype, as downregulation of the MsSPL9 transcript levels caused a nitrate-tolerant nodulation phenotype. Taken together, our results suggest that MsSPL9 regulates nodulation in alfalfa in response to nitrate.

Plants Recruit Peptides and Micro RNAs to Regulate Nutrient Acquisition from Soil and Symbiosis

Plants engage in symbiotic relationships with soil microorganisms to overcome nutrient limitations in their environment. Among the best studied endosymbiotic interactions in plants are those with arbuscular mycorrhizal (AM) fungi and N-fixing bacteria called rhizobia. The mechanisms regulating plant nutrient homeostasis and acquisition involve small mobile molecules such as peptides and micro RNAs (miRNAs). A large number of CLE (CLAVATA3/EMBRYO SURROUNDING REGION-RELATED) and CEP (C-TERMINALLY ENCODED PEPTIDE) peptide hormones as well as certain miRNAs have been reported to differentially respond to the availability of essential nutrients such as nitrogen (N) and phosphorus (P). Interestingly, a partially overlapping pool of these molecules is involved in plant responses to root colonization by rhizobia and AM fungi, as well as mineral nutrition. The crosstalk between root endosymbiosis and nutrient availability has been subject of intense investigations, and new insights in locally or systemically mobile molecules in nutrient- as well as symbiosis-related signaling continue to arise. Focusing on the key roles of peptides and miRNAs, we review the mechanisms that shape plant responses to nutrient limitation and regulate the establishment of symbiotic associations with beneficial soil microorganisms.

GmYSL7 controls iron uptake, allocation, and cellular response of nodules in soybean

Iron (Fe) is essential for DNA synthesis, photosynthesis and respiration of plants. The demand for Fe substantially increases during legumes-rhizobia symbiotic nitrogen fixation because of the synthesis of leghemoglobin in the host and Fe-containing proteins in bacteroids. However, the mechanism by which plant controls iron transport to nodules remains largely unknown. Here we demonstrate that GmYSL7 serves as a key regulator controlling Fe uptake from root to nodule and distribution in soybean nodules. GmYSL7 is Fe responsive and GmYSL7 transports iron across the membrane and into the infected cells of nodules. Alterations of GmYSL7 substantially affect iron distribution between root and nodule, resulting in defective growth of nodules and reduced nitrogenase activity. GmYSL7 knockout increases the expression of GmbHLH300, a transcription factor required for Fe response of nodules. Overexpression of GmbHLH300 decreases nodule number, nitrogenase activity and Fe content in nodules. Remarkably, GmbHLH300 directly binds to the promoters of ENOD93 and GmLbs, which regulate nodule number and nitrogenase activity, and represses their transcription. Our data reveal a new role of GmYSL7 in controlling Fe transport from host root to nodule and Fe distribution in nodule cells, and uncover a molecular mechanism by which Fe affects nodule number and nitrogenase activity.

Roles of non-coding RNAs in the hormonal and nutritional regulation in nodulation and nitrogen fixation

Symbiotic nitrogen fixation is an important component in the nitrogen cycle and is a potential solution for sustainable agriculture. It is the result of the interactions between the plant host, mostly restricted to legume species, and the rhizobial symbiont. From the first encounter between the host and the symbiont to eventual successful nitrogen fixation, there are delicate processes involved, such as nodule organogenesis, rhizobial infection thread progression, differentiation of the bacteroid, deregulation of the host defense systems, and reallocation of resources. All these processes are tightly regulated at different levels. Recent evidence revealed that non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs), participate in these processes by controlling the transcription and translation of effector genes. In general, ncRNAs are functional transcripts without translation potential and are important gene regulators. MiRNAs, negative gene regulators, bind to the target mRNAs and repress protein production by causing the cleavage of mRNA and translational silencing. LncRNAs affect the formation of chromosomal loops, DNA methylation, histone modification, and alternative splicing to modulate gene expression. Both lncRNAs and circRNAs could serve as target mimics of miRNA to inhibit miRNA functions. In this review, we summarized and discussed the current understanding of the roles of ncRNAs in legume nodulation and nitrogen fixation in the root nodule, mainly focusing on their regulation of hormone signal transduction, the autoregulation of nodulation (AON) pathway and nutrient homeostasis in nodules. Unraveling the mediation of legume nodulation by ncRNAs will give us new insights into designing higher-performance leguminous crops for sustainable agriculture.

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.

Large-Scale Integrative Analysis of Soybean Transcriptome Using an Unsupervised Autoencoder Model

Plant tissues are distinguished by their gene expression patterns, which can help identify tissue-specific highly expressed genes and their differential functional modules. For this purpose, large-scale soybean transcriptome samples were collected and processed starting from raw sequencing reads in a uniform analysis pipeline. To address the gene expression heterogeneity in different tissues, we utilized an adversarial deconfounding autoencoder (AD-AE) model to map gene expressions into a latent space and adapted a standard unsupervised autoencoder (AE) model to help effectively extract meaningful biological signals from the noisy data. As a result, four groups of 1,743, 914, 2,107, and 1,451 genes were found highly expressed specifically in leaf, root, seed and nodule tissues, respectively. To obtain key transcription factors (TFs), hub genes and their functional modules in each tissue, we constructed tissue-specific gene regulatory networks (GRNs), and differential correlation networks by using corrected and compressed gene expression data. We validated our results from the literature and gene enrichment analysis, which confirmed many identified tissue-specific genes. Our study represents the largest gene expression analysis in soybean tissues to date. It provides valuable targets for tissue-specific research and helps uncover broader biological patterns. Code is publicly available with open source at https://github.com/LingtaoSu/SoyMeta.

Progress in soybean functional genomics over the past decade

Soybean is one of the most important oilseed and fodder crops. Benefiting from the efforts of soybean breeders and the development of breeding technology, large number of germplasm has been generated over the last 100 years. Nevertheless, soybean breeding needs to be accelerated to meet the needs of a growing world population, to promote sustainable agriculture and to address future environmental changes. The acceleration is highly reliant on the discoveries in gene functional studies. The release of the reference soybean genome in 2010 has significantly facilitated the advance in soybean functional genomics. Here, we review the research progress in soybean omics (genomics, transcriptomics, epigenomics and proteomics), germplasm development (germplasm resources and databases), gene discovery (genes that are responsible for important soybean traits including yield, flowering and maturity, seed quality, stress resistance, nodulation and domestication) and transformation technology during the past decade. At the end, we also briefly discuss current challenges and future directions.

Non-Coding RNAs in Legumes: Their Emerging Roles in Regulating Biotic/Abiotic Stress Responses and Plant Growth and Development

Noncoding RNAs, including microRNAs (miRNAs), small interference RNAs (siRNAs), circular RNA (circRNA), and long noncoding RNAs (lncRNAs), control gene expression at the transcription, post-transcription, and translation levels. Apart from protein-coding genes, accumulating evidence supports ncRNAs playing a critical role in shaping plant growth and development and biotic and abiotic stress responses in various species, including legume crops. Noncoding RNAs (ncRNAs) interact with DNA, RNA, and proteins, modulating their target genes. However, the regulatory mechanisms controlling these cellular processes are not well understood. Here, we discuss the features of various ncRNAs, including their emerging role in contributing to biotic/abiotic stress response and plant growth and development, in addition to the molecular mechanisms involved, focusing on legume crops. Unravelling the underlying molecular mechanisms and functional implications of ncRNAs will enhance our understanding of the coordinated regulation of plant defences against various biotic and abiotic stresses and for key growth and development processes to better design various legume crops for global food security.

Differentially expressed microRNAs that target functional genes in mature soybean nodules

MicroRNAs (miRNAs) are important regulators of biological functions in plants. To find out what roles miRNAs play in regulating symbiotic nitrogen fixation (SNF) in soybean [Glycine max (L.) Merr.], we identified high-confidence differentially expressed (DE) miRNAs from uninoculated roots (UR), rhizobium-inoculated roots (IR), and nodules (NODs) of soybean by robust small RNA sequencing (sRNA-seq). Based on their predicted target messenger RNAs (mRNAs), the expression profiles of some of these DE miRNAs could be linked to nodule functions. In particular, several miRNAs associated with nutrient transportation genes were differentially expressed in IRs and mature NODs. MiR399b, specifically, was highly induced in IRs and NODs, as well as by inorganic phosphate (Pi) starvation. In composite soybean plants overexpressing miR399b, PHOSPHATE2 (PHO2), a known target of miR399b that inhibits the activities of high-affinity Pi transporters, was strongly repressed. In addition, the overexpression of miR399b in the roots of transgenic composite plants significantly improved whole-plant Pi and ureide concentrations and the overall growth in terms of leaf node numbers and whole-plant dry weight. Our findings suggest that the induction of miR399b in NODs could enhance nitrogen fixation and soybean growth, possibly via improving Pi uptake to achieve a better Pi-nitrogen balance to promote SNF in nodules.

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.

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.

miR169c-NFYA-C-ENOD40 modulates nitrogen inhibitory effects in soybean nodulation

Legume crops contribute a great portion of clean nitrogen (N) to agro-ecosystems through symbiotic N2 fixation in the nodule; however, the nodulation is always inhibited by high N availability which is known as the N inhibitory effect through largely unknown mechanisms. We functionally investigated miR169c-GmNFYA-C-GmENOD40 under multiple N conditions in soybean (Glycine max) (ENOD, Early Nodulin; NFYA, Nuclear Factor-Y Subunit A). We elucidated their regulatory roles in soybean nodulation through analyzing expression patterns, micro-messenger RNA (miRNA-mRNA) interactions, phenotypes of transgenic soybean plants and genetic interactions. We found that miR169c expression was induced by high N, whereas its target GmNFYA-C was preferentially expressed in nodules and induced by rhizobium inoculation. Overexpression of miR169c inhibited nodulation through targeting 3'-UTR of GmNFYA-C, whereas knockout miR169c through CRISPR-cas9 promoted nodulation. However, overexpression of GmNFYA-C promoted soybean nodulation through facilitating rhizobium infection and increasing the expression of symbiotic signaling gene GmENOD40. Besides, GmNFYA-C directly induced the expression of GmENOD40. In addition, overexpression of GmNFYA-C without the target site of miR169c partially attenuated the inhibitory effect of high N on soybean nodulation. We discovered a new regulatory pathway involving the miR169c-NFYA-C-ENOD40 module that regulates soybean nodulation in response to N availability. This pathway provides substantial new insights into the mechanisms underlying the N inhibitory effect on nodulation.

Origin and Evolutionary Dynamics of the miR2119 and ADH1 Regulatory Module in Legumes

MicroRNAs are important regulators of gene expression in eukaryotes. Previously, we reported that in Phaseolus vulgaris, the precursor for miR2119 is located in the same gene as miR398a, conceiving a dicistronic MIR gene. Both miRNA precursors are transcribed and processed from a single transcript resulting in two mature microRNAs that regulate the mRNAs encoding ALCOHOL DEHYDROGENASE 1 (ADH1) and COPPER-ZINC SUPEROXIDE DISMUTASE 1 (CSD1). Genes for miR398 are distributed throughout the spermatophytes; however, miR2119 is only found in Leguminosae species, indicating its recent emergence. Here, we used public databases to explore the presence of the miR2119 sequence in several plant species. We found that miR2119 is present only in specific clades within the Papilionoideae subfamily, including important crops used for human consumption and forage. Within this subfamily, MIR2119 and MIR398a are found together as a single gene in the genomes of the Millettioids and Hologalegina. In contrast, in the Dalbergioids MIR2119 is located in a different locus from MIR398a, suggesting this as the ancestral genomic organization. To our knowledge, this is a unique example where two separate MIRNA genes have merged to generate a single polycistronic gene. Phylogenetic analysis of ADH1 gene sequences in the Papilionoideae subfamily revealed duplication events resulting in up to four ADH1 genes in certain species. Notably, the presence of MIR2119 correlates with the conservation of target sites in particular ADH1 genes in each clade. Our results suggest that post-transcriptional regulation of ADH1 genes by miR2119 has contributed to shaping the expansion and divergence of this gene family in the Papilionoideae. Future experimental work on ADH1 regulation by miR2119 in more legume species will help to further understand the evolutionary history of the ADH1 gene family and the relevance of miRNA regulation in this process.

Knock-down of δ-aminolevulinic acid dehydratase via virus-induced gene silencing alters the microRNA biogenesis and causes stress-related reactions in citrus plants

The δ-aminolevulinic acid (δ-ALA) is an intermediate in the biosynthetic pathway of tetrapyrroles. Tetrapyrroles play vital roles in many biological processes such as photosynthesis, respiration, and light-sensing. ALA-dehydratase (ALAD) combines two molecules of δ-ALA to form porphobilinogen. In citrus, the silencing of ALAD caused discrete yellow spots and necrosis in leaves and stems. Additionally, it caused rapid death in developing new shoots. Herein, we hypothesize that the accumulation of δ-ALA results in severe stress and reduced meristem development. For that reason, we investigated the dynamic changes in the expression profiles of 23 microRNA (miRNA) identified through small RNA sequencing, from CTV-tALAD plants in comparison with healthy C. macrophylla and C. macrophylla infiltrated with CTV-wt. Furthermore, we reported the effect of ALAD silencing on the total phenolics, H₂O_2, and reactive oxygen species (ROS) levels, to examine the possibilities of miRNAs involving the regulation of these pathways. Our results showed that the total phenolics content, H₂O₂, and O₂^- levels were increased in CTV-tALAD plants. Moreover, 63 conserved miRNA members belonging to 23 different miRNA families were differentially expressed in CTV-tALAD plants compared to controls. The identified miRNAs are implicated in auxin biosynthesis and signaling, axillary shoot meristem formation and leaf morphology, starch metabolism, and oxidative stress. Collectively, our findings suggested that ALAD silencing initiates stress on citrus plants. As a result, CTV-tALAD plants exhibit reduced metabolic rate, growth, and development in order to cope with the stress that resulted from the accumulation of δ-ALA. This cascade of events led to leaf, stem, and meristem necrosis and failure of new shoot development.

Interactions of gene expression, alternative splicing, and DNA methylation in determining nodule identity

Soybean nodulation is a highly controlled process that involves complex gene regulation at both transcriptional and post-transcriptional levels. In the present study, we profiled gene expression changes, alternative splicing events, and DNA methylation patterns during nodule formation, development, and senescence. The transcriptome data uncovered key transcription patterns of nodule development that included 9669 core genes and 7302 stage-specific genes. Alternative splicing analysis uncovered a total of 2323 genes that undergo alternative splicing events in at least one nodule developmental stage, with activation of exon skipping and repression of intron retention being the most common splicing events in nodules compared to roots. Approximately 40% of the differentially spliced genes were also differentially expressed at the same nodule developmental stage, implying a substantial association between gene expression and alternative splicing. Genome-wide-DNA methylation analysis revealed dynamic changes in nodule methylomes that were specific to each nodule stage, occurred in a sequence-specific manner, and impacted the expression of 1864 genes. An attractive hypothesis raised by our data is that increased DNA methylation may contribute to the efficiency of alternative splicing. Together, our results provide intriguing insights into the associations between gene expression, alternative splicing, and DNA methylation that may shape transcriptome complexity and proteome specificity in developing soybean nodules.

The role of microRNAs in the legume-Rhizobium nitrogen-fixing symbiosis

Under nitrogen starvation, most legume plants form a nitrogen-fixing symbiosis with Rhizobium bacteria. The bacteria induce the formation of a novel organ called the nodule in which rhizobia reside as intracellular symbionts and convert atmospheric nitrogen into ammonia. During this symbiosis, miRNAs are essential for coordinating the various plant processes required for nodule formation and function. miRNAs are non-coding, endogenous RNA molecules, typically 20-24 nucleotides long, that negatively regulate the expression of their target mRNAs. Some miRNAs can move systemically within plant tissues through the vascular system, which mediates, for example, communication between the stem/leaf tissues and the roots. In this review, we summarize the growing number of miRNAs that function during legume nodulation focusing on two model legumes, Lotus japonicus and Medicago truncatula, and two important legume crops, soybean (Glycine max) and common bean (Phaseolus vulgaris). This regulation impacts a variety of physiological processes including hormone signaling and spatial regulation of gene expression. The role of mobile miRNAs in regulating legume nodule number is also highlighted.

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.

MicroRNAs and Their Regulatory Roles in Plant-Environment Interactions

MicroRNAs (miRNAs) are 20-24 nucleotide noncoding RNAs abundant in plants and animals. The biogenesis of plant miRNAs involves transcription of miRNA genes, processing of primary miRNA transcripts by DICER-LIKE proteins into mature miRNAs, and loading of mature miRNAs into ARGONAUTE proteins to form miRNA-induced silencing complex (miRISC). By targeting complementary sequences, miRISC negatively regulates gene expression, thereby coordinating plant development and plant-environment interactions. In this review, we present and discuss recent updates on the mechanisms and regulation of miRNA biogenesis, miRISC assembly and actions as well as the regulatory roles of miRNAs in plant developmental plasticity, abiotic/biotic responses, and symbiotic/parasitic interactions. Finally, we suggest future directions for plant miRNA research.

Characterization of the Spatial and Temporal Expression of Two Soybean miRNAs Identifies SCL6 as a Novel Regulator of Soybean Nodulation

MicroRNAs (miRNAs) control expression of endogenous target genes through transcript cleavage or translational inhibition. Legume plants can form a specialized organ, the nodule, through interaction with nitrogen fixing soil bacteria. To understand the regulatory roles of miRNAs in the nodulation process, we functionally validated gma-miR171o and gma-miR171q and their target genes in soybean. These two miRNA sequences are identical in sequence but their miRNA genes are divergent and show unique, tissue-specific expression patterns. The expression levels of the two miRNAs are negatively correlated with that of their target genes. Ectopic expression of these miRNAs in transgenic hairy roots resulted in a significant reduction in nodule formation. Both gma-miR171o and gma-miR171q target members of the GRAS transcription factor superfamily, namely GmSCL-6 and GmNSP2. Transient interaction of miRNAs and their target genes in tobacco cells further confirmed their cleavage activity. The results suggest that gma-miR171o and gma-miR171q regulate GmSCL-6 and GmNSP2, which in turn, influence expression of the Nodule inception (NIN), Early Nodulin 40 (ENOD40), and Ethylene Response Factor Required for Nodulation (ERN) genes during the Bradyrhizobium japonicum-soybean nodulation process. Collectively, our data suggest a role for two miRNAs, gma-miR171o and gma-miR171q, in regulating the spatial and temporal aspects of soybean nodulation.

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

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

Transcriptional Reprogramming of Legume Genomes: Perspective and Challenges Associated With Single-Cell and Single Cell-Type Approaches During Nodule Development

Transcriptomic approaches revealed thousands of genes differentially or specifically expressed during nodulation, a biological process resulting from the symbiosis between leguminous plant roots and rhizobia, atmospheric nitrogen-fixing symbiotic bacteria. Ultimately, nodulation will lead to the development of a new root organ, the nodule. Through functional genomic studies, plant transcriptomes have been used by scientists to reveal plant genes potentially controlling nodulation. However, it is important to acknowledge that the physiology, transcriptomic programs, and biochemical properties of the plant cells involved in nodulation are continuously regulated. They also differ between the different cell-types composing the nodules. To generate a more accurate picture of the transcriptome, epigenome, proteome, and metabolome of the cells infected by rhizobia and cells composing the nodule, there is a need to implement plant single-cell and single cell-types strategies and methods. Accessing such information would allow a better understanding of the infection of plant cells by rhizobia and will help understanding the complex interactions existing between rhizobia and the plant cells. In this mini-review, we are reporting the current knowledge on legume nodulation gained by plant scientists at the level of single cell-types, and provide perspectives on single cell/single cell-type approaches when applied to legume nodulation.

Conserved Proteins of the RNA Interference System in the Arbuscular Mycorrhizal Fungus Rhizoglomus irregulare Provide New Insight into the Evolutionary History of Glomeromycota

Horizontal gene transfer (HGT) is an important mechanism in the evolution of many living organisms particularly in Prokaryotes where genes are frequently dispersed between taxa. Although, HGT has been reported in Eukaryotes, its accumulative effect and its frequency has been questioned. Arbuscular mycorrhizal fungi (AMF) are an early diverged fungal lineage belonging to phylum Glomeromycota, whose phylogenetic position is still under debate. The history of AMF and land plant symbiosis dates back to at least 460 Ma. However, Glomeromycota are estimated to have emerged much earlier than land plants. In this study, we surveyed genomic and transcriptomic data of the model arbuscular mycorrhizal fungus Rhizoglomus irregulare (synonym Rhizophagus irregularis) and its relatives to search for evidence of HGT that occurred during AMF evolution. Surprisingly, we found a signature of putative HGT of class I ribonuclease III protein-coding genes that occurred from autotrophic cyanobacteria genomes to R. irregulare. At least one of two HGTs was conserved among AMF species with high levels of sequence similarity. Previously, an example of intimate symbiosis between AM fungus and cyanobacteria was reported in the literature. Ribonuclease III family enzymes are important in small RNA regulation in Fungi together with two additional core proteins (Argonaute/piwi and RdRP). The eukaryotic RNA interference system found in AMF was conserved and showed homology with high sequence similarity in Mucoromycotina, a group of fungi closely related to Glomeromycota. Prior to this analysis, class I ribonuclease III has not been identified in any eukaryotes. Our results indicate that a unique acquisition of class I ribonuclease III in AMF is due to a HGT event that occurred from cyanobacteria to Glomeromycota, at the latest before the divergence of the two Glomeromycota orders Diversisporales and Glomerales.

CRISPR/Cas9 and TALENs generate heritable mutations for genes involved in small RNA processing of Glycine max and Medicago truncatula

Processing of double-stranded RNA precursors into small RNAs is an essential regulator of gene expression in plant development and stress response. Small RNA processing requires the combined activity of a functionally diverse group of molecular components. However, in most of the plant species, there are insufficient mutant resources to functionally characterize each encoding gene. Here, mutations in loci encoding protein machinery involved in small RNA processing in soya bean and Medicago truncatula were generated using the CRISPR/Cas9 and TAL-effector nuclease (TALEN) mutagenesis platforms. An efficient CRISPR/Cas9 reagent was used to create a bi-allelic double mutant for the two soya bean paralogous Double-stranded RNA-binding2 (GmDrb2a and GmDrb2b) genes. These mutations, along with a CRISPR/Cas9-generated mutation of the M. truncatula Hua enhancer1 (MtHen1) gene, were determined to be germ-line transmissible. Furthermore, TALENs were used to generate a mutation within the soya bean Dicer-like2 gene. CRISPR/Cas9 mutagenesis of the soya bean Dicer-like3 gene and the GmHen1a gene was observed in the T0 generation, but these mutations failed to transmit to the T1 generation. The irregular transmission of induced mutations and the corresponding transgenes was investigated by whole-genome sequencing to reveal a spectrum of non-germ-line-targeted mutations and multiple transgene insertion events. Finally, a suite of combinatorial mutant plants were generated by combining the previously reported Gmdcl1a, Gmdcl1b and Gmdcl4b mutants with the Gmdrb2ab double mutant. Altogether, this study demonstrates the synergistic use of different genome engineering platforms to generate a collection of useful mutant plant lines for future study of small RNA processing in legume crops.

Transcription factors network in root endosymbiosis establishment and development

Root endosymbioses are mutualistic interactions between plants and the soil microorganisms (Fungus, Frankia or Rhizobium) that lead to the formation of nitrogen-fixing root nodules and/or arbuscular mycorrhiza. These interactions enable many species to survive in different marginal lands to overcome the nitrogen-and/or phosphorus deficient environment and can potentially reduce the chemical fertilizers used in agriculture which gives them an economic, social and environmental importance. The formation and the development of these structures require the mediation of specific gene products among which the transcription factors play a key role. Three of these transcription factors, viz., CYCLOPS, NSP1 and NSP2 are well conserved between actinorhizal, legume, non-legume and mycorrhizal symbioses. They interact with DELLA proteins to induce the expression of NIN in nitrogen fixing symbiosis or RAM1 in mycorrhizal symbiosis. Recently, the small non coding RNA including micro RNAs (miRNAs) have emerged as major regulators of root endosymbioses. Among them, miRNA171 targets NSP2, a TF conserved in actinorhizal, legume, non-legume and mycorrhizal symbioses. This review will also focus on the recent advances carried out on the biological function of others transcription factors during the root pre-infection/pre-contact, infection or colonization. Their role in nodule formation and AM development will also be described.

Salicylic Acid Perturbs sRNA-Gibberellin Regulatory Network in Immune Response of Potato to Potato virus Y Infection

Potato virus Y is the most economically important potato viral pathogen. We aimed at unraveling the roles of small RNAs (sRNAs) in the complex immune signaling network controlling the establishment of tolerant response of potato cv. Désirée to the virus. We constructed a sRNA regulatory network connecting sRNAs and their targets to link sRNA level responses to physiological processes. We discovered an interesting novel sRNAs-gibberellin regulatory circuit being activated as early as 3 days post inoculation (dpi) before viral multiplication can be detected. Two endogenous sRNAs, miR167 and phasiRNA931 were predicted to regulate gibberellin biosynthesis genes GA20-oxidase and GA3-oxidase. The increased expression of phasiRNA931 was also reflected in decreased levels of GA3-oxidase transcripts. Moreover, decreased concentration of gibberellin confirmed this regulation. The functional relation between lower activity of gibberellin signaling and reduced disease severity was previously confirmed in Arabidopsis-virus interaction using knockout mutants. We further showed that this regulation is salicylic acid-dependent as the response of sRNA network was attenuated in salicylic acid-depleted transgenic counterpart NahG-Désirée expressing severe disease symptoms. Besides downregulation of gibberellin signaling, regulation of immune receptor transcripts by miR6022 as well as upregulation of miR164, miR167, miR169, miR171, miR319, miR390, and miR393 in tolerant Désirée, revealed striking similarities to responses observed in mutualistic symbiotic interactions. The intertwining of different regulatory networks revealed, shows how developmental signaling, disease symptom development, and stress signaling can be balanced.

Plant response to biotic stress: Is there a common epigenetic response during plant-pathogenic and symbiotic interactions?

Plants constantly interact with pathogenic and symbiotic microorganisms. Recent studies have revealed several regulatory mechanisms controlling these interactions. Among them, the plant defense system is activated not only in response to pathogenic, but also in response to symbiotic microbes. Interestingly, shortly after symbiotic microbial recognition, the plant defense system is suppressed to promote plant infection by symbionts. Research studies have demonstrated the influence of the plant epigenome in modulating both pathogenic and symbiotic plant-microbe interactions, thereby influencing plant survival, adaptation and evolution of the plant response to microbial infections. It is however unclear if plant pathogenic and symbiotic responses share similar epigenomic profiles or if epigenomic changes differentially regulate plant-microbe symbiosis and pathogenesis. In this mini-review, we provide an update of the current knowledge of epigenomic control on plant immune responses and symbiosis, with a special attention being paid to knowledge gap and potential strategies to fill-in the missing links.

GmTIR1/GmAFB3-based auxin perception regulated by miR393 modulates soybean nodulation

Auxins play important roles in the nodulation of legumes. However, the mechanism by which auxin signaling regulates root nodulation is largely unknown. In particular, the role of auxin receptors and their regulation in determinate nodule development remains elusive. We checked the expression pattern of the auxin receptor GmTIR1/GmAFB3 genes in soybean. We analyzed the functions of GmTIR1/AFB3 in the regulation of rhizobial infection and nodule number, and also tested the functions of miR393 during nodulation and its relationship with GmTIR1/AFB3. The results showed that GmTIR1 and GmAFB3 genes exhibit diverse expression patterns during nodulation and overexpression of GmTIR1 genes significantly increased inflection foci and eventual nodule number. GmTIR1/AFB3 genes were post-transcriptionally cleaved by miR393 family and knock-down of the miR393 family members significantly increased rhizobial infection and the nodule number. Overexpression of the mutated form of GmTIR1C at the miR393 cleavage site that is resistant to miR393 cleavage led to a further increase in the number of infection foci and nodules, suggesting that miR393s modulate nodulation by directly targeting GmTIR1C. This study demonstrated that GmTIR1- and GmAFB3-mediated auxin signaling, that is spatio-temporally regulated by miR393, plays a crucial role in determinate nodule development in soybean.

Improvement of Soybean Products Through the Response Mechanism Analysis Using Proteomic Technique

Soybean is rich in protein/vegetable oil and contains several phytochemicals such as isoflavones and phenolic compounds. Because of the predominated nutritional values, soybean is considered as traditional health benefit food. Soybean is a widely cultivated crop; however, its growth and yield are markedly affected by adverse environmental conditions. Proteomic techniques make it feasible to map protein profiles both during soybean growth and under unfavorable conditions. The stress-responsive mechanisms during soybean growth have been uncovered with the help of proteomic studies. In this review, the history of soybean as food and the morphology/physiology of soybean are described. The utilization of proteomics during soybean germination and development is summarized. In addition, the stress-responsive mechanisms explored using proteomic techniques are reviewed in soybean.

Plant microRNAs: key regulators of root architecture and biotic interactions

Contents 22 I. 22 II. 24 III. 25 IV. 27 V. 29 VI. 10 31 References 32 SUMMARY: Plants have evolved a remarkable faculty of adaptation to deal with various and changing environmental conditions. In this context, the roots have taken over nutritional aspects and the root system architecture can be modulated in response to nutrient availability or biotic interactions with soil microorganisms. This adaptability requires a fine tuning of gene expression. Indeed, root specification and development are highly complex processes requiring gene regulatory networks involved in hormonal regulations and cell identity. Among the different molecular partners governing root development, microRNAs (miRNAs) are key players for the fast regulation of gene expression. miRNAs are small RNAs involved in most developmental processes and are required for the normal growth of organisms, by the negative regulation of key genes, such as transcription factors and hormone receptors. Here, we review the known roles of miRNAs in root specification and development, from the embryonic roots to the establishment of root symbioses, highlighting the major roles of miRNAs in these processes.

Early signaling, synthesis, transport and metabolism of ureides

The symbiosis between α nitrogen (N2)-fixing Proteobacteria (family Rhizobiaceae) and legumes belonging to the Fabaceae (a single phylogenetic group comprising three subfamilies: Caesalpinioideae, Mimosoideae and Papilionoideae) results in the formation of a novel root structure called a nodule, where atmospheric N2 is fixed into NH3(+). In the determinate type of nodules harbored by Rhizobium-nodulated Fabaceae species, newly synthesized NH3(+) is finally converted into allantoin (C4H6N4O3) and allantoic acid (C4H8N4O4) (ureides) through complex pathways involving at least 20 different enzymes that act synchronously in two types of nodule cells with contrasting ultrastructure, including the tree nodule cell organelles. Newly synthesized ureides are loaded into the network of nodule-root xylem vessels and transported to aerial organs by the transpirational water current. Once inside the leaves, ureides undergo an enzymatically driven reverse process to yield NH4(+) that is used for growth. This supports the role of ureides as key nitrogen (N)-compounds for the growth and yield of legumes nodulated by Rhizobium that grow in soils with a low N content. Thus, a concrete understanding of the mechanisms underlying ureide biogenesis and catabolism in legumes may help agrobiologists to achieve greater agricultural discoveries. In this review we focus on the transmembranal and transorganellar symplastic and apoplastic movement of N-precursors within the nodules, as well as on the occurrence, localization and properties of enzymes and genes involved in the biogenesis and catabolism of ureides. The synthesis and transport of ureides are not unique events in Rhizobium-nodulated N2-fixing legumes. Thus, a brief description of the synthesis and catabolism of ureides in non-legumes was included for comparison. The establishment of the symbiosis, nodule organogenesis and the plant's control of nodule number, synthesis and translocation of ureides via feed-back inhibition mechanisms are also reviewed.

Noncoding RNAs, Emerging Regulators in Root Endosymbioses

Endosymbiosis interactions allow plants to grow in nutrient-deficient soil environments. The arbuscular mycorrhizal (AM) symbiosis is an ancestral interaction between land plants and fungi, whereas nitrogen-fixing symbioses are highly specific for certain plants, notably major crop legumes. The signaling pathways triggered by specific lipochitooligosaccharide molecules involved in these interactions have common components that also overlap with plant root development. These pathways include receptor-like kinases, transcription factors (TFs), and various intermediate signaling effectors, including noncoding (nc)RNAs. These latter molecules have emerged as major regulators of gene expression and small ncRNAs, composed of micro (mi)RNAs and small interfering (si)RNAs, are known to control gene expression at transcriptional (chromatin) or posttranscriptional levels. In this review, we describe exciting recent data connecting variants of conserved si/miRNAs with the regulation of TFs, such as NSP2, NFY-A1, auxin-response factors, and AP2-like proteins, known to be involved in symbiosis. The link between hormonal regulations and these si- and miRNA-TF nodes is proposed in a model in which different feedback loops or regulations controlling endosymbiosis signaling are integrated. The diversity and emerging regulatory networks of young legume miRNAs are also highlighted.

System approaches to study root hairs as a single cell plant model: current status and future perspectives

Our current understanding of plant functional genomics derives primarily from measurements of gene, protein and/or metabolite levels averaged over the whole plant or multicellular tissues. These approaches risk diluting the response of specific cells that might respond strongly to the treatment but whose signal is diluted by the larger proportion of non-responding cells. For example, if a gene is expressed at a low level, does this mean that it is indeed lowly expressed or is it highly expressed, but only in a few cells? In order to avoid these issues, we adopted the soybean root hair cell, derived from a single, differentiated root epidermal cell, as a single-cell model for functional genomics. Root hair cells are intrinsically interesting since they are major conduits for root water and nutrient uptake and are also the preferred site of infection by nitrogen-fixing rhizobium bacteria. Although a variety of other approaches have been used to study single plant cells or single cell types, the root hair system is perhaps unique in allowing application of the full repertoire of functional genomic and biochemical approaches. In this mini review, we summarize our published work and place this within the broader context of root biology, with a significant focus on understanding the initial events in the soybean-rhizobium interaction.