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.

A Select and Resequence Approach Reveals Strain-Specific Effects of Medicago Nodule-Specific PLAT-Domain Genes.

Genetic studies of legume symbiosis with nitrogen-fixing rhizobial bacteria have traditionally focused on nodule and nitrogen-fixation phenotypes when hosts are inoculated with a single rhizobial strain. These approaches overlook the potential effect of host genes on rhizobial fitness (i.e. how many rhizobia are released from host nodules) and strain-specific effects of host genes (i.e. genome × genome interactions). Using Medicago truncatula mutants in the recently described nodule-specific PLAT domain (NPD) gene family, we show how inoculating plants with a mixed inoculum of 68 rhizobial strains (Ensifer meliloti) via a select-and-resequence approach can be used to efficiently assay host mutants for strain-specific effects of late-acting host genes on interacting bacteria. The deletion of a single NPD gene (npd2) or all five members of the NPD gene family (npd1-5) differentially altered the frequency of rhizobial strains in nodules even though npd2 mutants had no visible nodule morphology or N-fixation phenotype. Also, npd1-5 nodules were less diverse and had larger populations of colony-forming rhizobia despite their smaller size. Lastly, NPD mutations disrupt a positive correlation between strain fitness and wild-type host biomass. These changes indicate that the effects of NPD proteins are strain dependent and that NPD family members are not redundant with regard to their effects on rhizobial strains. Association analyses of the rhizobial strains in the mixed inoculation indicate that rhizobial genes involved in chromosome segregation, cell division, GABA metabolism, efflux systems, and stress tolerance play an important role in the strain-specific effects of NPD genes.

Nodule-specific PLAT domain proteins are expanded in the Medicago lineage and required for nodulation.

Symbiotic nitrogen fixation in legumes is mediated by an interplay of signaling processes between plant hosts and rhizobial symbionts. In legumes, several secreted protein families have undergone expansions and play key roles in nodulation. Thus, identifying lineage-specific expansions (LSEs) of nodulation-associated genes can be a strategy to discover candidate gene families. Using bioinformatic tools, we identified 13 LSEs of nodulation-related secreted protein families, each unique to either Glycine, Arachis or Medicago lineages. In the Medicago lineage, nodule-specific Polycystin-1, Lipoxygenase, Alpha Toxin (PLAT) domain proteins (NPDs) expanded to five members. We examined NPD function using CRISPR/Cas9 multiplex genome editing to create Medicago truncatula NPD knockout lines, targeting one to five NPD genes. Mutant lines with differing combinations of NPD gene inactivations had progressively smaller nodules, earlier onset of nodule senescence, or ineffective nodules compared to the wild-type control. Double- and triple-knockout lines showed dissimilar nodulation phenotypes but coincided in upregulation of a DHHC-type zinc finger and an aspartyl protease gene, possible candidates for the observed disturbance of proper nodule function. By postulating that gene family expansions can be used to detect candidate genes, we identified a family of nodule-specific PLAT domain proteins and confirmed that they play a role in successful nodule formation.

Validating Genome-Wide Association Candidates Controlling Quantitative Variation in Nodulation.

Genome-wide association (GWA) studies offer the opportunity to identify genes that contribute to naturally occurring variation in quantitative traits. However, GWA relies exclusively on statistical association, so functional validation is necessary to make strong claims about gene function. We used a combination of gene-disruption platforms (Tnt1 retrotransposons, hairpin RNA-interference constructs, and CRISPR/Cas9 nucleases) together with randomized, well-replicated experiments to evaluate the function of genes that an earlier GWA study in Medicago truncatula had identified as candidates contributing to variation in the symbiosis between legumes and rhizobia. We evaluated ten candidate genes found in six clusters of strongly associated single nucleotide polymorphisms, selected on the basis of their strength of statistical association, proximity to annotated gene models, and root or nodule expression. We found statistically significant effects on nodule production for three candidate genes, each validated in two independent mutants. Annotated functions of these three genes suggest their contributions to quantitative variation in nodule production occur through processes not previously connected to nodulation, including phosphorous supply and salicylic acid-related defense response. These results demonstrate the utility of GWA combined with reverse mutagenesis technologies to discover and validate genes contributing to naturally occurring variation in quantitative traits. The results highlight the potential for GWA to complement forward genetics in identifying the genetic basis of ecologically and economically important traits.

Genomic characterization of the LEED..PEEDs, a gene family unique to the medicago lineage.

The LEED..PEED (LP) gene family in Medicago truncatula (A17) is composed of 13 genes coding small putatively secreted peptides with one to two conserved domains of negatively charged residues. This family is not present in the genomes of Glycine max, Lotus japonicus, or the IRLC species Cicer arietinum. LP genes were also not detected in a Trifolium pratense draft genome or Pisum sativum nodule transcriptome, which were sequenced de novo in this study, suggesting that the LP gene family arose within the past 25 million years. M. truncatula accession HM056 has 13 LP genes with high similarity to those in A17, whereas M. truncatula ssp. tricycla (R108) and M. sativa have 11 and 10 LP gene copies, respectively. In M. truncatula A17, 12 LP genes are located on chromosome 7 within a 93-kb window, whereas one LP gene copy is located on chromosome 4. A phylogenetic analysis of the gene family is consistent with most gene duplications occurring prior to Medicago speciation events, mainly through local tandem duplications and one distant duplication across chromosomes. Synteny comparisons between R108 and A17 confirm that gene order is conserved between the two subspecies, although a further duplication occurred solely in A17. In M. truncatula A17, all 13 LPs are exclusively transcribed in nodules and absent from other plant tissues, including roots, leaves, flowers, seeds, seed shells, and pods. The recent expansion of LP genes in Medicago spp. and their timing and location of expression suggest a novel function in nodulation, possibly as an aftermath of the evolution of bacteroid terminal differentiation or potentially associated with rhizobial-host specificity.