Skip to main content

Oak displays common local but specific distant gene regulation responses to different mycorrhizal fungi

Abstract

Background

Associations of tree roots with diverse symbiotic mycorrhizal fungi have distinct effects on whole plant functioning. An untested explanation might be that such effect variability is associated with distinct impacts of different fungi on gene expression in local and distant plant organs. Using a large scale transcriptome sequencing approach, we compared the impact of three ectomycorrhizal (EMF) and one orchid mycorrhizal fungi (OMF) on gene regulation in colonized roots (local), non-colonized roots (short distance) and leaves (long distance) of the Quercus robur clone DF159 with reference to the recently published oak genome. Since different mycorrhizal fungi form symbiosis in a different time span and variable extents of apposition structure development, we sampled inoculated but non-mycorrhizal plants, for which however markedly symbiotic effects have been reported. Local root colonization by the fungi was assessed by fungal transcript analysis.

Results

The EMF induced marked and species specific effects on plant development in the analysed association stage, but the OMF did not. At local level, a common set of plant differentially expressed genes (DEG) was identified with similar patterns of responses to the three EMF, but not to the OMF. Most of these core DEG were down-regulated and correspond to already described but also new functions related to establishment of EMF symbiosis. Analysis of the fungal transcripts of two EMF in highly colonized roots also revealed onset of a symbiosis establishment. In contrast, in the OMF, the DEG were mainly related to plant defence. Already at short distances, high specificities in transcriptomic responses to the four fungi were detected, which were further enhanced at long distance in leaves, where almost no common DEG were found between the treatments. Notably, no correlation between phylogeny of the EMF and gene expression patterns was observed.

Conclusions

Use of clonal oaks allowed us to identify a core transcriptional program in roots colonized by three different EMF, supporting the existence of a common EMF symbiotic pathway. Conversely, the specific responses in non-colonized organs were more closely related to the specific impacts of the different of EMF on plant performance.

Background

The roots of numerous ecologically and economically important forest trees of the Pinaceae, Fagaceae, Betulaceae, Dipterocarpaceae, Myrtaceae and Salicaceae live in symbiosis with highly diverse ectomycorrhizal fungi (EMF) of the Basidiomycota and Ascomycota. In these mutualistic associations, the EMF transfer nutrients to the plant and receive photosynthetically derived sugar [1]. It is thought that EMF evolved repeatedly from saprophytic fungi [2] and that this evolutionary pattern is reflected by important ecological and genetic diversity among the EMF [3]. In the last decades, EMF were classified into morphotypes, some traits of which have been used to predict their ecological role on tree performance [4]. However, many EMF cannot be identified at the species level by anatomical description, and their effects on tree growth are sometimes even strain specific [5]. More importantly, root colonization levels and extent of symbiotic apposition structures (i.e. hyphal mantle and Hartig net) are not always correlated with effects on plant growth [6]. Consequently the anatomy of mycorrhiza or the fungal colonization pattern can be considered as a poor predictor of the functional effect of EMF on plants, and comparative studies should rather focus on the functional significance of the mycorrhizal fungal diversity (e.g. [7]).

Most previous transcriptomic analyses of EM symbioses have focused on genome-sequenced plant and fungal species, e.g. Populus trichocarpa interacting with Laccaria bicolor or Tuber melanosporum [8, 9]. Rapid advances in genomics have prompted new projects to sequence genomes of additional mycorrhizal fungi [3] and host trees, e.g., eucalyptus, spruce, oak, chestnut and pine [10,11,12,13,14]. These efforts are generating important resources for comparing the genomic regulation of plants interacting with different EMF partners, and discovering genes of ecological interest.

Establishment of EM symbiosis involves modifications in the development of both partners, including stimulation of roots formation and growth of fungal mycelium [1, 15, 16]. Some impacts such as those on root formation are induced prior to EM formation (early mycorrhizal stage) [17]. Even at these early stages of interaction, different mycobionts can have highly contrasting effects on growth of a same host plant species [18, 19].

Ectomycorrhizal symbiosis leads to huge modifications of plant and fungal gene expression, as shown by transcriptomic analyses of the partners in various plant/fungus associations [20,21,22,23]. Functional annotation of differentially expressed genes (DEG) has revealed similar plant gene expression patterns in roots colonized by different EMF, particularly for genes involved in plant cell wall modifications and nutrient transport. This suggests that similar plant metabolic pathways were activated during the parallel evolution of EM symbioses involving different fungal clades, although there is no evidence of a common symbiotic signalling pathway in EM associations [24]. However, most of such gene expression studies have focused on local effects within the mycorrhizal roots and few have considered short distance effects (on gene regulation in non-colonized roots) or long distance effects (in leaves) [25,26,27,28], although mycorrhizas influence the physiology of whole plants [29, 30].

The aim of this study was to elucidate how EMF with different evolutionary histories influence gene regulation in Quercus robur L. locally and in plant organs at both short and long distances to the colonized roots. We formulated three hypotheses. First, a core set of plant genes is locally regulated by all inoculated fungi, but a larger set of common core genes directly involved in the mycorrhizal symbiosis is locally regulated in a similar manner in associations with the different EMF. Second, the short and long distance plant responses is more specific to each inoculated EMF species, reflecting the variations in their specific effects on plant growth. Third, the phylogeny of the fungi is, at least locally correlated with changes in gene expression patterns induced in the host plant, i.e. phylogenetically related fungal taxa may induce a more similar plant regulation pattern than more distantly related EMF. To test these hypotheses, we inoculated genetically identical saplings of Q. robur clone DF159 (the TrophinOak platform: www.trophinoak.de) separately with four basidiomycetes: three EMF and one orchid mycorrhizal fungus (OMF). Two of the EMF, Paxillus involutus ATCC200175 and Pisolithus microcarpus 441, are members of the Boletales formally used for molecular studies on mycorrhiza formation with Betula and Eucalyptus, respectively [21, 31]. The other EMF, Laccaria bicolor S238N, a member of the Agaricales, has been widely investigated in EM associations with the model tree species Populus trichocarpa [32]. The three EMF belong to taxa shown to form fully developed mycorrhiza on DF159 [18]. The OMF, Serendipita vermifera MAFF 305830, belongs to the Sebacinales clade B, whose members were originally described as orchid mycorrhizae, but recent DNA studies have shown that they are able to form a broader spectrum of mycorrhiza [33,34,35]. We have found no reports of this OMF forming ectomycorrhizas with oaks [36], but it stimulates growth of several plant species, such as Arabidopsis thaliana, Panicum virgatum and Nicotiana attenuata without forming typical mycorrhizal structures [37,38,39]. Herrmann and Buscot [40] and Frettinger et al. [41] have shown that gene regulation patterns in the oak clone DF159 are largely similar in pre-mycorrhizal roots and fully developed EM. Hence, we compared responses of DF159 oak saplings to inoculation with the four fungi under high humidity conditions to avoid the formation of full mycorrhizas as described by Herrmann et al. [42]. This experimental strategyavoided biases due to comparing mycorrhizas at different stages of differentiation with a non-mycorrhiza forming OMF.

Transcriptomes of the saplings were analysed by Illumina sequencing, reads were aligned against the recently sequenced oak genome [12], and 12 of the differentially expressed transcripts were validated by qRT- PCR. Local, short distance and long distance responses to the fungi were then compared by analysing transcriptomic changes in colonized roots, non-colonized roots and leaves, respectively.

Results

Effects of fungal inoculation on oak growth

After 13 weeks of co-culture, total plant fresh weight and total root length were significantly enhanced by the three EMF P. microcarpus, P. involutus and L. bicolor (Additional file 1: Fig. S1a,c), but not affected by the OMF S. vermifera. The root/shoot ratio (R/S) was increased under the P. microcarpus and P. involutus treatments but not affected by L. bicolor and S. vermifera inoculations (Additional file 1: Fig. S1b). Total leaf area was significantly increased by treatments with P. microcarpus and L. bicolor, but not by P. involutus or S. vermifera treatments (Additional file 1: Fig. S1d).

Oak and fungal read alignments

Illumina RNAseq was performed on colonized roots, non-colonized roots and leaves of control plants (not inoculated) and plants inoculated with the four fungi. Between 26 and 32 million reads were successfully aligned on the oak genome and assigned to genes for colonized root, non-colonized root and leaf samples (Additional file 2: Figure S2).

