Forest Management Type Influences Diversity and Community Composition of Soil Fungi across Temperate Forest Ecosystems

Fungal communities have been shown to be highly sensitive toward shifts in plant diversity and species composition in forest ecosystems. However, little is known about the impact of forest management on fungal diversity and community composition of geographically separated sites. This study examined the effects of four different forest management types on soil fungal communities. These forest management types include age class forests of young managed beech (Fagus sylvatica L.), with beech stands age of approximately 30 years, age class beech stands with an age of approximately 70 years, unmanaged beech stands, and coniferous stands dominated by either pine (Pinus sylvestris L.) or spruce (Picea abies Karst.) which are located in three study sites across Germany. Soil were sampled from 48 study plots and we employed fungal ITS rDNA pyrotag sequencing to assess the soil fungal diversity and community structure. We found that forest management type significantly affects the Shannon diversity of soil fungi and a significant interaction effect of study site and forest management on the fungal operational taxonomic units richness. Consequently distinct fungal communities were detected in the three study sites and within the four forest management types, which were mainly related to the main tree species. Further analysis of the contribution of soil properties revealed that C/N ratio being the most important factor in all the three study sites whereas soil pH was significantly related to the fungal community in two study sites. Functional assignment of the fungal communities indicated that 38% of the observed communities were Ectomycorrhizal fungi (ECM) and their distribution is significantly influenced by the forest management. Soil pH and C/N ratio were found to be the main drivers of the ECM fungal community composition. Additional fungal community similarity analysis revealed the presence of study site and management type specific ECM genera. This study extends our knowledge on the impact of forest management type on general and ectomycorrhizal fungal diversity and community structure in temperate forests. High plasticity across management types but also study site specific spatial distribution revealed new insights in the ECM fungal distribution patterns.

Fungal communities have been shown to be highly sensitive toward shifts in plant diversity and species composition in forest ecosystems. However, little is known about the impact of forest management on fungal diversity and community composition of geographically separated sites. This study examined the effects of four different forest management types on soil fungal communities. These forest management types include age class forests of young managed beech (Fagus sylvatica L.), with beech stands age of approximately 30 years, age class beech stands with an age of approximately 70 years, unmanaged beech stands, and coniferous stands dominated by either pine (Pinus sylvestris L.) or spruce (Picea abies Karst.) which are located in three study sites across Germany. Soil were sampled from 48 study plots and we employed fungal ITS rDNA pyrotag sequencing to assess the soil fungal diversity and community structure. We found that forest management type significantly affects the Shannon diversity of soil fungi and a significant interaction effect of study site and forest management on the fungal operational taxonomic units richness. Consequently distinct fungal communities were detected in the three study sites and within the four forest management types, which were mainly related to the main tree species. Further analysis of the contribution of soil properties revealed that C/N ratio being the most important factor in all the three study sites whereas soil pH was significantly related to the fungal community in two study sites. Functional assignment of the fungal communities indicated that 38% of the observed communities were Ectomycorrhizal fungi (ECM) and their distribution is significantly influenced by the forest management. Soil pH and C/N ratio were found to be the main drivers of the ECM fungal community composition. Additional fungal community similarity analysis revealed the presence of study site and management type specific ECM genera. This study extends our knowledge on the impact of forest management type on general and ectomycorrhizal fungal diversity and community structure in temperate forests. High plasticity across management types but also study site specific spatial distribution revealed new insights in the ECM fungal distribution patterns.

