Next Article in Journal
Comparative Investigation of the Spectroscopic Behavior Based on High-Concentrated Solution in Nitrogen and Air Atmospheres
Next Article in Special Issue
Microbial Community Response to Various Types of Exogenous Organic Matter Applied to Soil
Previous Article in Journal
Enhancement of Arc Erosion Resistance in AgCuO Electrical Contact Materials through Rare Earth Element Doping: First-Principles and Experimental Studies
Previous Article in Special Issue
Tree Root-Associated Microbial Communities Depend on Various Floor Management Systems in an Intensive Apple (Malus × domestica Borkh.) Orchard
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 16
10.3390/ijms241612628
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Jolanta Behnke-Borowczyk
1,*,
Robert Korzeniewicz
1,
Adrian Łukowski
1,
Marlena Baranowska
1,
Radosław Jagiełło
1,
Bartosz Bułaj
1,
Maria Hauke-Kowalska
1,
Janusz Szmyt
1,
Jerzy M. Behnke
2,
Piotr Robakowski
1 and
Wojciech Kowalkowski
1
1
Faculty of Forestry and Wood Technology, Poznan University of Life Sciences, Wojska Polskiego 71C, 60-625 Poznań, Poland
2
School of Life Sciences, University Park Nottingham, Nottingham NG7 2RD, UK
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(16), 12628; https://doi.org/10.3390/ijms241612628
Submission received: 30 June 2023 / Revised: 28 July 2023 / Accepted: 4 August 2023 / Published: 10 August 2023
(This article belongs to the Special Issue Interaction between Plants, Microorganisms, and Soils in Different Ecosystems: 2nd Edition)
Browse Figures
Review Reports Versions Notes

Abstract

:
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.

1. 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.

2. 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.
Out of 339 Ascomycota taxa, we identified 298 in the rhizosphere of the seedling, 278 in small saplings, 289 in large sampling, and 282 in mature trees, and of these 243 were common in all stages of development. For Basidiomycota, 424 taxa were identified, of which 397 were in the rhizosphere of the seedings, 367 in small sampling, 390 in large sampling group, and 370 in mature trees, and 323 were common in all stages. The average number of taxa in the rhizosphere of the seedling was 783, in small saplings 726, in large saplings 765, and in mature trees 739, including 639 common taxa for all stages (Figure 2). The number of taxa found only in seedlings was 15. Five were found only in saplings up to 50 cm, nine in saplings above 50 cm, and 9 in mature trees (Figure 2).
It was observed that the abundance of genera, such as Amphinema (1.20–0.18%), Tylospora (0.34–0.006%), Thelephora (0.87–0.03%), (belonging to mycorrhizal), Armillaria (0.50–0.13%), Bryochiton (4.14–0.02%), Coniothyrium (1.96–0.20%), Drechslera (1.51–0.02%), Heterobasidion (1.13–0.52%), Scolecobasidium (0.21–0.00%), Zeloasperisporium (0.15–0.00%) (belonging to pathogens), Archaerhizomyces (0.03–0.49%), Basidiobolus (0.24–0.04%), Entoloma (1.70–0.47%), Exophiala (2.43–0.57%), Guttulispora (0.70–0.01%), Ischonoderma (0.70–0.09%), Kavina (1.13–0.09%), Kondoa (0.54–0.04%), Mycena (1.26–0.45%), Pyrenochaeta (11.78–1.39%), Stereum (1.20–0.18%), Thysanophora (0.36–0.06%), Tulasnella (1.57–0.29%), (belonging to saprotrophs) was highest in the seedlings or small saplings. The Amanita (0.33–0.04%), Inocybe (2.59–0.78%), Hygrophorus (0.56–0.00%), Paxillus (0.43–0.04%), Sceroderma (0.15–0.04%), Tomentella (15.94–4.54%) (belonging to mycorrhizal), Diplodia (2.59–1.57%), Foliporia (1.04–0.33%) (belonging to pathogens), Aureobasidium (0.36–0.08%), Cyanosporus (0.22–0.01%), Fibroporia (0.24–0.01%), Lycoperdon (2.11–0.13%), Sclerogaster (1.87–0.15%) (belonging to saprotrophs) was highest in large saplings or mature tree (Figure 3).
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.

3. 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 mycorrhiza-formers, 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.

4. Materials and Methods

