Piriformospora indica Stimulates Root Metabolism of Arabidopsis thaliana

Piriformospora indica is a root-colonizing fungus, which interacts with a variety of plants including Arabidopsis thaliana. This interaction has been considered as mutualistic leading to growth promotion of the host. So far, only indolic glucosinolates and phytohormones have been identified as key players. In a comprehensive non-targeted metabolite profiling study, we analyzed Arabidopsis thaliana’s roots, root exudates, and leaves of inoculated and non-inoculated plants by ultra performance liquid chromatography/electrospray ionization quadrupole-time-of-flight mass spectrometry (UPLC/(ESI)-QTOFMS) and gas chromatography/electron ionization quadrupole mass spectrometry (GC/EI-QMS), and identified further biomarkers. Among them, the concentration of nucleosides, dipeptides, oligolignols, and glucosinolate degradation products was affected in the exudates. In the root profiles, nearly all metabolite levels increased upon co-cultivation, like carbohydrates, organic acids, amino acids, glucosinolates, oligolignols, and flavonoids. In the leaf profiles, we detected by far less significant changes. We only observed an increased concentration of organic acids, carbohydrates, ascorbate, glucosinolates and hydroxycinnamic acids, and a decreased concentration of nitrogen-rich amino acids in inoculated plants. These findings contribute to the understanding of symbiotic interactions between plant roots and fungi of the order of Sebacinales and are a valid source for follow-up mechanistic studies, because these symbioses are particular and clearly different from interactions of roots with mycorrhizal fungi or dark septate endophytes


Introduction
Piriformospora indica is a root-interacting endophytic fungus and has been found in the Indian Thar-Dessert [1]. It belongs to the order of Sebacinaceous (Basidiomycota) [2] and yields a growth promotion effect with various crop plants such as barley, tobacco, maize, and tomato, but also with the model plant Arabidopsis thaliana [3]. Previous studies showed that P. indica promotes nutrient uptake and helps plants to survive under biotic (pathogenic organisms) [4,5] and abiotic (water, temperature, salt, toxins, heavy metal ions) stress conditions [6,7]. Furthermore, it stimulates plant growth, biomass, and seed production [8,9]. The fungus colonizes the epidermal and rhizodermal part of the roots and forms pearshaped spores, which accumulate within the roots and on the root surface. P. indica grows inter-and intracellularly [10] but does not invade the endodermis and aerial parts of the plant.

Results and Discussion
To study the interaction of P. indica with A. thaliana, a sterile hydroponic cultivation system was developed, which allows for the simultaneous profiling of roots and their exudates (Supplementary Figure S1). For this purpose, P. indica was precultivated on agar plates and A. thaliana on agar-filled, bottom-cut PCR tubes protruding into a liquid culture medium. After two weeks, both organisms were brought together in half-strength Murashige-Skoog (MS) medium supplemented with 0.5% sucrose (w/v) and Gamborg B5 vitamins such as myo-inositol, nicotinic acid, pyridoxin, and thiamine. According to our preliminary studies both components are essential for this symbiosis and hence the growth promotion effect of A. thaliana.

P. indica Promotes Shoot Growth of A. thaliana under Specific Culturing Conditions in a Hydroponic System after Root Colonization
If both components (sucrose and Gamborg B5 vitamins) were supplied for the co-cultivation studies, the shoot biomass increased by 22% (p = 4.2ˆ10´5) in inoculated samples compared to the control confirming the previously reported growth promotion effect in soil [31]. Although P. indica colonizes the roots, the root biomass did not change after two weeks of co-cultivation ( Figure 1) leading to the assumption that P. indica might provoke a systemic effect in A. thaliana. Although previous studies have shown a growth promotion effect in roots [23,31], we anticipated slight deviations in a hydroponic system compared to soil due to different physicochemical properties.

