Identifying functional multi-host shuttle plasmids to advance synthetic biology applications in Mesorhizobium and Bradyrhizobium

Ammonia availability has a crucial role in agriculture as it ensures healthy plant growth and increased crop yields. Since diazotrophs are the only organisms capable of reducing dinitrogen to ammonia, they have great ecological importance and potential to mitigate the environmental and economic costs of synthetic fertilizer use. Rhizobia are especially valuable being that they can engage in nitrogen-fixing symbiotic relationships with legumes, and they demonstrate great diversity and plasticity in genomic and phenotypic traits. However, few rhizobial species have sufficient genetic tractability for synthetic biology applications. This study established a basic genetic toolbox with antibiotic resistance markers, multi-host shuttle plasmids and a streamlined protocol for biparental conjugation with Mesorhizobium and Bradyrhizobium species. We identified two repABC origins of replication from Sinorhizobium meliloti (pSymB) and Rhizobium etli (p42d) that were stable across all three strains of interest. Furthermore, the NZP2235 genome was sequenced and phylogenetic analysis determined its reclassification to Mesorhizobium huakuii. These tools will enable the use of plasmid-based strategies for more advanced genetic engineering projects and ultimately contribute towards the development of more sustainable agriculture practices by means of novel nitrogen-fixing organelles, elite bioinoculants, or symbiotic association with nonlegumes.

Refining the taxonomy of the order Hyphomicrobiales (Rhizobiales) based on whole genome comparisons of over 130 type strains

The alphaproteobacterial order Hyphomicrobiales consists of 38 families comprising at least 152 validly published genera as of January 2024. The order Hyphomicrobiales was first described in 1957 and underwent important revisions in 2020. However, we show that several inconsistencies in the taxonomy of this order remain and we argue that there is a need for a consistent framework for defining families within the order. We propose a common genome-based framework for defining families within the order Hyphomicrobiales, suggesting that families represent monophyletic groups in core-genome phylogenies that share pairwise average amino acid identity values above ~75 % when calculated from a core set of 59 proteins. Applying this framework, we propose the formation of four new families and to reassign the genera Salaquimonas, Rhodoblastus, and Rhodoligotrophos into Salaquimonadaceae fam. nov., Rhodoblastaceae fam. nov., and Rhodoligotrophaceae fam. nov., respectively, and the genera Albibacter, Chenggangzhangella, Hansschlegelia, and Methylopila into Methylopilaceae fam. nov. We further propose to unify the families Bartonellaceae, Brucellaceae, Phyllobacteriaceae, and Notoacmeibacteraceae as Bartonellaceae; the families Segnochrobactraceae and Pseudoxanthobacteraceae as Segnochrobactraceae; the families Lichenihabitantaceae and Lichenibacteriaceae as Lichenihabitantaceae; and the families Breoghaniaceae and Stappiaceae as Stappiaceae. Lastly, we propose to reassign several genera to existing families. Specifically, we propose to reassign the genus Pseudohoeflea to the family Rhizobiaceae; the genera Oricola, Roseitalea, and Oceaniradius to the family Ahrensiaceae; the genus Limoniibacter to the emended family Bartonellaceae; the genus Faunimonas to the family Afifellaceae; and the genus Pseudochelatococcus to the family Chelatococcaceae. Our data also support the recent proposal to reassign the genus Prosthecomicrobium to the family Kaistiaceae.

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.

Diversity and Evolutionary History of Ti Plasmids of "tumorigenes" Clade of Rhizobium spp. and Their Differentiation from Other Ti and Ri Plasmids

Agrobacteria are important plant pathogens responsible for crown/cane gall and hairy root diseases. Crown/cane gall disease is associated with strains carrying tumor-inducing (Ti) plasmids, while hairy root disease is caused by strains harboring root-inducing (Ri) plasmids. In this study, we analyzed the sequences of Ti plasmids of the novel "tumorigenes" clade of the family Rhizobiaceae ("tumorigenes" Ti plasmids), which includes two species, Rhizobium tumorigenes and Rhizobium rhododendri. The sequences of reference Ti/Ri plasmids were also included, which was followed by a comparative analysis of their backbone and accessory regions. The "tumorigenes" Ti plasmids have novel opine signatures compared with other Ti/Ri plasmids characterized so far. The first group exemplified by pTi1078 is associated with production of agrocinopine, nopaline, and ridéopine in plant tumors, while the second group comprising pTi6.2 is responsible for synthesis of leucinopine. Bioinformatic and chemical analyses, including opine utilization assays, indicated that leucinopine associated with pTi6.2 most likely has D,L stereochemistry, unlike the L,L-leucinopine produced in tumors induced by reference strains Chry5 and Bo542. Most of the "tumorigenes" Ti plasmids have conjugative transfer system genes that are unusual for Ti plasmids, composed of avhD4/avhB and traA/mobC/parA regions. Next, our results suggested that "tumorigenes" Ti plasmids have a common origin, but they diverged through large-scale recombination events, through recombination with single or multiple distinct Ti/Ri plasmids. Lastly, we showed that Ti/Ri plasmids could be differentiated based on pairwise Mash or average amino-acid identity distance clustering, and we supply a script to facilitate application of the former approach by other researchers.

