Heterologous Complementation Reveals a Specialized Activity for BacA in the Medicago-Sinorhizobium meliloti Symbiosis

Symbiotic Nodule Development and Efficiency in the Medicago truncatula Mtefd-1 Mutant Is Highly Dependent on Sinorhizobium Strains

Symbiotic nitrogen fixation (SNF) can play a key role in agroecosystems to reduce the negative impact of nitrogen fertilizers. Its efficiency is strongly affected by the combination of bacterial and plant genotypes, but the mechanisms responsible for the differences in the efficiency of rhizobium strains are not well documented. In Medicago truncatula, SNF has been mostly studied using model systems, such as M. truncatula A17 in interaction with Sinorhizobium meliloti Sm2011. Here we analyzed both the wild-type (wt) A17 and the Mtefd-1 mutant in interaction with five S. meliloti and two Sinorhizobium medicae strains. ETHYLENE RESPONSE FACTOR REQUIRED FOR NODULE DIFFERENTIATION (MtEFD) encodes a transcription factor, which contributes to the control of nodule number and differentiation in M. truncatula. We found that, in contrast to Sm2011, four strains induce functional (Fix+) nodules in Mtefd-1, although less efficient for SNF than in wt A17. In contrast, the Mtefd-1 hypernodulation phenotype is not strain-dependent. We compared the plant nodule transcriptomes in response to SmBL225C, a highly efficient strain with A17, versus Sm2011, in wt and Mtefd-1 backgrounds. This revealed faster nodule development with SmBL225C and early nodule senescence with Sm2011. These RNA sequencing analyses allowed us to identify candidate plant factors that could drive the differential nodule phenotype. In conclusion, this work shows the value of having a set of rhizobium strains to fully evaluate the biological importance of a plant symbiotic gene.

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.

Incompatibility and Interchangeability in Molecular Evolution

There is remarkable variation in the rate at which genetic incompatibilities in molecular interactions accumulate. In some cases, minor changes-even single-nucleotide substitutions-create major incompatibilities when hybridization forces new variants to function in a novel genetic background from an isolated population. In other cases, genes or even entire functional pathways can be horizontally transferred between anciently divergent evolutionary lineages that span the tree of life with little evidence of incompatibilities. In this review, we explore whether there are general principles that can explain why certain genes are prone to incompatibilities while others maintain interchangeability. We summarize evidence pointing to four genetic features that may contribute to greater resistance to functional replacement: (1) function in multisubunit enzyme complexes and protein-protein interactions, (2) sensitivity to changes in gene dosage, (3) rapid rate of sequence evolution, and (4) overall importance to cell viability, which creates sensitivity to small perturbations in molecular function. We discuss the relative levels of support for these different hypotheses and lay out future directions that may help explain the striking contrasts in patterns of incompatibility and interchangeability throughout the history of molecular evolution.

Reference nodule transcriptomes for Melilotus officinalis and Medicago sativa cv. Algonquin

Host/symbiont compatibility is a hallmark of the symbiotic nitrogen-fixing interaction between rhizobia and legumes, mediated in part by plant-produced nodule-specific cysteine-rich (NCR) peptides and the bacterial BacA membrane protein that can act as a NCR peptide transporter. In addition, the genetic and metabolic properties supporting symbiotic nitrogen fixation often differ between compatible partners, including those sharing a common partner, highlighting the need for multiple study systems. Here, we report high-quality nodule transcriptome assemblies for Medicago sativa cv. Algonquin and Melilotus officinalis, two legumes able to form compatible symbioses with Sinorhizobium meliloti. The compressed M. sativa and M. officinalis assemblies consisted of 79,978 and 64,593 contigs, respectively, of which 33,341 and 28,278 were assigned putative annotations, respectively. As expected, the two transcriptomes showed broad similarity at a global level. We were particularly interested in the NCR peptide profiles of these plants, as these peptides drive bacterial differentiation during the symbiosis. A total of 412 and 308 NCR peptides were predicted from the M. sativa and M. officinalis transcriptomes, respectively, with approximately 9% of the transcriptome of both species consisting of NCR transcripts. Notably, transcripts encoding highly cationic NCR peptides (isoelectric point > 9.5), which are known to have antimicrobial properties, were ∼2-fold more abundant in M. sativa than in M. officinalis, and ∼27-fold more abundant when considering only NCR peptides in the six-cysteine class. We hypothesize that the difference in abundance of highly cationic NCR peptides explains our previous observation that some rhizobial bacA alleles which can support symbiosis with M. officinalis are unable to support symbiosis with M. sativa.

