Petra Veselá
Martina Vašutová
Magda Edwards-Jonášová3
Pavel Cudlín- 1Department of Carbon Storage in the Landscape, Global Change Research Institute of the Czech Academy of Sciences, Brno, Czechia
- 2Department of Forest Protection and Wildlife Management, Faculty of Forestry and Wood Technology, Mendel University in Brno, Brno, Czechia
- 3Department of Carbon Storage in the Landscape, Global Change Research Institute of the Czech Academy of Sciences, České Budějovice, Czechia
- 4Department of Botany, Faculty of Science, University of South Bohemia, České Budějovice, Czechia
- 5Department of Integrated Forest and Landscape Protection, Faculty of Forestry, Technical University in Zvolen, Zvolen, Slovakia
Due to ongoing climate change, forests are expected to face significant disturbances more frequently than in the past. Appropriate management is intended to facilitate forest regeneration. Because European temperate forests mostly consist of trees associated with ectomycorrhizal (ECM) fungi, understanding their role in these disturbances is important to develop strategies to minimize their consequences and effectively restore forests. Our aim was to determine how traditional (EXT) and nonintervention (NEX) management in originally Norway spruce (Picea abies) forests with an admixture of European larch (Larix decidua) affect ECM fungal communities and the potential to interconnect different tree species via ECM networks 15 years after a windstorm. Ten plots in NEX and 10 plots in EXT with the co-occurrences of Norway spruce, European larch, and silver birch (Betula pendula) were selected, and a total of 57 ECM taxa were identified using ITS sequencing from ECM root tips. In both treatments, five ECM species associated with all the studied tree species dominated, with a total abundance of approximately 50% in the examined root samples. Because there were no significant differences between treatments in the number of ECM species associated with different tree species combinations in individual plots, we concluded that the management type did not have a significant effect on networking. However, management significantly affected the compositions of ECM symbionts of Norway spruce and European larch but not those of silver birch. Although this result is explained by the occurrence of seedlings and ECM propagules that were present in the original forest, the consequences are difficult to assess without knowledge of the ecology of different ECM symbionts.
Introduction
Windstorms represent one of the major factors influencing European temperate forests (Fischer et al., 2013), and their intensity has been increasing during the last few decades (Gregow et al., 2017). Depending on wind severity, the influence of windstorms may extend from single trees to large areas of forests. Large windstorms, causing forest breakdown, affect the forest ecosystem in a complex way. Changes in forest structure lead to overall changes in microhabitat conditions, such as increased solar radiation connected to increased temperatures and changes in water availability (Gömöryová et al., 2008), nutrient cycling (Don et al., 2012), and vegetation cover (Dietz et al., 2020). Because European temperate forests mainly consist of trees living in symbiosis with ectomycorrhizal (ECM) fungi (Martin et al., 2016), these fungal symbionts also face these changes.
Severe windstorms generally lead to decrease in ECM fungal diversity and changes in their community composition (Egli et al., 2002; Ford et al., 2018; Vašutová et al., 2018). Subsequent management practices, such as clear-cut logging, usually enhance these changes as a possible consequence of soil mechanical damage (Ford et al., 2018; Vašutová et al., 2018). The main driver of fungal community shifts is the death of mature trees (Karst et al., 2014; Treu et al., 2014), because ECM fungi are obligate symbionts unable to live without their hosts (Baldrian, 2009; Lofgren et al., 2018). This results in an extinction of ECM fungal species associated with mature trees (= late-successional ECM fungi), while early-successional ECM fungi, with a low demand on photosynthates, survive on seedlings (Egli et al., 2002; Vašutová et al., 2018; Veselá et al., 2019b). Fungal inoculum may further persist in soil for some time as spores or hyphae emanating from dying or recently dead root tips (Simard, 2009; Guignabert et al., 2018). ECM species known to be able to survive or rapidly re-establish after disturbances are, e.g., Cenococcum sp., Laccaria sp., Piloderma sp., Rhizopogon sp., Thelephora sp., or Wilcoxina sp. (Visser, 1995; Jones et al., 2003; Kranabetter et al., 2004; Twieg et al., 2007).
The ability to succeed in young forest stages depends on various functional traits of the ECM fungal species, i.e., on their enzymatic activity, which enables them to decompose various substrates and mobilize nutrients (Velmala et al., 2014; Haas et al., 2018); spore traits which improve their resistance in soil (Jones et al., 2003); and morphology of mycorrhizas represented by exploration type which influences nutrient exploitation from soil and interconnection of trees through the mycorrhizal network (Agerer, 2001; Rosinger et al., 2018). For example, some Piloderma species secrete hemicellulases and hydrolytic enzymes involved in litter degradation, while Wilcoxina sp. is effective in degrading chitin (Velmala et al., 2014). Higher nitrate uptake and retranslocation to the plant host was also observed in Laccaria laccata compared with Suillus bovinus (Hobbie et al., 2008). The ascomycete species Cenococcum geophilum is supposed to be low carbon demanding due to having a thin mantel and short-distance exploration type (Guignabert et al., 2018; Defrenne et al., 2019), while long-distance explorers belonging to the genus Rhizopogon were found to dominate in dry, nutrient-poor areas, since they are able to transport resources from longer distances (Bakker et al., 2006).
Each individual tree can be associated with many ECM fungal species (Bahram et al., 2011), and each ECM fungal species differs in its specificity (Molina and Horton, 2015). ECM fungi with narrow host range (e.g., some Leccinum, Rhizopogon, and Suillus species) are usually restricted to a specific host genus, while ECM fungi with a broad host range (e.g., Boletus, Paxillus, and Laccaria species) are known to be associated with nearly every tree species that form ectomycorrhizae (Molina et al., 1992). ECM fungi with broad host ranges dominate in ECM fungal communities (Molina et al., 1992; Horton and Bruns, 1998; Taudiere et al., 2015; van der Linde et al., 2018; Defrenne et al., 2019; Rog et al., 2020). These fungi can interconnect different tree species and influence nutrient exchange among them (Simard et al., 1997; van der Heiden and Horton, 2009). On the other hand, almost all host tree species harbor in their symbiotic fungal assembly genus-restricted ECM fungi, which provide them exclusive access to resources (Molina et al., 1992). Nevertheless, ECM fungi with a broad host range may also show host preference depending on the ecological context (Lang et al., 2011), while genus-restricted ECM fungi may expand their host range in the absence of their primary host (Lofgren et al., 2018).
The ECM fungal community is a dynamic system interconnecting trees of the same species, but also different tree species through the mycorrhizal network (Simard et al., 1997; Klein et al., 2016; Rog et al., 2020). This phenomenon is not negligible when considering that up to approximately 40% of fine-root carbon can originate from mycorrhizal network exchange (Klein et al., 2016). Not only carbon can be transferred through the mycorrhizal network, but also water and other nutrients, or signaling molecules (Gorzelak et al., 2015). The mycorrhizal network is influenced by tree physiology, e.g., photosynthetic and growth rates, nutrient content, forest age, or infestation (Simard et al., 1997, 2012; van der Heiden and Horton, 2009; Song et al., 2015; Simard, 2018); the phylogenetic relationships among trees involved in the network (Rog et al., 2020); environmental conditions (Rosinger et al., 2018; Defrenne et al., 2019); and also by the ECM fungi themselves which protect their carbon donors (Song et al., 2015). This system is expected to be very dynamic and vulnerable to large disturbances, as demonstrated by Pec et al. (2020).
Severe windstorm damaged approximately 12,000 ha of mountain spruce forest in the Tatra Mts. (Slovakia) in 2004 (Fleischer, 2008). In most of the affected area, the fallen trees were extracted (EXT), while a small part was left for natural succession (NEX). To bring new perspectives to the discussion on the suitability of subsequent management, we were interested in how the windstorm and following timber harvest affected the ECM fungal community and its potential to form mycorrhizal networks among regenerating tree species, Norway spruce, European larch, and silver birch. We hypothesized that ECM fungi with a broad host range will prevail in plots with traditional management (e.g., fallen tree extraction) and that different tree species will therefore share a high percentage of ECM fungal species. On the other hand, plots left for natural succession will harbor more ECM fungi with narrow host range associated with target tree species, with fewer ECM fungal species shared among them. We also hypothesized that the number of ECM fungal species would decrease in managed plots due to soil compaction and mechanical damage to seedlings that survive and the ECM fungal extramatrical mycelium.
Materials and Methods
Study Site
The research plots are situated in the Tatra Mts. (Slovakia) in the Tatra National Park (Supplementary Figure 1). This area was affected by severe windstorm in 2004 during which about 12,000 ha of forest was destroyed. The area was originally covered by a Norway spruce (Picea abies) forest with an admixture of European larch (Larix decidua), and Scotch pine (Pinus sylvestris). Currently, this area is dominated by young trees of Norway spruce, birches (Betula sp.), European larch, and European rowan (Sorbus aucuparia) in variable proportions depending on the microhabitat conditions. Locally, Scots pine, goat willow (Salix caprea), and silver fir (Abies alba) are present. The area is characterized by a cold climate and has a mean annual temperature of 5.3°C and a mean annual precipitation of 833 mm (Fleischer, 2008). The soils in the plots are dystric cambisols on glacial moraine deposits with a loamy sand texture, an acidic pH of approximately 3.0 (Don et al., 2012), and carbon and nitrogen percentages of approximately 11.9 and 0.6%, respectively (Vašutová et al., 2018).
Following this disturbance event, permanent research sites with different management regimes were established in 2005, and each research site is about 100 ha wide (Fleischer, 2008). Within the current research, we studied two treatments: (i) site with traditional management (EXT), where fallen trees were extracted and the locality was consequently reforested, and (ii) a site without management (NEX) left for natural succession, in which fallen trees were left. Ten circular plots (d = 12.6 m and P = 125 m2) per treatment were selected so that there were at least two individuals of each tree species of interest, Norway spruce, European larch, and silver birch (Betula pendula), and the presence of other ECM tree species was avoided to minimize influencing the ECM fungal community.
Ectomycorrhizae on Tree Roots
ECM root samples (each consists of about 100 root tips) from 20 plots were collected during 2 years: three plots per treatment in September 2018 and seven plots per treatment in September 2019 (Supplementary Table 1). In each plot, two individuals for each tree species were selected, and three root samples from each individual tree were collected. The cover of trees, shrubs, herbs, mosses, and dead wood was visually estimated (Supplementary Table 1). Dead wood was classified into five categories (lying logs, piles of branches, root collars, stubs, and stumps).
The height and perimeter of each selected tree were measured (Supplementary Table 1). For each root sample, its depth, soil horizon, distance from the tree stem, and direction to other tree species were recorded (Supplementary Table 1). In total, 60 root samples from each tree species per treatment (2 trees × 3 root samples per plot × 10 plots) were sampled, with the exception of Norway spruce (due to the misidentification of one of the root samples, only 59 samples were analyzed).
The identity of the roots was guaranteed by direct tracing from a tree stem; in disputable cases, the identity was confirmed using chloroplast DNA sequence comparison against the GenBank database (Benson et al., 2007) using MegaBlast algorithm, with 97% similarity threshold as a molecular species criterion.
Processing of Root Samples
Each root sample was carefully washed under tap water, and all attached soil particles were removed using tweezers in a Petri dish with water. The root samples were examined under a stereomicroscope and documented (Canon 1,000D, Tokyo, Japan). The ECM root tips of each root sample were sorted into morphotypes according to their morphological features: color, presence of emanating hyphae and rhizomorphs, cystidia, and ramification (Agerer, 1991; Agerer and Rambold, 2004–2015). Each morphotype was excised and stored in a tube with silica gel. At least one morphotype per tree species and treatment was then used for molecular identification.
DNA Extraction and Molecular Identification
Genomic DNA was extracted from the selected ECM root tips using the DNeasy Plant Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. The extracted DNA was then used as a template for PCR amplification of the rDNA ITS region with the primer pairs ITS1F/ITS4 (White et al., 1990; Gardes and Bruns, 1993) for fungal species identification and of the noncoding region of chloroplast DNA with the primers trnL/trnF (Taberlet et al., 1991) for plant species identification. In case that we obtained multiple fungal amplicons from a root tip, we used primer pairs specific for basidiomycetous (ITS1F/ITS4B; Gardes and Bruns, 1993; Nikolcheva and Bärlocher, 2004) or ascomycetous fungi (ITS1/ITS4A; White et al., 1990; Nikolcheva and Bärlocher, 2004). When the amplicon separation was unsuccessful, we cut the bands directly from gel using QIAquick Gel Extraction Kit (Qiagen). PCR was performed in a 25-μl reaction mixture containing 1× reaction buffer (Bioline, London, United Kingdom), 500 nM of each primer, 0.75 U of MyTaq DNA polymerase (Bioline), and ∼20 ng of template DNA filled with dH2O up to 25 μl. The amplification conditions were 94°C for 2.5 min, followed by 37 cycles at 94°C for 30 s, 55°C for 40 s, 72°C for 30 s, and a final extension at 72°C for 4.5 min for the fungal amplicons; 94°C for 2 min, followed by 35 cycles at 94°C for 30 s, 56°C for 30 s, 72°C for 1.2 min, and a final extension at 72°C for 10 min for the plant amplicons. Sequencing was performed by Macrogen Inc., Amsterdam, Netherlands, on an ABI 3730 XL automated sequencer (Applied Biosystems, Foster City, CA, United States).
Molecular Data Analysis
Sequence chromatograms were analyzed using Chromas Lite v.2.6.5 (Technelysium Pty Ltd, Brisbane, Australia), and the sequences were blasted against the UNITE (Abarenkov et al., 2010) and GenBank (Benson et al., 2007) public databases. A 97% similarity threshold was used as a molecular species criterion and 90% similarity threshold as a molecular genus criterion (Nilsson et al., 2019; Vu et al., 2019). The sequences were deposited in GenBank (accession numbers: MT908274–MT908330). Based on the molecularly identified morphotypes, species names were given to all the documented morphotypes.
The assignment of ECM taxa into exploration types was made according to Agerer (2001), Agerer and Rambold (2004–2015), and our observations (Supplementary Table 2).
Statistical Analysis
The numbers of ECM fungal species shared by all tree species, a combination of two tree species, or unique for a single tree species were calculated for the entire dataset, EXT and NEX treatments and separately for each plot. The same approach was used in assessing the number of species with the exploration types with rhizomorphs (long-distance exploration type, medium-distance fringe exploration type, and medium-distance smooth exploration type), which could in particular provide networking.
The tree heights and perimeters per tree/tree species for the NEX and EXT treatments and the numbers of ECM fungal species per tree/tree species for the NEX and EXT treatments were tested by analysis of variance with Tukey’s post hoc contrasts in R 3.6.2 (R Development Core Team, 2019).
Rarefaction curves, using R package mobr (McGlinn et al., 2018), were prepared for all tree species within the NEX and EXT treatments (Supplementary Figure 2).
Three datasets were prepared for relating the plot variables (n = 20), trees (n = 120), and root samples (n = 359). Presence/absence data (roots) and number of root samples with each morphotype (trees, plots) were used to avoid difficulties in the proper identification of not well-developed morphotypes. The datasets were separately analyzed by multivariate methods in Canoco 5 program (ter Braak and Šmilauer, 2012), using detrended correspondence analysis (DCA) and canonical correspondence analysis (CCA). DCA was used to compare the ECM fungal species compositions in plots, trees, and roots. The year of sampling was used as a covariate. CCA was used to identify (i) the effects of management, vegetation cover, and geographical position on ECM fungal species composition in plots; (ii) the effects of management, tree species, height, and perimeter on ECM fungal species composition in trees; and (iii) the effects of management, tree species, distance from tree stem, and soil horizon on the ECM fungal species composition in roots. The dataset trees was divided based on tree species identity, and the resulting datasets were further analyzed separately. Various environmental variables were tested as explanatory variables in CCA analyses. Year (for all datasets) and plot identifier in interaction with year (only for datasets roots and trees) were used as covariates. We used hierarchical split-plot design for permutations and permutated whole plots (i.e., keeping the data from one plot together) for testing the effect of management. These whole plots were permuted at random within the covariate year, and the two split plots (trees) were not permuted. The effect of tree species was tested by split-plot permutations (within the whole plot, i.e., the data from different tree species were permutated within the plot). The effects of explanatory variables were summarized in partial analyses after removing the effects of covariates. The Monte Carlo permutation test was used and only significant variables were selected.
The effects of tree species and management on the number of ECM fungal species for tree were analyzed by linear mixed-effects model fit by REML using the lme function in nlme package in R 3.6.2., with plot as a random effect. The effect of plot was tested by comparing the model without plot (made by gls function) and the model with plot (made by lme function) by function ANOVA.
The numbers of shared and specific ECM fungal species in NEX and EXT treatments on P. abies, L. decidua, and B. pendula were presented by Venn diagrams made in web application BioVenn (Hulsen et al., 2008).
Two root samples (19Pa_2b and 2Ld_1a) and one tree sample (9Pa) were excluded from analysis as outliers.
Results
ECM Fungal Species and Potential ECM Networks
In total, 57 ECM taxa were detected (Supplementary Table 2): 51 taxa were assigned at the species level, two at the genus level, and four at the order level. There were no significant differences in the number of ECM fungal species per tree between the NEX and EXT treatments (p > 0.05, Table 1). However, differences among tree species were detected; silver birch had significantly higher number of ECM fungal species per tree than Norway spruce and European larch (p < 0.005), and Norway spruce and European larch did not differ significantly (p > 0.05, Table 1).
Table 1. Mean number of ectomycorrhizal (ECM) fungal species per tree, mean tree height, and mean tree perimeter in NEX (plots without management) and EXT (plots with traditional management) treatments and in Norway spruce, European larch, and silver birch.
The tree heights and perimeters within a treatment were significantly greater for silver birch than for the other tree species (p < 0.001, Table 1). In EXT, tree heights (p < 0.001, Table 1) and perimeters (p < 0.005, Table 1) were significantly greater than in NEX.
Fourteen ECM taxa were found in all tree species, five species were detected in Norway spruce and European larch, seven in Norway spruce and silver birch, five in European larch and silver birch, and eight were unique to Norway spruce, 10 to European larch, and 14 to silver birch. When comparing the numbers of ECM fungal species in single tree species between the NEX and EXT treatments (Figure 1), only Norway spruce had a lower number of ECM fungal species in EXT (six species) compared with NEX (11 species). When focusing on the potential mycorrhizal networking among trees in individual plots, we found no significant differences between NEX and EXT in the numbers of ECM fungal species that were shared by all host trees, between individual pairs of host trees, or associated with a single host tree only.
Figure 1. Venn diagram displaying the numbers of shared and specific ectomycorrhizal (ECM) fungal species in NEX (plots without management) and EXT (plots with traditional management) treatments on (A) Picea abies, (B) Larix decidua, and (C) Betula pendula.
Overall, 12 ECM species belonged to contact exploration type, 15 species to short-distance type, eight species to medium-distance smooth type, 12 species to medium-distance fringe type, and 11 species to long-distance type (Supplementary Table 2). When considering only the numbers of rhizomorph-forming ECM fungal species, i.e., medium-distance and long-distance exploration types, more ECM species shared by all tree species were found in NEX compared with EXT, but these differences were not significant (Supplementary Table 3). In case of NEX treatment, the potential networkers were Cortinarius bataillei, C. croceus, Paxillus involutus, Piloderma olivaceum, and Thelephora terrestris, and in case of EXT, the potential networkers were P. olivaceum and T. terrestris.
There were also no significant differences in the numbers of shared species among combinations of two tree species or in single tree species between NEX and EXT treatments.
Five ECM fungal species dominated in both treatments: Cenococcum geophilum, Lactarius rufus, P. involutus, P. olivaceum, and T. terrestris (Table 2). Most of these species were detected in all tree species and together accounted for 50.1% in the NEX plots and 51.7% in the EXT plots. The potential of these species to form mycorrhizal networks differed only for the case of P. involutus, which was absent in EXT on Norway spruce.
Table 2. Dominant ECM fungal species present in NEX (plots without management) and EXT (plots with traditional management) treatments.
ECM fungal species that were unique to individual tree species were generally rare. For Norway spruce, Piloderma bicolor (0.5%, both NEX and EXT) was present in both treatments; in NEX, Amphinema byssoides (1.8%), Helotiales1 (0.8%), and Imleria badia (0.5%) were present; in EXT, Tylopilus felleus (0.8%) was present. For European larch, Suillus gravillei (in NEX 0.8% and in EXT 2.6%) was present in both treatments; in NEX, Cortinarius cinnamomeus (0.8%) was present; in EXT, C. parvanulatus (0.5%), S. cavipes (0.3%), and S. viscidus (1.8%) were present. For silver birch, Leccinum holopus (1.5% in both NEX and EXT), L. variicolor (in NEX 3.5% and in EXT 3.3%), and Meliniomyces bicolor (in NEX 0.5% and in EXT 0.8%) were present in both treatments; in NEX, C. porphyropus (0.5%), Inocybe soluta (0.8%), Laccaria proxima (0.5%), and L. glyciosmus (2.3%) were present; in EXT, Amanita muscaria (0.5%) and Clavulina sp. (0.5%) were present.
Management Effect on ECM Fungal Species Composition
Management significantly influenced the ECM fungal species composition in the studied plots and explained 11.7% of the detected variability (Table 3). The other studied environmental variables (vegetation and dead wood cover, geographic position, morphometric parameters of host trees—height, perimeter, and root placement in soil) had no significant effect.
Table 3. Results of canonical correspondence analyses (CCA) of ECM fungal community.
Regarding the community composition of ECM fungal species based on the tree dataset, the only significant variables were tree species and management type (Table 3 and Figure 2). The morphometric parameters of trees (height and perimeter) had no effect on the ECM fungal species composition. When the tree species were analyzed separately, only the ECM fungal species compositions of Norway spruce and European larch were affected by the management type (Table 3 and Figures 3A,B). In Norway spruce trees, C. croceus, Meliniomyces vraolstadiae, P. involutus, T. terrestris, and T. fibrillosa were more abundant in NEX, whereas A. byssoides, L. rufus, and T. asterophora were more abundant in EXT. In European larch trees, C. bataillei, L. laccata, P. olivaceum, and T. terrestris were more abundant in NEX and Suillus spp. in EXT. The ECM fungal species composition of silver birch was not influenced by the studied variables (Table 3 and Figure 3C). The ECM fungal species composition, based on the root sample dataset, was significantly affected by tree species, management type, and soil horizon (Table 3 and Figure 4).
Figure 2. Ordination diagram from canonical correspondence analysis (CCA) of ECM fungal species presence on tree species (P—Picea abies, L—Larix decidua, and B—Betula pendula) in NEX (plots without management) and EXT (plots with traditional management) treatments. Species: AmaMus—Amanita muscaria, AmpBys—Amphinema byssoides, CenGeo—Cenococcum geophilum, ClaSp—Clavulina sp., CorAno—Cortinarius anomalu, CorBat—Cortinarius bataillei, CorCin—Cortinarius cinnamomeus, CorCro—Cortinarius croceus, CorPar—Cortinarius parvannulatus, CorPor—Cortinarius porphyropus, CorRap—Cortinarius raphanoides, CorVit—Cortinarius vitiosus, Hel1—Helotiales1, Hel2—Helotiales2, Hel3—Helotiales3, ImlBad—Imleria badia, InoSol—Inocybe soluta, LacLac—Laccaria laccata, LacPro—Laccaria proxima, LacGly—Lactarius glyciosmus, LacHel—Lactarius helvus, LacLig—Lactarius lignyotus, LacNec—Lactarius necator, LacRuf—Lactarius rufus, LacTab—Lactarius tabidus, LecHol—Leccinum holopus, LecVar—Leccinum variicolor, MelBic—Meliniomyces bicolor, MelVar—Meliniomyces variabilis, MelVra—Meliniomyces vraolstadiae, NeoLur—Neoboletus luridiformis, PaxInv—Paxillus involutus, PhiFor—Phialocephala fortinii, PhiSph—Phialocephala sphaeroides, PilBic—Piloderma bicolor, PilOli—Piloderma olivaceum, PilSph—Piloderma sphaerosporum, PseFla—Pseudotomentella flavovirens, PseVep—Pseudotomentella vepallidospora, RhiSal—Rhizopogon salebrosus, RhoBut—Rhodocollybia butyracea, RusAtr—Russula atrorubens, RusBet—Russula betularum, RusEme—Russula emetica, RusOch—Russula ochroleuca, RusSph—Russula sphagnophila, Seb—Sebacinales, SuiCav—Suillus cavipes, SuiGra—Suillus grevillei, SuiVis—Suillus viscidus, TheSp—Thelephora sp., TheTer—Thelephora terrestris, TomBad—Tomentella badia, TomSub—Tomentella sublilacina, TylFel—Tylopilus felleus, TylAst—Tylospora asterophora, and TylFib—Tylospora fibrillosa.
Figure 3. Ordination diagram from detrended correspondence analysis (DCA) with the ECM fungal species composition in NEX (plots without management) and EXT (plots with traditional management) treatments for single tree species: (A) Picea abies, (B) Larix decidua, and (C) Betula pendula.
Figure 4. Ordination diagram from canonical correspondence analysis (CCA) of ECM fungal species presence on tree roots in NEX (plots without management) and EXT (plots with traditional management) treatments for single tree species (P—Picea abies, L—Larix decidua, and B—Betula pendula) differentiated according to organic (F, H) and organo-mineral (A) soil horizons. Species: AmaMus—Amanita muscaria, AmpBys—Amphinema byssoides, CenGeo—Cenococcum geophilum, ClaSp—Clavulina sp., CorAno—Cortinarius anomalu, CorBat—Cortinarius bataillei, CorCin—Cortinarius cinnamomeus, CorCro—Cortinarius croceus, CorPar—Cortinarius parvannulatus, CorPor—Cortinarius porphyropus, CorRap—Cortinarius raphanoides, CorVit—Cortinarius vitiosus, Hel1—Helotiales1, Hel2—Helotiales2, Hel3—Helotiales3, ImlBad—Imleria badia, InoSol—Inocybe soluta, LacLac—Laccaria laccata, LacPro—Laccaria proxima, LacGly—Lactarius glyciosmus, LacHel—Lactarius helvus, LacLig—Lactarius lignyotus, LacNec—Lactarius necator, LacRuf—Lactarius rufus, LacTab—Lactarius tabidus, LecHol—Leccinum holopus, LecVar—Leccinum variicolor, MelBic—Meliniomyces bicolor, MelVar—Meliniomyces variabilis, MelVra—Meliniomyces vraolstadiae, NeoLur—Neoboletus luridiformis, PaxInv—Paxillus involutus, PhiFor—Phialocephala fortinii, PhiSph—Phialocephala sphaeroides, PilBic—Piloderma bicolor, PilOli—Piloderma olivaceum, PilSph—Piloderma sphaerosporum, PseFla—Pseudotomentella flavovirens, PseVep—Pseudotomentella vepallidospora, RhiSal—Rhizopogon salebrosus, RhoBut—Rhodocollybia butyracea, RusAtr—Russula atrorubens, RusBet—Russula betularum, RusEme—Russula emetica, RusOch—Russula ochroleuca, RusSph—Russula sphagnophila, Seb—Sebacinales, SuiCav—Suillus cavipes, SuiGra—Suillus grevillei, SuiVis—Suillus viscidus, TheSp—Thelephora sp., TheTer—Thelephora terrestris, TomBad—Tomentella badia, TomSub—Tomentella sublilacina, TylFel—Tylopilus felleus, TylAst—Tylospora asterophora, and TylFib—Tylospora fibrillosa.
Discussion
Potential Mycorrhizal Networks and ECM Fungal Diversity
Our study on the effect of different management practices on ECM fungal communities of regenerating trees in the post-windstorm area in the Tatra Mts. revealed no effect on dominant ECM fungi and their potential to form mycorrhizal networks. For both treatment types, five ECM fungal species with broad host range dominated and included both early-colonizers (C. geophilum, P. involutus, and T. terrestris) and late-successional ECM fungi (L. rufus and P. olivaceum; Visser, 1995; Twieg et al., 2007; Kim et al., 2021). The windstorm apparently affected forest ecosystem in such an extreme manner that the management type produced only insignificant adjustments and had only limited influence on the potential to form ECM networks.
In contrast to our hypothesis, we also found almost no difference in the total number of potentially shared ECM fungal species among different tree species between the traditionally managed plots (EXT) and plots without management (NEX). Even if we focus on ECM species with exploration types with rhizomorphs, we found no significant differences, although the numbers of species shared by all tree species slightly differ in favor of NEX. Exploration types with rhizomorphs are more carbon demanding compared with contact and short-distance types (Defrenne et al., 2019; Wasyliw and Karst, 2020); on the other hand, they are able to form common mycelial networks and exploit water and nutrients from a larger area (Simard et al., 2003; Bakker et al., 2006). So, a higher number of these species per plot in NEX could be explained by lower damage to the mycelium caused by management practices, since hyphae are able to emanate from dying or recently dead root tips for a short time (Simard, 2009).
Norway spruce was the only tree species, which showed a higher number of specific ECM fungal species in NEX than in EXT (Figure 1). This difference is also evident from rarefaction curves (Supplementary Figure 2). In contrast to other tree species in our study, Norway spruce is a shade-tolerant species, which is able to regenerate under mature forest canopies (Metslaid et al., 2007). Moreover, Norway spruce requires an adequate water supply. Gömöryová et al. (2008) detected lower soil moisture levels 2 years after the windstorm in EXT compared with NEX, which could, along with the increased insolation following tree removal, negatively affect spruce fitness. The better Norway spruce seedling fitness in NEX (although not supported by morphological characteristics) combined with the connectivity of the original forest could explain the differences observed. In contrast, both silver birch and European larch possessed the same numbers of specific ECM symbionts in both treatments. The highest ECM fungal species richness in silver birch can be the consequence of its larger dimensions related to its pioneer strategy and higher adaptability to clear-cut sites (Dubois et al., 2020). Similar findings were also reported by Twieg et al. (2007), who detected more ECM fungal species on paper birch than Douglas fir in young stands after clear-cut logging. However, in that case, the higher species richness of paper birch was also attributed to the presence of birch stumps.
No or small differences in quantitative parameters between treatments could be also related to a common problem in evaluating the responses of fungal communities to environmental factors or management. Generally, most fungal communities are composed of a small number of dominant species with a wide ecological valence that only respond to very pronounced gradients. These species are accompanied by few subdominant species and numerous rare species. Due to their rarity and the patchiness of ECM fungi dropping after 3–4 m (Lilleskov et al., 2004; Pickles et al., 2012), it is impossible to disentangle the random occurrence of these species and their response to environmental factors without including many more research plots and more intensive root sampling. Rosinger et al. (2018) found the majority of ECM taxa in one or two plots out of approximately 100 plots. While the overall fungal diversity in temperate forests is high, the species diversity at small scales is low and does not usually exceed 10–20 species (Tedersoo et al., 2003; Anderson et al., 2014). Thus, the theoretical absence of a few species for a total area of 10 times 125 m2 could indicate larger losses at larger scales. However, this hypothesis must be verified by further research focusing on the ecology of rare ECM fungal species, in which all the above mentioned limitations would be taken into consideration (i.e., rarity and patchiness of ECM fungal species and sampling effort). The combination of fruit body monitoring together with environmental soil sequencing could be a promising nondestructive method, but unfortunately without data on the associations with host trees.
Effect of Management on ECM Fungal Species Composition 15 Years After a Windstorm
Interestingly, management is the only significant variable that explains ECM fungal species composition 15 years after a windstorm on the studied plots and for Norway spruce and European larch trees. The similarity of the ECM fungal communities did not reflect the phylogenetic relationships of the studied trees, which contradicts to the observations of Rog et al. (2020).
The observed differences in ECM fungi composition between NEX and EXT in Norway spruce are difficult to explain because there are no clear patterns in exploration types (Agerer, 2001) or spore characteristics (Bernicchia and Gorjón, 2010; Knudsen and Vesterholt, 2012). We also do not expect differences in ability to spread, because both treatments contain species with pileate and resupinate fruit bodies (Lilleskov and Bruns, 2005; Halbwachs et al., 2016). However, the affinity of T. fibrillosa and T. terrestris for NEX and of A. byssoides and T. asterophora for EXT has been previously observed in 4–9-year-old seedlings from a similar area in the Tatra Mts. (Vašutová et al., 2018), and these differences persist in Norway spruce between treatments. These differences are difficult to evaluate without detailed knowledge of the ECM fungal ecology and physiology. European larch trees are often occupied by C. bataillei, L. laccata, P. olivaceum, and T. terrestris in NEX and by Suillus spp. in EXT. Coleman et al. (1989) found several Suillus spp. to be drought tolerant, and in contrast, L. laccata was not drought tolerant, which could partly explain the observed differences.
Although the studied sites were chosen to be the same as each other and the lack of differences in silver birch-associated ECM fungal community supported a good experimental design, we could not fully exclude the differences in ECM fungal community of original Norway spruce forest with European larch and Scotch pine admixtures (Lariceto–Picetum) that persisted in spores and as mycorrhizae on surviving seedlings. It would be interesting to determine whether these differences will persist 10–15 years later when ECM fungal diversity should be saturated after canopy closure (Visser, 1995; Twieg et al., 2007). More replications under different local conditions and with other tree species are necessary to generalize our results.
Regarding the effects of vertical distribution on the ECM fungal species composition, Tylospora spp., Amphinema spp., and Meliniomyces spp. showed an affinity for an organic layer, and Tomentella and Pseudotomentella spp. showed an affinity for an organo-mineral layer (Figure 4), as was previously observed by Tedersoo et al. (2003). They detected athelioid and thelephoroid species as indicators of organic layers, but Tomentella spp. were present in all horizons depending on the species. These findings support the importance of vertical niche partitioning, as was already discussed by Bahram et al. (2015).
In conclusion, different post-windstorm management practices do not significantly influence the potential for mycorrhizal networking among Norway spruce, European larch, and silver birch, which is mostly provided by common ECM fungal species with broad ecological niches. Management type affects the ECM fungal species compositions on Norway spruce and European larch; however, the consequences are difficult to assess without knowledge of the ecology and physiology of the affected ECM fungi. Norway spruce in NEX harbored a higher number of unique species than in EXT, but extremely intensive sampling is required to evaluate the possible loss of rare species due to traditional management.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
PV, MV, FH, PC, and PF designed the research. PV, MV, and FH performed the research and collected the data. ME-J analyzed the data. PV, MV, and ME-J interpreted the data. PV and MV wrote the manuscript with the assistance from all co-authors. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Ministry of Education, Youth and Sports of CR within the CzeCOS program, grant number LM2018123.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Zdeňka Kovaříková, Vladana Procházková, Jan Purkyt, and students from the Faculty of Forestry and Wood Technology, Mendel University in Brno for field and lab assistance.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.641232/full#supplementary-material
Supplementary Figure 1 | Map of the research area, Tatra Mountains National Park (TANAP), Slovakia (SK). Two management treatments, site without management (NEX) and site with traditional management (EXT), each with ten research plots.
Supplementary Figure 2 | Rarefaction curves showing the ECM species richness on tree species (P – Picea abies, L – Larix decidua, and B – Betula pendula) in NEX (plots without management) and EXT (plots with traditional management) treatments.
Supplementary Table 1 | Characteristics of the plots in NEX (plots without management) and EXT (plots with traditional management) treatments (geographic position, cover of trees, cover of dominant herbs and shrubs, and dead wood cover); characteristics of the selected trees (height and perimeter); and characteristics of selected roots (depth, soil horizon, distance from the tree stem, and direction to other tree species).
Supplementary Table 2 | List of ECM taxa, including exploration of each taxa, detected on root samples of the studied tree species (Picea abies, Larix decidua, and Betula pendula), including the numbers of root samples with taxon presence divided into NEX (plots without management) and EXT (plots with traditional management) treatments. For each taxon, the GenBank Accession Number, closest hit and percentage of similarity are shown.
Supplementary Table 3 | Presence/absence matrix showing rhizomorph-forming ECM fungal species in the plots in NEX (plots without management) and EXT (plots with traditional management) treatments.
References
Abarenkov, K., Nilsson, H. R., Larsson, K. H., Alexander, I. J., Eberhardt, U., Erland, S., et al. (2010). The UNITE database for molecular identification of fungi – recent updates and future perspectives. New Phytol. 186, 281–285. doi: 10.1111/j.1469-8137.2009.03160.x
Agerer, R. (1991). “Characeterisation of ectomycorrhizae,” in Techniques for the study of mycorrhiza, eds J. R. Norris, D. J. Read, and A. K. Varma, (London: Academic Press), 25–73.
Agerer, R., and Rambold, G. (2004–2015). DEEMY – an Information System for Characterization and Determination of Ectomycorrhizae. Available online at: http://www.deemy.de (accessed November 20, 2019).
Anderson, I. C., Genney, D. R., and Alexander, I. J. (2014). Fine-scale diversity and distribution of ectomycorrhizal fungal mycelium in a scots pine forest. New Phytol. 201, 1423–1430. doi: 10.1111/nph.12637
Bahram, M., Peay, K. G., and Tedersoo, L. (2015). Local-scale biogeography and spatiotemporal variability in communities of mycorrhizal fungi. New Phytol. 205, 1454–1463. doi: 10.1111/nph.13206
Bahram, M., Põlme, S., Kõljalg, U., and Tedersoo, L. (2011). A single European aspen (Populus tremula) tree individual may potentially harbour dozens of Cenococcum geophilum ITS genotypes and hundreds of species of ectomycorrhizal fungi. FEMS Microbiol. Ecol. 75, 313–320. doi: 10.1111/j.1574-6941.2010.01000.x
Bakker, M. R., Augusto, L., and Achat, D. L. (2006). Fine root distribution of trees and understory in mature stands of maritime pine (Pinus pinaster) on dry and humid sites. Plant Soil 286, 37–51. doi: 10.1007/s11104-006-9024-4
Baldrian, P. (2009). Ectomycorrhizal fungi and their enzymes in soils: is there enough evidence for their role as facultative soil saprotrophs? Oecologia 161, 657–660. doi: 10.1007/s00442-009-1433-7
Benson, D. A., Cavanaugh, M., Clark, K., Karsch-Mizrachi, I., Lipman, D. J., Ostell, J., et al. (2007). GenBank. Nucleic Acids Res. 41, 36–41.
Bernicchia, A., and Gorjón, S. P. (2010). Fungi Europaei, 12: Corticiaceae s.l. Alassio: Alassio Candusso, 1008.
Coleman, M. D., Bledsoe, C. S., and Lopushinsky, W. (1989). Pure culture response of ectomycorrhizal fungi to imposed water stress. Can. J. Bot. 67, 29–39. doi: 10.1139/b89-005
Defrenne, C. E., Philpott, T. J., Guichon, S. H. A., Roach, W. J., Pickles, B. J., and Simard, S. W. (2019). Shifts in ectomycorrhizal fungal communities and exploration types relate to the environment and fine-root traits across interior Douglas-fir forests of western Canada. Front. Plant Sci. 10:643. doi: 10.3389/fpls.2019.00643
Dietz, L., Collet, C., Dupouey, J., Lacombe, E., Laurent, L., and Gégout, J. (2020). Windstorm-induced canopy oppenings accelerate temperate forest adaptation to global warming. Glob. Ecol. Biogeogr. 29, 2067–2077.doi10.1111/geb.13177
Don, A., Bärwolff, M., Kalbitz, K., Andruschkewitsch, R., Jungkunst, H. F., and Schulze, E. (2012). No rapid soil carbon loss after a windthrow event in the High Tatra. For. Ecol. Manag. 276, 239–246. doi: 10.1016/j.foreco.2012.04.010
Dubois, H., Verkasalo, E., and Claessens, H. (2020). Potential of birch (Betula pendula roth and B. pubescens Ehrh.) for forestry and forest-based industry sector within the changing climatic and socio-economic context of Western Europe. Forests 11:336. doi: 10.3390/f11030336
Egli, S., Peter, M., and Falcato, S. (2002). Dynamics of ectomycorrhizal fungi after windthrow. For. Snow Landsc. Res. 77, 81–88.
Fischer, A., Marshall, P., and Camp, A. (2013). Disturbances in deciduous forest ecosystems of the northern hemisphere: their effects on both recent and future forest development. Biodivers. Conserv. 22, 1863–1893. doi: 10.1007/s10531-013-0525-1
Fleischer, P. (2008). Windfall research and monitoring in the High Tatra Mts, objectives, principles, methods and current status. Contribut. Geophys. Geod. 38, 233–248.
Ford, S. A., Kleinman, J. S., and Hart, J. L. (2018). Effects of wind disturbance and salvage harvesting on macrofungal communities in a Pinus Woodland. Forest Ecol. Manag. 407, 31–46. doi: 10.1016/j.foreco.2017.10.010
Gardes, M., and Bruns, T. D. (1993). ITS primers with enhanced specificity for basidiomycetes: application to the identification of mycorrhizae and rusts. Mol. Ecol. 2, 113–118. doi: 10.1111/j.1365-294x.1993.tb00005.x
Gömöryová, E., Střelcová, K., Škvarenina, J., Bebej, J., and Gömöry, D. (2008). The impact of windthrow and fire disturbances on selected soil properties in the Tatranational park. Soil Water Res. 3, S74–S80. doi: 10.17221/9/2008-SWR
Gorzelak, M. A., Asay, A. K., Pickles, B. J., and Simard, S. W. (2015). Inter-plant communication through mycorrhizal networks mediates complexadaptive behaviour in plant communities. AoB Plants 7:lv050. doi: 10.1093/aobpla/plv050
Gregow, H., Laaksonen, A., and Alper, M. E. (2017). Increasing large scale windstorm damage in Western, Central and Northern European forets, 1951-2010. Sci. Rep. 7:46397.
Guignabert, A., Delerue, F., Gonzalez, M., Augusto, L., and Bakker, M. R. (2018). Effects of management practices and topography on ectomycorrhizal fungi of maritime pine during seedling recruitment. Forests 9:245. doi: 10.3390/f9050245
Haas, J. C., Street, N. R., Sjödin, A., Lee, N. M., Högberg, M. N., Näsholm, A., et al. (2018). Microbial community response to growing season and plant nutrient optimisation in a boreal Norway spruce forest. Soil Biol. Biochem. 125, 197–209. doi: 10.1016/j.soilbio.2018.07.005
Halbwachs, H., Simmel, J., and Bässler, C. (2016). Tales and mysteries of fungal fruiting: How morphological and physiological traits affect a pileate lifestyle. Fungal Biol. Rev. 30, 36–61. doi: 10.1016/j.fbr.2016.04.002
Hobbie, E. A., Colpeart, J. V., White, M. W., Ouimette, A. P., and Macko, S. A. (2008). Nitrogen form, availability and mycorrhizal colonization affect biomass and nitrogen isotope patterns in Pinus sylvestris. Plant Soil 310, 121–136. doi: 10.1007/s11104-008-9637-x
Horton, T. R., and Bruns, T. D. (1998). Multiple-host fungi are the most frequent and abundant ectomycorrhizal types in a mixed stand of Douglas fir (Pseudotsuga menziesii) and bishop pine (Pinus muricata). New Phytol. 139, 331–339. doi: 10.1046/j.1469-8137.1998.00185.x
Hulsen, T., de Vlieg, J., and Alkema, W. (2008). BioVenn – a web application for the comparison and visualization of biological lists using area-proportional Venn diagrams. BMC Genomics 9:488. doi: 10.1186/1471-2164-9-488
Jones, M. D., Durall, D. M., and Cairney, J. W. G. (2003). Ectomycorrhizal fungal communities in young forest stands regenerating after clearcut logging. New Phytol. 157, 399–422. doi: 10.1046/j.1469-8137.2003.00698.x
Karst, J., Randall, M. J., and Gehring, C. A. (2014). Consequences for ectomycorrhizal fungi of the selective loss or gain of pine across landscape. Botany 92, 855–865. doi: 10.1139/cjb-2014-0063
Kim, S., Axelsson, E. P., Girona, M. M., and Senior, J. K. (2021). Continuous-cover forestry maintains soil fungal communities in Norway spruce dominated boreal forests. For. Ecol. Manag. 480:118659. doi: 10.1016/j.foreco.2020.118659
Klein, T., Siegwolf, R. T. W., and Körner, C. (2016). Belowground carbon trade among tall trees in a temperate forest. Forest Ecol. 352, 342–344. doi: 10.1126/science.aad6188
Knudsen, H., and Vesterholt, J. (2012). Funga Nordica: Agaricoid, Boletoid, Clavarioid, Cyphelloid and Gastroid Genera. Copenhagen: Nordsvamp, 1083.
Kranabetter, J. M., Friesen, J., Gamiet, S., and Kroeger, P. (2004). Ectomycorrhizal mushroom distribution by stand age in western hemlock – lodgepole pine forests of northwestern British Columbia. Can. J. For. Res. 35, 1527–1539. doi: 10.1139/x05-095
Lang, C., Seven, J., and Polle, A. (2011). Host preferences and differential contributions of deciduous tree species shape Mycorrhizal species richness in a mixed Centra European forest. Mycorrhiza 21, 297–308. doi: 10.1007/s00572-010-0338-y
Lilleskov, E. A., and Bruns, T. D. (2005). Spore dispesal of a resupinate ectomycorrhizal fungus, Tomentella sublilacina, via soil food webs. Mycologia 97, 762–769. doi: 10.3852/mycologia.97.4.762
Lilleskov, E. A., Bruns, T. D., Horton, T. R., Taylor, D. L., and Grogan, P. (2004). Detection of forest stand-level spatial structure in ectomycorrhizal fungal communities. FEMS Microbiol. Ecol. 49, 319–332. doi: 10.1016/j.femsec.2004.04.004
Lofgren, L., Nguyen, N. H., and Kennedy, P. G. (2018). Ectomycorrhizal host specifity in a changing world: can legacy effects explain anomalous current associations? New Phytol. 220, 1273–1284. doi: 10.1111/nph.15008
Martin, F., Kohler, A., Murat, C., Veneault-Fourrey, C., and Hibbett, D. S. (2016). Unearthing the roots of ectomycorrhizal symbioses. Nat. Rev. Microbiol. 14, 760–773. doi: 10.1038/nrmicro.2016.149
McGlinn, D. J., Xiao, X., May, F., Gotelli, N. J., Engel, T., Blowes, S. A., et al. (2018). Measurement of biodiversity (MoB): a method to separate the scale-dependent effects of species abundance distribution, density, and aggregation on diversity change. Methods Ecol. Evol. 10, 258–269. doi: 10.1111/2041-210x.13102
Metslaid, M., Jõgiste, K., Nikinmaa, E., Moser, W. K., and Porcar-Castell, A. (2007). Tree variables related to growth response and acclimation of advance regeneration of Norway spruce and other coniferous species after release. Forest Ecol. Manag. 250, 56–63. doi: 10.1016/j.foreco.2007.03.009
Molina, R., and Horton, T. (2015). “Mycorrhiza specificity: its role in the development and function of common mycelial networks,” in Mycorrhizal Networks. Ecological Studies (Analysis and Synthesis), ed. T. Horton, (Dordrech: Springer), 1–39. doi: 10.1007/978-94-017-7395-9_1
Molina, R., Massicotte, H. B., and Trappe, J. M. (1992). “Specificity phenomena in mycorrhizal symbioses: community-ecological consequences and practical implication,” in Mycorrhizal Functioning: An Integrated Plant Fungal Process, ed. M. F. Allen, (London: Chapman and Hall), 357–423.
Nikolcheva, L. G., and Bärlocher, F. (2004). Taxon-specific fungal primers reveal unexpectedly high diversity during leaf decomposition in a stream. Mycol. Prog. 3, 41–49. doi: 10.1007/s11557-006-0075-y
Nilsson, R. H., Anslan, S., Bahram, M., Wurzbacher, C., Baldrian, P., and Tedersoo, L. (2019). Mycobiome diversity: high-throughput sequencing and identification of fungi. Nat. Rev. Microbiol. 17, 95–109. doi: 10.1038/s41579-018-0116-y
Pec, G. J., Simard, S. W., Cahill, J. F., and Karst, J. (2020). The effects of ectomycorrhizal fungal networks on seedling establishment are contingent on species and severity of overstorey mortality. Mycorrhiza 30, 173–183. doi: 10.1007/s00572-020-00940-4
Pickles, B. J., Genney, D. R., Anderson, I. C., and Alexander, I. J. (2012). Spatial analysis of ectomycorrhizal fungi reveals that root tip communities are structured by competitive interactions. Mol. Ecol. 21, 5110–5123. doi: 10.1111/j.1365-294X.2012.05739.x
R Development Core Team, (2019). R: A Language and Environment for Statistical Computing. Vienna: The R Foundation for Statistical Computing.
Rog, I., Rosenstock, N. P., Körner, C., and Klein, T. (2020). Share the wealth: trees with greater ectomycorrhizal species overlap share more carbon. Mol. Ecol. 29, 2321–2333. doi: 10.1111/mec.15351
Rosinger, C., Sandén, H., Matthews, D., Mayer, M., and Godbold, D. L. (2018). Patterns in ectomycorrhizal diversity, community composition, and exploration types in European beech, pine, and spruce forests. Forests 9:445. doi: 10.3390/f9080445
Simard, S. (2018). “Mycorrhizal networks facilitate tree communication, learning, and memory,” in Memory and Learning in Plant. Signaling and Communication in Plants, eds F. Baluska, M. Gagliano, and G. Witzany, (Cham: Springer), 191–213. doi: 10.1007/978-3-319-75596-0_10
Simard, S. W. (2009). The foundaitional role of mycorrhizal networks in self-organization of interioir Douglaska-fir forests. For. Ecol. Manag. 258, 95–107. doi: 10.1002/9780470724460.ch8
Simard, S. W., Beiler, K. J., Bingham, M. A., Deslippe, J. R., Philip, L. J., and Teste, F. J. (2012). Mycorrhizal networks: Mechanisms, ecology and modelling. Fungal Biol. Rev. 26, 39–60. doi: 10.1016/j.fbr.2012.01.001
Simard, S. W., Jones, M. D., and Durall, D. M. (2003). “Carbon and nutrient fluxes within and between mycorrhizal plants,” in Mycorrhizal Ecology. Ecological Stusies (Analysis and Synthesis), eds M. G. A. van der Heijden, I and R. Sanders, (Berlin: Springer), 33–74. doi: 10.1007/978-3-540-38364-2_2
Simard, S. W., Perry, D. A., Jones, M. D., Myrold, D. D., Durall, D. M., and Molina, R. (1997). Net transfer of carbon between ectomyycorrhizal tree species in the field. Nature 388, 579–582. doi: 10.1038/41557
Song, Y. Y., Simard, S. W., Carroll, A., Mohn, W. W., and Zeng, R. S. (2015). Defoliation of interior Douglaska-fir elicits carbon transfer and stress signalling to ponderosa pine neighbors through ectomycorrhizal networks. Sci. Rep. 5:8495.
Taberlet, P., Gielly, L., Pautou, G., and Bouvet, J. (1991). Universal primer for amplification of three non-coding regions chloroplast DNA. Plant Mol. Biol. 17, 1105–1109. doi: 10.1007/bf00037152
Taudiere, A., Munoz, F., Lesne, A., Monnet, A., Bellanger, J., Selosse, M., et al. (2015). Beyond ectomycorrhizal bipartite networks: projected networks demonstrate contrasted patterns between early- and late-successional plants in Corsica. Front. Plant Sci. 6:881. doi: 10.3389/fpls.2015.00881
Tedersoo, L., Kõljalg, U., Hallenberg, N., and Larsson, K. (2003). Fine scale distribution of ectomycorrhizal fungi and roots across substrate layers including coarse woody debris in a mixed forest. New Phytol. 159, 153–165. doi: 10.1046/j.0028-646x.2003.00792.x
ter Braak, C. J. F., and Šmilauer, P. (2012). Canoco Reference Manual and User’s Guide: Software for Ordination (version 5.0). New York, NY: Microcomputer Power.
Treu, R., Karst, J., Randall, M., Pec, G. J., Cigan, P. W., Simard, S. W., et al. (2014). Decline of ectomycorrhizal fungi following a mountain pine beetle epidemic. Ecology 95, 1096–1103. doi: 10.1890/13-1233.1
Twieg, B. D., Durall, D. M., and Simard, S. W. (2007). Ectomycorrhizal fungal succession in mixed temperate forests. New Phytol. 176, 437–447. doi: 10.1111/j.1469-8137.2007.02173.x
van der Heiden, M. G. A., and Horton, T. R. (2009). Socialism in soil? The importance of mycorrhizal fungal networks for facilitation in natural ecosystems. J. Ecol. 97, 1139–1150. doi: 10.1111/j.1365-2745.2009.01570.x
van der Linde, S., Suz, L. M., Orme, C. D. L., Cox, F., Andreae, H., Asi, E., et al. (2018). Environment and host as large-scale controls of ectomycorrhizal fungi. Nature 558, 243–248. doi: 10.1038/s41586-018-0189-9
Vašutová, M., Edwars-Jonášová, M., Veselá, P., Effenberková, L., Fleischer, P., and Cudlín, P. (2018). Management regime is the most important factor influencing ectomycorrhizal species community in Norway spruce forests after windthrow. Mycorrhiza 28, 221–233. doi: 10.1007/s00572-018-0820-5
Velmala, S. M., Rajala, T., Heinonsalo, J., Taylor, A. F. S., and Pennanen, T. (2014). Profiling functions of ectomycorrhizal diversity and root structuring in seedling of Norway spruce (Picea abies) with fast- and slow-growing phenotypes. New Phytol. 201, 610–622. doi: 10.1111/nph.12542
Veselá, P., Vašutová, M., Hofmannová, K., Edwards-jonášová, M., and Cudlín, P. (2019b). Ectomycorrhizal community on Norway spruce seedlings following bark beetle infestation. Forest 10:740. doi: 10.3390/f10090740
Visser, S. (1995). Ectomycorrhizal fungal succession in jack pine stands following wildfire. New Phytol. 129, 389–401. doi: 10.1111/j.1469-8137.1995.tb04309.x
Vu, D., Groenewald, M., De Vries, M., Gehrmann, T., Stielow, B., Eberhardt, U., et al. (2019). Large-scale generation and analysis of filamentous fungal DNA barcodes boosts coverage for kingdom fungi and reveals thresholds for fungal species and higher taxon delimitation. Stud. Mycol. 92, 135–154. doi: 10.1016/j.simyco.2018.05.001
Wasyliw, J., and Karst, J. (2020). Shifts in ectomycorrhizal exploration types parallel leaf and fine root area with forest age. J. Ecol 108, 2270–2282. doi: 10.1111/1365-2745.13484
White, T. J., Bruns, T., Lee, S., and Taylor, J. (1990). “Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics,” in PCR Protocols: A Guide to Methods and Applications, eds N. Innis, D. Gelfand, J. Sninsky, and T. White, (New York, NY: Academic Press, Inc), 315–322. doi: 10.1016/b978-0-12-372180-8.50042-1
Keywords: disturbances, diversity, ectomycorrhizal fungi, exploration types, mycorrhizal networks
Citation: Veselá P, Vašutová M, Edwards-Jonášová M, Holub F, Fleischer P and Cudlín P (2021) Management After Windstorm Affects the Composition of Ectomycorrhizal Symbionts of Regenerating Trees but Not Their Mycorrhizal Networks. Front. Plant Sci. 12:641232. doi: 10.3389/fpls.2021.641232
Received: 13 December 2020; Accepted: 06 April 2021;
Published: 14 May 2021.
Edited by:
Douglas Godbold, University of Natural Resources and Life Sciences Vienna, AustriaReviewed by:
Miguel Montoro Girona, Université du Québec en Abitibi Témiscamingue, CanadaRodica Pena, University of Reading, United Kingdom
Copyright © 2021 Veselá, Vašutová, Edwards-Jonášová, Holub, Fleischer and Cudlín. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Petra Veselá, vesela.p@czechglobe.cz