Next Article in Journal
Description of Four Novel Species in Pleosporales Associated with Coffee in Yunnan, China
Previous Article in Journal
Roles of BrlA and AbaA in Mediating Asexual and Insect Pathogenic Lifecycles of Metarhizium robertsii
Previous Article in Special Issue
Diversity and Distribution of Macrofungi in Protected Mountain Forest Habitats in Serbia and Its Relation to Abiotic Factors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Principal Drivers of Fungal Communities Associated with Needles, Shoots, Roots and Adjacent Soil of Pinus sylvestris

by
Diana Marčiulynienė
1,*,
Adas Marčiulynas
1,
Valeriia Mishcherikova
1,
Jūratė Lynikienė
1,
Artūras Gedminas
1,
Iva Franic
2 and
Audrius Menkis
3
1
Institute of Forestry, Lithuanian Research Centre for Agriculture and Forestry, Liepų Str. 1, Girionys, 53101 Kaunas District, Lithuania
2
Department of Southern Swedish Forest Research Centre, Swedish University of Agricultural Sciences, P.O. Box 190, SE-23422 Lomma, Sweden
3
Department of Forest Mycology and Plant Pathology, Uppsala BioCenter, Swedish University of Agricultural Sciences, P.O. Box 7026, SE-75007 Uppsala, Sweden
*
Author to whom correspondence should be addressed.
J. Fungi 2022, 8(10), 1112; https://doi.org/10.3390/jof8101112
Submission received: 6 September 2022 / Revised: 12 October 2022 / Accepted: 18 October 2022 / Published: 21 October 2022
(This article belongs to the Special Issue Fungal Diversity in Europe)

Abstract

:
The plant- and soil-associated microbial communities are critical to plant health and their resilience to stressors, such as drought, pathogens, and pest outbreaks. A better understanding of the structure of microbial communities and how they are affected by different environmental factors is needed to predict and manage ecosystem responses to climate change. In this study, we carried out a country-wide analysis of fungal communities associated with Pinus sylvestris growing under different environmental conditions. Needle, shoot, root, mineral, and organic soil samples were collected at 30 sites. By interconnecting the high-throughput sequencing data, environmental variables, and soil chemical properties, we were able to identify key factors that drive the diversity and composition of fungal communities associated with P. sylvestris. The fungal species richness and community composition were also found to be highly dependent on the site and the substrate they colonize. The results demonstrated that different functional tissues and the rhizosphere soil of P. sylvestris are associated with diverse fungal communities, which are driven by a combination of climatic (temperature and precipitation) and edaphic factors (soil pH), and stand characteristics.

1. Introduction

Changing climate induces multiple abiotic and biotic risks to forests and forestry worldwide. Indeed, climate change affects the location, composition, structure and function of forests in many parts of the world, including high-latitude forests dominated by boreal species [1]. Nevertheless, forests may respond to climate change in various ways driven by local environmental conditions and the adaptive potential of the trees. Climate change is also changing the local environments, affecting local microbial communities and, thus, their metabolic activity and capacity [2]. Yet, climate, tree composition, or the distribution of the host tree do not always explain spatial patterns in the above- and below-ground microbial communities. This is important to consider because different parts of terrestrial plants, including forest trees, harbor several communities of microorganisms that form complex associations and play important roles that can determine the functioning and health of their hosts. There are many ways to tackle the challenges of climate change, including adaptation and mitigation strategies. However, little attention has been given to microbial adaptation [3]. Although microorganisms tend to adapt rapidly to changing environmental conditions [4], little is known about how these changes will feed back into plant–microbe interactions [5]. The overall effects of plant-associated microbes on host health and fitness are determined by a number of factors, including host and microbial genotypes, interactions within microbiota, and various abiotic factors [6]. The main abiotic factors that determine the diversity of tree microbiome include temperature, light, nutrient and water availability, and soil pH [7]. However, the role of abiotic factors in determining the microbiome structure of forest trees, especially in different habitats, is not well understood [8]. Tree-associated microbes occupy specific habitats, which are referred to as rhizosphere, rhizoplane, phyllosphere, or endosphere [9]. Microbes in these habitats can establish beneficial, neutral, or detrimental associations of varying intimacy with their hosts [10]. Fungi are among the dominant groups of plant-associated microorganisms as they play key roles in regulating plant health, maintaining interactions between plants and other organisms, and maintaining the functioning of entire ecosystems [11,12]. Maintaining functional tree-microbial interactions has been shown to be important for tree growth and development and for adaptation to the changing environment [13]. Given that different fungal species and functional groups are driven differently by climatic, nutritional, or biotic factors, it is expected that different fungal guilds may also react differently to changes in environmental conditions [14].
Some of the most complex microbiomes are found in soils [15]. Fungi are of particular interest as they represent a significant fraction of the soil microbial community and influence several ecosystem processes [16]. The abundance of fungal functional groups in the soil changes in response to forest disturbance and indicates a close interaction between the above-ground plant community and the below-ground soil biological community [14]. Moreover, structural and functional modifications in the soil microbiota have a crucial impact on above-ground ecosystems [17]. A global meta-analysis found that forest degradation reduces soil C and N content, increases soil pH, and increases C decomposition rates. The study also found a decrease in soil fungal biomass in disturbed sites but increased species diversity [18]. Changes in soil pH were also shown to be significantly correlated with changes in the fungal community composition in the soil [19]. Besides, the loss of host species and forest disturbance can also cause the direct loss of mycorrhizal and other rhizosphere fungi [20].
Several previous studies have examined the effects of individual disturbances on soil fungal composition, focusing on ectomycorrhizal (ECM) fungi [21,22,23,24] but largely excluding other fungal guilds with which they interact and synergize to perform multiple ecosystem processes [25,26]. ECM fungi are widely associated with trees in natural forests and plantations [27]. They are a key component of the forest soil microbiome, forming symbioses with ca. 60% of trees on Earth [28]. Their importance increases with the decrease of soil fertility and under harsh environmental conditions as they improve the host trees’ nutrition and stress resistance [20,29]. As previous studies have shown, fungal communities associated with tree roots show pronounced differences compared with the communities of the surrounding soil [7]. Differences between root and rhizosphere microbial communities have been found to be influenced by the seasonal production of roots and a variation in their activity [30] and stand age [31]. Root-associated fungal communities are also known to depend on environmental parameters [32,33,34], such as climatic conditions [35], soil chemical composition, especially pH, and the content of organic matter [36]. Soil and/or root-associated fungi can also be pathogenic [37] or interact as saprotrophs involved in nutrient recycling [38]. Pathogenic fungi can interfere with tree growth and alter the diversity and composition of the plant community [39].
Different above-ground functional tissues of forest trees can also host diverse communities of fungi [7]. For example, living needles of coniferous trees are associated with a diverse fungal community, and these fungi may have various effects on their host [40]. Such fungi commonly include endophytes and epiphytes that colonize the interior and exterior surfaces of living needles, respectively [41]. Besides, their abundances in needles vary depending on seasonal precipitation and temperature [42]. There are several factors that influence the composition of the foliar fungal community. These factors include the host species, nutrient content in the needles, needle age, air temperature, precipitation, and air pollution [43,44]. Recent studies also show that the fungal community of pine trees growing in a common environment show host species-specific structures [45]. The importance of endophytic and other wood-inhabiting fungi has also been emphasized. Endophytes are common in living tree tissues [46,47,48] and may influence many biotic processes in trees, but under certain conditions, they may shift to pathogenic or saprophytic lifestyles [7]. Therefore, they can also be involved in the early decomposition of deadwood, which is an important carbon stock in forest ecosystems [49,50] and an important habitat for many forest living organisms [51].
As fungi play a vital role in the structural and functional dynamics of terrestrial ecosystems [52], it is important to understand how the above- and below-ground fungal communities will respond to changing environmental conditions. Furthermore, exploring the impacts of local variables on the structure of fungal communities can help to better understand the ecological functions of fungi in different habitats. The focus of the present study was a Scots pine (Pinus sylvestris L.), which is one of the most economically and ecologically important tree species in Europe [53]. Still, climate change is likely to affect its survival and growth differently in different parts of its distribution range [54]. The aim of the present study was to study the diversity and composition of above- (needles, shoots) and below-ground (roots, soil) fungal communities associated with P. sylvestris, and to evaluate the impact of abiotic factors on these fungal communities. We hypothesized that above- and below-ground fungal communities will change in quantity and composition along with the changes in environmental factors (temperature, precipitation, soil chemical properties) and stand age. Furthermore, due to the importance of soil fungi for the establishment and function (nutrient uptake) of P. sylvestris, we also hypothesized that changes in soil fungal communities will have a greater impact on aboveground fungal communities than changes in aboveground fungal communities on fungi in the soil. Understanding these structural associations between above- and below-ground fungal communities will provide a relevant information about important functions linked to resilience and ecosystem functioning of P. sylvestris forests in the future. The results were also expected to provide new knowledge on overall fungal biodiversity associated with P. sylvestris that can contribute to the development of sustainable forest management strategies to effectively maintain forest biodiversity and different ecosystem services.

2. Materials and Methods

2.1. Study Sites and Sampling

The study included thirty sampling sites, each of which was 500 m2 in size, situated at least 100 m away from the forest edge, at least 15 km apart from each other, and followed the principal distribution of P. sylvestris in Lithuania (Figure 1). Information on the stand characteristics for each site is in Table 1. In each site, the health status of 30 trees was assessed by evaluating the proportion of dry branches, dechromation, and defoliation [55,56]. The sampling of needles, shoots, roots, and the rhizosphere soil was carried out between April and May 2019. At each site, five P. sylvestris trees were randomly selected, and from each tree, five branches that were ca. 17 m above the ground and growing out from the main stem were cut using a telescopic pruner. These branches were used to collone-year-old old, i.e., from the previous growing season (2018), needles and shoots. Four random needles were taken from each branch and put together, making one representative sample per tree consisting of 20 needles. For a sampling of shoots, needles were removed, and ca. 5 cm-long segment of a one-year-old shoot was randomly taken from each branch, making one representative sample per tree that consisted of five shoots. For a sampling of soil, the litter layer was removed, and three individual samples were taken down to 25 cm depth in the vicinity of five P. sylvestris trees using a 2.5 cm diameter soil core, which was carefully cleaned between individual samples. The positions of sampled soil cores were at the northern, northeastern, and southwestern sides of each tree and within a distance of 0.25 m of the tree trunk. After the collection, organic and mineral layers were separated and sieved (mesh size 2 mm × 2 mm) to remove larger particles and roots, and each layer from different samples pooled together. For soil chemical analyses, five additional sub-samples per site were taken in the vicinity of P. sylvestris trees and processed as described above. Soil chemical analyses were carried out at the Agrochemical Research Laboratory of the Lithuanian Research Centre for Agriculture and Forestry, Kaunas, Lithuania. Fine roots (lateral with root tips) were excavated within a distance of 0.5 m from the stem of five P. sylvestris trees and separated from the soil and other particles. All equipment was thoroughly cleaned between individual samples. After the collection, individual needle, shoot, root, and soil samples were placed in sterile plastic bags, labeled, placed on dry ice, the same day transported to the laboratory, and stored at −20 °C until further processing. In total, there were 150 needles, 150 shoots, 150 roots, and 180 soil samples collected. Climate data were obtained from the nearest meteorological stations.