The Genome of the Acid Soil-Adapted Strain Rhizobium favelukesii OR191 Encodes Determinants for Effective Symbiotic Interaction With Both an Inverted Repeat Lacking Clade and a Phaseoloid Legume Host

Although Medicago sativa forms highly effective symbioses with the comparatively acid-sensitive genus Ensifer, its introduction into acid soils appears to have selected for symbiotic interactions with acid-tolerant R. favelukesii strains. Rhizobium favelukesii has the unusual ability of being able to nodulate and fix nitrogen, albeit sub-optimally, not only with M. sativa but also with the promiscuous host Phaseolus vulgaris. Here we describe the genome of R. favelukesii OR191 and genomic features important for the symbiotic interaction with both of these hosts. The OR191 draft genome contained acid adaptation loci, including the highly acid-inducible lpiA/acvB operon and olsC, required for production of lysine- and ornithine-containing membrane lipids, respectively. The olsC gene was also present in other acid-tolerant Rhizobium strains but absent from the more acid-sensitive Ensifer microsymbionts. The OR191 symbiotic genes were in general more closely related to those found in Medicago microsymbionts. OR191 contained the nodA, nodEF, nodHPQ, and nodL genes for synthesis of polyunsaturated, sulfated and acetylated Nod factors that are important for symbiosis with Medicago, but contained a truncated nodG, which may decrease nodulation efficiency with M. sativa. OR191 contained an E. meliloti type BacA, which has been shown to specifically protect Ensifer microsymbionts from Medicago nodule-specific cysteine-rich peptides. The nitrogen fixation genes nifQWZS were present in OR191 and P. vulgaris microsymbionts but absent from E. meliloti-Medicago microsymbionts. The ability of OR191 to nodulate and fix nitrogen symbiotically with P. vulgaris indicates that this host has less stringent requirements for nodulation than M. sativa but may need rhizobial strains that possess nifQWZS for N₂-fixation to occur. OR191 possessed the exo genes required for the biosynthesis of succinoglycan, which is required for the Ensifer-Medicago symbiosis. However, ¹H-NMR spectra revealed that, in the conditions tested, OR191 exopolysaccharide did not contain a succinyl substituent but instead contained a 3-hydroxybutyrate moiety, which may affect its symbiotic performance with Medicago hosts. These findings provide a foundation for the genetic basis of nodulation requirements and symbiotic effectiveness with different hosts.

Varietas Delectat: Exploring Natural Variations in Nitrogen-Fixing Symbiosis Research

The nitrogen-fixing symbiosis between leguminous plants and soil bacteria collectively called rhizobia plays an important role in the global nitrogen cycle and is an essential component of sustainable agriculture. Genetic determinants directing the development and functioning of the interaction have been identified with the help of a very limited number of model plants and bacterial strains. Most of the information obtained from the study of model systems could be validated on crop plants and their partners. The investigation of soybean cultivars and different rhizobia, however, has revealed the existence of ineffective interactions between otherwise effective partners that resemble gene-for-gene interactions described for pathogenic systems. Since then, incompatible interactions between natural isolates of model plants, called ecotypes, and different bacterial partner strains have been reported. Moreover, diverse phenotypes of both bacterial mutants on different host plants and plant mutants with different bacterial strains have been described. Identification of the genetic factors behind the phenotypic differences did already and will reveal novel functions of known genes/proteins, the role of certain proteins in some interactions, and the fine regulation of the steps during nodule development.

Scent of a Symbiont: The Personalized Genetic Relationships of Rhizobium—Plant Interaction

Many molecular signals are exchanged between rhizobia and host legume plants, some of which are crucial for symbiosis to take place, while others are modifiers of the interaction, which have great importance in the competition with the soil microbiota and in the genotype-specific perception of host plants. Here, we review recent findings on strain-specific and host genotype-specific interactions between rhizobia and legumes, discussing the molecular actors (genes, gene products and metabolites) which play a role in the establishment of symbiosis, and highlighting the need for research including the other components of the soil (micro)biota, which could be crucial in developing rational-based strategies for bioinoculants and synthetic communities’ assemblage.

How rhizobia adapt to the nodule environment

