Expression profiling in Medicago truncatula identifies more than 750 genes differentially expressed during nodulation, including many potential regulators of the symbiotic program

Genome and Transcriptome Analysis of the Torreya grandis WRKY Gene Family during Seed Development

Torreya grandis, an economically significant evergreen tree species exclusive to subtropical China, is highly valued for its seeds. However, the seed development process of T. grandis remains relatively unexplored. Given the pivotal role WRKY transcription factors (TFs) play in coordinating diverse cellular and biological activities, as well as crucial signaling pathways essential for plant growth and development, and the lack of comprehensive investigation into their specific functions in T. grandis, our study investigated its genome and successfully isolated 78 WRKY genes and categorized them into three distinct clades. A conserved motif analysis unveiled the presence of the characteristic WRKY domain in each identified TgWRKY protein. The examination of gene structures revealed variable numbers of introns (ranging from zero to eight) and exons (ranging from one to nine) among TgWRKY genes. A chromosomal distribution analysis demonstrated the presence of TgWRKY across eight chromosomes in T. grandis. Tissue-specific expression profiling unveiled distinctive patterns of these 78 TgWRKY genes across various tissues. Remarkably, a co-expression analysis integrating RNA-seq data and morphological assessments pinpointed the pronounced expression of TgWRKY25 during the developmental stages of T. grandis seeds. Moreover, a KEGG enrichment analysis, focusing on genes correlated with TgWRKY25 expression, suggested its potential involvement in processes such as protein processing in the endoplasmic reticulum, starch, and sucrose metabolism, thereby modulating seed development in T. grandis. These findings not only underscore the pivotal role of WRKY genes in T. grandis seed development but also pave the way for innovative breeding strategies.

The Rpf107 gene, a homolog of LOR, is required for the symbiotic nodulation of Robinia pseudoacacia

MAIN CONCLUSION

Rpf107 is involved in the infection process of rhizobia and the maintenance of symbiotic nitrogen fixation in black locust root nodules. The LURP-one related (LOR) protein family plays a pivotal role in mediating plant defense responses against both biotic and abiotic stresses. However, our understanding of its function in the symbiotic interaction between legumes and rhizobia remains limited. Here, Rpf107, a homolog of LOR, was identified in Robinia pseudoacacia (black locust). The subcellular localization of Rpf107 was analyzed, and its function was investigated using RNA interference (RNAi) and overexpression techniques. The subcellular localization assay revealed that Rpf107 was mainly distributed in the plasma membrane and nucleus. Rpf107 silencing prevented rhizobial infection and hampered plant growth. The number of infected cells in the nitrogen fixation zone of the Rpf107-RNAi nodules was also noticeably lower than that in the control nodules. Notably, Rpf107 silencing resulted in bacteroid degradation and the premature aging of nodules. In contrast, the overexpression of Rpf107 delayed the senescence of nodules and prolonged the nitrogen-fixing ability of nodules. These results demonstrate that Rpf107 was involved in the infection of rhizobia and the maintenance of symbiotic nitrogen fixation in black locust root nodules. The findings reveal that a member of the LOR protein family plays a role in leguminous root nodule symbiosis, which is helpful to clarify the functions of plant LOR protein family and fully understand the molecular mechanisms underlying legume-rhizobium symbiosis.

Genomic and Metabolic Characterization of Plant Growth-Promoting Rhizobacteria Isolated from Nodules of Clovers Grown in Non-Farmed Soil

The rhizosphere microbiota, which includes plant growth-promoting rhizobacteria (PGPR), is essential for nutrient acquisition, protection against pathogens, and abiotic stress tolerance in plants. However, agricultural practices affect the composition and functions of microbiota, reducing their beneficial effects on plant growth and health. Among PGPR, rhizobia form mutually beneficial symbiosis with legumes. In this study, we characterized 16 clover nodule isolates from non-farmed soil to explore their plant growth-promoting (PGP) potential, hypothesizing that these bacteria may possess unique, unaltered PGP traits, compared to those affected by common agricultural practices. Biolog profiling revealed their versatile metabolic capabilities, enabling them to utilize a wide range of carbon and energy sources. All isolates were effective phosphate solubilizers, and individual strains exhibited 1-aminocyclopropane-1-carboxylate deaminase and metal ion chelation activities. Metabolically active strains showed improved performance in symbiotic interactions with plants. Comparative genomics revealed that the genomes of five nodule isolates contained a significantly enriched fraction of unique genes associated with quorum sensing and aromatic compound degradation. As the potential of PGPR in agriculture grows, we emphasize the importance of the molecular and metabolic characterization of PGP traits as a fundamental step towards their subsequent application in the field as an alternative to chemical fertilizers and supplements.

RNA-Seq analysis of mung bean (Vigna radiata L.) roots shows differential gene expression and predicts regulatory pathways responding to taxonomically different rhizobia

Symbiotic interaction among legume and rhizobia is a complex phenomenon which results in the formation of nitrogen-fixing nodules. Mung bean is promiscuous host however expression profile of this important legume plant in response to rhizobial infection was particularly lacking and urgently needed. We have demonstrated the pattern of gene expression of mung bean roots inoculated with two symbionts Bradyrhizobium yuanmingense Vr50 and Sinorhizobium (Ensifer) aridi Vr33 and non-inoculated control (CK). The RNA-Seq data analyzed at two growth stages i.e., 1-3 h and 10-16 days post inoculation revealed significantly higher number of differentially expressed genes (DEGs) at nodulation stage. The DEGs encoding receptor kinases identified at early stage might be involved in perception of Nod factors produced by different rhizobia. At nodulation stage important genes involved in plant hormone signal transduction, nitrogen and sulfur metabolism were identified. KEGG pathway enrichment analysis showed that metabolic pathways were most prominent in both groups (Group 1: Vr33 vs CK; Group 2: Vr50 vs CK), followed by biosynthesis of secondary metabolites, plant hormone signal transduction and biosynthesis of amino acids. Furthermore, DEGs involved in cell communication and plant hormone signal transduction were found to be different among two symbiotic systems while DEGs involved in carbon, nitrogen and sulfur metabolism were similar but their expression varied in response to two rhizobial strains. This study provides the first insight into the mechanisms underlying interactions of mung bean host with two taxonomically different symbionts (Bradyrhizobium and Sinorhizobium) and the candidate genes for better understanding the mechanisms of symbiotic host-specificity.

RNAseq-Based Working Model for Transcriptional Regulation of Crosstalk between Simultaneous Abiotic UV-B and Biotic Stresses in Plants

Plants adjust their secondary metabolism by altering the expression of corresponding genes to cope with both abiotic and biotic stresses. In the case of UV-B radiation, plants produce protective flavonoids; however, this reaction is impeded during pattern-triggered immunity (PTI) induced by pathogens. Pathogen attack can be mimicked by the application of microbial associated molecular patterns (e.g., flg22) to study crosstalk between PTI and UV-B-induced signaling pathways. Switching from Arabidopsis cell cultures to in planta studies, we analyzed whole transcriptome changes to gain a deeper insight into crosstalk regulation. We performed a comparative transcriptomic analysis by RNAseq with four distinct mRNA libraries and identified 10778, 13620, and 11294 genes, which were differentially expressed after flg22, UV-B, and stress co-treatment, respectively. Focusing on genes being either co-regulated with the UV-B inducible marker gene chalcone synthase CHS or the flg22 inducible marker gene FRK1 identified a large set of transcription factors from diverse families, such as MYB, WRKY, or NAC. These data provide a global view of transcriptomic reprogramming during this crosstalk and constitute a valuable dataset for further deciphering the underlying regulatory mechanism(s), which appear to be much more complex than previously anticipated. The possible involvement of MBW complexes in this context is discussed.

A Nonspecific Lipid Transfer Protein with Potential Functions in Infection and Nodulation

