Symbiosis genes show a unique pattern of introgression and selection within a Rhizobium leguminosarum species complex

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.

Rhizobium brockwellii sp. nov., Rhizobium johnstonii sp. nov. and Rhizobium beringeri sp. nov., three genospecies within the Rhizobium leguminosarum species complex

Genomic evidence indicates that the Rhizobium leguminosarum species complex comprises multiple distinct species, perhaps 18 or more. Of the five earliest genospecies (gs) to be described, only two have formal names: R. leguminosarum sensu stricto (gsE) and Rhizobium ruizarguesonis (gsC). Here, we provide formal descriptions and names for the other three genospecies, based on the publicly available genome sequences for multiple strains of each species: Rhizobium brockwellii sp. nov. (gsA, 37 strains, type strain CC275eT=LMG 6122T = ICMP 2163T=NZP 561T = PDDCC 2163T=HAMBI 13T), Rhizobium johnstonii sp. nov. (gsB, 54 strains, type strain 3841T = LMG 32736T=DSM 114642T) and Rhizobium beringeri sp. nov. (gsD, 8 strains, type strain SM51T = LMG 32895T = DSM 115206T). Each species forms a well-supported clade in a phylogeny based on 120 concatenated core genes. All strains have average nucleotide identity (ANI) above 96 % with the relevant type strain and below 96 % with all other type strains. Each species is characterised by a number of genes that are absent or rare in other species.

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.

Quantitative genetic analysis reveals potential to breed for improved white clover growth in symbiosis with nitrogen-fixing Rhizobium bacteria

White clover (Trifolium repens) is integral to mixed pastures in New Zealand and temperate agriculture globally. It provides quality feed and a sustainable source of plant-available nitrogen (N) via N-fixation through symbiosis with soil-dwelling Rhizobium bacteria. Improvement of N-fixation in white clover is a route to enhancing sustainability of temperate pasture production. Focussing on seedling growth critical for crop establishment and performance, a population of 120 half-sibling white clover families was assessed with either N-supplementation or N-fixation via inoculation with a commercial Rhizobium strain (TA1). Quantitative genetic analysis identified significant (p < 0.05) family additive genetic variance for Shoot and Root Dry Matter (DM) and Symbiotic Potential (SP), and Root to Shoot ratio. Estimated narrow-sense heritabilities for above-ground symbiotic traits were moderate (0.24-0.33), and the strong (r ≥ 0.97) genetic correlation between Shoot and Root DM indicated strong pleiotropy or close linkage. The moderate (r = 0.47) phenotypic correlation between Shoot DM under symbiosis vs. under N-supplementation suggested plant growth with mineral-N was not a strong predictor of symbiotic performance. At 5% among-family selection pressure, predicted genetic gains per selection cycle of 19 and 17% for symbiotic traits Shoot DM and Shoot SP, respectively, highlighted opportunities for improved early seedling establishment and growth under symbiosis. Single and multi-trait selection methods, including a Smith-Hazel index focussing on an ideotype of high Shoot DM and Shoot SP, showed commonality of top-ranked families among traits. This study provides a platform for proof-of-concept crosses to breed for enhanced seedling growth under Rhizobium symbiosis and is informative for other legume crops.

Genomes within genomes: nested symbiosis and its implications for plant evolution

Many important plant traits are products of nested symbiosis: mobile genetic elements (MGEs) are nested within microbes, which in turn, are nested within plants. Plant trait variation is therefore not only determined by the plant's genome, but also by loci within microbes and MGEs. Yet it remains unclear how interactions and coevolution within nested symbiosis impacts the evolution of plant traits. Despite the complexities of nested symbiosis, including nonadditive interactions, understanding the evolution of plant traits is facilitated by combining quantitative genetic and functional genomic approaches that explicitly consider sources of nested genetic variation (from loci in MGEs to microbiomes). Additionally, understanding coevolution within nested symbiosis enables us to design or select for MGEs that promote plant health.

