Metabolic modelling reveals the specialization of secondary replicons for niche adaptation in Sinorhizobium meliloti

Reconstruction of the genome-scale metabolic network model of Sinorhizobium fredii CCBAU45436 for free-living and symbiotic states

Sinorhizobium fredii CCBAU45436 is an excellent rhizobium that plays an important role in agricultural production. However, there still needs more comprehensive understanding of the metabolic system of S. fredii CCBAU45436, which hinders its application in agriculture. Therefore, based on the first-generation metabolic model iCC541 we developed a new genome-scale metabolic model iAQY970, which contains 970 genes, 1,052 reactions, 942 metabolites and is scored 89% in the MEMOTE test. Cell growth phenotype predicted by iAQY970 is 81.7% consistent with the experimental data. The results of mapping the proteome data under free-living and symbiosis conditions to the model showed that the biomass production rate in the logarithmic phase was faster than that in the stable phase, and the nitrogen fixation efficiency of rhizobia parasitized in cultivated soybean was higher than that in wild-type soybean, which was consistent with the actual situation. In the symbiotic condition, there are 184 genes that would affect growth, of which 94 are essential; In the free-living condition, there are 143 genes that influence growth, of which 78 are essential. Among them, 86 of the 94 essential genes in the symbiotic condition were consistent with the prediction of iCC541, and 44 essential genes were confirmed by literature information; meanwhile, 30 genes were identified by DEG and 33 genes were identified by Geptop. In addition, we extracted four key nitrogen fixation modules from the model and predicted that sulfite reductase (EC 1.8.7.1) and nitrogenase (EC 1.18.6.1) as the target enzymes to enhance nitrogen fixation by MOMA, which provided a potential focus for strain optimization. Through the comprehensive metabolic model, we can better understand the metabolic capabilities of S. fredii CCBAU45436 and make full use of it in the future.

A genome-scale metabolic reconstruction of soybean and Bradyrhizobium diazoefficiens reveals the cost-benefit of nitrogen fixation

Nitrogen-fixing symbioses allow legumes to thrive in nitrogen-poor soils at the cost of diverting some photoassimilate to their microsymbionts. Effort is being made to bioengineer nitrogen fixation into nonleguminous crops. This requires a quantitative understanding of its energetic costs and the links between metabolic variations and symbiotic efficiency. A whole-plant metabolic model for soybean (Glycine max) with its associated microsymbiont Bradyrhizobium diazoefficiens was developed and applied to predict the cost-benefit of nitrogen fixation with varying soil nitrogen availability. The model predicted a nitrogen-fixation cost of c. 4.13 g C g-1 N, which when implemented into a crop scale model, translated to a grain yield reduction of 27% compared with a non-nodulating plant receiving its nitrogen from the soil. Considering the lower nitrogen content of cereals, the yield cost to a hypothetical N-fixing cereal is predicted to be less than half that of soybean. Soybean growth was predicted to be c. 5% greater when the nodule nitrogen export products were amides versus ureides. This is the first metabolic reconstruction in a tropical crop species that simulates the entire plant and nodule metabolism. Going forward, this model will serve as a tool to investigate carbon use efficiency and key mechanisms within N-fixing symbiosis in a tropical species forming determinate nodules.

Independent origins and evolution of the secondary replicons of the class Gammaproteobacteria

Multipartite genomes, consisting of more than one replicon, have been found in approximately 10 % of bacteria, many of which belong to the phylum Proteobacteria. Many aspects of their origin and evolution, and the possible advantages related to this type of genome structure, remain to be elucidated. Here, we performed a systematic analysis of the presence and distribution of multipartite genomes in the class Gammaproteobacteria, which includes several genera with diverse lifestyles. Within this class, multipartite genomes are mainly found in the order Alteromonadales (mostly in the genus Pseudoalteromonas) and in the family Vibrionaceae. Our data suggest that the emergence of secondary replicons in Gammaproteobacteria is rare and that they derive from plasmids. Despite their multiple origins, we highlighted the presence of evolutionary trends such as the inverse proportionality of the genome to chromosome size ratio, which appears to be a general feature of bacteria with multipartite genomes irrespective of taxonomic group. We also highlighted some functional trends. The core gene set of the secondary replicons is extremely small, probably limited to essential genes or genes that favour their maintenance in the genome, while the other genes are less conserved. This hypothesis agrees with the idea that the primary advantage of secondary replicons could be to facilitate gene acquisition through horizontal gene transfer, resulting in replicons enriched in genes associated with adaptation to different ecological niches. Indeed, secondary replicons are enriched both in genes that could promote adaptation to harsh environments, such as those involved in antibiotic, biocide and metal resistance, and in functional categories related to the exploitation of environmental resources (e.g. carbohydrates), which can complement chromosomal functions.

Control of the rhizobium-legume symbiosis by the plant nitrogen demand is tightly integrated at the whole plant level and requires inter-organ systemic signaling

Symbiotic nodules formed on legume roots with rhizobia fix atmospheric N₂. Bacteria reduce N₂ to NH₄ ⁺ that is assimilated into amino acids by the plant. In return, the plant provides photosynthates to fuel the symbiotic nitrogen fixation. Symbiosis is tightly adjusted to the whole plant nutritional demand and to the plant photosynthetic capacities, but regulatory circuits behind this control remain poorly understood. The use of split-root systems combined with biochemical, physiological, metabolomic, transcriptomic, and genetic approaches revealed that multiple pathways are acting in parallel. Systemic signaling mechanisms of the plant N demand are required for the control of nodule organogenesis, mature nodule functioning, and nodule senescence. N-satiety/N-deficit systemic signaling correlates with rapid variations of the nodules' sugar levels, tuning symbiosis by C resources allocation. These mechanisms are responsible for the adjustment of plant symbiotic capacities to the mineral N resources. On the one hand, if mineral N can satisfy the plant N demand, nodule formation is inhibited, and nodule senescence is activated. On the other hand, local conditions (abiotic stresses) may impair symbiotic activity resulting in plant N limitation. In these conditions, systemic signaling may compensate the N deficit by stimulating symbiotic root N foraging. In the past decade, several molecular components of the systemic signaling pathways controlling nodule formation have been identified, but a major challenge remains, that is, to understand their specificity as compared to the mechanisms of non-symbiotic plants that control root development and how they contribute to the whole plant phenotypes. Less is known about the control of mature nodule development and functioning by N and C nutritional status of the plant, but a hypothetical model involving the sucrose allocation to the nodule as a systemic signaling process, the oxidative pentose phosphate pathway, and the redox status as potential effectors of this signaling is emerging. This work highlights the importance of organism integration in plant biology.

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.

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.

New Insights into the Taxonomy of Bacteria in the Genomic Era and a Case Study with Rhizobia