INTRODUCTION
Soils, which are habitats for a wide range of microorganisms including fungi, are known to provide many ecosystem functions (Deacon, 2009). Saprotrophic and ectomycorrhizal (ECM) fungi play especially important roles in decomposition and nutrient cycling (Cairney and Meharg, 2002). Saprotrophic fungi take major part in decay processes (Hobbie et al., 1999), whereas ECM fungi are both decomposers and mutualists with plant roots and can constitute up to 80% of the fungal biomass in forest soils (Nehls, 2008). Soil fungal communities are shaped by a number of biotic and abiotic factors. Because of their close associations with plants, fungi show particularly high sensitivity to shifts in vegetation (Lauber et al., 2008). Trees, as the main plants in forest ecosystems, affect soil fungi passively by shading the ground and by regulating soil temperature and moisture as well as by influencing understory vegetation (Gömöryová et al., 2013). Differences in root architecture and exudation patterns (Aleklett and Hart, 2013;Bakker et al., 2014) are also known to impact soil properties such as nutrient concentrations (Jones et al., 2004) and pH (Hartmann et al., 2008), which directly affect the soil fungal community composition.
Forest management, i.e., any anthropogenic actions linked to timber production, can change the original composition of a forest through the removal and/or replacement of tree species, and by altering the age class structure, exporting biomass and modifying the proportion of dead wood (Paillet et al., 2010). Such changes initiate cascades of consequences that lead to modifications of abiotic and biotic factors both above-and belowground (Paillet et al., 2010). In the natural state, Central Europe would be a deciduous forest (Bengtsson et al., 2000). However, the native vegetation in this area has been altered through management practices (Vanbergen et al., 2005). Currently less than 3% of the area of Central European forests is natural (FAO, 2010). For instance, the introduction of non-native coniferous species like Norway spruce (Picea abies Karst.) and Scots pine (Pinus sylvestris L.) represent one of the major changes in German forest within the last century (BMELV, 2011). The tree monocultures or human-induced fire events are known worldwide to change understory vegetation there by altering forest structures (Wingfield et al., 2001;Royo and Carson, 2006;Grant et al., 2007). In response to management practices, new types of microhabitats emerged which tend to result in a shift in the soil fungal and microbial communities. However, little is known about the effect of forest management regimes on soil fungal communities in Central European forests (Lang and Polle, 2011;Teste et al., 2012;Urbanová et al., 2015).
The German Biodiversity Exploratories, a long-term research platform located at three study sites in Germany, were implemented to assess the impact of forest management (with a focus on timber production) on the biodiversity and functions of forest ecosystems (Fischer et al., 2010). In this frame, previous surveys in the Biodiversity Exploratories investigating the effects of forest management on biodiversity were focusing on invertebrates in soil (Ferlian and Scheu, 2014;Klarner et al., 2014) or litter (Lange et al., 2011;Ott et al., 2014). At the microbial scale these effects of management were studied on soil bacterial communities (Nacke et al., 2011), wood-inhabiting bacteria (Hoppe et al., 2015), and wood-inhabiting fungi (Purahong et al., 2014) as well as soil fungi (Wubet et al., 2012) and yeasts (Yurkov et al., 2011) in beech-dominated forests. However, soil fungal communities under different forest management types within the Exploratories have not been studied so far.
To analyze the impact of four different forest management types within the three Exploratories (study sites ;Fischer et al., 2010;Solly et al., 2014) on soil fungal diversity and community composition, we used the fungal ITS rDNA pyrotag sequencing approach. The main objectives of this study were: (i) to evaluate the impacts of study site and soil environmental conditions on fungal diversity and community composition, (ii) to assess the influence of forest management types on fungal diversity and community structure, and (iii) identify the effects of study site and forest management type on the relative distribution of ECM fungal taxa. We hypothesized that; (1) In accordance with the principle that 'everything is everywhere but the environment selects' (Bass Becking, 1934), different sites and management types will exhibit different levels of diversity and community composition of general and ECM fungi. (2) Forest management, changes in the main tree species and soil properties are important factors shaping the fungal community structure.

Study Sites
The study was performed as part of the Biodiversity Exploratories project, which uses three geographically separated study sites (Fischer et al., 2010). These sites are located in the south west (Swabian Alb), the centre (Hainich-Dün) and the north east (Schorfheide-Chorin) of Germany. Besides the variation in their topogeography they also differ with respect to their geology and climate ( Table 1). We selected four forest management types according to Klarner et al. (2014). Briefly, the types are: age class forests of young managed beech (Fagus sylvatica L.), with beech stands age of approximately 30 years (B30); age class beech stands with an age of approximately 70 years (B70); unmanaged beech stands (unm B), and coniferous stands (Conif) consisting of pine (Pinus sylvestris L.) in Schorfheide-Chorin and spruce (Picea abies Karst.) in Swabian Alb and Hainich-Dün.

Soil Sample Collection
Soil samples were collected from four replicate experimental plots for each of the four forest management types in early May 2011 in parallel soil sampling campaigns conducted in all three study sites. On each of these 48 experimental plots (100 m × 100 m; Supplementary Table S1), 14 soil cores with a diameter of 5 cm were taken along two transects of 40 m length from north to south and from west to east at 1, 7, 13, 19, 31, and 37 m each. Organic layers were removed before taking the soil cores. For our study the upper 10 cm of the 14 soil cores from each plot were mixed. The composite samples were sieved at a mesh size of 2 mm. From the pooled soil samples, aliquots of 50 g were stored at −80 • C for further molecular analysis. Additionally, a 500 g aliquot of soil was dried at 40 • C for basic soil analyses.  (Fischer et al., 2010;Solly et al., 2014); soil properties display the mean values (see section "Basic soil analyses").

Basic Soil Analyses
All soil analyses were performed with air dried <2 mm samples. The pH was determined in duplicate, using a glass electrode in the supernatant of soil suspensions prepared using 1:2.5 mixtures of soil and 0.01 M CaCl 2 . Soil samples were ground to a size of 100 µm and analyzed for total carbon (C) and nitrogen (N) by dry combustion with a CN analyzer "Vario Max" (Elementar Analysensysteme GmbH, Hanau, Germany). After removal of organic C by combustion of samples for 16 h at 450 • C, inorganic C was determined using the same method. Organic C concentrations were then calculated from the differences between total and inorganic C concentrations (Solly et al., 2014).

DNA Extraction, Amplicon Library Preparation and Pyrosequencing
Microbial genomic DNA was extracted from two independent 0.5 g frozen subsamples of each soil sample collected during the sampling campaigns using a MO BIO Power Soil DNA isolation kit (MO BIO Laboratories, Carlsbad, CA, USA) following the manufacturer's protocol. The two soil DNA extracts from each sample were pooled and DNA concentrations were quantified using a NanoDrop UV-Vis spectrophotometer (Peqlab Biotechnologie GmbH, Erlangen, Germany). The fungal ITS rDNA barcode region was amplified using custom ITS1F primers (Gardes and Bruns, 1993) containing Roche 454 pyrosequencing adaptor B and the universal primer ITS4 (White et al., 1990) containing Roche 454 pyrosequencing adaptor A and a samplespecific multiplex identifier (MID). The PCR reactions were performed as described previously (Wubet et al., 2012) in a total volume of 50 µl reaction mix containing 1 µl DNA template (7-15 ng), 25 µl Go Taq Green Master mix (Promega, Mannheim, Germany) and 1 µl of a solution containing 25 pmol of each of the ITS region-specific primers.
All samples were amplified in triplicate and purified using a Qiagen gel extraction kit (Qiagen, Hilden, Germany); DNA concentrations were then measured using a fluorescence spectrophotometer (Cary Eclipse, Agilent Technologies, Waldbronn, Germany) and the samples were pooled to give equimolar representation of each. Unidirectional pyrosequencing from the ITS4 end of the amplicons was performed using a 454 Titanium amplicon sequencing kit and the Roche GS-FLX + 454 pyrosequencer (Roche, Mannheim, Germany) at the Department of Soil Ecology, Helmholtz Centre of Environmental Research (UFZ, Halle, Germany). The raw ITS rDNA sequences were deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under study accession number SRP049544.

Bioinformatic Analyses and Ecological Grouping
Quality filtering and analysis of the 454 ITS sequences was performed in a sequential analysis using mainly MOTHUR (Schloss et al., 2009). In the initial filtering step, sequences with ambiguous bases, homo-polymers and primer differences of more than eight bases were removed. Simultaneously all primer and barcode sequences were discarded. At the same time, sequence reads with a quality score lower than 20 and a read length of less than 300 bp were removed, using the keepfirst 300 bp command and thereby chopping at least 50 bp of the sequence end to remove sequencing noise. This resulted in a sequence read fragment of 300 bp length covering the ITS2 region. All samples were normalized to the smallest sample size (2040 reads per sample) by random removal using the subsample command as implemented in MOTHUR. Sequences were checked for chimeric sequences using the UCHIME algorithm (Edgar et al., 2011) implemented in MOTHUR. The remaining, non-chimeric, sequences were clustered into operational taxonomic units (OTUs) using cd-hit-est (Li and Godzik, 2006) at a threshold of 97% pairwise identity. Taxonomic assignment of the representative sequences for the OTUs was done with the classify.seq command of MOTHUR applied to the UNITE fungal ITS reference database version 6 ( Kõljalg et al., 2013). To improve the taxonomical resolution, those OTUs that had been identified only down to the family level were then subjected to a BLASTn search (e.g., Johnson et al., 2008) against the NCBI GenBank database (Benson et al., 2015). Finally, the fungal OTUs that had been assigned at the genus level were put into ecological groups on the basis of information from literature, in order to link taxonomic information to potential functions.

Statistical Analyses
Statistical analyses were performed using the softwares R version 3.1.1 (R Development Core Team, 2008) and PAST version 2.17b (Hammer et al., 2001).
In order to define the data matrix for our statistical analyses addressing our objectives and verify our hypothesis, we first tested the effect of removing rare fungal taxa on community composition. To assess the influence of rare fungal OTUs (OTUs represented by ≤3 reads), we calculated the nonmetric multidimensional scaling (NMDS) ordination with 20 random starts from the dataset both with all OTUs and with only the abundant fungal OTUs (OTUs represented by >3 reads). The congruence between the two ordination sets was tested by Procrustes correlation analysis using the protest function (Peres-Neto et al., 2006) of the R package vegan (Oksanen et al., 2015) with 999 permutations. We found that fungal community composition was not significantly affected by the presence or absence of rare fungal OTUs (Procrustes correlation coefficient = 0.9985; p < 0.001, suggesting nearly identical ordination). We also tested the need for re-normalization of the abundant fungal OTU data matrix and compared the congruence of the NMDS plots based on the dominant fungal OTU data matrix and on a re-normalized abundant OTU data matrix using Procrustes correlation analysis. We found that this normalization step did not affect the fungal community composition (Procrustes correlation coefficient = 0.9995; p < 0.001). Hence, all subsequent analyses were performed using the relative abundance fungal community matrix excluding singletons, doubletons and tripletons.
Fungal OTU diversity was assessed by calculating the Shannon-Wiener diversity index (Shannon, 1948) using the diversity function in vegan (Oksanen et al., 2015). Effects of study site and forest management type on general and ectomycorrhizal fungal OTU richness and Shannon diversity as well as soil chemical properties were tested using two-way analysis of variances (ANOVA) followed by a Tukey post hoc test. Correlation analyses were performed to test the association between general and ectomycorrhizal fungal OTU richness, respectively, Shannon diversity and soil chemical properties. Whether the correlations were positive or negative was revealed by linear modeling. To assess similarities in the general and ectomycorrhizal fungal community structure among the three study sites, we performed two-way non-parametric multivariate analysis of variances (NPMANOVA) using PAST (Hammer et al., 2001). Relationships between fungal communities and forest management type within the three study sites were visualized using NMDS on the basis of a Bray-Curtis distance matrix and 30 random starts using the metaMDS and ordihull functions of the vegan package (Oksanen et al., 2015). The function envfit was used to test the goodness of fit of environmental parameters on the fungal community NMDS ordination plot based on 999 random permutations. Significant correlations were plotted as vectors.
The relative contributions made by the ecological group of ectomycorrhizal fungal taxa (ECM) to community composition under the different forest management types and different study sites was assessed by means of similarity percentage analysis (SIMPER) based on the relative abundance of the fungal OTUs using PAST (Hammer et al., 2001). This method compares average relative abundance and examines the contribution of each ECM genus to the observed overall dissimilarities between groups or similarities within a given group (Clarke and Warwick, 2001;Gosling et al., 2013). The percentage contributions of the ECM genera, which add up to a total of 90% of the observed dissimilarities between the four different management types, across the three study regions were additionally visualized on heatmaps plotted in R using the package ggplot2 (Wickham, 2009). Significance of differences in the relative abundance of ECM in the respective forest management types were tested by a one-way ANOVA, comparing the mean relative abundance of each genus. To test individual effect of soil properties on ECM genera within the respective forest management types, we performed correlation analysis between the mean abundances with soil pH and C/N.

