Pangenome Evolution Reconciles Robustness and Instability of Rhizobial Symbiosis

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.

Host-imposed control mechanisms in legume-rhizobia symbiosis

Legumes are ecologically and economically important plants that contribute to nutrient cycling and agricultural sustainability, features tied to their intimate symbiosis with nitrogen-fixing rhizobia. Rhizobia vary dramatically in quality, ranging from highly growth-promoting to non-beneficial; therefore, legumes must optimize their symbiosis with rhizobia through host mechanisms that select for beneficial rhizobia and limit losses to non-beneficial strains. In this Perspective, we examine the considerable scientific progress made in decoding host control over rhizobia, empirically examining both molecular and cellular mechanisms and their effects on rhizobia symbiosis and its benefits. We consider pre-infection controls, which require the production and detection of precise molecular signals by the legume to attract and select for compatible rhizobia strains. We also discuss post-infection mechanisms that leverage the nodule-level and cell-level compartmentalization of symbionts to enable host control over rhizobia development and proliferation in planta. These layers of host control each contribute to legume fitness by directing host resources towards a narrowing subset of more-beneficial rhizobia.

Comparative genomic analysis of Bradyrhizobium strains with natural variability in the efficiency of nitrogen fixation, competitiveness, and adaptation to stressful edaphoclimatic conditions

Bradyrhizobium is known for fixing atmospheric nitrogen in symbiosis with agronomically important crops. This study focused on two groups of strains, each containing eight natural variants of the parental strains, Bradyrhizobium japonicum SEMIA 586 (=CNPSo 17) or Bradyrhizobium diazoefficiens SEMIA 566 (=CNPSo 10). CNPSo 17 and CNPSo 10 were used as commercial inoculants for soybean crops in Brazil at the beginning of the crop expansion in the southern region in the 1960s-1970s. Variants derived from these parental strains were obtained in the late 1980s through a strain selection program aimed at identifying elite strains adapted to a new cropping frontier in the central-western Cerrado region, with a higher capacity of biological nitrogen fixation (BNF) and competitiveness. Here, we aimed to detect genetic variations possibly related to BNF, competitiveness for nodule occupancy, and adaptation to the stressful conditions of the Brazilian Cerrado soils. High-quality genome assemblies were produced for all strains. The core genome phylogeny revealed that strains of each group are closely related, as confirmed by high average nucleotide identity values. However, variants accumulated divergences resulting from horizontal gene transfer, genomic rearrangements, and nucleotide polymorphisms. The B. japonicum group presented a larger pangenome and a higher number of nucleotide polymorphisms than the B. diazoefficiens group, possibly due to its longer adaptation time to the Cerrado soil. Interestingly, five strains of the B. japonicum group carry two plasmids. The genetic variability found in both groups is discussed considering the observed differences in their BNF capacity, competitiveness for nodule occupancy, and environmental adaptation.IMPORTANCEToday, Brazil is a global leader in the study and use of biological nitrogen fixation with soybean crops. As Brazilian soils are naturally void of soybean-compatible bradyrhizobia, strain selection programs were established, starting with foreign isolates. Selection searched for adaptation to the local edaphoclimatic conditions, higher efficiency of nitrogen fixation, and strong competitiveness for nodule occupancy. We analyzed the genomes of two parental strains of Bradyrhizobium japonicum and Bradyrhizobium diazoefficiens and eight variant strains derived from each parental strain. We detected two plasmids in five strains and several genetic differences that might be related to adaptation to the stressful conditions of the soils of the Brazilian Cerrado biome. We also detected genetic variations in specific regions that may impact symbiotic nitrogen fixation. Our analysis contributes to new insights into the evolution of Bradyrhizobium, and some of the identified differences may be applied as genetic markers to assist strain selection programs.

Programming Dynamic Division of Labor Using Horizontal Gene Transfer

The metabolic engineering of microbes has broad applications, including biomanufacturing, bioprocessing, and environmental remediation. The introduction of a complex, multistep pathway often imposes a substantial metabolic burden on the host cell, restraining the accumulation of productive biomass and limiting pathway efficiency. One strategy to alleviate metabolic burden is the division of labor (DOL) in which different subpopulations carry out different parts of the pathway and work together to convert a substrate into a final product. However, the maintenance of different engineered subpopulations is challenging due to competition and convoluted interstrain population dynamics. Through modeling, we show that dynamic division of labor (DDOL), which we define as the DOL between indiscrete populations capable of dynamic and reversible interchange, can overcome these limitations and enable the robust maintenance of burdensome, multistep pathways. We propose that DDOL can be mediated by horizontal gene transfer (HGT) and use plasmid genomics to uncover evidence that DDOL is a strategy utilized by natural microbial communities. Our work suggests that bioengineers can harness HGT to stabilize synthetic metabolic pathways in microbial communities, enabling the development of robust engineered systems for deployment in a variety of contexts.