Since early studies, the history of prokaryotes taxonomy has dealt with many changes driven by the development of new and more robust technologies. As a result, the number of new taxa descriptions is exponentially increasing, while an increasing number of others has been subject of reclassification, demanding from the taxonomists more effort to maintain an organized hierarchical system. However, expectations are that the taxonomy of prokaryotes will acquire a more stable status with the genomic era. Other analyses may continue to be necessary to determine microbial features, but the use of genomic data might be sufficient to provide reliable taxa delineation, helping taxonomy to reach the goal of correct classification and identification. Here we describe the evolution of prokaryotes' taxonomy until the genomic era, emphasizing bacteria and taking as an example the history of rhizobia taxonomy. This example was chosen because of the importance of the symbiotic nitrogen fixation of legumes with rhizobia to the nitrogen input to both natural ecosystems and agricultural crops. This case study reports the technological advances and the methodologies used to classify and identify bacterial species and indicates the actual rules required for an accurate description of new taxa.

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.

Genome-Scale Metabolic Modelling of Lifestyle Changes in Rhizobium leguminosarum

Biological nitrogen fixation in rhizobium-legume symbioses is of major importance for sustainable agricultural practices. To establish a mutualistic relationship with their plant host, rhizobia transition from free-living bacteria in soil to growth down infection threads inside plant roots and finally differentiate into nitrogen-fixing bacteroids. We reconstructed a genome-scale metabolic model for Rhizobium leguminosarum and integrated the model with transcriptome, proteome, metabolome, and gene essentiality data to investigate nutrient uptake and metabolic fluxes characteristic of these different lifestyles. Synthesis of leucine, polyphosphate, and AICAR is predicted to be important in the rhizosphere, while myo-inositol catabolism is active in undifferentiated nodule bacteria in agreement with experimental evidence. The model indicates that bacteroids utilize xylose and glycolate in addition to dicarboxylates, which could explain previously described gene expression patterns. Histidine is predicted to be actively synthesized in bacteroids, consistent with transcriptome and proteome data for several rhizobial species. These results provide the basis for targeted experimental investigation of metabolic processes specific to the different stages of the rhizobium-legume symbioses. IMPORTANCE Rhizobia are soil bacteria that induce nodule formation on plant roots and differentiate into nitrogen-fixing bacteroids. A detailed understanding of this complex symbiosis is essential for advancing ongoing efforts to engineer novel symbioses with cereal crops for sustainable agriculture. Here, we reconstruct and validate a genome-scale metabolic model for Rhizobium leguminosarum bv. viciae 3841. By integrating the model with various experimental data sets specific to different stages of symbiosis formation, we elucidate the metabolic characteristics of rhizosphere bacteria, undifferentiated bacteria inside root nodules, and nitrogen-fixing bacteroids. Our model predicts metabolic flux patterns for these three distinct lifestyles, thus providing a framework for the interpretation of genome-scale experimental data sets and identifying targets for future experimental studies.

Plant developmental stage drives the differentiation in ecological role of the maize microbiome

BACKGROUND

Plants live with diverse microbial communities which profoundly affect multiple facets of host performance, but if and how host development impacts the assembly, functions and microbial interactions of crop microbiomes are poorly understood. Here we examined both bacterial and fungal communities across soils, epiphytic and endophytic niches of leaf and root, and plastic leaf of fake plant (representing environment-originating microbes) at three developmental stages of maize at two contrasting sites, and further explored the potential function of phylloplane microbiomes based on metagenomics.

RESULTS

Our results suggested that plant developmental stage had a much stronger influence on the microbial diversity, composition and interkingdom networks in plant compartments than in soils, with the strongest effect in the phylloplane. Phylloplane microbiomes were co-shaped by both plant growth and seasonal environmental factors, with the air (represented by fake plants) as its important source. Further, we found that bacterial communities in plant compartments were more strongly driven by deterministic processes at the early stage but a similar pattern was for fungal communities at the late stage. Moreover, bacterial taxa played a more important role in microbial interkingdom network and crop yield prediction at the early stage, while fungal taxa did so at the late stage. Metagenomic analyses further indicated that phylloplane microbiomes possessed higher functional diversity at the early stage than the late stage, with functional genes related to nutrient provision enriched at the early stage and N assimilation and C degradation enriched at the late stage. Coincidently, more abundant beneficial bacterial taxa like Actinobacteria, Burkholderiaceae and Rhizobiaceae in plant microbiomes were observed at the early stage, but more saprophytic fungi at the late stage.

CONCLUSIONS

Our results suggest that host developmental stage profoundly influences plant microbiome assembly and functions, and the bacterial and fungal microbiomes take a differentiated ecological role at different stages of plant development. This study provides empirical evidence for host exerting strong effect on plant microbiomes by deterministic selection during plant growth and development. These findings have implications for the development of future tools to manipulate microbiome for sustainable increase in primary productivity. Video Abstract.

Metabolic control of nitrogen fixation in rhizobium-legume symbioses

Rhizobia induce nodule formation on legume roots and differentiate into bacteroids, which catabolize plant-derived dicarboxylates to reduce atmospheric N₂ into ammonia. Despite the agricultural importance of this symbiosis, the mechanisms that govern carbon and nitrogen allocation in bacteroids and promote ammonia secretion to the plant are largely unknown. Using a metabolic model derived from genome-scale datasets, we show that carbon polymer synthesis and alanine secretion by bacteroids facilitate redox balance in microaerobic nodules. Catabolism of dicarboxylates induces not only a higher oxygen demand but also a higher NADH/NAD⁺ ratio than sugars. Modeling and ¹³C metabolic flux analysis indicate that oxygen limitation restricts the decarboxylating arm of the tricarboxylic acid cycle, which limits ammonia assimilation into glutamate. By tightly controlling oxygen supply and providing dicarboxylates as the energy and electron source donors for N₂ fixation, legumes promote ammonia secretion by bacteroids. This is a defining feature of rhizobium-legume symbioses.

Rhizobial Chemoattractants, the Taste and Preferences of Legume Symbionts

The development of host-microbe interactions between legumes and their cognate rhizobia requires localization of the bacteria to productive sites of initiation on the plant roots. This end is achieved by the motility apparatus that propels the bacterium and the chemotaxis system that guides it. Motility and chemotaxis aid rhizobia in their competitiveness for space, resources, and nodulation opportunities. Here, we examine studies on chemotaxis of three major model rhizobia, namely Sinorhizobium meliloti, Rhizobium leguminosarum, and Bradyrhizobium japonicum, cataloging their range of attractant molecules and correlating this in the context of root and seed exudate compositions. Current research areas will be summarized, gaps in knowledge discussed, and future directions described.

Characterization of the sorbitol dehydrogenase SmoS from Sinorhizobium meliloti 1021

