Do Differences in Chemical Composition of Stem and Cap of Amanita muscaria Fruiting Bodies Correlate with Topsoil Type?

Fly agaric (Amanita muscaria) was investigated using a 1H NMR-based metabolomics approach. The caps and stems were studied separately, revealing different metabolic compositions. Additionally, multivariate data analyses of the fungal basidiomata and the type of soil were performed. Compared to the stems, A. muscaria caps exhibited higher concentrations of isoleucine, leucine, valine, alanine, aspartate, asparagine, threonine, lipids (mainly free fatty acids), choline, glycerophosphocholine (GPC), acetate, adenosine, uridine, 4-aminobutyrate, 6-hydroxynicotinate, quinolinate, UDP-carbohydrate and glycerol. Conversely, they exhibited lower concentrations of formate, fumarate, trehalose, α- and β-glucose. Six metabolites, malate, succinate, gluconate, N-acetylated compounds (NAC), tyrosine and phenylalanine, were detected in whole A. muscaria fruiting bodies but did not show significant differences in their levels between caps and stems (P value>0.05 and/or OPLS-DA loading correlation coefficient <0.4). This methodology allowed for the differentiation between the fruiting bodies of A. muscaria from mineral and mineral-organic topsoil. Moreover, the metabolomic approach and multivariate tools enabled to ascribe the basidiomata of fly agaric to the type of topsoil. Obtained results revealed that stems metabolome is more dependent on the topsoil type than caps. The correlation between metabolites and topsoil contents together with its properties exhibited mutual dependences.


Introduction
Furthermore, the possible influence of the collection site and the collection date on the toxin composition, their time stability and the tissue distribution, as the example of A. phalloides carpophores have been also analysed [44,45]. These investigations indicated that environmental factors and mainly the soil type can clearly affect the amatoxin and phallotoxin composition of A. phalloides carpophores, including their histologically different parts.
According to our knowledge, no substantial attempts have been made to relate metabolic fingerprints of A. muscaria to the main structural parts of the mature fruiting body or to determine their correlation with the physico-chemical properties of the microhabitat. Therefore, a metabolomic approach was used to study fly agaric collected in southwestern Poland, with the purposes of improving the knowledge of its chemical composition and determining whether the compounds correlate quantitatively or qualitatively with the topsoil properties or the structural parts of the mature fruiting bodies. The following hypotheses were examined: a) the levels of metabolites in A. muscaria vary with the structure (cap or stem) of the mature fruiting body, b) the soil at the collection site, due to microhabitat factors (chemical element content, acidity, humus and total carbon content), influences the metabolic fingerprint of the A. muscaria fruiting body.

Materials and Methods
Sample sites, mushrooms and soil sampling A. muscaria fruiting bodies were collected from June-October (20010-2011) from nine sites representing various macrohabitats in southwestern Poland (Table 1, Figure 1). Within these sites, representative areas of 30630 m constituted the sampling plots for soil and basidiomata. The collected basidiomata varied in age and size, but very old or damaged specimens were rejected. Fresh basidiomata were cleaned using a plastic knife to remove adherent plant and substrate debris, and the bottom parts of the stems were cut away. Then, they were air-dried with an electric desiccator at 40˚C for 24 h and stored in plastic bags for chemical analysis. Additionally, 27 soil samples (9 sites x 3 replicates each) of the upper forest soil horizon (0-20 cm) were collected after removing the surface layer. The soil substrate samples were stored in plastic bags and air-dried at room temperature in clean conditions for a few weeks for further chemical analysis.

Topsoil analysis
All metals were quantified using ASA 1100B (Perkin-Elmer), except for Mg, which was analyzed using ASA Avanta g (GBC), and Hg, which was analyzed using AMA 254 (ALTEC). The measurements were performed in a commercial certified laboratory according to Polish standards: PB 1 ed.

Sample preparation for NMR spectroscopy
A total of 22 fungi samples were used in this study. The analyzed carpophores represented mature, fully-developed specimens -comparable by their development stage and sizes. Each fruiting body was divided into cap and stem sections, resulting in 44 paired samples. Weighed tissue (Table S1) was transferred into a steal bead homogenizer (Tissuelyser LT; QIAGEN, Hilden, Germany) and disrupted for 10 minutes with shaking at 50 Hz. The resulting powder was suspended in 1 mL of 0.25 M phosphate buffer solution (PBS), pH57.4, containing 10% D 2 O, 3 mM sodium azide and 1 mM sodium salt of trimethylsilyl-2,29,3,39-tetradeuteropropionic acid (TSP), which was used as an internal standard. The solutions were again shaken for 10 minutes at the same frequency and placed in an ultrasound bath at ambient temperature for 10 minutes. To remove cellular debris and substances that are insoluble in PBS,