Complete Genome Sequences of Rhizobium leguminosarum bv. phaseoli BT01, Rhizobium sp. Strain BT03, and Rhizobium sp. Strain BT04, Isolated from Black Turtle Bean Nodules in Ontario, Canada

We report the complete genome sequences of three Rhizobium strains isolated from nodules of heritage black turtle bean (Phaseolus vulgaris) plants grown in a community garden in Ontario, Canada. The genomes are between 6.91 Mb and 7.98 Mb long and consist of five to seven DNA replicons.

Whole genome assemblies of Zophobas morio and Tenebrio molitor

Zophobas morio (=Zophobas atratus) and Tenebrio molitor are darkling beetles with industrial importance due to their use as feeder insects and their apparent ability to biodegrade plastics. High quality genome assemblies were recently reported for both species. Here, we report additional independent Z. morio and T. molitor genome assemblies generated from Nanopore and Illumina data. Following scaffolding against the published genomes, haploid assemblies of 462 Mb (scaffold N90 of 16.8 Mb) and 258 Mb (scaffold N90 of 5.9 Mb) were produced for Z. morio and T. molitor, respectively. Gene prediction led to the prediction of 28,544 and 19,830 genes for Z. morio and T. molitor, respectively. Benchmarking Universal Single Copy Orthologs (BUSCO) analyses suggested that both assemblies have a high level of completeness; 91.5 and 89.0% of the BUSCO endopterygota marker genes were complete in the Z. morio assembly and proteome, respectively, while 99.1 and 92.8% were complete in the T. molitor assembly and proteome, respectively. Phylogenomic analyses of four genera from the family Tenebrionidae yielded phylogenies consistent with those previously constructed based on mitochondrial genomes. Synteny analyses revealed large stretches of macrosynteny across the family Tenebrionidae, as well as numerous within-chromosome rearrangements. Finally, orthogroup analysis identified ∼28,000 gene families across the family Tenebrionidae, of which 8,185 were identified in all five of the analyzed species, and 10,837 were conserved between Z. morio and T. molitor. We expect that the availability of multiple whole genome sequences for Z. morio and T. molitor will facilitate population genetics studies to identify genetic variation associated with industrially relevant phenotypes.

Complete Genome Sequences of the Species Type Strains Sinorhizobium garamanticum LMG 24692 and Sinorhizobium numidicum LMG 27395 and CIP 109850

The genus Sinorhizobium comprises rhizobia that fix nitrogen in symbiosis with legumes. To support taxonomic studies of this genus and of rhizobia more broadly, we report complete genome sequences and annotations for the species type strains Sinorhizobium garamanticum LMG 24692 and Sinorhizobium numidicum LMG 27395 and CIP 109850.

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.

The evolutionary ecology of rhizobia: multiple facets of competition before, during, and after symbiosis with legumes

Rhizobial bacteria have complex lifestyles that involve growth and survival in bulk soil, plant rhizospheres and rhizoplanes, legume infection threads, and mature and senescing legume nodules. In nature, rhizobia coexist and compete with many other rhizobial strains and species to form host associations. We review recent work defining competitive interactions across these environments. We highlight the use of sophisticated measurement tools and sequencing technologies to examine competition mechanisms in planta, and highlight environments (e.g. soil and senescing nodules) where we still know exceedingly little. We argue that moving toward an explicitly ecological framework (types of competition, resources, and genetic differentiation) will clarify the evolutionary ecology of these foundational organisms and open doors for engineering sustainable, beneficial associations with hosts.

Genomics of the “tumorigenes” clade of the family Rhizobiaceae and description of Rhizobium rhododendri sp. nov

Tumorigenic members of the family Rhizobiaceae, known as agrobacteria, are responsible for crown and cane gall diseases of various crops worldwide. Tumorigenic agrobacteria are commonly found in the genera Agrobacterium, Allorhizobium, and Rhizobium. In this study, we analyzed a distinct “tumorigenes” clade of the genus Rhizobium, which includes the tumorigenic species Rhizobium tumorigenes, as well as strains causing crown gall disease on rhododendron. Here, high‐quality, closed genomes of representatives of the “tumorigenes” clade were generated, followed by comparative genomic and phylogenomic analyses. Additionally, the phenotypic characteristics of representatives of the “tumorigenes” clade were analyzed. Our results showed that the tumorigenic strains isolated from rhododendron represent a novel species of the genus Rhizobium for which the name Rhizobium rhododendri sp. nov. is proposed. This species also includes additional strains originating from blueberry and Himalayan blackberry in the United States, whose genome sequences were retrieved from GenBank. Both R. tumorigenes and R. rhododendri contain multipartite genomes, including a chromosome, putative chromids, and megaplasmids. Synteny and phylogenetic analyses indicated that a large putative chromid of R. rhododendri resulted from the cointegration of an ancestral megaplasmid and two putative chromids, following its divergence from R. tumorigenes. Moreover, gene clusters specific for both species of the “tumorigenes” clade were identified, and their biological functions and roles in the ecological diversification of R. rhododendri and R. tumorigenes were predicted and discussed.