The response of Alnus glutinosa to Frankia alni ACN14a is driven by several sequential physiological events from calcium spiking and root-hair deformation to the development of the nodule. Early stages of actinorhizal symbiosis were monitored at the transcriptional level to observe plant host responses to Frankia alni. Forty-two genes were significantly upregulated in inoculated compared with noninoculated roots. Most of these genes encode proteins involved in biological processes induced during microbial infection, such as oxidative stress or response to stimuli, but a large number of them are not differentially modulated or downregulated later in the process of nodulation. In contrast, several of them remained upregulated in mature nodules, and this included the gene most upregulated, which encodes a nonspecific lipid transfer protein (nsLTP). Classified as an antimicrobial peptide, this nsLTP was immunolocalized on the deformed root-hair surfaces that are points of contact for Frankia spp. during infection. Later in nodules, it binds to the surface of F. alni ACN14a vesicles, which are the specialized cells for nitrogen fixation. This nsLTP, named AgLTP24, was biologically produced in a heterologous host and purified for assay on F. alni ACN14a to identify physiological effects. Thus, the activation of the plant immunity response occurs upon first contact, while the recognition of F. alni ACN14a genes switches off part of the defense system during nodulation. AgLTP24 constitutes a part of the defense system that is maintained all along the symbiosis, with potential functions such as the formation of infection threads or nodule primordia to the control of F. alni proliferation. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Genes Associated with Biological Nitrogen Fixation Efficiency Identified Using RNA Sequencing in Red Clover (Trifolium pratense L.)

Commonly studied in the context of legume-rhizobia symbiosis, biological nitrogen fixation (BNF) is a key component of the nitrogen cycle in nature. Despite its potential in plant breeding and many years of research, information is still lacking as to the regulation of hundreds of genes connected with plant-bacteria interaction, nodulation, and nitrogen fixation. Here, we compared root nodule transcriptomes of red clover (Trifolium pratense L.) genotypes with contrasting nitrogen fixation efficiency, and we found 491 differentially expressed genes (DEGs) between plants with high and low BNF efficiency. The annotation of genes expressed in nodules revealed more than 800 genes not yet experimentally confirmed. Among genes mediating nodule development, four nod-ule-specific cysteine-rich (NCR) peptides were confirmed in the nodule transcriptome. Gene duplication analyses revealed that genes originating from tandem and dispersed duplication are significantly over-represented among DEGs. Weighted correlation network analysis (WGCNA) organized expression profiles of the transcripts into 16 modules linked to the analyzed traits, such as nitrogen fixation efficiency or sample-specific modules. Overall, the results obtained broaden our knowledge about transcriptomic landscapes of red clover's root nodules and shift the phenotypic description of BNF efficiency on the level of gene expression in situ.

Proteomic Analysis of Alfalfa (Medicago sativa L.) Roots in Response to Rhizobium Nodulation and Salt Stress

(1) Background: Alfalfa is an important legume forage throughout the world. Although alfalfa is considered moderately tolerant to salinity, its production and nitrogen-fixing activity are greatly limited by salt stress. (2) Methods: We examined the physiological changes and proteomic profiles of alfalfa with active nodules (NA) and without nodules (NN) under NaCl treatment. (3) Results: Our data suggested that NA roots showed upregulation of the pathways of abiotic and biotic stress responses (e.g., heat shock proteins and pathogenesis-related proteins), antioxidant enzyme synthesis, protein synthesis and degradation, cell wall degradation and modification, acid phosphatases, and porin transport when compared with NN plants under salt stress conditions. NA roots also upregulated the processes or proteins of lipid metabolism, heat shock proteins, protein degradation and folding, and cell cytoskeleton, downregulated the DNA and protein synthesis process, and vacuolar H⁺-ATPase proteins under salt stress. Besides, NA roots displayed a net H⁺ influx and low level of K⁺ efflux under salt stress, which may enhance the salt tolerance of NA plants. (4) Conclusions: The rhizobium symbiosis conferred the host plant salt tolerance by regulating a series of physiological processes to enhance stress response, improve antioxidant ability and energy use efficiency, and maintain ion homeostasis.

Time Series Transcriptome Analysis in Medicago truncatula Shoot and Root Tissue During Early Nodulation

In response to colonization by rhizobia bacteria, legumes are able to form nitrogen-fixing nodules in their roots, allowing the plants to grow efficiently in nitrogen-depleted environments. Legumes utilize a complex, long-distance signaling pathway to regulate nodulation that involves signals in both roots and shoots. We measured the transcriptional response to treatment with rhizobia in both the shoots and roots of Medicago truncatula over a 72-h time course. To detect temporal shifts in gene expression, we developed GeneShift, a novel computational statistics and machine learning workflow that addresses the time series replicate the averaging issue for detecting gene expression pattern shifts under different conditions. We identified both known and novel genes that are regulated dynamically in both tissues during early nodulation including leginsulin, defensins, root transporters, nodulin-related, and circadian clock genes. We validated over 70% of the expression patterns that GeneShift discovered using an independent M. truncatula RNA-Seq study. GeneShift facilitated the discovery of condition-specific temporally differentially expressed genes in the symbiotic nodulation biological system. In principle, GeneShift should work for time-series gene expression profiling studies from other systems.

Environmental Cues Contribute to Dynamic Plasma Membrane Organization of Nanodomains Containing Flotillin-1 and Hypersensitive Induced Reaction-1 Proteins in Arabidopsis thaliana

Plasma membranes are heterogeneous and contain multiple functional nanodomains. Although several signaling proteins have been shown to function by moving into or out of nanodomains, little is known regarding the effects of environmental cues on nanodomain organization. In this study, we investigated the heterogeneity and organization of distinct nanodomains, including those containing Arabidopsis thaliana flotillin-1 (AtFlot1) and hypersensitive induced reaction-1 proteins (AtHIR1), in response to biotic and abiotic stress. Variable-angle total internal reflection fluorescence microscopy coupled with single-particle tracking (SPT) revealed that AtFlot1 and AtHIR1 exhibit different lateral dynamics and inhabit different types of nanodomains. Furthermore, via SPT and fluorescence correlation spectroscopy, we observed lower density and intensity of AtFlot1 fluorescence in the plasma membrane after biotic stress. In contrast, the density and intensity of signal indicating AtHIR1 markedly increased in response to biotic stress. In response to abiotic stress, the density and intensity of both AtFlot1 and AtHIR1 signals decreased significantly. Importantly, SPT coupled with fluorescence recovery after photobleaching revealed that biotic and abiotic stress can regulate the dynamics of AtFlot1; however, only the abiotic stress can regulate AtHIR1 dynamics. Taken together, these findings suggest that a plethora of highly distinct nanodomains coexist in the plasma membrane (PM) and that different nanodomains may perform distinct functions in response to biotic and abiotic stresses. These phenomena may be explained by the spatial clustering of plasma membrane proteins with their associated signaling components within dedicated PM nanodomains.

Distinct genetic basis for root responses to lipo-chitooligosaccharide signal molecules from different microbial origins

Lipo-chitooligosaccharides (LCOs) were originally found as symbiotic signals called Nod Factors (Nod-LCOs) controlling the nodulation of legumes by rhizobia. More recently, LCOs were also found in symbiotic fungi and, more surprisingly, very widely in the kingdom Fungi, including in saprophytic and pathogenic fungi. The LCO-V(C18:1, fucosylated/methyl fucosylated), hereafter called Fung-LCOs, are the LCO structures most commonly found in fungi. This raises the question of how legume plants such as Medicago truncatula can discriminate between Nod-LCOs and Fung-LCOs. To address this question, we performed a genome-wide association study on 173 natural accessions of M. truncatula, using a root branching phenotype and a newly developed local score approach. Both Nod-LCOs and Fung-LCOs stimulated root branching in most accessions, but the root responses to these two types of LCO molecules were not correlated. In addition, the heritability of the root response was higher for Nod-LCOs than for Fung-LCOs. We identified 123 loci for Nod-LCO and 71 for Fung-LCO responses, of which only one was common. This suggests that Nod-LCOs and Fung-LCOs both control root branching but use different molecular mechanisms. The tighter genetic constraint of the root response to Fung-LCOs possibly reflects the ancestral origin of the biological activity of these molecules.

Interaction of Symbiotic Rhizobia and Parasitic Root-Knot Nematodes in Legume Roots: From Molecular Regulation to Field Application

