Nodulating strains of Rhizobium loti arise through chromosomal symbiotic gene transfer in the environment

Beav: a bacterial genome and mobile element annotation pipeline

Comprehensive and accurate genome annotation is crucial for inferring the predicted functions of an organism. Numerous tools exist to annotate genes, gene clusters, mobile genetic elements, and other diverse features. However, these tools and pipelines can be difficult to install and run, be specialized for a particular element or feature, or lack annotations for larger elements that provide important genomic context. Integrating results across analyses is also important for understanding gene function. To address these challenges, we present the Beav annotation pipeline. Beav is a command-line tool that automates the annotation of bacterial genome sequences, mobile genetic elements, molecular systems and gene clusters, key regulatory features, and other elements. Beav uses existing tools in addition to custom models, scripts, and databases to annotate diverse elements, systems, and sequence features. Custom databases for plant-associated microbes are incorporated to improve annotation of key virulence and symbiosis genes in agriculturally important pathogens and mutualists. Beav includes an optional Agrobacterium-specific pipeline that identifies and classifies oncogenic plasmids and annotates plasmid-specific features. Following the completion of all analyses, annotations are consolidated to produce a single comprehensive output. Finally, Beav generates publication-quality genome and plasmid maps. Beav is on Bioconda and is available for download at https://github.com/weisberglab/beav.

IMPORTANCE

Annotation of genome features, such as the presence of genes and their predicted function, or larger loci encoding secretion systems or biosynthetic gene clusters, is necessary for understanding the functions encoded by an organism. Genomes can also host diverse mobile genetic elements, such as integrative and conjugative elements and/or phages, that are often not annotated by existing pipelines. These elements can horizontally mobilize genes encoding for virulence, antimicrobial resistance, or other adaptive functions and alter the phenotype of an organism. We developed a software pipeline, called Beav, that combines new and existing tools for the comprehensive annotation of these and other major features. Existing pipelines often misannotate loci important for virulence or mutualism in plant-associated bacteria. Beav includes custom databases and optional workflows for the improved annotation of plant-associated bacteria. Beav is designed to be easy to install and run, making comprehensive genome annotation broadly available to the research community.

The role of microbial interactions on rhizobial fitness

Rhizobia are soil bacteria that can establish a nitrogen-fixing symbiosis with legume plants. As horizontally transmitted symbionts, the life cycle of rhizobia includes a free-living phase in the soil and a plant-associated symbiotic phase. Throughout this life cycle, rhizobia are exposed to a myriad of other microorganisms that interact with them, modulating their fitness and symbiotic performance. In this review, we describe the diversity of interactions between rhizobia and other microorganisms that can occur in the rhizosphere, during the initiation of nodulation, and within nodules. Some of these rhizobia-microbe interactions are indirect, and occur when the presence of some microbes modifies plant physiology in a way that feeds back on rhizobial fitness. We further describe how these interactions can impose significant selective pressures on rhizobia and modify their evolutionary trajectories. More extensive investigations on the eco-evolutionary dynamics of rhizobia in complex biotic environments will likely reveal fascinating new aspects of this well-studied symbiotic interaction and provide critical knowledge for future agronomical applications.

Mobile Genetic Element Flexibility as an Underlying Principle to Bacterial Evolution

Mobile genetic elements are key to the evolution of bacteria and traits that affect host and ecosystem health. Here, we use a framework of a hierarchical and modular system that scales from genes to populations to synthesize recent findings on mobile genetic elements (MGEs) of bacteria. Doing so highlights the role that emergent properties of flexibility, robustness, and genetic capacitance of MGEs have on the evolution of bacteria. Some of their traits can be stored, shared, and diversified across different MGEs, taxa of bacteria, and time. Collectively, these properties contribute to maintaining functionality against perturbations while allowing changes to accumulate in order to diversify and give rise to new traits. These properties of MGEs have long challenged our abilities to study them. Implementation of new technologies and strategies allows for MGEs to be analyzed in new and powerful ways.

Identifying and tracking mobile elements in evolving compost communities yields insights into the nanobiome

Microbial evolution is driven by rapid changes in gene content mediated by horizontal gene transfer (HGT). While mobile genetic elements (MGEs) are important drivers of gene flux, the nanobiome-the zoo of Darwinian replicators that depend on microbial hosts-remains poorly characterised. New approaches are necessary to increase our understanding beyond MGEs shaping individual populations, towards their impacts on complex microbial communities. A bioinformatic pipeline (xenoseq) was developed to cross-compare metagenomic samples from microbial consortia evolving in parallel, aimed at identifying MGE dissemination, which was applied to compost communities which underwent periodic mixing of MGEs. We show that xenoseq can distinguish movement of MGEs from demographic changes in community composition that otherwise confounds identification, and furthermore demonstrate the discovery of various unexpected entities. Of particular interest was a nanobacterium of the candidate phylum radiation (CPR) which is closely related to a species identified in groundwater ecosystems (Candidatus Saccharibacterium), and appears to have a parasitic lifestyle. We also highlight another prolific mobile element, a 313 kb plasmid hosted by a Cellvibrio lineage. The host was predicted to be capable of nitrogen fixation, and acquisition of the plasmid coincides with increased ammonia production. Taken together, our data show that new experimental strategies combined with bioinformatic analyses of metagenomic data stand to provide insight into the nanobiome as a driver of microbial community evolution.

DUF2285 is a novel helix-turn-helix domain variant that orchestrates both activation and antiactivation of conjugative element transfer in proteobacteria

Horizontal gene transfer is tightly regulated in bacteria. Often only a fraction of cells become donors even when regulation of horizontal transfer is coordinated at the cell population level by quorum sensing. Here, we reveal the widespread 'domain of unknown function' DUF2285 represents an 'extended-turn' variant of the helix-turn-helix domain that participates in both transcriptional activation and antiactivation to initiate or inhibit horizontal gene transfer. Transfer of the integrative and conjugative element ICEMlSymR7A is controlled by the DUF2285-containing transcriptional activator FseA. One side of the DUF2285 domain of FseA has a positively charged surface which is required for DNA binding, while the opposite side makes critical interdomain contacts with the N-terminal FseA DUF6499 domain. The QseM protein is an antiactivator of FseA and is composed of a DUF2285 domain with a negative surface charge. While QseM lacks the DUF6499 domain, it can bind the FseA DUF6499 domain and prevent transcriptional activation by FseA. DUF2285-domain proteins are encoded on mobile elements throughout the proteobacteria, suggesting regulation of gene transfer by DUF2285 domains is a widespread phenomenon. These findings provide a striking example of how antagonistic domain paralogues have evolved to provide robust molecular control over the initiation of horizontal gene transfer.

Discordant population structure among rhizobium divided genomes and their legume hosts

Symbiosis often occurs between partners with distinct life history characteristics and dispersal mechanisms. Many bacterial symbionts have genomes comprising multiple replicons with distinct rates of evolution and horizontal transmission. Such differences might drive differences in population structure between hosts and symbionts and among the elements of the divided genomes of bacterial symbionts. These differences might, in turn, shape the evolution of symbiotic interactions and bacterial evolution. Here we use whole genome resequencing of a hierarchically structured sample of 191 strains of Sinorhizobium meliloti collected from 21 locations in southern Europe to characterize population structures of this bacterial symbiont, which forms a root nodule symbiosis with the host plant Medicago truncatula. S. meliloti genomes showed high local (within-site) variation and little isolation by distance. This was particularly true for the two symbiosis elements, pSymA and pSymB, which have population structures that are similar to each other, but distinct from both the bacterial chromosome and the host plant. Given limited recombination on the chromosome, compared to the symbiosis elements, distinct population structures may result from differences in effective gene flow. Alternatively, positive or purifying selection, with little recombination, may explain distinct geographical patterns at the chromosome. Discordant population structure between hosts and symbionts indicates that geographically and genetically distinct host populations in different parts of the range might interact with genetically similar symbionts, potentially minimizing local specialization.