Cryptic Prophages Contribution for Campylobacter jejuni and Campylobacter coli Introgression

Campylobacter coli and C. jejuni, the causing agents of campylobacteriosis, are described to be undergoing introgression events, i.e., the transference of genetic material between different species, with some isolates sharing almost a quarter of its genome. The participation of phages in introgression events and consequent impact on host ecology and evolution remain elusive. Three distinct prophages, named C. jejuni integrated elements 1, 2, and 4 (CJIE1, CJIE2, and CJIE4), are described in C. jejuni. Here, we identified two unreported prophages, Campylobacter coli integrated elements 1 and 2 (CCIE1 and CCIE2 prophages), which are C. coli homologues of CJIE1 and CJIE2, respectively. No induction was achieved for both prophages. Conversely, induction assays on CJIE1 and CJIE2 point towards the inducibility of these prophages. CCIE2-, CJIE1-, and CJIE4-like prophages were identified in a Campylobacter spp. population of 840 genomes, and phylogenetic analysis revealed clustering in three major groups: CJIE1-CCIE1, CJIE2-CCIE2, and CJIE4, clearly segregating prophages from C. jejuni and C. coli, but not from human- and nonhuman-derived isolates, corroborating the flowing between animals and humans in the agricultural context. Punctual bacteriophage host-jumps were observed in the context of C. jejuni and C. coli, and although random chance cannot be fully discarded, these observations seem to implicate prophages in evolutionary introgression events that are modulating the hybridization of C. jejuni and C. coli species.

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'.

Major effect loci for plant size before onset of nitrogen fixation allow accurate prediction of yield in white clover

Accurate genomic prediction of yield within and across generations was achieved by estimating the genetic merit of individual white clover genotypes based on extensive genetic replication using cloned material. White clover is an agriculturally important forage legume grown throughout temperate regions as a mixed clover-grass crop. It is typically cultivated with low nitrogen input, making yield dependent on nitrogen fixation by rhizobia in root nodules. Here, we investigate the effects of clover and rhizobium genetic variation by monitoring plant growth and quantifying dry matter yield of 704 combinations of 145 clover genotypes and 170 rhizobium inocula. We find no significant effect of rhizobium variation. In contrast, we can predict yield based on a few white clover markers strongly associated with plant size prior to nitrogen fixation, and the prediction accuracy for polycross offspring yield is remarkably high. Several of the markers are located near a homolog of Arabidopsis thaliana GIGANTUS 1, which regulates growth rate and biomass accumulation. Our work provides fundamental insight into the genetics of white clover yield and identifies specific candidate genes as breeding targets.

Introducing a Novel, Broad Host Range Temperate Phage Family Infecting Rhizobium leguminosarum and Beyond

Temperate phages play important roles in bacterial communities but have been largely overlooked, particularly in non-pathogenic bacteria. In rhizobia the presence of temperate phages has the potential to have significant ecological impacts but few examples have been described. Here we characterize a novel group of 5 Rhizobium leguminosarum prophages, capable of sustaining infections across a broad host range within their host genus. Genome comparisons identified further putative prophages infecting multiple Rhizobium species isolated globally, revealing a wider family of 10 temperate phages including one previously described lytic phage, RHEph01, which appears to have lost the ability to form lysogens. Phylogenetic discordance between prophage and host phylogenies suggests a history of active mobilization between Rhizobium lineages. Genome comparisons revealed conservation of gene content and order, with the notable exception of an approximately 5 kb region of hypervariability, containing almost exclusively hypothetical genes. Additionally, several horizontally acquired genes are present across the group, including a putative antirepressor present only in the RHEph01 genome, which may explain its apparent inability to form lysogens. In summary, both phenotypic and genomic comparisons between members of this group of phages reveals a clade of viruses with a long history of mobilization within and between Rhizobium species.