Legumes form two types of root organs in response to signals from microbes, namely, nodules and root galls. In the field, these interactions occur concurrently and often interact with each other. The outcomes of these interactions vary and can depend on natural variation in rhizobia and nematode populations in the soil as well as abiotic conditions. While rhizobia are symbionts that contribute fixed nitrogen to their hosts, parasitic root-knot nematodes (RKN) cause galls as feeding structures that consume plant resources without a contribution to the plant. Yet, the two interactions share similarities, including rhizosphere signaling, repression of host defense responses, activation of host cell division, and differentiation, nutrient exchange, and alteration of root architecture. Rhizobia activate changes in defense and development through Nod factor signaling, with additional functions of effector proteins and exopolysaccharides. RKN inject large numbers of protein effectors into plant cells that directly suppress immune signaling and manipulate developmental pathways. This review examines the molecular control of legume interactions with rhizobia and RKN to elucidate shared and distinct mechanisms of these root-microbe interactions. Many of the molecular pathways targeted by both organisms overlap, yet recent discoveries have singled out differences in the spatial control of expression of developmental regulators that may have enabled activation of cortical cell division during nodulation in legumes. The interaction of legumes with symbionts and parasites highlights the importance of a comprehensive view of root-microbe interactions for future crop management and breeding strategies.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Fulvic acid increases forage legume growth inducing preferential up-regulation of nodulation and signalling-related genes

The use of potential biostimulants is of broad interest in plant science for improving yields. The application of a humic derivative called fulvic acid (FA) may improve forage crop production. FA is an uncharacterized mixture of chemicals and, although it has been reported to increase growth parameters in many species including legumes, its mode of action remains unclear. Previous studies of the action of FA have lacked appropriate controls, and few have included field trials. Here we report yield increases due to FA application in three European Medicago sativa cultivars, in studies which include the appropriate nutritional controls which hitherto have not been used. No significant growth stimulation was seen after FA treatment in grass species in this study at the treatment rate tested. Direct application to bacteria increased Rhizobium growth and, in M. sativa trials, root nodulation was stimulated. RNA transcriptional analysis of FA-treated plants revealed up-regulation of many important early nodulation signalling genes after only 3 d. Experiments in plate, glasshouse, and field environments showed yield increases, providing substantial evidence for the use of FA to benefit M. sativa forage production.

Responses of mature symbiotic nodules to the whole-plant systemic nitrogen signaling

In symbiotic root nodules of legumes, terminally differentiated rhizobia fix atmospheric N2 producing an NH4+ influx that is assimilated by the plant. The plant, in return, provides photosynthates that fuel the symbiotic nitrogen acquisition. Mechanisms responsible for the adjustment of the symbiotic capacity to the plant N demand remain poorly understood. We have investigated the role of systemic signaling of whole-plant N demand on the mature N2-fixing nodules of the model symbiotic association Medicago truncatula/Sinorhizobium using split-root systems. The whole-plant N-satiety signaling rapidly triggers reductions of both N2 fixation and allocation of sugars to the nodule. These responses are associated with the induction of nodule senescence and the activation of plant defenses against microbes, as well as variations in sugars transport and nodule metabolism. The whole-plant N-deficit responses mirror these changes: a rapid increase of sucrose allocation in response to N-deficit is associated with a stimulation of nodule functioning and development resulting in nodule expansion in the long term. Physiological, transcriptomic, and metabolomic data together provide evidence for strong integration of symbiotic nodules into whole-plant nitrogen demand by systemic signaling and suggest roles for sugar allocation and hormones in the signaling mechanisms.

Plant-Microbe Symbiosis: What Has Proteomics Taught Us?

Beneficial microbes have a positive impact on the productivity and fitness of the host plant. A better understanding of the biological impacts and underlying mechanisms by which the host derives these benefits will help to address concerns around global food production and security. The recent development of omics-based technologies has broadened our understanding of the molecular aspects of beneficial plant-microbe symbiosis. Specifically, proteomics has led to the identification and characterization of several novel symbiosis-specific and symbiosis-related proteins and post-translational modifications that play a critical role in mediating symbiotic plant-microbe interactions and have helped assess the underlying molecular aspects of the symbiotic relationship. Integration of proteomic data with other "omics" data can provide valuable information to assess hypotheses regarding the underlying mechanism of symbiosis and help define the factors affecting the outcome of symbiosis. Herein, an update is provided on the current and potential applications of symbiosis-based "omic" approaches to dissect different aspects of symbiotic plant interactions. The application of proteomics, metaproteomics, and secretomics as enabling approaches for the functional analysis of plant-associated microbial communities is also discussed.

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.

Reprogramming of Root Cells during Nitrogen-Fixing Symbiosis Involves Dynamic Polysome Association of Coding and Noncoding RNAs

Translational control is a widespread mechanism that allows the cell to rapidly modulate gene expression in order to provide flexibility and adaptability to eukaryotic organisms. We applied translating ribosome affinity purification combined with RNA sequencing to characterize translational regulation of mRNAs at early stages of the nitrogen-fixing symbiosis established between Medicago truncatula and Sinorhizobium meliloti Our analysis revealed a poor correlation between transcriptional and translational changes and identified hundreds of regulated protein-coding and long noncoding RNAs (lncRNAs), some of which are regulated in specific cell types. We demonstrated that a short variant of the lncRNA Trans-acting small interference RNA3 (TAS3) increased its association to the translational machinery in response to rhizobia. Functional analysis revealed that this short variant of TAS3 might act as a target mimic that captures microRNA390, contributing to reduce trans acting small interference Auxin Response Factor production and modulating nodule formation and rhizobial infection. The analysis of alternative transcript variants identified a translationally upregulated mRNA encoding subunit 3 of the SUPERKILLER complex (SKI3), which participates in mRNA decay. Knockdown of SKI3 decreased nodule initiation and development, as well as the survival of bacteria within nodules. Our results highlight the importance of translational control and mRNA decay pathways for the successful establishment of the nitrogen-fixing symbiosis.

Gene Expression in Nitrogen-Fixing Symbiotic Nodule Cells in Medicago truncatula and Other Nodulating Plants

Root nodules formed by plants of the nitrogen-fixing clade (NFC) are symbiotic organs that function in the maintenance and metabolic integration of large populations of nitrogen-fixing bacteria. These organs feature unique characteristics and processes, including their tissue organization, the presence of specific infection structures called infection threads, endocytotic uptake of bacteria, symbiotic cells carrying thousands of intracellular bacteria without signs of immune responses, and the integration of symbiont and host metabolism. The early stages of nodulation are governed by a few well-defined functions, which together constitute the common symbiosis-signaling pathway (CSSP). The CSSP activates a set of transcription factors (TFs) that orchestrate nodule organogenesis and infection. The later stages of nodule development require the activation of hundreds to thousands of genes, mostly expressed in symbiotic cells. Many of these genes are only active in symbiotic cells, reflecting the unique nature of nodules as plant structures. Although how the nodule-specific transcriptome is activated and connected to early CSSP-signaling is poorly understood, candidate TFs have been identified using transcriptomic approaches, and the importance of epigenetic and chromatin-based regulation has been demonstrated. We discuss how gene regulation analyses have advanced our understanding of nodule organogenesis, the functioning of symbiotic cells, and the evolution of symbiosis in the NFC.

Organization and dynamics of functional plant membrane microdomains

Plasma membranes are heterogeneous and laterally compartmentalized into distinct microdomains. These membrane microdomains consist of special lipids and proteins and are thought to act as signaling platforms. In plants, membrane microdomains have been detected by super-resolution microscopy, and there is evidence that they play roles in several biological processes. Here, we review current knowledge about the lipid and protein components of membrane microdomains. Furthermore, we summarize the dynamics of membrane microdomains in response to different stimuli. We also explore the biological functions associated with membrane microdomains as signal integration hubs. Finally, we outline challenges and questions for further studies.

Identifying Temporally Regulated Root Nodulation Biomarkers Using Time Series Gene Co-Expression Network Analysis

