Population genomics of the symbiotic plasmids of sympatric nitrogen-fixing Rhizobium species associated with Phaseolus vulgaris

How might bacteriophages shape biological invasions?

Invasions by eukaryotes dependent on environmentally acquired bacterial mutualists are often limited by the ability of bacterial partners to survive and establish free-living populations. Focusing on the model legume-rhizobium mutualism, we apply invasion biology hypotheses to explain how bacteriophages can impact the competitiveness of introduced bacterial mutualists. Predicting how phage-bacteria interactions affect invading eukaryotic hosts requires knowing the eco-evolutionary constraints of introduced and native microbial communities, as well as their differences in abundance and diversity. By synthesizing research from invasion biology, as well as bacterial, viral, and community ecology, we create a conceptual framework for understanding and predicting how phages can affect biological invasions through their effects on bacterial mutualists.

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.

Optimistic contributions of plant growth-promoting bacteria for sustainable agriculture and climate stress alleviation

Global climate change is the major cause of abiotic and biotic stresses that have adverse effects on agricultural productivity to an irreversible level, thus threatening to limit gains in production and imperil sustainable agriculture. These climate change-induced abiotic stresses, especially saline, drought, extreme temperature, and so on affect plant morphological, physiological, biochemical, and metabolic characteristics through various pathways and mechanisms, ultimately hindering plant growth, development, and productivity. However, overuse and other inappropriate uses of agrochemicals are not conducive to the protection of natural resources and the environment, thus hampering sustainable agricultural development. With the vigorous development of modern agriculture, the application of plant growth-promoting bacteria (PGPB) can better ensure sustainable agriculture, due to their ability to improve soil properties and confer stress tolerance in plants. This review deciphered the underlying mechanisms of PGPB involved in enhancing plant stress tolerance and performance under various abiotic and biotic stresses. Moreover, the recent advancements in PGPB inoculation techniques, the commercialization of PGPB-based technology and the current applications of PGPB in sustainable agriculture were extensively discussed. Finally, an outlook on the future directions of microbe-aided agriculture was pointed out. Providing insights into plant-PGPB interactions under biotic and abiotic stresses and offering evidence and strategies for PGPB better commercialization and implementation can inspire the development of innovative solutions exploiting PGPB under climatological conditions.

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.

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.

The changing paradigm of rhizobial taxonomy and its systematic growth upto postgenomic technologies

Rhizobia are a diazotrophic group of bacteria that are usually isolated form the nodules in roots, stem of leguminous plants and are able to form nodules in the host plant owing to the presence of symbiotic genes. The rhizobial community is highly diverse, and therefore, the taxonomy and genera-wise classification of rhizobia has been constantly changing since the last three decades. This is mainly due to technical advancements, and shifts in definitions, resulting in a changing paradigm of rhizobia taxonomy. Initially, the taxonomic definitions at the species and sub species level were based on phylogenetic analysis of 16S rRNA sequence, followed by polyphasic approach to have phenotypic, biochemical, and genetic analysis including multilocus sequence analysis. Rhizobia mainly belonging to α- and β-proteobacteria, and recently new additions from γ-proteobacteria had been classified. Nowadays rhizobial taxonomy has been replaced by genome-based taxonomy that allows gaining more insights of genomic characteristics. These omics-technologies provide genome specific information that considers nodulation and symbiotic genes, along with molecular markers as taxonomic traits. Taxonomy based on complete genome sequence (genotaxonomy), average nucleotide identity, is now being considered as primary approach, resulting in an ongoing paradigm shift in rhizobial taxonomy. Also, pairwise whole-genome comparisons, phylogenomic analyses offer correlations between DNA and DNA re-association values that have delineated biologically important species. This review elaborates the present classification and taxonomy of rhizobia, vis-a-vis development of technical advancements, parameters and controversies associated with it, and describe the updated information on evolutionary lineages of rhizobia.

Beyond the ABCs—Discovery of Three New Plasmid Types in Rhodobacterales (RepQ, RepY, RepW)