Sinorhizobium meliloti 1021 is a Gram-negative alphaproteobacterium with a robust capacity for carbohydrate metabolism. The enzymes that facilitate these reactions assist in the survival of the bacterium across a range of environmental niches, and they may also be suitable for use in industrial processes. SmoS is a dehydrogenase that catalyzes the oxidation of the commonly occurring sugar alcohols sorbitol and galactitol to fructose and tagatose, respectively, using NAD⁺ as a cofactor. The main objective of this study was to evaluate SmoS using biochemical techniques. The nucleotide sequence was codon-optimized for heterologous expression in Escherichia coli BL21 (DE3) Gold cells and the protein was subsequently overexpressed and purified. Size-exclusion chromatography and X-ray diffraction experiments suggest that SmoS is a tetramer. SmoS was crystallized, and crystals obtained in the absence of substrate diffracted to 2.1 Å resolution and those of a complex with sorbitol diffracted to 2.0 Å resolution. SmoS was characterized kinetically and shown to have a preference for sorbitol despite having a higher affinity for galactitol. Computational ligand-docking experiments suggest that tagatose binds the protein in a more energetically favourable complex than fructose, which is retained in the active site over a longer time frame following oxidation and reduces the rate of the reaction. These results supplement the inventory of biomolecules with potential for industrial applications and enhance the understanding of metabolism in the model organism S. meliloti.

Reconstruction and analysis of a genome-scale metabolic model for Agrobacterium tumefaciens

The plant pathogen Agrobacterium tumefaciens causes crown gall disease and is a widely used tool for generating transgenic plants owing to its virulence. The pathogenic process involves a shift from an independent to a living form within a host plant. However, comprehensive analyses of metabolites, genes, and reactions contributing to this complex process are lacking. To gain new insights about the pathogenicity from the viewpoints of physiology and cellular metabolism, a genome-scale metabolic model (GSMM) was reconstructed for A. tumefaciens. The model, referred to as iNX1344, contained 1,344 genes, 1,441 reactions, and 1,106 metabolites. It was validated by analyses of in silico cell growth on 39 unique carbon or nitrogen sources and the flux distribution of carbon metabolism. A. tumefaciens metabolic characteristics under three ecological niches were modelled. A high capacity to access and metabolize nutrients is more important for rhizosphere colonization than in the soil, and substantial metabolic changes were detected during the shift from the rhizosphere to tumour environments. Furthermore, by integrating transcriptome data for tumour conditions, significant alterations in central metabolic pathways and secondary metabolite metabolism were identified. Overall, the GSMM and constraint-based analysis could decode the physiological and metabolic features of A. tumefaciens as well as interspecific interactions with hosts, thereby improving our understanding of host adaptation and infection mechanisms.

Pervasive RNA Regulation of Metabolism Enhances the Root Colonization Ability of Nitrogen-Fixing Symbiotic α-Rhizobia

The rhizosphere and rhizoplane are nutrient-rich but selective environments for the root microbiome. Here, we deciphered a posttranscriptional network regulated by the homologous trans-small RNAs (sRNAs) AbcR1 and AbcR2, which rewire the metabolism of the nitrogen-fixing α-rhizobium Sinorhizobium meliloti during preinfection stages of symbiosis with its legume host alfalfa. The LysR-type regulator LsrB, which transduces the cell redox state, is indispensable for AbcR1 expression in actively dividing bacteria, whereas the stress-induced transcription of AbcR2 depends on the alternative σ factor RpoH1. MS2 affinity purification coupled with RNA sequencing unveiled exceptionally large and overlapping AbcR1/2 mRNA interactomes, jointly representing ⁓6% of the S. meliloti protein-coding genes. Most mRNAs encode transport/metabolic proteins whose translation is silenced by base pairing to two distinct anti-Shine Dalgarno motifs that function independently in both sRNAs. A metabolic model-aided analysis of the targetomes predicted changes in AbcR1/2 expression driven by shifts in carbon/nitrogen sources, which were confirmed experimentally. Low AbcR1/2 levels in some defined media anticipated overexpression growth phenotypes linked to the silencing of specific mRNAs. As a proof of principle, we confirmed AbcR1/2-mediated downregulation of the l-amino acid AapQ permease. AbcR1/2 interactomes are well represented in rhizosphere-related S. meliloti transcriptomic signatures. Remarkably, a lack of AbcR1 specifically compromised the ability of S. meliloti to colonize the root rhizoplane. The AbcR1 regulon likely ranks the utilization of available substrates to optimize metabolism, thus conferring on S. meliloti an advantage for efficient rhizosphere/rhizoplane colonization. AbcR1 regulation is predicted to be conserved in related α-rhizobia, which opens unprecedented possibilities for engineering highly competitive biofertilizers. IMPORTANCE Nitrogen-fixing root nodule symbioses between rhizobia and legume plants provide more than half of the combined nitrogen incorporated annually into terrestrial ecosystems, rendering plant growth independent of environmentally unfriendly chemical fertilizers. The success of symbiosis depends primarily on the capacity of rhizobia to establish competitive populations in soil and rhizosphere environments. Here, we provide insights into the regulation and architecture of an extensive RNA posttranscriptional network that fine-tunes the metabolism of the alfalfa symbiont S. meliloti, thereby enhancing the ability of this beneficial bacterium to colonize nutrient-rich but extremely selective niches, such as the rhizosphere of its host plant. This pervasive RNA regulation of metabolism is a major adaptive mechanism, predicted to operate in diverse rhizobial species. Because RNA regulation relies on modifiable base-pairing interactions, our findings open unexplored avenues for engineering the legumes rhizobiome within sustainable agricultural practices.

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.

The ModelSEED Biochemistry Database for the integration of metabolic annotations and the reconstruction, comparison and analysis of metabolic models for plants, fungi and microbes

For over 10 years, ModelSEED has been a primary resource for the construction of draft genome-scale metabolic models based on annotated microbial or plant genomes. Now being released, the biochemistry database serves as the foundation of biochemical data underlying ModelSEED and KBase. The biochemistry database embodies several properties that, taken together, distinguish it from other published biochemistry resources by: (i) including compartmentalization, transport reactions, charged molecules and proton balancing on reactions; (ii) being extensible by the user community, with all data stored in GitHub; and (iii) design as a biochemical 'Rosetta Stone' to facilitate comparison and integration of annotations from many different tools and databases. The database was constructed by combining chemical data from many resources, applying standard transformations, identifying redundancies and computing thermodynamic properties. The ModelSEED biochemistry is continually tested using flux balance analysis to ensure the biochemical network is modeling-ready and capable of simulating diverse phenotypes. Ontologies can be designed to aid in comparing and reconciling metabolic reconstructions that differ in how they represent various metabolic pathways. ModelSEED now includes 33,978 compounds and 36,645 reactions, available as a set of extensible files on GitHub, and available to search at https://modelseed.org/biochem and KBase.

