MucR Is Required for Transcriptional Activation of Conserved Ion Transporters to Support Nitrogen Fixation of Sinorhizobium fredii in Soybean Nodules

Metal nutrition and transport in the process of symbiotic nitrogen fixation

Symbiotic nitrogen fixation (SNF) facilitated by the interaction between legumes and rhizobia is a well-documented and eco-friendly alternative to chemical nitrogen fertilizers. Host plants obtain fixed nitrogen from rhizobia by providing carbon and mineral nutrients. These mineral nutrients, which are mostly in the form of metal ions, are implicated in various stages of the SNF process. This review describes the functional roles played by metal ions in nodule formation and nitrogen fixation and specifically addresses their transport mechanisms and associated transporters within root nodules. Future research directions and potential strategies for enhancing SNF efficiency are also discussed.

Brucella MucR acts as an H-NS-like protein to silence virulence genes and structure the nucleoid

Histone-like nucleoid structuring (H-NS) and H-NS-like proteins serve as global gene silencers and work with antagonistic transcriptional activators (counter-silencers) to properly coordinate the expression of virulence genes in pathogenic bacteria. In Brucella, MucR has been proposed as a novel H-NS-like gene silencer, but direct experimental evidence is lacking. Here, we show that MucR serves as an H-NS-like silencer of the Brucella abortus genes encoding the polar autotransporter adhesins BtaE and BmaC, the c-di-GMP-specific phosphodiesterase BpdB, and the quorum-sensing regulator BabR. We also demonstrate that the MarR-type transcriptional activator MdrA can displace MucR from the btaE promoter, supporting the existence of MucR counter-silencers in Brucella. Moreover, our chromatin immunoprecipitation (ChIP)-seq analysis identified 546 MucR enrichment peaks along the genome, including in the promoters of the genes encoding the Type IV secretion machinery and effectors and the quorum-sensing regulator VjbR. Importantly, MucR ChIP-seq peaks overlap with the previously described binding sites for the transcriptional activators VjbR, BvrR, and CtrA suggesting that these regulators serve as MucR counter-silencers and work in concert with MucR to coordinate virulence gene expression in Brucella. In addition, using chromosome conformation capture (Hi-C), we show that like H-NS in Escherichia coli, MucR alters the global structure of the Brucella nucleoid. Finally, a copy of the E. coli hns rescues the distinctive growth defect and elevated btaE expression of a B. abortus mucR mutant. Together, these findings solidify the role of MucR as a novel type of H-NS-like protein and suggest that MucR's gene-silencing properties play a key role in virulence in Brucella. IMPORTANCE Histone-like nucleoid structuring (H-NS) and H-NS-like proteins coordinate host-associated behaviors in many pathogenic bacteria, often through forming silencer/counter-silencer pairs with signal-responsive transcriptional activators to tightly control gene expression. Brucella and related bacteria do not encode H-NS or homologs of known H-NS-like proteins, and it is unclear if they have other proteins that perform analogous functions during pathogenesis. In this work, we provide compelling evidence for the role of MucR as a novel H-NS-like protein in Brucella. We show that MucR possesses many of the known functions attributed to H-NS and H-NS-like proteins, including the formation of silencer/counter-silencer pairs to control virulence gene expression and global structuring of the nucleoid. These results uncover a new role for MucR as a nucleoid structuring protein and support the importance of temporal control of gene expression in Brucella and related bacteria.

Rhizobial migration toward roots mediated by FadL-ExoFQP modulation of extracellular long-chain AHLs

Migration from rhizosphere to rhizoplane is a key selecting process in root microbiome assembly, but not fully understood. Rhizobiales members are overrepresented in the core root microbiome of terrestrial plants, and here we report a genome-wide transposon-sequencing of rhizoplane fitness genes of beneficial Sinorhizobium fredii on wild soybean, cultivated soybean, rice, and maize. There were few genes involved in broad-host-range rhizoplane colonization. The fadL mutant lacking a fatty acid transporter exhibited high colonization rates, while mutations in exoFQP (encoding membrane proteins directing exopolysaccharide polymerization and secretion), but not those in exo genes essential for exopolysaccharide biosynthesis, led to severely impaired colonization rates. This variation was not explainable by their rhizosphere and rhizoplane survivability, and associated biofilm and exopolysaccharide production, but consistent with their migration ability toward rhizoplane, and associated surface motility and the mixture of quorum-sensing AHLs (N-acylated-L-homoserine lactones). Genetics and physiology evidences suggested that FadL mediated long-chain AHL uptake while ExoF mediated the secretion of short-chain AHLs which negatively affected long-chain AHL biosynthesis. The fadL and exoF mutants had elevated and depleted extracellular long-chain AHLs, respectively. A synthetic mixture of long-chain AHLs mimicking that of the fadL mutant can improve rhizobial surface motility. When this AHL mixture was spotted into rhizosphere, the migration toward roots and rhizoplane colonization of S. fredii were enhanced in a diffusible way. This work adds novel parts managing extracellular AHLs, which modulate bacterial migration toward rhizoplane. The FadL-ExoFQP system is conserved in Alphaproteobacteria and may shape the "home life" of diverse keystone rhizobacteria.

Intracellular common gardens reveal niche differentiation in transposable element community during bacterial adaptive evolution

The distribution and abundance of transposable elements across the tree of life have significantly shaped the evolution of cellular organisms, but the underlying mechanisms shaping these ecological patterns remain elusive. Here we establish a "common garden" approach to study causal ecological interactions between a xenogeneic conditional lethal sacB gene and the community of transposable insertion sequences (ISs) in a multipartite prokaryote genome. Xenogeneic sacB of low, medium, or high GC content was individually inserted into three replicons of a model bacterium Sinorhizobium fredii, and exhibited replicon- and GC-dependent variation in genetic stability. This variation was largely attributable to multidimensional niche differentiation for IS community members. The transposition efficiency of major active ISs depended on the nucleoid-associated xenogeneic silencer MucR. Experimentally eliminating insertion activity of specific ISs by deleting MucR strongly demonstrated a dominant role of niche differentiation among ISs. This intracellular common garden approach in the experimental evolution context allows not only for evaluating genetic stability of natural and synthetic xenogeneic genes of different sequence signatures in host cells but also for tracking and testing causal relationships in unifying ecological principles in genome ecology.

Adaptive Evolution of Rhizobial Symbiosis beyond Horizontal Gene Transfer: From Genome Innovation to Regulation Reconstruction