Rhizobia are a phylogenetically diverse group of soil bacteria that engage in mutualistic interactions with legume plants. Although specifics of the symbioses differ between strains and plants, all symbioses ultimately result in the formation of specialized root nodule organs which host the nitrogen-fixing microsymbionts called bacteroids. Inside nodules, bacteroids encounter unique conditions that necessitate global reprogramming of physiological processes and rerouting of their metabolism. Decades of research have addressed these questions using genetics, omics approaches, and more recently computational modelling. Here we discuss the common adaptations of rhizobia to the nodule environment that define the core principles of bacteroid functioning. All bacteroids are growth-arrested and perform energy-intensive nitrogen fixation fueled by plant-provided C₄-dicarboxylates at nanomolar oxygen levels. At the same time, bacteroids are subject to host control and sanctioning that ultimately determine their fitness and have fundamental importance for the evolution of a stable mutualistic relationship.

Determinants of Host Range Specificity in Legume-Rhizobia Symbiosis

Leguminous plants possess the almost unique ability to enter symbiosis with soil-resident, nitrogen fixing bacteria called rhizobia. During this symbiosis, the bacteria physically colonize specialized organs on the roots of the host plant called nodules, where they reduce atmospheric nitrogen into forms that can be assimilated by the host plant and receive photosynthates in return. In order for nodule development to occur, there is extensive chemical cross-talk between both parties during the formative stages of the symbiosis. The vast majority of the legume family are capable of forming root nodules and typically rhizobia are only able to fix nitrogen within the context of this symbiotic association. However, many legume species only enter productive symbiosis with a few, or even single rhizobial species or strains, and vice-versa. Permitting symbiosis with only rhizobial strains that will be able to fix nitrogen with high efficiency is a crucial strategy for the host plant to prevent cheating by rhizobia. This selectivity is enforced at all stages of the symbiosis, with partner choice beginning during the initial communication between the plant and rhizobia. However, it can also be influenced even once nitrogen-fixing nodules have developed on the root. This review sets out current knowledge about the molecular mechanisms employed by both parties to influence host range during legume-rhizobia symbiosis.

Symbiotic and Nonsymbiotic Members of the Genus Ensifer (syn. Sinorhizobium) Are Separated into Two Clades Based on Comparative Genomics and High-Throughput Phenotyping

Rhizobium–legume symbioses serve as paradigmatic examples for the study of mutualism evolution. The genus Ensifer (syn. Sinorhizobium) contains diverse plant-associated bacteria, a subset of which can fix nitrogen in symbiosis with legumes. To gain insights into the evolution of symbiotic nitrogen fixation (SNF), and interkingdom mutualisms more generally, we performed extensive phenotypic, genomic, and phylogenetic analyses of the genus Ensifer. The data suggest that SNF likely emerged several times within the genus Ensifer through independent horizontal gene transfer events. Yet, the majority (105 of 106) of the Ensifer strains with the nodABC and nifHDK nodulation and nitrogen fixation genes were found within a single, monophyletic clade. Comparative genomics highlighted several differences between the “symbiotic” and “nonsymbiotic” clades, including divergences in their pangenome content. Additionally, strains of the symbiotic clade carried 325 fewer genes, on average, and appeared to have fewer rRNA operons than strains of the nonsymbiotic clade. Initial characterization of a subset of ten Ensifer strains identified several putative phenotypic differences between the clades. Tested strains of the nonsymbiotic clade could catabolize 25% more carbon sources, on average, than strains of the symbiotic clade, and they were better able to grow in LB medium and tolerate alkaline conditions. On the other hand, the tested strains of the symbiotic clade were better able to tolerate heat stress and acidic conditions. We suggest that these data support the division of the genus Ensifer into two main subgroups, as well as the hypothesis that pre-existing genetic features are required to facilitate the evolution of SNF in bacteria.

Effectiveness of nitrogen fixation in rhizobia

Biological nitrogen fixation in rhizobia occurs primarily in root or stem nodules and is induced by the bacteria present in legume plants. This symbiotic process has fascinated researchers for over a century, and the positive effects of legumes on soils and their food and feed value have been recognized for thousands of years. Symbiotic nitrogen fixation uses solar energy to reduce the inert N2 gas to ammonia at normal temperature and pressure, and is thus today, especially, important for sustainable food production. Increased productivity through improved effectiveness of the process is seen as a major research and development goal. The interaction between rhizobia and their legume hosts has thus been dissected at agronomic, plant physiological, microbiological and molecular levels to produce ample information about processes involved, but identification of major bottlenecks regarding efficiency of nitrogen fixation has proven to be complex. We review processes and results that contributed to the current understanding of this fascinating system, with focus on effectiveness of nitrogen fixation in rhizobia.