The ModelSEED Biochemistry Database for the integration of metabolic annotations and the reconstruction, comparison and analysis of metabolic models for plants, fungi and microbes

For over ten years, ModelSEED has been a primary resource for the construction of draft genome-scale metabolic models based on annotated microbial or plant genomes. Now being released, the biochemistry database serves as the foundation of biochemical data underlying ModelSEED and KBase. The biochemistry database embodies several properties that, taken together, distinguish it from other published biochemistry resources by: (i) including compartmentalization, transport reactions, charged molecules and proton balancing on reactions;; (ii) being extensible by the user community, with all data stored in GitHub; and (iii) design as a biochemical “Rosetta Stone” to facilitate comparison and integration of annotations from many different tools and databases. The database was constructed by combining chemical data from many resources, applying standard transformations, identifying redundancies, and computing thermodynamic properties. The ModelSEED biochemistry is continually tested using flux balance analysis to ensure the biochemical network is modeling-ready and capable of simulating diverse phenotypes. Ontologies can be designed to aid in comparing and reconciling metabolic reconstructions that differ in how they represent various metabolic pathways. ModelSEED now includes 33,978 compounds and 36,645 reactions, available as a set of extensible files on GitHub, and available to search at https://modelseed.org and KBase.

Proteome Analysis Reveals a Significant Host-Specific Response in Rhizobium leguminosarum bv. viciae Endosymbiotic Cells

The Rhizobium-legume symbiosis is a beneficial interaction in which the bacterium converts atmospheric nitrogen into ammonia and delivers it to the plant in exchange for carbon compounds. This symbiosis implies the adaptation of bacteria to live inside host plant cells. In this work, we apply RP-LC-MS/MS and isobaric tags as relative and absolute quantitation techniques to study the proteomic profile of endosymbiotic cells (bacteroids) induced by Rhizobium leguminosarum bv viciae strain UPM791 in legume nodules. Nitrogenase subunits, tricarboxylic acid cycle enzymes, and stress-response proteins are among the most abundant from over 1000 rhizobial proteins identified in pea (Pisum sativum) bacteroids. Comparative analysis of bacteroids induced in pea and in lentil (Lens culinaris) nodules revealed the existence of a significant host-specific differential response affecting dozens of bacterial proteins, including stress-related proteins, transcriptional regulators, and proteins involved in the carbon and nitrogen metabolisms. A mutant affected in one of these proteins, homologous to a GntR-like transcriptional regulator, showed a symbiotic performance significantly impaired in symbiosis with pea but not with lentil plants. Analysis of the proteomes of bacteroids isolated from both hosts also revealed the presence of different sets of plant-derived nodule-specific cysteine-rich peptides, indicating that the endosymbiotic bacteria find a host-specific cocktail of chemical stressors inside the nodule. By studying variations of the bacterial response to different plant cell environments, we will be able to identify specific limitations imposed by the host that might give us clues for the improvement of rhizobial performance.

Plasmids Related to the Symbiotic Nitrogen Fixation Are Not Only Cooperated Functionally but Also May Have Evolved over a Time Span in Family Rhizobiaceae

Rhizobia are soil bacteria capable of forming symbiotic nitrogen-fixing nodules associated with leguminous plants. In fast-growing legume-nodulating rhizobia, such as the species in the family Rhizobiaceae, the symbiotic plasmid is the main genetic basis for nitrogen-fixing symbiosis, and is susceptible to horizontal gene transfer. To further understand the symbioses evolution in Rhizobiaceae, we analyzed the pan-genome of this family based on 92 genomes of type/reference strains and reconstructed its phylogeny using a phylogenomics approach. Intriguingly, although the genetic expansion that occurred in chromosomal regions was the main reason for the high proportion of low-frequency flexible gene families in the pan-genome, gene gain events associated with accessory plasmids introduced more genes into the genomes of nitrogen-fixing species. For symbiotic plasmids, although horizontal gene transfer frequently occurred, transfer may be impeded by, such as, the host’s physical isolation and soil conditions, even among phylogenetically close species. During coevolution with leguminous hosts, the plasmid system, including accessory and symbiotic plasmids, may have evolved over a time span, and provided rhizobial species with the ability to adapt to various environmental conditions and helped them achieve nitrogen fixation. These findings provide new insights into the phylogeny of Rhizobiaceae and advance our understanding of the evolution of symbiotic nitrogen fixation.

The Impacts of Domestication and Agricultural Practices on Legume Nutrient Acquisition Through Symbiosis With Rhizobia and Arbuscular Mycorrhizal Fungi

Legumes are unique among plants as they can obtain nitrogen through symbiosis with nitrogen-fixing rhizobia that form root nodules in the host plants. Therefore they are valuable crops for sustainable agriculture. Increasing nitrogen fixation efficiency is not only important for achieving better plant growth and yield, but it is also crucial for reducing the use of nitrogen fertilizer. Arbuscular mycorrhizal fungi (AMF) are another group of important beneficial microorganisms that form symbiotic relationships with legumes. AMF can promote host plant growth by providing mineral nutrients and improving the soil ecosystem. The trilateral legume-rhizobia-AMF symbiotic relationships also enhance plant development and tolerance against biotic and abiotic stresses. It is known that domestication and agricultural activities have led to the reduced genetic diversity of cultivated germplasms and higher sensitivity to nutrient deficiencies in crop plants, but how domestication has impacted the capability of legumes to establish beneficial associations with rhizospheric microbes (including rhizobia and fungi) is not well-studied. In this review, we will discuss the impacts of domestication and agricultural practices on the interactions between legumes and soil microbes, focusing on the effects on AMF and rhizobial symbioses and hence nutrient acquisition by host legumes. In addition, we will summarize the genes involved in legume-microbe interactions and studies that have contributed to a better understanding of legume symbiotic associations using metabolic modeling.

Soil origin and plant genotype structure distinct microbiome compartments in the model legume Medicago truncatula

BACKGROUND

Understanding the genetic and environmental factors that structure plant microbiomes is necessary for leveraging these interactions to address critical needs in agriculture, conservation, and sustainability. Legumes, which form root nodule symbioses with nitrogen-fixing rhizobia, have served as model plants for understanding the genetics and evolution of beneficial plant-microbe interactions for decades, and thus have added value as models of plant-microbiome interactions. Here we use a common garden experiment with 16S rRNA gene amplicon and shotgun metagenomic sequencing to study the drivers of microbiome diversity and composition in three genotypes of the model legume Medicago truncatula grown in two native soil communities.

RESULTS