Sequence Quality Control and Characterization of Soil Fungi
From the 48 soil samples, consisting of four replicates of four different forest management types collected at each of three study sites, a total of 201,381 reads were obtained. Subsequent sequence quality filtering and normalization resulted in 97,920 sequences representing 2040 reads per sample. Further removal of a total of 3113 potential chimeric and non-fungal sequences resulted in 5333 fungal OTUs including 3710 rare OTUs. As described in detail in the section "Material and Methods", neither removal of rare taxa nor re-normalization of the dominant data matrix had a significant effect on the composition of the fungal community (the number of sequences in each dataset is given in Supplementary Table S2). Consequently only the 1623 abundant fungal OTUs were used for subsequent analysis.

Effects of Study Site and Forest Management Type on Fungal OTU Richness, Shannon Diversity and Soil Properties
Analysis of the impact of study site and forest management type on general fungal OTU richness and Shannon diversity using a two-way ANOVA showed that study site had no impact Frontiers in Microbiology | www.frontiersin.org on fungal diversity and fungal richness. In contrast, Shannon diversity was significantly influenced by forest management type. Furthermore, Shannon diversity and OTU richness were affected by the interaction of study site and forest management type ( Table 2). A Tukey post hoc test indicated that the lowest diversity was found in the unmanaged beech forests of the Swabian Alb which was significantly smaller than in young beech forests of the same study site, unmanaged beech stands in Hainich-Dün, and coniferous forests in Schorfheide-Chorin (Supplementary Figure S2). However, we found no significant impact of the forest management type within the study sites on fungal OTU richness. In contrast to total fungal OTU richness, two-way ANOVA showed the effect of study site and forest management type on ECM fungi, the major ecological group in this study ( Table 2). The unmanaged beech and coniferous forest management types in Schorfheide-Chorin showed significantly lower OTU richness than in old and unmanaged beech forests in the two study sites (Swabian Alb and Hainich-Dün, Supplementary Figure S3).
Further two-way ANOVA to assess the effects of study site and forest management type on the soil properties revealed that study site had a significant impact on all six tested soil properties, whereas forest management type affected only C/N ratio ( Table 2). For instance, pH was in general lower in Schorfheide-Chorin. Particularly, the soil in old beech forest in this study site was more acidic than in young and unmanaged beech forest in Swabian Alb and coniferous forest in Hainich-Dün (Supplementary Figure S4). In contrast, in Schorfheide-Chorin all observed forest management types showed a significant higher C/N ratio. An exception are the unmanaged beech forests in this study site because they were not significantly different from the coniferous forests in the other two study sites (Swabian Alb and Hainich-Dün, Supplementary Figure S4). Soil properties did not correlate significantly with either general fungal OTU richness or with the Shannon diversity. Ectomycorrhizal fungal OTU richness and Shannon diversity, however, showed positive correlations with pH and negative correlations with C/N ratio (Supplementary Table S3).

Fungal Community Similarity and Factors Correlating with Community Composition
Comparisons of fungal community similarity among the three study sites using two-way NPMANOVA revealed differences in community composition under the different forest management types (F = 1.885, p < 0.05), study sites (F = 5.774, p < 0.05) and a significant interaction between the two factors (F = 1.481, p < 0.05). Consequently, we performed oneway NPMANOVA to assess the pairwise effect of both study sites and forest management ( Table 3). The results indicated that the fungal community composition at Schorfheide-Chorin differed significantly from those in Swabian Alb and Hainich-Dün (p < 0.05; see Table 3) and that the community composition under conifer differs always from the communities under beech (p < 0.05, Table 3). Similarly, NMDS based ordination of the fungal community displayed the correlation of study site (p < 0.05; Supplementary Figure S5A) and forest management  Figure S5B). Study site based NMDS ordination plot analysis also revealed the correlation with forest management type in the respective study sites (Figure 1). In addition, soil chemical properties also contributed FIGURE 1 | Non-metric multidimensional scaling (NMDS) ordination displaying fungal community composition across the three study sites, (A) Swabian Alb (Alb), (B) Hainich-Dün (Hai), and (C) Schorfheide-Chorin (Sch), in relation to forest management type (o B30-young beech forest, B70-old beech forest, + Conif-coniferous forest, x unm_B-unmanaged beech forest) and soil properties (C tot -total carbon, C org -organic carbon, N tot -total nitrogen, C/N ratio-ratio of organic carbon to total nitrogen and pH). Stress values represent percentage.
significantly to the observed differences in fungal community composition among the three study sites (displayed as vectors in Figures 1A-C and summarized in Supplementary Table S4). Although C/N ratio showed a consistent correlation with the fungal community composition at all three study sites, the effect of pH was detected only at Swabian Alb and Hainich-Dün (Figures 1A,B).
The comparison of the similarity matrices of the whole fungal community with the one of the ECM revealed that the ordination was up to 88% identical (Procrustes correlation coefficient = 0.883, p < 0.05), suggesting that the community composition of ECM fungi followed the pattern of the whole fungal community. The ECM fungal communities in Schorfheide-Chorin differed significantly from those in Swabian Alb and Hainich-Dün (p < 0.05, Supplementary Figure  S5C). A similar trend was found for the effects of forest management with distinct communities under conifers as compared to beech-dominated stands (p < 0.05; Supplementary Figure S5D). These findings were supported by a two-way NPMANOVA, which revealed the effect of study site (F = 5.602, p < 0.05), forest management type (F = 1.744, p < 0.05) and their interaction effect (F = 1.440, p < 0.05) on the ECM fungal community. Indeed, pairwise effects of study site and forest management type followed the same trends than for the general soil fungi ( Table 3). The correlation with soil chemical properties was absent in Schorfheide-Chorin. However, in Swabian Alb pH and C/N ratio corresponded significantly to the ECM fungal community composition. In contrast, at Hainich-Dün all tested parameters except C/N ratio significantly affected the ECM community (Supplementary  Table S5).

Diversity and Community Composition of Ectomycorrhizal Fungal Taxa
Analysis of percentage similarity (SIMPER) and heatmapbased visualization of the ECM community showed the relative importance of the different fungal genera across the three study sites in relation to the four forest management types (overall dissimilarity = 63.07%; Figure 2). We detected 16 ectomycorrhizal genera that contributed to a cumulative percentage of 90% of the overall variation at all three study sites. All the 16 ECM genera appeared in Swabian Alb and Hainich-Dün, while the genera Sebacina and Hydnum were completely absent in Schorfheide-Chorin. The relative abundance of the ECM fungal genera detected in this study also differed among the three study sites. Furthermore, as a result of forest management types, those ECM fungal genera that are of high abundance and/or importance under beech forest management types are less so in coniferous forests, e.g., Genea in Schorfheide-Chorin, and vice versa, e.g., Tylospora in Swabian Alb and Hainich-Dün. Correlation analysis of soil pH and C/N ratio with the ECM genera revealed very few significant relationships (Supplementary Table S6). Nonetheless, we found no clear pattern between ECM features like exploration type or hydrophobicity and soil pH or C/N ratio.
FIGURE 2 | Distribution of mean abundance of the most important ectomycorrhizal genera across the forest management types (young beech forest-B30, old beech forest-B70, coniferous forest-Conif, and unmanaged beech forest-unm B) at the study sites Alb -Swabian Alb, Hai -Hainich-Dün and Sch-Schorfheide-Chorin visualized on a heatmap. Mean abundances based on SIMPER using Bray-Curtis distance, plotting ectomycorrhizal genera contributing to a cumulative percentage of 90%. Different letters within squares indicate significant differences (p < 0.05) based on Tukey post hoc pairwise comparison.

Fungal OTU Characterization
In our study we found that 1263 OTUs (≈78% of all abundant OTUs) belonged to the Asco-or Basidiomycota. Although the phylum of the Ascomycota is the largest among the Fungi (James et al., 2006) we detected more basidiomycetous OTUs in the temperate forest in this study. Of the OTUs identified down to the genus level and classified by ecological function, 38% were ECM fungi. The fungal partners in most ECM symbioses are dominant members of the Basidiomycota followed by Ascomycota (Smith and Read, 2008). The importance of ECM fungi in temperate forest ecosystems may therefore explain the relatively high proportion of basidiomycetous OTUs in our dataset. The primer pairs used in this study have been widely used to amplify the whole fungal ITS rDNA region (Amend et al., 2010;Ihrmark et al., 2012;Wubet et al., 2012). Nevertheless, it is well know that these primer pairs tend to favor Asco-and Basidiomycota over other fungal taxa.

