Variability of Functional Groups of Rhizosphere Fungi of Norway Spruce (Picea abies (L.) H.Karst.) in the Boreal Range: The Wigry National Park, Poland

Rhizosphere microbial communities can influence plant growth and development. Natural regeneration processes take place in the tree stands of protected areas, which makes it possible to observe the natural changes taking place in the rhizosphere along with the development of the plants. This study aimed to determine the diversity (taxonomic and functional) of the rhizosphere fungal communities of Norway spruce growing in one of four developmental stages. Our research was based on the ITS region using Illumina system sequencing. Saprotrophs dominated in the studied rhizospheres, but their percentage share decreased with the age of the development group (for 51.91 from 43.13%). However, in the case of mycorrhizal fungi, an opposite trend was observed (16.96–26.75%). The most numerous genera were: saprotrophic Aspergillus (2.54–3.83%), Penicillium (6.47–12.86%), Pyrenochaeta (1.39–11.78%), pathogenic Curvularia (0.53–4.39%), and mycorrhizal Cortinarius (1.80–5.46%), Pseudotomentella (2.94–5.64%) and Tomentella (4.54–15.94%). The species composition of rhizosphere fungal communities was favorable for the regeneration of natural spruce and the development of multi-generational Norway spruce stands. The ratio of the abundance of saprotrophic and mycorrhizal fungi to the abundance of pathogens was high and promising for the durability of the large proportion of spruce in the Wigry National Park and for forest ecosystems in general.


Introduction
The rhizosphere is where interactions between microorganisms and plants take place.Roots in the rhizosphere have an intimate association with the microbes that inhabit the soil [1].This rhizosphere microbiome includes microorganisms associated with the surface of the roots and the soil layer adjacent to them [2].For forest ecosystems, fungi are of particular importance, as they have a comprehensive effect on the operation of both the tree and the entire woodland community [3].Changes in taxonomic and functional diversity, and the abundance of fungal communities, depend on environmental conditions and the health, vitality, and developmental stages of host plants [4][5][6][7][8].Many plant-associated microorganisms affect seed germination, seedling vigour, growth and development, nutrition, disease susceptibility, and plant productivity [9].The community consists of species that perform several functions in the soil, ranging from saprotrophs through endotrophs or mycorrhizal fungi to pathogens.Each functional group affects, directly and indirectly, the dynamics of the forest ecosystems.
Ectomycorrhizal fungi combine with plant roots in symbiotic systems that provide plants with water and mineral salts and oxidize organic matter to acquire nitrogen, dissolve phosphate, and reduce pathogens [10].The vast majority of mycorrhizal fungi are generalists.Once a fungus has colonized a host plant, its mycelium can grow and extend for a long distance in the soil and reach and colonize the roots of neighboring plants of the same or different species [11].Mycorrhizal fungi connect individual plants to form common mycorrhizal networks (CMN) [12][13][14].This can modulate the survival and behaviour of connected plants, their competitiveness and cooperation (symbiosis), and consequently affect plant diversity [15,16].For example, improvement in seedling survival has been reported due to plant-fungal interactions [17,18].Saprotrophic fungi break down soil organic matter to release carbon, thereby contributing to the nutrient cycle [19].Elements used to build new cells are released as a result of their action.Endophytes colonize plants' tissues without visible disease symptoms.They inhabit plants by vertical transfer (from seeds) or penetrate from the external environments [3,6].Suppression of diseases by beneficial microorganisms in the rhizosphere can occur using several mechanisms of action, such as antagonism associated with the production of antibiotics, antifungal metabolites, competition for space and nutrients with phytopathogens and other microorganisms harmful to the rhizosphere, and induction of resistance in plants.The pathogenic fungi cause plant diseases, impairing growth and reproduction and they may even cause the death of host plants.They require conditions conducive to infection, such as drought or wounds [20].
The rhizosphere microbiome changes and adapts to current environmental conditions, i.e., biotic and abiotic stress [21][22][23].In the early stages of plant growth, when there is a poor root density, the fungal community is dominated by species with relatively long hyphae and fast growth (e.g., the genera Cenococcum, Cortinarius, Inocybe, Laccaria, Scleroderma, Sebacina, Thelephora i Tomentella, and the order Pezizales) [24].As the plants in the community develop, the ectomycorrhizal species, in particular, become more diverse [24,25].
Following disturbances arising from storms, fires, or pest outbreaks, forest stands in protected areas such as national parks typically comprise some surviving large living trees, significant deadwood, gaps in the canopy, and coexistence of different developmental stages, and are then succeeded by natural regeneration dynamics.As a result, the newly established stands are diverse in terms of age and biometric characteristics [26].Understanding the mechanisms and conditions that influence natural regeneration is a critical topic in forest ecology and management, especially in the context of the ongoing climate crisis [27].In recent years, the role of microorganisms potentially involved in natural regeneration has attracted much attention.The development of molecular methods, including next-generation sequencing, has enabled the full spectrum of microorganisms from tested samples to be identified and quantified.Consequently, this progress in technology has changed the possibilities for research in this field, from just observing a single taxon to monitoring the whole community and facilitating the identification of the interactions between functional and phylogenetic groups [28].
Norway spruce (Picea abies L. (H.Karst.)) is one of the dominant species of coniferous trees in the northern temperate and boreal forests of Europe.It grows in various climatic and edaphic conditions, but it is best adapted to cooler climates.It most often occurs on acidic and nutrient-rich soils with good moisture availability [29,30].There are two distributional domains of this species in Poland-eastern (boreal) and southern (lowland and mountainous).We are increasingly observing the weakening of the health of spruce stands throughout Europe due to a variety of abiotic and biotic factors, which are reflected in the mass die-off of populations [31].Therefore, it is crucially important to identify the factors that could promote the natural regeneration and development, and enhance the resilience of this species in its natural environment.
The current study aimed to determine the relationship between the developmental stages of spruce and the diversity of rhizosphere fungal communities, with a focus on the function performed by individual taxa.The following hypotheses were formulated: (i) the species composition and abundance of particular functional groups of rhizosphere fungal communities will differ between the development stages of Norway spruce, from the seedling to the mature old trees in the forest; (ii) along with the change in the developmental stages, there will be a concurrent change in the share of saprotrophs and mycorrhizal fungi, and (iii) a high abundance of saprotrophs and mycorrhizal fungi will be associated with low abundance of pathogenic fungi.