Bacterial diversity decreased between external (rhizosphere) and internal plant compartments (root endosphere, nodule endosphere, and leaf endosphere). Community composition was shaped by strong compartment × soil origin and compartment × plant genotype interactions, driven by significant soil origin effects in the rhizosphere and significant plant genotype effects in the root endosphere. Nevertheless, all compartments were dominated by Ensifer, the genus of rhizobia that forms root nodule symbiosis with M. truncatula, and additional shotgun metagenomic sequencing suggests that the nodulating Ensifer were not genetically distinguishable from those elsewhere in the plant. We also identify a handful of OTUs that are common in nodule tissues, which are likely colonized from the root endosphere.

CONCLUSIONS

Our results demonstrate strong host filtering effects, with rhizospheres driven by soil origin and internal plant compartments driven by host genetics, and identify several key nodule-inhabiting taxa that coexist with rhizobia in the native range. Our results set the stage for future functional genetic experiments aimed at expanding our pairwise understanding of legume-rhizobium symbiosis toward a more mechanistic understanding of plant microbiomes. Video Abstract.

The Effects of Domestication on Secondary Metabolite Composition in Legumes

Legumes are rich in secondary metabolites, such as polyphenols, alkaloids, and saponins, which are important defense compounds to protect the plant against herbivores and pathogens, and act as signaling molecules between the plant and its biotic environment. Legume-sourced secondary metabolites are well known for their potential benefits to human health as pharmaceuticals and nutraceuticals. During domestication, the color, smell, and taste of crop plants have been the focus of artificial selection by breeders. Since these agronomic traits are regulated by secondary metabolites, the basis behind the genomic evolution was the selection of the secondary metabolite composition. In this review, we will discuss the classification, occurrence, and health benefits of secondary metabolites in legumes. The differences in their profiles between wild legumes and their cultivated counterparts will be investigated to trace the possible effects of domestication on secondary metabolite compositions, and the advantages and drawbacks of such modifications. The changes in secondary metabolite contents will also be discussed at the genetic level to examine the genes responsible for determining the secondary metabolite composition that might have been lost due to domestication. Understanding these genes would enable breeding programs and metabolic engineering to produce legume varieties with favorable secondary metabolite profiles for facilitating adaptations to a changing climate, promoting beneficial interactions with biotic factors, and enhancing health-beneficial secondary metabolite contents for human consumption.

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.

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.

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.

Deciphering the Symbiotic Plant Microbiome: Translating the Most Recent Discoveries on Rhizobia for the Improvement of Agricultural Practices in Metal-Contaminated and High Saline Lands

Rhizosphere and plant-associated microorganisms have been intensely studied for their beneficial effects on plant growth and health. These mainly include nitrogen-fixing bacteria (NFB) and plant-growth promoting rhizobacteria (PGPR). This beneficial fraction is involved in major functions such as plant nutrition and plant resistance to biotic and abiotic stresses, which include water deficiency and heavy-metal contamination. Consequently, crop yield emerges as the net result of the interactions between the plant genome and its associated microbiome. Here, we provide a review covering recent studies on PGP rhizobia as effective inoculants for agricultural practices in harsh soil, and we propose models for inoculant combinations and genomic manipulation strategies to improve crop yield.

Exploring the evolutionary dynamics of Rhizobium plasmids through bipartite network analysis

The genus Rhizobium usually has a multipartite genome architecture with a chromosome and several plasmids, making these bacteria a perfect candidate for plasmid biology studies. As there are no universally shared genes among typical plasmids, network analyses can complement traditional phylogenetics in a broad-scale study of plasmid evolution. Here, we present an exhaustive analysis of 216 plasmids from 49 complete genomes of Rhizobium by constructing a bipartite network that consists of two classes of nodes, the plasmids and homologous protein families that connect them. Dissection of the network using a hierarchical clustering strategy reveals extensive variety, with 34 homologous plasmid clusters. Four large clusters including one cluster of symbiotic plasmids and two clusters of chromids carrying some truly essential genes are widely distributed among Rhizobium. In contrast, the other clusters are quite small and rare. Symbiotic clusters and rare accessory clusters are exogenetic and do not appear to have co-evolved with the common accessory clusters; the latter ones have a large coding potential and functional complementarity for different lifestyles in Rhizobium. The bipartite network also provides preliminary evidence of Rhizobium plasmid variation and formation including genetic exchange, plasmid fusion and fission, exogenetic plasmid transfer, host plant selection, and environmental adaptation.

Metabolic Modeling of Pectobacterium parmentieri SCC3193 Provides Insights into Metabolic Pathways of Plant Pathogenic Bacteria

Understanding plant⁻microbe interactions is crucial for improving plants' productivity and protection. Constraint-based metabolic modeling is one of the possible ways to investigate the bacterial adaptation to different ecological niches and may give insights into the metabolic versatility of plant pathogenic bacteria. We reconstructed a raw metabolic model of the emerging plant pathogenic bacterium Pectobacterium parmentieri SCC3193 with the use of KBase. The model was curated by using inParanoind and phenotypic data generated with the use of the OmniLog system. Metabolic modeling was performed through COBRApy Toolbox v. 0.10.1. The curated metabolic model of P. parmentieri SCC3193 is highly reliable, as in silico obtained results overlapped up to 91% with experimental data on carbon utilization phenotypes. By mean of flux balance analysis (FBA), we predicted the metabolic adaptation of P. parmentieri SCC3193 to two different ecological niches, relevant for the persistence and plant colonization by this bacterium: soil and the rhizosphere. We performed in silico gene deletions to predict the set of essential core genes for this bacterium to grow in such environments. We anticipate that our metabolic model will be a valuable element for defining a set of metabolic targets to control infection and spreading of this plant pathogen.

Metabolic niches in the rhizosphere microbiome: new tools and approaches to analyse metabolic mechanisms of plant-microbe nutrient exchange

Plants nourish rhizospheric microbes via provision of carbon substrates, and the composition of the microbiome is strongly influenced by metabolic phenomena such as niche differentiation, competitive exclusion, and cross-feeding. Despite intensive investigations of the taxonomic structure in root microbiomes, there is relatively little biochemical knowledge of the metabolic niches occupied by microbial strains in the rhizosphere. Here, we review new tools and approaches that are boosting our knowledge of the metabolic mechanisms that shape the composition of the root microbiome. New studies have elucidated biochemical pathways that mediate root colonisation and pathogen suppression, and synthetic communities are emerging as a powerful tool to understand microbe-microbe interactions. Knowledge of root exudate composition is being advanced by new metabolomics methodologies, which have highlighted that specific exudate components can inhibit pathogen growth, and that certain metabolites can recruit mutualistic strains according to substrate uptake preferences. Microbial genomics is rapidly advancing, with large collections of isolated rhizosphere strains and mutant libraries giving new insights into the metabolic mechanisms of root colonisation. Exometabolomics is emerging as a powerful methodology for directly observing microbial uptake of root metabolites, and also for profiling microbial cross-feeding. Integrative studies using these resources should enable rapid advances, particularly when applied to standardised experimental set-ups and model synthetic communities.