Copiotrophic marine bacteria of the Roseobacter group (Rhodobacterales, Alphaproteobacteria) are characterized by a multipartite genome organization. We sequenced the genomes of Sulfitobacter indolifex DSM 14862^T and four related plasmid-rich isolates in order to investigate the composition, distribution, and evolution of their extrachromosomal replicons (ECRs). A combination of long-read PacBio and short-read Illumina sequencing was required to establish complete closed genomes that comprised up to twelve ECRs. The ECRs were differentiated in stably evolving chromids and genuine plasmids. Among the chromids, a diagnostic RepABC-8 replicon was detected in four Sulfitobacter species that likely reflects an evolutionary innovation that originated in their common ancestor. Classification of the ECRs showed that the most abundant plasmid system is RepABC, followed by RepA, DnaA-like, and RepB. However, the strains also contained three novel plasmid types that were designated RepQ, RepY, and RepW. We confirmed the functionality of their replicases, investigated the genetic inventory of the mostly cryptic plasmids, and retraced their evolutionary origin. Remarkably, the RepY plasmid of S. pontiacus DSM 110277 is the first high copy-number plasmid discovered in Rhodobacterales.

Spatial patterns in phage-Rhizobium coevolutionary interactions across regions of common bean domestication

Bacteriophages play significant roles in the composition, diversity, and evolution of bacterial communities. Despite their importance, it remains unclear how phage diversity and phage-host interactions are spatially structured. Local adaptation may play a key role. Nitrogen-fixing symbiotic bacteria, known as rhizobia, have been shown to locally adapt to domesticated common bean at its Mesoamerican and Andean sites of origin. This may affect phage-rhizobium interactions. However, knowledge about the diversity and coevolution of phages with their respective Rhizobium populations is lacking. Here, through the study of four phage-Rhizobium communities in Mexico and Argentina, we show that both phage and host diversity is spatially structured. Cross-infection experiments demonstrated that phage infection rates were higher overall in sympatric rhizobia than in allopatric rhizobia except for one Argentinean community, indicating phage local adaptation and host maladaptation. Phage-host interactions were shaped by the genetic identity and geographic origin of both the phage and the host. The phages ranged from specialists to generalists, revealing a nested network of interactions. Our results suggest a key role of local adaptation to resident host bacterial communities in shaping the phage genetic and phenotypic composition, following a similar spatial pattern of diversity and coevolution to that in the host.

Comparative genomics reveals high rates of horizontal transfer and strong purifying selection on rhizobial symbiosis genes

Horizontal transfer (HT) alters the repertoire of symbiosis genes in rhizobial genomes and may play an important role in the on-going evolution of the rhizobia-legume symbiosis. To gain insight into the extent of HT of symbiosis genes with different functional roles (nodulation, N-fixation, host benefit and rhizobial fitness), we conducted comparative genomic and selection analyses of the full-genome sequences from 27 rhizobial genomes. We find that symbiosis genes experience high rates of HT among rhizobial lineages but also bear signatures of purifying selection (low Ka : Ks). HT and purifying selection appear to be particularly strong in genes involved in initiating the symbiosis (e.g. nodulation) and in genome-wide association candidates for mediating benefits provided to the host. These patterns are consistent with rhizobia adapting to the host environment through the loss and gain of symbiosis genes, but not with host-imposed positive selection driving divergence of symbiosis genes through recurring bouts of positive selection.

Phylogenomics reveals the basis of adaptation of Pseudorhizobium species to extreme environments and supports a taxonomic revision of the genus

The family Rhizobiaceae includes many genera of soil bacteria, often isolated for their association with plants. Herein, we investigate the genomic diversity of a group of Rhizobium species and unclassified strains isolated from atypical environments, including seawater, rock matrix or polluted soil. Based on whole-genome similarity and core genome phylogeny, we show that this group corresponds to the genus Pseudorhizobium. We thus reclassify Rhizobium halotolerans, R. marinum, R. flavum and R. endolithicum as P. halotolerans sp. nov., P. marinum comb. nov., P. flavum comb. nov. and P. endolithicum comb. nov., respectively, and show that P. pelagicum is a synonym of P. marinum. We also delineate a new chemolithoautotroph species, P. banfieldiae sp. nov., whose type strain is NT-26^T (=DSM 106348^T=CFBP 8663^T). This genome-based classification was supported by a chemotaxonomic comparison, with increasing taxonomic resolution provided by fatty acid, protein and metabolic profiles. In addition, we used a phylogenetic approach to infer scenarios of duplication, horizontal transfer and loss for all genes in the Pseudorhizobium pangenome. We thus identify the key functions associated with the diversification of each species and higher clades, shedding light on the mechanisms of adaptation to their respective ecological niches. Respiratory proteins acquired at the origin of Pseudorhizobium were combined with clade-specific genes to enable different strategies for detoxification and nutrition in harsh, nutrient-poor environments.

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

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