Results
Following normalization, a library of 412,560 Operational Taxonomic Units (OTUs) was obtained, representing 599 taxa.The identified taxa belonged to 6 phyla: Ascomycota, Basidiomycota, Chytridiomycota, Entomophthoromycota, Mucoromycota, and Zoopagomycota.In the four stages of development, the range of individual taxa ranged from 32.84 to 55.88% for Ascomycota and Basidiomycota from 31.38 to 67.03%.Figure 1 shows a clear trend reflecting a decreasing share of Ascomycota taxa and an increasing share of those from Basidiomycota with consecutive stages in tree development.Statistical tests revealed a significant difference between these two phyla (Table 1), but no substantial differences were found between tree stands.
progressed from seedlings to mature trees.In the case of mycorrhizal fungi, the average of OTUs in the rhizosphere of seedlings was 16.97%, saplings up to 50 cm (small saplings) 20.35% saplings above 50 cm (large saplings) 28.45%, and mature trees (mature trees) 26.76%.For pathogens, the average relative abundance in the rhizosphere in seedling ranged from 14.77% OTUs, in small saplings from 16.70% in large saplings from 6.81%, and in mature trees from 11.57% (Figure 4).Similarly, the developmental stage of trees was an important factor in explaining the variability of the abundance of functional groups: ectomycorrhizas and pathogens (Table 2).Figure 3a shows a general trend reflecting an increasing share by the former, and Figure 3b a falling trend in the latter, with consecutive three developmental stages.The abundance of saprotrophs, which dominate quantitatively, showed a decreasing trend from the seedling stage up to the mature tree stage (Figure 3c), which was just marginally the wrong side of our cut-off for significance (p ~ 0.06 level).
As a result of the alpha diversity analysis, communities differing in terms of species composition and abundance of individual taxa were identified.However, statistical significance was only found in the post-hoc test of the ChaoI index.The highest mean value of Chao was for small saplings (Table 3).Non-significant differences were found within fixed factors for the Shannon and Simpson indices.1.  1.

