Coevolutionary genetic variation in the legume-rhizobium transcriptome

Coevolutionary change requires reciprocal selection between interacting species, where the partner genotypes that are favoured in one species depend on the genetic composition of the interacting species. Coevolutionary genetic variation is manifested as genotype x genotype (G x G) interactions for fitness in interspecific interactions. Although quantitative genetic approaches have revealed abundant evidence for G x G interactions in symbioses, the molecular basis of this variation remains unclear. Here we study the molecular basis of G x G interactions in a model legume-rhizobium mutualism using gene expression microarrays. We find that, like quantitative traits such as fitness, variation in the symbiotic transcriptome may be partitioned into additive and interactive genetic components. Our results suggest that plant genetic variation had the largest influence on nodule gene expression and that plant genotype and the plant genotype x rhizobium genotype interaction determine global shifts in rhizobium gene expression that in turn feedback to influence plant fitness benefits. Moreover, the transcriptomic variation we uncover implicates regulatory changes in both species as drivers of symbiotic gene expression variation. Our study is the first to partition genetic variation in a symbiotic transcriptome and illuminates potential molecular routes of coevolutionary change.