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

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.

Herbivory modifies plant symbiont number and impact on host plant performance in the field

Species interactions are a unifying theme in ecology and evolution. Both fields are currently moving beyond their historical focus on isolated pairwise relationships to understand how ecological communities affect focal interactions. Additional species can modify both the number of interactions and the fitness consequences of each interaction (i.e., selection). Although only selection affects the evolution of the focal interaction, the two are often conflated, limiting our understanding of the evolution of multispecies interactions. We manipulated aboveground herbivory on the legume Medicago lupulina in the field and quantified its effect on number of symbionts and the per-symbiont impact on plant performance in two belowground symbioses: mutualistic rhizobia bacteria (Ensifer meliloti) and parasitic root-knot nematodes (Meloidogyne hapla). We found that herbivores modified the number of rhizobia nodules, as well as the benefit per nodule. However, each effect was specific to a distinct herbivory regime: natural herbivory affected nodule number, whereas leafhoppers (Cicadellidae) weakened the per nodule benefit. We did not detect any effect of herbivory on nematode gall number or the cost of infection. Our data demonstrate that distinguishing between symbiont number from the fitness consequences of symbiosis is crucial to accurately infer how pairwise interactions will evolve in a community.

Host-Associated Rhizobial Fitness: Dependence on Nitrogen, Density, Community Complexity, and Legume Genotype

The environmental context of the nitrogen-fixing mutualism between leguminous plants and rhizobial bacteria varies over space and time. Variation in resource availability, population density, and composition likely affect the ecology and evolution of rhizobia and their symbiotic interactions with hosts. We examined how host genotype, nitrogen addition, rhizobial density, and community complexity affected selection on 68 rhizobial strains in the Sinorhizobium meliloti–Medicago truncatula mutualism. As expected, host genotype had a substantial effect on the size, number, and strain composition of root nodules (the symbiotic organ). The understudied environmental variable of rhizobial density had a stronger effect on nodule strain frequency than the addition of low nitrogen levels. Higher inoculum density resulted in a nodule community that was less diverse and more beneficial but only in the context of the more selective host genotype. Higher density resulted in more diverse and less beneficial nodule communities with the less selective host. Density effects on strain composition deserve additional scrutiny as they can create feedback between ecological and evolutionary processes. Finally, we found that relative strain rankings were stable across increasing community complexity (2, 3, 8, or 68 strains). This unexpected result suggests that higher-order interactions between strains are rare in the context of nodule formation and development. Our work highlights the importance of examining mechanisms of density-dependent strain fitness and developing theoretical predictions that incorporate density dependence. Furthermore, our results have translational relevance for overcoming establishment barriers in bioinoculants and motivating breeding programs that maintain beneficial plant-microbe interactions across diverse agroecological contexts.

IMPORTANCE Legume crops establish beneficial associations with rhizobial bacteria that perform biological nitrogen fixation, providing nitrogen to plants without the economic and greenhouse gas emission costs of chemical nitrogen inputs. Here, we examine the influence of three environmental factors that vary in agricultural fields on strain relative fitness in nodules. In addition to manipulating nitrogen, we also use two biotic variables that have rarely been examined: the rhizobial community's density and complexity. Taken together, our results suggest that (i) breeding legume varieties that select beneficial strains despite environmental variation is possible, (ii) changes in rhizobial population densities that occur routinely in agricultural fields could drive evolutionary changes in rhizobial populations, and (iii) the lack of higher-order interactions between strains will allow the high-throughput assessments of rhizobia winners and losers during plant interactions.

Phenotypic and genomic signatures of interspecies cooperation and conflict in naturally occurring isolates of a model plant symbiont

Given the need to predict the outcomes of (co)evolution in host-associated microbiomes, whether microbial and host fitnesses tend to trade-off, generating conflict, remains a pressing question. Examining the relationships between host and microbe fitness proxies at both the phenotypic and genomic levels can illuminate the mechanisms underlying interspecies cooperation and conflict. We examined naturally occurring genetic variation in 191 strains of the model microbial symbiont Sinorhizobium meliloti, paired with each of two host Medicago truncatula genotypes in single- or multi-strain experiments to determine how multiple proxies of microbial and host fitness were related to one another and test key predictions about mutualism evolution at the genomic scale, while also addressing the challenge of measuring microbial fitness. We found little evidence for interspecies fitness conflict; loci tended to have concordant effects on both microbe and host fitnesses, even in environments with multiple co-occurring strains. Our results emphasize the importance of quantifying microbial relative fitness for understanding microbiome evolution and thus harnessing microbiomes to improve host fitness. Additionally, we find that mutualistic coevolution between hosts and microbes acts to maintain, rather than erode, genetic diversity, potentially explaining why variation in mutualism traits persists in nature.

Combining GWAS and population genomic analyses to characterize coevolution in a legume-rhizobia symbiosis

The mutualism between legumes and rhizobia is clearly the product of past coevolution. However, the nature of ongoing evolution between these partners is less clear. To characterize the nature of recent coevolution between legumes and rhizobia, we used population genomic analysis to characterize selection on functionally annotated symbiosis genes as well as on symbiosis gene candidates identified through a two-species association analysis. For the association analysis, we inoculated each of 202 accessions of the legume host Medicago truncatula with a community of 88 Sinorhizobia (Ensifer) meliloti strains. Multistrain inoculation, which better reflects the ecological reality of rhizobial selection in nature than single-strain inoculation, allows strains to compete for nodulation opportunities and host resources and for hosts to preferentially form nodules and provide resources to some strains. We found extensive host by symbiont, that is, genotype-by-genotype, effects on rhizobial fitness and some annotated rhizobial genes bear signatures of recent positive selection. However, neither genes responsible for this variation nor annotated host symbiosis genes are enriched for signatures of either positive or balancing selection. This result suggests that stabilizing selection dominates selection acting on symbiotic traits and that variation in these traits is under mutation-selection balance. Consistent with the lack of positive selection acting on host genes, we found that among-host variation in growth was similar whether plants were grown with rhizobia or N-fertilizer, suggesting that the symbiosis may not be a major driver of variation in plant growth in multistrain contexts.

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.

Prioritizing host phenotype to understand microbiome heritability in plants

Breeders and evolutionary geneticists have grappled with the complexity of the 'genotype-to-phenotype map' for decades. Now, recent studies highlight the relevance of this concept for understanding heritability of plant microbiomes. Because host phenotype is a more proximate cause of microbiome variation than host genotype, microbiome heritability varies across plant anatomy and development. Fine-scale variation of plant traits within organs suggests that the well-established concept of 'microbiome compartment' should be refined. Additionally, recent work shows that the balance of deterministic processes (including host genetic effects) vs stochastic processes also varies over time and space. Together, these findings suggest that re-centering plant phenotype - both as a predictor and a readout of microbiome function - will accelerate new insights into microbiome heritability.

Toward Unifying Evolutionary Ecology and Genomics to Understand Positive Plant-Plant Interactions Within Wild Species

In a local environment, plant networks include interactions among individuals of different species and among genotypes of the same species. While interspecific interactions are recognized as main drivers of plant community patterns, intraspecific interactions have recently gained attention in explaining plant community dynamics. However, an overview of intraspecific genotype-by-genotype interaction patterns within wild plant species is still missing. From the literature, we identified 91 experiments that were mainly designed to investigate the presence of positive interactions based on two contrasting hypotheses. Kin selection theory predicts partisan help given to a genealogical relative. The rationale behind this hypothesis relies on kin/non-kin recognition, with the positive outcome of kin cooperation substantiating it. On the other hand, the elbow-room hypothesis supports intraspecific niche partitioning leading to positive outcome when genetically distant genotypes interact. Positive diversity-productivity relationship rationalizes this hypothesis, notably with the outcome of overyielding. We found that both these hypotheses have been highly supported in experimental studies despite their opposite predictions between the extent of genetic relatedness among neighbors and the level of positive interactions. Interestingly, we identified a highly significant effect of breeding system, with a high proportion of selfing species associated with the presence of kin cooperation. Nonetheless, we identified several shortcomings regardless of the species considered, such as the lack of a reliable estimate of genetic relatedness among genotypes and ecological characterization of the natural habitats from which genotypes were collected, thereby impeding the identification of selective drivers of positive interactions. We therefore propose a framework combining evolutionary ecology and genomics to establish the eco-genomic landscape of positive GxG interactions in wild plant species.

Whose trait is it anyways? Coevolution of joint phenotypes and genetic architecture in mutualisms

Evolutionary biologists typically envision a trait's genetic basis and fitness effects occurring within a single species. However, traits can be determined by and have fitness consequences for interacting species, thus evolving in multiple genomes. This is especially likely in mutualisms, where species exchange fitness benefits and can associate over long periods of time. Partners may experience evolutionary conflict over the value of a multi-genomic trait, but such conflicts may be ameliorated by mutualism's positive fitness feedbacks. Here, we develop a simulation model of a host-microbe mutualism to explore the evolution of a multi-genomic trait. Coevolutionary outcomes depend on whether hosts and microbes have similar or different optimal trait values, strengths of selection and fitness feedbacks. We show that genome-wide association studies can map joint traits to loci in multiple genomes and describe how fitness conflict and fitness feedback generate different multi-genomic architectures with distinct signals around segregating loci. Partner fitnesses can be positively correlated even when partners are in conflict over the value of a multi-genomic trait, and conflict can generate strong mutualistic dependency. While fitness alignment facilitates rapid adaptation to a new optimum, conflict maintains genetic variation and evolvability, with implications for applied microbiome science.

Evolving together, evolving apart: measuring the fitness of rhizobial bacteria in and out of symbiosis with leguminous plants

Most plant-microbe interactions are facultative, with microbes experiencing temporally and spatially variable selection. How this variation affects microbial evolution is poorly understood. Given its tractability and ecological and agricultural importance, the legume-rhizobia nitrogen-fixing symbiosis is a powerful model for identifying traits and genes underlying bacterial fitness. New technologies allow high-throughput measurement of the relative fitness of bacterial mutants, strains and species in mixed inocula in the host, rhizosphere and soil environments. I consider how host genetic variation (G × G), other environmental factors (G × E), and host life-cycle variation may contribute to the maintenance of genetic variation and adaptive trajectories of rhizobia - and, potentially, other facultative symbionts. Lastly, I place these findings in the context of developing beneficial inoculants in a changing climate.