Lupin, a Unique Legume That Is Nodulated by Multiple Microsymbionts: The Role of Horizontal Gene Transfer

Lupin is a high-protein legume crop that grows in a wide range of edaphoclimatic conditions where other crops are not viable. Its unique seed nutrient profile can promote health benefits, and it has been proposed as a phytoremediation plant. Most rhizobia nodulating Lupinus species belong to the genus Bradyrhizobium, comprising strains that are phylogenetically related to B. cytisi, B. hipponenese, B. rifense, B. iriomotense/B. stylosanthis, B. diazoefficiens, B. japonicum, B. canariense/B. lupini, and B. retamae/B. valentinum. Lupins are also nodulated by fast-growing bacteria within the genera Microvirga, Ochrobactrum, Devosia, Phyllobacterium, Agrobacterium, Rhizobium, and Neorhizobium. Phylogenetic analyses of the nod and nif genes, involved in microbial colonization and symbiotic nitrogen fixation, respectively, suggest that fast-growing lupin-nodulating bacteria have acquired their symbiotic genes from rhizobial genera other than Bradyrhizobium. Horizontal transfer represents a key mechanism allowing lupin to form symbioses with bacteria that were previously considered as non-symbiotic or unable to nodulate lupin, which might favor lupin’s adaptation to specific habitats. The characterization of yet-unstudied Lupinus species, including microsymbiont whole genome analyses, will most likely expand and modify the current lupin microsymbiont taxonomy, and provide additional knowledge that might help to further increase lupin’s adaptability to marginal soils and climates.

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.

Population genomics of Australian indigenous Mesorhizobium reveals diverse nonsymbiotic genospecies capable of nitrogen-fixing symbioses following horizontal gene transfer

Mesorhizobia are soil bacteria that establish nitrogen-fixing symbioses with various legumes. Novel symbiotic mesorhizobia frequently evolve following horizontal transfer of symbiosis-gene-carrying integrative and conjugative elements (ICESyms) to indigenous mesorhizobia in soils. Evolved symbionts exhibit a wide range in symbiotic effectiveness, with some fixing nitrogen poorly or not at all. Little is known about the genetic diversity and symbiotic potential of indigenous soil mesorhizobia prior to ICESym acquisition. Here we sequenced genomes of 144 Mesorhizobium spp. strains cultured directly from cultivated and uncultivated Australian soils. Of these, 126 lacked symbiosis genes. The only isolated symbiotic strains were either exotic strains used previously as legume inoculants, or indigenous mesorhizobia that had acquired exotic ICESyms. No native symbiotic strains were identified. Indigenous nonsymbiotic strains formed 22 genospecies with phylogenomic diversity overlapping the diversity of internationally isolated symbiotic Mesorhizobium spp. The genomes of indigenous mesorhizobia exhibited no evidence of prior involvement in nitrogen-fixing symbiosis, yet their core genomes were similar to symbiotic strains and they generally lacked genes for synthesis of biotin, nicotinate and thiamine. Genomes of nonsymbiotic mesorhizobia harboured similar mobile elements to those of symbiotic mesorhizobia, including ICESym-like elements carrying aforementioned vitamin-synthesis genes but lacking symbiosis genes. Diverse indigenous isolates receiving ICESyms through horizontal gene transfer formed effective symbioses with Lotus and Biserrula legumes, indicating most nonsymbiotic mesorhizobia have an innate capacity for nitrogen-fixing symbiosis following ICESym acquisition. Non-fixing ICESym-harbouring strains were isolated sporadically within species alongside effective symbionts, indicating chromosomal lineage does not predict symbiotic potential. Our observations suggest previously observed genomic diversity amongst symbiotic Mesorhizobium spp. represents a fraction of the extant diversity of nonsymbiotic strains. The overlapping phylogeny of symbiotic and nonsymbiotic clades suggests major clades of Mesorhizobium diverged prior to introduction of symbiosis genes and therefore chromosomal genes involved in symbiosis have evolved largely independent of nitrogen-fixing symbiosis.

Two New Rhizobiales Species Isolated from Root Nodules of Common Sainfoin (Onobrychis viciifolia) Show Different Plant Colonization Strategies

Root nodules of legume plants are primarily inhabited by rhizobial nitrogen-fixing bacteria. Here, we propose two new Rhizobiales species isolated from root nodules of common sainfoin (Onobrychis viciifolia), as shown by core-gene phylogeny, overall genome relatedness indices, and pan-genome analysis. Mesorhizobium onobrychidis sp. nov. actively induces nodules and achieves atmospheric nitrogen and carbon dioxide fixation. This species appears to be depleted in motility genes and is enriched in genes for direct effects on plant growth performance. Its genome reveals functional and plant growth-promoting signatures, like a large unique chromosomal genomic island with high density of symbiotic genetic traits. Onobrychidicola muellerharveyae gen. nov. sp. nov. is described as a type species of the new genus Onobrychidicola in Rhizobiaceae. This species comprises unique genetic features and plant growth-promoting traits (PGPTs), which strongly indicate its function in biotic stress reduction and motility. We applied a newly developed bioinformatics approach for in silico prediction of PGPTs (PGPT-Pred), which supports the different lifestyles of the two new species and the plant growth-promoting performance of M. onobrychidis in the greenhouse trial. IMPORTANCE The intensive use of chemical fertilizers has a variety of negative effects on the environment. Increased utilization of biological nitrogen fixation (BNF) is one way to mitigate those negative impacts. In order to optimize BNF, suitable candidates for different legume species are required. Despite intensive search for new rhizobial bacteria associated with legumes, no new rhizobia have recently been identified from sainfoin (Onobrychis viciifolia). Here, we report on the discovery of two new rhizobial species associated with sainfoin, which are of high importance for the host and may help to increase sustainability in agricultural practices. We employed the combination of in silico prediction and in planta experiments, which is an effective way to detect promising plant growth-promoting bacteria.

Microevolution, speciation and macroevolution in rhizobia: Genomic mechanisms and selective patterns

Nodule bacteria (rhizobia), N₂-fixing symbionts of leguminous plants, represent an excellent model to study the fundamental issues of evolutionary biology, including the tradeoff between microevolution, speciation, and macroevolution, which remains poorly understood for free-living organisms. Taxonomically, rhizobia are extremely diverse: they are represented by nearly a dozen families of α-proteobacteria (Rhizobiales) and by some β-proteobacteria. Their genomes are composed of core parts, including house-keeping genes (hkg), and of accessory parts, including symbiotically specialized (sym) genes. In multipartite genomes of evolutionary advanced fast-growing species (Rhizobiaceae), sym genes are clustered on extra-chromosomal replicons (megaplasmids, chromids), facilitating gene transfer in plant-associated microbial communities. In this review, we demonstrate that in rhizobia, microevolution and speciation involve different genomic and ecological mechanisms: the first one is based on the diversification of sym genes occurring under the impacts of host-induced natural selection (including its disruptive, frequency-dependent and group forms); the second one-on the diversification of hkgs under the impacts of unknown factors. By contrast, macroevolution represents the polyphyletic origin of super-species taxa, which are dependent on the transfer of sym genes from rhizobia to various soil-borne bacteria. Since the expression of newly acquired sym genes on foreign genomic backgrounds is usually restricted, conversion of resulted recombinants into the novel rhizobia species involves post-transfer genetic changes. They are presumably supported by host-induced selective processes resulting in the sequential derepression of nod genes responsible for nodulation and of nif/fix genes responsible for symbiotic N₂ fixation.