2.2. DNA Work

Prior to DNA extraction, root samples were carefully washed in sterile water to remove any of the remaining soil and cut into ca. 5 mm-long segments. All needle, shoot, root, and soil samples were freeze-dried for 48h using Labconco FreeZone Benchtop Freeze Dryer (Cole-Parmer, Vernon Hills, IL, USA). Then, ca. 0.5 g of freeze-dried material was taken from each sample and individually ground to a fine powder using a Fast prep shaker (Montigny-le-Bretonneux, France). Approximately 30 mg of this powder per sample was used for DNA extraction. The total DNA was extracted using the CTAB method as described by Marčiulynas et al. [6]. Extracted DNA was quantified using a NanoDrop™ One spectrophotometer (Thermo Scientific, Rochester, NY, USA). PCR amplification of ITS rRNA region using barcoded primers gITS7 [57] and barcoded primers ITS4 [58] was performed according to the protocol of [59]. Within the same site, samples of the same substrate (needles, shoots, roots, and the soil) were amplified using primers with the same barcode, resulting in these samples pooling together following PCR. Individual PRRs were performed to increase the representativeness of each site. PCR amplification was performed in 50 μL reactions using an Applied Biosystems 2720 thermal cycler (Foster City, CA, USA). The PCR program started with an initial denaturation step at 94 °C for 5 min, followed by 30 cycles of 94 °C for 30 s, and annealing at 56 °C for 30 s and 72 °C for 30 s, followed by a final extension step at 72 °C for 7 min. The PCR products were assessed using gel electrophoresis on 1% agarose gels stained with GelRed (Biotium, Fremont, CA, USA). The PCR products were purified using 3 M sodium acetate (pH 5.2) (Applichem GmbH, Darmstadt, Germany) and 96% ethanol mixture (1:25). After quantification of PCR products using a Qubit fluorometer 4.0 (Life Technologies, Stockholm, Sweden), samples were pooled in an equimolar mix and used for PacBio sequencing using two SMRT cells at the SciLifeLab in Uppsala, Sweden.

2.3. Bioinformatics