Root nodulation results from a symbiotic relationship between a plant host and Rhizobium bacteria. Synchronized gene expression patterns over the course of rhizobial infection result in activation of pathways that are unique but overlapping with the highly conserved pathways that enable mycorrhizal symbiosis. We performed RNA sequencing of 30 Medicago truncatula root maturation zone samples at five distinct time points. These samples included plants inoculated with Sinorhizobium medicae and control plants that did not receive any Rhizobium. Following gene expression quantification, we identified 1,758 differentially expressed genes at various time points. We constructed a gene co-expression network (GCN) from the same data and identified link community modules (LCMs) that were comprised entirely of differentially expressed genes at specific time points post-inoculation. One LCM included genes that were up-regulated at 24 h following inoculation, suggesting an activation of allergen family genes and carbohydrate-binding gene products in response to Rhizobium. We also identified two LCMs that were comprised entirely of genes that were down regulated at 24 and 48 h post-inoculation. The identity of the genes in these modules suggest that down-regulating specific genes at 24 h may result in decreased jasmonic acid production with an increase in cytokinin production. At 48 h, coordinated down-regulation of a specific set of genes involved in lipid biosynthesis may play a role in nodulation. We show that GCN-LCM analysis is an effective method to preliminarily identify polygenic candidate biomarkers of root nodulation and develop hypotheses for future discovery.

The Ectomycorrhizal Fungus Laccaria bicolor Produces Lipochitooligosaccharides and Uses the Common Symbiosis Pathway to Colonize Populus Roots

Mycorrhizal fungi form mutualistic associations with the roots of most land plants and provide them with mineral nutrients from the soil in exchange for fixed carbon derived from photosynthesis. The common symbiosis pathway (CSP) is a conserved molecular signaling pathway in all plants capable of associating with arbuscular mycorrhizal fungi. It is required not only for arbuscular mycorrhizal symbiosis but also for rhizobia-legume and actinorhizal symbioses. Given its role in such diverse symbiotic associations, we hypothesized that the CSP also plays a role in ectomycorrhizal associations. We showed that the ectomycorrhizal fungus Laccaria bicolor produces an array of lipochitooligosaccharides (LCOs) that can trigger both root hair branching in legumes and, most importantly, calcium spiking in the host plant Populus in a CASTOR/POLLUX-dependent manner. Nonsulfated LCOs enhanced lateral root development in Populus in a calcium/calmodulin-dependent protein kinase (CCaMK)-dependent manner, and sulfated LCOs enhanced the colonization of Populus by L. bicolor Compared with the wild-type Populus, the colonization of CASTOR/POLLUX and CCaMK RNA interference lines by L. bicolor was reduced. Our work demonstrates that similar to other root symbioses, L. bicolor uses the CSP for the full establishment of its mutualistic association with Populus.

Nodulation Induces Systemic Resistance of Medicago truncatula and Pisum sativum Against Erysiphe pisi and Primes for Powdery Mildew-Triggered Salicylic Acid Accumulation

Plants encounter beneficial and detrimental microorganisms both above- and belowground and the health status of the plant depends on the composition of this pan-microbiome. Beneficial microorganisms contribute to plant nutrition or systemically or locally protect plants against pathogens, thus facilitating adaptation to a variety of environments. Induced systemic resistance, caused by root-associated microbes, manifests as aboveground resistance against necrotrophic pathogens and is mediated by jasmonic acid/ethylene-dependent signaling. By contrast, systemic acquired resistance relies on salicylic acid (SA) signaling and confers resistance against secondary infection by (hemi)biotrophic pathogens. To investigate whether symbiotic rhizobia that are ubiquitously found in natural ecosystems are able to modulate resistance against biotrophs, we tested the impact of preestablished nodulation of Medicago truncatula and pea (Pisum sativum) plants against infection by the powdery mildew fungus Erysiphe pisi. We found that root symbiosis interfered with fungal penetration of M. truncatula and reduced asexual spore formation on pea leaves independently of symbiotic nitrogen fixation. Improved resistance of nodulated plants correlated with elevated levels of free SA and SA-dependent marker gene expression upon powdery mildew infection. Our results suggest that nodulation primes the plants systemically for E. pisi-triggered SA accumulation and defense gene expression, resulting in increased resistance.

Transcriptome analysis of soybean (Glycine max) root genes differentially expressed in rhizobial, arbuscular mycorrhizal, and dual symbiosis

Soybean (Glycine max) roots establish associations with nodule-inducing rhizobia and arbuscular mycorrhizal (AM) fungi. Both rhizobia and AM fungi have been shown to affect the activity of and colonization by the other, and their interactions can be detected within host plants. Here, we report the transcription profiles of genes differentially expressed in soybean roots in the presence of rhizobial, AM, or rhizobial-AM dual symbiosis, compared with those in control (uninoculated) roots. Following inoculation, soybean plants were grown in a glasshouse for 6 weeks; thereafter their root transcriptomes were analyzed using an oligo DNA microarray. Among the four treatments, the root nodule number and host plant growth were highest in plants with dual symbiosis. We observed that the expression of 187, 441, and 548 host genes was up-regulated and 119, 1,439, and 1,298 host genes were down-regulated during rhizobial, AM, and dual symbiosis, respectively. The expression of 34 host genes was up-regulated in each of the three symbioses. These 34 genes encoded several membrane transporters, type 1 metallothionein, and transcription factors in the MYB and bHLH families. We identified 56 host genes that were specifically up-regulated during dual symbiosis. These genes encoded several nodulin proteins, phenylpropanoid metabolism-related proteins, and carbonic anhydrase. The nodulin genes up-regulated by the AM fungal colonization probably led to the observed increases in root nodule number and host plant growth. Some other nodulin genes were down-regulated specifically during AM symbiosis. Based on the results above, we suggest that the contribution of AM fungal colonization is crucial to biological N2-fixation and host growth in soybean with rhizobial-AM dual symbiosis.

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

MAIN CONCLUSION

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

Identification and characterization of aquaporin genes in Arachis duranensis and Arachis ipaensis genomes, the diploid progenitors of peanut

BACKGROUND

Aquaporins (AQPs) facilitate transport of water and small solutes across cell membranes and play an important role in different physiological processes in plants. Despite their importance, limited data is available about AQP distribution and function in the economically important oilseed crop peanut, Arachis hypogea (AABB). The present study reports the identification and structural and expression analysis of the AQPs found in the diploid progenitor genomes of A. hypogea i.e. Arachis duranensis (AA) and Arachis ipaensis (BB).

RESULTS

Genome-wide analysis revealed the presence of 32 and 36 AQPs in A. duranensis and A. ipaensis, respectively. Phylogenetic analysis showed similar numbers of AQPs clustered in five distinct subfamilies including the plasma membrane intrinsic proteins (PIPs), the tonoplast intrinsic proteins (TIPs), the nodulin 26-like intrinsic proteins (NIPs), the small basic intrinsic proteins (SIPs), and the uncharacterized intrinsic proteins (XIPs). A notable exception was the XIP subfamily where XIP1 group was observed only in A. ipaensis genome. Protein structure evaluation showed a hydrophilic aromatic/arginine (ar/R) selectivity filter (SF) in PIPs whereas other subfamilies mostly contained a hydrophobic ar/R SF. Both genomes contained one NIP2 with a GSGR SF indicating a conserved ability within the genus to uptake silicon. Analysis of RNA-seq data from A. hypogea revealed a similar expression pattern for the different AQP paralogs of AA and BB genomes. The TIP3s showed seed-specific expression while the NIP1s' expression was confined to roots and root nodules.

CONCLUSIONS

The identification and the phylogenetic analysis of AQPs in both Arachis species revealed the presence of all five sub-families of AQPs. Within the NIP subfamily, the presence of a NIP2 in both genomes supports a conserved ability to absorb Si within plants of the genus. The global expression profile of AQPs in A. hypogea revealed a similar pattern of AQP expression regardless of the subfamilies or the genomes. The tissue-specific expression of AQPs suggests an important role in the development and function of the respective organs. The AQPs identified in the present study will serve as a resource for further characterization and possible exploitation of AQPs to understand their physiological role in A. hypogea.