Cold-inducible promoter-driven knockdown of Brachypodium antifreeze proteins confers freezing and phytopathogen susceptibility

The model forage crop, Brachypodium distachyon, has a cluster of ice recrystallization inhibition (BdIRI) genes, which encode antifreeze proteins that function by adsorbing to ice crystals and inhibiting their growth. The genes were targeted for knockdown using a cold-induced promoter from rice (prOsMYB1R35) to drive miRNA. The transgenic lines showed no apparent pleiotropic developmental defects but had reduced antifreeze activity as assessed by assays for ice-recrystallization inhibition, thermal hysteresis, electrolyte leakage, and leaf infrared thermography. Strikingly, the number of cold-acclimated transgenic plants that survived freezing at -8°C was reduced by half or killed entirely, depending on the line, compared with cold-acclimated wild type plants. In addition, more leaf damage was apparent at subzero temperatures in knockdowns after infection with an ice nucleating pathogen, Pseudomonas syringae. Although antifreeze proteins have been studied for almost 60 years, this is the first unequivocal demonstration of their function by knockdown in any organism, and their dual contribution to freeze protection as well as pathogen susceptibility, independent of obvious developmental defects. These proteins are thus of potential interest in a wide range of biotechnological applications from cryopreservation, to frozen product additives, to the engineering of transgenic crops with enhanced pathogen and freezing tolerance.

Whole genome sequencing of mesorhizobia isolated from northern Canada

Rhizobia are soil-dwelling bacteria that can form N2-fixing symbioses with legume plant species (Fabaceae). These bacteria are globally distributed; however, few studies have examined the genomics of rhizobia that live in cold environments. Here, we isolated and characterized three Mesorhizobium strains from legume nodules collected at a pair of distant low Arctic tundra and boreal forest sites in northern Canada. Phylogenetic and average nucleotide identity measurements suggested that the three strains are members of the genus Mesorhizobium, and that each strain represents a novel genospecies. Intriguingly, whereas most mesorhizobia contain the classical determinants of nodulation and nitrogen fixation on their chromosome, whole genome sequencing revealed that all three strains carry these genes on large symbiotic megaplasmids of ~750 to ~1,000 kb. Phylogenetic and sequence analyses of the common nodulation genes revealed highly conserved alleles amongst these northern mesorhizobial strains, leading us to propose that they belong to a novel symbiovar that we termed symbiovar oxytropis. Interestingly, these nod gene alleles are uncommon in mesorhizobia isolated from similar plant hosts in other climatic regions, suggesting potential functional adaptive differences.

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.

Taxonomy of Rhizobiaceae revisited: proposal of a new framework for genus delimitation

The alphaproteobacterial family Rhizobiaceae is highly diverse, with 168 species with validly published names classified into 17 genera with validly published names. Most named genera in this family are delineated based on genomic relatedness and phylogenetic relationships, but some historically named genera show inconsistent distribution and phylogenetic breadth. The most problematic is Rhizobium , which is notorious for being highly paraphyletic, as most newly described species in the family are assigned to this genus without consideration of their proximity to existing genera, or the need to create novel genera. Moreover, many Rhizobiaceae genera lack synapomorphic traits that would give them biological and ecological significance. We propose a common framework for genus delimitation within the family Rhizobiaceae , wherein genera are defined as monophyletic groups in a core-genome gene phylogeny, that are separated from related species using a pairwise core-proteome average amino acid identity (cpAAI) threshold of approximately 86 %. We further propose that additional genomic or phenotypic evidence can justify division of species into separate genera even if they share greater than 86 % cpAAI. Applying this framework, we propose to reclassify Rhizobium rhizosphaerae and Rhizobium oryzae into Xaviernesmea gen. nov. Data is also provided to support the formation of Peteryoungia aggregata comb. nov., Endobacterium yantingense comb. nov., Neorhizobium petrolearium comb. nov., Pararhizobium arenae comb. nov., Pseudorhizobium tarimense comb. nov. and Mycoplana azooxidifex comb. nov. Lastly, we present arguments that the unification of the genera Ensifer and Sinorhizobium in Opinion 84 of the Judicial Commission is no longer justified by current genomic and phenotypic data. Despite pairwise cpAAI values for all Ensifer species and all Sinorhizobium species being >86 %, additional genomic and phenotypic data suggest that they significantly differ in their biology and ecology. We therefore propose emended descriptions of Ensifer and Sinorhizobium , which we argue should be considered as separate genera.