To assess the degree of fungal mycelium association in colonized (local) and optically non-colonized (short distance) roots of oak, the RNA-seq reads from the two root types were aligned on the genomes of the interacting fungi. In colonized roots, 11% of reads aligned to L. bicolor and 13.5% to P. involutus genome, but only 3.8% to P. microcarpus and 1.2% to S. vermifera. As expected, the alignment rates were markedly lower in non-colonized roots, with 2.9, 1.2, 0.6 and 0.9% of reads for L. bicolor, P. involutus, P. microcarpus, and S. vermifera, respectively.

Fungal symbiosis-related genes are expressed in colonized roots

Ectomycorrhizas were not detected in the densely colonized roots, but the levels of read alignment on L. bicolor and P. involutus genomes represented 11 and 13%, of total reads, respectively. Since this corresponds to more than four million reads per sample, fungal gene expression between colonized roots and free living mycelium was compared. In total 962 genes were up- and 1132 down-regulated in L. bicolor and 499 genes were up- and 348 down-regulated in P. involutus (Additional file 3: Table S1).

Several highly up-regulated fungal genes encoded mycorrhiza-induced small secreted proteins (MiSSPs) that have been implicated in ECM formation (Additional file 3: Table S1 and Additional file 4: Fig. S3). The 40 most up-regulated genes for L. bicolor encoded six MiSSPs: LbMiSSP7, LbMiSSP7.61, LbMiSSP8, LbMiSSP11.4 and LbMiSSP17, with up-regulation levels of 80-fold (LbMiSSP7.61) to more than 7000-fold (LbMiSSP8). Other highly up-regulated genes in L. bicolor included cysteine proteinase inhibitors, ammonium transporters, glutamate dehydrogenase, zinc dependent metalloprotease, major facilitator superfamily transporters and carbohydrate-binding protein. In addition, L. bicolor transcripts encoding carbohydrate active enzymes were up-regulated, including endoglucanases from glycoside hydrolase family 5 (Lb319772) and family 12 (Lb385634, Lb320398 and Lb477020), as well as glycoside hydrolase family 28 polygalacturonase (Lb612983). As observed with L. bicolor, the most highly root contact-induced genes in P. involutus encoded small secreted polypeptides, such as Pi167671 with 105-fold and Pi20703 with 46-fold up-regulation. Genes specifically and strongly up-regulated in P. involutus encoded hydrophobin, protein kinase, phosphatidylserine decarboxylase, acetate transporter, glutaredoxin, thioredoxin disulfide reductase and cytochrome P450 monooxygenase. Carbohydrate active enzyme genes were not up-regulated in P. involutus.

Core DEG in EMF-colonized oak roots

In total, 2252 genes were differentially expressed, relative to controls, under the four fungal treatments: 1081, 1156, 777 and 481 in roots colonized by P. microcarpus, P. involutus, L. bicolor and S. vermifera, respectively (Fig. 1a). A core set of 31 genes was differentially expressed in colonized roots, relative to controls, under all four inoculation treatments (Fig. 1a, Table 1). Most (23) of these core genes were up-regulated, and related to cell-wall organization, including inositol oxygenase and trehalose-phosphate synthase encoding genes. Seven genes of the core set were down-regulated transporter and photosynthesis-related genes. Finally one gene of the core set was a polygalacturonase-encoding gene (Qrob_P0350390.2), which was down-regulated by the EMF, but up-regulated by the OMF (Table 1).

Fig. 1
figure1

Venn diagrams showing numbers of differentially expressed genes in oak after inoculation with mycorrhizal fungi. Differential expression thresholds of RNA sequencing data were at least 2-fold difference to control (no fungus) and a Benjamini-Hochberg adjusted P < 0.01. Numbers in brackets indicate the total number of differentially expressed genes under each treatment. (a) colonized oak roots, (b) non-colonized oak roots, and (c) leaves after inoculation of P. microcarpus, P. involutus, L. bicolor and S. vermifera

Table 1 Thirty-one core genes showing differential expression in colonized roots after inoculation by the mycorrhizal fungi P. microcarpus, P. involutus, L. bicolor and S. vermifera. Genes with at least 2-fold difference relative to fungus-free controls and a Benjamini-Hochberg adjusted P<0.01 were selected

In accordance with our first hypothesis, a stronger common response of 196 core DEG was observed in colonized oak roots when the comparison was restricted to the three EMF treatments (Fig. 1a, Table 2 and Additional file 5: Table S2). Clearly corresponding to a general plant response to EMF, all of these 196 core DEG were regulated in the same direction by the three EMF. Most of them (87%) were down-regulated, including genes encoding proteins involved in carbon metabolism, defense responses, phenolic pathways and transport (Table 2). Due to the well-known impact of biotic interaction on plant hormone equilibrium, phytohormone-related genes were checked for all interactions. Genes involved in auxin biosynthesis, transport and responses were globally affected, relative to controls, by interactions with the three EMF, but not by the OMF S. vermifera. In contrast, ethylene- and gibberellic acid-related genes were differentially regulated in both the EMF- and OMF-treated roots (Additional file 6: Table S3).

Table 2 Forty genes showing the largest differential expression in colonized roots after inoculation by the EMF P. microcarpus, P. involutus and L. bicolor. Genes with at least 2-fold difference relative to fungus-free controls and a Benjamini-Hochberg adjusted P < 0.01 were selected

Fungus-specific gene expression patterns identified in non-colonized oak roots

In non-colonized roots, 3320 genes were identified as differentially expressed, relative to controls, under the four fungal treatments. Notably, the number of DEG strongly varied, depending on the inoculated fungus. 1984 and 1599 genes were differentially expressed in non-colonized roots of plants inoculated with P. involutus and S. vermifera, respectively, but only 87 and 23 in non-colonized roots inoculated with P. microcarpus and L. bicolor, respectively (Fig. 1b). Moreover, only one core gene regulated in non-colonized roots under all four treatments was detected. Transcripts of this gene, Qrob_P0673920.2, encoding a leucine-rich receptor serine threonine kinase, were among the most strongly enriched (log2 fold change between 7.58 and 8.45) in all four interactions (Additional file 7: Table S4). In non-colonized roots, most of the DEG were linked to the interaction with just one fungus, few overlaps between three or even two inoculation treatments were observed. Regarding the EMF, no gene was differentially expressed in non-colonized roots under treatments by all three EMF. In accordance with the second hypothesis, that the fungi would induce species-specific ‘short-distance’ responses in non-colonized roots of host plants, the genes with the highest fold-changes under each of the four treatments had different functions. For example, the treatments with P. microcarpus and P. involutus induced down-regulation of genes related to phytohormone (cytokinin and ethylene) biosynthesis and reduction of ferric iron, respectively (Additional file 7: Table S4). Less specifically, a global response concerning the plant hormone-related genes, particularly genes involved in auxin pathways, was detected in non-colonized roots under both P. involutus and S. vermifera treatments (Additional file 6: Table S3). However, in accordance with the specificity hypothesized, this auxin response was coupled with modifications in gibberellic acid (Gibberellic acid methyltransferase, Gibberellin 20 oxidase, Gibberellin 2-beta-dioxygenase and Ent-kaurene) and ethylene (1-aminocyclopropane-1-carboxylate oxidases) pathways only in the presence of the OMF S. vermifera.

The specificity of responses observed in DEG in non-colonized roots were also reflected in GO-term enrichment patterns (Additional file 8: Table S5), in both the numbers of enriched GO-terms and the directions of regulation (up or down) of genes associated with given GO-terms. No GO-term was enriched among the few DEG associated with the L. bicolor treatment, and only two (photosynthesis, light harvesting and galactose metabolic process) were associated with the P. microcarpus treatment. In contrast, 25 GO-terms were enriched in non-colonized roots after P. involutus inoculation. These included GO-terms related to the regulation of transcription among the up-regulated genes, and GO-terms related to cell-wall organization and defense responses among the down-regulated genes.

In contrast, the S. vermifera treatment led to up-regulation of defense response genes in non-colonized roots, as revealed by enrichment in defense response and response to biotic stimulus GO-terms. Strong expansion of genes related to biotic interactions has been detected in the oak genome, particularly those encoding nucleotide-binding site leucine-rich repeat (NB–LRR) proteins and receptor-like kinases (RLK), for which 95 and 55% expansion rates have been estimated, respectively [12]. Thus, we focused on these groups of potential defense-related genes. Many of these plant disease resistance genes were differentially expressed in non-colonized roots under the fungal treatments. They included 165 of the 1091 (14.8%) NB-LRR and 256 of the 1247 (20.5%) RLK genes (Additional file 9: Fig. S4 and Additional file 10: Table S6). Most of these plant genes responded in a fungus-specific manner. Interestingly, even though the proportion of RLK genes was similar with P. involutus and S. vermifera (5.5 and 5.6%, respectively), under the S. vermifera treatment, genes of LRR-RLK subgroup XIIa, which include receptors known to participate in responses to microbial infection, as well as lectin receptor kinases known for their role in plant defense and immunity, were specifically differentially expressed in non-colonized roots (Additional file 10: Table S6).

Confirmation of fungus-specific long distance responses in leaves

In leaves, 1408 DEG were identified, relative to controls, under the four fungal treatments (Additional file 11: Table S7). For 12 leaf DEG in plants inoculated with L. bicolor, we found a good agreement between estimates of the differential expression levels obtained by qRT-PCR and transcriptome analysis, confirming the overall value of the RNA-seq data (Additional file 12: Fig. S5). Numbers of DEG of the same order of magnitude (235 to 417) were obtained after inoculation with all three EMFs, while the OMF (S. vermifera) treatment resulted in the highest number (653). Only two genes were differentially expressed in leaves under all four fungal treatments (Fig. 1c), one encoding a metacaspase 1 (Qrob_P0407070.2) and the other a superoxide dismutase copper chaperone (Qrob_P0502610.2). Even under the three EMF treatments, only two more common DEG were identified in leaves, one encoding a sucrose synthase (Qrob_P0403240.2) and the other an F-box kelch-repeat protein (Qrob_P0245400.2). Similarly, the only consistent pattern of changes in expression of phytohormone-related genes in leaves under all three EMF treatments was the down-regulation of genes involved in abscisic acid degradation. However higher specificity in responses to the four fungi was observed in expression patterns of the other hormone-related genes (Additional file 6: Table S3). In leaves, the high specificity of changes induced by the four fungi was confirmed by GO-term enrichment analyses, which revealed specific modifications to plant functions induced by each fungus at long distance (Table 3). Laccaria bicolor inoculation induced enrichment of GO-terms associated with plant cell wall organization including lignin catabolic process, microtubule-based process and plant-type cell wall organization. In contrast, the two other EMF inoculations led to enrichment of only one GO-term, photosynthesis and light harvesting with P. involutus and oxidation-reduction process with P. microcarpus. Finally, S. vermifera inoculation induced enrichment of GO-terms related to cell wall biogenesis. In conclusion, the oak displayed highly fungus-specific responses at long distance in the leaves, demonstrating the specificity of the distant impacts of the three EMF on the host tree.

Table 3 GO-term enrichment analysis (Biological processes) performed on leaves for the different treatments. Values correspond to Benjamini-Hochberg adjusted p-values obtained using GOseq [43], in red for up-regulated GO-terms and in blue for down-regulated GO-terms. n.s.: not significant

As characteristic distant plant responses were detected for each fungus, we compared the DEG to mobile mRNAs in A. thaliana listed by [44,45,46,47] to identify whether the long distance effects were correlated with potential mobile mRNAs. Closest homologs in A. thaliana of about 14–23% of the genes differentially expressed in both roots and leaves under each inoculation treatment correspond to putative mobile mRNAs (Additional file 13: Table S8). Up- or down-regulation of most of these putative mobile mRNAs was associated with a specific mycorrhizal fungus, and contributed to the previously described specificity of distant plant responses. Under P. microcarpus and L. bicolor treatments, the DEG among the putative mobile mRNAs were related to cell wall processes, whereas under P. involutus and S. vermifera treatments they related to photosynthesis and defense responses, respectively.

Is fungal phylogeny a good predictor of induced plant regulation patterns?

The phylogeny of the four studied fungi was compared to the changes they induced in the expression patterns. The analysis was based on the DEG data and expression levels of the 5000 most highly expressed genes in each of the three considered plant parts. In colonized roots, the changes induced by the three EMF were more similar to each other than those induced by the OMF S. vermifera, according to heatmap-based clustering of data for both the DEG and 5000 most highly expressed genes (Fig. 2). However, the EMF phylogeny was not congruent with the observed clustering, as heatmaps of changes induced by L. bicolor (Agaricales) and P. involutus (Boletales) clustered together, while the heatmap of changes induced by the other EMF member of the Boletales, P. microcarpus, was more distant. In non-colonized roots, the three EMF induced closer DEG profiles than the OMF, but again the similarity of changes induced by the EMF did not correlate with their phylogenetic closeness. Clustering based on the 5000 most highly expressed genes did not confirm this pattern of responses in non-colonized roots, as the changes in gene expression induced by P. involutus differed most strongly from those induced by the other fungi (two other EMF and an OMF). In leaves, the clustering based on DEG and the 5000 most highly expressed genes showed that the three EMF did not induce more similar changes than the OMF (Fig. 2), confirming thereby the specificity of the distant plant responses to the four tested species of mycorrhizal fungi.

Fig. 2
figure2

Heatmaps of DEG and expression levels of the 5000 most highly expressed genes. Genes are shown in colonized roots, non-colonized roots and leaves

Discussion

Formation of ectomycorrhizas on roots by different EMF requires variable timespans after inoculation, and extent of apposition structures is heterogeneous and their development is not synchronous even on root systems developed in microcosms [18]. Consequently, it cannot be expected that different fungi inoculated on a same plant species will reach similar colonization intensities and rates at a given sampling time, even in a standardized culture system as the one with clonal oak cuttings used in this work [17]. In addition the presented experiment included a Serendipita non able to differentiate EM. For these reasons our experiment deliberately analyzed plants in interactions without EM differentiation, knowing that in the used oak system, inoculation alone exerts typical symbiotic effects on development and physiology [42] as well as on gene regulation [41] of plants, which are similar to the one observed after complete formation of ectomycorrhizas. This functional impact on the plant in absence of symbiotic structures was confirmed for the EMF tested here by their influence on plant growth and fresh weight. Additionally, in this experiment, the different proportions of fungal and plant transcripts found in colonized vs. non-colonized roots allowed us to confirm that at least for two out of the three EMF (P. involutus and L. bicolor), the fungi markedly colonized host roots locally, close to the inoculation points in the microcosms. Despite a weaker colonization revealed by its transcript analysis, the third EMF P. microcarpus had significant effects on plant growth, which have already been described for other fungi [1, 48].

Fungal response to oak roots includes MiSSP and plant cell wall degrading enzyme gene expression

Although sampling of the plant material was performed before mycorrhiza differentiation, L. bicolor genes related to the full establishment of the symbiosis were up-regulated in the colonized roots. Of the mycorrhiza induced small secreted polypeptide (MiSSP) gene products of L. bicolor, MiSSP7 induction leads to alterations of the plant transcriptomic profile and modification of plant immunity during colonization [32], LbMiSSP8 is involved at early stages of ectomycorrhiza formation [49], and MiSSP7.6 has been implicated in the modulation of plant gene expression during root colonization [50]. Interestingly, MiSSP7 is induced by poplar root exudates prior to root colonization process, and the abundance of MiSSP7.6 transcripts increases similarly in ectomycorrhizas and in the extramatricial mycelium close to the roots, indicating that direct physical contact with roots is not required to trigger MiSSP7 or MiSSP7.6 induction [50, 51]. These observations are in line with ours on MiSSP gene induction prior to ectomycorrhiza formation. The abundance of plant cell wall degrading enzyme transcripts of L. bicolor also increased on the colonized roots. The survey of ectomycorrhizal fungal genomes by Kohler et al. [3] demonstrated substantial losses in such genes as compared to saprophytic fungi, but the cell-wall-degrading enzyme genes are expressed in a limited manner during ECM establishment. Zhang et al. [52] showed that Lb319772 endoglucanase, one of the induced genes on colonized oak roots, is involved in plant cell wall modifications during symbiosis formation in poplar.

Different EMF induced a common local core transcriptional program in colonized roots