Genetic Variation in Host-Specific Competitiveness of the Symbiont Rhizobium leguminosarum Symbiovar viciae

Legumes of the Fabeae tribe form nitrogen-fixing root nodules resulting from symbiotic interaction with the soil bacteria Rhizobium leguminosarum symbiovar viciae (Rlv). These bacteria are all potential symbionts of the Fabeae hosts but display variable partner choice when co-inoculated in mixture. Because partner choice and symbiotic nitrogen fixation mostly behave as genetically independent traits, the efficiency of symbiosis is often suboptimal when Fabeae legumes are exposed to natural Rlv populations present in soil. A core collection of 32 Rlv bacteria was constituted based on the genomic comparison of a collection of 121 genome sequences, representative of known worldwide diversity of Rlv. A variable part of the nodD gene sequence was used as a DNA barcode to discriminate and quantify each of the 32 bacteria in mixture. This core collection was co-inoculated on a panel of nine genetically diverse Pisum sativum, Vicia faba, and Lens culinaris genotypes. We estimated the relative Early Partner Choice (EPC) of the bacteria with the Fabeae hosts by DNA metabarcoding on the nodulated root systems. Comparative genomic analyses within the bacterial core collection identified molecular markers associated with host-dependent symbiotic partner choice. The results revealed emergent properties of rhizobial populations. They pave the way to identify genes related to important symbiotic traits operating at this level.

Recombination facilitates adaptive evolution in rhizobial soil bacteria

Homologous recombination is expected to increase natural selection efficacy by decoupling the fate of beneficial and deleterious mutations and by readily creating new combinations of beneficial alleles. Here, we investigate how the proportion of amino acid substitutions fixed by adaptive evolution (α) depends on the recombination rate in bacteria. We analyze 3,086 core protein-coding sequences from 196 genomes belonging to five closely related species of the genus Rhizobium. These genes are found in all species and do not display any signs of introgression between species. We estimate α using the site frequency spectrum (SFS) and divergence data for all pairs of species. We evaluate the impact of recombination within each species by dividing genes into three equally sized recombination classes based on their average level of intragenic linkage disequilibrium. We find that α varies from 0.07 to 0.39 across species and is positively correlated with the level of recombination. This is both due to a higher estimated rate of adaptive evolution and a lower estimated rate of nonadaptive evolution, suggesting that recombination both increases the fixation probability of advantageous variants and decreases the probability of fixation of deleterious variants. Our results demonstrate that homologous recombination facilitates adaptive evolution measured by α in the core genome of prokaryote species in agreement with studies in eukaryotes.

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.

Genetic variation is associated with differences in facilitative and competitive interactions in the Rhizobium leguminosarum species complex

Competitive and facilitative interactions influence bacterial community composition, diversity and functioning. However, the role of genetic diversity for determining interactions between coexisting strains of the same, or closely related, species remains poorly understood. Here, we investigated the type (facilitative/inhibitory) and potential underlying mechanisms of pairwise interactions between 24 genetically diverse bacterial strains belonging to three genospecies (gsA,C,E) of the Rhizobium leguminosarum species complex. Interactions were determined indirectly, based on secreted compounds in cell-free supernatants, and directly, as growth inhibition in cocultures. We found supernatants mediated both facilitative and inhibitory interactions that varied greatly between strains and genospecies. Overall, gsE strains indirectly suppressed growth of gsA strains, while their own growth was facilitated by other genospecies' supernatants. Similar genospecies-level patterns were observed in direct competition, where gsA showed the highest susceptibility and gsE the highest inhibition capacity. At the genetic level, increased gsA susceptibility was associated with a non-random distribution of quorum sensing and secondary metabolite genes across genospecies. Together, our results suggest that genetic variation is associated with facilitative and competitive interactions, which could be important ecological mechanisms explaining R. leguminosarum diversity.

Impact of a bacterial consortium on the soil bacterial community structure and maize (Zea mays L.) cultivation