Genomic sequencing, genome-scale metabolic network reconstruction, and in silico flux analysis of the grape endophytic fungus Alternaria sp. MG1

BACKGROUND

Alternaria sp. MG1, an endophytic fungus isolated from grape, is a native producer of resveratrol, which has important application potential. However, the metabolic characteristics and physiological behavior of MG1 still remains mostly unraveled. In addition, the resveratrol production of the strain is low. Thus, the whole-genome sequencing is highly required for elucidating the resveratrol biosynthesis pathway. Furthermore, the metabolic network model of MG1 was constructed to provide a computational guided approach for improving the yield of resveratrol.

RESULTS

Firstly, a draft genomic sequence of MG1 was generated with a size of 34.7 Mbp and a GC content of 50.96%. Genome annotation indicated that MG1 possessed complete biosynthesis pathways for stilbenoids, flavonoids, and lignins. Eight secondary metabolites involved in these pathways were detected by GC-MS analysis, confirming the metabolic diversity of MG1. Furthermore, the first genome-scale metabolic network of Alternaria sp. MG1 (named iYL1539) was reconstructed, accounting for 1539 genes, 2231 metabolites, and 2255 reactions. The model was validated qualitatively and quantitatively by comparing the in silico simulation with experimental data, and the results showed a high consistency. In iYL1539, 56 genes were identified as growth essential in rich medium. According to constraint-based analysis, the importance of cofactors for the resveratrol biosynthesis was successfully demonstrated. Ethanol addition was predicted in silico to be an effective method to improve resveratrol production by strengthening acetyl-CoA synthesis and pentose phosphate pathway, and was verified experimentally with a 26.31% increase of resveratrol. Finally, 6 genes were identified as potential targets for resveratrol over-production by the recently developed methodology. The target-genes were validated using salicylic acid as elicitor, leading to an increase of resveratrol yield by 33.32% and the expression of gene 4CL and CHS by 1.8- and 1.6-fold, respectively.

CONCLUSIONS

This study details the diverse capability and key genes of Alternaria sp. MG1 to produce multiple secondary metabolites. The first model of the species Alternaria was constructed, providing an overall understanding of the physiological behavior and metabolic characteristics of MG1. The model is a highly useful tool for enhancing productivity by rational design of the metabolic pathway for resveratrol and other secondary metabolites.

Tn-Core: A Toolbox for Integrating Tn-seq Gene Essentiality Data and Constraint-Based Metabolic Modeling

The design of synthetic cells requires a detailed understanding of the relevance of genes and gene networks underlying complex cellular phenotypes. Transposon-sequencing (Tn-seq) and constraint-based metabolic modeling can be used to probe the core genetic and metabolic networks underlying a biological process. Integrating these highly complementary experimental and in silico approaches has the potential to yield a highly comprehensive understanding of the core networks of a cell. Specifically, it can facilitate the interpretation of Tn-seq data sets and identify gaps in the data that could hinder the engineering of the cellular system, while also providing refined models for the accurate predictions of cellular metabolism. Here, we present Tn-Core, the first easy-to-use computational pipeline specifically designed for integrating Tn-seq data with metabolic modeling, prepared for use by both experimental and computational biologists. Tn-Core is a MATLAB toolbox that contains several custom functions, and it is built upon existing functions within the COBRA Toolbox and the TIGER Toolbox. Tn-Core takes as input a genome-scale metabolic model, Tn-seq data, and optionally RNA-seq data, and returns: (i) a context-specific core metabolic model; (ii) an evaluation of redundancies within core metabolic pathways, and optionally (iii) a refined genome-scale metabolic model. A simple, user-friendly workflow, requiring limited knowledge of metabolic modeling, is provided that allows users to run the analyses and export the data as easy-to-explore files of value to both experimental and computational biologists. We demonstrate the utility of Tn-Core using Sinorhizobium meliloti, Pseudomonas aeruginosa, and Rhodobacter sphaeroides genome-scale metabolic reconstructions as case studies.

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.

Ancestral Genome Estimation Reveals the History of Ecological Diversification in Agrobacterium

Horizontal gene transfer (HGT) is considered as a major source of innovation in bacteria, and as such is expected to drive adaptation to new ecological niches. However, among the many genes acquired through HGT along the diversification history of genomes, only a fraction may have actively contributed to sustained ecological adaptation. We used a phylogenetic approach accounting for the transfer of genes (or groups of genes) to estimate the history of genomes in Agrobacterium biovar 1, a diverse group of soil and plant-dwelling bacterial species. We identified clade-specific blocks of cotransferred genes encoding coherent biochemical pathways that may have contributed to the evolutionary success of key Agrobacterium clades. This pattern of gene coevolution rejects a neutral model of transfer, in which neighboring genes would be transferred independently of their function and rather suggests purifying selection on collectively coded acquired pathways. The acquisition of these synapomorphic blocks of cofunctioning genes probably drove the ecological diversification of Agrobacterium and defined features of ancestral ecological niches, which consistently hint at a strong selective role of host plant rhizospheres.

Harnessing Rhizobia to Improve Heavy-Metal Phytoremediation by Legumes

Rhizobia are bacteria that can form symbiotic associations with plants of the Fabaceae family, during which they reduce atmospheric di-nitrogen to ammonia. The symbiosis between rhizobia and leguminous plants is a fundamental contributor to nitrogen cycling in natural and agricultural ecosystems. Rhizobial microsymbionts are a major reason why legumes can colonize marginal lands and nitrogen-deficient soils. Several leguminous species have been found in metal-contaminated areas, and they often harbor metal-tolerant rhizobia. In recent years, there have been numerous efforts and discoveries related to the genetic determinants of metal resistance by rhizobia, and on the effectiveness of such rhizobia to increase the metal tolerance of host plants. Here, we review the main findings on the metal resistance of rhizobia: the physiological role, evolution, and genetic determinants, and the potential to use native and genetically-manipulated rhizobia as inoculants for legumes in phytoremediation practices.

Creation and Characterization of a Genomically Hybrid Strain in the Nitrogen-Fixing Symbiotic Bacterium Sinorhizobium meliloti