H NMR measurements
All NMR spectra were recorded with an NMR spectrometer (Bruker Biospin Avance II; Bruker, GmBH, Germany) operating at a proton frequency of 600.58 MHz. 1 H NMR Carr-Purcell-Meiboom-Gill (CPMG) pulse sequence with water suppression was utilized to filter out broad spectral resonances that could have arisen from protein content (number of loops 580, spin echo delay 5400 ms). The spectra were acquired at 300 K after the samples were kept inside the spectrometer for at least 3 minutes. For each sample, 128 scans were collected, resulting in 64 k data points and a spectral width of 20.01 ppm. All spectra were manually phased and baseline-corrected using Topspin 1.3 software (Bruker, GmBH, Germany) with reference to the TSP signal (d50 ppm). Signal assignments were performed based on in-house and online (HMDB, BMRB) databases and confirmed using two-dimensional NMR experiments, 1 H-1 H COSY and 1 H-1 H TOCSY.

Multivariate data analysis
Prior multivariate data analysis spectra were preprocessed. Data reduction was performed by binning the spectra into 0.001 ppm-width buckets (region 0.5-  10.0 ppm), leading to 9,251 variables after excluding the residual water signal region (4.70-4.95 ppm). The obtained dataset (X matrix) was normalized by dividing each spectral bucket by the corresponding sample mass, and therefore, signal areas should reflect metabolite abundances relative to the fungal weight. Such a transformed data matrix was the input for SIMCA-P+ software (v 13.0, Umetrics, Umeå, Sweden), where principal component analysis (PCA), partial least-squares discriminant analysis (PLS-DA) and orthogonal partial least-squares discriminant analysis (OPLS-DA) were performed. The data were scaled using Pareto scaling.

Statistical analysis
Univariate statistical analysis was performed using a Mann-Whitney-Wilcoxon (MWW) test. All statistical results were evaluated at a50.05. In addition to the P value from the MWW test, an absolute value of the OPLS-DA correlation loading (p(corr).0.4) was used as a criterion for significance. Statistical computations employed the Statistica software package (v 10, StatSoft, Tulsa, USA). The Spearman rank correlation coefficient was used to observe relationships between soil factors and metabolite contents in each structure of the basidiomata (stems and caps). The heat maps were generated using MATLAB.

Ethics statement
No specific permits were required for the field studies apart from that from The Main Pharmaceutical Inspector of the Republic of Poland for the collecting and processing of psychotropic chemicals. This permit was obtained for the described field and laboratory studies (permission number: GIF-N-P/4420/6/09). The field of this studies did not involve endangered or protected species.