P. indica Promotes Shoot Growth of A. thaliana under Specific Culturing Conditions in a Hydroponic System after Root Colonization
If both components (sucrose and Gamborg B5 vitamins) were supplied for the co-cultivation studies, the shoot biomass increased by 22% (p = 4.2 × 10 −5 ) in inoculated samples compared to the control confirming the previously reported growth promotion effect in soil [31]. Although P. indica colonizes the roots, the root biomass did not change after two weeks of co-cultivation ( Figure 1) leading to the assumption that P. indica might provoke a systemic effect in A. thaliana. Although previous studies have shown a growth promotion effect in roots [23,31], we anticipated slight deviations in a hydroponic system compared to soil due to different physicochemical properties. To investigate how P. indica interacts with the host in a hydroponic system, fluorescence microscopy images were recorded using green fluorescent protein (GFP)-labeled P. indica and the interaction monitored at 14 dpi. P. indica grows inter-and intracellularly in root cells of A. thaliana when co-cultivated in soil [31]. In order to test if P. indica still forms fungal hyphae at the root surface under hydroponic conditions, a GFP-labeled P. indica strain was used to visualize colonization. Only weak autofluorescence signals were detected in the non-inoculated roots and roots inoculated with the ) after co-cultivation with P. indica (P.i.) in a hydroponic system. A. thaliana was co-cultivated for two weeks with an agar plug containing mycelia of P. indica. For control A. thaliana was solely cultivated with an agar plug in 0.5ˆMurashige & Skoog (MS) medium supplemented with 0.5% sucrose (w/v) and vitamins: (A) shoot fresh weight (FW); (B) root fresh weight (FW). Values represent the mean˘SD (standard deviation) of three independent experiments (control samples: n = 3ˆ(3´5) and co-cultivated samples: n = 3ˆ5). Each replicate n comprises a pool of 24 plants. Significance analysis of differences between control and co-cultivated samples was performed by t-test: ***, p ď 0.001.
To investigate how P. indica interacts with the host in a hydroponic system, fluorescence microscopy images were recorded using green fluorescent protein (GFP)-labeled P. indica and the interaction monitored at 14 dpi.

P. indica Promotes Shoot Growth of A. thaliana under Specific Culturing Conditions in a Hydroponic System after Root Colonization
If both components (sucrose and Gamborg B5 vitamins) were supplied for the co-cultivation studies, the shoot biomass increased by 22% (p = 4.2 × 10 −5 ) in inoculated samples compared to the control confirming the previously reported growth promotion effect in soil [31]. Although P. indica colonizes the roots, the root biomass did not change after two weeks of co-cultivation ( Figure 1) leading to the assumption that P. indica might provoke a systemic effect in A. thaliana. Although previous studies have shown a growth promotion effect in roots [23,31], we anticipated slight deviations in a hydroponic system compared to soil due to different physicochemical properties. To investigate how P. indica interacts with the host in a hydroponic system, fluorescence microscopy images were recorded using green fluorescent protein (GFP)-labeled P. indica and the interaction monitored at 14 dpi. P. indica grows inter-and intracellularly in root cells of A. thaliana when co-cultivated in soil [31]. In order to test if P. indica still forms fungal hyphae at the root surface under hydroponic conditions, a GFP-labeled P. indica strain was used to visualize colonization. Only weak autofluorescence signals were detected in the non-inoculated roots and roots inoculated with the P. indica grows inter-and intracellularly in root cells of A. thaliana when co-cultivated in soil [31]. In order to test if P. indica still forms fungal hyphae at the root surface under hydroponic conditions, a GFP-labeled P. indica strain was used to visualize colonization. Only weak autofluorescence signals were detected in the non-inoculated roots and roots inoculated with the non-labeled P. indica strain (Supplementary Figure S2). In contrast, roots inoculated with the GFP-labeled strain exhibited a very strong fluorescence already after a 3 s exposure time showing that P. indica colonizes the root surface and penetrates the root of A. thaliana ( Figure 2). Interestingly, P. indica was predominantly detected in lateral roots. According to these observations, we concluded that P. indica colonizes roots of A. thaliana and as a consequence likely leads to changes in root and shoot metabolism. So far, only indolic glucosinolates and hormones have been discussed in depth [24,26].