In accordance with our first hypothesis that a common set of genes related to symbiosis are regulated similarly by different EMF, a much broader range of plant core genes were found to be induced by all three EMF species within colonized roots compared to the number of genes common to the four mycorrhizal treatments. The proportion of DEG that accounted for these ‘core genes’ is similar to proportions observed in chickpea colonized by rhizobia, pathogenic oomycetes or arbuscular mycorrhizal fungi [53]. However, our study revealed several core genes that were not previously identified in other plant-mycorrhizal fungus interaction studies. These include genes encoding inositol oxygenase, an enzyme involved in cell wall biosynthesis [54], which is consistent with its likely participation in establishment of the apposition structure between plant and fungal walls in the Hartig net of mycorrhizae. Most of the core genes that were differentially expressed under all three EMF treatments (including genes encoding components of defense response and phenolic pathways) were down-regulated. This is consistent with the attenuation of cortical host cells’ defenses in roots that promotes fungal penetration [22, 23, 55]. At this local site of interaction, a modification of the auxin pathways was observed in the host plant during interaction with the three EMF but not with the OMF S. vermifera. This is congruent with auxin’s major role in EMF-induced root differentiation [56,57,58]. It is also consistent with increases in root/shoot ratios previously recorded following inoculation of the same oak clone (DF159) with the EMF Piloderma croceum [42], and may therefore be responsible for the increases we observed in root system lengths of plants inoculated with each of the three EMF. In addition, comparison of changes induced by the different EMF partners during interaction with genetically identical oak saplings enabled identification of further locally regulated core genes. These include genes previously reported to be differentially expressed in oak mycorrhizal roots [23] such as genes encoding for laccase, cellulose synthase, polygalacturonase and xyloglucan endotransglucosylase hydrolase, but also other genes with potential role in EM establishment such as a plant disease resistance response protein and a receptor protein kinase encoding gene.

The detected plant core response can reflect a common symbiotic pathway [59] activated by different EMF. This is supported by our finding of fungal symbiosis-related genes activated in L. bicolor and P. involutus hyphae during root colonization. Formation of ECM symbiosis requires lineage-specific genes encoding small secreted proteins that were evidenced in the fungal transcripts of both L. bicolor and P. involutus. Further found genes like hydrophobin, phosphatidylserine decarboxylases, glutaredoxin, and cytochrome P450 monooxygenase have also previously been reported to be up-regulated in ECM tissues [60, 61], confirming that even without differentiation of mycorrhiza a symbiotic plant-fungus dialog was already set up in the colonized roots.

Distant plant responses are highly fungus-specific

The second hypothesis, that short and long distance responses of plants to specific EMF may differ was validated in terms of plant growth responses, numbers of DEG and the metabolic pathways affected. Notably, we observed two patterns of short-distance responses, as few DEG were associated with the P. microcarpus and L. bicolor treatments, but high numbers were associated with the P. involutus and S. vermifera treatments. The smaller short-distance effects of P. microcarpus and L. bicolor on gene expression was not due to a general lack of host plant responses, as these interactors strongly affected gene expression in leaves. Moreover, despite only changing the expression of a few plant genes in non-colonized roots, P. microcarpus and L. bicolor had strong morphological effects, including more than 400 and 250% increases in root length and total plant fresh weight, respectively. Inoculation with P. involutus had strong effects on both host plant growth and gene expression in non-colonized roots. The DEG included genes involved in auxin and jasmonic acid pathways, which are known to play roles in both root growth [62] and formation of ectomycorrhiza and arbuscular mycorrhiza [16, 32, 58, 63]. In contrast, the short-distance responses induced by the non-EMF S. vermifera concern the regulation of numerous genes, particularly up-regulation of genes involved in defense responses, which is consistent with the probably non-symbiotic relationship between this OMF and the oak. As the high humidity conditions of our experiment hamper formation of ectomycorrhiza [42], it was not suitable to explore the potential of the OMF to form EMF. However, contrarily to the three tested EMF species, the OMF had no impact on the plant growth.

Large-scale comparative molecular studies on long-distance plant responses to EMF are scarce. A study focusing on defense-related genes in shoots of Medicago truncatula in presence of three arbuscular mycorrhizal species has been reported [28]. Here we show that more general long-distance plant responses, in leaves, are clearly induced by inoculation with EMF and even with an orchid mycobiont, and that the responses are strongly dependent on the fungal partner. The specific changes observed in expression patterns in leaves included up-regulation of genes related to cell wall modification-related genes following inoculation with L. bicolor, and the identification of putative mobile mRNAs related to cell wall organization among products of DEG detected under treatment with this fungus is consistent with the enhancement of aboveground plant growth. In contrast, no significant increase in leaf area occurred after P. involutus inoculation, but here modification in the regulation concerned genes involved in carbohydrate allocation to roots. Thus, P. involutus seems to influence energy storage, putatively changing the balance between soluble sugar and starch. Induction of an increase in sucrose contents of Populus euphratica leaves by P. involutus has been previously described, even in the apparent absence of mycorrhiza or changes in leaf growth [64]. This suggests that P. involutus has evolved a more specific strategy to benefit from plant carbohydrates than the two other studied EMF.

Particularity of the non-host interaction with S. vermifera

In contrast to the EMF, S. vermifera induced up-regulation of genes encoding chitinases, and the induction of several genes encoding beta-1.3-glucanases, ethylene responsive transcription factors, caffeic acid methyltransferases and flavonol synthases, some of which were also identified as putative mobile mRNAs. The defense responses were detected both at short distance in non-colonized roots and long distance in leaves, and these observations suggest that the OMF induces general defense responses and secondary metabolism pathways, perhaps in a similar manner to the reported induction of oak defenses against the biotrophic pathogen Microsphaera alphitoides induced by a mycorrhiza helper bacterium [65]. The OMF generally form friendly associations with their hosts without eliciting defense responses [66]. However, an interesting outcome of our comparison is that S. vermifera seems to elicit responses in the oak saplings, which are different from their responses to true EMF. The up-regulation of defense-related genes suggests that S. vermifera may induce a non-host response in oak, including up-regulation of large proportions of the NB-LRR and RLKs genes identified in Plomion et al. [12] in non-colonized roots, particularly lectin receptor kinases and LRR-RLKs known to induce defense responses and mediate plant immunity [67, 68]. Such defense response might be related to oak being a non-host of S. vermifera. Although other Sebacina species are ectomycorrhizal in oaks [36], S. vermifera does not colonize oaks in nature, and thus the S. vermifera isolate used in the present study might show an incompatibility reaction to the oak clone.

Fungal phylogeny is not correlated to induced plant responses

Contrary to the third hypothesis, fungal phylogeny was not a good predictor of oak gene expression patterns. As expected, the transcriptomic changes induced by the OMF in the host plants’ colonized roots differed from those induced by the three EMF. However, the differences did not agree with the phylogenic distances between the three EMF, as the patterns of responses to P. involutus (Boletales) and L. bicolor (Agaricales) were more similar to each other than to those induced by the other member of the Boletales, P. microcarpus. Interestingly, we found that the phylogenic distances between the tested fungal taxa were less strongly connected to transcriptomic changes in the distant plant parts than in the infected roots. In this respect, there was a clear distinction between the three EMF and the OMF in infected roots, but not in non-colonized roots and leaves. These results are consistent with the common traits of EMF apposition tissues (i.e. hyphal mantle and Hartig net within EM roots), but sometimes markedly different effects on the growth of the whole host plant, due to specific patterns of distant gene activation. EMF have evolved independently from multiple fungal phyla, and it has been estimated that more than 80 lineages exist [69]. Phylogenetic and genomic analyses suggest that they convergently evolved from white- and brown-rot fungi and from soil saprotrophs [2, 70], and even within the Boletales, the ectomycorrhiza formation capacity of Paxillus and Pisolithus spp. may have emerged independently [3]. If so, clear associations between evolutionary relationships of EMF and the changes they induce in plant expression patterns would only exist among taxa of the same ectomycorrhizal lineage or which are derived from the same ancestor [70].

Conclusion

Our study allowed the first molecular-level comparison of local and distant effects of three EMF and one OMF on different parts of a genetic identical host tree (Q. robur clone DF159), with reference to the recently released oak genome. Through comparative transcriptome analysis, we identified a core set of EMF-induced genes in the oak at the local site of interaction. This supports the existence of common EMF symbiotic pathways [59], despite dissimilar “symbiotic toolkits” of the EMF previously reported [3]. Conversely, in non-colonized roots and leaves, effects were strongly fungus-specific and did not correlate with their phylogenetic relationships. Mycorrhizal inoculation may trigger diverse distant mechanisms like signaling or transfer of metabolites to target organs, which is also supported by the identification of affiliated putative mobile mRNAs known to migrate through the plants. These distant mechanisms vary according to the EMF partner, leading to high heterogeneity of the distant responses, in accordance with evidence that the ectomycorrhizal functional trait evolved independently in distinct fungal lineages [2], resulting in different EMF taxa having different regulatory effects on the same host tree.

Methods

Plant and fungal growth conditions

The pedunculate oak clone DF159 (Quercus robur L.) was micro-propagated and rooted as described by Herrmann et al. [17]. The four used mycorrhizal fungal strains included the ectomycorrhiza forming fungi Pisolithus microcarpus 441, Paxillus involutus ATCC200175 and Laccaria bicolor S238N, as well as the orchid mycorrhiza forming fungus Serendipita vermifera MAFF305830. They were pre-cultivated in Petri dishes on Melin Norkrans Modified agar medium (MMN agar medium [71];) supplemented with 0.1% (w/v) casein hydrolysate (MMNC [17];) in darkness at 20 °C. Fungal inoculum was produced by inoculating a mixture of 750 ml of perlite with 100 ml of MMN liquid medium supplemented with 0.5% (w/v) glucose and 100 ml of a 2–3 weeks old liquid fungal culture produced in MMNC medium, in darkness at 20 °C with shaking at 100 rpm. Rooted micro-cuttings were transferred into closed glass vessels containing 150 ml of the two-week-old inoculum, 150 ml of sterilized perlite and 30 ml of MMN liquid medium with no carbohydrate and tenth-strength P and N (MMN1/10 [17];). Non-inoculated control plants were grown in similar glass vessels with 300 ml of sterilized perlite and 50 ml of MMN1/10 liquid medium. Plants were grown at 23 °C with 16 h photoperiod, with a photosynthetic photon flux density of 180 μmol m− 2 s− 1 during light phases, for 13 weeks.

Morphological analyses of the oak saplings

At harvest, the oaks at bud rest stage A of the endogenous rhythmic growth [25] were gently lifted from the vessels, perlite was carefully removed with forceps and placed on filter paper for photography. Plants measurements were performed on these photos using WinRhizo and WinFolia pro V2005b (Regent Instruments). Plant growth parameters (fresh weight) were evaluated on the complete plant using a minimum of eight replicates (excepted for saplings inoculated with L. bicolor, 6 replicates) and compared using the Kruskal-Wallis test implemented in R software [72]. Between-treatment differences in these parameters were considered significant if P < 0.05.

RNA extractions and Illumina RNA-Seq

For RNA extractions, three types of plant material were harvested: colonized roots represented by short roots with visible, dense (L. bicolor, P. involutus, P. microcarpus) or loose (S. vermifera) fungal mycelium covering their surfaces; non-colonized roots characterized by short roots without any visible fungal mycelium covering their surfaces, and fully developed source leaves belonging to the last developed shoot flush. An absence of mycorrhizas was controlled by de visu observation with dissection microscope on a small set of lateral roots. At harvest, tissues were immediately submerged in liquid nitrogen and the material was stored at − 80 °C. Total RNA was extracted from 50 mg portions of leaf material and 100 mg portions of colonized and non-colonized root material using a MasterPure Plant RNA Purification Kit (Epicentre, Hessisch Oldendorf, Germany). Three biological replicates from single plants were used for non-colonized roots and leaves. For the colonized roots, two (L.bicolor, P. microcarpus and S. vermifera) to three (P. involutus) biological replicates from four to six plants were considered. The extracted RNA was treated with DNase I (Fermentas, St Leon-Rot, Germany), and quantified using a NanoDrop spectrophotometer (Thermo Scientific, Passau, Germany) and a Quant-iT RiboGreen RNA Assay Kit (Invitrogen, Darmstadt, Germany). Its quality was checked using a Nano Chip and Bioanalyzer 2100 (Agilent, Böblingen, Germany). cDNA libraries were prepared from the total RNA using Illumina TruSeq RNA Sample Preparation kit v2 and sequenced using a 2 × 100 bp Illumina HiSeq 2000 platform at the Beijing Genomics Institute, Hong Kong, China. Raw sequence data are available in the Sequence Read Archive at the National Center for Biotechnology Information (BioProject accession PRJNA516042).

Read processing and statistical analysis

Low complexity and low-quality reads were removed with SeqClean http://sourceforge.net/projects/seqclean/files/). Nucleotides with quality score < 20 were removed from the ends of the reads using a custom Java script. Sequences < 50 bp were discarded, as well as sequences lacking paired-end information. Illumina reads were aligned against the haploid oak genome [12] by TopHat2 v.2.1.0 [73] and gene expression levels quantified by featureCounts in the subread package v.1.4.6 [74]. The significance of differences in gene expression was assessed using edgeR v.3.10.2 [75] implemented in the Bioconductor package. Genes were considered differentially expressed when the Benjamini-Hochberg adjusted P-value was less than 1% and the log2-fold change of gene expression was greater than 2. Blast2GO [76] was used to classify the genes in Gene Ontology (GO) terms [77]), and GO enrichment analysis was performed with the Bioconductor package GOseq v.1.20.0 [43], and GO-terms were considered enriched when the Benjamini-Hochberg adjusted P-value was less than 5%. Differentially expressed genes (DEG) were subjected to blastx search against the NCBI nr and Arabidopsis thaliana TAIR databases. Only hits with an E-value of at least <1e-5 were taken into account. Heatmapping was performed using the DEG data and the number of reads of the 5000 most strongly expressed genes using the heatmap function in the R software. The list of the nucleotide-binding site leucine-rich repeat (NB–LRR)-related protein genes and receptor-like kinase (RLK)-encoding genes were extracted from the data published by Plomion et al. [12] and compared to the DEGs identified. For analysis of mobile mRNAs, DEG detected in both roots (colonized and non-colonized) and leaves were compared to the list of 2249 known mobile mRNAs of A. thaliana identified by [44,45,46,47].

To assess how the fungal mycelium contributes to gene expression in colonized and non-colonized roots, reads from the colonized and non-colonized roots samples were aligned using TopHat2 v.2.1.0 [73] on the genomes of the interacting fungi: L. bicolor, P. involutus, P. microcarpus and S. vermifera. Fungal genomes information was retrieved from Joint Genome Institute’s Mycocosm database [78]. Alignment results of colonized roots showed that from the 40 million reads, the highest mapping rates were for L. bicolor (average 4,440,291 reads comprising 11% of input), followed by P. involutus (3,717,004, 9.3%), P. microspermum (1,355,623, 3.3%) and S. vermifera (430,732, 1%). The respective alignments in lateral roots were considerably lower, L. bicolor with 1,001,193 aligned reads comprising 2.5%, P. involutus (424,155, 1%), P. microcarpum (213,166, 0.46%) and S. vermifera (318,539, 0.53%). Fungal gene expression analysis was thus possible only for the colonised root samples with L. bicolor and P. involutus, and for this purpose, we considered the control samples of free living mycelium grown on MMN agar medium by Shah et al. [79] and with following accession numbers: L. bicolor SRR1752511, SRR1752510 and SRR1752509, and P. involutus SRR1752470, SRR1752472 and SRR1752473. The RNA seq analysis pipeline was according to the oak RNA-seq analysis. Illumina reads from colonized roots and from vegetative mycelium were aligned against L. bicolor or P. involutus genome by TopHat2, gene expression levels were quantified by featureCounts, and the significance of differences in gene expression was assessed using edgeR. The thresholds for differential expression were fold change > 2 and Benjamini-Hochberg adjusted P < 0.01. Gene expression analysis was not performed for P. microcarpus and S. vermifera due to the low levels of fungal reads in colonized roots.

Primer design and validation by qRT-PCR

To validate the differential expression of genes revealed by Illumina RNA sequencing, the expression of 12 of those genes was quantified by qRT-PCR in leaf samples from control and L. bicolor-inoculated plants. The same original RNA samples were used for cDNA synthesis for sequencing experiments as for qRT-PCR analyses. Among these 12 genes some were affected by all four fungi (Qrob_P0407070.2 and Qrob_P0502610.2), others by the three EMF (Qrob_P0403240.2), or by three fungi (L. bicolor, S. vermifera and P. involutus or P. microcarpus; Qrob_P0412360.2, Qrob_P0286850.2, Qrob_P0758050.2, Qrob_P0112730.2, Qrob_P0457340.2, Qrob_P0656320.2 and Qrob_P0491740.2), and two specifically by L. bicolor (Qrob_P0135910.2 and Qrob_P0273540.2). Primer pairs were constructed using the oak genome [12] as a reference and tested for functionality, amplicon size, specificity (single peak produced for each primer after melting curve analysis) and efficiency (using standard curve with dilution series of cDNA as template), and the qRT-PCR amplifications were performed according to Tarkka et al. [23]. Sequences of constructed primer pairs are listed in Additional file 14: Table S9. Briefly, using an iScript One-Step RT-PCR Kit with SYBR Green (Bio-Rad) and 18S rRNA as the reference gene, transcript abundances in the leaf samples were determined based on their Ct values using the Relative Expression Software Tool (REST, [80]). The coefficient of variation was used as a reproducibility indicator, with a maximal value of 6.0. Differential gene expression was determined by a randomization test implemented in REST.

Availability of data and materials

Raw sequence data are available in the Sequence Read Archive at the National Center for Biotechnology Information (BioProject accession PRJNA516042).

Abbreviations

DEG:

Differentially expressed gene

ECM:

Ectomycorrhiza

EMF:

Ectomycorrhizal fungus

GO-term:

Gene ontology term

MMN medium:

Melin Norkrans modified medium

MMN1/10:

Melin Norkrans modified medium with 1/10 P and N.

MMNC medium:

Melin Norkrans modified medium supplemented with 0.1% (w/v) casein hydrolysate

NB-LRR protein:

Nucleotide-binding site leucine-rich repeat protein

OMF:

Orchid mycorrhizal fungus

RLK:

Receptor-like kinase

References

  1. 1.

    Smith SE, Read DJ. Mycorrhizal symbiosis. London: Academic Press; 2008.

  2. 2.

    Hibbett DS, Gilbert L-B, Donoghue MJ. Evolutionary instability of ectomycorrhizal symbioses in basidiomycetes. Nature. 2000;407:506–8.

    CAS  PubMed  Article  Google Scholar 

  3. 3.

    Kohler A, Kuo A, Nagy LG, Morin E, Barry KW, Buscot F, et al. Convergent losses of decay mechanisms and rapid turnover of symbiosis genes in mycorrhizal mutualists. Nat Genet. 2015;47:410–5.

    CAS  PubMed  Article  Google Scholar 

  4. 4.

    Agerer R. Exploration types of ectomycorrhizae. Mycorrhiza. 2001;11:107–14.

    Article  Google Scholar 

  5. 5.

    Burgess T, Dell B, Malajczuk N. Variation in mycorrhizal development and growth stimulation by 20 Pisolithus isolates inoculated on to Eucalyptus grandis W. Hill ex Maiden New Phytologist. 1994;127:731-9.

  6. 6.

    Lu X, Malajczuk N, Dell B. Mycorrhiza formation and growth of Eucalyptus globulus seedlings inoculated with spores of various ectomycorrhizal fungi. Mycorrhiza. 1998;8:81–6.

    Article  Google Scholar 

  7. 7.

    Pena R, Polle A. Attributing functions to ectomycorrhizal fungal identities in assemblages for nitrogen acquisition under stress. ISME J. 2014;8:321–30.

    CAS  PubMed  Article  Google Scholar 

  8. 8.

    Plett JM, Tisserant E, Brun A, Morin E, Grigoriev IV, Kuo A, et al. The mutualist Laccaria bicolor expresses a core gene regulon during the colonization of diverse host plants and a variable regulon to counteract host-specific defenses. MPMI. 2014;28:261–73.

    Article  CAS  Google Scholar 

  9. 9.

    Tisserant E, Silva CD, Kohler A, Morin E, Wincker P, Martin F. Deep RNA sequencing improved the structural annotation of the Tuber melanosporum transcriptome. New Phytol. 2011;189:883–91.

    CAS  PubMed  Article  Google Scholar 

  10. 10.

    Myburg AA, Grattapaglia D, Tuskan GA, Hellsten U, Hayes RD, Grimwood J, et al. The genome of Eucalyptus grandis. Nature. 2014;510:356–62.

    CAS  PubMed  Article  Google Scholar 

  11. 11.

    Nystedt B, Street NR, Wetterbom A, Zuccolo A, Lin Y-C, Scofield DG, et al. The Norway spruce genome sequence and conifer genome evolution. Nature. 2013;497:579–84.

    CAS  PubMed  Article  Google Scholar 

  12. 12.

    Plomion C, Aury J-M, Amselem J, Leroy T, Murat F, Duplessis S, et al. Oak genome reveals facets of long lifespan. Nature Plants. 2018;4:440.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Staton M, Zhebentyayeva T, Olukolu B, Fang GC, Nelson D, Carlson JE, et al. Substantial genome synteny preservation among woody angiosperm species: comparative genomics of Chinese chestnut (Castanea mollissima) and plant reference genomes. BMC Genomics. 2015;16:744.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  14. 14.

    Zimin A, Stevens KA, Crepeau MW, Holtz-Morris A, Koriabine M, Marçais G, et al. Sequencing and assembly of the 22-gb loblolly pine genome. Genetics. 2014;196:875–90.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    Barker SJ, Tagu D, Delp G. Regulation of root and fungal morphogenesis in mycorrhizal symbioses. Plant Physiol. 1998;116:1201–7.

    CAS  Article  Google Scholar 

  16. 16.

    Laurans F, Pepin R, Gay G. Fungal auxin overproduction affects the anatomy of Hebeloma cylindrosporumPinus pinaster ectomycorrhizas. Tree Physiol. 2001;21:533–40.

    CAS  PubMed  Article  Google Scholar 

  17. 17.

    Herrmann S, Munch J-C, Buscot F. A gnotobiotic culture system with oak microcuttings to study specific effects of mycobionts on plant morphology before, and in the early phase of, ectomycorrhiza formation by Paxillus involutus and Piloderma croceum. New Phytol. 1998;138:203–12.

    Article  Google Scholar 

  18. 18.

    Buscot F, Herrmann S. At the frontier between basiodiomycotes and plants: reciprocal interactions between mycorrhiza formation and root development in an in vitro system with oaks and hymenomycetes. In: Agerer R, Piepenbring M, Blanz P (eds), Frontiers in Basidiomycote mycology. Eching IHW-Verlag; 2004. pp. 361-76.

  19. 19.

    Sarjala T, Niemi K, Häggman H. Mycorrhiza formation is not needed for early growth induction and growth-related changes in polyamines in scots pine seedlings in vitro. Plant Physiol Biochem. 2010;48:596–601.

    CAS  PubMed  Article  Google Scholar 

  20. 20.

    Duplessis S, Courty P-E, Tagu D, Martin F. Transcript patterns associated with ectomycorrhiza development in Eucalyptus globulus and Pisolithus microcarpus. New Phytol. 2005;165:599–611.

    CAS  PubMed  Article  Google Scholar 

  21. 21.

    Plett JM, Kohler A, Khachane A, Keniry K, Plett KL, Martin F, et al. The effect of elevated carbon dioxide on the interaction between Eucalyptus grandis and diverse isolates of Pisolithus sp. is associated with a complex shift in the root transcriptome. New Phytol. 2015;206:1423–36.

    CAS  PubMed  Article  Google Scholar 

  22. 22.

    Sebastiana M, Vieira B, Lino-Neto T, Monteiro F, Figueiredo A, Sousa L, et al. Oak root response to ectomycorrhizal symbiosis establishment: RNA-Seq derived transcript identification and expression profiling. PLoS One. 2014;9:e98376.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  23. 23.

    Tarkka MT, Herrmann S, Wubet T, Feldhahn L, Recht S, Kurth F, et al. OakContigDF159.1, a reference library for studying differential gene expression in Quercus robur during controlled biotic interactions: use for quantitative transcriptomic profiling of oak roots in ectomycorrhizal symbiosis. New Phytol. 2013;199:529–40.

    CAS  PubMed  Article  Google Scholar 

  24. 24.

    Garcia K, Delaux P-M, Cope KR, Ané J-M. Molecular signals required for the establishment and maintenance of ectomycorrhizal symbioses. New Phytol. 2015;208:79–87.

    PubMed  Article  Google Scholar 

  25. 25.

    Herrmann S, Recht S, Boenn M, Feldhahn L, Angay O, Fleischmann F, et al. Endogenous rhythmic growth in oak trees is regulated by internal clocks rather than resource availability. J Exp Bot. 2015;66:7113-27.

  26. 26.

    Kaling M, Schmidt A, Moritz F, Rosenkranz M, Witting M, Kasper K, et al. Mycorrhiza-triggered transcriptomic and metabolomic networks impinge on herbivore fitness. Plant Physiol. 2018;176:2639–56.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  27. 27.

    Kurth F, Feldhahn L, Bönn M, Herrmann S, Buscot F, Tarkka MT. Large scale transcriptome analysis reveals interplay between development of forest trees and a beneficial mycorrhiza helper bacterium. BMC Genomics. 2015;16:658.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  28. 28.

    Liu J, Maldonado-Mendoza I, Lopez-Meyer M, Cheung F, Town CD, Harrison MJ. Arbuscular mycorrhizal symbiosis is accompanied by local and systemic alterations in gene expression and an increase in disease resistance in the shoots. Plant J. 2007;50:529–44.

    CAS  PubMed  Article  Google Scholar 

  29. 29.

    Loewe A, Einig W, Shi L, Dizengremel P, Hampp R. Mycorrhiza formation and elevated CO2 both increase the capacity for sucrose synthesis in source leaves of spruce and aspen. New Phytol. 2000;145:565–74.

    CAS  Article  Google Scholar 

  30. 30.

    Luo Z-B, Li K, Gai Y, Göbel C, Wildhagen H, Jiang X, et al. The ectomycorrhizal fungus (Paxillus involutus) modulates leaf physiology of poplar towards improved salt tolerance. Environ Exp Bot. 2011;72:304–11.

    CAS  Article  Google Scholar 

  31. 31.

    Wright DP, Johansson T, Le Quéré A, Söderström B, Tunlid A. Spatial patterns of gene expression in the extramatrical mycelium and mycorrhizal root tips formed by the ectomycorrhizal fungus Paxillus involutus in association with birch (Betula pendula) seedlings in soil microcosms. New Phytol. 2005;167:579–96.

    CAS  PubMed  Article  Google Scholar 

  32. 32.

    Plett JM, Daguerre Y, Wittulsky S, Vayssières A, Deveau A, Melton SJ, et al. Effector MiSSP7 of the mutualistic fungus Laccaria bicolor stabilizes the Populus JAZ6 protein and represses jasmonic acid (JA) responsive genes. PNAS. 2014:201322671.

  33. 33.

    Ray P, Craven KD. Sebacina vermifera: a unique root symbiont with vast agronomic potential. World J Microbiol Biotechnol. 2016;32:16.

    PubMed  Article  Google Scholar 

  34. 34.

    Weiss M, Selosse M-A, Rexer K-H, Urban A, Oberwinkler F. Sebacinales: a hitherto overlooked cosm of heterobasidiomycetes with a broad mycorrhizal potential. Mycol Res. 2004;108:1003–10.

    PubMed  Article  Google Scholar 

  35. 35.

    Oberwinkler F, Riess K, Bauer R, Selosse M-A, Weiß M, Garnica S, et al. Enigmatic Sebacinales. Mycol Progress. 2013;12:1–27.

    Article  Google Scholar 

  36. 36.

    Wang Q, He XH, Guo L-D. Ectomycorrhizal fungus communities of Quercus liaotungensis Koidz of different ages in a northern China temperate forest. Mycorrhiza. 2012;22:461–70.

    PubMed  Article  Google Scholar 

  37. 37.

    Barazani O, von Dahl CC, Baldwin IT. Sebacina vermifera promotes the growth and fitness of Nicotiana attenuata by inhibiting ethylene signaling. Plant Physiol. 2007;144:1223–32.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  38. 38.

    Ghimire SR, Charlton ND, Craven KD. The mycorrhizal fungus, Sebacina vermifera, enhances seed germination and biomass production in switchgrass (Panicum virgatum L). Bioenerg Res. 2009;2:51–8.

    Article  Google Scholar 

  39. 39.

    Lahrmann U, Strehmel N, Langen G, Frerigmann H, Leson L, Ding Y, et al. Mutualistic root endophytism is not associated with the reduction of saprotrophic traits and requires a noncompromised plant innate immunity. New Phytol. 2015;207:841–57.

    CAS  PubMed  Article  Google Scholar 

  40. 40.

    Herrmann S, Buscot F. Cross talks at the morphogenetic, physiological and gene regulation levels between the mycobiont Piloderma croceum and oak microcuttings (Quercus robur) during formation of ectomycorrhizas. Phytochemistry. 2007;68:52–67.

    CAS  PubMed  Article  Google Scholar 

  41. 41.

    Frettinger P, Derory J, Herrmann S, Plomion C, Lapeyrie F, Oelmüller R, et al. Transcriptional changes in two types of pre-mycorrhizal roots and in ectomycorrhizas of oak microcuttings inoculated with Piloderma croceum. Planta. 2007;225:331–40.

    CAS  PubMed  Article  Google Scholar 

  42. 42.

    Herrmann S, Oelmüller R, Buscot F. Manipulation of the onset of ectomycorrhiza formation by indole-3-acetic acid, activated charcoal or relative humidity in the association between oak microcuttings and Piloderma croceum: influence on plant development and photosynthesis. J Plant Physiol. 2004;161:509–17.

    CAS  PubMed  Article  Google Scholar 

  43. 43.

    Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11:R14.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  44. 44.

    Deeken R, Ache P, Kajahn I, Klinkenberg J, Bringmann G, Hedrich R. Identification of Arabidopsis thaliana phloem RNAs provides a search criterion for phloem-based transcripts hidden in complex datasets of microarray experiments. Plant J. 2008;55:746–59.

    CAS  PubMed  Article  Google Scholar 

  45. 45.

    Doering-Saad C, Newbury HJ, Couldridge CE, Bale JS, Pritchard J. A phloem-enriched cDNA library from Ricinus: insights into phloem function. J Exp Bot. 2006;57:3183–93.

    CAS  PubMed  Article  Google Scholar 

  46. 46.

    Kanehira A, Yamada K, Iwaya T, Tsuwamoto R, Kasai A, Nakazono M, et al. Apple phloem cells contain some mRNAs transported over long distances. Tree Genet Genomes. 2010;6:635–42.

    Article  Google Scholar 

  47. 47.

    Thieme CJ, Rojas-Triana M, Stecyk E, Schudoma C, Zhang W, Yang L, et al. Endogenous Arabidopsis messenger RNAs transported to distant tissues. Nature Plants. 2015;1:15025.

    CAS  PubMed  Article  Google Scholar 

  48. 48.

    Parladé J, Pera J, Alvarez IF. Inoculation of containerized Pseudotsuga menziesii and Pinus pinaster seedlings with spores of five species of ectomycorrhizal fungi. Mycorrhiza. 1996;6:237–45.

    Article  Google Scholar 

  49. 49.

    Pellegrin C, Daguerre Y, Ruytinx J, Guinet F, Kemppainen M, Frey NF dit, et al. Laccaria bicolor MiSSP8 is a small-secreted protein decisive for the establishment of the ectomycorrhizal symbiosis. Environ Microbiol. 2019;21:3765–3779.

  50. 50.

    Kang H, Chen X, Kemppainen M, Pardo AG, Veneault-Fourrey C, Kohler A, et al. The small secreted effector protein MiSSP7.6 of Laccaria bicolor is required for the establishment of ectomycorrhizal symbiosis. Environ Microbiol. 2020;22:1435–46.

    CAS  PubMed  Article  Google Scholar 

  51. 51.

    Plett JM, Martin F. Poplar root exudates contain compounds that induce the expression of MiSSP7 in Laccaria bicolor. Plant Signal Behav. 2012;7:12–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  52. 52.

    Zhang F, Anasontzis GE, Labourel A, Champion C, Haon M, Kemppainen M, et al. The ectomycorrhizal basidiomycete Laccaria bicolor releases a secreted β-1,4 endoglucanase that plays a key role in symbiosis development. New Phytol. 2018;220:1309–21.

    CAS  PubMed  Article  Google Scholar 

  53. 53.

    Plett JM, Plett KL, Bithell SL, Mitchell C, Moore K, Powell JR, et al. Improved Phytophthora resistance in commercial chickpea (Cicer arietinum) varieties negatively impacts symbiotic gene signalling and symbiotic potential in some varieties. Plant Cell Environ. 2016;39:1858–69.

    CAS  PubMed  Article  Google Scholar 

  54. 54.

    Kanter U, Usadel B, Guerineau F, Li Y, Pauly M, Tenhaken R. The inositol oxygenase gene family of Arabidopsis is involved in the biosynthesis of nucleotide sugar precursors for cell-wall matrix polysaccharides. Planta. 2005;221:243–54.

    CAS  PubMed  Article  Google Scholar 

  55. 55.

    Le Quéré A, Wright DP, Söderström B, Tunlid A, Johansson T. Global patterns of gene regulation associated with the development of ectomycorrhiza between birch (Betula pendula Roth.) and Paxillus involutus (Batsch) Fr. Mol Plant Microbe Interact. 2005;18:659–73.

    PubMed  Article  CAS  Google Scholar 

  56. 56.

    Felten J, Kohler A, Morin E, Bhalerao RP, Palme K, Martin F, et al. The ectomycorrhizal fungus Laccaria bicolor stimulates lateral root formation in poplar and Arabidopsis through auxin transport and signaling. Plant Physiol. 2009;151:1991–2005.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  57. 57.

    Sukumar P, Legué V, Vayssières A, Martin F, Tuskan GA, Kalluri UC. Involvement of auxin pathways in modulating root architecture during beneficial plant–microorganism interactions. Plant Cell Environ. 2013;36:909–19.

    CAS  PubMed  Article  Google Scholar 

  58. 58.

    Vayssières A, Pěnčík A, Felten J, Kohler A, Ljung K, Martin F, et al. Development of the poplar-Laccaria bicolor ectomycorrhiza modifies root auxin metabolism, signaling, and response. Plant Physiol. 2015;169:890–902.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  59. 59.

    Cope KR, Bascaules A, Irving TB, Venkateshwaran M, Maeda J, Garcia K, et al. The ectomycorrhizal fungus Laccaria bicolor produces lipochitooligosaccharides and uses the common symbiosis pathway to colonize Populus roots. Plant Cell. 2019;31:2386–410.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  60. 60.

    Plett JM, Gibon J, Kohler A, Duffy K, Hoegger PJ, Velagapudi R, et al. Phylogenetic, genomic organization and expression analysis of hydrophobin genes in the ectomycorrhizal basidiomycete Laccaria bicolor. Fungal Genet Biol. 2012;49:199–209.

    CAS  PubMed  Article  Google Scholar 

  61. 61.

    Morel M, Jacob C, Kohler A, Johansson T, Martin F, Chalot M, et al. Identification of genes differentially expressed in extraradical mycelium and ectomycorrhizal roots during Paxillus involutus-Betula pendula ectomycorrhizal symbiosis. Appl Environ Microbiol. 2005;71:382–91.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  62. 62.

    Santner A, Estelle M. Recent advances and emerging trends in plant hormone signalling. Nature. 2009;459:1071–8.

    CAS  PubMed  Article  Google Scholar 

  63. 63.

    Hause B, Mrosk C, Isayenkov S, Strack D. Jasmonates in arbuscular mycorrhizal interactions. Phytochemistry. 2007;68:101–10.

    CAS  PubMed  Article  Google Scholar 

  64. 64.

    Luo Z-B, Li K, Jiang X, Polle A. Ectomycorrhizal fungus (Paxillus involutus) and hydrogels affect performance of Populus euphratica exposed to drought stress. Ann For Sci. 2009;66:106p1–10.

    CAS  Article  Google Scholar 

  65. 65.

    Kurth F, Mailänder S, Bönn M, Feldhahn L, Herrmann S, Große I, et al. Streptomyces-induced resistance against oak powdery mildew involves host plant responses in defense, photosynthesis, and secondary metabolism pathways. MPMI. 2014;27:891–900.

    CAS  PubMed  Article  Google Scholar 

  66. 66.

    Perotto S, Rodda M, Benetti A, Sillo F, Ercole E, Rodda M, et al. Gene expression in mycorrhizal orchid protocorms suggests a friendly plant–fungus relationship. Planta. 2014;239:1337–49.

    CAS  PubMed  Article  Google Scholar 

  67. 67.

    Chinchilla D, Zipfel C, Robatzek S, Kemmerling B, Nürnberger T, Jones JDG, et al. A flagellin-induced complex of the receptor FLS2 and BAK1 initiates plant defence. Nature. 2007;448:497–500.

    CAS  PubMed  Article  Google Scholar 

  68. 68.

    Saijo Y, Loo EP, Yasuda S. Pattern recognition receptors and signaling in plant–microbe interactions. Plant J. 2018;93:592–613.

    CAS  PubMed  Article  Google Scholar 

  69. 69.

    Strullu-Derrien C, Selosse M-A, Kenrick P, Martin FM. The origin and evolution of mycorrhizal symbioses: from palaeomycology to phylogenomics. New Phytol. 2018;220:1012–30.

    PubMed  Article  Google Scholar 

  70. 70.

    Martin F, Kohler A, Murat C, Veneault-Fourrey C, Hibbett DS. Unearthing the roots of ectomycorrhizal symbioses. Nat Rev Microbiol. 2016;14:760–73.

    CAS  PubMed  Article  Google Scholar 

  71. 71.

    Marx DH. The influence of ectotrophic mycorrhizal fungi on the resistance of pine roots to pathogenic infections. I. Antagonism of mycorrhizal fungi to root pathogenic fungi and soil bacteria. Phytopathology. 1969;59:153–63.

    Google Scholar 

  72. 72.

    R Core Team. R: A language and environment for statistical computing. 2013.

    Google Scholar 

  73. 73.

    Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  74. 74.

    Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30.

    CAS  PubMed  Article  Google Scholar 

  75. 75.

    Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  76. 76.

    Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, et al. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008;36:3420–35.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  77. 77.

    Harris MA, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, et al. The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res. 2004;32(suppl_1):D258–61.

    CAS  PubMed  Google Scholar 

  78. 78.

    Grigoriev IV, Nikitin R, Haridas S, Kuo A, Ohm R, Otillar R, et al. MycoCosm portal: gearing up for 1000 fungal genomes. Nucleic Acids Res. 2014;42:D699–704.

    CAS  PubMed  Article  Google Scholar 

  79. 79.

    Shah F, Nicolás C, Bentzer J, Ellström M, Smits M, Rineau F, et al. Ectomycorrhizal fungi decompose soil organic matter using oxidative mechanisms adapted from saprotrophic ancestors. New Phytol. 2016;209:1705–19.

    CAS  PubMed  Article  Google Scholar 

  80. 80.

    Pfaffl MW, Horgan GW, Dempfle L. Relative expression software tool (REST©) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002;30:e36.

    PubMed  PubMed Central  Article  Google Scholar 

Download references

Acknowledgements

We thank the Helmholtz Center for Environmental Research (UFZ) for giving us the opportunity to use the infrastructure of the TrophinOak platform (www.trophinoak.de). We also thank Ines Krieg and Barbara Krause for the establishment of the microcosms, and Kerstin Hommel for qRT-PCR analyses.

Funding

We gratefully acknowledge the German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig for the support of the Post Doc position of Marie-Lara Bouffaud in the frame of the Flexpool project 50170649_#3.

Author information

Affiliations

Authors

Contributions

FB and SH conceived the study; SH and MT performed the experiments; MLB, MB and LF analyzed the data; MLB, SH, MT and FB wrote the manuscript. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Sylvie Herrmann.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1:Figure S1. Influence of mycorrhizal fungi on oak growth.

Additional file 2:Figure S2. Numbers of reads obtained from samples under indicated treatments after sequence processing.

Additional file 3:Table S1. List of the differentially expressed genes in colonized roots for L. bicolor and P. involutus.

Additional file 4:Figure S3. Differential L. bicolor and P. involutus gene expression between colonized oak roots and free-living mycelium.

Additional file 5:Table S2. Plant genes differentially expressed after inoculation with the three EMF in colonized roots.

Additional file 6:Table S3. Differentially expressed hormone-related genes identified in colonized roots, non-colonized roots and leaves.

Additional file 7:Table S4. Most strongly differentially expressed plant genes after inoculation with all four fungi in non-colonized roots.

Additional file 8:Table S5. GO-terms identified among differentially expressed genes in non-colonized roots.

Additional file 9:Figure S4. Differential expression of disease-resistance gene families.

Additional file 10:Table S6. Differentially expressed defense-related genes.

Additional file 11:Table S7. Differentially expressed genes in leaves.

Additional file 12:Figure S5. Comparison of Illumina RNA-sequencing and qRT-PCR results for 12 selected genes differentially expressed in leaf samples of L. bicolor.

Additional file 13:Table S8. Putative mobile mRNAs corresponding to genes differentially expressed after inoculation with each of the fungi.

Additional file 14:Table S9. List of the primers used for real-time qRT-PCR analysis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Bouffaud, M., Herrmann, S., Tarkka, M.T. et al. Oak displays common local but specific distant gene regulation responses to different mycorrhizal fungi. BMC Genomics 21, 399 (2020). https://doi.org/10.1186/s12864-020-06806-5

Download citation

Keywords

  • Quercus robur
  • Local and distant effects
  • EMF interaction
  • OMF interaction
  • RNA-seq