Cold Acclimation in Brachypodium Is Accompanied by Changes in Above-Ground Bacterial and Fungal Communities

Shifts in microbiota undoubtedly support host plants faced with abiotic stress, including low temperatures. Cold-resistant perennials prepare for freeze stress during a period of cold acclimation that can be mimicked by transfer from growing conditions to a reduced photoperiod and a temperature of 4 °C for 2-6 days. After cold acclimation, the model cereal, Brachypodium distachyon, was characterized using metagenomics supplemented with amplicon sequencing (16S ribosomal RNA gene fragments and an internal transcribed spacer region). The bacterial and fungal rhizosphere remained largely unchanged from that of non-acclimated plants. However, leaf samples representing bacterial and fungal communities of the endo- and phyllospheres significantly changed. For example, a plant-beneficial bacterium, Streptomyces sp. M2, increased more than 200-fold in relative abundance in cold-acclimated leaves, and this increase correlated with a striking decrease in the abundance of Pseudomonas syringae (from 8% to zero). This change is of consequence to the host, since P. syringae is a ubiquitous ice-nucleating phytopathogen responsible for devastating frost events in crops. We posit that a responsive above-ground bacterial and fungal community interacts with Brachypodium's low temperature and anti-pathogen signalling networks to help ensure survival in subsequent freeze events, underscoring the importance of inter-kingdom partnerships in the response to cold stress.

The Brachypodium distachyon cold-acclimated plasma membrane proteome is primed for stress resistance

In order to survive subzero temperatures, some plants undergo cold acclimation (CA) where low, nonfreezing temperatures, and/or shortened day lengths allow cold-hardening and survival during subsequent freeze events. Central to this response is the plasma membrane (PM), where low temperature is perceived and cellular homeostasis must be preserved by maintaining membrane integrity. Here, we present the first PM proteome of cold-acclimated Brachypodium distachyon, a model species for the study of monocot crops. A time-course experiment investigated CA-induced changes in the proteome following two-phase partitioning PM enrichment and label-free quantification by nano-liquid chromatography-mass spectrophotometry. Two days of CA were sufficient for membrane protection as well as an initial increase in sugar levels and coincided with a significant change in the abundance of 154 proteins. Prolonged CA resulted in further increases in soluble sugars and abundance changes in more than 680 proteins, suggesting both a necessary early response to low-temperature treatment, as well as a sustained CA response elicited over several days. A meta-analysis revealed that the identified PM proteins have known roles in low-temperature tolerance, metabolism, transport, and pathogen defense as well as drought, osmotic stress, and salt resistance suggesting crosstalk between stress responses, such that CA may prime plants for other abiotic and biotic stresses. The PM proteins identified here present keys to an understanding of cold tolerance in monocot crops and the hope of addressing economic losses associated with modern climate-mediated increases in frost events.

Proteobacteria Contain Diverse flg22 Epitopes That Elicit Varying Immune Responses in Arabidopsis thaliana

Bacterial flagellin protein is a potent microbe-associated molecular pattern. Immune responses are triggered by a 22-amino-acid epitope derived from flagellin, known as flg22, upon detection by the pattern recognition receptor FLAGELLIN-SENSING2 (FLS2) in multiple plant species. However, increasing evidence suggests that flg22 epitopes of several bacterial species are not universally immunogenic to plants. We investigated whether flg22 immunogenicity systematically differs between classes of the phylum Proteobacteria, using a dataset of 2,470 flg22 sequences. To predict which species encode highly immunogenic flg22 epitopes, we queried a custom motif (¹¹[ST]xx[DN][DN]xAGxxI²¹) in the flg22 sequences, followed by sequence conservation analysis and protein structural modeling. These data led us to hypothesize that most flg22 epitopes of the γ- and β-Proteobacteria are highly immunogenic, whereas most flg22 epitopes of the α-, δ-, and ε-Proteobacteria are weakly to moderately immunogenic. To test this hypothesis, we generated synthetic peptides representative of the flg22 epitopes of each proteobacterial class, and we monitored their ability to elicit an immune response in Arabidopsis thaliana. The flg22 peptides of γ- and β-Proteobacteria triggered strong oxidative bursts, whereas peptides from the ε-, δ-, and α-Proteobacteria triggered moderate, weak, or no response, respectively. These data suggest flg22 immunogenicity is not highly conserved across the phylum Proteobacteria. We postulate that sequence divergence of each taxonomic class was present prior to the evolution of FLS2, and that the ligand specificity of A. thaliana FLS2 was driven by the flg22 epitopes of the γ- and β-Proteobacteria, a monophyletic group containing many common phytopathogens.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

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.

Rhizobium indicum sp. nov., isolated from root nodules of pea (Pisum sativum) cultivated in the Indian trans-Himalayas