4.1. 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).

4.2. 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.

4.3. 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.

4.4. 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.

4.5. 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].

5. 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.

Author Contributions

Conceptualization, J.B.-B., M.B. and W.K.; methodology, J.B.-B. and M.B.; software, J.B.-B. and R.J.; formal analysis, J.B.-B., R.J., R.K. and A.Ł.; investigation, J.B.-B., R.K., A.Ł., M.B., R.J., B.B., M.H.-K., J.S., P.R. and W.K.; resources, J.B.-B.; writing—original draft preparation, J.B.-B.; writing—review and editing, J.M.B., R.K., A.Ł., M.B., R.J., B.B., M.H.-K., J.S., P.R. and W.K.; visualization, R.J. and J.B.-B.; supervision, P.R., W.K. and J.M.B.; project administration, W.K. and J.B.-B. funding acquisition, J.B.-B.; R.K., A.Ł., M.B., R.J., B.B., M.H.-K., J.S., P.R. and W.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The State Forests (contract no. EZ.0290.1.21.2022 of 23 August 2022) under the research grant entitled “Genetic structure and population dynamics of Norway spruce in the changing environment of the Wigry National Park”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data supporting reported results can be found at https://doi.org/10.6084/m9.figshare.23733399.v1 (accessed on 23 July 2023).

Acknowledgments