Ectomycorrhizas
An accurate identification up to the species or genus level was conducted to enablae the functional characteristics of taxa to be grouped, as most functional elements of fungi are conserved at the generic level and sometimes at higher taxonomic levels [32][33][34].Analyzing the communities with respect to functional groups, saprotrophs were the most numerous group.Across four development stages, the average abundance of individual functional groups varied and in seedling, we noted 51.90% OTUs, in small saplings 44.02%, in large saplings 40.29%, and in mature trees 43.12%A decreasing trend in the proportion of saprotrophic fungi was observed in small saplings and large saplings as development progressed from seedlings to mature trees.In the case of mycorrhizal fungi, the average of OTUs in the rhizosphere of seedlings was 16.97%, saplings up to 50 cm (small saplings) 20.35% saplings above 50 cm (large saplings) 28.45%, and mature trees (mature trees) 26.76%.For pathogens, the average relative abundance in the rhizosphere in seedling ranged from 14.77% OTUs, in small saplings from 16.70% in large saplings from 6.81%, and in mature trees from 11.57% (Figure 4).Similarly, the developmental stage of trees was an important factor in explaining the variability of the abundance of functional groups: ectomycorrhizas and pathogens (Table 2).Figure 3a shows a general trend reflecting an increasing share by the former, and Figure 3b a falling trend in the latter, with consecutive three developmental stages.The abundance of saprotrophs, which dominate quantitatively, showed a decreasing trend from the seedling stage up to the mature tree stage (Figure 3c), which was just marginally the wrong side of our cut-off for significance (p ~0.06 level).
As a result of the alpha diversity analysis, communities differing in terms of species composition and abundance of individual taxa were identified.However, statistical significance was only found in the post-hoc test of the ChaoI index.The highest mean value of Chao was for small saplings (Table 3).Non-significant differences were found within fixed factors for the Shannon and Simpson indices.