Three strains of rhizobia isolated from effective root nodules of pea (Pisum sativum L.) collected from the Indian trans-Himalayas were characterized using 16S rRNA, atpD and recA genes. Phylogeny of the 16S rRNA genes revealed that the newly isolated strains were members of the genus Rhizobium with ≥99.9% sequence similarity to the members within the "Rhizobium leguminosarum" group. Phylogenetic analyses based on the concatenated sequences of atpD and recA gene, and 92 core genes extracted from the genome sequences indicated that strains JKLM 12A2T and JKLM 13E are grouped as a separate clade closely related to R. laguerreae FB206T. In contrast, the strain JKLM 19E was placed with "R. hidalgonense" FH14T. Whole-genome average nucleotide identity (ANI) values were 97.6% within strains JKLM 12A2T and JKLM 13E, and less than 94% with closely related species. The digital DNA-DNA hybridization (dDDH) values were 81.45 within the two strains and less than 54.8% to closely related species. The major cellular fatty acids were C18:1w7c in summed feature 8, C14:0 3OH/C16:1 iso I in summed feature 2, and C18:0. The DNA G+C content of JKLM 12A2T and JKLM 13E was 60.8mol%. The data on genomic, chemotaxonomic, and phenotypic characteristics indicates that the strains JKLM 12A2T and JKLM 13E represent a novel species, Rhizobium indicum sp. nov. The type strain is JKLM 12A2T (= MCC 3961T=KACC 21380T=JCM 33658T). However, the strain JKLM 19E represents a member of "R. hidalgonense" and the symbiovar viciae.

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.

Chromids Aid Genome Expansion and Functional Diversification in the Family Burkholderiaceae

Multipartite genomes, containing at least two large replicons, are found in diverse bacteria; however, the advantage of this genome structure remains incompletely understood. Here, we perform comparative genomics of hundreds of finished β-proteobacterial genomes to gain insights into the role and emergence of multipartite genomes. Almost all essential secondary replicons (chromids) of the β-proteobacteria are found in the family Burkholderiaceae. These replicons arose from just two plasmid acquisition events, and they were likely stabilized early in their evolution by the presence of core genes. On average, Burkholderiaceae genera with multipartite genomes had a larger total genome size, but smaller chromosome, than genera without secondary replicons. Pangenome-level functional enrichment analyses suggested that interreplicon functional biases are partially driven by the enrichment of secondary replicons in the accessory pangenome fraction. Nevertheless, the small overlap in orthologous groups present in each replicon's pangenome indicated a clear functional separation of the replicons. Chromids appeared biased to environmental adaptation, as the functional categories enriched on chromids were also overrepresented on the chromosomes of the environmental genera (Paraburkholderia and Cupriavidus) compared with the pathogenic genera (Burkholderia and Ralstonia). Using ancestral state reconstruction, it was predicted that the rate of accumulation of modern-day genes by chromids was more rapid than the rate of gene accumulation by the chromosomes. Overall, the data are consistent with a model where the primary advantage of secondary replicons is in facilitating increased rates of gene acquisition through horizontal gene transfer, consequently resulting in replicons enriched in genes associated with adaptation to novel environments.

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.

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.

Pulse-Modulated Radio-Frequency Alternating-Current-Driven Atmospheric-Pressure Glow Discharge for Continuous-Flow Synthesis of Silver Nanoparticles and Evaluation of Their Cytotoxicity toward Human Melanoma Cells

An innovative and environmentally friendly method for the synthesis of size-controlled silver nanoparticles (AgNPs) is presented. Pectin-stabilized AgNPs were synthesized in a plasma-reaction system in which pulse-modulated radio-frequency atmospheric-pressure glow discharge (pm-rf-APGD) was operated in contact with a flowing liquid electrode. The use of pm-rf-APGD allows for better control of the size of AgNPs and their stability and monodispersity. AgNPs synthesized under defined operating conditions exhibited average sizes of 41.62 ± 12.08 nm and 10.38 ± 4.56 nm, as determined by dynamic light scattering and transmission electron microscopy (TEM), respectively. Energy-dispersive X-ray spectroscopy (EDS) confirmed that the nanoparticles were composed of metallic Ag. Furthermore, the ξ-potential of the AgNPs was shown to be -43.11 ± 0.96 mV, which will facilitate their application in biological systems. Between 70% and 90% of the cancerous cells of the human melanoma Hs 294T cell line underwent necrosis following treatment with the synthesized AgNPs. Furthermore, optical emission spectrometry (OES) identified reactive species, such as NO, NH, N₂, O, and H, as pm-rf-APGD produced compounds that may be involved in the reduction of the Ag(I) ions.

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.

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.

The Divided Bacterial Genome: Structure, Function, and Evolution