Genome-scale metabolic reconstruction of the symbiosis between a leguminous plant and a nitrogen-fixing bacterium

The mutualistic association between leguminous plants and endosymbiotic rhizobial bacteria is a paradigmatic example of a symbiosis driven by metabolic exchanges. Here, we report the reconstruction and modelling of a genome-scale metabolic network of Medicago truncatula (plant) nodulated by Sinorhizobium meliloti (bacterium). The reconstructed nodule tissue contains five spatially distinct developmental zones and encompasses the metabolism of both the plant and the bacterium. Flux balance analysis (FBA) suggests that the metabolic costs associated with symbiotic nitrogen fixation are primarily related to supporting nitrogenase activity, and increasing N2-fixation efficiency is associated with diminishing returns in terms of plant growth. Our analyses support that differentiating bacteroids have access to sugars as major carbon sources, ammonium is the main nitrogen export product of N2-fixing bacteria, and N2 fixation depends on proton transfer from the plant cytoplasm to the bacteria through acidification of the peribacteroid space. We expect that our model, called 'Virtual Nodule Environment' (ViNE), will contribute to a better understanding of the functioning of legume nodules, and may guide experimental studies and engineering of symbiotic nitrogen fixation.

Multidisciplinary approaches for studying rhizobium–legume symbioses

The rhizobium–legume symbiosis is a major source of fixed nitrogen (ammonia) in the biosphere. The potential for this process to increase agricultural yield while reducing the reliance on nitrogen-based fertilizers has generated interest in understanding and manipulating this process. For decades, rhizobium research has benefited from the use of leading techniques from a very broad set of fields, including population genetics, molecular genetics, genomics, and systems biology. In this review, we summarize many of the research strategies that have been employed in the study of rhizobia and the unique knowledge gained from these diverse tools, with a focus on genome- and systems-level approaches. We then describe ongoing synthetic biology approaches aimed at improving existing symbioses or engineering completely new symbiotic interactions. The review concludes with our perspective of the future directions and challenges of the field, with an emphasis on how the application of a multidisciplinary approach and the development of new methods will be necessary to ensure successful biotechnological manipulation of the symbiosis.

Robustness encoded across essential and accessory replicons of the ecologically versatile bacterium Sinorhizobium meliloti

Bacterial genome evolution is characterized by gains, losses, and rearrangements of functional genetic segments. The extent to which large-scale genomic alterations influence genotype-phenotype relationships has not been investigated in a high-throughput manner. In the symbiotic soil bacterium Sinorhizobium meliloti, the genome is composed of a chromosome and two large extrachromosomal replicons (pSymA and pSymB, which together constitute 45% of the genome). Massively parallel transposon insertion sequencing (Tn-seq) was employed to evaluate the contributions of chromosomal genes to growth fitness in both the presence and absence of these extrachromosomal replicons. Ten percent of chromosomal genes from diverse functional categories are shown to genetically interact with pSymA and pSymB. These results demonstrate the pervasive robustness provided by the extrachromosomal replicons, which is further supported by constraint-based metabolic modeling. A comprehensive picture of core S. meliloti metabolism was generated through a Tn-seq-guided in silico metabolic network reconstruction, producing a core network encompassing 726 genes. This integrated approach facilitated functional assignments for previously uncharacterized genes, while also revealing that Tn-seq alone missed over a quarter of wild-type metabolism. This work highlights the many functional dependencies and epistatic relationships that may arise between bacterial replicons and across a genome, while also demonstrating how Tn-seq and metabolic modeling can be used together to yield insights not obtainable by either method alone.

S. meliloti, which has traditionally facilitated ground-breaking insights into symbiotic communication, is also emerging as an excellent model for studying the evolution of functional relationships between bacterial chromosomes and anciently acquired accessory replicons. Multi-replicon genome architecture is present in ~ 10% of presently sequenced bacterial genomes. The S. meliloti genome is composed of three circular replicons, two of which are dispensable even though they encompass nearly half of the protein-coding genes in this organism. The construction of strains lacking these replicons has enabled a straightforward, genome-wide analysis of interactions between the chromosome and the non-essential replicons, revealing extensive functional cooperation between these genomic components. This analysis enabled a substantial refinement of a metabolic network model for S. meliloti. The integration of massively parallel genotype-phenotype screening with in silico metabolic reconstruction has enhanced our understanding of metabolic network structure as it relates to genome evolution in S. meliloti, and exemplifies an approach that may be productively applied to other taxa. The combined experimental and computational approach employed here further provides unique insights into the pervasive genetic interactions that may exist within large bacterial genomes.