P. indica Alters the Exudation of Secondary Metabolites by A. thaliana Roots
Hormonal regulation has been described to accompany the colonization of P. indica on A. thaliana roots [22][23][24][32][33][34][35]. An enrichment analysis (Table 1 and Supplementary Table S1) of the upregulated root transcripts 14 dpi as published in Lahrmann et al. [24] revealed an overrepresentation of genes involved in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway "Biosynthesis of plant hormones" (ath01070). As shown in Supplementary Figure S3, P. indica significantly affects phytohormone levels in root exudates and roots, respectively. A higher concentration of hormones was found in exudates of co-cultivated samples as compared to control samples. This effect was in particular pronounced for jasmonate (JA), and jasmonyl-isoleucine (JA-Ile), both showing a more than 10-fold increase in the exudates and its potential role was discussed in reference [24]. In roots, the hormone content was also increased, but to a lower extent for JA, and JA-Ile, for which only a two-to four-fold increase was observed. 12-oxo-phytodienoic acid (OPDA), the precursor of JA and JA-Ile, also accumulated in roots but could not be detected in exudates irrespective of the conditions.
Besides the transcriptional regulation of hormone biosynthesis, hormone responses were overrepresented biological processes in Gene Ontology (GO:0009753 response to jasmonic acid stimulus, GO:0009751 response to salicylic acid stimulus). The analysis of Gene Ontology (GO) terms (Supplementary Table S1) further pointed to the involvement of secondary metabolic processes as the top two enriched processes (GO:0019748) ranked after defense response (GO:0006952). Consequently, roots and their exudates were comprehensively profiled for changes in primary and secondary metabolism upon P. indica colonization. Root exudates were only profiled for changes in semipolar metabolism, as in a screen for primary metabolites (GC/MS) all blank samples (samples without plant and/or fungus) already exhibited a considerable number of primary metabolites. Representative base peak chromatograms are depicted in Figure 3 and reveal a unique metabolic fingerprint for both ionization modes, ESI(+) and ESI(´). A principal component analysis (PCA) could verify this assumption. For both ionization modes 89% of the total variance was explained by the first principal component (PC1) and 3% ESI(+) as well as 4% ESI(´) by PC2 (Supplementary Figure S4).  Non-targeted UPLC/ESI(+/−)-QTOFMS-based metabolite profiling revealed that the concentration of 200 out of 341 detected ESI(+) components as well as 271 out of 377 ESI(−) components was significantly affected (p < 0.01) due to the presence of P. indica. A total of 28 (ESI(+)) as well as 24 (ESI(−)) components were down-and 172 (ESI(+)) as well as 247 (ESI(−)) components were upregulated due to the inoculation implying that P. indica stimulates root exudation of A. thaliana.
As already observed by Lahrmann et al. [24], the amount of compounds associated with nucleoside and aromatic amino acid metabolism was reduced in concentration by the inoculation, while that of aliphatic and indolic glucosinolate metabolism (except for 4-hydroxy-indole-3carbaldehyde), dihydroxybenzoic acid (DHBA) conjugates, JA metabolism as well as fatty acid and pantothenic acid metabolism was increased. A number of phenylpropanoids including coumarins and oligolignols (except for scopoletin and G(8-O-4)FA sulfate) showed reduced levels in the exudates of inoculated samples ( Figure 4) leading to the assumption that these oligomers are further metabolized inside the cell and not exuded, very likely to oligolignols or to lignin [36], a main constituent of the cell wall. Both, glucosinolate (ath00966) and phenylpropanoid biosynthesis (ath00940), were among the overrepresented KEGG pathways of root transcripts (Table 1).  Non-targeted UPLC/ESI(+/´)-QTOFMS-based metabolite profiling revealed that the concentration of 200 out of 341 detected ESI(+) components as well as 271 out of 377 ESI(´) components was significantly affected (p < 0.01) due to the presence of P. indica. A total of 28 (ESI(+)) as well as 24 (ESI(´)) components were down-and 172 (ESI(+)) as well as 247 (ESI(´)) components were upregulated due to the inoculation implying that P. indica stimulates root exudation of A. thaliana.
As already observed by Lahrmann et al. [24], the amount of compounds associated with nucleoside and aromatic amino acid metabolism was reduced in concentration by the inoculation, while that of aliphatic and indolic glucosinolate metabolism (except for 4-hydroxy-indole-3-carbaldehyde), dihydroxybenzoic acid (DHBA) conjugates, JA metabolism as well as fatty acid and pantothenic acid metabolism was increased. A number of phenylpropanoids including coumarins and oligolignols (except for scopoletin and G(8-O-4)FA sulfate) showed reduced levels in the exudates of inoculated samples ( Figure 4) leading to the assumption that these oligomers are further metabolized inside the cell and not exuded, very likely to oligolignols or to lignin [36], a main constituent of the cell wall. Both, glucosinolate (ath00966) and phenylpropanoid biosynthesis (ath00940), were among the overrepresented KEGG pathways of root transcripts (Table 1).
Nicotinic acid, an important precursor for vitamin B6, and thus, key player in the photoprotection of plants [37], also decreased in concentration upon co-cultivation in the root exudates ( Figure 4). Obviously, nicotinic acid is required by P. indica. If this compound was not supplemented, no growth-promoting effect was observed of the host.
In the exudates we also detected differences in the dipeptide pool, namely the concentration of Phe-Gly, Ile-Leu, Phe-Ile and Ile-Phe was enhanced, while that of Ile-Val, Leu-Val, Val-Leu, Leu-Pro and Leu-Tyr was reduced in the co-cultivated samples ( Figure 4). These differences might originate from different functionalities of the respective dipeptides and require further investigation. So far, Komarova et al. [38] showed that peptide transporters (AtPTR5 and AtPTR5) facilitate the uptake of nitrogen from the rhizosphere.

