High Fungal Diversity but Low Seasonal Dynamics and Ectomycorrhizal Abundance in a Mountain Beech Forest

Forests on steep slopes constitute a significant proportion of European mountain areas and are important as production and protection forests. This study describes the soil fungal community structure in a European beech-dominated mountain forest stands in the Northern Calcareous Alps and investigates how it is determined by season and soil properties. Samples were collected at high spatial resolution in an area of ca. 100 m × 700 m in May (spring) and August (summer). Illumina MiSeq high-throughput sequencing of the ITS2-region revealed distinct patterns for the soil fungal communities. In contrast to other studies from temperate European beech forest stands, Ascomycota dominated the highly diverse fungal community, while ectomycorrhizal fungi were of lower abundance. Russulaceae, which are often among the dominant ectomycorrhizal fungi associated with European beech, were absent from all samples. Potentially plant pathogenic fungi were more prevalent than previously reported. Only subtle seasonal differences were found between fungal communities in spring and summer. Especially, dominant saprotrophic taxa were largely unaffected by season, while slightly stronger effects were observed for ectomycorrhizal fungi. Soil characteristics like pH and organic carbon content, on the other hand, strongly shaped abundant taxa among the saprotrophic fungal community. Supplementary Information The online version contains supplementary material available at 10.1007/s00248-021-01736-5.


Introduction
Approximately 40% of European mountain areas are covered by forests [1]. In comparison to low-elevation forests, mountain forests are generally exposed to more extreme environmental conditions, such as higher solar radiation and precipitation, lower annual temperatures, harsher winters, and avalanches [2]. Although European mountain forest ecosystems are in a relatively natural state, they currently undergo major compositional and structural shifts-partly as a result of changes in land use, climate, and natural forest disturbance regimes [3]. European beech is the most abundant broadleaved tree species in Central Europe [4] and is a major constituent of Alpine forests. A third of the current potential beech forest area in Europe is located on calcareous soils and considered especially vulnerable to future climate scenarios [5,6]. Mountain forest soils developed on calcareous substrates, in particular from limestone and dolomite, are usually shallow, have high organic carbon and Ca contents, and nearneutral pH values [7,8]. While there is an increasing number of studies on soil processes in spruce or mixed mountain forests (on calcareous bedrock) [9,10], the knowledge base on shallow, calcareous beech forest soils in mountainous areas is still limited.
The steep slopes and varying topographies of mountain forests result in heterogeneity in microsite and microclimatic conditions, as well as interlocking patterns of different humus and soil types [11]. This causes also heterogeneity in ecosystem processes such as N mineralization, nitrification, and consequently N availability [12]. High heterogeneity in site properties and stand characteristics provides a variety of ecological niches and promotes a high (aboveground) biodiversity [13,14]. Calcareous forests in the European Alps are hotspots for vascular plant diversity [15]. In contrast to vascular plant diversity, which has been studied across all European forest types [15], studies of soil microbial diversity have been carried out mainly on homogeneous low-elevation forests [16][17][18] (but see [10]). These studies have primarily considered fertility gradients or the effects of understory vegetation [16][17][18]. Although the importance of spatial heterogeneity for microbial community composition has been widely recognized [19], it is not fully understood whether observations at large scales apply to standlevels [20], particularly for calcareous mountain forests with high soil spatial heterogeneity [21].
Soil fungi are a key component of microbial communities in forest ecosystems, especially as integral constituents of carbon (C) and nutrient cycling [22,23]. DNA sequence-based methods triggered manifold developments in fungal taxonomy [e.g., 24] and enabled recent surveys of soil fungi from a range of forests-revealing an extremely high soil fungal diversity [e.g., 25,26]. Three major functional groups ("guilds") of soil fungi can be distinguished according to their ecological lifestyle [27], i.e., saprotrophic, symbiotic, and plant pathogenic fungi. Symbiotic ectomycorrhizal (ECM) fungi acquire carbohydrates from living host plants [28] and often have a strong ecosystem-and host-specificity [25,26]. Their abundance, diversity, and activity can be strongly related to seasonal changes in photosynthetic activity in accordance with the phenological stages of the host plant [29][30][31][32]. Soil calcium concentration has been reported to influence the taxonomic richness of ECM symbionts on Alnus spp. [33]. Saprotrophic fungi rely for their nutrition entirely on dead organic matter and are therefore especially sensitive to quality and quantity of organic matter [34,35], but see [36] but were also reported to be affected by potential evapotranspiration [37]. In saprotrophic fungi, mycelia growth peaked in autumn and winter [10]. Plant pathogenic fungi in soil make up the bulk of microbes attacking roots [38] and are dominantly nonhost-specific necrotrophs, i.e., killing host tissue and thereby providing conditions favorable to their growth and sporulation. Štursová et al. recently found no indication of seasonal patterns of plant pathogenic fungal activity below ground in a mountainous, Picea abies forest [10].
Relative to the functions that soil fungi perform, there is a very limited understanding of the spatial and temporal variations of soil fungal communities and their driving forces [32]. Climatic conditions, plant species composition, and soil factors have been identified as key parameters shaping the functional and taxonomic composition of soil fungal communities. The main soil factors include soil C, nutrient availability, and pH [e.g., 25,26,37,39,40]. This study describes the soil fungal community structure in a mountain beech forest and investigates how it is determined by plant phenological changes concomitant with season and soil properties. It is hypothesized that (1) soil conditions on calcareous bedrock in conjunction with a highly heterogeneous topology of mountain forests harbor a distinct and diverse fungal community; during the short growing season in mountain beech forests, the compositions of (2) symbiotic and plant pathogenic fungal communities are subject to pronounced seasonal changes, while (3) saprotrophic fungal community composition is primarily determined by soil properties.