Many bacteria, often associated with eukaryotic hosts and of relevance for biotechnological applications, harbor a multipartite genome composed of more than one replicon. Biotechnologically relevant phenotypes are often encoded by genes residing on the secondary replicons. A synthetic biology approach to developing enhanced strains for biotechnological purposes could therefore involve merging pieces or entire replicons from multiple strains into a single genome. Here we report the creation of a genomic hybrid strain in a model multipartite genome species, the plant-symbiotic bacterium Sinorhizobium meliloti. We term this strain as cis-hybrid, since it is produced by genomic material coming from the same species' pangenome. In particular, we moved the secondary replicon pSymA (accounting for nearly 20% of total genome content) from a donor S. meliloti strain to an acceptor strain. The cis-hybrid strain was screened for a panel of complex phenotypes (carbon/nitrogen utilization phenotypes, intra- and extracellular metabolomes, symbiosis, and various microbiological tests). Additionally, metabolic network reconstruction and constraint-based modeling were employed for in silico prediction of metabolic flux reorganization. Phenotypes of the cis-hybrid strain were in good agreement with those of both parental strains. Interestingly, the symbiotic phenotype showed a marked cultivar-specific improvement with the cis-hybrid strains compared to both parental strains. These results provide a proof-of-principle for the feasibility of genome-wide replicon-based remodelling of bacterial strains for improved biotechnological applications in precision agriculture.

Genomic and Biotechnological Characterization of the Heavy-Metal Resistant, Arsenic-Oxidizing Bacterium Ensifer sp. M14

Ensifer (Sinorhizobium) sp. M14 is an efficient arsenic-oxidizing bacterium (AOB) that displays high resistance to numerous metals and various stressors. Here, we report the draft genome sequence and genome-guided characterization of Ensifer sp. M14, and we describe a pilot-scale installation applying the M14 strain for remediation of arsenic-contaminated waters. The M14 genome contains 6874 protein coding sequences, including hundreds not found in related strains. Nearly all unique genes that are associated with metal resistance and arsenic oxidation are localized within the pSinA and pSinB megaplasmids. Comparative genomics revealed that multiple copies of high-affinity phosphate transport systems are common in AOBs, possibly as an As-resistance mechanism. Genome and antibiotic sensitivity analyses further suggested that the use of Ensifer sp. M14 in biotechnology does not pose serious biosafety risks. Therefore, a novel two-stage installation for remediation of arsenic-contaminated waters was developed. It consists of a microbiological module, where M14 oxidizes As(III) to As(V) ion, followed by an adsorption module for As(V) removal using granulated bog iron ores. During a 40-day pilot-scale test in an abandoned gold mine in Zloty Stok (Poland), water leaving the microbiological module generally contained trace amounts of As(III), and dramatic decreases in total arsenic concentrations were observed after passage through the adsorption module. These results demonstrate the usefulness of Ensifer sp. M14 in arsenic removal performed in environmental settings.

Plasmids of Psychrotolerant Polaromonas spp. Isolated From Arctic and Antarctic Glaciers - Diversity and Role in Adaptation to Polar Environments

Cold-active bacteria of the genus Polaromonas (class Betaproteobacteria) are important components of glacial microbiomes. In this study, extrachromosomal replicons of 26 psychrotolerant Polaromonas strains, isolated from Arctic and Antarctic glaciers, were identified, sequenced, and characterized. The plasmidome of these strains consists of 13 replicons, ranging in size from 3,378 to 101,077 bp. In silico sequence analyses identified the conserved backbones of these plasmids, composed of genes required for plasmid replication, stable maintenance, and conjugal transfer. Host range analysis revealed that all of the identified plasmids are narrow-host-range replicons, only able to replicate in bacteria of closely related genera (Polaromonas and Variovorax) of the Comamonadaceae family. Special attention was paid to the identification of plasmid auxiliary genetic information, which may contribute to the adaptation of bacteria to environmental conditions occurring in glaciers. Detailed analysis revealed the presence of genes encoding proteins potentially involved in (i) protection against reactive oxygen species, ultraviolet radiation, and low temperatures; (ii) transport and metabolism of organic compounds; (iii) transport of metal ions; and (iv) resistance to heavy metals. Some of the plasmids also carry genes required for the molecular assembly of iron-sulfur [Fe-S] clusters. Functional analysis of the predicted heavy metal resistance determinants demonstrated that their activity varies, depending on the host strain. This study provides the first molecular insight into the mobile DNA of Polaromonas spp. inhabiting polar glaciers. It has generated valuable data on the structure and properties of a pool of plasmids and highlighted their role in the biology of psychrotolerant Polaromonas strains and their adaptation to the environmental conditions of Arctic and Antarctic glaciers.

Using genome-scale metabolic models to compare serovars of the foodborne pathogen Listeria monocytogenes

Listeria monocytogenes is a microorganism of great concern for the food industry and the cause of human foodborne disease. Therefore, novel methods of control are needed, and systems biology is one such approach to identify them. Using a combination of computational techniques and laboratory methods, genome-scale metabolic models (GEMs) can be created, validated, and used to simulate growth environments and discern metabolic capabilities of microbes of interest, including L. monocytogenes. The objective of the work presented here was to generate GEMs for six different strains of L. monocytogenes, and to both qualitatively and quantitatively validate these GEMs with experimental data to examine the diversity of metabolic capabilities of numerous strains from the three different serovar groups most associated with foodborne outbreaks and human disease. Following qualitative validation, 57 of the 95 carbon sources tested experimentally were present in the GEMs, and; therefore, these were the compounds from which comparisons could be drawn. Of these 57 compounds, agreement between in silico predictions and in vitro results for carbon source utilization ranged from 80.7% to 91.2% between strains. Nutrient utilization agreement between in silico predictions and in vitro results were also conducted for numerous nitrogen, phosphorous, and sulfur sources. Additionally, quantitative validation showed that the L. monocytogenes GEMs were able to generate in silico predictions for growth rate and growth yield that were strongly and significantly (p < 0.0013 and p < 0.0015, respectively) correlated with experimental results. These findings are significant because they show that these GEMs for L. monocytogenes are comparable to published GEMs of other organisms for agreement between in silico predictions and in vitro results. Therefore, as with the other GEMs, namely those for Escherichia coli, Staphylococcus aureus, Vibrio vulnificus, and Salmonella spp., they can be used to determine new methods of growth control and disease treatment.

Inter-replicon Gene Flow Contributes to Transcriptional Integration in the Sinorhizobium meliloti Multipartite Genome