Changes in the Root Metabolism of A. thaliana
The secondary metabolic changes detected in root exudates, especially that of glucosinolate biosynthesis, phenylpropanoid biosynthesis, and phenylalanine metabolism should also be reflected in root metabolism. In addition, transcriptionally enriched KEGG pathways of primary metabolism (Table 1), such as nitrogen metabolism (ath00910), glycine, serine and threonine metabolism (ath00260), and cyanoamino acid metabolism (ath00460) were expected in the GC/MS-based metabolite profiles.  Nicotinic acid, an important precursor for vitamin B6, and thus, key player in the photoprotection of plants [37], also decreased in concentration upon co-cultivation in the root exudates ( Figure 4). Obviously, nicotinic acid is required by P. indica. If this compound was not supplemented, no growth-promoting effect was observed of the host.
In the exudates we also detected differences in the dipeptide pool, namely the concentration of Phe-Gly, Ile-Leu, Phe-Ile and Ile-Phe was enhanced, while that of Ile-Val, Leu-Val, Val-Leu, Leu-Pro and Leu-Tyr was reduced in the co-cultivated samples ( Figure 4). These differences might originate from different functionalities of the respective dipeptides and require further investigation. So far, Komarova et al. [38] showed that peptide transporters (AtPTR5 and AtPTR5) facilitate the uptake of nitrogen from the rhizosphere.

Changes in the Root Metabolism of A. thaliana
The secondary metabolic changes detected in root exudates, especially that of glucosinolate biosynthesis, phenylpropanoid biosynthesis, and phenylalanine metabolism should also be reflected in root metabolism. In addition, transcriptionally enriched KEGG pathways of primary metabolism (Table 1), such as nitrogen metabolism (ath00910), glycine, serine and threonine metabolism (ath00260), and cyanoamino acid metabolism (ath00460) were expected in the GC/MS-based metabolite profiles. Figure 5 shows two representative extracted ion chromatograms of m/z 73 (equals C 3 H 9 Si + and is a typical fragment for trimethylsilylated compounds) obtained from a pool of inoculated and non-inoculated roots. Again, the inoculated profile is distinct compared to the non-inoculated root metabolic profile. Forty-eight percent of the total variance was explained by PC1 and 13% by PC2 and are plotted in Supplementary Figure S5. Non-targeted GC/EI-Q-MS based metabolite profiling revealed 287 out of 801 differentially accumulated components. Among them, we detected amino acids (e.g., Asn, Thr, Leu, 3-Cyano-Ala, beta-Ala, Val, Ala, Gln, ornithine, Pro, pyro-Glu, and GABA), organic acids (e.g., citrate, 2-oxoglutarate, fumarate, malate, oxalate, glycerate, fumarate, and 3-hydroxy-3-methylglutaric acid), carbohydrates (e.g., 1-O-methylglucopyranoside, 1-O-methylgalactopyraoside, maltose, raffinose, trehalose, xylose, ribose), polyols (erythritol, myo-inositol), phosphates (e.g., glycerol-3phosphate, phosphate, glycerophosphoglycerol), and sulfates (e.g., sulfate, thiamine, thiamine-hex) belonging to the starch and sucrose metabolism, glycolysis, tricarboxylic acid (TCA) cycle, amino acid metabolism, and urea metabolism. All compounds showed increased levels in the inoculated roots ( Figure 6) except for pyruvate, erythritol, allantoin, and 4-methyl-5thiazoleethanolglycopyranoside (for spectrum see Supplementary Figure S6) indicating that the initially applied amount of sucrose and thiamine is metabolized by P. indica. We observed an increase in the concentration of Asn, Gln, Ser, Thr, and Ala at 14 dpi. Serine was also increased in its Non-targeted GC/EI-Q-MS based metabolite profiling revealed 287 out of 801 differentially accumulated components. Among them, we detected amino acids (e.g., Asn, Thr, Leu, 3-Cyano-Ala, beta-Ala, Val, Ala, Gln, ornithine, Pro, pyro-Glu, and GABA), organic acids (e.g., citrate, 2-oxoglutarate, fumarate, malate, oxalate, glycerate, fumarate, and 3-hydroxy-3-methylglutaric acid), carbohydrates (e.g., 1-O-methylglucopyranoside, 1-O-methylgalactopyraoside, maltose, raffinose, trehalose, xylose, ribose), polyols (erythritol, myo-inositol), phosphates (e.g., glycerol-3-phosphate, phosphate, glycerophosphoglycerol), and sulfates (e.g., sulfate, thiamine, thiamine-hex) belonging to the starch and sucrose metabolism, glycolysis, tricarboxylic acid (TCA) cycle, amino acid metabolism, and urea metabolism. All compounds showed increased levels in the inoculated roots ( Figure 6) except for pyruvate, erythritol, allantoin, and 4-methyl-5-thiazoleethanolglycopyranoside (for spectrum see Supplementary Figure S6) indicating that the initially applied amount of sucrose and thiamine is metabolized by P. indica. We observed an increase in the concentration of Asn, Gln, Ser, Thr, and Ala at 14 dpi. Serine was also increased in its levels as described by Lahrmann et al. [15], but did not pass the defined significance level (p = 0.051). In general, the data collected are in good accordance with the transcriptional changes and lead us to the hypothesis that A. thaliana provides nitrogen to the fungus so that P. indica can maintain a biotrophic nutritional state [15]. In the leaf profiles, the N-rich amino acids (Gln, Arg, Asn, 3-Cyano-Ala, and ornithine) were among the few differentially accumulated compounds decreasing in concentration upon colonization and consequently showed the opposite trend (Supplementary Figure S7) compared to the roots. This raises the question if these amino acids are transported to the root to feed P. indica. Most likely, these amino acids are required to balance the nutritional state of P. indica. To trace the flow of nutrients, further investigations are required. The change in the concentration of organic acids and carbohydrates was comparable for roots and leaves except that less differential changes were observed in the leaf profiles. These results again show that P. indica activates primary root metabolism of A. thaliana. According to our data, both partners appear to offer each other nutrients to maintain a mutualistic interaction, since an enhanced amount of P, S, and N (in the form of amino acids) was observed in roots of A. thaliana colonized with P. indica.