Comprehensive screening of genomic and metagenomic data reveals a large diversity of tetracycline resistance genes

Tetracyclines are broad-spectrum antibiotics used to prevent or treat a variety of bacterial infections. Resistance is often mediated through mobile resistance genes, which encode one of the three main mechanisms: active efflux, ribosomal target protection or enzymatic degradation. In the last few decades, a large number of new tetracycline-resistance genes have been discovered in clinical settings. These genes are hypothesized to originate from environmental and commensal bacteria, but the diversity of tetracycline-resistance determinants that have not yet been mobilized into pathogens is unknown. In this study, we aimed to characterize the potential tetracycline resistome by screening genomic and metagenomic data for novel resistance genes. By using probabilistic models, we predicted 1254 unique putative tetracycline resistance genes, representing 195 gene families (<70 % amino acid sequence identity), whereof 164 families had not been described previously. Out of 17 predicted genes selected for experimental verification, 7 induced a resistance phenotype in an Escherichia coli host. Several of the predicted genes were located on mobile genetic elements or in regions that indicated mobility, suggesting that they easily can be shared between bacteria. Furthermore, phylogenetic analysis indicated several events of horizontal gene transfer between bacterial phyla. Our results also suggested that acquired efflux pumps originate from proteobacterial species, while ribosomal protection genes have been mobilized from Firmicutes and Actinobacteria. This study significantly expands the knowledge of known and putatively novel tetracycline resistance genes, their mobility and evolutionary history. The study also provides insights into the unknown resistome and genes that may be encountered in clinical settings in the future.

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.

Bacillus velezensis 83 a bacterial strain from mango phyllosphere, useful for biological control and plant growth promotion

Bacillus velezensis 83 was isolated from mango tree phyllosphere of orchards located in El Rosario, Sinaloa, México. The assessment of this strain as BCA (biological control agent), as well as PGPB (plant growth-promoting bacteria), were demonstrated through in vivo and in vitro assays. In vivo assays showed that B. velezensis 83 was able to control anthracnose (Kent mangoes) as efficiently as chemical treatment with Captan 50 PH™ or Cupravit hidro™. The inoculation of B. velezensis 83 to the roots of maize seedlings yielded an increase of 12% in height and 45% of root biomass, as compared with uninoculated seedlings. In vitro co-culture assays showed that B. velezensis 83 promoted Arabidopsis thaliana growth (root and shoot biomass) while, under the same experimental conditions, B. velezensis FZB42 (reference strain) had a suppressive effect on plant growth. In order to characterize the isolated strain, the complete genome sequence of B. velezensis 83 is reported. Its circular genome consists of 3,997,902 bp coding to 3949 predicted genes. The assembly and annotation of this genome revealed gene clusters related with plant-bacteria interaction and sporulation, as well as ten secondary metabolites biosynthetic gene clusters implicated in the biological control of phytopathogens. Despite the high genomic identity (> 98%) between B. velezensis 83 and B. velezensis FZB42, they are phenotypically different. Indeed, in vitro production of compounds such as surfactin and bacillomycin D (biocontrol activity) and γ-PGA (biofilm component) is significantly different between both strains.

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.

Genomic insight into the origins and evolution of symbiosis genes in Phaseolus vulgaris microsymbionts

Background

Phaseolus vulgaris (common bean) microsymbionts belonging to the bacterial genera Rhizobium, Bradyrhizobium, and Ensifer (Sinorhizobium) have been isolated across the globe. Individual symbiosis genes (e.g., nodC) of these rhizobia can be different within each genus and among distinct genera. Little information is available about the symbiotic structure of indigenous Rhizobium strains nodulating introduced bean plants or the emergence of a symbiotic ability to associate with bean plants in Bradyrhizobium and Ensifer strains. Here, we sequenced the genomes of 29 representative bean microsymbionts (21 Rhizobium, four Ensifer, and four Bradyrhizobium) and compared them with closely related reference strains to estimate the origins of symbiosis genes among these Chinese bean microsymbionts.

Results