Study Site and Soil Sampling
The "Molln" experimental site is located in the Reichraminger Hintergebirge, a mountain range located in the Northern Calcareous Alps of Austria (47°49′ 08″ N, 14°23′ 34″ E). The steeply sloping site (35°) is south exposed, ranging from 1000 to 1100 m a.s.l.; average annual air temperature and precipitation are 7.8°C and 1645 mm. The site is dominated by European beech (Fagus sylvatica L.) and sparsely intermixed by Norway spruce (Picea abies (L.) H. Karst), sycamore (Acer pseudoplatanus L.), European ash (Fraxinus excelsior L.), and silver fir (Abies alba Mill.). Stand age was 146 years in 2015 [41]. Sparse layers of understory tree seedlings, grasses, and herbs were present. The calcareous parent bedrock is limestone and the dominant soil types are Rendzic Leptosol and Chromic Cambisol [42]. Soils are shallow (meañ 30 cm, max.~50 cm) with a high rock content; the predominant humus form is Moder. See Online Resource 1 (Supplementary Materials and Methods, Figure S1, and Table S1) for details on (understory) vegetation, microclimate (i.e., temperature and precipitation from January to August 2015), and soil parameters.
Soil sampling took place during two sampling campaigns in May and August 2015, representing two key phenological growth stages of beech trees in the studied mountainous ecosystem coming along with seasonal changes. In May (i.e., "spring"),~20% of leaves were unfolded (but not yet at full size) on most trees (BBCH12); in August ("late summer"), the foliage was still green and terminal buds were developed (BBCH91) before signs of senescence started occurring in September (i.e., beginning of leaf discoloration and isolated leaf fall) [43]. In late April 2015, 16 plots (~25 m × 25 m) were established at the experimental site; at each plot, four subplots were selected by choosing triangular areas between three mature beech trees (Online Resource 1, Figure S2). In May, three out of four subplots per plot were sampled, resulting in 48 samples. In August, all 64 subplots were sampled. A distance of ≥ 1.5 m was kept to beech trees; no other tree species were in close vicinity to sampling points. For soil sampling, we focused on the mineral topsoil (A-horizon), as the organic horizons LF (litter/fragmented litter) and H (humus) were often rather shallow and highly heterogeneous in extend (data not shown), and the mineral topsoil layer featured the highest carbon accumulation compared to deeper mineral soil layers (Online Resource 1, Table S1). After the organic horizons were removed carefully, approximately 1 l of soil was collected from the upper 10 cm of the mineral soil within each subplot. To mark sampled locations, pits were backfilled with quartz sand. Sampling tools were wiped clean to minimize cross-contamination between (sub-)plots. Soil samples were sieved (2 mm) and homogenized in the field; sieves were thoroughly rinsed in tap water and air-dried between sample processing. For microbial community analysis, 0.5 g fresh soil was weighed immediately into 1.5 ml LifeGuard Soil Preservation Solution (MO BIO, Carlsbad, CA, USA) on site. Soil samples for chemical analyses were stored fresh in sealed plastic bags. All samples were kept at 4°C until further processing.

Soil Parameter and Molecular Fungal Community Analysis
In the laboratory, gravimetric soil moisture content, pH in CaCl 2 , total carbon (C tot ), inorganic C, and total nitrogen (N tot ) contents were determined according to standard procedures as outlined in Online Resource 1, Supplementary Materials and Methods. Organic carbon (C org ) was calculated as the difference between C tot and inorganic C. In early June 2015, detailed physio-chemical analyses were conducted on soil samples derived from additional soil pits (n = 4) to characterize the habitat. The following parameters were determined: fine and coarse root biomass, stone fraction, bulk density, pH, C org , N tot , soil nutrients, and exchangeable ions (Online Resource 1, Table S1).
For DNA isolation from soil samples, 800 μl of the soil suspension in LifeGuard Soil Preservation Solution (see above) was transferred to the wells of a Bead Plate from the PowerSoil-htp 96-Well Soil DNA Isolation Kit (MO BIO, Carlsbad, CA, USA). After centrifugation (4000g, 15 min) and removal of the supernatant, the protocol of the manufacturer was followed with modifications (see Online Resource 1, Supplementary Materials and Methods). The fungal ITS2-region, which has been suggested as preferred barcoding region for fungal community analyses [44], was amplified with primer pair ITS3Mix/ITS4Mix (adapted from [37], for details see Online Resource 1, Supplementary Materials and Methods). Library preparation and Illumina MiSeq sequencing of fungal amplicons were conducted as described [45]. Illumina MiSeq PE250 sequencing was performed at the NGS Unit of the Vienna Biocenter Core Facility GmbH (Vienna, Austria). For details, see Online Resource 1.

Data Evaluation and Statistical Analysis
Initial quality filtering was done with Trimmomatic v. 0.36 [46]. USEARCH (v. 9.0.2132) program suite [47] was used for merging the forward and reverse reads with a minimal overlap of 30 bp with fastq_mergepairs. Sequences < 280 bp were all of non-fungal origin and thus filtered out. Further steps essentially followed [45]. In brief, FASTX toolkit script fastx_barcode_splitter.pl was used to sort out project-specific fungal sequences; USEARCH scripts were used for chimera detection and filtering underrepresented sequences (< 10). VSEARCH [48] was used for clustering and counting sequences per cluster, using a 97% sequence similarity, which is a widely used threshold for the ITS region [e.g., 18] and lies between generally accepted limits for discrimination of species and genera [49]. The results were a sequence file of Operational Taxonomic Units (OTUs) holding one representative sequence per cluster and a table with counts of each OTU per sample. The OTU sequences were aligned using Clustal Omega, and a PhyML tree [50] was calculated. Taxonomic affiliation of OTUs was done with the UTAX script against the UNITE database [51]; manual editing increased phylogenetic accuracy [52]. Given the set sequence similarity threshold, species names provided for OTUs must be considered sensu lato. Non-fungal sequences were removed from further analyses. Ecological guild mapping of OTUs of the total fungal community (TOT) into three functional groups, i.e., saprotrophic (SAP), symbiotic (SYM), and potentially plant pathogenic (PAT) fungi, was based on scientific literature as described previously [53] (Online Resource 2, Table S2). The mapping was primarily conducted at the genus level and refined manually were necessary. Additionally, ectomycorrhizal morphotyping (A. Bittner, unpublished results) allowed mapping selected OTUs affiliated to the Hyaloscyphaceae as SYM. OTUs which could not be categorized were thus labeled "not assigned" (NA). Sequencing and associated data have been deposited at N C B I B i o P r o j e c t P R J N A 5 2 1 6 7 7 , B i o S a m p l e s SAMN12582230-SAMN12582341, and GenBank accession numbers MK626959-MK627467 (see also Online Resource 2, Table S2).
On average 67,875 fungal reads were obtained per sample. For all further calculations, samples were rarefied to 10,022 reads by the function Rarefy from R package "GUniFrac" [54]. No statistically significant differences were found for raw read numbers between seasons (p = 0.71) or plots (p = 0.12); observed OTU richness was independent of raw read numbers before rarefaction (p = 0.32). The observed richness in the number of fungal OTUs, the Berger-Parker index for the most abundant OTU, Simpson's inverse diversity index (1/D), and Shannon's diversity index (only shown for comparison with other studies, where no other diversity index is shown) were calculated for each single sample using R package "vegan" [55] (Online Resource 1, Figure S4). Additionally, to mimic pooling of individual soil samples, as practiced frequently by other studies, rarefied reads from 3 to 4 samples per plot were pooled by season and 1/D of the composite dataset was calculated.
Generalized UniFrac Distances (UF), taking phylogenetic information into account [56], were calculated using R package "vegan" [55]. Two additional indices for β-diversity, Morisita-Horn (MH) and Bray-Curtis (BC) [57], were calculated with the software EstimateS v.9.1 [58] for total fungal communities and functional guilds. Classical multidimensional scaling of the resulting data matrix was performed as principal coordinate analysis (PCoA) to obtain a 2-dimensional representation of dissimilarities [59]. For PCoA, the distances between fungal communities from separate samples were calculated as generalized UniFrac distances using the R package "GUniFrac" [54] with log-scaled data and an alpha value of 0.5 to avoid domination by overabundant species [56]. The "adonis" function [60] from R package "vegan" [55] was used to describe the significances of seasonal differences among the fungal communities.
Environmental distance (ED) for correlation analyses was calculated as the Euclidean distance from differences in space (i.e., horizontal distance of coordinates), soil pH, and soil C org similar to Goldmann et al. [61]. Details are provided in Online Resource 1, Supplementary Materials and Methods.
Random forest modeling was performed using the R package "randomForest" [62] as previously described [45]. Models were calculated with parameters season (i.e., sampling date), soil pH, and soil C org as dependent and OTUs as independent variables. Models with the parameter "season" were calculated separately for TOT and SYM. A total of 2000 trees per model were calculated. During model building, the "out of the bag error" was recorded and the error rates (classification) or mean squared errors (regression) were used to estimate the importance of each OTU. Importance of single variables for correct classification is given as a mean decrease in accuracy (Online Resource 2, Table S2) (see package documentation for details [62]). The most important OTUs for inferring the sampling date from TOT and SYM are listed in Table S3 and S4, respectively (Online Resource 1). OTUs resulting in ≥ 7.5% mean decrease in accuracy for predicting soil pH or C org are listed in Table S5.
Within the manuscript, means and standard deviation (mean ± SD) are given if not otherwise denoted.

Fungal Diversity
In the studied mountain beech forest, 509 fungal OTUs (Online Resource 2, Table S2) were retrieved from 112 soil samples collected in May and August 2015; on average, 231 ±25 OTUs were found per sample (Online Resource 1, Figure S4a). The most abundant OTU (OTU_3) accounted for 5.8% of the total community. Maximum abundance was 19±12% per sample. Certain samples were, however, remarkably uneven with one highly dominant OTU, which could account for up to 75% of the total community ( Figure S4b). The fungal community was highly diverse with an Inverse Simpson's Diversity Index of 20±11 per sample ( Figure S4c) and 32±12 for plots (i.e., pooled data of three to four single samples).
Both in May and August, the total fungal community (TOT) was dominated by Ascomycota followed by Basidiomycota, Mortierellomycota, and Mucoromycota (Fig.  1a). All other fungal phyla (e.g., Zoopagomycota and Chytridiomycota) possessed relative abundances <1% and/or were only occasionally detected. Interestingly, Glomeromycota, which can form arbuscular mycorrhizae with the roots of many plant species [28], were totally absent from the dataset, although some potential host plants were (sparsely) present at the study site as understorey vegetation (see Online Resource 1, Supplementary Materials and Methods). A list of detected fungal OTUs is provided in Online Resource 2, Table S2.

Seasonal Changes
Comparing May (spring) and August (late summer) 2015, no major differences were observed in relative abundances of the phyla Ascomycota, Basidiomycota, Mortierellomycota, or Mucoromycota (Fig. 1a) and of the ecological guilds (Fig.  1b). Only a minor shift in the total fungal composition (TOT) in response to season was observed (3.07% of variance explained by season in Adonis; Fig. 3a); the most important OTUs from TOT for inferring the season are listed in Online Resource 1, Table S3. The majority of OTUs with a high importance for explanation of seasonal variations belongs to the so-called basal fungal lineages, more specifically to the Mortierellomycota and the Chytridiomycota. Relative abundances of single OTUs with significant seasonal changes did, however, not exceed 1%, whereas dominant groups in the fungal community showed less pronounced seasonal shifts.
No major effect of season on the relative abundances of ecological guilds was observed (Fig. 1b). Separate analyses of the ecological guilds indicated that the community compositions of SYM and PAT showed slightly more pronounced seasonal shifts (4.39% and 3.96%, respectively) than those of SAP (3.01%) (Fig. 3b-d); underlying taxa could be identified (Online Resource 1, Table S4). Hygrophorus spp. (i.e., OTU_96 and 125) and ectomycorrhizal Hyaloscyphaceae (i.e., OTU_11 and 27) were among key taxa explaining the seasonal shift in SYM composition-both taxa being less abundant in August compared to May 2015 (Fig. 4a, b). Most OTUs affiliated to the Sebacinales had a greater abundance in August, rendering the seasonal difference between combined sebacinalean OTUs significant (Fig. 4c).

Spatial Heterogeneity and Environmental Influences
Soils at the Molln experimental site showed differences in C org , N tot , soil moisture, and pH (Online Resource 1, Figure S5). Contents of C org and N tot were highly correlated, with a relatively constant C:N ratio of 14.8±1.3. Therefore, the effects of N tot could not be separated from the effects of C org . Gravimetric soil moisture was generally greater in May compared to August and was strongly correlated to C org within each season (Online Resource 1, Figure S5b).
Organic carbon had a minor effect on fungal β-diversity (6.5% of the variation of the Bray-Curtis distance -BC), a  Table S2 for a list of OTUs). SAP, saprotrophic fungi; SYM, symbiotic fungi; PAT, potentially plant pathogenic fungi; NA, fungal guild not assigned  Figure S6). Geographic distance alone could explain 16.4% of the variation of BC. Combining geographic distance and differences in pH and C org into an environmental distance (ED) index explained 26.9% of total soil fungal BC β-diversity (Fig. 5a). ED was significantly related to BC β-diversity per guild (SAP, SYM, and PAT; Fig. 5b). SYM had the greatest β-diversity, PAT the least; in both guilds, BC slightly increased with ED; β-diversity of SAP, on the other hand, increased more strongly with ED. A substantially higher fraction of the variance was explained by ED for the SAP guild (28.5%) than for SYM (8.1%) or PAT (1.8%).
Although no significant influence of soil pH or C org on the relative abundances of the phyla (Ascomycota or Basidiomycota) or the fungal guilds occurred, both parameters  Table S5). A group of fungi in the Leotiomycetes (affiliated to the families Pseudeurotiaceae and Myxotrichaceae) increased in relative abundance from 1.5% at pH 4.7 to 29.7% at pH 7.2 (Fig. 6a); rel. abundances >20% were found at C org contents of 10-20% (Fig. 6b), which occurred at pH >6.5 (Fig.  6c). Similar, the basidiomycetous yeast Saitozyma podzolica (OTU_6) varied in relative abundances between <1.5 at pH <6.5 and C org <6.5% (Fig. 6d, e) to a relative abundance of 16.6% at pH = 7.2 and C org = 12.9% (Fig. 6f). In contrast, greater abundances in soils with a pH of roughly <6.5 and C org <10% were observed for three OTUs from the genus M o r t i e r e l l a , i . e . , M . a f f . e l o n g a t a ( O T U _ 4 7 ) , M. pseudozygospora (OTU_55), and Mortierella sp. (OTU_135; Fig. 6g h, i). The ECM fungus Suillellus luridus tended to occur in samples with a pH >6.8 (data not shown); while holding a maximum abundance of 74% in a single sample, S. luridus was generally a rare species (Online Resources 2, Table S2).

Fungal Community and Diversity
Several studies focused on soil fungal diversity in beech forest stands, often with a specific emphasis on mycorrhizal fungi [e.g 18, 64, 65]. The sites used in these studies were mainly characterized by acidic soils and situated in the lowlands. In contrast, the Molln experimental site is situated on a steep south-facing slope in the Austrian Northern Calcareous Alps, featuring shallow soils with pH values ranging predominantly from 6.5 to 7.3 (total range 4.7-7.3), and high Ca contents (~13 mg g -1 ; Table S1) in the top mineral soil [8].
In agreement with our first hypothesis, we found a highly diverse fungal community in the studied mountain forest. At the experimental site, the median of observed richness was 231 OTUs per sample ( Figure S4a). In contrast, beech forest stands in Northern, Central, and Southern Germany were reported to hold an observed soil fungal richness of <180 OTUs per sample and featured inversed Simpson's indices <20 (for pooled samples) [18]. Pooling of sequence data from three (May 2015) or four (August 2015) soil samples at our mountainous beech forest stand resulted in inverse Simpson's indices of 26 and 34 in May and August, respectively. This difference might be partially attributed to the use of different barcoding regions used for fungal community profiling, i.e., the whole ITS region by Wubet et al. [18] vs. the ITS2 region in our study. However, a greater fungal diversity at near neutral soil pH is in agreement with recent findings, where fungal diversity responded in a unimodal relationship to soil pH with a peak around pH 6 [66]. Similarly, a strong positive influence of soil Ca content on the richness of fungi was reported earlier [37]. Both point toward a more diverse soil fungal community in beech stands on calcareous bedrock compared to beech forests with more acidic soil conditions. Other potential effects of soil parameters are discussed below.
The horizontal spatial heterogeneity of the fungal community within a maximum geographic distance of~700 m was very high; three different indices of β-diversity were applied to examine the dissimilarities (Online Resource 1, Fig. 5 Environmental influences on fungal β-diversity at the Molln site. a Indices of β-diversity, Morisita-Horn (MH, blue), Bray-Curtis (BC, red), and UniFrac (UF, grey), were calculated from the fungal community data and plotted against the environmental distance. b BC β-diversity of symbiotic (SYM, green), saprotrophic (SAP, orange), and potentially plant pathogenic (PAT, purple) fungal guilds. Regression coefficients (r 2 ) and significance levels (p < 0.001) are shown in corresponding colors Figure S6a). The Bray-Curtis index (BC) possessed intermediate, and the UniFrac distance (UF) low dissimilarities, pointing towards a highly conserved phylogenetic composition of the fungal community. In contrast, the high Morisita-Horn index (MH) values indicate pronounced changes in abundance of dominant OTUs at short distances (<10 m) potentially related to the stands' heterogeneous soil topography.
The fungal community at the Molln experimental site was remarkably distinct compared to those reported earlier. Specifically, Ascomycota was more frequent than Basidiomycota at the studied site (Fig. 1a). Most previous studies, including those conducted in lowland beech forests, reported higher frequencies of Basidiomycota or at least a balanced ratio of both phyla [18,37,61,64,[67][68][69]. While the prevalence of Ascomycota was thus considered untypical for soil of temperate forest in general and beech forest in specific, our study indicates that this assumption might not hold true for mountainous beech stands on calcareous bedrock.
As many basidiomycetes are ECM fungi, the relative abundance of ECM fungi was also low, accounting for <25% of the total soil fungal community in the vast majority of samples (Fig. 1b)-lower than reported for other European beech forests. Temperate forest soils are generally dominated by ECM fungi [e.g., 70] and published data on relative ECM abundances in (lowland) beech forest soils range from 32 [67] to 90% [68]. The observed differences might be related to distinct environmental factors of the studied beech stand (as discussed below). A combination of soil factors, which favor saprotrophic fungi over symbiotic fungi, has probably contributed to low ECM abundance in soil. Beech root colonization by ECM fungi has previously been described to be higher [71] or lower [72] on calcareous soils with near-neutral pH compared to siliceous soils with acidic pH. The use of different primer target regions for fungal community profiling in different studies could have partially contributed to the observed differences. Op De Beeck et al., however, reported that using e.g. either ITS1 vs. ITS2 regions to determine the relative abundances of Ascomycota and Basidiomycota had a relatively minor effect [73]. Among ECM fungi, the genus Russula is one of the most broadly distributed worldwide [74]. As Russulaceae were common among the dominant members of the ECM community of beech stands [18,61,64,[67][68][69], the complete absence of Russulaceae from the ECM community of the studied mature mountainous beech stand was highly surprising. A strong bias by the amplification protocol cannot be expected, as no mismatches of the primers were found for representative sequences of Russulaceae in the NCBI GenBank. Furthermore, Russulaceae were regularly found in soil samples from other temperate deciduous forests with the original primers [37] (for visualization see Fig. 77 in [24]) and in a parallel study with the same amplification and sequencing approach as applied herein (M. Gorfer, unpublished results). As a decline in Russula sp. abundance with increasing pH and decreasing C:N ratio was previously reported by others [18,66,75,76], we speculated that the absence of Russula spp. at the site may be (partially) related to the near-neutral pH and a low C:N ratio [77], as discussed below in detail. Negative priority effects, i.e., pioneer fungal species reducing the subsequent colonization of a host by additional species, could have also contributed to the observed patterns-as previously shown for both ECM and AM colonization [78] and references within]. In addition, it has been suggested that long-distance spore dispersal is generally limited for Russula spp. [74]-potentially hampering spore dispersal particularly in a landscape dominated by mountain-valley systems and complex wind systems. Apart from the absence of Russulaceae, the ECM fungal community was similar to other beech-dominated forest stands, i.e., exhibiting relatively high abundances of Inocybe spp. and Sebacinales (Figure 2b, see [68] for comparison). A preference for high pH and low C:N ratios was previously described for Inocybe spp. and Sebacina spp. [18,79] but also for Suillellus luridus [80] which however, was a rare species at the studied site.
In contrast to symbiotic fungi, high abundances of potentially plant pathogenic fungi were observed at the study site (Fig. 1b). On a global scale, soil fungal communities in temperate deciduous forest harbor <5% phytopathogenic fungi [37], while relative abundances of >15% were found at our study site. The potentially plant pathogenic community consisted mostly of taxa from the nectriaceous genera Dactylonectria, Ilyonectria, and Neonectria-causing root rot and other diseases in a wide variety of plants [81]. Protection of roots from pathogens by ECM fungi is well described; thus, an unfavorable ratio of symbiotic to plant pathogenic soil fungi could make mountain beech forests especially vulnerable to biotic stress. We can only speculate if the environmental conditions in harsher mountain forest habitats related to temperature extremes and restricted resource availability (i.e., shorter growing season and poor soil conditions) deem roots more susceptible to pathogenic fungi [82]. The topic of linkages between tolerance to both abiotic stress and potentially plant pathogenic fungi warrants further investigation, particularly addressing the rhizosphere of forest ecosystems on marginal sites.

Seasonal Variation of Fungal Community Composition
A pronounced seasonality was reported recently for the fungal soil communities of temperate deciduous forests in China, with the diversity of the total fungal community being highest in July when trees possessed the most vigorous growthpotentially providing more C assimilate for fungal growth [32]. In our study, however, only subtle seasonal changes in the soil fungal community composition of the mountainous beech forest were found between spring (May) and late summer (August) 2015 (Figs. 1, 2, and 3). Species counts and βdiversity indices did not indicate major shifts in overall fungal community diversity (Online Resource 1, Figure S4). The abundance of fungal phyla (Fig. 1a) remained largely similar although some changes in specific taxa occurred. Those changes were generally more distinct not only in rare taxa, especially from the Mortierellomycota and the Chytridiomycota (both SAP; Table S3), but also in SYM-classified members of Hyaloscyphaceae and Sebacinales (Fig. 4c).
Based on previous findings, we hypothesized that taxa of the SYM guild, including Sebacinales with ectomycorrhizalike root interactions, would be strongly affected by seasonal changes between spring and late summer [10,32]. Recently fixed carbon is an important driver of soil biological processes [83], and C allocation belowground differs strongly with plant phenological and physiological states; thus, seasonal changes in the ECM community were previously linked to seasonal patterns in photosynthetic products and nutrient availability [84,85]. In contrast, the richness of SAP fungi was earlier reported to be directly sensitive to soil moisture availability [37]. In partial agreement with our hypotheses, the SYM guild showed indeed a slightly more pronounced community response to season compared to the SAP guild, although overall seasonal differences in both guilds remained small (Fig. 3b-d). However, SYM-associated taxa with significantly greater abundances either in spring (Hygrophorus spp., Hyaloscyphaceae) or summer (Sebacinales) were identified (Fig. 4), indicating specific effects of different symbiotic taxa to seasonal patterns [84] but not generally greater abundances of SYM-associated taxa in summer (Fig. 1b). As storage products can play an important role in seasonal carbon allocation (e.g., spring vs. autumn) [86] and up to 12% of C in ECM fungi may originate from soil [87], we speculate that ECM taxa with a greater abundance in spring may rely less on recently assimilated C and are replaced by other symbiotic taxa such as Sebacinales later in the growing season. In addition, a different availability of soil nutrients in spring vs. late summer, driven by an increased uptake by roots and decreasing soil moisture availabilities until August, may underlie the pattern. In contrast to our findings, Voříšková et al. [31] and Štursová et al. [10] reported significantly greater relative abundances of several ECM fungal species in summer than in spring. We can only speculate that the less productive, deciduous mountainous stand at the Molln site provides less pronounced differences in seasonal C and nutrient availability belowground compared to previously studied (coniferous) forests with longer growth periods and greater resource availabilities.
Here, seasonal changes in relative abundance were generally more distinct in rare taxa, especially in taxa from the SAPassociated group Mortierellomycota and Chytridiomycota ( Fig. 2a; Online Resource 1, Table S3). Three rare SAPclassified OTUs (relative abundance <1%) affiliated to the genus Mortierella showed greater relative abundances in (late) summer than in spring (Fig. 4b). Mortierella spp. are often classified as sugar fungi that as r-selected species can rapidly respond to transient nutrient peaks [88]. It is thus speculated that the observed higher relative abundance of selected Mortierella spp. in August might be caused by an increase in available nutrients by ongoing decomposition at sufficiently high soil moisture levels ( Figure S5b; see [53,88] for Mortierella spp., and [89] for Chytridiomycota's soil moisture dependencies). Our results are in accordance with recent findings in a mountainous Picea abies forest, where the abundance of SAP-associated fungi in soil generally varied little between summer and spring, but Mortierella's abundance increased in summer [10]. However, extreme weather events like extended drought spells, which were absent in the year of our study ( Figure S1), are expected to induce more pronounced seasonal effects on soil fungal communities in general [90] and saprotrophs in particular [37]. Soil moisture at the Molln experimental site was lower in August than in May 2015, but was well above levels where drought stress would occur ( Figure S5b). Seasonal differences in C availability to the soil fungal community thus warrant further investigation, particularly considering more parallel measurements of parameters underlying strong seasonal variation.

Soil pH and Organic C Effects on Fungal Community
Our third hypothesis stated that the saprotrophic fungal community composition is primarily determined by soil properties. Soil chemistry in general, and soil pH in particular, has often been identified as a major parameter shaping soil fungal communities [37,40,66,91], but see [92]. In accordance, the fungal β-diversity increased with environmental distance (ED) between sampling points, i.e., with combined changes in geographic distance, pH, and C org contents (Fig. 5a). A substantial fraction of the variation remained, however, unexplained by the variables pH, C org , and geographic distance. Factors that were highly correlated to C org , i.e., N tot and soil water content (Online Resource 1, Fig.  S5a,b), could not be analyzed separately; soil Ca and P concentrations-of potentially predictive value [37]were unfortunately not measured at the sample level. As hypothesized, ED explained a substantially higher fraction of the variability in BC diversity for SAP than for PAT or SAM guilds (Fig. 5b). Differences in soil pH and geographical distance alone could partially explain the increase in fungal β-diversity across the site, while C org was a weaker predictor of β-diversity (Online Resource 1, Figure S6).
The gradients in pH and C org allowed an identification of taxa, predominantly of the SAP guild and NA taxa, that showed strong responses to soil parameters (Table S5, Fig.  6). An abundant group of OTUs classified as Leotiomycetes i.s were preferentially found at mineral soil locations with a pH >6.5 and C org =10-20% (Fig. 6a-c). Interestingly, the genus Oidiodendron (Myxotrichaceae, Leotiomycetes i.s.) was shown by others to have a preference for acidic soils [93] but occurred at low abundance in more acidic parts of the studied site. We can only speculate that the abundance of Oidiodendron species is co-determined by the presence of intermediate to high C org contents as lower pH values were correlated with lower C org values at our study site ( Fig. 6c; Online Resource 1, Figure S5c). Selected OTUs from Mortierella showed a preference for soil pH <6.5 and an associated C org <10% (Fig. 6 g, h). Mortierellaceae were previously shown to rapidly colonize organic materials with preferred utilization of non-cellulosic C sources [88]. The preference for lower C org contents was thus surprising, but might be driven by pH preferences. Soil pH and C org had pronounced effects on rather abundant taxa, which were mostly classified as saprotrophic. Response curves to changes in pH and C org were linear, exponential or polynomic (Fig. 6), which reflects different subsections of unimodal response curves often described for environmental factors [e.g., 94]. A strong influence of soil pH, which directly affects availability of major nutrients and toxic elements, on soil microbial community composition has been repeatedly reported [e.g., 37,40]. Saprotrophic fungi, which largely depend on breakdown of complex organic matter through the activity of extracellular enzymes, are additionally affected by the pH optima of these enzymes.
The initially stated hypothesis that soil parameters primarily shape the SAP fungal community was thus supported by our results. A substantial fraction of the observed β-diversity remains, however, unexplained by the factors soil pH, C org , and geographic distance. It is thus assumed that additional factors including understorey vegetation, topographic heterogeneity, and further soil parameters such as Ca contents contribute to the observed high diversity.

Conclusion
The studied mountain beech forest stand in the Northern Calcareous Alps of Austria showed a distinct and highly diverse soil fungal community with comparatively low relative abundances of basidiomycetes and of ECM fungi, while potentially plant pathogenic fungi were more prevalent than in previous studies (in lowland beech forests). Seasonal differences between May and August were minor and influenced mainly rare taxa. Soil properties like pH and C org affected fungal community composition, affecting particularly the distribution pattern of dominant taxa and shaping the saprotrophic community.
Acknowledgements The research site was kindly provided by Austrian Federal Forests (ÖBf AG). The Melitta Group | Cofresco Frischhalteprodukte GmbH & Co. KG is acknowledged for kindly providing soil sampling bags free of charge. Alexander Urban is acknowledged for the provision of fruiting bodies of Russula spp. and Lactarius spp. for testing of primers. Illumina MiSeq PE250 sequencing was performed at the NGS Unit of the Vienna Biocenter Core Facility GmbH (Vienna, Austria; www.viennabiocenter.org/facilities/next-generationsequencing). Three anonymous reviewers provided constructive criticism that improved a previous version of the MS.
Code availability Software applications are fully disclosed in the "Materials and Methods" section. Data availability Sequencing and associated data have been deposited at NCBI BioProject PRJNA521677, BioSamples SAMN12582230-SAMN12582341, and GenBank accession numbers MK626959-MK627467. Environmental data will be made available upon request.

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