Dual RNA-seq reveals large-scale non-conserved genotype × genotype-specific genetic reprograming and molecular crosstalk in the mycorrhizal symbiosis

Arbuscular mycorrhizal fungi (AMF) impact plant growth and are a major driver of plant diversity and productivity. We quantified the contribution of intra-specific genetic variability in cassava (Manihot esculenta) and Rhizophagus irregularis to gene reprogramming in symbioses using dual RNA-sequencing. A large number of cassava genes exhibited altered transcriptional responses to the fungus but transcription of most of these plant genes (72%) responded in a different direction or magnitude depending on the plant genotype. Two AMF isolates displayed large differences in their transcription, but the direction and magnitude of the transcriptional responses for a large number of these genes was also strongly influenced by the genotype of the plant host. This indicates that unlike the highly conserved plant genes necessary for the symbiosis establishment, most of the plant and fungal gene transcriptional responses are not conserved and are greatly influenced by plant and fungal genetic differences, even at the within-species level. The transcriptional variability detected allowed us to identify an extensive gene network showing the interplay in plant-fungal reprogramming in the symbiosis. Key genes illustrated that the two organisms jointly program their cytoskeleton organization during growth of the fungus inside roots. Our study reveals that plant and fungal genetic variation has a strong role in shaping the genetic reprograming in response to symbiosis, indicating considerable genotype × genotype interactions in the mycorrhizal symbiosis. Such variation needs to be considered in order to understand the molecular mechanisms between AMF and their plant hosts in natural communities.