Approximately 10% of bacterial genomes are split between two or more large DNA fragments, a genome architecture referred to as a multipartite genome. This multipartite organization is found in many important organisms, including plant symbionts, such as the nitrogen-fixing rhizobia, and plant, animal, and human pathogens, including the genera Brucella, Vibrio, and Burkholderia. The availability of many complete bacterial genome sequences means that we can now examine on a broad scale the characteristics of the different types of DNA molecules in a genome. Recent work has begun to shed light on the unique properties of each class of replicon, the unique functional role of chromosomal and nonchromosomal DNA molecules, and how the exploitation of novel niches may have driven the evolution of the multipartite genome. The aims of this review are to (i) outline the literature regarding bacterial genomes that are divided into multiple fragments, (ii) provide a meta-analysis of completed bacterial genomes from 1,708 species as a way of reviewing the abundant information present in these genome sequences, and (iii) provide an encompassing model to explain the evolution and function of the multipartite genome structure. This review covers, among other topics, salient genome terminology; mechanisms of multipartite genome formation; the phylogenetic distribution of multipartite genomes; how each part of a genome differs with respect to genomic signatures, genetic variability, and gene functional annotation; how each DNA molecule may interact; as well as the costs and benefits of this genome structure.

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

The bacterium Sinorhizobium meliloti Rm2011 forms N2-fixing root nodules on alfalfa and other leguminous plants. The pSymB chromid contains a 110-kb region (the ETR region) showing high synteny to a chromosomally located region in Sinorhizobium fredii NGR234 and related rhizobia. We recently introduced the ETR region from S. fredii NGR234 into the S. meliloti chromosome. Here, we report that, unexpectedly, the S. fredii NGR234 ETR region did not complement deletion of the S. meliloti ETR region in symbiosis with Medicago sativa. This phenotype was due to the bacA gene of NGR234 not being functionally interchangeable with the S. meliloti bacA gene during M. sativa symbiosis. Further analysis revealed that, whereas bacA genes from S. fredii or Rhizobium leguminosarum bv. viciae 3841 failed to complement the Fix- phenotype of a S. meliloti bacA mutant with M. sativa, they allowed for further developmental progression prior to a loss of viability. In contrast, with Melilotus alba, bacA from S. fredii and R. leguminosarum supported N2 fixation by a S. meliloti bacA mutant. Additionally, the S. meliloti bacA gene can support N2 fixation of a R. leguminosarum bacA mutant during symbiosis with Pisum sativum. A phylogeny of BacA proteins illustrated that S. meliloti BacA has rapidly diverged from most rhizobia and has converged toward the sequence of pathogenic genera Brucella and Escherichia. These data suggest that the S. meliloti BacA has evolved toward a specific interaction with Medicago and highlights the limitations of using a single model system for the study of complex biological topics.

A putative 3-hydroxyisobutyryl-CoA hydrolase is required for efficient symbiotic nitrogen fixation in Sinorhizobium meliloti and Sinorhizobium fredii NGR234

We report that the smb20752 gene of the alfalfa symbiont Sinorhizobium meliloti is a novel symbiotic gene required for full N₂ -fixation. Deletion of smb20752 resulted in lower nitrogenase activity and smaller nodules without impacting overall nodule morphology. Orthologs of smb20752 were present in all alpha and beta rhizobia, including the ngr_b20860 gene of Sinorhizobium fredii NGR234. A ngr_b20860 mutant formed Fix^- determinate nodules that developed normally to a late stage of the symbiosis on the host plants Macroptilium atropurpureum and Vigna unguiculata. However an early symbiotic defect was evident during symbiosis with Leucaena leucocephala, producing Fix^- indeterminate nodules. The smb20752 and ngr_b20860 genes encode putative 3-hydroxyisobutyryl-CoA (HIB-CoA) hydrolases. HIB-CoA hydrolases are required for l-valine catabolism and appear to prevent the accumulation of toxic metabolic intermediates, particularly methacrylyl-CoA. Evidence presented here and elsewhere (Curson et al., , PLoS ONE 9:e97660) demonstrated that Smb20752 and NGR_b20860 can also prevent metabolic toxicity, are required for l-valine metabolism, and play an undefined role in 3-hydroxybutyrate catabolism. We present evidence that the symbiotic defect of the HIB-CoA hydrolase mutants is independent of the inability to catabolize l-valine and suggest it relates to the toxicity resulting from metabolism of other compounds possibly related to 3-hydroxybutyric acid.

Loss of malic enzymes leads to metabolic imbalance and altered levels of trehalose and putrescine in the bacterium Sinorhizobium meliloti

BACKGROUND

Malic enzymes decarboxylate the tricarboxylic acid (TCA) cycle intermediate malate to the glycolytic end-product pyruvate and are well positioned to regulate metabolic flux in central carbon metabolism. Despite the wide distribution of these enzymes, their biological roles are unclear in part because the reaction catalyzed by these enzymes can be by-passed by other pathways. The N2-fixing alfalfa symbiont Sinorhizobium meliloti contains both a NAD(P)-malic enzyme (DME) and a separate NADP-malic enzyme (TME) and to help understand the role of these enzymes, we investigated growth, metabolomic, and transcriptional consequences resulting from loss of these enzymes in free-living cells.

RESULTS

Loss of DME, TME, or both enzymes had no effect on growth with the glycolytic substrate, glucose. In contrast, the dme mutants, but not tme, grew slowly on the gluconeogenic substrate succinate and this slow growth was further reduced upon the addition of glucose. The dme mutant strains incubated with succinate accumulated trehalose and hexose sugar phosphates, secreted malate, and relative to wild-type, these cells had moderately increased transcription of genes involved in gluconeogenesis and pathways that divert metabolites away from the TCA cycle. While tme mutant cells grew at the same rate as wild-type on succinate, they accumulated the compatible solute putrescine.

CONCLUSIONS

NAD(P)-malic enzyme (DME) of S. meliloti is required for efficient metabolism of succinate via the TCA cycle. In dme mutants utilizing succinate, malate accumulates and is excreted and these cells appear to increase metabolite flow via gluconeogenesis with a resulting increase in the levels of hexose-6-phosphates and trehalose. For cells utilizing succinate, TME activity alone appeared to be insufficient to produce the levels of pyruvate required for efficient TCA cycle metabolism. Putrescine was found to accumulate in tme cells growing with succinate, and whether this is related to altered levels of NADPH requires further investigation.

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

The genome of about 10% of bacterial species is divided among two or more large chromosome-sized replicons. The contribution of each replicon to the microbial life cycle (for example, environmental adaptations and/or niche switching) remains unclear. Here we report a genome-scale metabolic model of the legume symbiont Sinorhizobium meliloti that is integrated with carbon utilization data for 1,500 genes with 192 carbon substrates. Growth of S. meliloti is modelled in three ecological niches (bulk soil, rhizosphere and nodule) with a focus on the role of each of its three replicons. We observe clear metabolic differences during growth in the tested ecological niches and an overall reprogramming following niche switching. In silico examination of the inferred fitness of gene deletion mutants suggests that secondary replicons evolved to fulfil a specialized function, particularly host-associated niche adaptation. Thus, genes on secondary replicons might potentially be manipulated to promote or suppress host interactions for biotechnological purposes.

The genome of some bacteria consists of two or more chromosomes or replicons. Here, diCenzo et al. integrate genome-scale metabolic modelling and growth data from a collection of mutants of the plant symbiont Sinorhizobium meliloti to estimate the fitness contribution of each replicon in three environments.

Genomic resources for identification of the minimal N2 -fixing symbiotic genome

The lack of an appropriate genomic platform has precluded the use of gain-of-function approaches to study the rhizobium-legume symbiosis, preventing the establishment of the genes necessary and sufficient for symbiotic nitrogen fixation (SNF) and potentially hindering synthetic biology approaches aimed at engineering this process. Here, we describe the development of an appropriate system by reverse engineering Sinorhizobium meliloti. Using a novel in vivo cloning procedure, the engA-tRNA-rmlC (ETR) region, essential for cell viability and symbiosis, was transferred from Sinorhizobium fredii to the ancestral location on the S. meliloti chromosome, rendering the ETR region on pSymB redundant. A derivative of this strain lacking both the large symbiotic replicons (pSymA and pSymB) was constructed. Transfer of pSymA and pSymB back into this strain restored symbiotic capabilities with alfalfa. To delineate the location of the single-copy genes essential for SNF on these replicons, we screened a S. meliloti deletion library, representing > 95% of the 2900 genes of the symbiotic replicons, for their phenotypes with alfalfa. Only four loci, accounting for < 12% of pSymA and pSymB, were essential for SNF. These regions will serve as our preliminary target of the minimal set of horizontally acquired genes necessary and sufficient for SNF.

Proline auxotrophy in Sinorhizobium meliloti results in a plant-specific symbiotic phenotype

In order to effectively manipulate rhizobium-legume symbioses for our benefit, it is crucial to first gain a complete understanding of the underlying genetics and metabolism. Studies with rhizobium auxotrophs have provided insight into the requirement for amino acid biosynthesis during the symbiosis; however, a paucity of available L-proline auxotrophs has limited our understanding of the role of L-proline biosynthesis. Here, we examined the symbiotic phenotypes of a recently described Sinorhizobium meliloti L-proline auxotroph. Proline auxotrophy was observed to result in a host-plant-specific phenotype. The S. meliloti auxotroph displayed reduced symbiotic capability with alfalfa (Medicago sativa) due to a decrease in nodule mass formed and therefore a reduction in nitrogen fixed per plant. However, the proline auxotroph formed nodules on white sweet clover (Melilotus alba) that failed to fix nitrogen. The rate of white sweet clover nodulation by the auxotroph was slightly delayed, but the final number of nodules per plant was not impacted. Examination of white sweet clover nodules by confocal microscopy and transmission electron microscopy revealed the presence of the S. meliloti proline auxotroph cells within the host legume cells, but few differentiated bacteroids were identified compared with the bacteroid-filled plant cells of WT nodules. Overall, these results indicated that L-proline biosynthesis is a general requirement for a fully effective nitrogen-fixing symbiosis, likely due to a transient requirement during bacteroid differentiation.

Genetic redundancy is prevalent within the 6.7 Mb Sinorhizobium meliloti genome

Biological pathways are frequently identified via a genetic loss-of-function approach. While this approach has proven to be powerful, it is imperfect as illustrated by well-studied pathways continuing to have missing steps. One potential limiting factor is the masking of phenotypes through genetic redundancy. The prevalence of genetic redundancy in bacterial species has received little attention, although isolated examples of functionally redundant gene pairs exist. Here, we made use of a strain of Sinorhizobium meliloti whose genome was reduced by 45 % through the complete removal of a megaplasmid and a chromid (3 Mb of the 6.7 Mb genome was removed) to begin quantifying the level of genetic redundancy within a large bacterial genome. A mutagenesis of the strain with the reduced genome identified a set of transposon insertions precluding growth of this strain on minimal medium. Transfer of these mutations to the wild-type background revealed that 10-15 % of these chromosomal mutations were located within duplicated genes, as they did not prevent growth of cells with the full genome. The functionally redundant genes were involved in a variety of metabolic pathways, including central carbon metabolism, transport, and amino acid biosynthesis. These results indicate that genetic redundancy may be prevalent within large bacterial genomes. Failing to account for redundantly encoded functions in loss-of-function studies will impair our understanding of a broad range of biological processes and limit our ability to use synthetic biology in the construction of designer cell factories.

Examination of prokaryotic multipartite genome evolution through experimental genome reduction

Many bacteria carry two or more chromosome-like replicons. This occurs in pathogens such as Vibrio cholerea and Brucella abortis as well as in many N₂-fixing plant symbionts including all isolates of the alfalfa root-nodule bacteria Sinorhizobium meliloti. Understanding the evolution and role of this multipartite genome organization will provide significant insight into these important organisms; yet this knowledge remains incomplete, in part, because technical challenges of large-scale genome manipulations have limited experimental analyses. The distinct evolutionary histories and characteristics of the three replicons that constitute the S. meliloti genome (the chromosome (3.65 Mb), pSymA megaplasmid (1.35 Mb), and pSymB chromid (1.68 Mb)) makes this a good model to examine this topic. We transferred essential genes from pSymB into the chromosome, and constructed strains that lack pSymB as well as both pSymA and pSymB. This is the largest reduction (45.4%, 3.04 megabases, 2866 genes) of a prokaryotic genome to date and the first removal of an essential chromid. Strikingly, strains lacking pSymA and pSymB (ΔpSymAB) lost the ability to utilize 55 of 74 carbon sources and various sources of nitrogen, phosphorous and sulfur, yet the ΔpSymAB strain grew well in minimal salts media and in sterile soil. This suggests that the core chromosome is sufficient for growth in a bulk soil environment and that the pSymA and pSymB replicons carry genes with more specialized functions such as growth in the rhizosphere and interaction with the plant. These experimental data support a generalized evolutionary model, in which non-chromosomal replicons primarily carry genes with more specialized functions. These large secondary replicons increase the organism's niche range, which offsets their metabolic burden on the cell (e.g. pSymA). Subsequent co-evolution with the chromosome then leads to the formation of a chromid through the acquisition of functions core to all niches (e.g. pSymB).

Rhizobia are free-living bacteria of agricultural and environmental importance that form root-nodules on leguminous plants and provide these plants with fixed nitrogen. Many of the rhizobia have a multipartite genome, as do several plant and animal pathogens. All isolates of the alfalfa symbiont, Sinorhizobium meliloti, carry three large replicons, the chromosome (∼3.7 Mb), pSymA megaplasmid (∼1.4 Mb), and pSymB chromid (∼1.7 Mb). To gain insight into the role and evolutionary history of these replicons, we have ‘reversed evolution’ by constructing a S. meliloti strain consisting solely of the chromosome and lacking the pSymB chromid and pSymA megaplasmid. As the resulting strain was viable, we could perform a detailed phenotypic analysis and these data provided significant insight into the biology and metabolism of S. meliloti. The data lend direct experimental evidence in understanding the evolution and role of the multipartite genome. Specifically the large secondary replicons increase the organism's niche range, and this advantage offsets the metabolic burden of these replicons on the cell. Additionally, the single-chromosome strain offers a useful platform to facilitate future forward genetic approaches to understanding and manipulating the symbiosis and plant-microbe interactions.