Fungal Diversity Relations with Study Site and Forest Management
We found a significant interaction effect of study site and forest management type on both general fungal Shannon diversity and OTU richness ( Table 2). This is mainly observed within the Swabian Alb young beech forests composed of more diverse fungi than the unmanaged beech forest. We assume that the recurrent disturbances that occur in managed stands (like B30 or Conif) could lead to higher fungal diversity, since ECM and other plant-associated fungi are known to be vulnerable to disturbance (Lazaruk et al., 2005;Gömöryová et al., 2013). In contrast, thinning improves the root growth of the remaining trees (Grant et al., 2007;Lin et al., 2011) and provides new habitats which could lead to an increase in fungal diversity. Thinning processes also alter light conditions, making it possible for other plants to fill the gaps (Brunet et al., 2010). However, plant diversity is not necessarily reflected in below-ground diversity (Gömöryová et al., 2013). We also found that unmanaged beech plots in Swabian Alb were less diverse than plots under highly managed stands (B30). Unmanaged beech forests include several features, such as decaying trees, different forms of deadwood, and uneven relief (mounds and pits; Bengtsson et al., 2000;Paillet et al., 2010), which lead to an expected greater fungal diversity within the soil environment. This could be only shown in Hainich-Dün where the Shannon diversity was among the highest in this forest management type. Recently Mölder et al. (2014) has reported a reduction of forest management leads to dominance of European beech and simultaneously to a decrease in understory vegetation. These authors suggested that dominance of one tree species could maintain an increase in fungal diversity (Paillet et al., 2010;Mölder et al., 2014), which is not thoroughly demonstrated by our results (Supplementary Figure S2). Moreover, in Schorfheide-Chorin the observed differences in the diversity and richness ECM fungal communities could be attributed to the soil properties. Soil pH was significantly and positively correlated to the observed OTU richness of ectomycorrhizal fungal while C/N was correlated significantly and negatively. This is in accordance with previous studies (Wubet et al., 2012;Tedersoo et al., 2014).

Role of the Main Tree Species in Shaping the Fungal Community Composition
The results of NMDS analysis revealed significant variation in the composition of the fungal community among the three study sites and management types within the study sites (Figure 1). In the respective study sites, the general and ectomycorrhizal fungal community under coniferous forest differed significantly from those found in beech-dominated stands. Fungi, in particular mutualists and decomposers, show high sensitivity to changing vegetation (Lauber et al., 2008), implying the effect of tree species as the "most important 'filter' " for shaping the fungal and general microbial community in root surrounding soils (Aleklett and Hart, 2013). This importance is due to various root traits. Deciduous tree species found in temperate forests have root architectures that are very different from those of the coniferous trees (specifically European beech, Norway spruce, and Scots pine) present in these forests (Kutschera and Lichtenegger, 2002). The development of shallow root systems, in comparison to deep rootedness, could explain differences in fungal community structures as seen in our results. Furthermore the presences of only one host species due to the plantation management practice (i.e., Norway spruce -P. abies Karst.), which provides only one kind of litter, could also lead to slight shifts in fungal communities. In general, our results indicated that fungal communities under beech management types did not differ significantly, suggesting no effect of the management types or age class among the beech forests.
Forest management types can influence the vertical stratification of forests. For example, regular timber harvests lead to more compact soil and affect the herb layer (Godefroid and Koedam, 2004). As a consequence, changes in the soil fungal community structure are to be expected. Soil pH, shaping soil fungal communities, can also affect patterns of herbal plants in forests (Brunet et al., 2010) and consequently the soil microbial community. Thus, the lower pH at Schorfheide-Chorin could additionally explain the distinctness of the fungal community as compared to those at the other two study sites. These results on the effect of soil pH are in line with a number of reports on the role of soil pH in shaping fungal communities in soil ecosystems (see review by Aleklett and Hart, 2013). Although pH depends mainly on soil type, texture, and the parental rock material, the plant community and its root deposits can also alter soil acidity (Jones et al., 2004;Lauber et al., 2008). Due to leaching and bioturbation processes, litter influences substrate availability in the topsoil (Kasel et al., 2008) and this may explain why we found a consistent effect of the C/N ratio on fungal community in all the three study sites (Figure 1).

Relationship between the Main Tree Species and ECM Fungal Community
Fungal community similarity analysis revealed the presence of study site and management type specific ECM genera. Due to comparable soil properties ( Table 1), all of the ECM fungal genera were shared between the Hainich-Dün and Swabian Alb study sites. In contrast, in Schorfheide-Chorin only 14 ECM genera were found as members of the ECM communities. Strong variations in the soil properties and organic matter substrates might account to the observed distinct ECM communities (Walker et al., 2014). The abundant ECM genera, Russula, Inocybe, and Lactarius, are known to be widely distributed ECM fungi (Kirk et al., 2008;Tedersoo et al., 2010). Although Russula has mainly been described as being distributed in coniferous forests (Geml et al., 2010;Tedersoo et al., 2010), we were able to show that this genus also occurs in beech dominated deciduous forests, particularly in the case of Schorfheide-Chorin, we found it to be abundant in beech-dominated forests. Consistent with the report of Twieg et al. (2007), who observed significant differences in fungal community composition between stands of different ages in both Douglas-fir and paper birch, we detected the trend of an increase in the genus Russula with increasing forest age in beech stands (Figure 2). In contrast to boreal forests, where the relative abundance of ECM DNA was higher in younger stands , our study showed lower ECM abundance in young temperate beech forests.
In addition, we observed the presence of ECM genera that were specific to coniferous or to beech forest. For instance, in Swabian Alb and Hainich-Dün clear differences were found between the beech and coniferous stands, with the genus Tylospora being significantly more abundant in the coniferous stands. This corresponds to recent findings (Miyamoto et al., 2015) and is in line with previous reports of strong host species preference by ECM fungi (e.g., Ishida et al., 2007;Gao et al., 2013). However, Lactarius, another major genus of the family Russulaceae (Kirk et al., 2008) and one which has already been described in coniferous forest (Pande et al., 2004;Wang, 2007;Ilyas et al., 2013) was almost absent under the coniferous management types in our study. However, this genus was very abundant under all beech management types, particularly at Hainich-Dün and Schorfheide-Chorin. Our results also showed that taxa like Cenococcum, which are commonly found in morphotyping studies in temperate forests (Jany et al., 2003;Goicoechea et al., 2009), were not among the prominent ECM genera at our study sites. Previous studies were able to display relationships between ECM exploration strategies (Agerer, 2001) and nutrient uptake under different soil chemical properties (Suz et al., 2014;Clemmensen et al., 2015;Sterkenburg et al., 2015). Although our results revealed the effect of soil pH on the short-range exploration ECM types Inocybe, Genea, and Hydrophorus in our analysis the significant correlations appear very scattered (Supplementary Table S6B) and do not support the concept of replacements of explorations types with changing soil fertility (Sterkenburg et al., 2015).

CONCLUSION
In general, this study was able to show that forest management types at more or less similar but geographically separated study sites induced a high level of distinctions in the composition of soil fungal communities in general and ECM fungi in particular. High aboveground heterogeneity, like in unmanaged beech forests, does not automatically result in high diversity of soil fungi (Gömöryová et al., 2013). More disturbed sites, such as young beech forests and highly managed coniferous stands, were found to be composed of diverse fungal communities. Dominant tree species in addition with soil properties are the main factors shaping fungal communities in temperate forests. Nevertheless, our analysis of the most important ECM fungi showed that the distribution pattern of certain taxa is more complicated and can only partly be explained by the effects of host tree species, soil pH or C/N ratio. Factors such as understory vegetation, climatic variation, rooting depths and root exudation profiles could also contribute in shaping the fungal communities and therefore should be considered in future studies. Temporal and spatial long-term monitoring of soil fungi in forest ecosystems also remains crucial to discover fungal diversity and ecosystem functioning.