Unified Transcriptomic Signature of Arbuscular Mycorrhiza Colonization in Roots of Medicago truncatula by Integration of Machine Learning, Promoter Analysis, and Direct Merging Meta-Analysis

Plant root symbiosis with Arbuscular mycorrhizal (AM) fungi improves uptake of water and mineral nutrients, improving plant development under stressful conditions. Unraveling the unified transcriptomic signature of a successful colonization provides a better understanding of symbiosis. We developed a framework for finding the transcriptomic signature of Arbuscular mycorrhiza colonization and its regulating transcription factors in roots of Medicago truncatula. Expression profiles of roots in response to AM species were collected from four separate studies and were combined by direct merging meta-analysis. Batch effect, the major concern in expression meta-analysis, was reduced by three normalization steps: Robust Multi-array Average algorithm, Z-standardization, and quartiling normalization. Then, expression profile of 33685 genes in 18 root samples of Medicago as numerical features, as well as study ID and Arbuscular mycorrhiza type as categorical features, were mined by seven models: RELIEF, UNCERTAINTY, GINI INDEX, Chi Squared, RULE, INFO GAIN, and INFO GAIN RATIO. In total, 73 genes selected by machine learning models were up-regulated in response to AM (Z-value difference > 0.5). Feature weighting models also documented that this signature is independent from study (batch) effect. The AM inoculation signature obtained was able to differentiate efficiently between AM inoculated and non-inoculated samples. The AP2 domain class transcription factor, GRAS family transcription factors, and cyclin-dependent kinase were among the highly expressed meta-genes identified in the signature. We found high correspondence between the AM colonization signature obtained in this study and independent RNA-seq experiments on AM colonization, validating the repeatability of the colonization signature. Promoter analysis of upregulated genes in the transcriptomic signature led to the key regulators of AM colonization, including the essential transcription factors for endosymbiosis establishment and development such as NF-YA factors. The approach developed in this study offers three distinct novel features: (I) it improves direct merging meta-analysis by integrating supervised machine learning models and normalization steps to reduce study-specific batch effects; (II) seven attribute weighting models assessed the suitability of each gene for the transcriptomic signature which contributes to robustness of the signature (III) the approach is justifiable, easy to apply, and useful in practice. Our integrative framework of meta-analysis, promoter analysis, and machine learning provides a foundation to reveal the transcriptomic signature and regulatory circuits governing Arbuscular mycorrhizal symbiosis and is transferable to the other biological settings.

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.

Compatibility between Legumes and Rhizobia for the Establishment of a Successful Nitrogen-Fixing Symbiosis

The root nodule symbiosis established between legumes and rhizobia is an exquisite biological interaction responsible for fixing a significant amount of nitrogen in terrestrial ecosystems. The success of this interaction depends on the recognition of the right partner by the plant within the richest microbial ecosystems on Earth, the soil. Recent metagenomic studies of the soil biome have revealed its complexity, which includes microorganisms that affect plant fitness and growth in a beneficial, harmful, or neutral manner. In this complex scenario, understanding the molecular mechanisms by which legumes recognize and discriminate rhizobia from pathogens, but also between distinct rhizobia species and strains that differ in their symbiotic performance, is a considerable challenge. In this work, we will review how plants are able to recognize and select symbiotic partners from a vast diversity of surrounding bacteria. We will also analyze recent advances that contribute to understand changes in plant gene expression associated with the outcome of the symbiotic interaction. These aspects of nitrogen-fixing symbiosis should contribute to translate the knowledge generated in basic laboratory research into biotechnological advances to improve the efficiency of the nitrogen-fixing symbiosis in agronomic systems.

The Medicago truncatula Genome: Genomic Data Availability

Medicago truncatula emerged in 1990 as a model for legumes, comprising the third largest land plant family. Most legumes form symbiotic nitrogen-fixing root nodules with compatible soil bacteria and thus are important contributors to the global nitrogen cycle and sustainable agriculture. Legumes and legume products are important sources for human and animal protein as well as for edible and industrial oils. In the years since M. truncatula was chosen as a legume model, many genetic, genomic, and molecular resources have become available, including reference quality genome sequences for two widely used genotypes. Accessibility of genomic data is important for many different types of studies with M. truncatula as well as for research involving crop and forage legumes. In this chapter, we discuss strategies to obtain archived M. truncatula genomic data originally deposited into custom databases that are no longer maintained but are now accessible in general databases. We also review key current genomic databases that are specific to M. truncatula as well as those that contain M. truncatula data in addition to data from other plants.

LysM Receptor-Like Kinase and LysM Receptor-Like Protein Families: An Update on Phylogeny and Functional Characterization

Members of plant specific families of receptor-like kinases (RLKs) and receptor-like proteins (RLPs), containing 3 extracellular LysMs have been shown to directly bind and/or to be involved in perception of lipo-chitooligosaccharides (LCO), chitooligosaccharides (CO), and peptidoglycan (PGN), three types of GlcNAc-containing molecules produced by microorganisms. These receptors are involved in microorganism perception by plants and can activate different plant responses leading either to symbiosis establishment or to defense responses against pathogens. LysM-RLK/Ps belong to multigenic families. Here, we provide a phylogeny of these families in eight plant species, including dicotyledons and monocotyledons, and we discuss known or putative biological roles of the members in each of the identified phylogenetic groups. We also report and discuss known biochemical properties of the LysM-RLK/Ps.

Effects of biosolids from a wastewater treatment plant receiving manufactured nanomaterials on Medicago truncatula and associated soil microbial communities at low nanomaterial concentrations

Concern has grown regarding engineered nanomaterials (ENMs) entering agricultural soils through the application of biosolids and their possible effects on agroecosystems, even though the ENMs are extensively transformed. The effects of exposure to biosolids containing transformation products of these ENMs at low concentrations remain largely unexplored. We examined the responses of Medicago truncatula and its symbiotic rhizobia Sinorhizobium meliloti exposed to soil amended with biosolids from WWTP containing low added concentrations of ENMs (ENM Low), bulk/dissolved metals (bulk/dissolved Low), or no metal additions (control). We targeted adding approximately 5mg/kg of Ag and 50mg/kg of Zn, and Ti. Measured endpoints included M. truncatula growth, nodulation, changes in the expression of stress response genes, uptake of metals (Ag, Zn and Ti) into shoots, and quantification of S. meliloti populations and soil microbial communities. After 30days exposure, no effects on root or shoot biomass were observed in ENM Low and bulk/dissolved Low treatments, whereas both treatments had a larger average number of nodules (5.7 and 5.57, respectively) compared to controls (0.33). There were no significant differences in either total accumulated metal or metal concentrations in shoots among the treatments. Expression of five stress-related genes (metal tolerance protein (MTP), metal transporter (MTR), peroxidase (PEROX), NADPH oxidase (NADPH) and 1-aminocyclopropane-1-carboxylate oxidase-like protein (ACC_Oxidase)) was significantly down-regulated in both bulk/dissolved Low and ENM Low treatments. However, a change in soil microbial community composition and a significant increase in total microbial biomass were observed in ENM Low relative to control. The ENM Low treatment had increased abundance of Gram-negative and anaerobic bacteria and reduced abundance of eukaryotes compared to control. The study demonstrated that although there were some subtle shifts in microbial community composition, plant health was minimally impacted by ENMs within the time frame and at the low exposure concentrations used in this study.

Transcriptome analysis of two recombinant inbred lines of common bean contrasting for symbiotic nitrogen fixation