Integration of newly acquired genes into existing regulatory networks is necessary for successful horizontal gene transfer (HGT). Ten percent of bacterial species contain at least two DNA replicons over 300 kilobases in size, with the secondary replicons derived predominately through HGT. The Sinorhizobium meliloti genome is split between a 3.7 Mb chromosome, a 1.7 Mb chromid consisting largely of genes acquired through ancient HGT, and a 1.4 Mb megaplasmid consisting primarily of recently acquired genes. Here, RNA-sequencing is used to examine the transcriptional consequences of massive, synthetic genome reduction produced through the removal of the megaplasmid and/or the chromid. Removal of the pSymA megaplasmid influenced the transcription of only six genes. In contrast, removal of the chromid influenced expression of ∼8% of chromosomal genes and ∼4% of megaplasmid genes. This was mediated in part by the loss of the ETR DNA region whose presence on pSymB is due to a translocation from the chromosome. No obvious functional bias among the up-regulated genes was detected, although genes with putative homologs on the chromid were enriched. Down-regulated genes were enriched in motility and sensory transduction pathways. Four transcripts were examined further, and in each case the transcriptional change could be traced to loss of specific pSymB regions. In particularly, a chromosomal transporter was induced due to deletion of bdhA likely mediated through 3-hydroxybutyrate accumulation. These data provide new insights into the evolution of the multipartite bacterial genome, and more generally into the integration of horizontally acquired genes into the transcriptome.

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.

The complete replicons of 16 Ensifer meliloti strains offer insights into intra- and inter-replicon gene transfer, transposon-associated loci, and repeat elements

Ensifer meliloti (formerly Rhizobium meliloti and Sinorhizobium meliloti) is a model bacterium for understanding legume–rhizobial symbioses. The tripartite genome of E. meliloti consists of a chromosome, pSymA and pSymB, and in some instances strain-specific accessory plasmids. The majority of previous sequencing studies have relied on the use of assemblies generated from short read sequencing, which leads to gaps and assembly errors. Here we used PacBio-based, long-read assemblies and were able to assemble, de novo, complete circular replicons. In this study, we sequenced, de novo-assembled and analysed 10 E. meliloti strains. Sequence comparisons were also done with data from six previously published genomes. We identified genome differences between the replicons, including mol% G+C and gene content, nucleotide repeats, and transposon-associated loci. Additionally, genomic rearrangements both within and between replicons were identified, providing insight into evolutionary processes at the structural level. There were few cases of inter-replicon gene transfer of core genes between the main replicons. Accessory plasmids were more similar to pSymA than to either pSymB or the chromosome, with respect to gene content, transposon content and G+C content. In our population, the accessory plasmids appeared to share an open genome with pSymA, which contains many nodulation- and nitrogen fixation-related genes. This may explain previous observations that horizontal gene transfer has a greater effect on the content of pSymA than pSymB, or the chromosome, and why some rhizobia show unstable nodulation phenotypes on legume hosts.

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.

Rhizobia: from saprophytes to endosymbionts

Rhizobia are some of the best-studied plant microbiota. These oligotrophic Alphaproteobacteria or Betaproteobacteria form symbioses with their legume hosts. Rhizobia must exist in soil and compete with other members of the microbiota before infecting legumes and forming N₂-fixing bacteroids. These dramatic lifestyle and developmental changes are underpinned by large genomes and even more complex pan-genomes, which encompass the whole population and are subject to rapid genetic exchange. The ability to respond to plant signals and chemoattractants and to colonize nutrient-rich roots are crucial for the competitive success of these bacteria. The availability of a large body of genomic, physiological, biochemical and ecological studies makes rhizobia unique models for investigating community interactions and plant colonization.

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.

Techniques for Large-Scale Bacterial Genome Manipulation and Characterization of the Mutants with Respect to In Silico Metabolic Reconstructions

The rate at which all genes within a bacterial genome can be identified far exceeds the ability to characterize these genes. To assist in associating genes with cellular functions, a large-scale bacterial genome deletion approach can be employed to rapidly screen tens to thousands of genes for desired phenotypes. Here, we provide a detailed protocol for the generation of deletions of large segments of bacterial genomes that relies on the activity of a site-specific recombinase. In this procedure, two recombinase recognition target sequences are introduced into known positions of a bacterial genome through single cross-over plasmid integration. Subsequent expression of the site-specific recombinase mediates recombination between the two target sequences, resulting in the excision of the intervening region and its loss from the genome. We further illustrate how this deletion system can be readily adapted to function as a large-scale in vivo cloning procedure, in which the region excised from the genome is captured as a replicative plasmid. We next provide a procedure for the metabolic analysis of bacterial large-scale genome deletion mutants using the Biolog Phenotype MicroArray™ system. Finally, a pipeline is described, and a sample Matlab script is provided, for the integration of the obtained data with a draft metabolic reconstruction for the refinement of the reactions and gene-protein-reaction relationships in a metabolic reconstruction.

Interaction and Regulation of Carbon, Nitrogen, and Phosphorus Metabolisms in Root Nodules of Legumes

Members of the plant family Leguminosae (Fabaceae) are unique in that they have evolved a symbiotic relationship with rhizobia (a group of soil bacteria that can fix atmospheric nitrogen). Rhizobia infect and form root nodules on their specific host plants before differentiating into bacteroids, the symbiotic form of rhizobia. This complex relationship involves the supply of C₄-dicarboxylate and phosphate by the host plants to the microsymbionts that utilize them in the energy-intensive process of fixing atmospheric nitrogen into ammonium, which is in turn made available to the host plants as a source of nitrogen, a macronutrient for growth. Although nitrogen-fixing bacteroids are no longer growing, they are metabolically active. The symbiotic process is complex and tightly regulated by both the host plants and the bacteroids. The metabolic pathways of carbon, nitrogen, and phosphate are heavily regulated in the host plants, as they need to strike a fine balance between satisfying their own needs as well as those of the microsymbionts. A network of transporters for the various metabolites are responsible for the trafficking of these essential molecules between the two partners through the symbiosome membrane (plant-derived membrane surrounding the bacteroid), and these are in turn regulated by various transcription factors that control their expressions under different environmental conditions. Understanding this complex process of symbiotic nitrogen fixation is vital in promoting sustainable agriculture and enhancing soil fertility.

Trade, Diplomacy, and Warfare: The Quest for Elite Rhizobia Inoculant Strains

Rhizobia form symbiotic nitrogen-fixing nodules on leguminous plants, which provides an important source of fixed nitrogen input into the soil ecosystem. The improvement of symbiotic nitrogen fixation is one of the main challenges facing agriculture research. Doing so will reduce the usage of chemical nitrogen fertilizer, contributing to the development of sustainable agriculture practices to deal with the increasing global human population. Sociomicrobiological studies of rhizobia have become a model for the study of the evolution of mutualistic interactions. The exploitation of the wide range of social interactions rhizobia establish among themselves, with the soil and root microbiota, and with the host plant, could constitute a great advantage in the development of a new generation of highly effective rhizobia inoculants. Here, we provide a brief overview of the current knowledge on three main aspects of rhizobia interaction: trade of fixed nitrogen with the plant; diplomacy in terms of communication and possible synergistic effects; and warfare, as antagonism and plant control over symbiosis. Then, we propose new areas of investigation and the selection of strains based on the combination of the genetic determinants for the relevant rhizobia symbiotic behavioral phenotypes.