Discrimination between stems and caps of A. muscaria basidiocarp
Data reduction and visualization were performed using PCA (two principal components; R 2 X50.725; Q 2 Y50.671). A score plot of the first two principal components (PCs) showed a clear grouping of the A. muscaria fruiting body parts, namely caps and stems ( Figure 2a). The caps were much more homogeneous group than stems. The first principal component (PC1) corresponded to 51.2% of the total variance and was responsible for the separation observed between stems and caps. In contrast, the second principal component (PC2) corresponded to 21.3% of the total variance and could be assigned to intra-group variability, mainly in stems. Loadings plot ( Figure 2b) revealed that lipids, choline and glycerophosphocholine (GPC) mainly contributed to PC1, while glucose and trehalose were greatest constituents of PC2.
To identify the regions of the fingerprint that are the most different between stems and caps, the discriminant variant of OPLS analysis was utilized. Again, as in the case of PCA, the model showed full separation between the analyzed groups and good model parameters (R 2 X50.774, R 2 Y50.938 and Q 2 Y50.904) ( Figure 3a). The OPLS-DA results were validated by checking cross-validated score plots (internal cross validation CV1) (Figure 3b), 999 permutation test ( Figure 3c) and crossvalidated Y predicted values (Figure 3d), which again showed full separation of the examined samples and provided evidence that the model is stable. Loading plots were used to evaluate the influence of particular metabolites on the model's discrimination ability. Generally, carbohydrates were present at higher concentrations in stems, while an aliphatic region revealed many metabolites of higher concentrations in caps (Figure 4a). Metabolites that were identified and further tested using univariate statistical analysis are listed in Table 2. A. muscaria caps were characterized by higher concentrations of several amino acids (AAs) including isoleucine, leucine, valine, alanine, aspartate, asparagine, and threonine, as well as lipids, glycerol, choline and GPC compared to stems. Elevated levels of acetate, adenosine, uridine, 4-aminobutyrate, 6-hydroxynicotinate, qinolinate and UDPcarbohydrate were also observed in caps. Conversely, the stems had higher concentrations of formate and fumarate, and carbohydrates, namely aand bglucose and trehalose. Six metabolites: malate, succinate, gluconate, N-acetylated compounds (NAC), tyrosine and phenylalanine, were detected in whole A. muscaria fruiting bodies but did not shown significant differences between caps and stems (P value.0.05) and/or OPLS-DA loading correlation coefficients (p(corr) ,0.4). A simplified representation of metabolite distributions in A. muscaria basidocarps is presented in Figure 4b.
The influence of topsoil on the metabolic composition of A. muscaria basidiocarp The samples used in this study were collected from two types of topsoil, mineral and organic-mineral. Sixteen mushroom samples were taken from mineral topsoil, while   six samples were collected from organic-mineral topsoil. Because A. muscaria stems and caps exhibit different metabolic compositions, the influence of the topsoil was studied separately in whole fruiting bodies, stems and caps. All three OPLS-DA models showed separation between samples originating from different topsoil types (mineral vs. organic-mineral) and were characterized with the following parameters: whole fruiting bodies (R 2 X50.924, R 2 Y50.944 and Q 2 Y50.790) (Figure 5a), stems (R 2 X50.814, R 2 Y50.830, Q 2 Y50.728) (Figure 5b), and caps (R 2 X50.608 R 2 Y50.732, Q 2 Y50.554) (Figure 5c). A general analysis of the corresponding loading plots ( Figure 5) showed that the topsoil type dependence is of a different origin in stems and caps. Therefore, the identified metabolites were tested using univariate statistics, and the results for all three comparisons are reported in Table 3. Nine metabolites were found in similar concentrations in the two types of topsoil namely: alanine, acetate, succinate, formate, UDP-carbohydrate, 6hydroxynicotinate, quinolinate, choline and lipids. A. muscaria basidocarps from mineral topsoil exhibited lower concentration of isoleucine, leucine, valine, asparagine, tyrosine, phenylalanine, NACs, uridine, adenosine and gluconate in comparison to basidocarps from organic-mineral topsoil. Only two metabolites, malate and trehalose, were present at higher concentrations in mineral topsoil compared to organic-mineral topsoil. There was only one compound -glycerolamong the identified metabolites that changed significantly in caps (P value ,0.05 and OPLS-DA correlation coefficient p(corr).0.4) but not in stems. In contrast, stems exhibited a strong dependence on the soil type, showing significantly different levels of seven metabolites that correlated with the topsoil type, whereas there were no differences in the levels of these metabolites in caps that correlated with the topsoil type. The metabolites at significantly lower levels in A. muscaria stems in mineral topsoil compared to organic-mineral topsoil were aspartate, threonine, 4-aminobutyrate, aand b-glucose and GPC, while only the level of fumarate was elevated.

Correlation analysis of topsoil composition and A. muscaria metabolites
To discover which topsoil properties influence the metabolic composition of A. muscaria fruiting bodies, a correlation analysis was performed. To better evaluate the biological meaning of the results, the selected topsoil properties were analyzed using the following order: toxic heavy metals, metalloids, non-metals, essential metals, pH and type of soil. The results were visualized using a heat map approach ( Figure 6). Clustering was applied only for columns (metabolites) since rows (topsoil properties) remained in established order.
In caps clustering revealed evident correlation patterns. Metabolites such as: bglucose, tyrosine, NAC, adenosine, gluconate, leucine, isoleucine, valine, phenylalanine, uridine, glycerol, 4-aminobutyrate, asparagine and succinate exhibited positive correlations with majority of soil parameters (with exception of pH and total Zn). Contrary, malate, UDP-carbohydrate, choline, trehalose, 6hydroxynicotinate, quinolinate and a-glucose showed negative correlations regarding to topsoil parameters with exception of P 2 O 5 , K 2 O, Mn contents and soil pH.  The differences in percentages show the higher or lower levels of metabolites in A. muscaria basidocarps growing on mineral topsoil compared with those growing on mineral-organic topsoil. On the contrary in stem the correlation patterns were less evident. For example, different correlation patterns were observed for caps and stems with respect to aand b-glucose. In stems both glucose forms were clustered together, while in caps presented different behavior. In stems band a-glucose, GPC, isoleucine, tyrosine, leucine, phenylalanine, uridine, adenosine, gluconate and asparagine had exhibited mostly positive correlations with heavy metals and boron. On the other hand, NAC, valine, threonine and lipids showed positive correlations with Cd, Pb and B. Moreover total carbon revealed the same pattern as humus in both parts of mushrooms. Those topsoil properties demonstrated more positive correlations in case of stems than caps.

Discussion
The different distributions of metabolites in A. muscaria basidocarp caps and stems reflect the diverse functions of these parts. In general, stems serve two purposes. First, they support the cap in spore dispersal, and second, they allow the flow of water and solutes from the mycelium growing underneath the ground to the developing fruiting body [47,48]. Our results shows that stems can fulfill storage functions, being reservoirs of energetic metabolites (trehalose, aand bglucose) and the accompanying TCA cycle-related compounds malate and fumarate. This result confirms the common knowledge that disaccharides and polysaccharides are especially good carbon sources for the production of fruiting bodies [47,49]. Formate, which is present in stems, can be used in many biochemical pathways, e.g., puryvate, glyoxylate and dicarboxylate metabolism.
Most of a mushroom's energy is directed into cap development. Therefore, the elevated levels in almost all detected amino acids (except tyrosine and phenylalanine) in caps are likely due to an increased protein biosynthesis (anabolism and catabolism). Lipids are necessary to form cell membranes. Higher amounts of choline, GPC and glycerol fulfill the demand for lipid biosynthesis. Cholines are efficient methyl group carriers that can be utilized for lipid chain elongation. Furthermore, glycerol is a backbone of various lipids; thus, GPC most likely serves as a source of both glycerol and phosphate for phospholipid biosynthesis. Finally, acetate is a basic compound that transfers methyl groups by the formation of acetyl-CoA; therefore, its increased level in caps suggests intensified lipid biosynthesis. In addition, the higher concentrations of uridine and adenosine may reflect the higher reproductive rates of cap cells compared to stems in mature A. muscaria basidocarps. All of these compounds are necessary to maintain the functions of the cap as well as the entire fruiting body [50]. The presence of quinolinate molecules is especially interesting. These compounds might be responsible for the properties of the defined hallucinogenic substances, i.e., ibotenic acids and muscimol, that are responsible for the neurotoxic effect of A. muscaria [51].
The stronger dispersion observed for stem samples in the PCA score plot (Figure 2a) may have occurred for several reasons. One possibility is that there are different PBS extraction efficiencies in stems and caps. However, it is more likely that stems are more sensitive to external (physiological) stimuli than caps. This result may be attributable to the environmental factors, i.e., the type of soil and metal ions, that influence stems.
Indeed, seven metabolites were found at significantly different levels between the two topsoil types only in stems. The concentrations of aand b-glucose in A. muscaria stems depended on the topsoil type, which is likely related to total carbon and humus content. Furthermore, those changes were related to the levels of aspartate, threonine, 4-aminobutyrate, fumarate and GPC.
Fungi cannot obtain energy directly from the sun but rather use the energy stored in decomposing plant and animal biomass to create their own mass. Interestingly, in topsoil with lower total carbon and humus contents, both stems and caps exhibited much higher trehalose concentrations (200% and 455%, respectively) than in A. muscaria from organic-mineral topsoil. This result may indicate that when basic nutrients are not available, mushrooms switch metabolism strategies and store energy in the form of abiotic stress modulatortrehalose.
Because fungi require certain metallic and non-metallic elements to meet their nutritional requirements [52,53], metal ions in the soil may influence their biochemical pathways [53,54]. In the present study, we observed such an indirect influence. Therefore the strongest effect on the metabolic composition of A. muscaria basidocarp must come from the topsoil humus and total carbon content, we observed many strong correlations between mushroom metabolites and inorganic topsoil elements. The concentrations of diverse metal ions in various mushrooms are described in the literature [55][56][57][58][59][60][61][62]. Mushrooms may bioaccumulate some toxic elements, which are considered as the main source of environmental pollution [63].
Our study cannot determine which biochemical pathways are influenced by the metal and non-metal ions present in topsoil. However, it shows that this effect most likely exists. Thus, it is important to examine metabolic studies in a wider context, including studies on nutrient sources, especially when studying mushrooms. Our results clearly show that in discovering new mushroom natural products or studying their metabolic pathways, topsoil properties should always be considered.

Conclusions
This metabolomic study on Amanita muscaria revealed different metabolite contents depending on the morphological part of the fruiting body. Moreover, it is possible to distinguish mushrooms depending on the topsoil type (mineral vs. organic-mineral). Stems were found to be more sensitive to topsoil composition than caps. Correlations between topsoil elements and particular metabolites were found, showing the usefulness of holistic metabolomic approaches in the delineation and characterization of mushroom microenvironments and suggesting that topsoil mineral compositions influence the metabolic pathways of A. muscaria. In conclusion, in the discovery of new mushroom natural products or the study of mushroom metabolic pathways, topsoil properties should always be considered. Further studies in different regions should be performed using this species to determine the metal ion content in the soil as well as in mushroom fruiting bodies and to examine these data and with a metabolomics-based comprehensive approach.