We warmly thank workers of the Wigry National Park for help in the organisation of our research and providing all facilites in the Wigry National Park.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Philippot, L.; Raaijmakers, J.M.; Lemanceau, P.; van der Putten, W.H. Going Back to the Roots: The Microbial Ecology of the Rhizosphere. Nat. Rev. Microbiol. 2013, 11, 789–799. [Google Scholar] [CrossRef] [PubMed]
  2. Turner, T.R.; Ramakrishnan, K.; Walshaw, J.; Heavens, D.; Alston, M.; Swarbreck, D.; Osbourn, A.; Grant, A.; Poole, P.S. Comparative Metatranscriptomics Reveals Kingdom Level Changes in the Rhizosphere Microbiome of Plants. ISME J. 2013, 7, 2248–2258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Carroll, G. Fungal Endophytes in Stems and Leaves: From Latent Pathogen to Mutualistic Symbiont. Ecology 1988, 69, 2–9. [Google Scholar] [CrossRef]
  4. Menkis, A.; Burokienė, D.; Stenlid, J.; Stenström, E. High-Throughput Sequencing Shows High Fungal Diversity and Community Segregation in the Rhizospheres of Container-Grown Conifer Seedlings. Forests 2016, 7, 44. [Google Scholar] [CrossRef] [Green Version]
  5. Otsing, E.; Barantal, S.; Anslan, S.; Koricheva, J.; Tedersoo, L. Litter Species Richness and Composition Effects on Fungal Richness and Community Structure in Decomposing Foliar and Root Litter. Soil. Biol. Biochem. 2018, 125, 328–339. [Google Scholar] [CrossRef]
  6. LeBlanc, N.; Kinkel, L.L.; Kistler, H.C. Soil Fungal Communities Respond to Grassland Plant Community Richness and Soil Edaphics. Microb. Ecol. 2015, 70, 188–195. [Google Scholar] [CrossRef]
  7. Peay, K.G.; Kennedy, P.G.; Talbot, J.M. Dimensions of Biodiversity in the Earth Mycobiome. Nat. Rev. Microbiol. 2016, 14, 434–447. [Google Scholar] [CrossRef]
  8. Tedersoo, L.; Bahram, M.; Zobel, M. How Mycorrhizal Associations Drive Plant Population and Community Biology. Science 2020, 367, eaba1223. [Google Scholar] [CrossRef]
  9. Khoshru, B.; Moharramnejad, S.; Gharajeh, N.H.; Asgari Lajayer, B.; Ghorbanpour, M. Plant Microbiome and Its Important in Stressful Agriculture. In Plant Microbiome Paradigm; Springer International Publishing: Cham, Switzerland, 2020; pp. 13–48. [Google Scholar]
  10. Lindahl, B.D.; Tunlid, A. Ectomycorrhizal Fungi—Potential Organic Matter Decomposers, yet Not Saprotrophs. New Phytol. 2015, 205, 1443–1447. [Google Scholar] [CrossRef]
  11. Van der Heijden, M.G.A.; Horton, T.R. Socialism in Soil? The Importance of Mycorrhizal Fungal Networks for Facilitation in Natural Ecosystems. J. Ecol. 2009, 97, 1139–1150. [Google Scholar] [CrossRef]
  12. Figueiredo, A.F.; Boy, J.; Guggenberger, G. Common Mycorrhizae Network: A Review of the Theories and Mechanisms Behind Underground Interactions. Front. Fungal Biol. 2021, 2, 48. [Google Scholar] [CrossRef]
  13. Rhodes, C.J. The Whispering World of Plants: ‘The Wood Wide Web’. Sci. Prog. 2017, 100, 331–337. [Google Scholar] [CrossRef] [PubMed]
  14. Wipf, D.; Krajinski, F.; van Tuinen, D.; Recorbet, G.; Courty, P. Trading on the Arbuscular Mycorrhiza Market: From Arbuscules to Common Mycorrhizal Networks. New Phytol. 2019, 223, 1127–1142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Deslippe, J.R.; Hartmann, M.; Simard, S.W.; Mohn, W.W. Long-Term Warming Alters the Composition of Arctic Soil Microbial Communities. FEMS Microbiol. Ecol. 2012, 82, 303–315. [Google Scholar] [CrossRef] [Green Version]
  16. Bücking, H.; Mensah, J.A.; Fellbaum, C.R. Common Mycorrhizal Networks and Their Effect on the Bargaining Power of the Fungal Partner in the Arbuscular Mycorrhizal Symbiosis. Commun. Integr. Biol. 2016, 9, e1107684. [Google Scholar] [CrossRef]
  17. Bingham, M.A.; Simard, S.W. Do Mycorrhizal Network Benefits to Survival and Growth of Interior Douglas-Fir Seedlings Increase with Soil Moisture Stress? Ecol. Evol. 2011, 1, 306–316. [Google Scholar] [CrossRef]
  18. Seiwa, K.; Negishi, Y.; Eto, Y.; Hishita, M.; Masaka, K.; Fukasawa, Y.; Matsukura, K.; Suzuki, M. Successful Seedling Establishment of Arbuscular Mycorrhizal-Compared to Ectomycorrhizal-Associated Hardwoods in Arbuscular Cedar Plantations. For. Ecol. Manag. 2020, 468, 118155. [Google Scholar] [CrossRef]
  19. Cairney, J.W.; Meharg, A.A. Interactions between Ectomycorrhizal Fungi and Soil Saprotrophs: Implications for Decomposition of Organic Matter in Soils and Degradation of Organic Pollutants in the Rhizosphere. Can. J. Bot. 2002, 80, 803–809. [Google Scholar] [CrossRef]
  20. Termorshuizen, A.J. Ecology of Fungal Plant Pathogens. Microbiol. Spectr. 2016, 4. [Google Scholar] [CrossRef]
  21. Berendsen, R.L.; Pieterse, C.M.J.; Bakker, P.A.H.M. The Rhizosphere Microbiome and Plant Health. Trends Plant Sci. 2012, 17, 478–486. [Google Scholar] [CrossRef]
  22. Sharma, S.; Ghoshal, C.; Arora, A.; Samar, W.; Nain, L.; Paul, D. Strain Improvement of Native Saccharomyces Cerevisiae LN ITCC 8246 Strain Through Protoplast Fusion To Enhance Its Xylose Uptake. Appl. Biochem. Biotechnol. 2021, 193, 2455–2469. [Google Scholar] [CrossRef] [PubMed]
  23. Korenblum, E.; Dong, Y.; Szymanski, J.; Panda, S.; Jozwiak, A.; Massalha, H.; Meir, S.; Rogachev, I.; Aharoni, A. Rhizosphere Microbiome Mediates Systemic Root Metabolite Exudation by Root-to-Root Signaling. Proc. Natl. Acad. Sci. USA 2020, 117, 3874–3883. [Google Scholar] [CrossRef] [PubMed]
  24. Peay, K.G.; Belisle, M.; Fukami, T. Phylogenetic Relatedness Predicts Priority Effects in Nectar Yeast Communities. Proc. R. Soc. B Biol. Sci. 2012, 279, 749–758. [Google Scholar] [CrossRef]
  25. Nara, K.; Nakaya, H.; Wu, B.; Zhou, Z.; Hogetsu, T. Underground Primary Succession of Ectomycorrhizal Fungi in a Volcanic Desert on Mount Fuji. New Phytol. 2003, 159, 743–756. [Google Scholar] [CrossRef]
  26. Gossner, M.M.; Getzin, S.; Lange, M.; Pašalić, E.; Türke, M.; Wiegand, K.; Weisser, W.W. The Importance of Heterogeneity Revisited from a Multiscale and Multitaxa Approach. Biol. Conserv. 2013, 166, 212–220. [Google Scholar] [CrossRef]
  27. Wright, J.S. Plant Diversity in Tropical Forests: A Review of Mechanisms of Species Coexistence. Oecologia 2002, 130, 1–14. [Google Scholar] [CrossRef] [PubMed]
  28. Behnke-Borowczyk, J.; Kowalkowski, W.; Kartawik, N.; Baranowska, M.; Barzdajn, W. The Soil Fungal Communities in Nurseries Producing Abies Alba. Balt. For. 2020, 26, 426. [Google Scholar] [CrossRef]
  29. Seppä, H.; Alenius, T.; Bradshaw, R.H.W.; Giesecke, T.; Heikkilä, M.; Muukkonen, P. Invasion of Norway Spruce (Picea Abies) and the Rise of the Boreal Ecosystem in Fennoscandia. J. Ecol. 2009, 97, 629–640. [Google Scholar] [CrossRef]
  30. Jevšenak, J.; Tychkov, I.; Gričar, J.; Levanič, T.; Tumajer, J.; Prislan, P.; Arnič, D.; Popkova, M.; Shishov, V.V. Growth-Limiting Factors and Climate Response Variability in Norway Spruce (Picea Abies L.) along an Elevation and Precipitation Gradients in Slovenia. Int. J. Biometeorol. 2021, 65, 311–324. [Google Scholar] [CrossRef]
  31. Hentschel, R.; Rosner, S.; Kayler, Z.E.; Andreassen, K.; Børja, I.; Solberg, S.; Tveito, O.E.; Priesack, E.; Gessler, A. Norway Spruce Physiological and Anatomical Predisposition to Dieback. For. Ecol. Manag. 2014, 322, 27–36. [Google Scholar] [CrossRef]
  32. Tedersoo, L.; Bahram, M.; Põlme, S.; Kõljalg, U.; Yorou, N.S.; Wijesundera, R.; Ruiz, L.V.; Vasco-Palacios, A.M.; Thu, P.Q.; Suija, A.; et al. Global Diversity and Geography of Soil Fungi. Science 2014, 346, 1256688. [Google Scholar] [CrossRef] [Green Version]
  33. Nguyen, N.H.; Song, Z.; Bates, S.T.; Branco, S.; Tedersoo, L.; Menke, J.; Schilling, J.S.; Kennedy, P.G. FUNGuild: An Open Annotation Tool for Parsing Fungal Community Datasets by Ecological Guild. Fungal Ecol. 2016, 20, 241–248. [Google Scholar] [CrossRef]
  34. Zanne, A.E.; Abarenkov, K.; Afkhami, M.E.; Aguilar-Trigueros, C.A.; Bates, S.; Bhatnagar, J.M.; Busby, P.E.; Christian, N.; Cornwell, W.K.; Crowther, T.W.; et al. Fungal Functional Ecology: Bringing a Trait-based Approach to Plant-associated Fungi. Biol. Rev. 2020, 95, 409–433. [Google Scholar] [CrossRef] [PubMed]
  35. Kohout, P.; Charvátová, M.; Štursová, M.; Mašínová, T.; Tomšovský, M.; Baldrian, P. Clearcutting Alters Decomposition Processes and Initiates Complex Restructuring of Fungal Communities in Soil and Tree Roots. ISME J. 2018, 12, 692–703. [Google Scholar] [CrossRef] [Green Version]
  36. Baldrian, P.; Kolařík, M.; Štursová, M.; Kopecký, J.; Valášková, V.; Větrovský, T.; Žifčáková, L.; Šnajdr, J.; Rídl, J.; Vlček, Č.; et al. Active and Total Microbial Communities in Forest Soil Are Largely Different and Highly Stratified during Decomposition. ISME J. 2012, 6, 248–258. [Google Scholar] [CrossRef] [Green Version]
  37. Kwaśna, H.; Behnke-Borowczyk, J.; Gornowicz, R.; Łakomy, P. Effects of Preparation of Clear-cut Forest Sites on the Soil Mycobiota with Consequences for Scots Pine Growth and Health. For. Pathol. 2019, 49, e12494. [Google Scholar] [CrossRef]
  38. Kwaśna, H.; Łakomy, P.; Gornowicz, R.; Borowczyk-Behnke, J.; Kuźmiński, R. Wpływ Sposobu Przygotowania Gleby Na Aktywność Biologiczną Gleby Względem Patogenów Korzeni w 40-Letnim Drzewostanie Sosnowym. Sylwan 2015, 159, 125–177. [Google Scholar]
  39. Hilszczańska, D.; Sierota, Z. Wplyw Inokulum Mikoryzowego Grzyba Thelephora Terrestris Na Wzrost Sadzonek Sosny Zwyczajnej Pinus Sylvestris L. II. Badania Polowe. Sylwan 2006, 150, 20–28. [Google Scholar]
  40. Hayward, J.; Horton, T.R.; Nuñez, M.A. Ectomycorrhizal Fungal Communities Coinvading with Pinaceae Host Plants in Argentina: Gringos Bajo El Bosque. New Phytol. 2015, 208, 497–506. [Google Scholar] [CrossRef] [PubMed]
  41. Guo, Q.; Duan, Y.; Meng, N.; Liu, Y.; Luo, G. The N-Terminus of Sec3 Is Required for Cell Wall Integrity in Yeast. Biochimie 2020, 177, 30–39. [Google Scholar] [CrossRef]
  42. Hibbett, D.S.; Binder, M. Evolution of Complex Fruiting–Body Morphologies in Homobasidiomycetes. Proc. R. Soc. Lond. B Biol. Sci. 2002, 269, 1963–1969. [Google Scholar] [CrossRef] [PubMed]
  43. Zmitrovich, I. Conspectus systematis Polyporacearum v. 1.0. Folia Cryptogam. Petropolitanfa 2018, 6, 3–145. [Google Scholar]
  44. Dettman, J.R.; Eggertson, Q. Phylogenomic Analyses of Alternaria Section Alternaria: A High-Resolution, Genome-Wide Study of Lineage Sorting and Gene Tree Discordance. Mycologia 2021, 113, 1218–1232. [Google Scholar] [CrossRef] [PubMed]
  45. Zhang, S.; Zhang, J.; Wang, C.; Chen, X.; Zhao, X.; Jing, H.; Liu, H.; Li, Z.; Wang, L.; Shi, J. COVID-19 and Ischemic Stroke: Mechanisms of Hypercoagulability (Review). Int. J. Mol. Med. 2021, 47, 21. [Google Scholar] [CrossRef] [PubMed]
  46. Simmons, E.G. Alternaria Themes and Variations (287–304): Species on Caryophyllaceae. Mycotaxon 2002, 82, 1–40. [Google Scholar]
  47. Díaz, G.A.; Mostert, L.; Halleen, F.; Lolas, M.; Gutierrez, M.; Ferrada, E.; Latorre, B.A. Diplodia Seriata Associated with Botryosphaeria Canker and Dieback in Apple Trees in Chile. Plant Dis. 2019, 103, 1025. [Google Scholar] [CrossRef]
  48. Elena, G.; Luque, J. Seasonal Susceptibility of Grapevine Pruning Wounds and Cane Colonization in Catalonia, Spain Following Artificial Infection with Diplodia Seriata and Phaeomoniella Chlamydospora. Plant Dis. 2016, 100, 1651–1659. [Google Scholar] [CrossRef] [Green Version]
  49. Swart, W.J. and Wingfield, M.J. Biology and Control of Sphaeropsis sapinea on Pinus species in South Africa. Plant Dis. 1991, 75, 761–766. [Google Scholar] [CrossRef]
  50. Oliva, J.; Rommel, S.; Fossdal, C.G.; Hietala, A.M.; Nemesio-Gorriz, M.; Solheim, H.; Elfstrand, M. Transcriptional Responses of Norway Spruce (Picea Abies) Inner Sapwood against Heterobasidion Parviporum. Tree Physiol. 2015, 35, 1007–1015. [Google Scholar] [CrossRef] [Green Version]
  51. Kovalchuk, A.; Wen, Z.; Sun, H.; Asiegbu, F.O. Heterobasidion Annosum s.l.: Biology, Genomics, and Pathogenicity Factors. In Forest Microbiology; Elsevier: Amsterdam, The Netherlands, 2022; pp. 345–359. [Google Scholar]
  52. Sierota, Z. Heterobasidion Root Rot in Forests on Former Agricultural Lands in Poland: Scale of Threat and Prevention. Sci. Res. Essays 2013, 8, 2298–2305. [Google Scholar]
  53. Misiek, M.; Hoffmeister, D. Sesquiterpene Aryl Ester Natural Products in North American Armillaria Species. Mycol. Prog. 2012, 11, 7–15. [Google Scholar] [CrossRef]
  54. Zeiner, C.A.; Purvine, S.O.; Zink, E.M.; Paša-Tolić, L.; Chaput, D.L.; Haridas, S.; Wu, S.; LaButti, K.; Grigoriev, I.V.; Henrissat, B.; et al. Comparative Analysis of Secretome Profiles of Manganese(II)-Oxidizing Ascomycete Fungi. PLoS ONE 2016, 11, e0157844. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Sieber, T.N. Pyrenochaeta Ligni-Putridi Sp. Nov., a New Coelomycete Associated with Butt Rot of Picea Abies in Switzerland. Mycol. Res. 1995, 99, 274–276. [Google Scholar] [CrossRef]
  56. Infantino, A.; Pucci, N. A PCR-Based Assay for the Detection and Identification of Pyrenochaeta Lycopersici. Eur. J. Plant Pathol. 2005, 112, 337–347. [Google Scholar] [CrossRef]
  57. Butin, H. Tree Diseases and Disorders: Causes, Biology, and Control in Forest and Amenity Trees; Oxford University Press: Oxford, UK, 1995. [Google Scholar]
  58. Sivanesan, A. Acantharia, Gibbera and Their Anamorphs. Trans. Br. Mycol. Soc. 1984, 82, 507–529. [Google Scholar] [CrossRef]
  59. Kowalski, T.; Andruch, K. Mycobiota in Needles of Abies Alba with and without Symptoms of Herpotrichia Needle Browning. For. Pathol. 2012, 42, 183–190. [Google Scholar] [CrossRef]
  60. Barr, M.E. Notes on Some “Dimeriaceous” Fungi. Mycotaxon 1997, 64, 149–171. [Google Scholar]
  61. Harman, G.; Khadka, R.; Doni, F.; Uphoff, N. Benefits to Plant Health and Productivity From Enhancing Plant Microbial Symbionts. Front. Plant Sci. 2021, 11, 610065. [Google Scholar] [CrossRef]
  62. Schmidt, S.; Moser, M.; Sperandio, M. The Molecular Basis of Leukocyte Recruitment and Its Deficiencies. Mol. Immunol. 2013, 55, 49–58. [Google Scholar] [CrossRef]
  63. Bálint, M.; Schmidt, P.; Sharma, R.; Thines, M.; Schmitt, I. An Illumina Metabarcoding Pipeline for Fungi. Ecol. Evol. 2014, 4, 2642–2653. [Google Scholar] [CrossRef]
  64. Vilgalys, R.; Hester, M. Rapid Genetic Identification and Mapping of Enzymatically Amplified Ribosomal DNA from Several Cryptococcus Species. J. Bacteriol. 1990, 172, 4238–4246. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Aronesty, E. Comparison of Sequencing Utility Programs. Open Bioinforma J. 2013, 7, 1–8. [Google Scholar] [CrossRef] [Green Version]
  66. Benson, D.A.; Cavanaugh, M.; Clark, K.; Karsch-Mizrachi, I.; Ostell, J.; Pruitt, K.D.; Sayers, E.W. GenBank. Nucleic Acids Res. 2018, 46, D41–D47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Kubartová, A.; Ottosson, E.; Dahlberg, A.; Stenlid, J. Patterns of Fungal Communities among and within Decaying Logs, Revealed by 454 Sequencing. Mol. Ecol. 2012, 21, 4514–4532. [Google Scholar] [CrossRef]
  68. Raja, H.A.; Miller, A.N.; Pearce, C.J.; Oberlies, N.H. Fungal Identification Using Molecular Tools: A Primer for the Natural Products Research Community. J. Nat. Prod. 2017, 80, 756–770. [Google Scholar] [CrossRef]
  69. Põlme, S.; Abarenkov, K.; Henrik Nilsson, R.; Lindahl, B.D.; Clemmensen, K.E.; Kauserud, H.; Nguyen, N.; Kjøller, R.; Bates, S.T.; Baldrian, P.; et al. FungalTraits: A User-Friendly Traits Database of Fungi and Fungus-like Stramenopiles. Fungal Divers. 2020, 105, 1–16. [Google Scholar] [CrossRef]
  70. Oksanen, J.; Friendly, M.; Kindt, R. Measurement of Biodiversity (MoB) View Project. Available online: https://cran.r-project.org/web/packages/vegan/vegan.pdf (accessed on 1 May 2023).
  71. R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. 2022. Available online: https://www.R-project.org (accessed on 3 May 2023).