Sequence quality control and clustering were performed using the SCATA NGS sequencing pipeline (http://scata.mykopat.slu.se, accessed on 14 April 2021). Quality filtering was done by removing short sequences (<200 bp), sequences with low read quality (Q < 20), primer dimers, and homopolymers, which were collapsed to 3 base pairs (bp) before clustering. Sequences that did not have a tag or primer were excluded, but information about linking the sequence to the sample was stored as metadata. The sequences were clustered into different OTUs using single-linker clustering based on 98% similarity. For each cluster, the most common genotype (real read) was used to represent each OTU. For clusters with only two sequences, a consensus sequence was created. Fungal OTUs were taxonomically identified using both the RDP classifier available at https://pyro.cme.msu.edu/index.jsp, accessed on 18 May 2021 (Centre for Microbial Ecology, Michigan State University, East Lansing, MI, USA) and the GenBank (NCBI) database using Blastn algorithm. The criteria used for identification were: sequence coverage >80%, similarity to species level 98–100%, and similarity to genus level 94–97%. Representative sequences of fungal non-singletons as the Targeted Locus Study project have been deposited in GenBank under accession number KFVY00000000. Taxonomical information was also associated with each cluster using the SH mapping feature using the UNITE database (https://unite.ut.ee/analysis.php, accessed on 26 July 2021). Fungal functional groups were assigned using the FUNGuild database (https://github.com/UMNFuN/FUNGuild, accessed on 26 July 2021) according to Nguyen et al. [60] and Tedersoo et al. [61].

2.4. Statistical Analyses

The effects of the substrate, environmental variables, and soil characteristics at a site level on OTU richness were assessed using generalized linear mixed-effect models using the glmmTMB function from the glmmTMB package [62]. Correlation between predictor variables was assessed using the cor function in R [63]. When a correlation coefficient between two variables was higher than 0.7, only one variable was selected to be used in the final model. A final model contained the following variables: tree composition, age, defoliation (stand characteristics), air temperature, precipitation (environmental variables), and soil pH, P2O5, K2O, and Ca (soil variables). Interactions between each variable and the substrate were included in the model to assess if the effects of a variable are consistent across all substrates. The site was included in the model as a random factor. All continuous variables were scaled using the function scale in R [63]. Truncated Poisson distribution for errors was assumed as all samples contained fungi. Model predictions were calculated using ggpredict from the package ggeffects [64] in R and plotted using the ggplot function from the ggplot2 package [65]. Non-metric multidimensional scaling (NMDS) based on Bray-Curtis dissimilarity was used to visualize fungal community structure among different substrates (needles, shoots, roots, and the mineral and organic soil). Prior to NMDS analysis, which was completed using a metaMDS function from the vegan package [66] in R, fungal reads were rarefied to 500 reads per sample to remove the effect of a differential sequencing depth. This eliminated 39 out of 148 samples from the data set. We used the permutational multivariate analysis of variance (PERMANOVA) with the Bray–Curtis distance metric to assess the significance of community similarity as a function of the substrate. PERMANOVA was done using the adonis2 function from the vegan package [66] in R. Pairwise comparisons between substrate levels were performed using pairwise.perm.manova function from the RVAideMemoire [67] package in R with 999 permutations. The Benjamini & Hochberg method was used for the adjustment of p-values for multiple comparisons [68]. The Shannon diversity index and qualitative Sørensen similarity index were used to characterize the diversity of fungal communities [69,70]. The nonparametric Mann–Whitney test in Minitab v. 18.1 (State College, PA, USA) was used to test if the Shannon diversity index among different samples was statistically similar or not.

3. Results

A total of 277,050 reads passing quality control were clustered into 2914 non-singleton OTUs (Table 2), while singletons were removed. Among the non-singletons, the most abundant were fungi, with 2602 (89.3%) OTUs (Table 2). Non-fungal OTU were removed from further analyses. The Shannon diversity index of fungal OTUs varied between 1.6 and 5.1 for different samples and sites, with the lowest found in root samples and the highest in soil samples (Table 2). The Shannon diversity index of fungal OTUs was significantly higher in needles than in roots, in shoots than in roots, and in the mineral and organic soil than in roots (p < 0.05) (Figure 2). In a similar comparison, other samples did not differ significantly from each other (Figure 2).
Ascomycota was the most abundant phylum in all substrates and sites, accounting for 75.4% of all sequences, followed by Basidiomycota (21.4%), Zygomycota (2.9%), Chytridiomycota (0.2%), and Glomeromycota (<0.1%) (Supplementary Table S1). Venn diagram revealed the common and unique fungal OTUs among different substrates (Figure 3). There were 109 (4.21%) OTUs that were shared among all samples. Soil harbored the highest percentage of unique OTUs (32.1%, 832 OTUs), followed by roots (8.9%, 232 OTUs), shoots (7.1%, 183), and needles (5.8%, 150 OTUs). A high percentage of unique OTUs in the soil indicates substrate-specific fungal communities. Species accumulation curves of different sample types (needles, shoots, roots, and the soil) were approaching the asymptote, showing the sequencing depth was largely sufficient and that nearly all OTUs were detected (Figure 4).
At the class level, 19 dominant classes were identified, accounting for 99.4% of all high-quality reads (Figure 5). The distribution and relative abundance of fungal classes varied among different substrates. The most dominant fungal classes in shoot samples were Eurotyomycetes (33.0%), Dothideomycetes (31.0%), and Leotiomycetes (10.9%), in needle samples—Dothideomycetes (49.4%) and Eurotyomycetes (15.1%), in roots samples—Agaricomycetes (31.6%), Leotiomycetes (27.1%), and Mucoromycotina_Incertae sedis (9.6%), while in soil samples—Agaricomycetes (21.1%), Dothideomycetes (19.5%), Leotiomycetes (11.4%), Archaeorhizomycetes (11.3%), and Eurotyomycetes (10.0%).
Identification at least to genus level was possible for 1396 (53.7%) out of 2602 fungal OTUs. Information on the 30 most common fungal OTUs representing 50.0% of all high-quality reads is in Table 3. Among these, nine OTUs representing 14.6% of all high-quality reads could not be identified to the species or genus level. The most common OTUs in all the samples were Dothideomycetes sp. 5208_5 (6.9%), Unidentified sp. 5208_1 (5.5%), Archaeorhizomyces sp. 5208_0 (4.3%), Unidentified sp. 5208_2 (3.9%) and Helotiales sp. 5208_17 (2.5%). The most common OTUs in the needle samples were Unidentified sp. 5208_5 (20.4%), Phacidium lacerum (4.5%), Unidentified sp. 5208_2 (4.4%), Unidentified sp. 5208_1 (4.3%), and Sydowia polyspora (3.4%). The most common OTUs in the shoot samples were Unidentified sp. 5208_1 (12.8%), Unidentified sp. 5208_2 (7.1%), Rhinocladiella sp. 5208_3 (6.4%), Helotiales sp. 5208_17 (6.1%). The most common OTUs in root samples were Phialocephala fortinii (11.1%), Mycena cinerella (8.9%), Archaeorhizomyces sp. 5208_0 (8.4%), and Penicillium camemberti (5.3%). The most common OTUs in the soil samples (organic and mineral soil combined) were Archaeorhizomyces sp. 5208_0 (9.0%), Unidentified sp. 5208_5 (2.6%), Umbelopsis nana (2.5%), and Malassezia restricta (2.4%).
Among the most common fungal OTUs, there were several for which the relative abundance varied significantly among the 30 sampling sites and substrates (Table 3). For example, in shoots, these were Leptosphaeria sp. 5208_64, Lecania naegelii, and Rhinocladiella sp. 5208_3. In roots, these were P. camemberti, Mucor abundans, M. hiemalis, and Armillaria ostoyae.
The assessment of fungal functional groups showed that most fungal OTUs could not be assigned to any functional group, primarily because these could not be identified to the genus or species level. Among the remaining OTUs, 7.6–10.7% per site constituted plant pathogenic fungi, 9.5–14.0%—saprotrophs, and 2.0–14.9% endophytic fungi (Figure 6). In terms of sequence reads, in different samples, fungi of the unknown functional group constituted 54.6–81.4% of reads, pathogens—6.5–12.4%, saprotrophs–6.1–11.7%, and endophytic fungi—2.0–20.4% (Figure 6). Consequently, the distribution of fungal functional groups substantially differed when compared between OTUs and sequence reads.
Generalized linear mixed effect models showed that species richness of fungal OTUs associated with roots increased with the increase of stand age while the opposite was the case for fungi in needles and organic soil (Chisq = 89, df = 4, p < 0.05). The stand age did not have an effect on the fungal species richness in shoots or the mineral soil (Figure 7A). In both organic and mineral soil, the species richness increased with the increase of defoliation, while the opposite relationship was in needles, shoots, and roots (Chisq = 133, df = 4, p < 0.05) (Figure 7B). Climatic variables also showed a significant effect on the species richness that varied depending on the substrate (Chisq = 95, df = 4, p < 0.05 and Chisq = 252, df = 4, p < 0.05 for temperature and precipitation, respectively) (Figure 7C,D). In needles, the fungal species richness increased with the increase in temperatures and precipitation. The opposite was for shoots and soil fungi, i.e., the species richness decreased with the increase of temperature and precipitation. In roots, the fungal species richness increased with the increase in temperature, but no effect of precipitation was found (Figure 7C,D).
The soil parameters significantly affected the species richness (Figure 8). The increasing concentration of P2O5 increased the fungal species richness in needles but decreased in the shoots and mineral soil (Chisq = 88, df = 4, p < 0.05) (Figure 8A). For the roots and organic soil, the concentration of P2O5 had no significant effect (Figure 8A). The increasing concentration of K2O increased the fungal species richness in shoots but decreased in other substrates, i.e., needles, roots, minerals, and organic soil (Chisq = 24, df = 4, p < 0.05) (Figure 8B). Soil pH had a similar effect on fungal species richness in the needles, minerals, and organic soil, i.e., species richness decreased with the increase in pH (Chisq = 127, df = 4, p < 0.05). The opposite was in the shoots, while the fungal species richness in roots was generally unaffected (Figure 8C). In all substrates, the fungal species richness increased with the increase of Ca concentration (Chisq = 23, df = 4, p < 0.05) (Figure 8D).
The type of substrate (needles, shoots, roots, or soil) had a strong effect on the composition of the fungal communities as the fungal communities associated with different substrates differed significantly from each other (PERMANOVA: df = 4, SumOfSqs = 12.554, R2 = 0.32793, F = 12.687, p < 0.001) (Table 4). An exception was the fungal communities in organic and mineral soil, which were similar (p > 0.05). Similarly, differences in the fungal communities among different substrates were also indicated by the NMDS of Bray–Curtis dissimilarity based on the OTU-level abundance (stress = 0.15) (Figure 9). PERMANOVA has shown that the fungal community composition varied among study sites, which explained almost as much variation in the fungal community composition as it was explained by the substrate (df = 29, SumOfSqs = 8.055, R2 = 0.210, F = 1.179, p < 0.01).

4. Discussion

As climate change is expected to affect the natural distribution of forest tree species, entire forest ecosystems are likely to be subjected to major changes to adapt to new conditions. However, it is less clear whether and how global fungal biodiversity will respond to changes in the distribution range of the host trees. In the present study, to get a deeper understanding of these questions, the impact of abiotic factors on the above- and below-ground diversity and composition of fungal communities were studied in 30 P. sylvestris sites. Using the high-throughput sequencing data, different environmental variables, and soil chemical properties, we identified key factors that drive the diversity and composition of fungal communities associated with P. sylvestris. Moreover, by linking the fungal community data and different substrates, we have demonstrated that the fungal species richness and community composition depend on the substrate they colonize. Such patterns of fungal specificity for the substrate have already been shown in Picea abies stands [71]. Besides, the fungal species richness can be different in the same tissues depending on their life stage, e.g., living or dead wood [49]. In addition, we have shown that the fungal species richness is highly dependent on the site. The site-specific composition of fungal communities was shown previously for soil, root, or foliage samples [72,73,74,75]. Overall, our results demonstrated that different functional tissues and the rhizosphere soil of P. sylvestris are associated with diverse fungal communities driven by a combination of climatic and edaphic factors, thereby providing insights into possible ecological responses of fungal communities to climate change in northern Europe.
Climate is one of the main factors influencing fungal development, either directly or indirectly, by triggering tree responses [12]. However, the diversity of fungal communities can be influenced by a variety of factors, including tree species composition, stand age, habitat conditions, and edaphic factors [76]. Despite many studies on interactions between soil fungi and their hosts, the interactions between fungi and above- and below-ground tree traits are less studied. It is emphasized that when assessing the interaction of above- and below-ground communities, it is important to identify the dynamics of soil microbial communities that reflect changes in forest age and soil properties [77]. Our findings suggest that the fungal species richness in needles, organic soil, and roots of P. sylvestris change with the stand age. In agreement with other studies [78], the greatest positive effect of stand age was on fungal species richness in roots. Such results could be due to the accumulation of fungal species over time as well as due to changes in nutrient supply by host trees [71,78]. Several studies have also reported a strong effect of tree age on fungal communities, as the litter quality and understory vegetation in young stands differ markedly from older stands [79]. It was shown that fungal communities in roots might be specifically shaped at the fine scale, but this may disappear when averaged across an entire landscape, obscuring the specific environmental conditions and species interactions that drive fungal diversity [80]. Indeed, the assembly of root-associated fungi has been shown to be regulated by a wide range of spatial and temporal variables [81].
Fungal communities are dynamic components of terrestrial ecosystems and exhibit temporal and spatial variation [77]. Across the range of a single host tree species, the fungal community may change in response to climatic factors [82]. It is also known that moisture and temperature influence the growth of fungi and that characteristic weather conditions favorable for fungi may be used to predict their abundance and richness in habitats with different climatic conditions [83]. Although climatic factors such as temperature and precipitation were found to promote fungal richness globally [84], our results suggest that the richness of fungal OTUs in the organic soil layer decreases with the increase in temperature and precipitation. Temperature is one of the main factors determining soil microorganisms’ activity [85]. It is known that temperature enhances microbial metabolic activity, potentially leading to an accelerated litter decomposition rate [86]. Low temperature is considered to be one of the dominant forces protecting soil C from decomposition [87]. However, it seems likely that soil organic C will decrease with increasing temperature due to climate change [88]. With an increase in temperature, evaporation is enhanced, and plant productivity is reduced, ultimately resulting in a decrease in soil organic carbon input [89]. Finally, shifts in the amount of above- and below-ground organic matter inputs to the soil may also shape the composition and activity of fungal communities [89,90]. Studies by Dang et al. [77] have also shown that the abundance of dominant fungal communities is significantly correlated with organic C, total N, C: N, available N, and available P, indicating the dependence of these microbes on soil nutrients. Therefore, it could be assumed that as the temperature increases and the C content decreases, the abundance of microbial communities is prone to decrease. Yet, large knowledge gaps remain that feed uncertainty around the temperature and sensitivity of soil microbial processes [91].
Recent studies have shown that the decrease in the abundance of saprotrophic fungi may have been a direct result of forest disturbance or an indirect result of changes in soil pH and soil P [92]. Furthermore, the dominant P-solubilizing saprotrophic fungi are replaced by diverse facultative pathogenic fungi with weaker C decomposition ability. These changes potentially indicate a shift from soil phosphate limitation to carbon limitation following deforestation [92].
Some belowground properties, including plant richness and plant diversity, total carbon, total nitrogen, soil pH, and nutrient content, change simultaneously with changes in the aboveground structure [77]. Furthermore, these indices eventually influence the function and structure of the soil microbial community [77]. Our study has shown that P. sylvestris fungal communities in the soil, roots, and shoots negatively respond to the increase in soil pH. The effects of soil pH on the soil microbial community assembly occur through the availability of nutrients and carbon and the solubility of metals, which are strongly influenced by changes in soil pH concentration [19]. Previous studies also confirm that soils with lower pH concentrations increase the rate of fungal growth in contrast to soils with higher pH [93]. Decreased fungal species richness in forest soils with increasing pH is also associated with decreased availability of carbon [93,94]. The increase in soil pH also reduces the rhizosphere priming effect, which leads to an increase in the accumulation of carbon in the soil, which may also be one of the reasons for the decrease in fungal richness in the soil [94]. It is also important to mention that soil pH is closely related to vegetation and strongly influences the degree of ectomycorrhizal colonization. In particular, it has been observed that a decrease in pH results in a shift in vegetation, leading to an increase in ectomycorrhizal fungi [95].
In the present study, the most common fungal phyla among different P. sylvestris substrates were Ascomycota and Basidiomycota, which are known to dominate among fungi colonizing soils and terrestrial plant tissues [7,11,96,97]. Interestingly, P. sylvestris needles showed a higher relative abundance of fungal OTUs belonging to Helotiales than Dothideales, confirming the hypothesis that Helotiales are dominated in gymnosperms and Dothideales in angiosperms [47]. The major differences between the above- and below-ground communities were the dominance of the Dothideomycetes in shoots and needles, while Archaeorhizomyces dominated in roots and soil samples. This is not surprising as representatives of both Dothideomycetes and Archaeorhizomyces were among the most commonly found fungal classes in pine forests [98,99,100,101].
Fungi colonizing different P. sylvestris tissues accounted for only a fraction of all fungal OTUs. As a result, communities of soil fungi had higher diversity than those inhabiting the aboveground parts, and in the organic soil, the abundance of fungi was higher than in the mineral soil. Numerous studies indicate that soil is the primary source of microbiomes in terrestrial habitats [102]. Since resource availability is a key factor regulating biodiversity and ecosystem functioning [103], it could be that soil harbors a greater diversity of microbes by providing more enriched or easily accessible resources and being an ecotone for both above- and below-ground communities. It has previously been shown that plant species can modify the soil environment and support rhizosphere microbes, which in turn can provide feedback that promotes plant health and growth [104]. The study by Bulgarelli et al. [105] on Pinus ponderosa sampled from several locations around the globe revealed minimal geographic differentiation of phyllosphere bacterial communities, which supports the concept that the host plant species is a determinant for the structure of the phyllosphere community. In comparison to other tissues and the soil, P. sylvestris needles were more often colonized by plant pathogens. It is important to emphasize that foliar pathogens, while infecting and causing the disease to the needles or leaves, can also, especially at high infection levels, cause growth reduction, increase susceptibility to biotic and abiotic stresses, and in severe cases, cause tree mortality [106,107]. This may also suggest that fungal communities in needles and roots, which are the primary functional tissues of trees, respond faster to biotic and abiotic stresses, while in the soil and shoots, these changes are slower [108].
Among soil and root-associated fungi, ectomycorrhizal fungi (ECM) represent a key component needed for the successful establishment and growth of P. sylvestris trees. Among these, Suilloid fungi are of key importance for pines, especially those of the Suillus and Rhizopogon genera [109]. Despite the key importance of these ECM fungi, they occupied only 0.5% of the total diversity in the soil and 0.1% in the roots of investigated samples. Suilloid fungi represent species of early succession and, with the age of the stand, are often replaced by other ECM fungi, which leads to the fact that they become less common [110,111]. Besides, a higher proportion of this group of fungi is usually found after forest disturbance [109]. For example, Rhizopogon species were found to dominate ECM communities associated with pine roots in post-fire sites [112]. Among other more frequently detected fungi was Umbelopsis nana, which accounted for 2.7% of the total fungal diversity. Members of the genus Umbelopsis are common soil fungi that are found in forest ecosystems around the world [113,114]. This species appears to be an indicator of the mineral soil habitats in many northern forests [115]. Studies in Swedish clear-cut P. sylvestris forests also revealed that U. nana was one of the most common species of soil fungi [116]. In previous studies, members of this genus were also detected in the roots of P. sylvestris [101], P. abies [117], Picea mariana [115], Pseudotsuga menziesii, Pinus ponderosa [114], and Quercus sp. [118]. Umbelopsis fungi are also known to synthesize polyunsaturated fatty acids, which are important in plant stress defense [119] and can be important in the process of climate change due to the increase of stressors. Phacidium lacerum, which was frequently detected in the soil and shoot samples, is widely distributed throughout Europe, where it commonly occurs on P. sylvestris [120], but it is not known to cause the disease on this host [121]. However, it was recently reported to be the causal agent of postharvest rot in apples and pears [122]. Therefore, the pathogenic nature of pines should not be excluded. Sydowia polyspora, often found in needles and shoots, is commonly associated with conifers worldwide and is considered a pathogen on several hosts [123]. Moreover, Cleary et al. [124] recently found S. polyspora in asymptomatic seeds of several Pinus spp. obtained across Europe and North America, including P. pinaster, P. radiata, P. strobus, P. sylvestris, P. mugo, and P. pinea. Fungus is also vectored by insects [125] and very often reported as an endophyte or a saprophyte [126]. The endophytic way of life of pathogenic fungi can be a starting point at an early stage of the infection process before transitioning to a pathogenic lifestyle [127].

5. Conclusions

Our study has shown that P. sylvestris forests in Lithuania are associated with a high diversity of above- and below-ground fungal communities, which are strongly influenced by interrelated environmental factors and stand characteristics. Patterns of fungal community composition were strongly associated with particular substrates, showing that they provide different ecological niches that are preferred by different fungal species. To maintain the availability and suitability of ecological niches for different fungal species, it is essential to maintain an appropriate soil pH and the availability of soil nutrients. Disturbances related to climate change are likely to affect these important factors and, thus, may indirectly affect the structure and diversity of fungal communities. On the practical side, forest management and planning in the future will have to adapt to changing climatic conditions to ensure ecosystem services through the introduction of new management approaches. This is particularly important for maintaining soil health. In the future, the stands with other, possibly non-native, tree species, which better tolerate climate change, may also become unavoidable to maintain the diversity of microorganisms and ensure the growth and sustainability of forest stands.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jof8101112/s1, Table S1: Relative abundance (%) of fungal taxa associated with needles, shoots, roots, and the rhizosphere soil of Pinus sylvestris in Lithuania.

Author Contributions

Conceptualization, A.M. (Audrius Menkis) and D.M.; methodology, D.M., A.M. (Adas Marčiulynas) and A.M. (Audrius Menkis); validation, A.M. (Audrius Menkis) and D.M.; formal analysis, D.M., I.F., A.M. (Adas Marčiulynas) and V.M.; investigation, D.M., A.M. (Adas Marčiulynas), A.G., J.L. and V.M.; data curation, D.M. and A.M. (Audrius Menkis); writing—original draft preparation, D.M. and A.M. (Audrius Menkis); writing—review and editing, D.M., A.M. (Audrius Menkis), A.M. (Adas Marčiulynas), J.L., V.M. and I.F.; supervision, A.M. (Audrius Menkis); project administration, D.M.; funding acquisition, A.M. (Audrius Menkis). All authors have read and agreed to the published version of the manuscript.

Funding

This project has received funding from the European Social Fund (project No. 09.3.3-LMT-K-712-01-0039) under a grant agreement with the Research Council of Lithuania (LMTLT).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank Dangutė Šalevičienė for her help with sample collection and preparation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Oleksyn, J.; Wyka, T.P.; Żytkowiak, R.; Zadworny, M.; Mucha, J.; Dering, M.; Ufnalski, K.; Nihlgård, B.; Reich, P.B. A fingerprint of climate change across pine forests of Sweden. Ecol. Lett. 2020, 23, 1739–1746. [Google Scholar] [CrossRef]
  2. Cavicchioli, R.; Ripple, W.J.; Timmis, K.N.; Azam, F.; Bakken, L.R.; Baylis, M.; Behrenfeld, M.J.; Boetius, A.; Boyd, P.W.; Classen, A.T.; et al. Scientists’ warning to humanity: Microorganisms and climate change. Nat. Rev. Microbiol. 2019, 17, 569–586. [Google Scholar] [CrossRef]
  3. Singh, D.P.; Prabha, R. Microbial Interventions in Agriculture and Environment: Volume 3: Soil and Crop Health Management; Springer: Berlin/Heidelberg, Germany, 2019; pp. 1–491. [Google Scholar] [CrossRef]
  4. de Habiyaremye, J.D.; Goldmann, K.; Reitz, T.; Herrmann, S.; Buscot, F. Tree Root Zone Microbiome: Exploring the Magnitude of Environmental Conditions and Host Tree Impact. Front. Microbiol. 2020, 11, 749. [Google Scholar] [CrossRef]
  5. Barnes, E.M.; Tringe, S.G. Exploring the roles of microbes in facilitating plant adaptation to climate change. Biochem. J. 2022, 479, 327–335. [Google Scholar] [CrossRef]
  6. Hardoim, P.R.; van Overbeek, L.S.; Berg, G.; Pirttilä, A.M.; Compant, S.; Campisano, A.; Döring, M.; Sessitsch, A. The Hidden World within Plants: Ecological and Evolutionary Considerations for Defining Functioning of Microbial Endophytes. Microbiol. Mol. Biol. Rev. 2015, 79, 293–320. [Google Scholar] [CrossRef]
  7. Terhonen, E.; Blumenstein, K.; Kovalchuk, A.; Asiegbu, F.O. Forest tree microbiomes and associated fungal endophytes: Functional roles and impact on forest health. Forests 2019, 10, 42. [Google Scholar] [CrossRef]
  8. Kasanen, R. Abiotic factors affecting the composition of forest tree microbiomes. In Forest Microbiology; Academic Press: Cambridge, MA, USA, 2021; pp. 59–73. [Google Scholar]
  9. Wu, Y.; Qu, M.; Pu, X.; Lin, J.; Shu, B. Distinct microbial communities among different tissues of citrus tree Citrus reticulata cv. Chachiensis. Sci. Rep. 2020, 10, 6068. [Google Scholar] [CrossRef]
  10. Turner, T.R.; James, E.K.; Poole, P.S. The plant microbiome. Genome Biol. 2013, 14, 209. [Google Scholar] [CrossRef]
  11. Trivedi, P.; Leach, J.E.; Tringe, S.G.; Sa, T.; Singh, B.K. Plant–microbiome interactions: From community assembly to plant health. Nat. Rev. Microbiol. 2020, 18, 607–621. [Google Scholar] [CrossRef]
  12. Faticov, M.; Abdelfattah, A.; Roslin, T.; Vacher, C.; Hambäck, P.; Blanchet, F.G.; Lindahl, B.D.; Tack, A.J.M. Climate warming dominates over plant genotype in shaping the seasonal trajectory of foliar fungal communities on oak. New Phytol. 2021, 231, 1770–1783. [Google Scholar] [CrossRef]
  13. Lau, J.A.; Lennon, J.T.; Heath, K.D. Trees harness the power of microbes to survive climate change. Proc. Natl. Acad. Sci. USA 2017, 114, 11009–11011. [Google Scholar] [CrossRef]
  14. Bennett, A.E.; Classen, A.T. Climate change influences mycorrhizal fungal–plant interactions, but conclusions are limited by geographical study bias. Ecology 2020, 101, e02978. [Google Scholar] [CrossRef]
  15. Bender, S.F.; Wagg, C.; van der Heijden, M.G.A. An Underground Revolution: Biodiversity and Soil Ecological Engineering for Agricultural Sustainability. Trends Ecol. Evol. 2016, 31, 440–452. [Google Scholar] [CrossRef]
  16. Frac, M.; Hannula, S.E.; Belka, M.; Jȩdryczka, M. Fungal biodiversity and their role in soil health. Front. Microbiol. 2018, 9, 707. [Google Scholar] [CrossRef]
  17. Mercado-Blanco, J.; Abrantes, I.; Caracciolo, A.B.; Bevivino, A.; Ciancio, A.; Grenni, P.; Hrynkiewicz, K.; Kredics, L.; Proença, D.N. Belowground microbiota and the health of tree crops. Front. Microbiol. 2018, 9, 1006. [Google Scholar] [CrossRef]
  18. Zhou, Z.; Wang, C.; Luo, Y. Effects of forest degradation on microbial communities and soil carbon cycling: A global meta-analysis. Glob. Ecol. Biogeogr. 2018, 27, 110–124. [Google Scholar] [CrossRef]
  19. Rousk, J.; Brookes, P.C.; Bååth, E. Contrasting soil pH effects on fungal and bacterial growth suggest functional redundancy in carbon mineralization. Appl. Environ. Microbiol. 2009, 75, 1589–1596. [Google Scholar] [CrossRef]
  20. van der Heijden, M.G.A.; Martin, F.M.; Selosse, M.A.; Sanders, I.R. Mycorrhizal ecology and evolution: The past, the present, and the future. New Phytol. 2015, 205, 1406–1423. [Google Scholar] [CrossRef]
  21. Saravesi, K.; Aikio, S.; Wäli, P.R.; Ruotsalainen, A.L.; Kaukonen, M.; Huusko, K.; Suokas, M.; Brown, S.P.; Jumpponen, A.; Tuomi, J.; et al. Moth Outbreaks Alter Root-Associated Fungal Communities in Subarctic Mountain Birch Forests. Microb. Ecol. 2015, 69, 788–797. [Google Scholar] [CrossRef]
  22. Hartmann, M.; Niklaus, P.A.; Zimmermann, S.; Schmutz, S.; Kremer, J.; Abarenkov, K.; Lüscher, P.; Widmer, F.; Frey, B. Resistance and resilience of the forest soil microbiome to logging-associated compaction. ISME J. 2014, 8, 226–244. [Google Scholar] [CrossRef]
  23. Treu, R.; Karst, J.; Randall, M.; Pec, G.J.; Cigan, P.W.; Simard, S.W.; Cooke, J.E.K.; Erbilgin, N.; Cahill, J.F. Decline of ectomycorrhizal fungi following a mountain pine beetle epidemic. Ecology 2014, 95, 1096–1103. [Google Scholar] [CrossRef] [PubMed]
  24. Karst, J.; Erbilgin, N.; Pec, G.J.; Cigan, P.W.; Najar, A.; Simard, S.W.; Cahill, J.F. Ectomycorrhizal fungi mediate indirect effects of a bark beetle outbreak on secondary chemistry and establishment of pine seedlings. New Phytol. 2015, 208, 904–914. [Google Scholar] [CrossRef] [PubMed]
  25. Bödeker, I.T.M.; Lindahl, B.D.; Olson, Å.; Clemmensen, K.E. Mycorrhizal and saprotrophic fungal guilds compete for the same organic substrates but affect decomposition differently. Funct. Ecol. 2016, 30, 1967–1978. [Google Scholar] [CrossRef]
  26. Crowther, T.W.; van den Hoogen, J.; Wan, J.; Mayes, M.A.; Keiser, A.D.; Mo, L.; Averill, C.; Maynard, D.S. The global soil community and its influence on biogeochemistry. Science 2019, 365, eaav0550. [Google Scholar] [CrossRef]
  27. Usman, M.; Ho-Plágaro, T.; Frank, H.E.R.; Calvo-Polanco, M.; Gaillard, I.; Garcia, K.; Zimmermann, S.D. Mycorrhizal Symbiosis for Better Adaptation of Trees to Abiotic Stress Caused by Climate Change in Temperate and Boreal Forests. Front. For. Glob. Chang. 2021, 4, 141. [Google Scholar] [CrossRef]
  28. Steidinger, B.S.; Crowther, T.W.; Liang, J.; Van Nuland, M.E.; Werner, G.D.A.; Reich, P.B.; Nabuurs, G.; de-Miguel, S.; Zhou, M.; Picard, N.; et al. Climatic controls of decomposition drive the global biogeography of forest-tree symbioses. Nature 2019, 569, 404–408. [Google Scholar] [CrossRef]
  29. Park, K.H.; Oh, S.Y.; Yoo, S.; Park, M.S.; Fong, J.J.; Lim, Y.W. Successional Change of the Fungal Microbiome Pine Seedling Roots Inoculated With Tricholoma matsutake. Front. Microbiol. 2020, 11, 574146. [Google Scholar] [CrossRef]
  30. Lebeis, S.L. Greater than the sum of their parts: Characterizing plant microbiomes at the community-level. Curr. Opin. Plant Biol. 2015, 24, 82–86. [Google Scholar] [CrossRef]
  31. Wallander, H.; Johansson, U.; Sterkenburg, E.; Brandström Durling, M.; Lindahl, B.D. Production of ectomycorrhizal mycelium peaks during canopy closure in Norway spruce forests. New Phytol. 2010, 187, 1124–1134. [Google Scholar] [CrossRef]
  32. Bulgarelli, D.; Rott, M.; Schlaeppi, K.; Ver Loren van Themaat, E.; Ahmadinejad, N.; Assenza, F.; Rauf, P.; Huettel, B.; Reinhardt, R.; Schmelzer, E.; et al. Revealing structure and assembly cues for Arabidopsis root-inhabiting bacterial microbiota. Nature 2012, 488, 91–95. [Google Scholar] [CrossRef]
  33. Lundberg, D.S.; Lebeis, S.L.; Paredes, S.H.; Yourstone, S.; Gehring, J.; Malfatti, S.; Tremblay, J.; Engelbrektson, A.; Kunin, V.; del Rio, T.G.; et al. Defining the core Arabidopsis thaliana root microbiome. Nature 2012, 488, 86–90. [Google Scholar] [CrossRef] [PubMed]
  34. Kohout, P.; Sudová, R.; Brabcová, V.; Vosolsobě, S.; Baldrian, P.; Albrechtová, J. Forest Microhabitat Affects Succession of Fungal Communities on Decomposing Fine Tree Roots. Front. Microbiol. 2021, 12, 541583. [Google Scholar] [CrossRef] [PubMed]
  35. Lladó, S.; López-Mondéjar, R.; Baldrian, P. Drivers of microbial community structure in forest soils. Appl. Microbiol. Biotechnol. 2018, 102, 4331–4338. [Google Scholar] [CrossRef] [PubMed]
  36. Lareen, A.; Burton, F.; Schäfer, P. Plant root-microbe communication in shaping root microbiomes. Plant Mol. Biol. 2016, 90, 575–587. [Google Scholar] [CrossRef]
  37. 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]
  38. Mäkipää, R.; Rajala, T.; Schigel, D.; Rinne, K.T.; Pennanen, T.; Abrego, N.; Ovaskainen, O. Interactions between soil- and dead wood-inhabiting fungal communities during the decay of Norway spruce logs. ISME J. 2017, 11, 1964–1974. [Google Scholar] [CrossRef]
  39. Bagchi, R.; Gallery, R.E.; Gripenberg, S.; Gurr, S.J.; Narayan, L.; Addis, C.E.; Freckleton, R.P.; Lewis, O.T. Pathogens and insect herbivores drive rainforest plant diversity and composition. Nature 2014, 506, 85–88. [Google Scholar] [CrossRef]
  40. Millberg, H.; Boberg, J.; Stenlid, J. Changes in fungal community of Scots pine (Pinus sylvestris) needles along a latitudinal gradient in Sweden. Fungal Ecol. 2015, 17, 126–139. [Google Scholar] [CrossRef]
  41. Hyde, K.D.; Soytong, K. Understanding microfungal diversity: A critique. Cryptogam. Mycol. 2007, 28, 281–289. [Google Scholar]
  42. Osorio, M.; Stephan, B.R. Life cycle of Lophodermium piceae in Norway spruce needles. Eur. J. For. Pathol. 1991, 21, 152–163. [Google Scholar] [CrossRef]
  43. Larkin, B.G.; Hunt, L.S.; Ramsey, P.W. Foliar nutrients shape fungal endophyte communities in Western white pine (Pinus monticola) with implications for white-tailed deer herbivory. Fungal Ecol. 2012, 5, 252–260. [Google Scholar] [CrossRef]
  44. Zimmerman, N.B.; Vitousek, P.M. Fungal endophyte communities reflect environmental structuring across a Hawaiian landscape. Proc. Natl. Acad. Sci. USA 2012, 109, 13022–13027. [Google Scholar] [CrossRef] [PubMed]
  45. Romeralo, C.; Martín-García, J.; Martínez-Álvarez, P.; Muñoz-Adalia, E.J.; Gonçalves, D.R.; Torres, E.; Witzell, J.; Diez, J.J. Pine species determine fungal microbiome composition in a common garden experiment. Fungal Ecol. 2022, 56, 101137. [Google Scholar] [CrossRef]
  46. Rabiey, M.; Hailey, L.E.; Roy, S.R.; Grenz, K.; Al-Zadjali, M.A.S.; Barrett, G.A.; Jackson, R.W. Endophytes vs. tree pathogens and pests: Can they be used as biological control agents to improve tree health? Eur. J. Plant Pathol. 2019, 155, 711–729. [Google Scholar] [CrossRef]
  47. Sieber, T.N. Endophytic fungi in forest trees: Are they mutualists? Fungal Biol. Rev. 2007, 21, 75–89. [Google Scholar] [CrossRef]
  48. Mengistu, A.A. Endophytes: Colonization, Behaviour, and Their Role in Defense Mechanism. Int. J. Microbiol. 2020, 2020, 6927219. [Google Scholar] [CrossRef]
  49. Rajala, T.; Peltoniemi, M.; Pennanen, T.; Mäkipää, R. Fungal community dynamics in relation to substrate quality of decaying Norway spruce (Picea abies [L.] Karst.) logs in boreal forests. FEMS Microbiol. Ecol. 2012, 81, 494–505. [Google Scholar] [CrossRef]
  50. Purahong, W.; Wubet, T.; Krüger, D.; Buscot, F. Molecular evidence strongly supports deadwood-inhabiting fungi exhibiting unexpected tree species preferences in temperate forests. ISME J. 2018, 12, 289–295. [Google Scholar] [CrossRef]
  51. Olsson, J. Colonization Patterns of Wood-Inhabiting Fungi in Boreal Forest; Ekologi, Miljö Och Geovetenskap: Umeå, Sweden, 2008; p. 35. [Google Scholar]
  52. Bahram, M.; Netherway, T. Fungi as mediators linking organisms and ecosystems. FEMS Microbiol. Rev. 2022, 46, fuab058. [Google Scholar] [CrossRef]
  53. Hallingbäck, H.R.; Burton, V.; Vizcaíno-Palomar, N.; Trotter, F.; Liziniewicz, M.; Marchi, M.; Berlin, M.; Ray, D.; Benito Garzón, M. Managing Uncertainty in Scots Pine Range-Wide Adaptation Under Climate Change. Front. Ecol. Evol. 2021, 9, 724051. [Google Scholar] [CrossRef]
  54. Reich, P.B.; Oleksyn, J. Climate warming will reduce growth and survival of Scots pine except in the far north. Ecol. Lett. 2008, 11, 588–597. [Google Scholar] [CrossRef] [PubMed]
  55. Vorontsov, A.I.; Mozolevskaya, E.G.; Sokolova, E.S. Tekhnologiya Zashchity Lesa—M; Ekologiya: Moscow, Russia, 1991. [Google Scholar]
  56. Ozolincius, R.; Stakenas, V. Tree crown defoliation: Influencing factors. Balt. For. 1996, 2, 48–55. [Google Scholar]
  57. Ihrmark, K.; Bödeker, I.T.M.; Cruz-Martinez, K.; Friberg, H.; Kubartova, A.; Schenck, J.; Strid, Y.; Stenlid, J.; Brandström-Durling, M.; Clemmensen, K.E.; et al. New primers to amplify the fungal ITS2 region—Evaluation by 454-sequencing of artificial and natural communities. FEMS Microbiol. Ecol. 2012, 82, 666–677. [Google Scholar] [CrossRef] [PubMed]
  58. White, T.J.; Bruns, T.; Lee, S.; Taylor, J. Amplification and Direct Sequencing of Fungal Ribosomal Rna Genes for Phylogenetics. PCR Protoc. 1990, 18, 315–322. [Google Scholar] [CrossRef]
  59. Clemmensen, K.E.; Ihrmark, K.; Brandström-Durling, M.; Lindahl, B.D. Sample preparation for fungal community analysis by high-throughput sequencing of barcode amplicons. In Microbial Environmental Genomics (MEG); Methods in Molecular Biology; Martin, F., Uroz, S., Eds.; Humana Press: New York, NY, USA, 2016; Volume 1399, pp. 61–88. [Google Scholar] [CrossRef]
  60. 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]
  61. 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]
  62. Brooks, M.E.; Kristensen, K.; van Benthem, K.J.; Magnusson, A.; Berg, C.W.; Nielsen, A.; Skaug, H.J.; Mächler, M.; Bolker, B.M. glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. R J. 2017, 9, 378–400. [Google Scholar] [CrossRef]
  63. R: A Language and Environment for Statistical Computing. In R Foundation for Statistical Computing; R Core Team: Vienna, Austria, 2020.
  64. Lüdecke, D. ggeffects: Tidy Data Frames of Marginal Effects from Regression Models. J. Open Source Softw. 2018, 3, 772. [Google Scholar] [CrossRef]
  65. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016; ISBN 978-3-319-24277-4. [Google Scholar]
  66. Oksanen, J.; Blanchet, F.G.; Kindt, R.; Legendre, P.; Minchin, P.; O’Hara, B.; Simpson, G.; Solymos, P.; Stevens, H.; Wagner, H. Vegan: Community Ecology Package, R Package Version 2.5–3. 2018. Available online: https://cran.r-project.org (accessed on 28 November 2021).
  67. Hervé, M. RVAideMemoire: Testing and Plotting Procedures for Biostatistics, R Packag. Version 0.9-81-2. 2022. Available online: https://cran.r-project.org/web/packages/RVAideMemoire/RVAideMemoire.pdf (accessed on 15 June 2022).
  68. Haynes, W. Benjamini–Hochberg Method BT—Encyclopedia of Systems Biology; Dubitzky, W., Wolkenhauer, O., Cho, K.-H., Yokota, H., Eds.; Springer: New York, NY, USA, 2013; p. 78. ISBN 978-1-4419-9863-7. [Google Scholar]
  69. Shannon, C.E. A Mathematical Theory of Communication. Bell Syst. Technol. J. 1948, 27, 379–423. [Google Scholar] [CrossRef]
  70. Magurran, A.E. Ecological Diversity and Its Measurement. J. Appl. Ecol. 1989, 26, 1101. [Google Scholar] [CrossRef]
  71. Marčiulynas, A.; Marčiulynienė, D.; Mishcherikova, V.; Franić, I.; Lynikienė, J.; Gedminas, A.; Menkis, A. High Variability of Fungal Communities Associated with the Functional Tissues and Rhizosphere Soil of Picea abies in the Southern Baltics. Forests 2022, 13, 1103. [Google Scholar] [CrossRef]
  72. Hoeber, S.; Baum, C.; Weih, M.; Manzoni, S.; Fransson, P. Site-Dependent Relationships between Fungal Community Composition, Plant Genotypic Diversity and Environmental Drivers in a Salix Biomass System. Front. Fungal Biol. 2021, 2, 31. [Google Scholar] [CrossRef]
  73. Alzarhani, A.K.; Clark, D.R.; Underwood, G.J.C.; Ford, H.; Cotton, T.E.A.; Dumbrell, A.J. Are drivers of root-associated fungal community structure context specific? ISME J. 2019, 13, 1330–1344. [Google Scholar] [CrossRef] [PubMed]
  74. Adamo, I.; Ortiz-Malavasi, E.; Chazdon, R.; Chaverri, P.; Ter Steege, H.; Geml, J. Soil fungal community composition correlates with site-specific abiotic factors, tree community structure, and forest age in regenerating tropical rainforests. Biology 2021, 10, 1120. [Google Scholar] [CrossRef]
  75. Lumibao, C.Y.; Borer, E.T.; Condon, B.; Kinkel, L.; May, G.; Seabloom, E.W. Site-specific responses of foliar fungal microbiomes to nutrient addition and herbivory at different spatial scales. Ecol. Evol. 2019, 9, 12231–12244. [Google Scholar] [CrossRef]
  76. Smith, J.E.; Molina, R.; Huso, M.M.P.; Luoma, D.L.; McKay, D.; Castellano, M.A.; Lebel, T.; Valachovic, Y. Species richness, abundance, and composition of hypogeous and epigeous ectomycorrhizal fungal sporocarps in young, rotation-age, and old-growth stands of Douglas-fir (Pseudotsuga menziesii) in the Cascade Range of Oregon, U.S.A. Can. J. Bot. 2002, 80, 186–204. [Google Scholar] [CrossRef]
  77. Dang, P.; Yu, X.; Le, H.; Liu, J.; Shen, Z.; Zhao, Z. Effects of stand age and soil properties on soil bacterial and fungal community composition in Chinese pine plantations on the Loess Plateau. PLoS ONE 2017, 12, e0186501. [Google Scholar] [CrossRef]
  78. Odriozola, I.; Martinovic, T.; Bahnmann, B.D.; Ryšánek, D.; Mašínová, T.; Sedlák, P.; Merunková, K.; Kohout, P.; Tomšovský, M.; Baldrian, P. Stand age affects fungal community composition in a Central European temperate forest. Fungal Ecol. 2020, 48, 100985. [Google Scholar] [CrossRef]
  79. Unuk, T.; Martinović, T.; Finžgar, D.; Šibanc, N.; Grebenc, T.; Kraigher, H. Root-Associated Fungal Communities from Two Phenologically Contrasting Silver Fir (Abies alba Mill.) Groups of Trees. Front. Plant Sci. 2019, 10, 214. [Google Scholar] [CrossRef]
  80. Izzo, A.; Agbowo, J.; Bruns, T.D. Detection of plot-level changes in ectomycorrhizal communities across years in an old-growth mixed-conifer forest. New Phytol. 2005, 166, 619–630. [Google Scholar] [CrossRef]
  81. Barnes, C.J.; van der Gast, C.J.; McNamara, N.P.; Rowe, R.; Bending, G.D. Extreme rainfall affects assembly of the root-associated fungal community. New Phytol. 2018, 220, 1172–1184. [Google Scholar] [CrossRef] [PubMed]
  82. Bui, A.; Orr, D.; Lepori-Bui, M.; Konicek, K.; Young, H.S.; Moeller, H.V. Soil fungal community composition and functional similarity shift across distinct climatic conditions. FEMS Microbiol. Ecol. 2020, 96, fiaa193. [Google Scholar] [CrossRef] [PubMed]
  83. Talley, S.M.; Coley, P.D.; Kursar, T.A. The effects of weather on fungal abundance and richness among 25 communities in the Intermountain West. BMC Ecol. 2002, 2, 7–10. [Google Scholar] [CrossRef] [PubMed]
  84. Swaty, R.L.; Gehring, C.A.; Van Ert, M.; Theimer, T.A.D.C.; Keim, P.; Whitham, T.G. Temporal variation in temperature and rainfall differentially affects ectomycorrhizal colonization at two contrasting sites. New Phytol. 1998, 139, 733–739. [Google Scholar] [CrossRef]
  85. Järvan, M.; Edesi, L.; Adamson, A.; Võsa, T. Soil microbial communities and dehydrogenase activity depending on farming systems. Plant Soil Environ. 2014, 60, 459–463. [Google Scholar] [CrossRef]
  86. Schindlbacher, A.; Rodler, A.; Kuffner, M.; Kitzler, B.; Sessitsch, A.; Zechmeister-Boltenstern, S. Experimental warming effects on the microbial community of a temperate mountain forest soil. Soil Biol. Biochem. 2011, 43, 1417–1425. [Google Scholar] [CrossRef]
  87. Davidson, E.A.; Janssens, I.A. Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature 2006, 440, 165–173. [Google Scholar] [CrossRef]
  88. Kirschbaum, M.U.F. The temperature dependence of soil organic matter decomposition, and the effect of global warming on soil organic C storage. Soil Biol. Biochem. 1995, 27, 753–760. [Google Scholar] [CrossRef]
  89. Martin, M.P.; Wattenbach, M.; Smith, P.; Meersmans, J.; Jolivet, C.; Boulonne, L.; Arrouays, D. Spatial distribution of soil organic carbon stocks in France. Biogeosciences 2011, 8, 1053–1065. [Google Scholar] [CrossRef]
  90. de Graaff, M.-A.; Classen, A.T.; Castro, H.F.; Schadt, C.W. Labile soil carbon inputs mediate the soil microbial community composition and plant residue decomposition rates. New Phytol. 2010, 188, 1055–1064. [Google Scholar] [CrossRef]
  91. Bradford, M.A.; Wieder, W.R.; Bonan, G.B.; Fierer, N.; Raymond, P.A.; Crowther, T.W. Managing uncertainty in soil carbon feedbacks to climate change. Nat. Clim. Chang. 2016, 6, 751–758. [Google Scholar] [CrossRef]
  92. Shi, L.; Dossa, G.G.O.; Paudel, E.; Zang, H.; Xu, J.; Harrison, R.D. Changes in fungal communities across a forest disturbance gradient. Appl. Environ. Microbiol. 2019, 85, e00080-19. [Google Scholar] [CrossRef] [PubMed]
  93. Pennanen, T.; Fitze, H.; Vanhala, P.; Kiikkilä, O.; Neuvonen, S.; Bååth, E. Structure of a microbial community in soil after prolonged addition of low levels of simulated acid rain. Appl. Environ. Microbiol. 1998, 64, 2173–2180. [Google Scholar] [CrossRef] [PubMed]
  94. Albright, M.B.; Johansen, R.; Thompson, J.; Lopez, D.; Gallegos-Graves, L.V.; Kroeger, M.E.; Runde, A.; Mueller, R.C.; Washburne, A.; Munsky, B.; et al. Soil bacterial and fungal richness forecast patterns of early pine litter decomposition. Front. Microbiol. 2020, 11, 542220. [Google Scholar] [CrossRef] [PubMed]
  95. Nilsson, L.O.; Bååth, E.; Falkengren-Grerup, U.; Wallander, H. Growth of ectomycorrhizal mycelia and composition of soil microbial communities in oak forest soils along a nitrogen deposition gradient. Oecologia 2007, 153, 375–384. [Google Scholar] [CrossRef] [PubMed]
  96. Coleman-Derr, D.; Desgarennes, D.; Fonseca-Garcia, C.; Gross, S.; Clingenpeel, S.; Woyke, T.; North, G.; Visel, A.; Partida-Martinez, L.P.; Tringe, S.G. Plant compartment and biogeography affect microbiome composition in cultivated and native Agave species. New Phytol. 2016, 209, 798–811. [Google Scholar] [CrossRef]
  97. Fonseca-García, C.; Coleman-Derr, D.; Garrido, E.; Visel, A.; Tringe, S.G.; Partida-Martínez, L.P. The Cacti Microbiome: Interplay between habitat-filtering and host-specificity. Front. Microbiol. 2016, 7, 150. [Google Scholar] [CrossRef] [PubMed]
  98. Rosling, A.; Cox, F.; Cruz-Martinez, K.; Ihrmark, K.; Grelet, G.A.; Lindahl, B.D.; Menkis, A.; James, T.Y. Archaeorhizomycetes: Unearthing an ancient class of ubiquitous soil fungi. Science 2011, 333, 876–879. [Google Scholar] [CrossRef]
  99. Qadri, M.; Rajput, R.; Abdin, M.Z.; Vishwakarma, R.A.; Riyaz-Ul-Hassan, S. Diversity, Molecular Phylogeny, and Bioactive Potential of Fungal Endophytes Associated with the Himalayan Blue Pine (Pinus wallichiana). Microb. Ecol. 2014, 67, 877–887. [Google Scholar] [CrossRef]
  100. Jeewon, R.; Yeung, Q.S.Y.; Wannasinghe, D.N.; Rampadarath, S.; Puchooa, D.; Wang, H.K.; Hyde, K.D. Hidden mycota of pine needles: Molecular signatures from PCR-DGGE and Ribosomal DNA phylogenetic characterization of novel phylotypes. Sci. Rep. 2018, 8, 18053. [Google Scholar] [CrossRef]
  101. Lynikienė, J.; Marčiulynienė, D.; Marčiulynas, A.; Gedminas, A.; Vaičiukynė, M.; Menkis, A. Managed and unmanaged Pinus sylvestris forest stands harbour similar diversity and composition of the phyllosphere and soil fungi. Microorganisms 2020, 8, 259. [Google Scholar] [CrossRef] [PubMed]
  102. Bahram, M.; Netherway, T.; Hildebrand, F.; Pritsch, K.; Drenkhan, R.; Loit, K.; Anslan, S.; Bork, P.; Tedersoo, L. Plant nutrient—Acquisition strategies drive topsoil microbiome structure and function. New Phytol. 2020, 227, 1189–1199. [Google Scholar] [CrossRef]
  103. Hou, J.; Feng, H.; Wu, M. Incorporating Effect Factors into the Relationship between Biodiversity and Ecosystem Functioning (BEF). Diversity 2022, 14, 274. [Google Scholar] [CrossRef]
  104. de la Fuente Cantó, C.; Simonin, M.; King, E.; Moulin, L.; Bennett, M.J.; Castrillo, G.; Laplaze, L. An extended root phenotype: The rhizosphere, its formation and impacts on plant fitness. Plant J. 2020, 103, 951–964. [Google Scholar] [CrossRef] [PubMed]
  105. Bulgarelli, D.; Schlaeppi, K.; Spaepen, S.; Van Themaat, E.V.L.; Schulze-Lefert, P. Structure and functions of the bacterial microbiota of plants. Annu. Rev. Plant Biol. 2013, 64, 807–838. [Google Scholar] [CrossRef] [PubMed]
  106. Hanso, M.; Drenkhan, R. Lophodermium needle cast, insect defoliation and growth responses of young Scots pines in Estonia. For. Pathol. 2012, 42, 124–135. [Google Scholar] [CrossRef]
  107. Millberg, H. Foliar Fungi of Scots Pine (Pinus sylvestris). Ph.D. Thesis, Department of Forest Mycology and Plant Pathology, Swedish University of Agricultural Sciences, Uppsala, Zweden, 2015; pp. 1652–6880. [Google Scholar]
  108. Haas, J.C.; Street, N.R.; Sjödin, A.; Lee, N.M.; Högberg, M.N.; Näsholm, T.; Hurry, V. Microbial community response to growing season and plant nutrient optimisation in a boreal Norway spruce forest. Soil Biol. Biochem. 2018, 125, 197–209. [Google Scholar] [CrossRef]
  109. Policelli, N.; Bruns, T.D.; Vilgalys, R.; Nuñez, M.A. Suilloid fungi as global drivers of pine invasions. New Phytol. 2019, 222, 714–725. [Google Scholar] [CrossRef]
  110. Peay, K.G.; Baraloto, C.; Fine, P.V.A. Strong coupling of plant and fungal community structure across western Amazonian rainforests. ISME J. 2013, 7, 1852–1861. [Google Scholar] [CrossRef]
  111. van der Linde, S.; Suz, L.M.; Orme, C.D.L.; Cox, F.; Andreae, H.; Asi, E.; Atkinson, B.; Benham, S.; Carroll, C.; Cools, N.; et al. Author Correction: Environment and host as large-scale controls of ectomycorrhizal fungi. Nature 2018, 561, E42. [Google Scholar] [CrossRef]
  112. Baar, J.; Horton, T.R.; Kretzer, A.; Bruns, T.D. From Resistant Propagules After a Stand- Replacing Wildfire. New Phytol. 1999, 143, 409–418. [Google Scholar] [CrossRef]
  113. Qu, Z.L.; Santalahti, M.; Köster, K.; Berninger, F.; Pumpanen, J.; Heinonsalo, J.; Sun, H. Soil Fungal Community Structure in Boreal Pine Forests: From Southern to Subarctic Areas of Finland. Front. Microbiol. 2021, 12, 653896. [Google Scholar] [CrossRef] [PubMed]
  114. Hoff, J.A.; Klopfenstein, N.B.; Mcdonald, G.I.; Tonn, J.R.; Kim, M.S.; Zambino, P.J.; Hessburg, P.F.; Rogers, J.D.; Peever, T.L.; Carris, L.M. Fungal endophytes in woody roots of Douglas-fir (Pseudotsuga menziesii) and ponderosa pine (Pinus ponderosa). For. Pathol. 2004, 34, 255–271. [Google Scholar] [CrossRef]
  115. Summerbell, R.C. Root endophyte and mycorrhizosphere fungi of black spruce, Picea mariana, in a boreal forest habitat: Influence of site factors on fungal distributions. Stud. Mycol. 2005, 53, 121–145. [Google Scholar] [CrossRef]
  116. Bååth, E. Microfungi in a clear-cut pine forest soil in central Sweden. Can. J. Bot. 1981, 59, 1331–1337. [Google Scholar] [CrossRef]
  117. 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]
  118. Halmschlager, E.; Kowalski, T. The mycobiota in nonmycorrhizal roots of healthy and declining oaks. Can. J. Bot. 2004, 82, 1446–1458. [Google Scholar] [CrossRef]
  119. Sheng, Y.; Cong, J.; Lu, H.; Yang, L.; Liu, Q.; Li, D.; Zhang, Y. Broad-leaved forest types affect soil fungal community structure and soil organic carbon contents. Microbiologyopen 2019, 8, e874. [Google Scholar] [CrossRef]
  120. Crous, P.W.; Quaedvlieg, W.; Hansen, K.; Hawksworth, D.L.; Groenewald, J.Z. Phacidium and Ceuthospora (Phacidiaceae) are congeneric: Taxonomic and nomenclatural implications. IMA Fungus 2014, 5, 173–193. [Google Scholar] [CrossRef]
  121. Nawrot-Chorabik, K.; Grad, B.; Kowalski, T. Interactions between callus cultures of Pinus sylvestris and pine fungi with different trophic properties. For. Pathol. 2016, 46, 179–186. [Google Scholar] [CrossRef]
  122. Ozturk, I.K.; Amiri, A. Pathogenicity and Control of Phacidium lacerum, an Emerging Pome Fruit Pathogen in Washington State. Plant Dis. 2020, 104, 3124–3130. [Google Scholar] [CrossRef] [PubMed]
  123. Silva, A.C.; Henriques, J.; Diogo, E.; Ramos, A.P.; Bragança, H. First report of Sydowia polyspora causing disease on Pinus pinea shoots. For. Pathol. 2020, 50, 27–30. [Google Scholar] [CrossRef]
  124. Cleary, M.; Oskay, F.; Doğmuş, H.T.; Lehtijärvi, A.; Woodward, S.; Vettraino, A.M. Cryptic risks to forest biosecurity associated with the global movement of commercial seed. Forests 2019, 10, 459. [Google Scholar] [CrossRef]
  125. Davydenko, K.; Vasaitis, R.; Meshkova, V.; Menkis, A. Fungi associated with the red-haired bark beetle, Hylurgus ligniperda (Coleoptera: Curculionidae) in the forest-steppe zone in eastern Ukraine. Eur. J. Entomol. 2014, 111, 561–565. [Google Scholar] [CrossRef]
  126. Lygis, V.; Vasiliauskaite, I.; Matelis, A.; Pliūra, A.; Vasaitis, R. Fungi in living and dead stems and stumps of pinus mugo on coastal dunes of the baltic sea. Plant Prot. Sci. 2014, 50, 221–226. [Google Scholar] [CrossRef]
  127. Behnke, J.; Hanna, B. Fungi associated with Cyclaneusma needle cast in Scots pine in the west of Poland. For. Pathol. 2019, 49, e12487. [Google Scholar] [CrossRef]