There are ubiquitous variations in symbiotic performance of different rhizobial strains associated with the same legume host in agricultural practices. This is due to polymorphisms of symbiosis genes and/or largely unexplored variations in integration efficiency of symbiotic function. Here, we reviewed cumulative evidence on integration mechanisms of symbiosis genes. Experimental evolution, in concert with reverse genetic studies based on pangenomics, suggests that gain of the same circuit of key symbiosis genes through horizontal gene transfer is necessary but sometimes insufficient for bacteria to establish an effective symbiosis with legumes. An intact genomic background of the recipient may not support the proper expression or functioning of newly acquired key symbiosis genes. Further adaptive evolution, through genome innovation and reconstruction of regulation networks, may confer the recipient of nascent nodulation and nitrogen fixation ability. Other accessory genes, either co-transferred with key symbiosis genes or stochastically transferred, may provide the recipient with additional adaptability in ever-fluctuating host and soil niches. Successful integrations of these accessory genes with the rewired core network, regarding both symbiotic and edaphic fitness, can optimize symbiotic efficiency in various natural and agricultural ecosystems. This progress also sheds light on the development of elite rhizobial inoculants using synthetic biology procedures.

The Ros/MucR Zinc-Finger Protein Family in Bacteria: Structure and Functions

Ros/MucR is a widespread family of bacterial zinc-finger-containing proteins that integrate multiple functions, such as symbiosis, virulence, transcription regulation, motility, production of surface components, and various other physiological processes in cells. This regulatory protein family is conserved in bacteria and is characterized by its zinc-finger motif, which has been proposed as the ancestral domain from which the eukaryotic C₂H₂ zinc-finger structure has evolved. The first prokaryotic zinc-finger domain found in the transcription regulator Ros was identified in Agrobacterium tumefaciens. In the past decades, a large body of evidence revealed Ros/MucR as pleiotropic transcriptional regulators that mainly act as repressors through oligomerization and binding to AT-rich target promoters. The N-terminal domain and the zinc-finger-bearing C-terminal region of these regulatory proteins are engaged in oligomerization and DNA binding, respectively. These properties of the Ros/MucR proteins are similar to those of xenogeneic silencers, such as H-NS, MvaT, and Lsr2, which are mainly found in other lineages. In fact, a novel functional model recently proposed for this protein family suggests that they act as H-NS-'like' gene silencers. The prokaryotic zinc-finger domain exhibits interesting structural and functional features that are different from that of its eukaryotic counterpart (a βββα topology), as it folds in a significantly larger zinc-binding globular domain (a βββαα topology). Phylogenetic analysis of Ros/MucR homologs suggests an ancestral origin of this type of protein in α-Proteobacteria. Furthermore, multiple duplications and lateral gene transfer events contributing to the diversity and phyletic distribution of these regulatory proteins were found in bacterial genomes.

Coordinated regulation of symbiotic adaptation by NodD proteins and NolA in the type I peanut bradyrhizobial strain Bradyrhizobium zhanjiangense CCBAU51778

Type I peanut bradyrhizobial strains can establish efficient symbiosis in contrast to symbiotic incompatibility induced by type II strains with mung bean. The notable distinction in the two kinds of key symbiosis-related regulators nolA and nodD close to the nodABCSUIJ operon region between these two types of peanut bradyrhizobia was found. Therefore, we determined whether NolA and NodD proteins regulate the symbiotic adaptations of type I strains to different hosts. We found that NodD1-NolA synergistically regulated the symbiosis between the type I strain Bradyrhizobium zhanjiangense CCBAU51778 and mung bean, and NodD1-NodD2 jointly regulated nodulation ability. In contrast, NodD1-NolA coordinately regulated nodulation ability in the CCBAU51778-peanut symbiosis. Meanwhile, NodD1 and NolA collectively contributes to competitive nodule colonization of CCBAU51778 on both hosts. The Fucosylated Nod factors and intact type 3 secretion system (T3SS), rather than extra nodD2 and full-length nolA, were critical for effective symbiosis with mung bean. Unexpectedly, T3SS-related genes were activated by NodD2 but not NodD1. Compared to NodD1 and NodD2, NolA predominantly inhibits exopolysaccharide production by promoting exoR expression. Importantly, this is the first report that NolA regulates rhizobial T3SS-related genes. The coordinated regulation and integration of different gene networks to fine-tune the expression of symbiosis-related genes and other accessory genes by NodD1-NolA might be required for CCBAU51778 to efficiently nodulate peanut. This study shed new light on our understanding of the regulatory roles of NolA and NodD proteins in symbiotic adaptation, highlighting the sophisticated gene networks dominated by NodD1-NolA.

The convergent xenogeneic silencer MucR predisposes α-proteobacteria to integrate AT-rich symbiosis genes

Bacterial adaptation is largely shaped by horizontal gene transfer, xenogeneic silencing mediated by lineage-specific DNA bridgers (H-NS, Lsr2, MvaT and Rok), and various anti-silencing mechanisms. No xenogeneic silencing DNA bridger is known for α-proteobacteria, from which mitochondria evolved. By investigating α-proteobacterium Sinorhizobium fredii, a facultative legume microsymbiont, here we report the conserved zinc-finger bearing MucR as a novel xenogeneic silencing DNA bridger. Self-association mediated by its N-terminal domain (NTD) is required for DNA–MucR–DNA bridging complex formation, maximizing MucR stability, transcriptional silencing, and efficient symbiosis in legume nodules. Essential roles of NTD, CTD (C-terminal DNA-binding domain), or full-length MucR in symbiosis can be replaced by non-homologous NTD, CTD, or full-length protein of H-NS from γ-proteobacterium Escherichia coli, while NTD rather than CTD of Lsr2 from Gram-positive Mycobacterium tuberculosis can replace the corresponding domain of MucR in symbiosis. Chromatin immunoprecipitation sequencing reveals similar recruitment profiles of H-NS, MucR and various functional chimeric xenogeneic silencers across the multipartite genome of S. fredii, i.e. preferring AT-rich genomic islands and symbiosis plasmid with key symbiosis genes as shared targets. Collectively, the convergently evolved DNA bridger MucR predisposed α-proteobacteria to integrate AT-rich foreign DNA including symbiosis genes, horizontal transfer of which is strongly selected in nature.

Surface Motility Regulation of Sinorhizobium fredii HH103 by Plant Flavonoids and the NodD1, TtsI, NolR, and MucR1 Symbiotic Bacterial Regulators

Bacteria can spread on surfaces to colonize new environments and access more resources. Rhizobia, a group of α- and β-Proteobacteria, establish nitrogen-fixing symbioses with legumes that rely on a complex signal interchange between the partners. Flavonoids exuded by plant roots and the bacterial transcriptional activator NodD control the transcription of different rhizobial genes (the so-called nod regulon) and, together with additional bacterial regulatory proteins (such as TtsI, MucR or NolR), influence the production of different rhizobial molecular signals. In Sinorhizobium fredii HH103, flavonoids and NodD have a negative effect on exopolysaccharide production and biofilm production. Since biofilm formation and motility are often inversely regulated, we have analysed whether flavonoids may influence the translocation of S. fredii HH103 on surfaces. We show that the presence of nod gene-inducing flavonoids does not affect swimming but promotes a mode of surface translocation, which involves both flagella-dependent and -independent mechanisms. This surface motility is regulated in a flavonoid-NodD1-TtsI-dependent manner, relies on the assembly of the symbiotic type 3 secretion system (T3SS), and involves the participation of additional modulators of the nod regulon (NolR and MucR1). To our knowledge, this is the first evidence indicating the participation of T3SS in surface motility in a plant-interacting bacterium. Interestingly, flavonoids acting as nod-gene inducers also participate in the inverse regulation of surface motility and biofilm formation, which could contribute to a more efficient plant colonisation.

The zinc-finger bearing xenogeneic silencer MucR in α-proteobacteria balances adaptation and regulatory integrity

Foreign AT-rich genes drive bacterial adaptation to new niches while challenging the existing regulation network. Here we report that MucR, a conserved regulator in α-proteobacteria, balances adaptation and regulatory integrity in Sinorhizobium fredii, a facultative microsymbiont of legumes. Chromatin immunoprecipitation sequencing coupled with transcriptomic data reveal that average transcription levels of both target and non-target genes, under free-living and symbiotic conditions, increase with their conservation levels. Targets involved in environmental adaptation and symbiosis belong to genus or species core and can be repressed or activated by MucR in a condition-dependent manner, implying regulatory integrations. However, most targets are enriched in strain-specific genes of lower expression levels and higher AT%. Within each conservation levels, targets have higher AT% and average transcription levels than non-target genes and can be further up-regulated in the mucR mutant. This is consistent with higher AT% of spacers between -35 and -10 elements of promoters for target genes, which enhances transcription. The MucR recruitment level linearly increases with AT% and the number of a flexible pattern (with periodic repeats of Ts) of target sequences. Collectively, MucR directly represses AT-rich foreign genes with predisposed high transcription potential while progressive erosions of its target sites facilitate regulatory integrations of foreign genes.

Global Transcriptional Repression of Diguanylate Cyclases by MucR1 Is Essential for Sinorhizobium-Soybean Symbiosis

The ubiquitous bacterial second messenger c-di-GMP is intensively studied in pathogens but less so in mutualistic bacteria. Here, we report a genome-wide investigation of functional diguanylate cyclases (DGCs) synthesizing c-di-GMP from two molecules of GTP in Sinorhizobium fredii CCBAU45436, a facultative microsymbiont fixing nitrogen in nodules of diverse legumes, including soybean. Among 25 proteins harboring a putative GGDEF domain catalyzing the biosynthesis of c-di-GMP, eight functional DGCs were identified by heterogenous expression in Escherichia coli in a Congo red binding assay. This screening result was further verified by in vitro enzymatic assay with purified full proteins or the GGDEF domains from representative functional and nonfunctional DGCs. In the same in vitro assay, a functional EAL domain catalyzing the degradation of c-di-GMP into pGpG was identified in a protein that has an inactive GGDEF domain but with an active phosphodiesterase (PDE) function. The identified functional DGCs generally exhibited low transcription levels in soybean nodules compared to free-living cultures, as revealed in transcriptomes. An engineered upregulation of a functional DGC in nodules led to a significant increase of c-di-GMP level and symbiotic defects, which were not observed when a functional EAL domain was upregulated at the same level. Further transcriptional analysis and gel shift assay demonstrated that these functional DGCs were all transcriptionally repressed in nodules by a global pleiotropic regulator, MucR1, that is essential in Sinorhizobium-soybean symbiosis. These findings shed novel insights onto the systematic regulation of c-di-GMP biosynthesis in mutualistic symbiosis. IMPORTANCE The ubiquitous second messenger c-di-GMP is well-known for its role in biofilm formation and host adaptation of pathogens, whereas it is less investigated in mutualistic symbioses. Here, we reveal a cocktail of eight functional diguanylate cyclases (DGCs) catalyzing the biosynthesis of c-di-GMP in a broad-host-range Sinorhizobium that can establish nitrogen-fixing nodules on soybean and many other legumes. These functional DGCs are generally transcribed at low levels in soybean nodules compared to free-living conditions. The engineered nodule-specific upregulation of DGC can elevate the c-di-GMP level and cause symbiotic defects, while the upregulation of a phosphodiesterase that quenches c-di-GMP has no detectable symbiotic defects. Moreover, eight functional DGCs located on two different replicons are all directly repressed in nodules by a global silencer, MucR1, that is essential for Sinorhizobium-soybean symbiosis. These findings represent a novel mechanism of a strategic regulation of the c-di-GMP biosynthesis arsenal in prokaryote-eukaryote interactions.

Rhizobial Exopolysaccharides: Genetic Regulation of Their Synthesis and Relevance in Symbiosis with Legumes

Rhizobia are soil proteobacteria able to engage in a nitrogen-fixing symbiotic interaction with legumes that involves the rhizobial infection of roots and the bacterial invasion of new organs formed by the plant in response to the presence of appropriate bacterial partners. This interaction relies on a complex molecular dialogue between both symbionts. Bacterial N-acetyl-glucosamine oligomers called Nod factors are indispensable in most cases for early steps of the symbiotic interaction. In addition, different rhizobial surface polysaccharides, such as exopolysaccharides (EPS), may also be symbiotically relevant. EPS are acidic polysaccharides located out of the cell with little or no cell association that carry out important roles both in free-life and in symbiosis. EPS production is very complexly modulated and, frequently, co-regulated with Nod factors, but the type of co-regulation varies depending on the rhizobial strain. Many studies point out a signalling role for EPS-derived oligosaccharides in root infection and nodule invasion but, in certain symbiotic couples, EPS can be dispensable for a successful interaction. In summary, the complex regulation of the production of rhizobial EPS varies in different rhizobia, and the relevance of this polysaccharide in symbiosis with legumes depends on the specific interacting couple.