Ijms 24 12628 g001
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 Table 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 Table 1.
Ijms 24 12628 g001
Ijms 24 12628 g002
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. 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.
Ijms 24 12628 g002
Ijms 24 12628 g003
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 Table 2.
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 Table 2.
Ijms 24 12628 g003
Ijms 24 12628 g004
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%).
Ijms 24 12628 g004
Table 1. Results of the two-way analysis of variance for OTU Ascomycota: Basidiomycota ratio from different stands (S) and tree development stage (TDS) Df—degrees of freedom; F—the F-test value; P(α)—the probability of type I error.
Table 1. Results of the two-way analysis of variance for OTU Ascomycota: Basidiomycota ratio from different stands (S) and tree development stage (TDS) Df—degrees of freedom; F—the F-test value; P(α)—the probability of type I error.
Source of VariationDfFP (α)
TDS35.4740.003
S22.0550.140
TDS × S60.4040.872
Table 2. Results of the analysis of variance for mean OTU of taxa classified to three functional groups (ectomycorrhizas, pathogens and saprotrophs) from different stands (S) and tree development stage (TDS). Df—degrees of freedom; F—the F-test value; P(α)—the probability of type I error.
Table 2. Results of the analysis of variance for mean OTU of taxa classified to three functional groups (ectomycorrhizas, pathogens and saprotrophs) from different stands (S) and tree development stage (TDS). Df—degrees of freedom; F—the F-test value; P(α)—the probability of type I error.
Source of VariationDfEctomycorrhizasPathogensSaprotrophs
FP (α)FP (α)FP (α)
TDS35.6160.0023.9420.0142.5980.063
S21.3630.2660.9840.3825.5220.007
TDS × S61.5380.1871.7660.1271.7920.122
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).
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).
Tree
Development Stage		Chao1ShannonSimpson
Mean (Standard Error)
Seedlings718 (61) b3.91 (0.230)0.986 (0.044)
Small saplings1213 (128) a4.06 (0.209)0.923 (0.023)
Large saplings845 (76) b3.92 (0.195)0.901 (0.036)
Mature trees891 (89) ab4.23 (0.219)0.932 (0.030)
Analysis of Variance
Source of
Variation		DfFP (α)FP (α)FP (α)
TDS35.3070.0032.2490.0970.2700.847
S21.2960.2832.4770.0960.1220.886
TDS × S61.0460.4082.1260.0712.0120.083
Table 4. Vertical structure and species composition of the stands involved in the research. Ap—Acer platanoides L., Bp—Betula pendula Roth, Ca—Corylus avellana L., Ev—Euonymus verrucosus Scop, Fa—Frangula alnus Mill., Jc—Juniperus communis L., Lx—Lonicera xylosteum L., Pa—Picea abies (L.) H.Karst., Ps—Pinus sylvestris L., Pt – Populus tremula L., Qr—Quercus robur L., Sa—Sorbus aucuparia L., Tc—Tilia cordata Mill.; ns—not specified; * Data source: Forest Data Bank (https://www.bdl.lasy.gov.pl/portal/en), accessed on 30 May 2023. Data in the table are derived from the last update for the year 2013. The species are indicated as a gradient of decreasing share within the stratum.
Table 4. Vertical structure and species composition of the stands involved in the research. Ap—Acer platanoides L., Bp—Betula pendula Roth, Ca—Corylus avellana L., Ev—Euonymus verrucosus Scop, Fa—Frangula alnus Mill., Jc—Juniperus communis L., Lx—Lonicera xylosteum L., Pa—Picea abies (L.) H.Karst., Ps—Pinus sylvestris L., Pt – Populus tremula L., Qr—Quercus robur L., Sa—Sorbus aucuparia L., Tc—Tilia cordata Mill.; ns—not specified; * Data source: Forest Data Bank (https://www.bdl.lasy.gov.pl/portal/en), accessed on 30 May 2023. Data in the table are derived from the last update for the year 2013. The species are indicated as a gradient of decreasing share within the stratum.
Protection Zone, Compartment, GPS LocationArea
[ha]		Forest
Stratum		Cover
[%]		Species Composition *
[Age in Years]
Lipniak,
15f
54°08′53.4″ N 23°05′15.5″ E		9.75Overstorey80Pa (116), Ps (116)
II stratum20Pa (75), Pa (50), Bp (50), Tc (50)
Large saplings10Pa (25)
Small saplingsnsPa (5), Qr (3), Ap (3), Tc (3)
Undergrowth50Lx, Ca, Fa, Sa, Ev
Mikołajewo,
210c
54°02′10.5″ N 23°08′55.0″ E		6.55Overstorey70Pa (141), Ps (141), Pa (101)
II stratum30Pa (76), Pa (46), Qr (46), Pt (46)
III stratum20Pa (31), Qr (31)
Large saplings30Pa (21), Pa (11), Qr (21)
Small saplings10Pa (4), Qr (4)
Undergrowth 30Ca, Sa, Lx, Fa, Ev
Powały,
390c
53°59′59.7″ N 23°02′12.9″ E		7.98Overstorey60Ps (121), Pa (85), Pa (121), Bp (121),
Bp (85), Qr (85)
II stratum30Pa (35), Pa (50), Bp (35), Bv (50),
Qr (50), Jc (35)
Large saplings10Pa (15), Bp (15), Qr (15)
Small saplingsnsPa (3), Bp (3), Qr (3)
Undergrowth30Ca, Sa
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Behnke-Borowczyk, J.; Korzeniewicz, R.; Łukowski, A.; Baranowska, M.; Jagiełło, R.; Bułaj, B.; Hauke-Kowalska, M.; Szmyt, J.; Behnke, J.M.; Robakowski, P.; et al. Variability of Functional Groups of Rhizosphere Fungi of Norway Spruce (Picea abies (L.) H.Karst.) in the Boreal Range: The Wigry National Park, Poland. Int. J. Mol. Sci. 2023, 24, 12628. https://doi.org/10.3390/ijms241612628

AMA Style

Behnke-Borowczyk J, Korzeniewicz R, Łukowski A, Baranowska M, Jagiełło R, Bułaj B, Hauke-Kowalska M, Szmyt J, Behnke JM, Robakowski P, et al. Variability of Functional Groups of Rhizosphere Fungi of Norway Spruce (Picea abies (L.) H.Karst.) in the Boreal Range: The Wigry National Park, Poland. International Journal of Molecular Sciences. 2023; 24(16):12628. https://doi.org/10.3390/ijms241612628

Chicago/Turabian Style

Behnke-Borowczyk, Jolanta, Robert Korzeniewicz, Adrian Łukowski, Marlena Baranowska, Radosław Jagiełło, Bartosz Bułaj, Maria Hauke-Kowalska, Janusz Szmyt, Jerzy M. Behnke, Piotr Robakowski, and et al. 2023. "Variability of Functional Groups of Rhizosphere Fungi of Norway Spruce (Picea abies (L.) H.Karst.) in the Boreal Range: The Wigry National Park, Poland" International Journal of Molecular Sciences 24, no. 16: 12628. https://doi.org/10.3390/ijms241612628

APA Style

Behnke-Borowczyk, J., Korzeniewicz, R., Łukowski, A., Baranowska, M., Jagiełło, R., Bułaj, B., Hauke-Kowalska, M., Szmyt, J., Behnke, J. M., Robakowski, P., & Kowalkowski, W. (2023). Variability of Functional Groups of Rhizosphere Fungi of Norway Spruce (Picea abies (L.) H.Karst.) in the Boreal Range: The Wigry National Park, Poland. International Journal of Molecular Sciences, 24(16), 12628. https://doi.org/10.3390/ijms241612628

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop