Search for: All records

Abstract Mutualism breakdown can be prevented if partner species preferentially select and reward partners that provide greater benefit. We examined these two components using the legume Medicago truncatula and its nitrogen-fixing symbiont Sinorhizobium meliloti. First, we re-analysed data from 202 accessions to show significant genetic variation in the capacity of Medicago to restrict strain diversity, finding that hosts with shorter nodules were more selective. A genome-wide association study on host selectivity identified genes including the hormone leginsulin, pectin degradation, multidrug and toxic compound efflux, zinc transport and DNA methylation. Second, we used two well-studied Medicago genotypes with contrasting nodule morphologies to assess the effectiveness of adaptive enrichment mechanisms by sampling the relative frequencies of rhizobial strains in pools of small nodules (indicating a lack of host investment) compared to large nodules (indicating increased host investment) and pairing these results with previous single-strain assessments of strain benefits to hosts. While both hosts enriched beneficial strains in large nodules, the host that formed larger and more variably sized nodules and thus had greater 'potential' to increase rhizobial populations was less effective. Our findings reveal that host genetic variation affects strain selectivity and suggest that nodule morphology traits warrant attention when exploring mutualism evolution. 

Abstract In the mutualism between leguminous plants and rhizobial bacteria, rhizobia live inside root nodules, creating potential for host genes to shape the rhizobial selective environment. Many host genes that affect symbiosis have been identified; however, the extent to which these genes affect selection acting on rhizobia is unknown. In this study, we inoculated 18 Medicago truncatula symbiotic mutants (including mutants that alter Nodule Cysteine-Rich (NCR) peptide production, plant defence, and nodule number regulation) with a mixture of 86 Sinorhizobium meliloti strains. Most mutations resulted in reduced host benefits, but the effects on rhizobial benefit (i.e. relative strain fitness) varied widely, revealing widespread host-by-strain fitness interactions. Genome-wide association analyses identified variants on rhizobial replicons pSymA and pSymB as important in mediating strain fitness responses to host mutations. Whereas most top variants affected rhizobial fitness with one host mutation (limited effect variants), nine affected fitness across six or more host mutations. These pervasive variants occurred primarily on pSymA, the symbiotic replicon, and include fixL and some metabolic genes. In contrast to the limited effect variants, variants with pervasive positive effects on strain fitness when host genes were mutated tended to adversely affect fitness in wild-type hosts. Competition assays across Medicago genotypes confirmed a pervasive role for one candidate (malonyl-CoA synthase), and AlphaFold multimer modelling suggests that many rhizobial top candidates could interact with host NCR peptides. Our results reveal how host genetic mutations alter strain fitness, setting the stage for improving rhizobial inoculants and breeding legume hosts better adapted to multi-strain environments. 

ABSTRACT Genome-wide association studies (GWAS) can identify genetic variants responsible for naturally occurring and quantitative phenotypic variation. Association studies therefore provide a powerful complement to approaches that rely on de novo mutations for characterizing gene function. Although bacteria should be amenable to GWAS, few GWAS have been conducted on bacteria, and the extent to which nonindependence among genomic variants (e.g., linkage disequilibrium [LD]) and the genetic architecture of phenotypic traits will affect GWAS performance is unclear. We apply association analyses to identify candidate genes underlying variation in 20 biochemical, growth, and symbiotic phenotypes among 153 strains of Ensifer meliloti . For 11 traits, we find genotype-phenotype associations that are stronger than expected by chance, with the candidates in relatively small linkage groups, indicating that LD does not preclude resolving association candidates to relatively small genomic regions. The significant candidates show an enrichment for nucleotide polymorphisms (SNPs) over gene presence-absence variation (PAV), and for five traits, candidates are enriched in large linkage groups, a possible signature of epistasis. Many of the variants most strongly associated with symbiosis phenotypes were in genes previously identified as being involved in nitrogen fixation or nodulation. For other traits, apparently strong associations were not stronger than the range of associations detected in permuted data. In sum, our data show that GWAS in bacteria may be a powerful tool for characterizing genetic architecture and identifying genes responsible for phenotypic variation. However, careful evaluation of candidates is necessary to avoid false signals of association. IMPORTANCE Genome-wide association analyses are a powerful approach for identifying gene function. These analyses are becoming commonplace in studies of humans, domesticated animals, and crop plants but have rarely been conducted in bacteria. We applied association analyses to 20 traits measured in Ensifer meliloti , an agriculturally and ecologically important bacterium because it fixes nitrogen when in symbiosis with leguminous plants. We identified candidate alleles and gene presence-absence variants underlying variation in symbiosis traits, antibiotic resistance, and use of various carbon sources; some of these candidates are in genes previously known to affect these traits whereas others were in genes that have not been well characterized. Our results point to the potential power of association analyses in bacteria, but also to the need to carefully evaluate the potential for false associations.