Common bean (Phaseolus vulgaris L.) fixes atmospheric nitrogen (N2) through symbiotic nitrogen fixation (SNF) at levels lower than other grain legume crops. An understanding of the genes and molecular mechanisms underlying SNF will enable more effective strategies for the genetic improvement of SNF traits in common bean. In this study, transcriptome profiling was used to identify genes and molecular mechanisms underlying SNF differences between two common bean recombinant inbred lines that differed in their N-fixing abilities. Differential gene expression and functional enrichment analyses were performed on leaves, nodules and roots of the two lines when grown under N-fixing and non-fixing conditions. Receptor kinases, transmembrane transporters, and transcription factors were among the differentially expressed genes identified under N-fixing conditions, but not under non-fixing conditions. Genes up-regulated in the stronger nitrogen fixer, SA36, included those involved in molecular functions such as purine nucleoside binding, oxidoreductase and transmembrane receptor activities in nodules, and transport activity in roots. Transcription factors identified in this study are candidates for future work aimed at understanding the functional role of these genes in SNF. Information generated in this study will support the development of gene-based markers to accelerate genetic improvement of SNF in common bean.

Plant Carbonic Anhydrases: Structures, Locations, Evolution, and Physiological Roles

Carbonic anhydrases (CAs) are zinc metalloenzymes that catalyze the interconversion of CO2 and HCO3- and are ubiquitous in nature. Higher plants contain three evolutionarily distinct CA families, αCAs, βCAs, and γCAs, where each family is represented by multiple isoforms in all species. Alternative splicing of CA transcripts appears common; consequently, the number of functional CA isoforms in a species may exceed the number of genes. CAs are expressed in numerous plant tissues and in different cellular locations. The most prevalent CAs are those in the chloroplast, cytosol, and mitochondria. This diversity in location is paralleled in the many physiological and biochemical roles that CAs play in plants. In this review, the number and types of CAs in C3, C4, and crassulacean acid metabolism (CAM) plants are considered, and the roles of the α and γCAs are briefly discussed. The remainder of the review focuses on plant βCAs and includes the identification of homologs between species using phylogenetic approaches, a consideration of the inter- and intracellular localization of the proteins, along with the evidence for alternative splice forms. Current understanding of βCA tissue-specific expression patterns and what controls them are reviewed, and the physiological roles for which βCAs have been implicated are presented.

Plant NF-Y transcription factors: Key players in plant-microbe interactions, root development and adaptation to stress

NF-Ys are heterotrimeric transcription factors composed by the NF-YA, NF-YB and NF-YC subunits. In plants, NF-Y subunits are encoded by multigene families whose members show structural and functional diversifications. An increasing number of NF-Y genes has been shown to play key roles during different stages of root nodule and arbuscular mycorrhizal symbiosis, as well as during the interaction of plants with pathogenic microorganisms. Individual members of the NF-YA and NF-YB families have also been implicated in the development of primary and lateral roots. In addition, different members of the NF-YA and NF-YB gene families from mono- and di-cotyledonous plants have been involved in plant responses to water and nutrient scarcity. This review presents the most relevant and striking results concerning these NF-Y subunits. A phylogenetic analysis of the functionally characterized NF-Y genes revealed that, across plant species, NF-Y proteins functioning in the same biological process tend to belong to common phylogenetic groups. Finally, we discuss the forthcoming challenges of plant NF-Y research, including the detailed dissection of expression patterns, the elucidation of functional specificities as well as the characterization of the potential NF-Y-mediated epigenetic mechanisms by which they control the expression of their target genes. This article is part of a Special Issue entitled: Nuclear Factor Y in Development and Disease, edited by Prof. Roberto Mantovani.

Root nodule symbiosis in Lotus japonicus drives the establishment of distinctive rhizosphere, root, and nodule bacterial communities

Lotus japonicus has been used for decades as a model legume to study the establishment of binary symbiotic relationships with nitrogen-fixing rhizobia that trigger root nodule organogenesis for bacterial accommodation. Using community profiling of 16S rRNA gene amplicons, we reveal that in Lotus, distinctive nodule- and root-inhabiting communities are established by parallel, rather than consecutive, selection of bacteria from the rhizosphere and root compartments. Comparative analyses of wild-type (WT) and symbiotic mutants in Nod factor receptor5 (nfr5), Nodule inception (nin) and Lotus histidine kinase1 (lhk1) genes identified a previously unsuspected role of the nodulation pathway in the establishment of different bacterial assemblages in the root and rhizosphere. We found that the loss of nitrogen-fixing symbiosis dramatically alters community structure in the latter two compartments, affecting at least 14 bacterial orders. The differential plant growth phenotypes seen between WT and the symbiotic mutants in nonsupplemented soil were retained under nitrogen-supplemented conditions that blocked the formation of functional nodules in WT, whereas the symbiosis-impaired mutants maintain an altered community structure in the nitrogen-supplemented soil. This finding provides strong evidence that the root-associated community shift in the symbiotic mutants is a direct consequence of the disabled symbiosis pathway rather than an indirect effect resulting from abolished symbiotic nitrogen fixation. Our findings imply a role of the legume host in selecting a broad taxonomic range of root-associated bacteria that, in addition to rhizobia, likely contribute to plant growth and ecological performance.

Dissecting the Root Nodule Transcriptome of Chickpea (Cicer arietinum L.)

A hallmark trait of chickpea (Cicer arietinum L.), like other legumes, is the capability to convert atmospheric nitrogen (N2) into ammonia (NH3) in symbiotic association with Mesorhizobium ciceri. However, the complexity of molecular networks associated with the dynamics of nodule development in chickpea need to be analyzed in depth. Hence, in order to gain insights into the chickpea nodule development, the transcriptomes of nodules at early, middle and late stages of development were sequenced using the Roche 454 platform. This generated 490.84 Mb sequence data comprising 1,360,251 reads which were assembled into 83,405 unigenes. Transcripts were annotated using Gene Ontology (GO), Cluster of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways analysis. Differential expression analysis revealed that a total of 3760 transcripts were differentially expressed in at least one of three stages, whereas 935, 117 and 2707 transcripts were found to be differentially expressed in the early, middle and late stages of nodule development respectively. MapMan analysis revealed enrichment of metabolic pathways such as transport, protein synthesis, signaling and carbohydrate metabolism during root nodulation. Transcription factors were predicted and analyzed for their differential expression during nodule development. Putative nodule specific transcripts were identified and enriched for GO categories using BiNGO which revealed many categories to be enriched during nodule development, including transcription regulators and transporters. Further, the assembled transcriptome was also used to mine for genic SSR markers. In conclusion, this study will help in enriching the transcriptomic resources implicated in understanding of root nodulation events in chickpea.

Genome-wide analysis identifies gain and loss/change of function within the small multigenic insecticidal Albumin 1 family of Medicago truncatula

BACKGROUND

Albumin 1b peptides (A1b) are small disulfide-knotted insecticidal peptides produced by Fabaceae (also called Leguminosae). To date, their diversity among this plant family has been essentially investigated through biochemical and PCR-based approaches. The availability of high-quality genomic resources for several fabaceae species, among which the model species Medicago truncatula (Mtr), allowed for a genomic analysis of this protein family aimed at i) deciphering the evolutionary history of A1b proteins and their links with A1b-nodulins that are short non-insecticidal disulfide-bonded peptides involved in root nodule signaling and ii) exploring the functional diversity of A1b for novel bioactive molecules.

RESULTS

Investigating the Mtr genome revealed a remarkable expansion, mainly through tandem duplications, of albumin1 (A1) genes, retaining nearly all of the same canonical structure at both gene and protein levels. Phylogenetic analysis revealed that the ancestral molecule was most probably insecticidal giving rise to, among others, A1b-nodulins. Expression meta-analysis revealed that many A1b coding genes are silent and a wide tissue distribution of the A1 transcripts/peptides within plant organs. Evolutionary rate analyses highlighted branches and sites with positive selection signatures, including two sites shown to be critical for insecticidal activity. Seven peptides were chemically synthesized and folded in vitro, then assayed for their biological activity. Among these, AG41 (aka MtrA1013 isoform, encoded by the orphan TA24778 contig.), showed an unexpectedly high insecticidal activity. The study highlights the unique burst of diversity of A1 peptides within the Medicago genus compared to the other taxa for which full-genomes are available: no A1 member in Lotus, or in red clover to date, while only a few are present in chick pea, soybean or pigeon pea genomes.

CONCLUSION