A Novel Module Promotes Horizontal Gene Transfer in Azorhizobium caulinodans ORS571

Azorhizobium caulinodans ORS571 contains an 87.6 kb integrative and conjugative element (ICE^Ac) that conjugatively transfers symbiosis genes to other rhizobia. Many hypothetical redundant gene fragments (rgfs) are abundant in ICE^Ac, but their potential function in horizontal gene transfer (HGT) is unknown. Molecular biological methods were employed to delete hypothetical rgfs, expecting to acquire a minimal ICE^Ac and consider non-functional rgfs as editable regions for inserting genes related to new symbiotic functions. We determined the significance of rgf4 in HGT and identified the physiological function of genes designated rihF1a (AZC_3879), rihF1b (AZC_RS26200), and rihR (AZC_3881). In-frame deletion and complementation assays revealed that rihF1a and rihF1b work as a unit (rihF1) that positively affects HGT frequency. The EMSA assay and lacZ-based reporter system showed that the XRE-family protein RihR is not a regulator of rihF1 but promotes the expression of the integrase (intC) that has been reported to be upregulated by the LysR-family protein, AhaR, through sensing host’s flavonoid. Overall, a conservative module containing rihF1 and rihR was characterized, eliminating the size of ICE^Ac by 18.5%. We propose the feasibility of constructing a minimal ICE^Ac element to facilitate the exchange of new genetic components essential for symbiosis or other metabolic functions between soil bacteria.

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

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

MGEs as the MVPs of Partner Quality Variation in Legume-Rhizobium Symbiosis

Despite decades of research, we are only just beginning to understand the forces maintaining variation in the nitrogen-fixing symbiosis between rhizobial bacteria and leguminous plants. In their recent work, Alexandra Weisberg and colleagues use genomics to document the breadth of mobile element diversity that carries the symbiosis genes of Bradyrhizobium in natural populations. Studying rhizobia from the perspective of their mobile genetic elements, which have their own transmission modes and fitness interests, reveals novel mechanisms for the generation and maintenance of diversity in natural populations of these ecologically and economically important mutualisms.

Ecosystem consequences of introducing plant growth promoting rhizobacteria to managed systems and potential legacy effects

The rapidly growing industry of crop biostimulants leverages the application of plant growth promoting rhizobacteria (PGPR) to promote plant growth and health. However, introducing nonnative rhizobacteria may impact other aspects of ecosystem functioning and have legacy effects; these potential consequences are largely unexplored. Nontarget consequences of PGPR may include changes in resident microbiomes, nutrient cycling, pollinator services, functioning of other herbivores, disease suppression, and organic matter persistence. Importantly, we lack knowledge of whether these ecosystem effects may manifest in adjacent ecosystems. The introduced PGPR can leave a functional legacy whether they persist in the community or not. Legacy effects include shifts in resident microbiomes and their temporal dynamics, horizontal transfer of genes from the PGPR to resident taxa, and changes in resident functional groups and interaction networks. Ecosystem functions may be affected by legacies PGPR leave following niche construction, such as when PGPR alter soil pH that in turn alters biogeochemical cycling rates. Here, we highlight new research directions to elucidate how introduced PGPR impact resident microbiomes and ecosystem functions and their capacity for legacy effects.

Prevalence of an Insect-Associated Genomic Region in Environmentally Acquired Burkholderiaceae Symbionts

Microbial symbionts are critical for the development and survival of many eukaryotes. Recent research suggests that the genes enabling these relationships can be localized in horizontally transferred regions of microbial genomes termed "symbiotic islands." Recently, a putative symbiotic island was found that may facilitate symbioses between true bugs and numerous Burkholderia species, based on analysis of five Burkholderia symbionts. We expanded on this work by exploring the putative island's prevalence, origin, and association with colonization across the bacterial family Burkholderiaceae. We performed a broad comparative analysis of 229 Burkholderiaceae genomes, including 8 new genomes of insect- or soil-associated Burkholderia sequenced for this study. We detected the region in 23% of the genomes; these were located solely within two Burkholderia clades. Our analyses suggested that the contiguous region arose at the common ancestor of plant- and insect-associated Burkholderia clades, but the genes themselves are ancestral. Although the region was initially discovered on plasmids and we did detect two likely instances of horizontal transfer within Burkholderia, we found that the region is almost always localized to a chromosome and does not possess any of the mobility elements that typify genomic islands. Finally, to attempt to deduce the region's function, we combined our data with information on several strains' abilities to colonize the insect's symbiotic organ. Although the region was associated with improved colonization of the host, this relationship was confounded with, and likely driven by, Burkholderia clade membership. These findings advance our understanding of the genomic underpinnings of a widespread insect-microbe symbiosis. IMPORTANCE Many plants and animals form intricate associations with bacteria. These pairings can be mediated by genomic islands, contiguous regions containing numerous genes with cohesive functionality. Pathogen-associated islands are well described, but recent evidence suggests that mutualistic islands, which benefit both host and symbiont, may also be common. Recently, a putative symbiosis island was found in Burkholderia symbionts of insects. We determined that this genomic region is located in only two clades of Burkholderia (the plant- and insect-associated species) and that although it has undergone horizontal transfer, it is most likely a symbiosis-associated region rather than a true island. This region is associated with improved host colonization, although this is may be due to specific Burkholderia clades' abilities to colonize rather than presence of the region. By studying the genomic basis of the insect-Burkholderia symbiosis, we can better understand how mutualisms evolve in animals.

Rhizobial diversity is associated with inoculation history at a two-continent scale

A total of 120 Mesorhizobium strains collected from the central dry zone of Myanmar were analyzed in a pot experiment to evaluate nodulation and symbiotic effectiveness (SE%) in chickpea plants. Phylogenetic analyses revealed all strains belonged to the genus Mesorhizobium according to 16–23S rDNA IGS and the majority of chickpea nodulating rhizobia in Myanmar soils were most closely related to M. gobiense, M. muleiense, M. silamurunense, M. tamadayense and M. temperatum. Around two-thirds of the Myanmar strains (68%) were most closely related to Indian strain IC-2058 (CA-181), which is also most closely related to M. gobiense. There were no strains that were closely related to the cognate rhizobial species to nodulate chickpea: M. ciceri and M. mediterraneum. Strains with diverse 16S–23S rDNA IGS shared similar nodC and nifH gene sequences with chickpea symbionts. Detailed sequence analysis of nodC and nifH found that the strains in Myanmar were somewhat divergent from the group including M. ciceri and were more closely related to M. muleiense and IC-2058. A cross-continent analysis between strains isolated in Australia compared with Myanmar found that there was little overlap in species, where Australian soils were dominated with M. ciceri, M. temperatum and M. huakuii. The only co-occurring species found in both Myanmar and Australia were M. tamadayense and M. silumurunense. Continued inoculation with CC1192 may have reduced diversity of chickpea strains in Australian soils. Isolated strains in Australian and Myanmar had similar adaptive traits, which in some cases were also phylogenetically related. The genetic discrepancy between chickpea nodulating strains in Australia and Myanmar is not only due to inoculation history but to adaptation to soil conditions and crop management over a long period, and there has been virtually no loss of symbiotic efficiency over this time in strains isolated from soils in Myanmar.

