RNA sequencing and analysis of three Lupinus nodulomes provide new insights into specific host-symbiont relationships with compatible and incompatible Bradyrhizobium strains

Nitrogen fixation in the legume root-nodule symbiosis has a critical importance in natural and agricultural ecosystems and depends on the proper choice of the symbiotic partners. However, the genetic determinism of symbiotic specificity remains unclear. To study this process, we inoculated three Lupinus species (L. albus, L. luteus, L. mariae-josephae), belonging to the under-investigated tribe of Genistoids, with two Bradyrhizobium strains (B. japonicum, B. valentinum) presenting contrasted degrees of symbiotic specificity depending on the host. We produced the first transcriptomes (RNA-Seq) from lupine nodules in a context of symbiotic specificity. For each lupine species, we compared gene expression between functional and non-functional interactions and determined differentially expressed (DE) genes. This revealed that L. luteus and L. mariae-josephae (nodulated by only one of the Bradyrhizobium strains) specific nodulomes were richest in DE genes than L. albus (nodulation with both microsymbionts, but non-functional with B. valentinum) and share a higher number of these genes between them than with L. albus. In addition, a functional analysis of DE genes highlighted the central role of the genetic pathways controlling infection and nodule organogenesis, hormones, secondary, carbon and nitrogen metabolisms, as well as the implication of plant defence in response to compatible or incompatible Bradyrhizobium strains.