Microorganisms are often applied as biofertilizer to crops to stimulate plant growth, increase yields and reduce inorganic N application. The survival and proliferation of these allochthonous microorganisms in soil is a necessary requisite for them to promote plant growth. We applied a sterilized or unsterilized not commercialized bacterial consortium mixed with cow manure leachate used by a farmer as biofertilizer to maize (Zea mays L.) in a greenhouse experiment, while maize development and the bacterial community structure was determined just before the biofertilizer was applied a first time (day 44), after three applications (day 89) and after six application at the end of the experiment (day 130). Application of sterilized or unsterilized biofertilizer with pH 4.3 and 864 mg NH₄⁺-N kg^-1 had no significant effect on maize growth. The application of the biofertilizer dominated by Lactobacillus (relative abundance 11.90%) or the sterilized biofertilizer changed the relative abundance of a limited number of bacterial groups, i.e. Delftia, Halomonas, Lactobacillus and Stenotrophomonas, without altering significantly the bacterial community structure. Cultivation of maize, however, affected significantly the bacterial community structure, which showed large significant variations over time in the cultivated and uncultivated soil. It was concluded that the bacteria applied as a biofertilizer had only a limited effect on the relative abundance of these groups in uncultivated or soil cultivated with maize.

The impact of intra-specific diversity in the rhizobia-legume symbiosis

Rhizobia - nitrogen-fixing, root-nodulating bacteria - play a critical role in both plant ecosystems and sustainable agriculture. Rhizobia form intracellular infections within legumes roots where they produce plant accessible nitrogen from atmospheric nitrogen and thus reduce the reliance on industrial inputs. The rhizobia-legume symbiosis is often treated as a pairwise relationship between single genotypes, both in research and in the production of rhizobial inoculants. However in nature individual plants are infected by a high diversity of rhizobia symbionts. How this diversity affects productivity within the symbiosis is unclear. Here, we use a powerful statistical approach to assess the impact of diversity within the Rhizobium leguminosarum - clover symbiosis using a biodiversity-ecosystem function framework. Statistically, we found no significant impact of rhizobium diversity. However this relationship was weakly positive - rather than negative - indicating that there is no significant cost to increasing inoculant diversity. Productivity was influenced by the identity of the strains within an inoculant; strains with the highest individual performance showed a significant positive contribution within mixed inoculants. Overall, inoculant effectiveness was best predicted by the individual performance of the best inoculant member, and only weakly predicted by the worst performing member. Collectively, our data suggest that the Rhizobium leguminosarum - clover symbiosis displays a weak diversity-function relationship, but that inoculant performance can be improved through the inclusion of high performing strains. Given the wide environmental dependence of rhizobial inoculant quality, multi-strain inoculants could be highly successful as they increase the likelihood of including a strain well adapted to local conditions across different environments.

Defining the Rhizobium leguminosarum Species Complex

Bacteria currently included in Rhizobium leguminosarum are too diverse to be considered a single species, so we can refer to this as a species complex (the Rlc). We have found 429 publicly available genome sequences that fall within the Rlc and these show that the Rlc is a distinct entity, well separated from other species in the genus. Its sister taxon is R. anhuiense. We constructed a phylogeny based on concatenated sequences of 120 universal (core) genes, and calculated pairwise average nucleotide identity (ANI) between all genomes. From these analyses, we concluded that the Rlc includes 18 distinct genospecies, plus 7 unique strains that are not placed in these genospecies. Each genospecies is separated by a distinct gap in ANI values, usually at approximately 96% ANI, implying that it is a 'natural' unit. Five of the genospecies include the type strains of named species: R. laguerreae, R. sophorae, R. ruizarguesonis, "R. indicum" and R. leguminosarum itself. The 16S ribosomal RNA sequence is remarkably diverse within the Rlc, but does not distinguish the genospecies. Partial sequences of housekeeping genes, which have frequently been used to characterize isolate collections, can mostly be assigned unambiguously to a genospecies, but alleles within a genospecies do not always form a clade, so single genes are not a reliable guide to the true phylogeny of the strains. We conclude that access to a large number of genome sequences is a powerful tool for characterizing the diversity of bacteria, and that taxonomic conclusions should be based on all available genome sequences, not just those of type strains.

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.