Comparative genomics demonstrated horizontal gene transfer exclusively at the plasmid level, leading to expanded diversity of bean-nodulating Rhizobium strains. Analysis of vertically transferred genes uncovered 191 (out of the 2654) single-copy core genes with phylogenies strictly consistent with the taxonomic status of bacterial species, but none were found on symbiosis plasmids. A common symbiotic region was wholly conserved within the Rhizobium genus yet different from those of the other two genera. A single strain of Ensifer and two Bradyrhizobium strains shared similar gene content with soybean microsymbionts in both chromosomes and symbiotic regions.

Conclusions

The 19 native bean Rhizobium microsymbionts were assigned to four defined species and six putative novel species. The symbiosis genes of R. phaseoli, R. sophoriradicis, and R. esperanzae strains that originated from Mexican bean-nodulating strains were possibly introduced alongside bean seeds. R. anhuiense strains displayed distinct host ranges, indicating transition into bean microsymbionts. Among the six putative novel species exclusive to China, horizontal transfer of symbiosis genes suggested symbiosis with other indigenous legumes and loss of originally symbiotic regions or non-symbionts before the introduction of common bean into China. Genome data for Ensifer and Bradyrhizobium strains indicated symbiotic compatibility between microsymbionts of common bean and other hosts such as soybean.

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.

Exploring the evolutionary dynamics of Rhizobium plasmids through bipartite network analysis

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

Phylogenomic Rhizobium Species Are Structured by a Continuum of Diversity and Genomic Clusters

The bacterial genus Rhizobium comprises diverse symbiotic nitrogen-fixing species associated with the roots of plants in the Leguminosae family. Multiple genomic clusters defined by whole genome comparisons occur within Rhizobium, but their equivalence to species is controversial. In this study we investigated such genomic clusters to ascertain their significance in a species phylogeny context. Phylogenomic inferences based on complete sets of ribosomal proteins and stringent core genome markers revealed the main lineages of Rhizobium. The clades corresponding to R. etli and R. leguminosarum species show several genomic clusters with average genomic nucleotide identities (ANI > 95%), and a continuum of divergent strains, respectively. They were found to be inversely correlated with the genetic distance estimated from concatenated ribosomal proteins. We uncovered evidence of a Rhizobium pangenome that was greatly expanded, both in its chromosomes and plasmids. Despite the variability of extra-chromosomal elements, our genomic comparisons revealed only a few chromid and plasmid families. The presence/absence profile of genes in the complete Rhizobium genomes agreed with the phylogenomic pattern of species divergence. Symbiotic genes were distributed according to the principal phylogenomic Rhizobium clades but did not resolve genome clusters within the clades. We distinguished some types of symbiotic plasmids within Rhizobium that displayed different rates of synonymous nucleotide substitutions in comparison to chromosomal genes. Symbiotic plasmids may have been repeatedly transferred horizontally between strains and species, in the process displacing and substituting pre-existing symbiotic plasmids. In summary, the results indicate that Rhizobium genomic clusters, as defined by whole genomic identities, might be part of a continuous process of evolutionary divergence that includes the core and the extrachromosomal elements leading to species formation.

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.

Nodule bacteria from the cultured legume Phaseolus dumosus (belonging to the Phaseolus vulgaris cross-inoculation group) with common tropici phenotypic characteristics and symbiovar but distinctive phylogenomic position and chromid

Phaseolus dumosus is an endemic species from mountain tops in Mexico that was found in traditional agriculture areas in Veracruz, Mexico. P. dumosus plants were identified by ITS sequences and their nodules were collected from agricultural fields or from trap plant experiments in the laboratory. Bacteria from P. dumosus nodules were identified as belonging to the phaseoli-etli-leguminosarum (PEL) or to the tropici group by 16S rRNA gene sequences. We obtained complete closed genomes from two P. dumosus isolates CCGE531 and CCGE532 that were phylogenetically placed within the tropici group but with a distinctive phylogenomic position and low average nucleotide identity (ANI). CCGE531 and CCGE532 had common phenotypic characteristics with tropici type B rhizobial symbionts. Genome synteny analysis and ANI showed that P. dumosus isolates had different chromids and our analysis suggests that chromids have independently evolved in different lineages of the Rhizobium genus. Finally, we considered that P. dumosus and Phaseolus vulgaris plants belong to the same cross-inoculation group since they have conserved symbiotic affinites for rhizobia.

Dynamic genomic architecture of mutualistic cooperation in a wild population of Mesorhizobium