Lineage-Specific Rewiring of Core Pathways Predating Innovation of Legume Nodules Shapes Symbiotic Efficiency

The interkingdom coevolution innovated the rhizobium-legume symbiosis. The application of this nitrogen-fixing system in sustainable agriculture is usually impeded by incompatible interactions between partners. However, the progressive evolution of rhizobium-legume compatibility remains elusive. In this work, deletions of rhcV encoding a structural component of the type three secretion system allow related Sinorhizobium strains to nodulate a previously incompatible soybean cultivar (Glycine max). These rhcV mutants show low to medium to high symbiotic efficiency on the same cultivated soybean while being indistinguishable on wild soybean plants (Glycine soja). The dual pantranscriptomics reveals nodule-specific activation of core symbiosis genes of Sinorhizobium and Glycine genes associated with genome duplication events along the chronogram. Unexpectedly, symbiotic efficiency is in line with lineage-dependent transcriptional profiles of core pathways which predate the diversification of Fabaceae and Sinorhizobium. This is supported by further physiological and biochemical experiments. Particularly, low-efficiency nodules show disordered antioxidant activity and low-energy status, which restrict nitrogen fixation activity. Collectively, the ancient core pathways play a crucial role in optimizing the function of later-evolved mutualistic arsenals in the rhizobium-legume coevolution.

IMPORTANCE Significant roles of complex extracellular microbiota in environmental adaptation of eukaryotes in ever-changing circumstances have been revealed. Given the intracellular infection ability, facultative endosymbionts can be considered pioneers within complex extracellular microbiota and are ideal organisms for understanding the early stage of interkingdom adaptation. This work reveals that the later innovation of key symbiotic arsenals and the lineage-specific network rewiring in ancient core pathways, predating the divergence of legumes and rhizobia, underline the progressive evolution of rhizobium-legume compatibility. This insight not only is significant for improving the application benefits of rhizobial inoculants in sustainable agriculture but also advances our general understanding of the interkingdom coevolution which is theoretically explored by all host-microbiota interactions.

PHO1 family members transport phosphate from infected nodule cells to bacteroids in Medicago truncatula

Legumes play an important role in the soil nitrogen availability via symbiotic nitrogen fixation (SNF). Phosphate (Pi) deficiency severely impacts SNF because of the high Pi requirement of symbiosis. Whereas PHT1 transporters are involved in Pi uptake into nodules, it is unknown how Pi is transferred from the plant infected cells to nitrogen-fixing bacteroids. We hypothesized that Medicago truncatula genes homologous to Arabidopsis PHO1, encoding a vascular apoplastic Pi exporter, are involved in Pi transfer to bacteroids. Among the seven MtPHO1 genes present in M. truncatula, we found that two genes, namely MtPHO1.1 and MtPHO1.2, were broadly expressed across the various nodule zones in addition to the root vascular system. Expressions of MtPHO1.1 and MtPHO1.2 in Nicotiana benthamiana mediated specific Pi export. Plants with nodule-specific downregulation of both MtPHO1.1 and MtPHO1.2 were generated by RNA interference (RNAi) to examine their roles in nodule Pi homeostasis. Nodules of RNAi plants had lower Pi content and a three-fold reduction in SNF, resulting in reduced shoot growth. Whereas the rate of 33Pi uptake into nodules of RNAi plants was similar to control, transfer of 33Pi from nodule cells into bacteroids was reduced and bacteroids activated their Pi-deficiency response. Our results implicate plant MtPHO1 genes in bacteroid Pi homeostasis and SNF via the transfer of Pi from nodule infected cells to bacteroids.

Ancestral zinc-finger bearing protein MucR in alpha-proteobacteria: A novel xenogeneic silencer?

The MucR/Ros family protein is conserved in alpha-proteobacteria and characterized by its zinc-finger motif that has been proposed as the ancestral domain from which the eukaryotic C2H2 zinc-finger structure evolved. In the past decades, accumulated evidences have revealed MucR as a pleiotropic transcriptional regulator that integrating multiple functions such as virulence, symbiosis, cell cycle and various physiological processes. Scattered reports indicate that MucR mainly acts as a repressor, through oligomerization and binding to multiple sites of AT-rich target promoters. The N-terminal region and zinc-finger bearing C-terminal region of MucR mediate oligomerization and DNA-binding, respectively. These features are convergent to those of xenogeneic silencers such as H-NS, MvaT, Lsr2 and Rok, which are mainly found in other lineages. Phylogenetic analysis of MucR homologs suggests an ancestral origin of MucR in alpha- and delta-proteobacteria. Multiple independent duplication and lateral gene transfer events contribute to the diversity and phyletic distribution of MucR. Finally, we posed questions which remain unexplored regarding the putative roles of MucR as a xenogeneic silencer and a general manager in balancing adaptation and regulatory integration in the pangenome context.

Structural Insight of the Full-Length Ros Protein: A Prototype of the Prokaryotic Zinc-Finger Family

Ros/MucR is a widespread family of bacterial zinc-finger (ZF) containing proteins that integrate multiple functions such as virulence, symbiosis and/or cell cycle transcription. NMR solution structure of Ros DNA-binding domain (region 56–142, i.e. Ros87) has been solved by our group and shows that the prokaryotic ZF domain shows interesting structural and functional features that differentiate it from its eukaryotic counterpart as it folds in a significantly larger zinc-binding globular domain. We have recently proposed a novel functional model for this family of proteins suggesting that they may act as H-NS-‘like’ gene silencers. Indeed, the N-terminal region of this family of proteins appears to be responsible for the formation of functional oligomers. No structural characterization of the Ros N-terminal domain (region 1–55) is available to date, mainly because of serious solubility problems of the full-length protein. Here we report the first structural characterization of the N-terminal domain of the prokaryotic ZF family examining by means of MD and NMR the structural preferences of the full-length Ros protein from Agrobacterium tumefaciens.

Transition metal transporters in rhizobia: tuning the inorganic micronutrient requirements to different living styles