Why are rhizobial symbiosis genes mobile?

Rhizobia are one of the most important and best studied groups of bacterial symbionts. They are defined by their ability to establish nitrogen-fixing intracellular infections within plant hosts. One surprising feature of this symbiosis is that the bacterial genes required for this complex trait are not fixed within the chromosome, but are encoded on mobile genetic elements (MGEs), namely plasmids or integrative and conjugative elements. Evidence suggests that many of these elements are actively mobilizing within rhizobial populations, suggesting that regular symbiosis gene transfer is part of the ecology of rhizobial symbionts. At first glance, this is counterintuitive. The symbiosis trait is highly complex, multipartite and tightly coevolved with the legume hosts, while transfer of genes can be costly and disrupt coadaptation between the chromosome and the symbiosis genes. However, horizontal gene transfer is a process driven not only by the interests of the host bacterium, but also, and perhaps predominantly, by the interests of the MGEs that facilitate it. Thus understanding the role of horizontal gene transfer in the rhizobium-legume symbiosis requires a 'mobile genetic element's-eye view' on the ecology and evolution of this important symbiosis. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Evolution of rhizobial symbiosis islands through insertion sequence-mediated deletion and duplication

Symbiosis between organisms influences their evolution via adaptive changes in genome architectures. Immunity of soybean carrying the Rj2 allele is triggered by NopP (type III secretion system [T3SS]-dependent effector), encoded by symbiosis island A (SymA) in B. diazoefficiens USDA122. This immunity was overcome by many mutants with large SymA deletions that encompassed T3SS (rhc) and N2 fixation (nif) genes and were bounded by insertion sequence (IS) copies in direct orientation, indicating homologous recombination between ISs. Similar deletion events were observed in B. diazoefficiens USDA110 and B. japonicum J5. When we cultured a USDA122 strain with a marker gene sacB inserted into the rhc gene cluster, most sucrose-resistant mutants had deletions in nif/rhc gene clusters, similar to the mutants above. Some deletion mutants were unique to the sacB system and showed lower competitive nodulation capability, indicating that IS-mediated deletions occurred during free-living growth and the host plants selected the mutants. Among 63 natural bradyrhizobial isolates, 2 possessed long duplications (261-357 kb) harboring nif/rhc gene clusters between IS copies in direct orientation via homologous recombination. Therefore, the structures of symbiosis islands are in a state of flux via IS-mediated duplications and deletions during rhizobial saprophytic growth, and host plants select mutualistic variants from the resultant pools of rhizobial populations. Our results demonstrate that homologous recombination between direct IS copies provides a natural mechanism generating deletions and duplications on symbiosis islands.

An epigenetic switch activates bacterial quorum sensing and horizontal transfer of an integrative and conjugative element

Horizontal transfer of the integrative and conjugative element ICEMlSym^R7A converts non-symbiotic Mesorhizobium spp. into nitrogen-fixing legume symbionts. Here, we discover subpopulations of Mesorhizobium japonicum R7A become epigenetically primed for quorum-sensing (QS) and QS-activated horizontal transfer. Isolated populations in this state termed R7A* maintained these phenotypes in laboratory culture but did not transfer the R7A* state to recipients of ICEMlSym^R7A following conjugation. We previously demonstrated ICEMlSym^R7A transfer and QS are repressed by the antiactivator QseM in R7A populations and that the adjacently-coded DNA-binding protein QseC represses qseM transcription. Here RNA-sequencing revealed qseM expression was repressed in R7A* cells and that RNA antisense to qseC was abundant in R7A but not R7A*. Deletion of the antisense-qseC promoter converted cells into an R7A*-like state. An adjacently coded QseC2 protein bound two operator sites and repressed antisense-qseC transcription. Plasmid overexpression of QseC2 stimulated the R7A* state, which persisted following curing of this plasmid. The epigenetic maintenance of the R7A* state required ICEMlSym^R7A-encoded copies of both qseC and qseC2. Therefore, QseC and QseC2, together with their DNA-binding sites and overlapping promoters, form a stable epigenetic switch that establishes binary control over qseM transcription and primes a subpopulation of R7A cells for QS and horizontal transfer.

Comparative analysis of integrative and conjugative mobile genetic elements in the genus Mesorhizobium

Members of the Mesorhizobium genus are soil bacteria that often form nitrogen-fixing symbioses with legumes. Most characterised Mesorhizobium spp. genomes are ~8 Mb in size and harbour extensive pangenomes including large integrative and conjugative elements (ICEs) carrying genes required for symbiosis (ICESyms). Here, we document and compare the conjugative mobilome of 41 complete Mesorhizobium genomes. We delineated 56 ICEs and 24 integrative and mobilizable elements (IMEs) collectively occupying 16 distinct integration sites, along with 24 plasmids. We also demonstrated horizontal transfer of the largest (853,775 bp) documented ICE, the tripartite ICEMspSym^AA22. The conjugation systems of all identified ICEs and several plasmids were related to those of the paradigm ICESym ICEMlSym^R7A, with each carrying conserved genes for conjugative pilus formation (trb), excision (rdfS), DNA transfer (rlxS) and regulation (fseA). ICESyms have likely evolved from a common ancestor, despite occupying a variety of distinct integration sites and specifying symbiosis with diverse legumes. We found extensive evidence for recombination between ICEs and particularly ICESyms, which all uniquely lack the conjugation entry-exclusion factor gene trbK. Frequent duplication, replacement and pseudogenization of genes for quorum-sensing-mediated activation and antiactivation of ICE transfer suggests ICE transfer regulation is constantly evolving. Pangenome-wide association analysis of the ICE identified genes potentially involved in symbiosis, rhizosphere colonisation and/or adaptation to distinct legume hosts. In summary, the Mesorhizobium genus has accumulated a large and dynamic pangenome that evolves through ongoing horizontal gene transfer of large conjugative elements related to ICEMlSym^R7A.

Rapid compensatory evolution can rescue low fitness symbioses following partner switching

Partner switching plays an important role in the evolution of symbiosis, enabling local adaptation and recovery from the breakdown of symbiosis. Because of intergenomic epistasis, partner-switched symbioses may possess novel combinations of phenotypes but may also exhibit low fitness due to their lack of recent coevolutionary history. Here, we examine the structure and mechanisms of intergenomic epistasis in the Paramecium-Chlorella symbiosis and test whether compensatory evolution can rescue initially low fitness partner-switched symbioses. Using partner-switch experiments coupled with metabolomics, we show evidence for intergenomic epistasis wherein low fitness is associated with elevated symbiont stress responses either in dark or high irradiance environments, potentially owing to mismatched light management traits between the host and symbiont genotypes. Experimental evolution under high light conditions revealed that an initially low fitness partner-switched non-native host-symbiont pairing rapidly adapted, gaining fitness equivalent to the native host-symbiont pairing in less than 50 host generations. Compensatory evolution took two alternative routes: either hosts evolved higher symbiont loads to mitigate for their new algal symbiont's poor performance, or the algal symbionts themselves evolved higher investment in photosynthesis and photoprotective traits to better mitigate light stress. These findings suggest that partner switching combined with rapid compensatory evolution can enable the recovery and local adaptation of symbioses in response to changing environments.

