Mobile gene clusters and coexpressed plant-rhizobium pathways drive partner quality variation in symbiosis

Plant-microbe symbioses such as the legume-rhizobium mutualism are vital in the web of ecological relationships within both natural and managed ecosystems, influencing primary productivity, crop yield, and ecosystem services. The outcome of these interactions for plant hosts varies quantitatively and can range from highly beneficial to even detrimental depending on natural genetic variation in microbial symbionts. Here, we take a systems genetics approach, harnessing the genetic diversity present in wild rhizobial populations to predict genes and molecular pathways crucial in determining partner quality, i.e., the benefits of symbiosis for legume hosts. We combine traits, dual-RNAseq of both partners from active nodules, pangenomics/pantranscriptomics, and Weighted Gene Co-expression Network Analysis (WGCNA) for a panel of 20 Sinorhizobium meliloti strains that vary in symbiotic partner quality. We find that genetic variation in the nodule transcriptome predicts host plant biomass, and WGCNA reveals networks of genes in plants and rhizobia that are coexpressed and associated with high-quality symbiosis. Presence-absence variation of gene clusters on the symbiosis plasmid (pSymA), validated in planta, is associated with high or low-quality symbiosis and is found within important coexpression modules. Functionally our results point to management of oxidative stress, amino acid and carbohydrate transport, and NCR peptide signaling mechanisms in driving symbiotic outcomes. Our integrative approach highlights the complex genetic architecture of microbial partner quality and raises hypotheses about the genetic mechanisms and evolutionary dynamics of symbiosis.