The expansion of the A1 family in the Medicago genus is reminiscent of the situation described for another disulfide-rich peptide family, the "Nodule-specific Cysteine-Rich" (NCR), discovered within the same species. The oldest insecticidal A1b toxin was described from the Sophorae, dating the birth of this seed-defense function to more than 58 million years, and making this model of plant/insect toxin/receptor (A1b/insect v-ATPase) one of the oldest known.

How legumes recognize rhizobia

Legume plants have developed the capacity to establish symbiotic interactions with soil bacteria (known as rhizobia) that can convert N2 to molecular forms that are incorporated into the plant metabolism. The first step of this relationship is the recognition of bacteria by the plant, which allows to distinguish potentially harmful species from symbiotic partners. The main molecular determinant of this symbiotic interaction is the Nod Factor, a diffusible lipochitooligosaccharide molecule produced by rhizobia and perceived by LysM receptor kinases; however, other important molecules involved in the specific recognition have emerged over the years. Secreted exopolysaccharides and the lipopolysaccharides present in the bacterial cell wall have been proposed to act as signaling molecules, triggering the expression of specific genes related to the symbiotic process. In this review we will briefly discuss how transcriptomic analysis are helping to understand how multiple signaling pathways, triggered by the perception of different molecules produced by rhizobia, control the genetic programs of root nodule organogenesis and bacterial infection. This knowledge can help to understand how legumes have evolved to recognize and establish complex ecological relationships with particular species and strains of rhizobia, adjusting gene expression in response to identity determinants of bacteria.

MtNF-YA1, A Central Transcriptional Regulator of Symbiotic Nodule Development, Is Also a Determinant of Medicago truncatula Susceptibility toward a Root Pathogen

Plant NF-Y transcription factors control a wide array of biological functions enabling appropriate reproductive and developmental processes as well as adaptation to various abiotic and biotic environments. In Medicago truncatula, MtNF-YA1 was previously identified as a key determinant for nodule development and establishment of rhizobial symbiosis. Here, we highlight a new role for this protein in compatibility to Aphanomyces euteiches, a root pathogenic oomycete. The Mtnf-ya1-1 mutant plants showed better survival rate, reduced symptoms, and increased development of their root apparatus as compared to their wild-type (WT) background A17. MtNF-YA-1 was specifically up-regulated by A. euteiches in F83005.5, a highly susceptible natural accession of M. truncatula while transcript level remained stable in A17, which is partially resistant. The role of MtNF-YA1 in F83005.5 susceptibility was further documented by reducing MtNF-YA1 expression either by overexpression of the miR169q, a microRNA targeting MtNF-YA1, or by RNAi approaches leading to a strong enhancement in the resistance of this susceptible line. Comparative analysis of the transcriptome of WT and Mtnf-ya1-1 led to the identification of 1509 differentially expressed genes. Among those, almost 36 defense-related genes were constitutively expressed in Mtnf-ya1-1, while 20 genes linked to hormonal pathways were repressed. In summary, we revealed an unexpected dual role for this symbiotic transcription factor as a key player in the compatibility mechanisms to a pathogen.

De novo transcriptome analysis of Medicago falcata reveals novel insights about the mechanisms underlying abiotic stress-responsive pathway

Background

The entire world is facing a deteriorating environment. Understanding the mechanisms underlying plant responses to external abiotic stresses is important for breeding stress-tolerant crops and herbages. Phytohormones play critical regulatory roles in plants in the response to external and internal cues to regulate growth and development. Medicago falcata is one of the stress-tolerant candidate leguminous species and is able to fix atmospheric nitrogen. This ability allows leguminous plants to grow in nitrogen deficient soils.

Methods

We performed Illumina sequencing of cDNA prepared from abiotic stress treated M. falcata. Sequencedreads were assembled to provide a transcriptome resource. Transcripts were annotated using BLASTsearches against the NCBI non-redundant database and gene ontology definitions were assigned. Acomparison among the three abiotic stress treated samples was carried out. The expression of transcriptswas confirmed with qRT-PCR.

Results

We present an abiotic stress-responsive M. falcata transcriptome using next-generation sequencing data from samples grown under standard, dehydration, high salinity, and cold conditions. We combined reads from all samples and de novo assembled 98,515 transcripts to build the M. falcata gene index. A comprehensive analysis of the transcriptome revealed abiotic stress-responsive mechanisms underlying the metabolism and core signalling components of major phytohormones. We identified nod factor signalling pathways during early symbiotic nodulation that are modified by abiotic stresses. Additionally, a global comparison of homology between the M. falcata and M. truncatula transcriptomes, along with five other leguminous species, revealed a high level of global sequence conservation within the family.

Conclusions