Discussion
In this paper, we have identified the rhizosphere fungal taxa of Norway spruce from the Wigry National Park, Northeastern Poland, and determined their functional groups, on four developmental stages of these conifers: on current year seedlings, saplings up to 50 cm (small saplings) and over 50 cm in height (large saplings) and mature trees.As expected, the greatest alpha diversity was in the rhizosphere of small spruces.The highest abundance of Ascomycota was in the rhizospheres of current-year seedlings and small saplings up to 50 cm in height, while the highest abundance of Basidiomycota was in the rhizosphere of saplings that were over 50 cm in height and in mature trees.The Ascomycota phylum dominates most ecosystems.This is related to their ecological requirements and functional features, which translate into the abundance of spores and how they spread [35].The ability to disperse is the main factor determining the composition of early succession fungal communities.In the case of the Ascomycota, saprotrophic taxa predominate [36], and therefore, as expected the share of saprotrophic species in our study was higher in the rhizospheres of younger trees compared with older trees.The more significant share of saprotrophic fungi in the rhizosphere of the juvenile development stages of spruce stands can be explained to some extent by the appearance of natural regeneration in the gaps in the tree canopy formed after the death of trees, next to decaying pieces of wood and even on rotten and only partially decomposed spruce trunks [37], and a relatively small root tissue volume in the sample.The abundance of wood lying in gaps in the forest canopy is conducive to the mass occurrence of saprotrophs, and at the same time, these spaces create favorable conditions for the emergence and development of a young generation of trees.Among the saprotrophic fungi, genera such as Aspergillus, Talaromyces, Phlebia, Phlebiopsis, and Penicillium prevailed in the Wigry National Park.These taxa are considered antagonists to the genera Armillaria and Heterobasidion [38][39][40].The mycelium and its metabolites reduce the extent of rot caused by Armillaria spp.and the extent of base necrosis caused by Heterobasidion in Scots pine and other tree species.These taxa compete with pathogens for living space and food.They have the ability to dissolve phosphates and produce antifungal agents, e.g., citrinin, griseofulvin, fellutanin A, and archidic acid [41].
The presence of the genus Mycena in the analyzed community indicates that decaying wood (like that of spruce roots) lies near sampling sites.Mycena spp. is the cause of white rot decay and has been previously identified as a significant component of fungal communities in later stages of root [35] or wood decay [36].Our results did not confirm the abundant occurrence of this taxon in the late developmental stages of trees, which may be due to the different (rhizosphere not soil) material collected for the current research project.The rhizosphere consists of small roots and the soil adjacent to it, and not of decomposing plant tissues.
It is known that saprotrophs and mycorrhizal fungi are essential in creating suitable conditions for future stands of trees to grow and develop from the juvenile stage [37][38][39].As these grow, the diversity and abundance of mycorrhizal species increase.Among the mycorrhizal fungi, the most numerous were: Cortinarius, Amanita, Inocybe, Pseudotomentella, and Tomentella.Species of these genera are common mycorrhizal species found in Pinaceae worldwide [40,41].Taxa from the genera Cenococcum, Cortinarius, Inocybe, and Tomentella, as well as the order Pezizales in the above work, were identified from the early stages of development [40], which confirms their importance as the first colonizers of roots.Some of these mycorrhizal fungi have features of saprotrophs and usually show a poor ability to decompose carbon compounds [10].Fungal communities make it possible for trees to adapt to the current environment [32,42], an ability associated especially with the genus Amanita, among others.This genus is a good indicator of the presence of mature trees.After the removal of the host plant, it disappears very quickly from the site [35].In our research, this genus was most abundant in the rhizosphere of mature trees.Representatives of this genus were scarce in the rhizosphere of the other developmental stages.Thelephora is an ectomycorrhizal basidiomycetes genus, of which all constituent species.form ectomycorrhizal relationships with diverse plants [39].As mycorrhizaformers, Thelephora plays a very important role in pioneering microhabitats in coniferous forests and is also involved in the decomposition of dead wood [43].Thelephora wakefieldiae is one of the most important fungi establishing symbiotic relationships with Norway spruce in the Wigierski National Park.It dominated in each of the samples, although it was most common in large saplings.
We observed the prevailing dominance of saprotrophs and mycorrhizal fungi in our study sites.Accordingly, the high abundance of mycorrhizal species constitutes a natural biological defense against pathogens, especially those causing diseases of root systems, including fungi belonging to the genera Fusarium, Verticillium, Typhula, Gaeumannomyces, Calonectria, Nectria, Rhizoctonia and fungus-like organisms (Oomycota): Phytium and Phytophthora genera, as well as plant-parasitic nematodes.These species are responsible for the most significant threat to the early stages of development of saplings; among other pathogenic consequences, they cause the blight of seedlings.In mature tree stands, mycorrhizal symbiosis is common and has a positive effect on the growth and development of trees.
We identified also taxa belonging to the genus Alternaria, which includes species of opportunistic plant pathogens contributing to losses in fruit and seed storage, found worldwide [44].Alternaria spp.cause 20-80% agricultural losses in field, horticultural, plantation, and forest plants, as well as post-harvest storage.The Alternaria species cause needle blight in Pinus bungeana in China [45].Infected seeds suffer from blight and death of the necrotic cotyledons.Infection occurs through the lower leaves of seedlings, which are infected by spores from the soil.Initially, there are slight changes on the leaves, which develop to a size of 5-15 mm in diameter with light and dark concentric rings.The spores are carried also by rain and wind and infect other tissues.Wet and dry weather favors the life cycle (germination and spore release) of Alternaria [46].
The genus Diplodia was identified in the analyzed community.Diplodia sapinea, a species common throughout Europe, although in our research, it was rarely cause dieback of shoots of coniferous trees [47][48][49].
Among the pathogens, we identified taxa responsible for diseases of the root system belonging to the genus Heterobasidion, including H. parviporum and H. annosum, and a small contribution of Armillaria spp.Heterobasidion parviporum mainly infects spruce [50], while H. annosum has a broader host spectrum and infects both Scots pine and Norway spruce [51].The highest abundance was in seedlings, although we did not observe any pathological symptoms of Heterobasidion infection.This pathogen quickly kills young trees, while in older trees, it causes long-term disease, consisting of the decay of wood and roots of the lower parts of the trunk.In older spruces decomposition only occurs in the heartwood, and as a result, their vital functions remain normal [52].These fungi are not present in smaller roots.
The genus Armillaria causes dieback of trees older than 2 years.Colonization of the root system of large and old trees is slow and sometimes takes decades.In younger trees, the dying process is rapid and lasts from several months to a year after infection.The presence of the rhizomorph Armillaria spp. was found in the soil in the stand in Mikołajewo, while in the analyzed community the abundance of the DNA of this pathogen was low.However, the abundance of pathogens and mycorrhizal species did not differ significantly between the plots.Statistical differences occurred in the case of saprotrophs for this area.In the case of Armillaria spp., the risk of infection increases as a result of drought.Tree roots weakened by drought are more susceptible to infection with rhizomorphs of the pathogen [53].However, it should be remembered that Armillaria has so far been identified in roots larger than 2-3 cm, while in this study we investigated only the fine root.
Pyrenochaeta inhabit soil and plant debris worldwide and are well-known as pathogens of various plants.Several species of Pyrenochaeta fungi infect plants, particularly plant roots [54].In the environment, numerous Pyrenochaeta species are found as saprotrophs in soil, plant debris, and wood [55].Some species have been identified as tree endophytes.Some of the Pyrenochaeta species cause serious plant diseases in agriculture and forestry [56].Herpotrichia juniperi and Neopeckia coulteri (Pyrenochaeta anamorphs) cause shoot and needle diseases of conifers [57].For example, Pyrenochaeta parasitica/Nematostoma parasiticum occurs on Picea and Tsuga (the Herpotrichia needle browning) [58].The species has been recorded in Austria, Switzerland, Denmark, Germany, Norway, Great Britain, and Poland [57,59], and sporadically in North America [60].In our results, we identified this taxon only at the genus level.Fungi belonging to this genus may reduce the young spruce generations in the Wigry National Park.
Our study provides new insights into the variation of fungal communities among the developmental stages of spruce stands, however, we are aware of the limitations of our findings.Even within the same genus, some species may be symbiotic and others pathogenic.There are controversies in concluding that a species is pathogenic based on its ITS sequence alone.Even for species that are clearly identified as plant pathogens, there are reports of complex symbiotic relationships in which they do not express pathogenicity but rather enhance the plant's immune system [61].
The fungal community consisting of a large abundance of saprotrophs and mycorrhizal species significantly reduces pathogen threats to the Norway spruce.In addition, plants have access to the elements necessary for life and for building new cells Therefore, the fungal community identified in Wigry National Park potentially benefits the natural regeneration, growth, and development of Norway spruce.