A group of bacteria known as rhizobia are key players in symbiotic nitrogen fixation (SNF) in partnership with legumes. After a molecular exchange, the bacteria end surrounded by a plant membrane forming symbiosomes, organelle-like structures, where they differentiate to bacteroids and fix nitrogen. This symbiotic process is highly dependent on dynamic nutrient exchanges between the partners. Among these are transition metals (TM) participating as inorganic and organic cofactors of fundamental enzymes. While the understanding of how plant transporters facilitate TMs to the very near environment of the bacteroid is expanding, our knowledge on how bacteroid transporters integrate to TM homeostasis mechanisms in the plant host is still limited. This is significantly relevant considering the low solubility and scarcity of TMs in soils, and the in crescendo gradient of TM bioavailability rhizobia faces during the infection and bacteroid differentiation processes. In the present work, we review the main metal transporter families found in rhizobia, their role in free-living conditions and, when known, in symbiosis. We focus on discussing those transporters which could play a significant role in TM-dependent biochemical and physiological processes in the bacteroid, thus paving the way towards an optimized SNF.

Experimental Evolution of Legume Symbionts: What Have We Learnt?

Rhizobia, the nitrogen-fixing symbionts of legumes, are polyphyletic bacteria distributed in many alpha- and beta-proteobacterial genera. They likely emerged and diversified through independent horizontal transfers of key symbiotic genes. To replay the evolution of a new rhizobium genus under laboratory conditions, the symbiotic plasmid of Cupriavidus taiwanensis was introduced in the plant pathogen Ralstonia solanacearum, and the generated proto-rhizobium was submitted to repeated inoculations to the C. taiwanensis host, Mimosa pudica L.. This experiment validated a two-step evolutionary scenario of key symbiotic gene acquisition followed by genome remodeling under plant selection. Nodulation and nodule cell infection were obtained and optimized mainly via the rewiring of regulatory circuits of the recipient bacterium. Symbiotic adaptation was shown to be accelerated by the activity of a mutagenesis cassette conserved in most rhizobia. Investigating mutated genes led us to identify new components of R. solanacearum virulence and C. taiwanensis symbiosis. Nitrogen fixation was not acquired in our short experiment. However, we showed that post-infection sanctions allowed the increase in frequency of nitrogen-fixing variants among a non-fixing population in the M. pudica-C. taiwanensis system and likely allowed the spread of this trait in natura. Experimental evolution thus provided new insights into rhizobium biology and evolution.

Modulation of Symbiotic Compatibility by Rhizobial Zinc Starvation Machinery

The rhizobium-legume symbiosis contributes around 65% of biological nitrogen fixation in agriculture systems and is critical for sustainable agriculture by reducing the amount of chemical nitrogen fertilizer being used. Rhizobial inocula have been commercialized for more than 100 years, but the efficiency of inoculation can vary among legume cultivars, field sites, and years. These long-lasting challenging problems impede the establishment of a sustainable agriculture, particularly in developing countries. Here, we report that rhizobial zinc starvation machinery containing a conserved high-affinity zinc transporter and accessory components makes cumulative contributions to modulating rhizobial symbiotic compatibility. This work highlights a critical role of largely unexplored nutritional immunity in the rhizobium-legume symbiosis, which makes zinc starvation machinery an attractive target for improving rhizobial symbiotic compatibility.

Pathogenic bacteria need high-affinity zinc uptake systems to counteract the nutritional immunity exerted by infected hosts. However, our understanding of zinc homeostasis in mutualistic systems such as the rhizobium-legume symbiosis is limited. Here, we show that the conserved high-affinity zinc transporter ZnuABC and accessory transporter proteins (Zip1, Zip2, and c06450) made cumulative contributions to nodulation of the broad-host-range strain Sinorhizobium fredii CCBAU45436. Zur acted as a zinc-dependent repressor for the znuC-znuB-zur operon, znuA, and c06450 by binding to the associated Zur box, but did not regulate zip1 and zip2. ZnuABC was the major zinc transporter. Combined mutants lacking znuA and one of the three accessory genes had more severe defects in nodulation and growth under zinc starvation conditions than the znuA mutant, though rhizoplane colonization by these mutants was not impaired. In contrast to the elite strain CCBAU45436, more drastic symbiotic defects were observed for the znuA mutants of other Sinorhizobium strains, which lack at least one of the three accessory genes in their genomes and are characterized by their limited host range and geographical distribution. The znu-derived mutants showed a higher expression level of nod genes involved in Nod factor biosynthesis and a reduced expression of genes encoding a type three secretion system and its effector NopP, which can interfere with the host immune system. Application of exogenous zinc restored the nodulation ability of these znu-derived mutants. Therefore, the conserved ZnuABC and accessory components in the zinc starvation machinery play an important role in modulating symbiotic compatibility.

Metabolic Analyses of Nitrogen Fixation in the Soybean Microsymbiont Sinorhizobium fredii Using Constraint-Based Modeling

Nitrogen is the most limiting macronutrient for plant growth, and rhizobia are important bacteria for agriculture because they can fix atmospheric nitrogen and make it available to legumes through the establishment of a symbiotic relationship with their host plants. In this work, we studied the nitrogen fixation process in the microsymbiont Sinorhizobium fredii at the genome level. A metabolic model was built using genome annotation and literature to reconstruct the symbiotic form of S. fredii. Genes controlling the nitrogen fixation process were identified by simulating gene knockouts. Additionally, the nitrogen-fixing capacities of S. fredii CCBAU45436 in symbiosis with cultivated and wild soybeans were evaluated. The predictions suggested an outperformance of S. fredii with cultivated soybean, consistent with published experimental evidence. The reconstruction presented here will help to understand and improve nitrogen fixation capabilities of S. fredii and will be beneficial for agriculture by reducing the reliance on fertilizer applications.

Rhizobia are soil bacteria able to establish symbiosis with diverse host plants. Specifically, Sinorhizobium fredii is a soil bacterium that forms nitrogen-fixing root nodules in diverse legumes, including soybean. The strain S. fredii CCBAU45436 is a dominant sublineage of S. fredii that nodulates soybeans in alkaline-saline soils in the Huang-Huai-Hai Plain region of China. Here, we present a manually curated metabolic model of the symbiotic form of Sinorhizobium fredii CCBAU45436. A symbiosis reaction was defined to describe the specific soybean-microsymbiont association. The performance and quality of the reconstruction had a 70% score when assessed using a standardized genome-scale metabolic model test suite. The model was used to evaluate in silico single-gene knockouts to determine the genes controlling the nitrogen fixation process. One hundred forty-one of 541 genes (26%) were found to influence the symbiotic process, wherein 121 genes were predicted as essential and 20 others as having a partial effect. Transcriptomic profiles of CCBAU45436 were used to evaluate the nitrogen fixation capacity in cultivated versus in wild soybean inoculated with the microsymbiont. The model quantified the nitrogen fixation activities of the strain in these two hosts and predicted a higher nitrogen fixation capacity in cultivated soybean. Our results are consistent with published data demonstrating larger amounts of ureides and total nitrogen in cultivated soybean than in wild soybean. This work presents the first metabolic network reconstruction of S. fredii as an example of a useful tool for exploring the potential benefits of microsymbionts to sustainable agriculture and the ecosystem.

IMPORTANCE Nitrogen is the most limiting macronutrient for plant growth, and rhizobia are important bacteria for agriculture because they can fix atmospheric nitrogen and make it available to legumes through the establishment of a symbiotic relationship with their host plants. In this work, we studied the nitrogen fixation process in the microsymbiont Sinorhizobium fredii at the genome level. A metabolic model was built using genome annotation and literature to reconstruct the symbiotic form of S. fredii. Genes controlling the nitrogen fixation process were identified by simulating gene knockouts. Additionally, the nitrogen-fixing capacities of S. fredii CCBAU45436 in symbiosis with cultivated and wild soybeans were evaluated. The predictions suggested an outperformance of S. fredii with cultivated soybean, consistent with published experimental evidence. The reconstruction presented here will help to understand and improve nitrogen fixation capabilities of S. fredii and will be beneficial for agriculture by reducing the reliance on fertilizer applications.

Sinorhizobium fredii HH103 RirA Is Required for Oxidative Stress Resistance and Efficient Symbiosis with Soybean

Members of Rhizobiaceae contain a homologue of the iron-responsive regulatory protein RirA. In different bacteria, RirA acts as a repressor of iron uptake systems under iron-replete conditions and contributes to ameliorate cell damage during oxidative stress. In Rhizobium leguminosarum and Sinorhizobium meliloti, mutations in rirA do not impair symbiotic nitrogen fixation. In this study, a rirA mutant of broad host range S. fredii HH103 has been constructed (SVQ780) and its free-living and symbiotic phenotypes evaluated. No production of siderophores could be detected in either the wild-type or SVQ780. The rirA mutant exhibited a growth advantage under iron-deficient conditions and hypersensitivity to hydrogen peroxide in iron-rich medium. Transcription of rirA in HH103 is subject to autoregulation and inactivation of the gene upregulates fbpA, a gene putatively involved in iron transport. The S. fredii rirA mutant was able to nodulate soybean plants, but symbiotic nitrogen fixation was impaired. Nodules induced by the mutant were poorly infected compared to those induced by the wild-type. Genetic complementation reversed the mutant's hypersensitivity to H₂O₂, expression of fbpA, and symbiotic deficiency in soybean plants. This is the first report that demonstrates a role for RirA in the Rhizobium-legume symbiosis.

Identifying the region responsible for Brucella abortus MucR higher-order oligomer formation and examining its role in gene regulation

MucR is a member of the Ros/MucR family of prokaryotic zinc-finger proteins found in the α-proteobacteria which regulate the expression of genes required for the successful pathogenic and symbiotic interactions of these bacteria with the eukaryotic hosts. The structure and function of their distinctive zinc-finger domain has been well-studied, but only recently the quaternary structure of the full length proteins was investigated demonstrating their ability to form higher-order oligomers. The aim of this study was to identify the region of MucR involved in higher-order oligomer formation by analysing deletion and point mutants of this protein by Light Scattering, and to determine the role that MucR oligomerization plays in the regulatory function of this protein. Here we demonstrate that a conserved hydrophobic region at the N-terminus of MucR is responsible for higher-order oligomer formation and that MucR oligomerization is essential for its regulatory function in Brucella. All these features of MucR are shared by the histone-like nucleoid structuring protein, (H-NS), leading us to propose that the prokaryotic zinc-finger proteins in the MucR/Ros family control gene expression employing a mechanism similar to that used by the H-NS proteins, rather than working as classical transcriptional regulators.

Evidence for Phosphate Starvation of Rhizobia without Terminal Differentiation in Legume Nodules

Phosphate homeostasis is tightly modulated in all organisms, including bacteria, which harbor both high- and low-affinity transporters acting under conditions of fluctuating phosphate levels. It was thought that nitrogen-fixing rhizobia, named bacteroids, inhabiting root nodules of legumes are not phosphate limited. Here, we show that the high-affinity phosphate transporter PstSCAB, rather than the low-affinity phosphate transporter Pit, is essential for effective nitrogen fixation of Sinorhizobium fredii in soybean nodules. Symbiotic and growth defects of the pst mutant can be effectively restored by knocking out PhoB, the transcriptional repressor of pit. The pst homologs of representative rhizobia were actively transcribed in bacteroids without terminal differentiation in nodules of diverse legumes (soybean, pigeonpea, cowpea, common bean, and Sophora flavescens) but exhibited a basal expression level in terminally differentiated bacteroids (alfalfa, pea, and peanut). Rhizobium leguminosarum bv. viciae Rlv3841 undergoes characteristic nonterminal and terminal differentiations in nodules of S. flavescens and pea, respectively. The pst mutant of Rlv3841 showed impaired adaptation to the nodule environment of S. flavescens but was indistinguishable from the wild-type strain in pea nodules. Taken together, root nodule rhizobia can be either phosphate limited or nonlimited regarding the rhizobial differentiation fate, which is a host-dependent feature.

Extracellular polysaccharide protects Rhizobium leguminosarum cells against zinc stress in vitro and during symbiosis with clover

Rhizobium leguminosarum bv. trifolii is a soil bacterium that establishes symbiosis with clover (Trifolium spp.) under nitrogen-limited conditions. This microorganism produces exopolysaccharide (EPS), which plays an important role in symbiotic interactions with the host plant. The aim of the current study was to establish the role of EPS in the response of R. leguminosarum bv. trifolii cells, free-living and during symbiosis, to zinc stress. We show that EPS-deficient mutants were more sensitive to Zn²⁺ exposure than EPS-producing strains, and that EPS overexpression conferred some protection onto the strains beyond that observed in the wild type. Exposure of the bacteria to Zn²⁺ ions stimulated EPS and biofilm production, and increased cell hydrophobicity. However, zinc stress negatively affected the motility and attachment of bacteria to clover roots, as well as the symbiosis with the host plant. In the presence of Zn²⁺ ions, cell viability, root attachment, biofilm formation and symbiotic efficiency of EPS-overproducing strains were significantly higher than those of the EPS-deficient mutants. We conclude that EPS plays an important role in the adaptation of rhizobia to zinc stress, in both the free-living stage and during symbiosis.

Coordinated regulation of core and accessory genes in the multipartite genome of Sinorhizobium fredii

Prokaryotes benefit from having accessory genes, but it is unclear how accessory genes can be linked with the core regulatory network when developing adaptations to new niches. Here we determined hierarchical core/accessory subsets in the multipartite pangenome (composed of genes from the chromosome, chromid and plasmids) of the soybean microsymbiont Sinorhizobium fredii by comparing twelve Sinorhizobium genomes. Transcriptomes of two S. fredii strains at mid-log and stationary growth phases and in symbiotic conditions were obtained. The average level of gene expression, variation of expression between different conditions, and gene connectivity within the co-expression network were positively correlated with the gene conservation level from strain-specific accessory genes to genus core. Condition-dependent transcriptomes exhibited adaptive transcriptional changes in pangenome subsets shared by the two strains, while strain-dependent transcriptomes were enriched with accessory genes on the chromid. Proportionally more chromid genes than plasmid genes were co-expressed with chromosomal genes, while plasmid genes had a higher within-replicon connectivity in expression than chromid ones. However, key nitrogen fixation genes on the symbiosis plasmid were characterized by high connectivity in both within- and between-replicon analyses. Among those genes with host-specific upregulation patterns, chromosomal znu and mdt operons, encoding a conserved high-affinity zinc transporter and an accessory multi-drug efflux system, respectively, were experimentally demonstrated to be involved in host-specific symbiotic adaptation. These findings highlight the importance of integrative regulation of hierarchical core/accessory components in the multipartite genome of bacteria during niche adaptation and in shaping the prokaryotic pangenome in the long run.

Adaptive evolution of rhizobial symbiotic compatibility mediated by co-evolved insertion sequences

Mutualism between bacteria and eukaryotes has essential roles in the history of life, but the evolution of their compatibility is poorly understood. Here we show that different Sinorhizobium strains can form either nitrogen-fixing nodules or uninfected pseudonodules on certain cultivated soybeans, while being all effective microsymbionts of some wild soybeans. However, a few well-infected nodules can be found on a commercial soybean using inocula containing a mixed pool of Tn5 insertion mutants derived from an incompatible strain. Reverse genetics and genome sequencing of compatible mutants demonstrated that inactivation of T3SS (type three secretion system) accounted for this phenotypic change. These mutations in the T3SS gene cluster were dominated by parallel transpositions of insertion sequences (ISs) other than the introduced Tn5. This genetic and phenotypic change can also be achieved in an experimental evolution scenario on a laboratory time scale using incompatible wild-type strains as inocula. The ISs acting in the adaptive evolution of Sinorhizobium strains exhibit broader phyletic and replicon distributions than other ISs, and prefer target sequences of low GC% content, a characteristic feature of symbiosis plasmid where T3SS genes are located. These findings suggest an important role of co-evolved ISs in the adaptive evolution of rhizobial compatibility.

Transcriptomic Studies of the Effect of nod Gene-Inducing Molecules in Rhizobia: Different Weapons, One Purpose

Simultaneous quantification of transcripts of the whole bacterial genome allows the analysis of the global transcriptional response under changing conditions. RNA-seq and microarrays are the most used techniques to measure these transcriptomic changes, and both complement each other in transcriptome profiling. In this review, we exhaustively compiled the symbiosis-related transcriptomic reports (microarrays and RNA sequencing) carried out hitherto in rhizobia. This review is specially focused on transcriptomic changes that takes place when five rhizobial species, Bradyrhizobium japonicum (=diazoefficiens) USDA 110, Rhizobium leguminosarum biovar viciae 3841, Rhizobium tropici CIAT 899, Sinorhizobium (=Ensifer) meliloti 1021 and S. fredii HH103, recognize inducing flavonoids, plant-exuded phenolic compounds that activate the biosynthesis and export of Nod factors (NF) in all analysed rhizobia. Interestingly, our global transcriptomic comparison also indicates that each rhizobial species possesses its own arsenal of molecular weapons accompanying the set of NF in order to establish a successful interaction with host legumes.

Regulatory Elements Located in the Upstream Region of the Rhizobium leguminosarum rosR Global Regulator Are Essential for Its Transcription and mRNA Stability

Rhizobium leguminosarum bv. trifolii is a soil bacterium capable of establishing a symbiotic relationship with clover (Trifolium spp.). Previously, the rosR gene, encoding a global regulatory protein involved in motility, synthesis of cell-surface components, and other cellular processes was identified and characterized in this bacterium. This gene possesses a long upstream region that contains several regulatory motifs, including inverted repeats (IRs) of different lengths. So far, the role of these motifs in the regulation of rosR transcription has not been elucidated in detail. In this study, we performed a functional analysis of these motifs using a set of transcriptional rosR-lacZ fusions that contain mutations in these regions. The levels of rosR transcription for different mutant variants were evaluated in R. leguminosarum using both quantitative real-time PCR and β-galactosidase activity assays. Moreover, the stability of wild type rosR transcripts and those with mutations in the regulatory motifs was determined using an RNA decay assay and plasmids with mutations in different IRs located in the 5′-untranslated region of the gene. The results show that transcription of rosR undergoes complex regulation, in which several regulatory elements located in the upstream region and some regulatory proteins are engaged. These include an upstream regulatory element, an extension of the -10 element containing three nucleotides TGn (TGn-extended -10 element), several IRs, and PraR repressor related to quorum sensing.

The nitrate-reduction gene cluster components exert lineage-dependent contributions to optimization of Sinorhizobium symbiosis with soybeans

Receiving nodulation and nitrogen fixation genes does not guarantee rhizobia an effective symbiosis with legumes. Here, variations in gene content were determined for three Sinorhizobium species showing contrasting symbiotic efficiency on soybeans. A nitrate-reduction gene cluster absent in S. sojae was found to be essential for symbiotic adaptations of S. fredii and S. sp. III. In S. fredii, the deletion mutation of the nap (nitrate reductase), instead of nir (nitrite reductase) and nor (nitric oxide reductase), led to defects in nitrogen-fixation (Fix- ). By contrast, none of these core nitrate-reduction genes were required for the symbiosis of S. sp. III. However, within the same gene cluster, the deletion of hemN1 (encoding oxygen-independent coproporphyrinogen III oxidase) in both S. fredii and S. sp. III led to the formation of nitrogen-fixing (Fix+ ) but ineffective (Eff- ) nodules. These Fix+ /Eff- nodules were characterized by significantly lower enzyme activity of glutamine synthetase indicating rhizobial modulation of nitrogen-assimilation by plants. A distant homologue of HemN1 from S. sojae can complement this defect in S. fredii and S. sp. III, but exhibited a more pleotropic role in symbiosis establishment. These findings highlighted the lineage-dependent optimization of symbiotic functions in different rhizobial species associated with the same host.