M. falcata is shown to be a model candidate for studying abiotic stress-responsive mechanisms in legumes. This global gene expression analysis provides new insights into the biochemical and molecular mechanisms involved in the acclimation to abiotic stresses. Our data provides many gene candidates that might be used for herbage and crop breeding. Additionally, FalcataBase (http://bioinformatics.cau.edu.cn/falcata/) was built for storing these data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2019-x) contains supplementary material, which is available to authorized users.

A JAZ Protein in Astragalus sinicus Interacts with a Leghemoglobin through the TIFY Domain and Is Involved in Nodule Development and Nitrogen Fixation

Leghemoglobins (Lbs) play an important role in legumes-rhizobia symbiosis. Lbs bind O2 and protect nitrogenase activity from damage by O2 in nodules, therefore, they are regarded as a marker of active nitrogen fixation in nodules. Additionally, Lbs are involved in the nitric oxide (NO) signaling pathway, acting as a NO scavenger during nodule development and nitrogen fixation. However, regulators responsible for Lb expression and modulation of Lb activity have not been characterized. In our previous work, a Jasmonate-Zim-domain (JAZ) protein interacting with a Lb (AsB2510) in Astragalus sinicus was identified and designated AsJAZ1. In this study, the interaction between AsJAZ1 and AsB2510 was verified using a yeast two-hybrid system and in vitro Glutathione S-transferase (GST) pull-down assays, resulting in identification of the interaction domain as a TIFY (previously known as zinc-finger protein expressed in inflorescence meristem, ZIM) domain. TIFY domain is named after the most conserved amino acids within the domain. Bimolecular fluorescence complementation (BiFC) was used to confirm the interaction between AsJAZ1 and AsB2510 in tobacco cells, demonstrating that AsJAZ1-AsB2510 interaction was localized to the cell membrane and cytoplasm. Furthermore, the expression patterns and the symbiotic phenotypes of AsJAZ1 were investigated. Knockdown of AsJAZ1 expression via RNA interference led to decreased number of nodules, abnormal development of bacteroids, accumulation of poly-x-hydroxybutyrate (PHB) and loss of nitrogenase activity. Taken together, our results suggest that AsJAZ1 interacts with AsB2510 and participates in nodule development and nitrogen fixation. Our results provide novel insights into the functions of Lbs or JAZ proteins during legume-rhizobia symbiosis.

Combined genetic and transcriptomic analysis reveals three major signalling pathways activated by Myc-LCOs in Medicago truncatula

Myc-LCOs are newly identified symbiotic signals produced by arbuscular mycorrhizal (AM) fungi. Like rhizobial Nod factors, they are lipo-chitooligosaccharides that activate the common symbiotic signalling pathway (CSSP) in plants. To increase our limited understanding of the roles of Myc-LCOs we aimed to analyse Myc-LCO-induced transcriptional changes and their genetic control. Whole genome RNA sequencing (RNA-seq) was performed on roots of Medicago truncatula wild-type plants, and dmi3 and nsp1 symbiotic mutants affected in nodulation and mycorrhizal signalling. Plants were treated separately with the two major types of Myc-LCOs, sulphated and nonsulphated. Generalized linear model analysis identified 2201 differentially expressed genes and classified them according to genotype and/or treatment effects. Three genetic pathways for Myc-LCO-regulation of transcriptomic reprogramming were highlighted: DMI3- and NSP1-dependent; DMI3-dependent and NSP1-independent; and DMI3- and NSP1-independent. Comprehensive analysis revealed overlaps with previous AM studies, and highlighted certain functions, especially signalling components and transcription factors. These data provide new insights into mycorrhizal signalling mechanisms, supporting a role for NSP1, and specialisation for NSP1-dependent and -independent pathways downstream of DMI3. Our data also indicate significant Myc-LCO-activated signalling upstream of DMI3 and/or parallel to the CSSP and some constitutive activity of the CSSP.

Deep Sequencing of the Medicago truncatula Root Transcriptome Reveals a Massive and Early Interaction between Nodulation Factor and Ethylene Signals

The legume-rhizobium symbiosis is initiated through the activation of the Nodulation (Nod) factor-signaling cascade, leading to a rapid reprogramming of host cell developmental pathways. In this work, we combine transcriptome sequencing with molecular genetics and network analysis to quantify and categorize the transcriptional changes occurring in roots of Medicago truncatula from minutes to days after inoculation with Sinorhizobium medicae. To identify the nature of the inductive and regulatory cues, we employed mutants with absent or decreased Nod factor sensitivities (i.e. Nodulation factor perception and Lysine motif domain-containing receptor-like kinase3, respectively) and an ethylene (ET)-insensitive, Nod factor-hypersensitive mutant (sickle). This unique data set encompasses nine time points, allowing observation of the symbiotic regulation of diverse biological processes with high temporal resolution. Among the many outputs of the study is the early Nod factor-induced, ET-regulated expression of ET signaling and biosynthesis genes. Coupled with the observation of massive transcriptional derepression in the ET-insensitive background, these results suggest that Nod factor signaling activates ET production to attenuate its own signal. Promoter:β-glucuronidase fusions report ET biosynthesis both in root hairs responding to rhizobium as well as in meristematic tissue during nodule organogenesis and growth, indicating that ET signaling functions at multiple developmental stages during symbiosis. In addition, we identified thousands of novel candidate genes undergoing Nod factor-dependent, ET-regulated expression. We leveraged the power of this large data set to model Nod factor- and ET-regulated signaling networks using MERLIN, a regulatory network inference algorithm. These analyses predict key nodes regulating the biological process impacted by Nod factor perception. We have made these results available to the research community through a searchable online resource.

Rhizobial symbiosis effect on the growth, metal uptake, and antioxidant responses of Medicago lupulina under copper stress

The effects of rhizobial symbiosis on the growth, metal uptake, and antioxidant responses of Medicago lupulina in the presence of 200 mg kg(-1) Cu(2+) throughout different stages of symbiosis development were studied. The symbiosis with Sinorhizobium meliloti CCNWSX0020 induced an increase in plant growth and nitrogen content irrespective of the presence of Cu(2+). The total amount of Cu uptake of inoculated plants significantly increased by 34.0 and 120.4% in shoots and roots, respectively, compared with non-inoculated plants. However, although the rhizobial symbiosis promoted Cu accumulation both in shoots and roots, the increase in roots was much higher than in shoots, thus decreasing the translocation factor and helping Cu phytostabilization. The rate of lipid peroxidation was significantly decreased in both shoots and roots of inoculated vs. non-inoculated plants when measured either 8, 13, or 18 days post-inoculation. In comparison with non-inoculated plants, the activities of superoxide dismutase and ascorbate peroxidase of shoots of inoculated plants exposed to excess Cu were significantly elevated at different stages of symbiosis development; similar increases occurred in the activities of superoxide dismutase, catalase, and glutathione reductase of inoculated roots. The symbiosis with S. meliloti CCNWSX0020 also upregulated the corresponding genes involved in antioxidant responses in the plants treated with excess Cu. The results indicated that the rhizobial symbiosis with S. meliloti CCNWSX0020 not only enhanced plant growth and metal uptake but also improved the responses of plant antioxidant defense to excess Cu stress.

Activation of symbiosis signaling by arbuscular mycorrhizal fungi in legumes and rice

Establishment of arbuscular mycorrhizal interactions involves plant recognition of diffusible signals from the fungus, including lipochitooligosaccharides (LCOs) and chitooligosaccharides (COs). Nitrogen-fixing rhizobial bacteria that associate with leguminous plants also signal to their hosts via LCOs, the so-called Nod factors. Here, we have assessed the induction of symbiotic signaling by the arbuscular mycorrhizal (Myc) fungal-produced LCOs and COs in legumes and rice (Oryza sativa). We show that Myc-LCOs and tetra-acetyl chitotetraose (CO4) activate the common symbiosis signaling pathway, with resultant calcium oscillations in root epidermal cells of Medicago truncatula and Lotus japonicus. The nature of the calcium oscillations is similar for LCOs produced by rhizobial bacteria and by mycorrhizal fungi; however, Myc-LCOs activate distinct gene expression. Calcium oscillations were activated in rice atrichoblasts by CO4, but not the Myc-LCOs, whereas a mix of CO4 and Myc-LCOs activated calcium oscillations in rice trichoblasts. In contrast, stimulation of lateral root emergence occurred following treatment with Myc-LCOs, but not CO4, in M. truncatula, whereas both Myc-LCOs and CO4 were active in rice. Our work indicates that legumes and non-legumes differ in their perception of Myc-LCO and CO signals, suggesting that different plant species respond to different components in the mix of signals produced by arbuscular mycorrhizal fungi.

Mutualism-parasitism paradigm synthesized from results of root-endophyte models

Plant tissues host a variety of fungi. One important group is the dark septate endophytes (DSEs) that colonize plant roots and form characteristic intracellular structures - melanized hyphae and microsclerotia. The DSE associations are common and frequently observed in various biomes and plant taxa. Reviews suggest that the proportion of plant species colonized by DSE equal that colonized by AM and microscopic studies show that the proportion of the root system colonized by fungi DSE can equal, or even exceed, the colonization by AM fungi. Despite the high frequency and suspected ecological importance, the effects of DSE colonization on plant growth and performance have remained unclear. Here, we draw from over a decade of experimentation with the obscure DSE symbiosis and synthesize across large bodies of published and unpublished data from Arabidopsis thaliana and Allium porrum model systems as well as from experiments that use native plants to better resolve the host responses to DSE colonization. The data indicate similar distribution of host responses in model and native plant studies, validating the use of model plants for tractable dissection of DSE symbioses. The available data also permit empirical testing of the environmental modulation of host responses to DSE colonization and refining the "mutualism-parasitism-continuum" paradigm for DSE symbioses. These data highlight the context dependency of the DSE symbioses: not only plant species but also ecotypes vary in their responses to populations of conspecific DSE fungi - environmental conditions further shift the host responses similar to those predicted based on the mutualism-parasitism-continuum paradigm. The model systems provide several established avenues of inquiry that permit more detailed molecular and functional dissection of fungal endophyte symbioses, identifying thus likely mechanisms that may underlie the observed host responses to endophyte colonization.

Two CCAAT-box-binding transcription factors redundantly regulate early steps of the legume-rhizobia endosymbiosis

During endosymbiotic interactions between legume plants and nitrogen-fixing rhizobia, successful root infection by bacteria and nodule organogenesis requires the perception and transduction of bacterial lipo-chitooligosaccharidic signal called Nod factor (NF). NF perception in legume roots leads to the activation of an early signaling pathway and of a set of symbiotic genes which is controlled by specific early transcription factors (TFs) including CYCLOPS/IPD3, NSP1, NSP2, ERN1 and NIN. In this study, we bring convincing evidence that the Medicago truncatula CCAAT-box-binding NF-YA1 TF, previously associated with later stages of rhizobial infection and nodule meristem formation is, together with its closest homolog NF-YA2, also an essential positive regulator of the NF-signaling pathway. Here we show that NF-YA1 and NF-YA2 are both expressed in epidermal cells responding to NFs and their knock-down by reverse genetic approaches severely affects the NF-induced expression of symbiotic genes and rhizobial infection. Further over-expression, transactivation and ChIP-PCR approaches indicate that NF-YA1 and NF-YA2 function, at least in part, via the direct activation of ERN1. We thus propose a model in which NF-YA1 and NF-YA2 appear as early symbiotic regulators acting downstream of DMI3 and NIN and possibly within the same regulatory complexes as NSP1/2 to directly activate the expression of ERN1.