Lineage-specific evolution of Aquibium, a close relative of Mesorhizobium, during habitat adaptation

The novel genus Aquibium that lacks nitrogenase was recently reclassified from the Mesorhizobium genus. The genomes of Aquibium species isolated from water were smaller and had higher GC contents than those of Mesorhizobium species. Six Mesorhizobium species lacking nitrogenase were found to exhibit low similarity in the average nucleotide identity values to the other 24 Mesorhizobium species. Therefore, they were classified as the non-N2-fixing Mesorhizobium lineage (N-ML), an evolutionary intermediate species. The results of our phylogenomic analyses and the loss of Rhizobiales-specific fur/mur indicated that Mesorhizobium species may have evolved from Aquibium species through an ecological transition. Halotolerant and alkali-resistant Aquibium and Mesorhizobium microcysteis belonging to N-ML possessed many tripartite ATP-independent periplasmic transporter and sodium/proton antiporter subunits composed of seven genes (mrpABCDEFG). These genes were not present in the N2-fixing Mesorhizobium lineage (ML), suggesting that genes acquired for adaptation to highly saline and alkaline environments were lost during the evolution of ML as the habitat changed to soil. Land-to-water habitat changes in Aquibium species, close relatives of Mesorhizobium species, could have influenced their genomic evolution by the gain and loss of genes. Our study indicated that lineage-specific evolution could have played a significant role in shaping their genome architecture and conferring their ability to thrive in different habitats.IMPORTANCEPhylogenetic analyses revealed that the Aquibium lineage (AL) and non-N2-fixing Mesorhizobium lineage (N-ML) were monophyletically grouped into distinct clusters separate from the N2-fixing Mesorhizobium lineage (ML). The N-ML, an evolutionary intermediate species having characteristics of both ancestral and descendant species, could provide a genomic snapshot of the genetic changes that occur during adaptation. Genomic analyses of AL, N-ML, and ML revealed that changes in the levels of genes related to transporters, chemotaxis, and nitrogen fixation likely reflect adaptations to different environmental conditions. Our study sheds light on the complex and dynamic nature of the evolution of rhizobia in response to changes in their environment and highlights the crucial role of genomic analysis in understanding these processes.

Comparative genomics analysis reveals genetic characteristics and nitrogen fixation profile of Bradyrhizobium

Bradyrhizobium is a genus of nitrogen-fixing bacteria, with some species producing nodules in leguminous plants. Investigations into Bradyrhizobium have recently revealed its substantial genetic resources and agricultural benefits, but a comprehensive survey of its genetic diversity and functional properties is lacking. Using a panel of various strains (N = 278), this study performed a comparative genomics analysis to anticipate genes linked with symbiotic nitrogen fixation. Bradyrhizobium's pan-genome consisted of 84,078 gene families, containing 824 core genes and 42,409 accessory genes. Core genes were mainly involved in crucial cell processes, while accessory genes served diverse functions, including nitrogen fixation and nodulation. Three distinct genetic profiles were identified based on the presence/absence of gene clusters related to nodulation, nitrogen fixation, and secretion systems. Most Bradyrhizobium strains from soil and non-leguminous plants lacked major nif/nod genes and were evolutionarily more closely related. These findings shed light on Bradyrhizobium's genetic features for symbiotic nitrogen fixation.

Diverse mobile genetic elements shaped the evolution of Streptomyces virulence

Mobile genetic elements can innovate bacteria with new traits. In plant pathogenic Streptomyces, frequent and recent acquisition of integrative and conjugative or mobilizable genetic elements is predicted to lead to the emergence of new lineages that gained the capacity to synthesize Thaxtomin, a phytotoxin neccesary for induction of common scab disease on tuber and root crops. Here, we identified components of the Streptomyces-potato pathosystem implicated in virulence and investigated them as a nested and interacting system to reevaluate evolutionary models. We sequenced and analysed genomes of 166 strains isolated from over six decades of sampling primarily from field-grown potatoes. Virulence genes were associated to multiple subtypes of genetic elements differing in mechanisms of transmission and evolutionary histories. Evidence is consistent with few ancient acquisition events followed by recurrent loss or swaps of elements carrying Thaxtomin A-associated genes. Subtypes of another genetic element implicated in virulence are more distributed across Streptomyces. However, neither the subtype classification of genetic elements containing virulence genes nor taxonomic identity was predictive of pathogenicity on potato. Last, findings suggested that phytopathogenic strains are generally endemic to potato fields and some lineages were established by historical spread and further dispersed by few recent transmission events. Results from a hierarchical and system-wide characterization refine our understanding by revealing multiple mechanisms that gene and bacterial dispersion have had on shaping the evolution of a Gram-positive pathogen in agricultural settings.

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.