Medicago truncatula Zinc-Iron Permease6 provides zinc to rhizobia-infected nodule cells

Zinc is a micronutrient required for symbiotic nitrogen fixation. It has been proposed that in model legume Medicago truncatula, zinc is delivered by the root vasculature into the nodule and released in the infection/differentiation zone. There, transporters must introduce this element into rhizobia-infected cells to metallate the apoproteins that use zinc as a cofactor. MtZIP6 (Medtr4g083570) is an M. truncatula Zinc-Iron Permease (ZIP) that is expressed only in roots and nodules, with the highest expression levels in the infection/differentiation zone. Immunolocalization studies indicate that it is located in the plasma membrane of nodule rhizobia-infected cells. Down-regulating MtZIP6 expression levels with RNAi does not result in any strong phenotype when plants are fed mineral nitrogen. However, these plants displayed severe growth defects when they depended on nitrogen fixed by their nodules, losing of 80% of their nitrogenase activity. The reduction of this activity was likely an indirect effect of zinc being retained in the infection/differentiation zone and not reaching the cytosol of rhizobia-infected cells. These data are consistent with a model in which MtZIP6 would be responsible for zinc uptake by rhizobia-infected nodule cells in the infection/differentiation zone.

The Sinorhizobium fredii HH103 MucR1 Global Regulator Is Connected With the nod Regulon and Is Required for Efficient Symbiosis With Lotus burttii and Glycine max cv. Williams

Sinorhizobium fredii HH103 is a rhizobial strain showing a broad host range of nodulation. In addition to the induction of bacterial nodulation genes, transition from a free-living to a symbiotic state requires complex genetic expression changes with the participation of global regulators. We have analyzed the role of the zinc-finger transcriptional regulator MucR1 from S. fredii HH103 under both free-living conditions and symbiosis with two HH103 host plants, Glycine max and Lotus burttii. Inactivation of HH103 mucR1 led to a severe decrease in exopolysaccharide (EPS) biosynthesis but enhanced production of external cyclic glucans (CG). This mutant also showed increased cell aggregation capacity as well as a drastic reduction in nitrogen-fixation capacity with G. max and L. burttii. However, in these two legumes, the number of nodules induced by the mucR1 mutant was significantly increased and decreased, respectively, with respect to the wild-type strain, indicating that MucR1 can differently affect nodulation depending on the host plant. RNA-Seq analysis carried out in the absence and the presence of flavonoids showed that MucR1 controls the expression of hundreds of genes (including some related to EPS production and CG transport), some of them being related to the nod regulon.

Exopolysaccharide Production by Sinorhizobium fredii HH103 Is Repressed by Genistein in a NodD1-Dependent Manner

In the rhizobia-legume symbiotic interaction, bacterial surface polysaccharides, such as exopolysaccharide (EPS), lipopolysaccharide (LPS), K-antigen polysaccharide (KPS) or cyclic glucans (CG), appear to play crucial roles either acting as signals required for the progression of the interaction and/or preventing host defence mechanisms. The symbiotic significance of each of these polysaccharides varies depending on the specific rhizobia-legume couple. In this work we show that the production of exopolysaccharide by Sinorhizobium fredii HH103, but not by other S. fredii strains such as USDA257 or NGR234, is repressed by nod gene inducing flavonoids such as genistein and that this repression is dependent on the presence of a functional NodD1 protein. In agreement with the importance of EPS for bacterial biofilms, this reduced EPS production upon treatment with flavonoids correlates with decreased biofilm formation ability. By using quantitative RT-PCR analysis we show that expression of the exoY2 and exoK genes is repressed in late stationary cultures of S. fredii HH103 upon treatment with genistein. Results presented in this work show that in S. fredii HH103 EPS production is regulated just in the opposite way than other bacterial signals such as Nod factors and type 3 secreted effectors: it is repressed by flavonoids and NodD1 and enhanced by the nod repressor NolR. These results are in agreement with our previous observations showing that lack of EPS production by S. fredii HH103 is not only non-detrimental but even beneficial for symbiosis with soybean.

RNA-Seq Analysis of Differential Gene Expression Responding to Different Rhizobium Strains in Soybean (Glycine max) Roots

The root nodule symbiosis (RNS) between legume plants and rhizobia is the most efficient and productive source of nitrogen fixation, and has critical importance in agriculture and mesology. Soybean (Glycine max), one of the most important legume crops in the world, establishes a nitrogen-fixing symbiosis with different types of rhizobia, and the efficiency of symbiotic nitrogen fixation in soybean greatly depends on the symbiotic host-specificity. Although, it has been reported that rhizobia use surface polysaccharides, secretion proteins of the type-three secretion systems and nod factors to modulate host range, the host control of nodulation specificity remains poorly understood. In this report, the soybean roots of two symbiotic systems (Bradyrhizobium japonicum strain 113-2-soybean and Sinorhizobium fredii USDA205-soybean)with notable different nodulation phenotypes and the control were studied at five different post-inoculation time points (0.5, 7-24 h, 5, 16, and 21 day) by RNA-seq (Quantification). The results of qPCR analysis of 11 randomly-selected genes agreed with transcriptional profile data for 136 out of 165 (82.42%) data points and quality assessment showed that the sequencing library is of quality and reliable. Three comparisons (control vs. 113-2, control vs. USDA205 and USDA205 vs. 113-2) were made and the differentially expressed genes (DEGs) between them were analyzed. The number of DEGs at 16 days post-inoculation (dpi) was the highest in the three comparisons, and most of the DEGs in USDA205 vs. 113-2 were found at 16 dpi and 21 dpi. 44 go function terms in USDA205 vs. 113-2 were analyzed to evaluate the potential functions of the DEGs, and 10 important KEGG pathway enrichment terms were analyzed in the three comparisons. Some important genes induced in response to different strains (113-2 and USDA205) were identified and analyzed, and these genes primarily encoded soybean resistance proteins, NF-related proteins, nodulins and immunity defense proteins, as well as proteins involving flavonoids/flavone/flavonol biosynthesis and plant-pathogen interaction. Besides, 189 candidate genes are largely expressed in roots and\or nodules. The DEGs uncovered in this study provides molecular candidates for better understanding the mechanisms of symbiotic host-specificity and explaining the different symbiotic effects between soybean roots inoculated with different strains (113-2 and USDA205).