Non-Targeted GC/MS Based Metabolite Profiling Reveals Perturbations in the Primary Root Metabolism
Int. J. Mol. Sci. 2016, 17, 1091 8 of 18 levels as described by Lahrmann et al. [15], but did not pass the defined significance level (p = 0.051).
In general, the data collected are in good accordance with the transcriptional changes and lead us to the hypothesis that A. thaliana provides nitrogen to the fungus so that P. indica can maintain a biotrophic nutritional state [15]. In the leaf profiles, the N-rich amino acids (Gln, Arg, Asn, 3-Cyano-Ala, and ornithine) were among the few differentially accumulated compounds decreasing in concentration upon colonization and consequently showed the opposite trend (Supplementary Figure S7) compared to the roots. This raises the question if these amino acids are transported to the root to feed P. indica. Most likely, these amino acids are required to balance the nutritional state of P. indica. To trace the flow of nutrients, further investigations are required. The change in the concentration of organic acids and carbohydrates was comparable for roots and leaves except that less differential changes were observed in the leaf profiles. These results again show that P. indica activates primary root metabolism of A. thaliana. According to our data, both partners appear to offer each other nutrients to maintain a mutualistic interaction, since an enhanced amount of P, S, and N (in the form of amino acids) was observed in roots of A. thaliana colonized with P. indica.

LC/MS Based Non-Targeted Metabolite Profiling Shows an Induction of Aliphatic and Indolic Glucosinolate Metabolism, Flavonoids, and Oligolignols in Roots
Besides primary metabolism, secondary root metabolism was investigated, since one category "secondary metabolic process" was a highly ranked candidate in the GO enrichment analysis. A unique fingerprint was observed in the root LC/MS profiles (Figure 7). According to Supplementary Figure S5, 76% of the entire variance was explained by PC1 and 0.07% by PC2 for the positive mode. These values were similar for the negative mode (PC1: 74%; PC2: 0.08%) and indicate that secondary metabolism is perturbed to a greater extent than primary metabolism.

LC/MS Based Non-Targeted Metabolite Profiling Shows an Induction of Aliphatic and Indolic Glucosinolate Metabolism, Flavonoids, and Oligolignols in Roots
Besides primary metabolism, secondary root metabolism was investigated, since one category "secondary metabolic process" was a highly ranked candidate in the GO enrichment analysis. A unique fingerprint was observed in the root LC/MS profiles (Figure 7). According to Supplementary Figure S5, 76% of the entire variance was explained by PC1 and 0.07% by PC2 for the positive mode. These values were similar for the negative mode (PC1: 74%; PC2: 0.08%) and indicate that secondary metabolism is perturbed to a greater extent than primary metabolism.