Sampling Site Description
Three uneven-aged stands with a differentiated vertical and species structure, and a significant share of the Norway spruce were selected for the study (Table 4).The Serraulo-Pinetum (W.Mat.1981) J.Mat. 1988 (subboreal mixed forest) community dominates all stands.They were chosen to represent three parts of the Park-north, south-western, and south-eastern, where forests dominate.Since 1989 these stands have been excluded from silvicultural management.They are protected from interference and are reserved for scientific study, including observation of natural processes as a main objective (Poland, Regulation of the Council of Ministers of 27 June 1988-Journal of Laws No. 25 of 21 July 1988, item 173).

Sampling
The collection of research material was carried out in June 2022.For the analysis of the communities, twenty trees were selected randomly in each stand as follows.Five vigorous individuals of Norway spruce, with no disease symptoms, from four forest strata were selected in each stand: current-year seedlings with visible cotyledons from the undergrowth, saplings up to 50 cm in height, and over 50 cm in height from the small and large saplings strata respectively, and mature trees with visible generative organs from the overstorey.Their roots were tracked, and surrounding soil at a depth of up to 15 cm was collected.Finally, sixty samples were collected (three stands × four forest strata × five samples = 60).Samples were placed in a plastic container and immediately taken to the laboratory.Comparably sized subsamples of fine root branches (2 cm long with soil) were placed in 2 mL screw-cap tubes.Samples were stored at −20 • C until molecular analyses.