Figure 1. Map of Lithuania showing the sampling sites in Pinus sylvestris forest stands. The green color on the map shows the distribution of P. sylvestris stands and its abundance (%) in different areas.
Figure 1. Map of Lithuania showing the sampling sites in Pinus sylvestris forest stands. The green color on the map shows the distribution of P. sylvestris stands and its abundance (%) in different areas.
Jof 08 01112 g001
Figure 2. Differences in the Shannon diversity index of fungal OTUs among different substrates. In the Mann–Whitney test, p-values were: roots vs. shoots p < 0.05, roots vs. needles p < 0.05, roots vs. organic soil p < 0.05, root vs. mineral soil p < 0.05, shoot vs. needles p = 0.97, shoots vs. organic soil p = 0.06, shoots vs. mineral soil p = 0.23, needles vs. organic soil p = 0.06, needles vs. mineral soil p = 0.06, and organic soil vs. mineral soil p = 0.99.
Figure 2. Differences in the Shannon diversity index of fungal OTUs among different substrates. In the Mann–Whitney test, p-values were: roots vs. shoots p < 0.05, roots vs. needles p < 0.05, roots vs. organic soil p < 0.05, root vs. mineral soil p < 0.05, shoot vs. needles p = 0.97, shoots vs. organic soil p = 0.06, shoots vs. mineral soil p = 0.23, needles vs. organic soil p = 0.06, needles vs. mineral soil p = 0.06, and organic soil vs. mineral soil p = 0.99.
Jof 08 01112 g002
Figure 3. Venn diagram showing the diversity and overlap of fungal taxa in different sample types from P. sylvestris stands. Different colors represent different substrates: Yellow—Soil (organic and mineral combined), Orange—Needles, Blue—Shoots, and Gray—Roots.
Figure 3. Venn diagram showing the diversity and overlap of fungal taxa in different sample types from P. sylvestris stands. Different colors represent different substrates: Yellow—Soil (organic and mineral combined), Orange—Needles, Blue—Shoots, and Gray—Roots.
Jof 08 01112 g003
Figure 4. Species accumulation curves showing the relationship between the cumulative number of fungal OTUs and the number of ITS rRNA sequences shown as Yellow—Soil, Orange—Needles, Blue—Shoots, and Gray—Roots.
Figure 4. Species accumulation curves showing the relationship between the cumulative number of fungal OTUs and the number of ITS rRNA sequences shown as Yellow—Soil, Orange—Needles, Blue—Shoots, and Gray—Roots.
Jof 08 01112 g004
Figure 5. Relative abundance of fungal classes in needles, shoots, roots, and the soil of Pinus sylvestris. Other represented fungal classes with a relative abundance of <1%. The data from the 30 different sample sites are combined.
Figure 5. Relative abundance of fungal classes in needles, shoots, roots, and the soil of Pinus sylvestris. Other represented fungal classes with a relative abundance of <1%. The data from the 30 different sample sites are combined.
Jof 08 01112 g005
Figure 6. Fungal functional groups in the needle, shoot, root, and soil samples from the 30 sampling sites of Pinus sylvestris are shown as a relative abundance (%) of the total number of fungal taxa.
Figure 6. Fungal functional groups in the needle, shoot, root, and soil samples from the 30 sampling sites of Pinus sylvestris are shown as a relative abundance (%) of the total number of fungal taxa.
Jof 08 01112 g006
Figure 7. The relationship between the species richness of fungal OTUs in different substrates (needles, shoots, roots, mineral, and organic soil) of Pinus sylvestris and different stand (A,B) and climatic (C,D) parameters. The semitransparent field around each curve denotes the size of the deviation from the mean value.
Figure 7. The relationship between the species richness of fungal OTUs in different substrates (needles, shoots, roots, mineral, and organic soil) of Pinus sylvestris and different stand (A,B) and climatic (C,D) parameters. The semitransparent field around each curve denotes the size of the deviation from the mean value.
Jof 08 01112 g007
Figure 8. The relationship between the species richness of fungal OTUs in the different substrates (needles, shoots, roots, mineral, and organic soil) of Pinus sylvestris and different soil parameters (AD). The semitransparent field around each curve denotes the size of the deviation from the mean value.
Figure 8. The relationship between the species richness of fungal OTUs in the different substrates (needles, shoots, roots, mineral, and organic soil) of Pinus sylvestris and different soil parameters (AD). The semitransparent field around each curve denotes the size of the deviation from the mean value.
Jof 08 01112 g008
Figure 9. Non-metric multidimensional scaling (NMDS) of fungal communities associated with different substrates (N—needles, R—roots, Sh—shoots, Sm—mineral soil, and So—organic soil) of Pinus sylvestris.
Figure 9. Non-metric multidimensional scaling (NMDS) of fungal communities associated with different substrates (N—needles, R—roots, Sh—shoots, Sm—mineral soil, and So—organic soil) of Pinus sylvestris.
Jof 08 01112 g009
Table 1. Principal information on the 30 sampling sites in Pinus sylvestris forest stands.
Table 1. Principal information on the 30 sampling sites in Pinus sylvestris forest stands.
SitePositionTree SpeciesAgeSoil Chemical ParametersClimate Data
No.Longitude (N)Latitude (E)Composition, % *(y)pH, mol/L (KCl)P2O5, mg/kgK2O, mg/kgCa, mg/kgMg, mg/kgCl, mg/kgSalts, ms/cmAverage Annual Temp. °CPrecipitation, mm/year
155°37′14″25°11′25″100P604.034525155483.62.558.0508.8
253°58′52″23°56′56″100P984.29331212363.61.618.0484.6
353°58′27″24°25′48″100P534.212325256743.61.898.0484.6
454°8′20″24°11′51″100P534.46423260603.61.448.0484.6
554°52′40″23°42′17″80P,20B674.16019205503.62.188.3614.9
655°10′27″24°26′43″100P574.013919234555.32.418.3614.9
755°1′10″24°12′56″100P713.51957181455.32.718.6519.8
854°52′15″23°26′41″100P493.91861298623.63.458.9491.7
956°11′40″21°28′13″100P593.834656431155.34.458.9491.7
1056°1′1″21°6′45″90P,10S474.09231246463.61.748.3643.6
1155°45′36″24°41′38″70P,20S,10B723.41336191503.63.078.1479.8
1255°45′36″24°41′38″80P,20B685.9716041006706.920.28.1505.5
1356°11′12″22°23′4″100P584.2258620812725.35.188.0480.4
1454°52′16″25°41′56″100P654.14157339683.62.818.0755.4
1554°34′28″23°57′2″100P474.011541243573.62.218.6595.2
1655°28′16″23°26′39″100P554.27039373665.32.087.8458.0
1755°29′54″21°57′21″70P,20S,10B652.813838111712717.115.97.7514.6
1856°3′26″25°42′28″90P,10S593.91051208413.63.028.4617.0
1955°1′46″22°42′14″70P,20J,10B363.97830178533.61.788.6519.8
2054°17′17″25°39′41″90P,10S503.94444191545.33.298.0530.6
2155°10′12″25°41′56″100P654.4571397283.62.627.4669.2
2255°19′38″22°27′14″70P,20B,10S603.5161748342813.658.6591.3
2355°45′34″21°43′12″70P,20S,10B653.43083383815.34.298.0598.7
2454°44′58″24°41′59″100P585.92716716242153.68.497.3614.9
2554°56′54″24°41′46″100P554.910429399685.33.198.0535.8
2655°10′47″24°42′10″100P503.913036183523.62.728.0506.9
2755°19′1″25°42′54″100P654.53748675975.34.867.7505.4
2854°25′23″24°57′33″90P,10B484.11433180445.32.58.0593.2
2954°25′39″24°27′15″90P,10B704.33848323663.63.888.0540.9
3054°7′57″23°41′50″100P695.3624010531153.64.038.4565.2
* P—Pinus sylvestris, S—Picea abies, B—Betula pendula, J—Alnus glutinosa. Tree species composition is based on volume.
Table 2. The number of high-quality sequences and fungal OTUs from each study site.
Table 2. The number of high-quality sequences and fungal OTUs from each study site.
SiteNo. of High-Quality Sequences/OTUsShannon Diversity Index (H)
No.RootsShootsNeedlesSoil O *Soil M *Soil O + MRootsShootsNeedlesSoil O *Soil M *Soil O + M
16569/74271/431117/1363165/283209/773374/3112.282.783.784.203.844.29
22558/117224/591529/134414/85639/1481053/1932.833.343.263.354.064.19
32189/91597/612518/1251092/145863/1931955/2793.212.612.993.584.574.34
41107/503779/134870/613205/2371962/1315167/3122.602.912.823.643.533.91
5907/67458/93770/853180/273742/873922/3053.043.883.063.203.053.40
6139/3349/301846/182236/101390/136626/1913.063.254.034.054.004.27
73026/108517/1041974/1741/121/822/92.933.813.54-1.591.70
82390/971006/125479/76708/1421609/1572317/2463.303.533.354.343.704.31
9969/88708/1192286/1801427/2341888/2083315/3223.183.843.754.474.024.42
101309/1261715/2224841/32539/292407/1842446/1983.044.194.243.243.803.85
11411/1043032/2481058/128382/97191/57573/1253.764.063.783.823.504.05
12885/893134/237689/971074/223731/801805/2533.184.032.584.662.934.37
1375/362018/1693590/1886088/2423955/26210043/3713.303.803.013.814.084.09
145076/1031726/156677/1021275/2671018/2022293/3672.473.683.234.864.655.11
15223/561252/1154097/2191508/1521339/1372847/2453.283.563.833.713.204.02
1679/144344/2213936/232168/7587/46255/1071.973.673.773.713.484.02
1775/303770/1943037/178883/123640/1381523/2313.043.853.293.734.104.42
18841/515417/2023649/209113/48285/78398/1061.973.683.873.153.543.75
19174/505273/1865146/20427/21152/79179/933.213.363.372.973.974.12
20139/382440/1503122/1401355/242574/1061929/2893.003.433.284.813.864.84
21121/42913/862819/189701/142234/92935/2103.183.183.194.184.174.61
22-/-2368/1534445/21687/5461/31148/83-3.373.713.703.354.04
23169/334620/2331885/1371047/1281752/2002799/2862.563.672.863.474.404.54
24219/412890/1582572/14868/16136/35204/452.653.893.391.841.651.88
252193/611490/155682/894009/3451721/1595730/3961.623.933.264.394.064.57
26149/354207/239397/891575/183553/1122128/2452.523.953.663.573.373.74
271395/884081/159417/713611/3212120/1485731/3922.713.363.084.403.484.44
281868/109652/851102/1191257/1661188/1952445/2983.123.323.183.624.344.39
29699/585143/217710/911960/2422598/2184558/3642.583.473.484.564.214.67
30-/-10667/233776/1001502/1661234/1832736/276-2.783.423.834.134.27
Total35954/73478760/90763031/95342117/144031254/123373371/1854
* O—organic layer, M—mineral layer.
Table 3. Relative abundance (%) of fungal OTUs associated with needles, shoots, roots, and organic and mineral soil of Pinus sylvestris. All study sites are combined.
Table 3. Relative abundance (%) of fungal OTUs associated with needles, shoots, roots, and organic and mineral soil of Pinus sylvestris. All study sites are combined.
Fungal OTUPhylum *Genbank/UNITE ReferenceSequence Similarity, %Needles, %Shoots, %Roots, %Soil O **, %Soil M **, %Soil all, %All, %
Dothideomycetes sp. 5208_5AKX9084729920.4453.3890.0393.2171.7092.5756.953
Unidentified sp. 5208_1AKP8913981004.31412.7860.0171.3631.6541.4875.530
Archaeorhizomyces sp. 5208_0AMH2480431001.2310.3828.3528.3939.8399.0094.257
Unidentified sp. 5208_2AMN9023671004.4067.0960.0112.2751.8142.0783.940
Helotiales sp. 5208_17AKY7425931001.5776.1360.0110.6170.5730.5982.497
Rhinocladiella sp. 5208_3AKM056296980.3986.4350.0080.7690.5180.6622.313
Phacidium lacerumAMN5881631004.5150.1130.0363.0391.2482.2761.839
Phialocephala fortiniiAMN9473951000.1600.04111.0560.4680.4580.4631.771
Scoliciosporum umbrinumAKX1330081000.8874.2470.0030.4770.6940.5701.722
Sydowia polysporaAMN9006301003.4281.8940.0310.5670.4670.5251.612
Mycena cinerellaBKT9001461000.0240.0098.8471.1350.2050.7391.491
Malassezia restrictaBLT8546971001.2280.6700.4561.8543.1802.4191.291
Microsphaeropsis olivaceaAMT5613961002.5880.8190.1110.8031.6411.1601.261
Penicillium camembertiAMT3555661000.1590.0245.3261.6481.0371.3871.215
Phaeomoniella pinifoliorumAMK7625951000.9821.7690.0060.5580.4830.5260.956
Unidentified sp. 5208_36AMG8283111000.6952.0510.0080.4800.3680.4320.945
Cladosporium herbarumAMT6352881001.8260.4320.6310.7840.8480.8110.921
Chaetothyriomycetidae sp. 5208_56AKX589170980.2272.400-0.0780.4830.2510.883
Leptosphaeria sp. 5208_64AJQ044439970.2522.417-0.0240.0220.0230.828
Unidentified sp. 5208_39AMT2420101001.5991.0230.0060.3800.2620.3300.820
Hyaloscypha variabilisAMT4699251000.057-3.2211.0991.1421.1180.802
Umbelopsis dimorphaZMT1386161000.1140.0040.0281.8853.2092.4490.749
Chaetothyriales sp. 5208_15AKP4005721000.8361.3460.0190.2970.3810.3330.732
Trechispora sp. 5208_19BJX392812990.002-3.9190.5820.4770.5370.718
Unidentified sp. 5208_12AFJ5535821000.0620.0040.0062.1232.7262.3800.713
Unidentified sp. 5208_32BMN9023631002.1890.325-0.2520.0700.1740.702
Cenococcum geophilumAHM1897241000.010-2.5170.5601.8171.0960.683
Unidentified sp. 5208_72AMN9023961001.6370.4580.0140.5250.2240.3970.673
Unidentified sp. 5208_23AMT2370781001.1090.6840.0060.5650.4610.5210.646
Unidentified sp. 5208_69AMN9023871001.4740.543-0.2160.0830.1590.587
* A—Ascomycota, B—Basidiomycota, Z—Zygomycota; ** O—organic layer, M—mineral layer.
Table 4. Pairwise comparison of the fungal communities among different substrates (needles, roots, shoots, mineral soil, and organic soil) of Pinus sylvestris.
Table 4. Pairwise comparison of the fungal communities among different substrates (needles, roots, shoots, mineral soil, and organic soil) of Pinus sylvestris.
SubstrateNeedlesRootsShootsSoilM
R2FpR2FpR2FpR2Fp
Roots0.3724.190.0011---------
Shoots0.2314.860.00110.3622.060.0011------
SoilM0.2615.640.00110.208.250.00110.2514.500.0011---
SoilO0.2111.880.00110.187.680.00110.2011.290.00110.020.900.591
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Marčiulynienė, D.; Marčiulynas, A.; Mishcherikova, V.; Lynikienė, J.; Gedminas, A.; Franic, I.; Menkis, A. Principal Drivers of Fungal Communities Associated with Needles, Shoots, Roots and Adjacent Soil of Pinus sylvestris. J. Fungi 2022, 8, 1112. https://doi.org/10.3390/jof8101112

AMA Style

Marčiulynienė D, Marčiulynas A, Mishcherikova V, Lynikienė J, Gedminas A, Franic I, Menkis A. Principal Drivers of Fungal Communities Associated with Needles, Shoots, Roots and Adjacent Soil of Pinus sylvestris. Journal of Fungi. 2022; 8(10):1112. https://doi.org/10.3390/jof8101112

Chicago/Turabian Style

Marčiulynienė, Diana, Adas Marčiulynas, Valeriia Mishcherikova, Jūratė Lynikienė, Artūras Gedminas, Iva Franic, and Audrius Menkis. 2022. "Principal Drivers of Fungal Communities Associated with Needles, Shoots, Roots and Adjacent Soil of Pinus sylvestris" Journal of Fungi 8, no. 10: 1112. https://doi.org/10.3390/jof8101112

APA Style

Marčiulynienė, D., Marčiulynas, A., Mishcherikova, V., Lynikienė, J., Gedminas, A., Franic, I., & Menkis, A. (2022). Principal Drivers of Fungal Communities Associated with Needles, Shoots, Roots and Adjacent Soil of Pinus sylvestris. Journal of Fungi, 8(10), 1112. https://doi.org/10.3390/jof8101112

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