LC/MS Based Non-Targeted Metabolite Profiling Shows an Induction of Aliphatic and Indolic Glucosinolate Metabolism, Flavonoids, and Oligolignols in Roots
Besides primary metabolism, secondary root metabolism was investigated, since one category "secondary metabolic process" was a highly ranked candidate in the GO enrichment analysis. A unique fingerprint was observed in the root LC/MS profiles (Figure 7). According to Supplementary Figure S5, 76% of the entire variance was explained by PC1 and 0.07% by PC2 for the positive mode. These values were similar for the negative mode (PC1: 74%; PC2: 0.08%) and indicate that secondary metabolism is perturbed to a greater extent than primary metabolism.   In these profiles, 167 out of 329 detected compounds (ESI(+)) were altered in abundance and 188 out of 359 for the negative ionization mode due to the presence of P. indica. Similarly to the exudates, a higher number of compounds displayed upregulated abundance in the inoculated samples compared to the non-inoculated samples. From these numbers one can once more conclude that P. indica stimulates secondary root metabolism as well.
In accordance with Lahrmann et al. [24], aliphatic and indolic glucosinolates as well as their breakdown products, aromatic amino acids, coumarins, oligolignols, and flavonoids accumulated in inoculated roots (Figure 8) confirming the transcript data (KEGG, Table 1: glucosinolate ath00966 and phenylpropanoid biosynthesis ath00940). Although the plant seems to be in a defensive stage, no camalexin was detected in these profiles. In the leaf profiles an increased amount of aliphatic and indolic glucosinolates as well as their breakdown products, JA conjugates, oligolignols, and hydroxycinnamic acid amides was detected (Supplementary Figure S8). Several flavonoids (glycosylated kaempferol and quercetin) were only detected as differential in the root profiles and not in the leaf profiles, leading to the conclusion that this substance class plays an important role in the mutualistic interaction of A. thaliana and P. indica. Recently, Lahrmann et al. [24] stated that it remains to be clarified if flavonoids are accumulating in roots of A. thaliana upon interaction with P. indica. Indeed, we show that flavonoids accumulate in roots of A. thaliana upon co-cultivation with P. indica. Most likely, enhanced flavonoid biosynthesis, in addition to JA signaling [39], may also function as a signal for P. indica. In these profiles, 167 out of 329 detected compounds (ESI(+)) were altered in abundance and 188 out of 359 for the negative ionization mode due to the presence of P. indica. Similarly to the exudates, a higher number of compounds displayed upregulated abundance in the inoculated samples compared to the non-inoculated samples. From these numbers one can once more conclude that P. indica stimulates secondary root metabolism as well.
In accordance with Lahrmann et al. [24], aliphatic and indolic glucosinolates as well as their breakdown products, aromatic amino acids, coumarins, oligolignols, and flavonoids accumulated in inoculated roots (Figure 8) confirming the transcript data (KEGG, Table 1: glucosinolate ath00966 and phenylpropanoid biosynthesis ath00940). Although the plant seems to be in a defensive stage, no camalexin was detected in these profiles. In the leaf profiles an increased amount of aliphatic and indolic glucosinolates as well as their breakdown products, JA conjugates, oligolignols, and hydroxycinnamic acid amides was detected (Supplementary Figure S8). Several flavonoids (glycosylated kaempferol and quercetin) were only detected as differential in the root profiles and not in the leaf profiles, leading to the conclusion that this substance class plays an important role in the mutualistic interaction of A. thaliana and P. indica. Recently, Lahrmann et al. [24] stated that it remains to be clarified if flavonoids are accumulating in roots of A. thaliana upon interaction with P. indica. Indeed, we show that flavonoids accumulate in roots of A. thaliana upon co-cultivation with P. indica. Most likely, enhanced flavonoid biosynthesis, in addition to JA signaling [39], may also function as a signal for P. indica.

3.2.
Pre-Cultivation of P. indica P. indica was cultured on agar plates (1.5% (w/v) agar) for 3 weeks at 28˝C in the dark using Aspergillus minimal medium [29]. For this purpose, a punched out agar block with mycelia of P. indica was placed in the center of a culture plate.

Conduction of Co-Cultivation Studies and Production of Plant Material
Co-cultivation studies were performed as previously described [25]. In short, two-week old A. thaliana plantlets were co-cultivated for two weeks with P. indica in a hydroponic system under short day conditions (23˝C, 8 h light, 180 µE¨m´2¨s´1 and 21˝C, 16 h dark). After two weeks of co-cultivation (four-week old plants), the medium containing the nutrient solution and the root exudates was filtered and stored at 4˝C in Schott flasks until further processing. At harvest, roots were cut below the bottom of the PCR tube and blotted dry with a paper towel before shock freezing in liquid nitrogen. Finally, they were stored at´80˝C until further processing. More technical details are visualized in Supplementary Figure S1. Media composition is summarized in Appendix A.

LC/MS-Based Metabolite Profiling
For LC/MS-based metabolite profiling (UPLC: Acquity, Waters, Eschborn, Germany; MS: MicrOTOF-Q I hybrid quadrupole time-of-flight mass spectrometer equipped with an Apollo II electrospray ion source, Bruker Daltonik GmbH, Bremen, Germany), the ground tissue material was processed by solid liquid extraction using methanol/water, 80/20 (v/v) (40 mg root fresh weight corresponds to 200 µL extraction solution and 50 mg leaf fresh weight corresponds to 400 µL extraction solution). Analytes of the nutrient solution were extracted by a reversed-phase solid phase extraction procedure (180 mL medium result in 120 µL analysis solution).