DNA Extraction and PCR Amplification
The collected material was ground in a mortar chilled to −70 • C. Then DNA was isolated using the DNeasy PowerSoil Kit (QIAGEN, Hilden, Germany).The DNA was then purified using the Anti-Inhibitor Kit (A&A Biotechnology, Gdynia, Poland).Total DNA was sent for sequencing and preparation of OTU (operational taxonomic unit) libraries at Genomed SA (Warsaw, Poland).Fungal community analysis was performed based on the Internal transcribed spacer 1 (ITS1) region using specific primers ITS1FI2 5 -GAA CCW GCG GAR GGA TCA-3 [62,63] and 5.8S 5 -CGC TGC GTT CTT CAT CG-3 [64].PCR was performed using Q5 Hot Start High-Fidelity 2X Master Mix (New England Biolabs, Ipswich, MA, USA).Reaction conditions were as recommended by the manufacturer.Sequencing was performed on a MiSeq sequencer in paired-end (PE) technology.Negative samples (without DNA) were also sequenced to remove artifacts.

Sequence Processing and Statistical Analysis
The obtained sequences were subjected to bioinformatics analysis to generate the OTU library.The first step was classifying reads to species level using the QIIME2 software package (ver.2017.6.0).The study consisted in removing the adapter sequences in the cutadapt software (v4.4) and then in the qualitative analysis of the reads.Low-quality sequences (quality < 25, length up to 30 bp) were removed.In the next group, paired sequences were combined (fastq-join algorithm) [65], chimeric sequences were removed (usearch algorithm), clustering was performed based on the selected reference sequence database (uclust algorithm) and taxonomy was assigned based on the NCBI (National Center for Biotechnology Information, https://www.ncbi.nlm.nih.gov/(blast algorithm ver.BLAST+ 2.14.0)) (accessed on 27 May 2023) [66].For identification, the percentage of sequence similarity with the reference sequence was assumed to be 98-100%, with a minimum coverage of 90% [67,68].Taxa were identified down to the lowest possible taxonomic level.OTU sequences not belonging to Fungi or Oomycota were removed from further analysis.Further, OTUs were assigned to ecophysiological categories using FungalTraits [69], with subsequent manual corrections.Finally, the taxa were divided into pathogens, saprotrophs, endotrophs, mycorrhizal fungi, and lichens.The applicable nomenclature was adopted from MycoBank (https://www.mycobank.org/,accessed on 13 May 2023).Rarefaction curves were determined for the obtained OTU library.Rarification samples were made to normalise the samples to one size.As a result of normalization, two negative samples were removed.The ChaoI species diversity index and Shannon diversity index and Simpson index were calculated using function richness (vegan) [70].VennDiagram (https://bioinfogp.cnb.csic.es/tools/venny/,accessed on 14 May 2023) has been used to illustrate dependencies between number of taxa in different tree development stages.

Statistical Analyses
A two-way analysis of variance model was used to compare fungal frequencies expressed in OTU units from different tree stages (seedling, small samplings, large saplings, and mature trees) and stands (Lipniak, Mikołajewo and Powały).This model was applied to analyse six dependent variables: Ascomycota: Basidiomycota ratio, fungi functional groups (ectomycorrhizas, pathogens, and saprotrophs), and two alpha diversity indexes (Chao, Shannon and Simpson).From a total of 58 independent observations a few outlying data points were removed for the purposes of the individual analyses.This information is emphasised in the tables containing the ANOVA results.Homogeneity of variance within fixed factors (tree groups, stands) was checked by the Levene F test using the function leveneTest(car ver.3.1-1).The goodness of fit of the data and residuals of ANOVA to a normal distribution was checked by Shapiro-Wilk test using the function shaprio.test(stats ver.4.2.2).To meet this requirement of parametric models two variables were transformed (by natural logarithm or square root), as stated in the tables providing the results of analysis by ANOVA.For statistically significant factors, least square means were calculated (function emmeans(emmeasn ver.1.8.2)) and the Tukey HSD post-hoc test was implemented to identify a significant difference between factor levels (function contrast(emmeans ver.1.8.2)).All statistical analyses were performed in R 4.2.2 [71].

Conclusions
The gradual decrease in abundance of saprotrophic species and increase of ectomycorrhizal in consecutive tree developmental stages indicates favorable conditions for spruce stand regeneration.The communities were dominated by saprotrophic and mycorrhizal species, with a small share of pathogenic species.The many species antagonistic to Heterobsasidion spp.and Armillaria spp.are beneficial for the growth and development of multi-generation Norway spruce stands in the Wigry National Park.

Figure 1 .
Figure 1.Mean and standard error of the mean for the relative abundance of taxa from the Ascomycota and Basidiomycota divisions in different tree developmental stages.The results of the Tukey HSD post hoc test are presented above the columns.Columns denoted by a different letter are significantly different (α = 0.05) for the Ascomycota/Basidiomycota.The results of the analysis of variance for the Ascomycota/Basidiomycota ratio are presented in Table1.

Figure 1 .
Figure 1.Mean and standard error of the mean for the relative abundance of taxa from the Ascomycota and Basidiomycota divisions in different tree developmental stages.The results of the Tukey HSD post hoc test are presented above the columns.Columns denoted by a different letter are significantly different (α = 0.05) for the Ascomycota/Basidiomycota.The results of the analysis of variance for the Ascomycota/Basidiomycota ratio are presented in Table1.

16 Figure 2 .
Figure 2. Venn diagram showing the number of taxa in different developmental stages.Treatments: (I) seedlings-current year plants with visible cotyledons; (II) small saplings up to 50 cm grown under the canopy of the mature stand, (III) large saplings over 50 cm in height, grown under the canopy of the mature stand, and (IV) mature trees.

Figure 2 .
Figure 2. Venn diagram showing the number of taxa in different developmental stages.Treatments: (I) seedlings-current year plants with visible cotyledons; (II) small saplings up to 50 cm grown under the canopy of the mature stand, (III) large saplings over 50 cm in height, grown under the canopy of the mature stand, and (IV) mature trees.

Figure 3 .
Figure 3. Mean and standard error of the mean for relative abundance of taxa classified into three functional groups: (a) ectomycorrhizas, (b) pathogens, and (c) saprotrophs.The results of the Tukey HSD post-hoc test are presented above the columns.Columns denoted by a different letter are significantly different (α = 0.05).The results of the analysis of variance are presented in Table2.

Figure 3 .
Figure 3. Mean and standard error of the mean for relative abundance of taxa classified into three functional groups: (a) ectomycorrhizas, (b) pathogens, and (c) saprotrophs.The results of the Tukey HSD post-hoc test are presented above the columns.Columns denoted by a different letter are significantly different (α = 0.05).The results of the analysis of variance are presented in Table2.

Figure 4 .
Figure 4. Heat map for the relative abundance of genera (the minimum threshold >1%).Figure 4. Heat map for the relative abundance of genera (the minimum threshold >1%).

Figure 4 .
Figure 4. Heat map for the relative abundance of genera (the minimum threshold >1%).Figure 4. Heat map for the relative abundance of genera (the minimum threshold >1%).

Table 3 .
Mean, standard error of the mean and ANOVA results for the Chao1 and Shannon diversity indices from different tree development stages (TDS) and stands (S).Df-degrees of freedom; F-the F-test value; P(α)-the probability of type I error.Levels not connected by the same letter are significantly different (Tukey's HSD test).