Competition, Nodule Occupancy, and Persistence of Inoculant Strains: Key Factors in the Rhizobium-Legume Symbioses

Biological nitrogen fixation by Rhizobium-legume symbioses represents an environmentally friendly and inexpensive alternative to the use of chemical nitrogen fertilizers in legume crops. Rhizobial inoculants, applied frequently as biofertilizers, play an important role in sustainable agriculture. However, inoculants often fail to compete for nodule occupancy against native rhizobia with inferior nitrogen-fixing abilities, resulting in low yields. Strains with excellent performance under controlled conditions are typically selected as inoculants, but the rates of nodule occupancy compared to native strains are rarely investigated. Lack of persistence in the field after agricultural cycles, usually due to the transfer of symbiotic genes from the inoculant strain to naturalized populations, also limits the suitability of commercial inoculants. When rhizobial inoculants are based on native strains with a high nitrogen fixation ability, they often have superior performance in the field due to their genetic adaptations to the local environment. Therefore, knowledge from laboratory studies assessing competition and understanding how diverse strains of rhizobia behave, together with assays done under field conditions, may allow us to exploit the effectiveness of native populations selected as elite strains and to breed specific host cultivar-rhizobial strain combinations. Here, we review current knowledge at the molecular level on competition for nodulation and the advances in molecular tools for assessing competitiveness. We then describe ongoing approaches for inoculant development based on native strains and emphasize future perspectives and applications using a multidisciplinary approach to ensure optimal performance of both symbiotic partners.

qPCR assay targeting Bradyrhizobium japonicum shows that row spacing and soybean density affects Bradyrhizobium population

The ability for a soybean plant to be efficiently nodulated when grown as a crop is dependent on the number of effective Bradyrhizobium japonicum that can be found in close proximity to the developing seedling shortly after planting. In Manitoba, the growing of soybean as a crop has increased from less than 500 000 acres in 2008 to over 2.3 million acres in 2017. Since the large increase in soybean production is relatively recent, populations of B. japonicum have not yet developed. In response to this, we developed a primer pair that can identify B. japonicum, and be used to determine the titre found in field soil. Their utility was demonstrated by being used to determine whether row spacing of soybean affects B. japonicum populations, as well as to follow B. japonicum populations in a soybean field over the course of a field season. The data show that plant density can affect B. japonicum populations. Moreover, evidence is presented that suggests plant development affects overall B. japonicum populations.

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

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

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

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

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

Evolution of diverse effective N2-fixing microsymbionts of Cicer arietinum following horizontal transfer of the Mesorhizobium ciceri CC1192 symbiosis integrative and conjugative element

Rhizobia are soil bacteria capable of forming N₂-fixing symbioses with legumes, with highly effective strains often selected in agriculture as inoculants to maximize symbiotic N₂ fixation. When rhizobia in the genus Mesorhizobium have been introduced with exotic legumes into farming systems, horizontal transfer of symbiosis Integrative and Conjugative Elements (ICEs) from the inoculant strain to soil bacteria has resulted in the evolution of ineffective N₂-fixing rhizobia that are competitive for nodulation with the target legume. In Australia, Cicer arietinum (chickpea) has been inoculated since the 1970's with Mesorhizobium ciceri sv. ciceri CC1192, a highly effective strain from Israel. Although the full genome sequence of this organism is available, little is known about the mobility of its symbiosis genes and the diversity of cultivated C. arietinum-nodulating organisms. Here, we show the CC1192 genome harbors a 419-kb symbiosis ICE (ICEMcSym¹¹⁹²) and a 648-kb repABC-type plasmid pMC1192 carrying putative fix genes. We sequenced the genomes of 11 C. arietinum nodule isolates from a field site exclusively inoculated with CC1192 and showed they were diverse unrelated Mesorhizobium carrying ICEMcSym¹¹⁹², indicating they had acquired the ICE by environmental transfer. No exconjugants harboured pMc1192 and the plasmid was not essential for N₂ fixation in CC1192. Laboratory conjugation experiments confirmed ICEMcSym¹¹⁹² is mobile, integrating site-specifically within the 3' end of one of the four ser-tRNA genes in the R7ANS recipient genome. Strikingly, all ICEMcSym¹¹⁹² exconjugants were as efficient at fixing N₂ with C. arietinum as CC1192, demonstrating ICE transfer does not necessarily yield ineffective microsymbionts as previously observed.Importance Symbiotic N₂ fixation is a key component of sustainable agriculture and in many parts of the world legumes are inoculated with highly efficient strains of rhizobia to maximise fixed N₂ inputs into farming systems. Symbiosis genes for Mesorhizobium spp. are often encoded chromosomally within mobile gene clusters called Integrative and Conjugative Elements or ICEs. In Australia, where all agricultural legumes and their rhizobia are exotic, horizontal transfer of ICEs from inoculant Mesorhizobium strains to native rhizobia has led to the evolution of inefficient strains that outcompete the original inoculant, with the potential to render it ineffective. However, the commercial inoculant strain for Cicer arietinum (chickpea), M. ciceri CC1192, has a mobile symbiosis ICE (ICEMcSym¹¹⁹²) which can support high rates of N₂ fixation following either environmental or laboratory transfer into diverse Mesorhizobium backgrounds, demonstrating ICE transfer does not necessarily yield ineffective microsymbionts as previously observed.

Exopolysaccharide Characterization of Rhizobium favelukesii LPU83 and Its Role in the Symbiosis With Alfalfa

One of the greatest inputs of available nitrogen into the biosphere occurs through the biological N2-fixation to ammonium as result of the symbiosis between rhizobia and leguminous plants. These interactions allow increased crop yields on nitrogen-poor soils. Exopolysaccharides (EPS) are key components for the establishment of an effective symbiosis between alfalfa and Ensifer meliloti, as bacteria that lack EPS are unable to infect the host plants. Rhizobium favelukesii LPU83 is an acid-tolerant rhizobia strain capable of nodulating alfalfa but inefficient to fix nitrogen. Aiming to identify the molecular determinants that allow R. favelukesii to infect plants, we studied its EPS biosynthesis. LPU83 produces an EPS I identical to the one present in E. meliloti, but the organization of the genes involved in its synthesis is different. The main gene cluster needed for the synthesis of EPS I in E. meliloti, is split into three different sections in R. favelukesii, which probably arose by a recent event of horizontal gene transfer. A R. favelukesii strain devoided of all the genes needed for the synthesis of EPS I is still able to infect and nodulate alfalfa, suggesting that attention should be directed to other molecules involved in the development of the symbiosis.

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.

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.

A bacteriophage infecting Mesorhizobium species has a prolate capsid and shows similarities to a family of Caulobacter crescentus phages

Mesorhizobium phage vB_MloS_Cp1R7A-A1 was isolated from soil planted with chickpea in Saskatchewan. It is dissimilar in sequence and morphology to previously described rhizobiophages. It is a B3 morphotype virus with a distinct prolate capsid and belongs to the tailed phage family Siphoviridae. Its genome has a GC content of 60.3% and 238 predicted genes. Putative functions were predicted for 57 genes, which include 27 tRNA genes with anticodons corresponding to 18 amino acids. This represents the highest number of tRNA genes reported yet in a rhizobiophage. The gene arrangement shows a partially modular organization. Most of the structural genes are found in one module, whereas tRNA genes are in another. Genes for replication, recombination, and nucleotide metabolism form the third module. The arrangement of the replication module resembles the replication module of Enterobacteria phage T5, raising the possibility that it uses a recombination-based replication mechanism, but there is also a suggestion that a T7-like replication mechanism could be used. Phage termini appear to be long direct repeats of just over 12 kb in length. Phylogenetic analysis revealed that Cp1R7A-A1 is more closely related to PhiCbK-like Caulobacter phages and other B3 morphotype phages than to other rhizobiophages sequenced thus far.