Research on mutualism seeks to explain how cooperation can be maintained when uncooperative mutants co-occur with cooperative kin. Gains and losses of the gene modules required for cooperation punctuate symbiont phylogenies and drive lifestyle transitions between cooperative symbionts and uncooperative free-living lineages over evolutionary time. Yet whether uncooperative symbionts commonly evolve from within cooperative symbiont populations or from within distantly related lineages with antagonistic or free-living lifestyles (i.e., third-party mutualism exploiters or parasites), remains controversial. We use genomic data to show that genotypes that differ in the presence or absence of large islands of symbiosis genes are common within a single wild recombining population of Mesorhizobium symbionts isolated from host tissues and are an important source of standing heritable variation in cooperation in this population. In a focal population of Mesorhizobium, uncooperative variants that lack a symbiosis island segregate at 16% frequency in nodules, and genome size and symbiosis gene number are positively correlated with cooperation. This finding contrasts with the genomic architecture of variation in cooperation in other symbiont populations isolated from host tissues in which the islands of genes underlying cooperation are ubiquitous and variation in cooperation is primarily driven by allelic substitution and individual gene gain and loss events. Our study demonstrates that uncooperative mutants within mutualist populations can comprise a significant component of genetic variation in nature, providing biological rationale for models and experiments that seek to explain the maintenance of mutualism in the face of non-cooperators.

Horizontal operon transfer, plasmids, and the evolution of photosynthesis in Rhodobacteraceae

The capacity for anoxygenic photosynthesis is scattered throughout the phylogeny of the Proteobacteria. Their photosynthesis genes are typically located in a so-called photosynthesis gene cluster (PGC). It is unclear (i) whether phototrophy is an ancestral trait that was frequently lost or (ii) whether it was acquired later by horizontal gene transfer. We investigated the evolution of phototrophy in 105 genome-sequenced Rhodobacteraceae and provide the first unequivocal evidence for the horizontal transfer of the PGC. The 33 concatenated core genes of the PGC formed a robust phylogenetic tree and the comparison with single-gene trees demonstrated the dominance of joint evolution. The PGC tree is, however, largely incongruent with the species tree and at least seven transfers of the PGC are required to reconcile both phylogenies. The origin of a derived branch containing the PGC of the model organism Rhodobacter capsulatus correlates with a diagnostic gene replacement of pufC by pufX. The PGC is located on plasmids in six of the analyzed genomes and its DnaA-like replication module was discovered at a conserved central position of the PGC. A scenario of plasmid-borne horizontal transfer of the PGC and its reintegration into the chromosome could explain the current distribution of phototrophy in Rhodobacteraceae.

A Novel Strategy for Detecting Recent Horizontal Gene Transfer and Its Application to Rhizobium Strains

Recent horizontal gene transfer (HGT) is crucial for enabling microbes to rapidly adapt to their novel environments without relying upon rare beneficial mutations that arise spontaneously. For several years now, computational approaches have been developed to detect HGT, but they typically lack the sensitivity and ability to detect recent HGT events. Here we introduce a novel strategy, named RecentHGT. The number of genes undergoing recent HGT between two bacterial genomes was estimated by a new algorithm derived from the expectation-maximization algorithm and is based on the theoretical sequence-similarity distribution of orthologous genes. We tested the proposed strategy by applying it to a set of 10 Rhizobium genomes, and detected several large-scale recent HGT events. We also found that our strategy was more sensitive than other available HGT detection methods. These HGT events were mainly mediated by symbiotic plasmids. Our new strategy can provide clear evidence of recent HGT events and thus it brings us closer to the goal of detecting these potentially adaptive evolution processes in rhizobia as well as pathogens.

Comparative analysis of rhizobial chromosomes and plasmids to estimate their evolutionary relationships

In the present study, complete genomic sequences retrieved from 57 rhizobial strains that covered four genera including 11 species were analyzed comprehensively. The four types of replicons: chromosomes, chromids, nonsymbiotic plasmids, and symbiotic plasmids were investigated and compared among these strains. Results showed that co-evolution occurred among these four replicons based on the similarities in average nucleotide identity. High correlation coefficient r values were observed between chromosomes and chromids, as well as between chromosomes and nonsymbiotic plasmids. Chromosomes and symbiotic plasmids showed different phylogenetic topology based on their core genes. Population structure analyses were performed to extrapolate the evolutionary histories of the test strains based on their chromosomal and symbiotic plasmid background. This resulted in seven ancestral types for chromosomal genes and three ancestral types for symbiotic plasmid genes. Rhizobial strains containing chromosome genes with ancestral type E tend to contain symbiotic plasmid genes with ancestral type II, while rhizobial strains containing chromosome genes with ancestral type G tend to contain symbiotic plasmid genes with ancestral type III. Seventeen strains associated with different host plant species which harbored the symbiotic genes with ancestral type I, exhibited high genetic diversity. In addition, Fu's test of the symbiotic plasmid genes with ancestral type III had undergone an expansion event, implying the influence of negative selection on these symbiotic plasmid genes.

Genomic insight into the taxonomy of Rhizobium genospecies that nodulate Phaseolus vulgaris

Due to the wide cultivation of bean (Phaseolus vulgaris L.), rhizobia associated with this plant have been isolated from many different geographical regions. In order to investigate the species diversity of bean rhizobia, comparative genome sequence analysis was performed in the present study for 69 Rhizobium strains mainly isolated from root nodules of bean and clover (Trifolium spp.). Based on genome average nucleotide identity, digital DNA:DNA hybridization, and phylogenetic analysis of 1,458 single-copy core genes, these strains were classified into 28 clusters, consistent with their species definition based on multilocus sequence analysis (MLSA) of atpD, glnII, and recA. The bean rhizobia were found in 16 defined species and nine putative novel species; in addition, 35 strains previously described as Rhizobium etli, Rhizobium phaseoli, Rhizobium vallis, Rhizobium gallicum, Rhizobium leguminosarum and Rhizobium spp. should be renamed. The phylogenetic patterns of symbiotic genes nodC and nifH were highly host-specific and inconsistent with the genomic phylogeny. Multiple symbiovars (sv.) within the Rhizobium species were found as a common feature: sv. phaseoli, sv. trifolii and sv. viciae in Rhizobium anhuiense; sv. phaseoli and sv. mimosae in Rhizobium sophoriradicis/R. etli/Rhizobium sp. III; sv. phaseoli and sv. trifolii in Rhizobium hidalgonense/Rhizobium acidisoli; sv. phaseoli and sv. viciae in R. leguminosarum/Rhizobium sp. IX; sv. trifolii and sv. viciae in Rhizobium laguerreae. Thus, genomic comparison revealed great species diversity in bean rhizobia, corrected the species definition of some previously misnamed strains, and demonstrated the MLSA a valuable and simple method for defining Rhizobium species.

Genetic diversity and symbiotic effectiveness of Phaseolus vulgaris-nodulating rhizobia in Kenya

Phaseolus vulgaris (common bean) was introduced to Kenya several centuries ago but the rhizobia that nodulate it in the country remain poorly characterised. To address this gap in knowledge, 178 isolates recovered from the root nodules of P. vulgaris cultivated in Kenya were genotyped stepwise by the analysis of genomic DNA fingerprints, PCR-RFLP and 16S rRNA, atpD, recA and nodC gene sequences. Results indicated that P. vulgaris in Kenya is nodulated by at least six Rhizobium genospecies, with most of the isolates belonging to Rhizobium phaseoli and a possibly novel Rhizobium species. Infrequently, isolates belonged to Rhizobium paranaense, Rhizobium leucaenae, Rhizobium sophoriradicis and Rhizobium aegyptiacum. Despite considerable core-gene heterogeneity among the isolates, only four nodC gene alleles were observed indicating conservation within this gene. Testing of the capacity of the isolates to fix nitrogen (N₂) in symbiosis with P. vulgaris revealed wide variations in effectiveness, with ten isolates comparable to Rhizobium tropici CIAT 899, a commercial inoculant strain for P. vulgaris. In addition to unveiling effective native rhizobial strains with potential as inoculants in Kenya, this study demonstrated that Kenyan soils harbour diverse P. vulgaris-nodulating rhizobia, some of which formed phylogenetic clusters distinct from known lineages. The native rhizobia differed by site, suggesting that field inoculation of P. vulgaris may need to be locally optimised.

Rhizobia: from saprophytes to endosymbionts

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

The Composite 259-kb Plasmid of Martelella mediterranea DSM 17316T-A Natural Replicon with Functional RepABC Modules from Rhodobacteraceae and Rhizobiaceae

A multipartite genome organization with a chromosome and many extrachromosomal replicons (ECRs) is characteristic for Alphaproteobacteria. The best investigated ECRs of terrestrial rhizobia are the symbiotic plasmids for legume root nodulation and the tumor-inducing (Ti) plasmid of Agrobacterium tumefaciens. RepABC plasmids represent the most abundant alphaproteobacterial replicon type. The currently known homologous replication modules of rhizobia and Rhodobacteraceae are phylogenetically distinct. In this study, we surveyed type-strain genomes from the One Thousand Microbial Genomes (KMG-I) project and identified a roseobacter-specific RepABC-type operon in the draft genome of the marine rhizobium Martelella mediterranea DSM 17316^T. PacBio genome sequencing demonstrated the presence of three circular ECRs with sizes of 593, 259, and 170-kb. The rhodobacteral RepABC module is located together with a rhizobial equivalent on the intermediate sized plasmid pMM259, which likely originated in the fusion of a pre-existing rhizobial ECR with a conjugated roseobacter plasmid. Further evidence for horizontal gene transfer (HGT) is given by the presence of a roseobacter-specific type IV secretion system on the 259-kb plasmid and the rhodobacteracean origin of 62% of the genes on this plasmid. Functionality tests documented that the genuine rhizobial RepABC module from the Martelella 259-kb plasmid is only maintained in A. tumefaciens C58 (Rhizobiaceae) but not in Phaeobacter inhibens DSM 17395 (Rhodobacteraceae). Unexpectedly, the roseobacter-like replication system is functional and stably maintained in both host strains, thus providing evidence for a broader host range than previously proposed. In conclusion, pMM259 is the first example of a natural plasmid that likely mediates genetic exchange between roseobacters and rhizobia.

Interspecies cross-talk between co-cultured Pseudomonas putida and Escherichia coli

Pseudomonas putida and Escherichia coli are ubiquitous microorganisms that can be isolated from soil rhizosphere, the surface of vegetables, fresh waters and wastewaters - environments in which they likely co-exist. Despite this, the potential interactions between these microbes have not been studied in detail. To analyse these interactions, we carried out RNA-seq transcriptomic analysis of these microbes as monocultures and as co-cultures. Our results show that co-culture of these microbes significantly alters transcriptional profiles. The most dramatic transcriptional changes in both microorganisms were involved in central carbon metabolism, as well as adhesion to surfaces and the activation of drug efflux pumps. We also found that acetate production was one of the mechanisms used by E. coli K-12 MG1655 in response to the presence of P. putida DOT-T1E.

Complete Genome Sequences of Eight Rhizobium Symbionts Associated with Common Bean (Phaseolus vulgaris)

We present here the high-quality complete genome sequences of eight strains of Rhizobium-nodulating Phaseolus vulgaris. Comparative analyses showed that some of them belonged to different genomic and evolutionary lineages with common symbiotic properties. Two novel symbiotic plasmids (pSyms) with P. vulgaris specificity are reported here.

Natural rhizobial diversity helps to reveal genes and QTLs associated with biological nitrogen fixation in common bean

Common bean is one of the most important crops for human feed, and the most important legume for direct consumption by millions of people, especially in developing countries. It is a promiscuous host legume in terms of nodulation, able to associate with a broad and diverse range of rhizobia, although the competitiveness for nodulation and the nitrogen fixation capacity of most of these strains is generally low. As a result, common bean is very inefficient for symbiotic nitrogen fixation, and nitrogen has to be supplied with chemical fertilizers. In the last years, symbiotic nitrogen fixation has received increasing attention as a sustainable alternative to nitrogen fertilizers, and also as a more economic and available one in poor countries. Therefore, optimization of nitrogen fixation of bean-rhizobia symbioses and selection of efficient rhizobial strains should be a priority, which begins with the study of the natural diversity of the symbioses and the rhizobial populations associated. Natural rhizobia biodiversity that nodulates common bean may be a source of adaptive alleles acting through phenotypic plasticity. Crosses between accessions differing for nitrogen fixation may combine alleles that never meet in nature. Another way to discover adaptive genes is to use association genetics to identify loci that common bean plants use for enhanced biological nitrogen fixation and, in consequence, for marker assisted selection for genetic improvement of symbiotic nitrogen fixation. In this review, rhizobial biodiversity resources will be discussed, together with what is known about the loci that underlie such genetic variation, and the potential candidate genes that may influence the symbiosis' fitness benefits, thus achieving an optimal nitrogen fixation capacity in order to help reduce reliance on nitrogen fertilizers in common bean.