Preparation of Nutrient Solutions for LC/MS Analysis
All exudate samples were prepared and analyzed by UPLC/ESI-QTOFMS as presented in Lahrmann et al. [24]. In short, the nutrient solution was spiked with 20 µM 2-(2,4-dichlorophenoxy) acetic acid, evaporated until dryness, reconstituted in 9 mL water/methanol 95/5 (v/v) and subjected to a Bond Elut PPL cartridge (200 mg, 3 mL, Agilent Technologies, Böblingen, Germany). Finally, the eluate was subjected to a solid phase extraction workup and reconstituted in 120 µL water/methanol 70/30 (v/v) prior to LC/MS analysis. Technical details of the solid phase extraction workup can be found in the Appendix B.

Sample Preparation and Profiling of Tissue Material for LC/MS Analysis
The plant material was processed according to Böttcher et al. [29]. As already described, the frozen material was extracted twice with methanol/water, 80/20 (v/v) and reconstituted in methanol/water, 30/70 (v/v) prior to LC/MS-analysis. More details of the extraction procedure can be found in the Appendix B.

Non-Targeted LC/MS-Based Profiling and Data Analysis
Changes in the secondary plant metabolism were analyzed by UPLC/ESI-QTOFMS. Samples were injected onto an Acquity UPLC system (Waters, Eschborn, Germany), equipped with an HSS T3 column (100ˆ1.0 mm, particle size 1.8 µm, Waters), and separated using a binary gradient (A: water/0.1% (v/v) formic acid; B: acetonitrile/0.1% (v/v) formic acid). Eluting compounds were detected in positive and negative ionization mode from m/z 100-1000 using a MicroTOF-Q I hybrid quadrupole time-of-flight mass spectrometer equipped with an Apollo II electrospray ion source (Bruker Daltonics, Billerica, MA, USA). All instrument parameters and further settings can be found in the Appendix B.
Raw data files were converted to mzData using CompassXPort version 1.3.10 (Bruker Daltonics). For feature detection, alignment, and filling of missing values the R package XCMS version 1.41.0 [40] was used. Settings are summarized in Appendix B.
The intensities of the resulting features (m/z-retention time pairs) were log 2 transformed and subjected to a two-sided Student's t-test. Relevant mass spectral features were extracted within a predefined range (isolation width:˘0.02 m/z) and elemental compositions were calculated applying a default error range (15 ppm). Putative elemental compositions were checked for consistency while analyzing elemental compositions of fragment ions and neutral losses of collision-induced dissociation (CID)-mass spectra. For acquisition of CID mass spectra quasi-molecular cluster ions were isolated at the Q1 (isolation width:˘3 m/z) and fragmented inside the collision cell using argon as collision gas. Product ions were detected as described above. All mass spectral data can be found in the MetaboLights repository (MTBLS341) [41].

Sample Preparation of Tissue Material
One hundred µL extract of the remainder from the LC/MS-based metabolite profiling studies was spiked with 100 µM succinic acid-2,2,3,3-d 4 , dried down in a vacuum concentrator, and stored at 20˝C until further processing.

Preparation of Samples for Non-Targeted Metabolite Profiling and Analysis of GC/MS Profiles
Dried down extracts were subjected to a two-step derivatization process using methoxyamine hydrochloride and N,O-bis(trimethylsilyl)trifluoroacetamide. Derivatized samples were injected splitless at 230˝C onto an Agilent 6890N GC equipped with a split/splitless inlet and a ZB-5 column (30 mˆ0.25 mm, 0.25 µm 95% dimethyl/5% diphenyl polysiloxane film, 10 m integrated guard column, Phenomenex, Aschaffenburg, Germany). Eluting components were detected from m/z 70-600 by using an Agilent 5975 Series Mass Selective Detector (Agilent Technologies, Waldbronn, Germany). For the generation of the metabolite profiles, chromatograms were baseline-corrected using Metalign [42]. Peak intensities above 500 arbitrary ion current units were imported into the TagFinder software [43], aligned using the retention index model of van den Dool and grouped according to their common retention time and mass spectral features. For statistical analysis, peak intensities of cluster (cluster size > 3) were normalized to the internal standard (succinic acid-2,2,3,3-d 4 ). Then, all data were log 2 -transformed and submitted to a two-sided Students t-test. Finally, resulting mass spectral features were identified via best mass spectral and retention index match using the Golm Metabolome Database [44] and the NIST2012 software (May 2011, National Institute of Standards and Technology, Gaithersburg, MD, USA). Details of the derivatization protocol and instrument parameters can be found in the Appendix C.
All statistical analysis was either performed with the R statistical language, the Bioconductor environment, the package pcaMethods or Microsoft Excel.

Hormone Analysis
Hormone profiling was conducted as described in Ziegler et al. [45] (for further information see Appendix D). Root material was homogenized, extracted in methanol, and processed firstly using a hydrophobic solid phase extraction cartridge (Chromabond Sorbent HR-XC, Macherey-Nagel, Düren, Germany) and secondly with an anion exchange solid phase extraction cartridge (Diethylaminoethyl Sephadex (DEAE-Sephadex)). For the root exudates the anion exchange step was omitted.
Analytes were separated by an Agilent 1290 Infinity HPLC system and detected on-line by ESI-MS/MS using an API 3200 triple-quadrupole LC-MS/MS system equipped with an ESI Turbo Ion Spray interface (AB Sciex, Darmstadt, Germany). Triple quadrupole scans were acquired in the multiple reaction monitoring mode (MRM) with Q1 and Q3 set at unit resolution. Scheduled MRM was performed with a window of 90 s and a target scan time of 0.1 s. Selected MRM transitions and compound specific parameters can be found in Ziegler et al. [45].

Microscope Images
Bright-field and fluorescence microscopic images were recorded with a Stemi 2000 Axio Imager stereomicroscope (Carl Zeiss MicroImaging GmbH, Göttingen, Germany). For bright-field images a Plan Apochromat 20ˆ/0.75 objective with 20ˆmagnification was used and for fluorescence images a Plan Apochromat 20ˆ/0.75 objective with 20ˆmagnification, a GFP-Filter 450-490 nm, filterset 9 and the Axio Imager camera.

Transcript Enrichment Analysis
Overrepresentation analysis of the overexpressed genes in Arabidopsis 14 dpi as published in Lahrmann et al. [24] was performed with DAVID [46,47] against the default background genes from TAIR using KEGG pathways [48] and Gene Ontology [49].

Data Availability
All data sets are available from the MetaboLights repository [41] under the accession number MTBLS341.

Conclusions
The mutualistic interaction of P. indica with A. thaliana resulted in an increased shoot biomass production, but not root biomass after a two-week co-cultivation. Interestingly, the presence of P. indica had an obvious effect on the root's primary and secondary metabolism and the exudation rate, but not on leaf metabolism of A. thaliana. Apparently, P. indica stimulates the belowground metabolism of A. thaliana, but not the shoot metabolism. The metabolic changes identified can be considered as potential biomarkers, which need to be tackled in the near future. Previous studies and this study have shown that indolic glucosinolates and hormones are important for the interaction. The induction of the defense response might indicate that the plant tries to balance fungal growth and maintain its mutualism. This assumption could be confirmed by the analysis of appropriate mutants. In the future, new mutants, especially of the flavonoid metabolism, need to be obtained to investigate the mutualistic interaction in more depth. It is possible that plant-growth promoting microorganisms can be valuable tools for crop improvement [7,50], as they promote the plant growth and help the plant to cope with abiotic and biotic stress factors. XCMS settings: Feature detection with the help of the centWave algorithm (sntresh: 3, prefilter: (3.100), ppm: 25, peak width: (5.12); feature alignment with the help of the XCMS function group.density (minfrac: 0.75, bw: 2, mzwid: 0.05); missing values replacement by the XCMS function fillPeaks.
Analysis of raw data: DataAnalysis 4.2 software (Bruker Daltonics) for deconvolution and generation of extracted ion chromatograms
HR-XC: The resin was conditioned with 1 mL methanol followed by 1 mL water (the liquid was passed through SPE 96 well plate (50 mg HR-XC resin per well) by centrifugation at 250ˆg for 5 min using a JS5.3 bucket rotor in an Avanti J-26XP centrifuge (Beckman). Samples were transferred to the SPE 96 well plate, the resin washed with 1 mL H 2 O. Analytes were eluted successively by adding 1 mL MeOH (for acidic hormones) and 1 mL methanolic ammonia (0.35 M) for zeatins. DEAE-Sephadex: The resin was washed with 1 mL methanol. The methanolic eluates from the HR-XC plates were loaded onto DEAE-Sephadex (acetate form, 50 mg¨well´1) filled. After washing with 1 mL methanol, the analytes were eluted with 1 mL of 3 M acetic acid in methanol.
Further processing: eluates were transferred to 2 mL Eppendorf tubes and evaporated to dryness; residues were dissolved in 40 µL of 20% (v/v) methanol, diluted with 40 µL of water and centrifuged at 10,000ˆg for 10 min.