MAUI-seq: Metabarcoding using amplicons with unique molecular identifiers to improve error correction

Sequencing and PCR errors are a major challenge when characterizing genetic diversity using high-throughput amplicon sequencing (HTAS). We have developed a multiplexed HTAS method, MAUI-seq, which uses unique molecular identifiers (UMIs) to improve error correction by exploiting variation among sequences associated with a single UMI. Erroneous sequences are recognized because, across the data set, they are over-represented among the minor sequences associated with UMIs. We show that two main advantages of this approach are efficient elimination of chimeric and other erroneous reads, outperforming dada2 and unoise3, and the ability to confidently recognize genuine alleles that are present at low abundance or resemble chimeras. The method provides sensitive and flexible profiling of diversity and is readily adaptable to most HTAS applications, including microbial 16S rRNA profiling and metabarcoding of environmental DNA.

Diverse Rhizobium strains isolated from root nodules of Trifolium alexandrinum in Egypt and symbiovars

Berseem clover (T. alexandrinum) is the main forage legume crop used as animal feed in Egypt. Here, eighty rhizobial isolates were isolated from root nodules of berseem clover grown in different regions in Egypt and were grouped by RFLP-16S rRNA ribotyping. Representative isolates were characterized using phylogenetic analyses of the 16S rRNA, rpoB, glnA, pgi, and nodC genes. We also investigated the performance of these isolates using phenotypic tests and nitrogen fixation efficiency assays. The majority of strains (<90%) were closely related to Rhizobium aegyptiacum and Rhizobium aethiopicum and of the remaining strains, six belonged to the Rhizobium leguminosarum genospecies complex and only one strain was assigned to Agrobacterium fabacearum. Despite their heterogeneous chromosomal background, most of the strains shared nodC gene alleles corresponding to symbiovar trifolii. Some of the strains closely affiliated to R. aegyptiacum and R. aethiopicum had superior nodulation and nitrogen fixation capabilities in berseem clover, compared to the commercial inoculant (Okadein®) and N-added treatments. R. leguminosarum strain NGB-CR 17 that harbored a nodC allele typical of symbiovar viciae, was also able to form an effective symbiosis with clover. Two strains with nodC alleles of symbiovar trifolii, R. aegyptiacum strains NGB-CR 129 and 136, were capable of forming effective nodules in Phaseolus vulgaris in axenic greenhouse conditions. This adds the symbiovar trifolii which is well-established in the Egyptian soils to the list of symbiovars that form nodules in P. vulgaris.

Greenotyper: Image-Based Plant Phenotyping Using Distributed Computing and Deep Learning

Image-based phenotype data with high temporal resolution offers advantages over end-point measurements in plant quantitative genetics experiments, because growth dynamics can be assessed and analysed for genotype-phenotype association. Recently, network-based camera systems have been deployed as customizable, low-cost phenotyping solutions. Here, we implemented a large, automated image-capture system based on distributed computing using 180 networked Raspberry Pi units that could simultaneously monitor 1,800 white clover (Trifolium repens) plants. The camera system proved stable with an average uptime of 96% across all 180 cameras. For analysis of the captured images, we developed the Greenotyper image analysis pipeline. It detected the location of the plants with a bounding box accuracy of 97.98%, and the U-net-based plant segmentation had an intersection over union accuracy of 0.84 and a pixel accuracy of 0.95. We used Greenotyper to analyze a total of 355,027 images, which required 24-36 h. Automated phenotyping using a large number of static cameras and plants thus proved a cost-effective alternative to systems relying on conveyor belts or mobile cameras.