Competitive interference among rhizobia reduces benefits to hosts

The capacity of beneficial microbes to compete for host infection-and the ability of hosts to discriminate among them-introduces evolutionary conflict that is predicted to destabilize mutualism. We investigated fitness outcomes in associations between legumes and their symbiotic rhizobia to characterize fitness impacts of microbial competition. Diverse Bradyrhizobium strains varying in their capacity to fix nitrogen symbiotically with a common host plant, Acmispon strigosus, were tested in full-factorial coinoculation experiments involving 28 pairwise strain combinations. We analyzed the effects of interstrain competition and host discrimination on symbiotic-interaction outcomes by relativizing fitness proxies to clonally infected and uninfected controls. More than one thousand root nodules of coinoculated plants were genotyped to quantify strain occupancy, and the Bradyrhizobium strain genome sequences were analyzed to uncover the genetic bases of interstrain competition outcomes. Strikingly, interstrain competition favored a fast-growing, minimally beneficial rhizobia strain. Host benefits were significantly diminished in coinoculation treatments relative to expectations from clonally inoculated controls, consistent with competitive interference among rhizobia that reduced both nodulation and plant growth. Competition traits appear polygenic, linked with inter-strain allelopathic interactions in the rhizosphere. This study confirms that competition among strains can destabilize mutualism by favoring microbes that are superior in colonizing host tissues but provide minimal benefits to host plants. Moreover, our findings help resolve the paradox that despite efficient host control post infection, legumes nonetheless encounter rhizobia that vary in their nitrogen fixation.

Whole Genome Sequencing-Based Tracing of a 2022 Introduction and Outbreak of Xanthomonas hortorum pv. pelargonii

Globalization has made agricultural commodities more accessible, available, and affordable. However, their global movement increases the potential for invasion by pathogens and necessitates development and implementation of sensitive, rapid, and scalable surveillance methods. Here, we used 35 strains, isolated by multiple diagnostic laboratories, as a case study for using whole genome sequence data in a plant disease diagnostic setting. Twenty-seven of the strains were isolated in 2022 and identified as Xanthomonas hortorum pv. pelargonii. Eighteen of these strains originated from material sold by a plant breeding company that had notified clients following a release of infected geranium cuttings. Analyses of whole genome sequences revealed epidemiological links among the 27 strains from different growers that confirmed a common source of the outbreak and uncovered likely secondary spread events within facilities that housed plants originating from different plant breeding companies. Whole genome sequencing data were also analyzed to reveal how preparatory and analytical methods can impact conclusions on outbreaks of clonal pathogenic strains. The results demonstrate the potential power of using whole genome sequencing among a network of diagnostic labs and highlight how sharing such data can help shorten response times to mitigate outbreaks more expediently and precisely than standard methods.

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.

Dynamic interactions between mega symbiosis ICEs and bacterial chromosomes maintain genome architecture

Acquisition of mobile genetic elements can confer novel traits to bacteria. Some integrative and conjugative elements confer upon members of Bradyrhizobium the capacity to fix nitrogen in symbiosis with legumes. These so-called symbiosis integrative conjugative elements (symICEs) can be extremely large and vary as monopartite and polypartite configurations within chromosomes of related strains. These features are predicted to impose fitness costs and have defied explanation. Here, we show that chromosome architecture is largely conserved despite diversity in genome composition, variations in locations of attachment sites recognized by integrases of symICEs, and differences in large-scale chromosomal changes that occur upon integration. Conversely, many simulated nonnative chromosome–symICE combinations are predicted to result in lethal deletions or disruptions to architecture. Findings suggest that there is compatibility between chromosomes and symICEs. We hypothesize that the size and structural flexibility of symICEs are important for generating combinations that maintain chromosome architecture across a genus of nitrogen-fixing bacteria with diverse and dynamic genomes.