Diversity and plant growth-promoting functions of diazotrophic/N-scavenging bacteria isolated from the soils and rhizospheres of two species of Solanum

Studies of the interactions between plants and their microbiome have been conducted worldwide in the search for growth-promoting representative strains for use as biological inputs for agriculture, aiming to achieve more sustainable agriculture practices. With a focus on the isolation of plant growth-promoting (PGP) bacteria with ability to alleviate N stress, representative strains that were found at population densities greater than 10⁴ cells g^-1 and that could grow in N-free semisolid media were isolated from soils under different management conditions and from the roots of tomato (Solanum lycopersicum) and lulo (Solanum quitoense) plants that were grown in those soils. A total of 101 bacterial strains were obtained, after which they were phylogenetically categorized and characterized for their basic PGP mechanisms. All strains belonged to the Proteobacteria phylum in the classes Alphaproteobacteria (61% of isolates), Betaproteobacteria (19% of isolates) and Gammaproteobacteria (20% of isolates), with distribution encompassing nine genera, with the predominant genus being Rhizobium (58.4% of isolates). Strains isolated from conventional horticulture (CH) soil composed three bacterial genera, suggesting a lower diversity for the diazotrophs/N scavenger bacterial community than that observed for soils under organic management (ORG) or secondary forest coverture (SF). Conversely, diazotrophs/N scavenger strains from tomato plants grown in CH soil comprised a higher number of bacterial genera than did strains isolated from tomato plants grown in ORG or SF soils. Furthermore, strains isolated from tomato were phylogenetically more diverse than those from lulo. BOX-PCR fingerprinting of all strains revealed a high genetic diversity for several clonal representatives (four Rhizobium species and one Pseudomonas species). Considering the potential PGP mechanisms, 49 strains (48.5% of the total) produced IAA (2.96–193.97 μg IAA mg protein^-1), 72 strains (71.3%) solubilized FePO₄ (0.40–56.00 mg l^-1), 44 strains (43.5%) solubilized AlPO₄ (0.62–17.05 mg l^-1), and 44 strains produced siderophores (1.06–3.23). Further, 91 isolates (90.1% of total) showed at least one PGP trait, and 68 isolates (67.3%) showed multiple PGP traits. Greenhouse trials using the bacterial collection to inoculate tomato or lulo plants revealed increases in plant biomass (roots, shoots or both plant tissues) elicited by 65 strains (54.5% of the bacterial collection), of which 36 were obtained from the tomato rhizosphere, 15 were obtained from the lulo rhizosphere, and 14 originated from samples of soil that lacked plants. In addition, 18 strains showed positive inoculation effects on both Solanum species, of which 12 were classified as Rhizobium spp. by partial 16S rRNA gene sequencing. Overall, the strategy adopted allowed us to identify the variability in the composition of culturable diazotroph/N-scavenger representatives from soils under different management conditions by using two Solanum species as trap plants. The present results suggest the ability of tomato and lulo plants to enrich their belowground microbiomes with rhizobia representatives and the potential of selected rhizobial strains to promote the growth of Solanum crops under limiting N supply.

Rhizobia Isolated from the Relict Legume Vavilovia formosa Represent a Genetically Specific Group within Rhizobium leguminosarum biovar viciae

Twenty-two rhizobia strains isolated from three distinct populations (North Ossetia, Dagestan, and Armenia) of a relict legume Vavilovia formosa were analysed to determine their position within Rhizobium leguminosarum biovar viciae (Rlv). These bacteria are described as symbionts of four plant genera Pisum, Vicia, Lathyrus, and Lens from the Fabeae tribe, of which Vavilovia is considered to be closest to its last common ancestor (LCA). In contrast to biovar viciae, bacteria from Rhizobium leguminosarum biovar trifolii (Rlt) inoculate plants from the Trifolieae tribe. Comparison of house-keeping (hkg: 16S rRNA, glnII, gltA, and dnaK) and symbiotic (sym: nodA, nodC, nodD, and nifH) genes of the symbionts of V. formosa with those of other Rlv and Rlt strains reveals a significant group separation, which was most pronounced for sym genes. A remarkable feature of the strains isolated from V. formosa was the presence of the nodX gene, which was commonly found in Rlv strains isolated from Afghanistan pea genotypes. Tube testing of different strains on nine plant species, including all genera from the Fabeae tribe, demonstrated that the strains from V. formosa nodulated the same cross inoculation group as the other Rlv strains. Comparison of nucleotide similarity in sym genes suggested that their diversification within sym-biotypes of Rlv was elicited by host plants. Contrariwise, that of hkg genes could be caused by either local adaptation to soil niches or by genetic drift. Long-term ecological isolation, genetic separation, and the ancestral position of V. formosa suggested that symbionts of V. formosa could be responsible for preserving ancestral genotypes of the Rlv biovar.