Rhizobium leguminosarum bv. trifolii NodD2 Enhances Competitive Nodule Colonization in the Clover-Rhizobium Symbiosis

Establishment of the symbiotic relationship that develops between rhizobia and their legume hosts is contingent upon an interkingdom signal exchange. In response to host legume flavonoids, NodD proteins from compatible rhizobia activate expression of nodulation genes that produce lipochitin oligosaccharide signaling molecules known as Nod factors. Root nodule formation commences upon legume recognition of compatible Nod factor. Rhizobium leguminosarum was previously considered to contain one copy of nodD; here, we show that some strains of the Trifolium (clover) microsymbiont R. leguminosarum bv. trifolii contain a second copy designated nodD2. nodD2 genes were present in 8 out of 13 strains of R. leguminosarum bv. trifolii, but were absent from the genomes of 16 R. leguminosarum bv. viciae strains. Analysis of single and double nodD1 and nodD2 mutants in R. leguminosarum bv. trifolii strain TA1 revealed that NodD2 was functional and enhanced nodule colonization competitiveness. However, NodD1 showed significantly greater capacity to induce nod gene expression and infection thread formation. Clover species are either annual or perennial and this phenological distinction is rarely crossed by individual R. leguminosarum bv. trifolii microsymbionts for effective symbiosis. Of 13 strains with genome sequences available, 7 of the 8 effective microsymbionts of perennial hosts contained nodD2, whereas the 3 microsymbionts of annual hosts did not. We hypothesize that NodD2 inducer recognition differs from NodD1, and NodD2 functions to enhance competition and effective symbiosis, which may discriminate in favor of perennial hosts.IMPORTANCE Establishment of the rhizobium-legume symbiosis requires a highly specific and complex signal exchange between both participants. Rhizobia perceive legume flavonoid compounds through LysR-type NodD regulators. Often, rhizobia encode multiple copies of nodD, which is one determinant of host specificity. In some species of rhizobia, the presence of multiple copies of NodD extends their symbiotic host-range. Here, we identified and characterized a second copy of nodD present in some strains of the clover microsymbiont Rhizobium leguminosarum bv. trifolii. The second nodD gene contributed to the competitive ability of the strain on white clover, an important forage legume. A screen for strains containing nodD2 could be utilized as one criterion to select strains with enhanced competitive ability for use as inoculants for pasture production.

Symbiosis islands of Loteae-nodulating Mesorhizobium comprise three radiating lineages with concordant nod gene complements and nodulation host-range groupings

Mesorhizobium is a genus of soil bacteria, some isolates of which form an endosymbiotic relationship with diverse legumes of the Loteae tribe. The symbiotic genes of these mesorhizobia are generally carried on integrative and conjugative elements termed symbiosis islands (ICESyms). Mesorhizobium strains that nodulate Lotus spp. have been divided into host-range groupings. Group I (GI) strains nodulate L. corniculatus and L. japonicus ecotype Gifu, while group II (GII) strains have a broader host range, which includes L. pedunculatus. To identify the basis of this extended host range, and better understand Mesorhizobium and ICESym genomics, the genomes of eight Mesorhizobium strains were completed using hybrid long- and short-read assembly. Bioinformatic comparison with previously sequenced mesorhizobia genomes indicated host range was not predicted by Mesorhizobium genospecies but rather by the evolutionary relationship between ICESym symbiotic regions. Three radiating lineages of Loteae ICESyms were identified on this basis, which correlate with Lotus spp. host-range grouping and have lineage-specific nod gene complements. Pangenomic analysis of the completed GI and GII ICESyms identified 155 core genes (on average 30.1 % of a given ICESym). Individual GI or GII ICESyms carried diverse accessory genes with an average of 34.6 % of genes unique to a given ICESym. Identification and comparative analysis of NodD symbiotic regulatory motifs – nod boxes – identified 21 branches across the NodD regulons. Four of these branches were associated with seven genes unique to the five GII ICESyms. The nod boxes preceding the host-range gene nodZ in GI and GII ICESyms were disparate, suggesting regulation of nodZ may differ between GI and GII ICESyms. The broad host-range determinant(s) of GII ICESyms that confer nodulation of L. pedunculatus are likely present amongst the 53 GII-unique genes identified.

Connecting signals and benefits through partner choice in plant-microbe interactions

Stabilizing mechanisms in plant-microbe symbioses are critical to maintaining beneficial functions, with two main classes: host sanctions and partner choice. Sanctions are currently presumed to be more effective and widespread, based on the idea that microbes rapidly evolve cheating while retaining signals matching cooperative strains. However, hosts that effectively discriminate among a pool of compatible symbionts would gain a significant fitness advantage. Using the well-characterized legume-rhizobium symbiosis as a model, we evaluate the evidence for partner choice in the context of the growing field of genomics. Empirical studies that rely upon bacteria varying only in nitrogen-fixation ability ignore host-symbiont signaling and frequently conclude that partner choice is not a robust stabilizing mechanism. Here, we argue that partner choice is an overlooked mechanism of mutualism stability and emphasize that plants need not use the microbial services provided a priori to discriminate among suitable partners. Additionally, we present a model that shows that partner choice signaling increases symbiont and host fitness in the absence of sanctions. Finally, we call for a renewed focus on elucidating the signaling mechanisms that are critical to partner choice while further aiming to understand their evolutionary dynamics in nature.

Novel putative rhizobial species with different symbiovars nodulate Lotus creticus and their differential preference to distinctive soil properties

Phylogenetically diverse rhizobial strains endemic to Tunisia were isolated from symbiotic nodules of Lotus creticus, growing in different arid extremophile geographical regions of Tunisia, and speciated using multiloci-phylogenetic analysis as Neorhizobium huautlense (LCK33, LCK35, LCO42 and LCO49), Ensifer numidicus (LCD22, LCD25, LCK22 and LCK25), Ensifer meliloti (LCK8, LCK9 and LCK12) and Mesorhizobium camelthorni (LCD11, LCD13, LCD31 and LCD33). In addition, phylogenetic analyses revealed eight additional strains with previously undescribed chromosomal lineages within the genera Ensifer (LCF5, LCF6 and LCF8),Rhizobium (LCF11, LCF12 and LCF14) and Mesorhizobium (LCF16 and LCF19). Analysis using the nodC gene identified five symbiovar groups, four of which were already known. The remaining group composed of two strains (LCD11 and LCD33) represented a new symbiovar of Mesorhizobium camelthorni, which we propose designating as sv. hedysari. Interestingly, we report that soil properties drive and structure the symbiosis of L. creticus and its rhizobia.

Young JP, Andersen SU. Symbiosis genes show a unique pattern of introgression and selection within a Rhizobium leguminosarum species complex

Rhizobia supply legumes with fixed nitrogen using a set of symbiosis genes. These can cross rhizobium species boundaries, but it is unclear how many other genes show similar mobility. Here, we investigate inter-species introgression using de novo assembly of 196 Rhizobium leguminosarum sv. trifolii genomes. The 196 strains constituted a five-species complex, and we calculated introgression scores based on gene-tree traversal to identify 171 genes that frequently cross species boundaries. Rather than relying on the gene order of a single reference strain, we clustered the introgressing genes into four blocks based on population structure-corrected linkage disequilibrium patterns. The two largest blocks comprised 125 genes and included the symbiosis genes, a smaller block contained 43 mainly chromosomal genes, and the last block consisted of three genes with variable genomic location. All introgression events were likely mediated by conjugation, but only the genes in the symbiosis linkage blocks displayed overrepresentation of distinct, high-frequency haplotypes. The three genes in the last block were core genes essential for symbiosis that had, in some cases, been mobilized on symbiosis plasmids. Inter-species introgression is thus not limited to symbiosis genes and plasmids, but other cases are infrequent and show distinct selection signatures.

Experimental Evolution of Legume Symbionts: What Have We Learnt?

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

Genetically diverse lentil- and faba bean-nodulating rhizobia are present in soils across Central and Southern Ethiopia

In total 196 bacterial isolates were obtained from root nodules of lentil (Lens culinaris) and faba bean (Vicia faba) grown on soil samples collected from 10 different sites in central and southern parts of Ethiopia. All isolates were identified as members of the genus Rhizobium by using recA gene sequence analysis. In the recA phylogenetic tree 195 rhizobial strains were classified into nine genospecies. The phylogeny of symbiotic genes nodC and nifH revealed five and six distinct groups respectively, largely dominated by symbiovar viciae. A multivariate analysis showed that environmental variables of the sampling sites considered in this study had more effect on the distribution and composition of the genospecies than the host legumes of the strains. Twenty representative strains, selected based on their isolation site, host plant and nodC group, were able to nodulate all lentil, faba bean, field pea (Pisum abyssinicum) and grass pea (Lathyrus sativus) plants in a greenhouse test in axenic conditions. The majority of the rhizobial strains were effective nitrogen-fixing symbionts for all tested legumes, indicating their potential to serve as broad host-range inoculants in agriculture. The present work suggests the presence of taxonomically and symbiotically diverse rhizobial species for legumes in the Viciae tribe in Ethiopia.

Recurrent mutualism breakdown events in a legume rhizobia metapopulation

Bacterial mutualists generate major fitness benefits for eukaryotes, reshaping the host phenotype and its interactions with the environment. Yet, microbial mutualist populations are predicted to generate mutants that defect from providing costly services to hosts while maintaining the capacity to exploit host resources. Here, we examined the mutualist service of symbiotic nitrogen fixation in a metapopulation of root-nodulating Bradyrhizobium spp. that associate with the native legume Acmispon strigosus. We quantified mutualism traits of 85 Bradyrhizobium isolates gathered from a 700 km transect in California spanning 10 sampled A. strigosus populations. We clonally inoculated each Bradyrhizobium isolate onto A. strigosus hosts and quantified nodulation capacity and net effects of infection, including host growth and isotopic nitrogen concentration. Six Bradyrhizobium isolates from five populations were categorized as ineffective because they formed nodules but did not enhance host growth via nitrogen fixation. Six additional isolates from three populations failed to form root nodules. Phylogenetic reconstruction inferred two types of mutualism breakdown, including three to four independent losses of effectiveness and five losses of nodulation capacity on A. strigosus. The evolutionary and genomic drivers of these mutualism breakdown events remain poorly understood.

Novel putative Mesorhizobium and Ensifer genomospecies together with a novel symbiovar psoraleae nodulate legumes of agronomic interest grown in Tunisia

Forty rhizobial strains were isolated from Lotus creticus, L. pusillus and Bituminaria bituminosa endemic to Tunisia, and they belonged to the Mesorhizobium and Ensifer genera based on 16S rDNA sequence phylogeny. According to the concatenated recA and glnII sequence-based phylogeny, four Bituminaria isolates Pb5, Pb12, Pb8 and Pb17 formed a monophyletic group with Mesorhizobium chacoense ICMP14587^T, whereas four other strains Pb1, Pb6, Pb13 and Pb15 formed two separate lineages within the Ensifer genus. Among the L. pusillus strains, Lpus9 and Lpus10 showed a 96% identical nucleotide with Ensifer meliloti CCBAU83493^T; whereas six other strains could belong to previously undescribed Mesorhizobium and Ensifer species. For L. creticus strains, Lcus37, Lcus39 and Lcus44 showed 98% sequence identity with Ensifer aridi JNVU TP6, and Lcus42 shared a 96% identical nucleotide with Ensifer meliloti CCBAU83493^T; whereas another four strains were divergent from all the described Ensifer and Mesorhizobium species. The analysis of the nodC gene-based phylogeny identified four symbiovar groups; Mesorhizobium sp. sv. anthyllidis (Lpus3 and Lpus11 from L. pusillus, Lcus43 from L. creticus), Ensifer medicae sv. meliloti (four strains from L. creticus and two strains from L. pusillus), E. meliloti sv. meliloti (four from L. creticus, four from L. pusillus and four from B. bituminosa). In addition, four B. bituminosa strains (Pb5, Pb8, Pb12, and Pb17) displayed a distinctive nodC sequence distant from those of other symbiovars described to date. According to their symbiotic gene sequences and host range, the B. bituminosa symbionts (Pb5, Pb8, Pb12 and Pb17) would represent a new symbiovar of M. chacoense for which sv. psoraleae is proposed.

noeM, a New Nodulation Gene Involved in the Biosynthesis of Nod Factors with an Open-Chain Oxidized Terminal Residue and in the Symbiosis with Mimosa pudica

The β-rhizobium Cupriavidus taiwanensis is a nitrogen-fixing symbiont of Mimosa pudica. Nod factors produced by this species were previously found to be pentameric chitin-oligomers carrying common C18:1 or C16:0 fatty acyl chains, N-methylated and C-6 carbamoylated on the nonreducing terminal N-acetylglucosamine and sulfated on the reducing terminal residue. Here, we report that, in addition, C. taiwanensis LMG19424 produces molecules where the reducing sugar is open and oxidized. We identified a novel nodulation gene located on the symbiotic plasmid pRalta, called noeM, which is involved in this atypical Nod factor structure. noeM encodes a transmembrane protein bearing a fatty acid hydroxylase domain. This gene is expressed during symbiosis with M. pudica and requires NodD and luteolin for optimal expression. The closest noeM homologs formed a separate phylogenetic clade containing rhizobial genes only, which are located on symbiosis plasmids downstream from a nod box. Corresponding proteins, referred to as NoeM, may have specialized in symbiosis via the connection to the nodulation pathway and the spread in rhizobia. noeM was mostly found in isolates of the Mimoseae tribe, and specifically detected in all tested strains able to nodulate M. pudica. A noeM deletion mutant of C. taiwanensis was affected for the nodulation of M. pudica, confirming the role of noeM in the symbiosis with this legume.

Exploring Genetic Diversity and Signatures of Horizontal Gene Transfer in Nodule Bacteria Associated with Lotus japonicus in Natural Environments

To investigate the genetic diversity and understand the process of horizontal gene transfer (HGT) in nodule bacteria associated with Lotus japonicus, we analyzed sequences of three housekeeping and five symbiotic genes using samples from a geographically wide range in Japan. A phylogenetic analysis of the housekeeping genes indicated that L. japonicus in natural environments was associated with diverse lineages of Mesorhizobium spp., whereas the sequences of symbiotic genes were highly similar between strains, resulting in remarkably low nucleotide diversity at both synonymous and nonsynonymous sites. Guanine-cytosine content values were lower in symbiotic genes, and relative frequencies of recombination between symbiotic genes were also lower than those between housekeeping genes. An analysis of molecular variance showed significant genetic differentiation among populations in both symbiotic and housekeeping genes. These results confirm that the Mesorhizobium genes required for symbiosis with L. japonicus behave as a genomic island (i.e., a symbiosis island) and suggest that this island has spread into diverse genomic backgrounds of Mesorhizobium via HGT events in natural environments. Furthermore, our data compilation revealed that the genetic diversity of symbiotic genes in L. japonicus-associated symbionts was among the lowest compared with reports of other species, which may be related to the recent population expansion proposed in Japanese populations of L. japonicus.

Conjugative transfer between Rhizobium etli endosymbionts inside the root nodule

Since the discovery that biological nitrogen fixation ensues in nodules resulting from the interaction of rhizobia with legumes, nodules were thought to be exclusive for hosting nitrogen-fixing and plant growth promoting bacteria. In this work, we uncover a novel function of nodules, as a niche permissive to acquisition of plasmids via conjugative transfer. We used Rhizobium etli CFN42, which nodulates Phaseolus vulgaris. The genome of R. etli CFN42 contains a chromosome and six plasmids. pRet42a is a conjugative plasmid regulated by Quorum-Sensing (QS), and pRet42d is the symbiotic plasmid. Here, using confocal microscopy and flow cytometry, we show that pRet42a transfers on the root's surface, and unexpectedly, inside the nodules. Conjugation still took place inside nodules, even when it was restricted on the plant surface by placing the QS traI regulator under the promoter of the nitrogenase gene, which is only expressed inside the nodules, or by inhibiting the QS transcriptional induction of transfer genes with a traM antiactivator on an unstable vector maintained on the plant surface and lost inside the nodules. These results conclusively confirm the occurrence of conjugation in these structures, defining them as a protected environment for bacterial diversification.

Legacy of prior host and soil selection on rhizobial fitness in planta

Measuring selection acting on microbial populations in natural or even seminatural environments is challenging because many microbial populations experience variable selection. The majority of rhizobial bacteria are found in the soil. However, they also live symbiotically inside nodules of legume hosts and each nodule can release thousands of daughter cells back into the soil. We tested how past selection (i.e., legacies) by two plant genotypes and by the soil alone affected selection and genetic diversity within a population of 101 strains of Ensifer meliloti. We also identified allelic variants most strongly associated with soil- and host-dependent fitness. In addition to imposing direct selection on rhizobia populations, soil and host environments had lasting effects across host generations. Host presence and genotype during the legacy period explained 22% and 12% of the variance in the strain composition of nodule communities in the second cohort, respectively. Although strains with high host fitness in the legacy cohort tended to be enriched in the second cohort, the diversity of the strain community was greater when the second cohort was preceded by host rather than soil legacies. Our results indicate the potential importance of soil selection driving the evolution of these plant-associated microbes.

Global-level population genomics reveals differential effects of geography and phylogeny on horizontal gene transfer in soil bacteria

Although microorganisms are known to dominate Earth's biospheres and drive biogeochemical cycling, little is known about the geographic distributions of microbial populations or the environmental factors that pattern those distributions. We used a global-level hierarchical sampling scheme to comprehensively characterize the evolutionary relationships and distributional limitations of the nitrogen-fixing bacterial symbionts of the crop chickpea, generating 1,027 draft whole-genome sequences at the level of bacterial populations, including 14 high-quality PacBio genomes from a phylogenetically representative subset. We find that diverse Mesorhizobium taxa perform symbiosis with chickpea and have largely overlapping global distributions. However, sampled locations cluster based on the phylogenetic diversity of Mesorhizobium populations, and diversity clusters correspond to edaphic and environmental factors, primarily soil type and latitude. Despite long-standing evolutionary divergence and geographic isolation, the diverse taxa observed to nodulate chickpea share a set of integrative conjugative elements (ICEs) that encode the major functions of the symbiosis. This symbiosis ICE takes 2 forms in the bacterial chromosome-tripartite and monopartite-with tripartite ICEs confined to a broadly distributed superspecies clade. The pairwise evolutionary relatedness of these elements is controlled as much by geographic distance as by the evolutionary relatedness of the background genome. In contrast, diversity in the broader gene content of Mesorhizobium genomes follows a tight linear relationship with core genome phylogenetic distance, with little detectable effect of geography. These results illustrate how geography and demography can operate differentially on the evolution of bacterial genomes and offer useful insights for the development of improved technologies for sustainable agriculture.

New Insight into the Evolution of Symbiotic Genes in Black Locust-Associated Rhizobia

Nitrogen fixation in legumes occurs via symbiosis with rhizobia. This process involves packages of symbiotic genes on mobile genetic elements that are readily transferred within or between rhizobial species, furnishing the recipient with the ability to interact with plant hosts. However, it remains elusive whether plant host migration has played a role in shaping the current distribution of genetic variation in symbiotic genes. Herein, we examined the genetic structure and phylogeographic pattern of symbiotic genes in 286 symbiotic strains of Mesorhizobium nodulating black locust (Robinia pseudoacacia), a cross-continental invasive legume species that is native to North America. We conducted detailed phylogeographic analysis and approximate Bayesian computation to unravel the complex demographic history of five key symbiotic genes. The sequencing results indicate an origin of symbiotic genes in Germany rather than North America. Our findings provide strong evidence of prehistoric lineage splitting and spatial expansion events resulting in multiple radiations of descendent clones from founding sequence types worldwide. Estimates of the timescale of divergence in North American and Chinese subclades suggest that black locust-specific symbiotic genes have been present in these continent many thousands of years before recent migration of plant host. Although numerous crop plants, including legumes, have found their centers of origin as centers of evolution and diversity, the number of legume-specific symbiotic genes with a known geographic origin is limited. This work sheds light on the coevolution of legumes and rhizobia.

Novel large-scale chromosomal transfer in Bacteroides fragilis contributes to its pan-genome and rapid environmental adaptation

Bacteroides fragilis, an important component of the human gastrointestinal microbiota, can cause lethal extra-intestinal infection upon escape from the gastrointestinal tract. We demonstrated transfer and recombination of large chromosomal segments from B. fragilis HMW615, a multidrug resistant clinical isolate, to B. fragilis 638R. In one example, the transfer of a segment of ~435 Kb/356 genes replaced ~413 Kb/326 genes of the B. fragilis 638R chromosome. In addition to transfer of antibiotic resistance genes, these transfers (1) replaced complete divergent polysaccharide biosynthesis loci; (2) replaced DNA inversion-controlled intergenic shufflons (that control expression of genes encoding starch utilization system outer membrane proteins) with more complex, divergent shufflons; and (3) introduced additional intergenic shufflons encoding divergent Type 1 restriction/modification systems. Conjugative transposon-like genes within a transferred segment and within a putative integrative conjugative element (ICE5) ~45 kb downstream from the transferred segment both encode proteins that may be involved in the observed transfer. These data indicate that chromosomal transfer is a driver of antigenic diversity and nutrient adaptation in Bacteroides that (1) contributes to the dissemination of the extensive B. fragilis pan-genome, (2) allows rapid adaptation to a changing environment and (3) can confer pathogenic characteristics to host symbionts.