Succinate Transport Is Not Essential for Symbiotic Nitrogen Fixation by Sinorhizobium meliloti or Rhizobium leguminosarum

Symbiotic nitrogen fixation (SNF) is an energetically expensive process performed by bacteria during endosymbiotic relationships with plants. The bacteria require the plant to provide a carbon source for the generation of reductant to power SNF. While C₄-dicarboxylates (succinate, fumarate, and malate) appear to be the primary, if not sole, carbon source provided to the bacteria, the contribution of each C₄-dicarboxylate is not known. We address this issue using genetic and systems-level analyses. Expression of a malate-specific transporter (MaeP) in Sinorhizobium meliloti Rm1021 dct mutants unable to transport C₄-dicarboxylates resulted in malate import rates of up to 30% that of the wild type. This was sufficient to support SNF with Medicago sativa, with acetylene reduction rates of up to 50% those of plants inoculated with wild-type S. melilotiRhizobium leguminosarum bv. viciae 3841 dct mutants unable to transport C₄-dicarboxylates but expressing the maeP transporter had strong symbiotic properties, with Pisum sativum plants inoculated with these strains appearing similar to plants inoculated with wild-type R. leguminosarum This was despite malate transport rates by the mutant bacteroids being 10% those of the wild type. An RNA-sequencing analysis of the combined P. sativum-R. leguminosarum nodule transcriptome was performed to identify systems-level adaptations in response to the inability of the bacteria to import succinate or fumarate. Few transcriptional changes, with no obvious pattern, were detected. Overall, these data illustrated that succinate and fumarate are not essential for SNF and that, at least in specific symbioses, l-malate is likely the primary C₄-dicarboxylate provided to the bacterium.IMPORTANCE Symbiotic nitrogen fixation (SNF) is an economically and ecologically important biological process that allows plants to grow in nitrogen-poor soils without the need to apply nitrogen-based fertilizers. Much research has been dedicated to this topic to understand this process and to eventually manipulate it for agricultural gains. The work presented in this article provides new insights into the metabolic integration of the plant and bacterial partners. It is shown that malate is the only carbon source that needs to be available to the bacterium to support SNF and that, at least in some symbioses, malate, and not other C₄-dicarboxylates, is likely the primary carbon provided to the bacterium. This work extends our knowledge of the minimal metabolic capabilities the bacterium requires to successfully perform SNF and may be useful in further studies aiming to optimize this process through synthetic biology approaches. The work describes an engineering approach to investigate a metabolic process that occurs between a eukaryotic host and its prokaryotic endosymbiont.

Integrated roles of BclA and DD-carboxypeptidase 1 in Bradyrhizobium differentiation within NCR-producing and NCR-lacking root nodules

Legumes harbor in their symbiotic nodule organs nitrogen fixing rhizobium bacteria called bacteroids. Some legumes produce Nodule-specific Cysteine-Rich (NCR) peptides in the nodule cells to control the intracellular bacterial population. NCR peptides have antimicrobial activity and drive bacteroids toward terminal differentiation. Other legumes do not produce NCR peptides and their bacteroids are not differentiated. Bradyrhizobia, infecting NCR-producing Aeschynomene plants, require the peptide uptake transporter BclA to cope with the NCR peptides as well as a specific peptidoglycan-modifying DD-carboxypeptidase, DD-CPase1. We show that Bradyrhizobium diazoefficiens strain USDA110 forms undifferentiated bacteroids in NCR-lacking soybean nodules. Unexpectedly, in Aeschynomene afraspera nodules the nitrogen fixing USDA110 bacteroids are hardly differentiated despite the fact that this host produces NCR peptides, suggesting that USDA110 is insensitive to the host peptide effectors and that nitrogen fixation can be uncoupled from differentiation. In agreement with the absence of bacteroid differentiation, USDA110 does not require its bclA gene for nitrogen fixing symbiosis with these two host plants. Furthermore, we show that the BclA and DD-CPase1 act independently in the NCR-induced morphological differentiation of bacteroids. Our results suggest that BclA is required to protect the rhizobia against the NCR stress but not to induce the terminal differentiation pathway.