Luis Daniel Prada-Salcedo1,2,3*
Juan Pablo Prada-Salcedo4
Anna Heintz-Buschart1,3,5
François Buscot1,3
Kezia Goldmann1- 1Department Soil Ecology, Helmholtz-Centre for Environmental Research (UFZ), Halle, Germany
- 2Department of Biology, University of Leipzig, Leipzig, Germany
- 3German Centre for Integrative Biodiversity Research (IDiv), Leipzig, Germany
- 4Department of Bioinformatics, Biocenter, University of Würzburg, Würzburg, Germany
- 5Swammerdam Institute for Life Sciences, University of Amsterdam, Amsterdam, Netherlands
Depending on their tree species composition, forests recruit different soil microbial communities. Likewise, the vertical nutrient gradient along soil profiles impacts these communities and their activities. In forest soils, bacteria and fungi commonly compete, coexist, and interact, which is challenging for understanding the complex mechanisms behind microbial structuring. Using amplicon sequencing, we analyzed bacterial and fungal diversity in relation to forest composition and soil depth. Moreover, employing random forest models, we identified microbial indicator taxa of forest plots composed of either deciduous or evergreen trees, or their mixtures, as well as of three soil depths. We expected that forest composition and soil depth affect bacterial and fungal diversity and community structure differently. Indeed, relative abundances of microbial communities changed more across soil depths than in relation to forest composition. The microbial Shannon diversity was particularly affected by soil depth and by the proportion of evergreen trees. Our results also reflected that bacterial communities are primarily shaped by soil depth, while fungi were influenced by forest tree species composition. An increasing proportion of evergreen trees did not provoke differences in main bacterial metabolic functions, e.g., carbon fixation, degradation, or photosynthesis. However, significant responses related to specialized bacterial metabolisms were detected. Saprotrophic, arbuscular mycorrhizal, and plant pathogenic fungi were related to the proportion of evergreen trees, particularly in topsoil. Prominent microbial indicator taxa in the deciduous forests were characterized to be r-strategists, whereas K-strategists dominated evergreen plots. Considering simultaneously forest composition and soil depth to unravel differences in microbial communities, metabolic pathways and functional guilds have the potential to enlighten mechanisms that maintain forest soil functionality and provide resistance against disturbances.
Introduction
Forest soil ecosystems hold high amounts of microbial biomass and diversity (Llado et al., 2018; He et al., 2020). Linked to the litter and exudate inputs from trees and understory vegetation, soil microbes sustain multiple ecosystem functions like decomposition, nutrient cycling, primary production, and multi-trophic interactions (Binkley and Giardina, 1998; Mori et al., 2017).
Soil microbial communities can be strongly affected by forest composition and are therefore known to react sensitively to forest conversion (Goldmann et al., 2015; Dukunde et al., 2019). Impacts of forest composition on soil microbes can be coupled to individual tree taxonomy and to differences at the functional level between evergreen and deciduous trees (Prescott and Grayston, 2013). Replacement of beech by spruce, for instance, is associated with changes in soil structure, including humus form and soil acidity, consequently impacting soil bacterial and fungal community composition (Berger and Berger, 2012; Nacke et al., 2016). Furthermore, divergent litter and root exudates of evergreen and deciduous trees can directly affect soil microbes (Eisenhauer et al., 2017). Indirectly, microbial communities are also shaped by interactions with other members of the soil food web (e.g., protists or nematodes), which in turn are also affected by forest tree composition (Bonkowski, 2004; Geisen et al., 2018). Litter fall often leads to an enriched microbial diversity and activity in the topsoil (Thoms et al., 2010; Uri et al., 2012; Ana et al., 2015). Consequently, and due to reduced vertical nutrient availability, it is assumed that bacterial and fungal diversity declines with increasing soil depth (Jobbágy and Jackson, 2001; Goebes et al., 2019). However, rhizodeposition affects microbes at deeper soil horizons (Certini et al., 2004; Lopez et al., 2020). Hence, microbial communities and their activities are differently affected by multiple factors in the topsoil compared to the subsoil (Blume et al., 2002; Loeppmann et al., 2016). However, frequently, only the top 10 cm of soil is studied, although microbial depth-dependence was recently reported: For instance, multiple bacterial and archaeal taxa in a poplar plantation exhibited significant changes across different soil depths, and indicated that main community transitions occurred at 10–20 cm (Feng et al., 2019). For soil fungi, a study in a Betula albosinensis forest detected that the abundance of the genus Paxillus decreased, whereas the abundances of Inocybe and Genabea increased with increasing soil depth (Du et al., 2017). Thus, ecological studies should consider at least the upper 30 cm (Richter and Markewitz, 1995; Goebes et al., 2019; Yost and Hartemink, 2020). While forest composition and soil depth are often considered independently, their interplay affects soil microbial communities, as for example, deciduous, evergreen, and mixed forests are related to different humus forms (Swift et al., 1979; Ponge, 2003). Likewise, soil profiles differ between deciduous and evergreen forests, where the former is characterized by strong mixing of organic material with the mineral soil components, and the latter by a distinct boundary between the mineral and organic horizons (Adams et al., 2019).
Although soil bacteria and fungi are affected by similar factors, and are closely interconnected, they differ in community structure, size, species richness, life history, and enzymatic capacity (Deveau et al., 2018). Therefore, bacterial and fungal assembly in relation to forest composition and soil depth varies from each other (Sun et al., 2017). In general, bacterial or fungal dominance depends on substrate degradation potential and nutrient availability (Van Der Heijden et al., 2008; Rinnan and Bååth, 2009). Fungi dominate soil systems with a high ratio of carbon to nitrogen (C:N), e.g., forest or no-till agricultural soils, while bacteria prefer systems with low C:N ratios, e.g., grasslands or conventional till agricultural soils (Moore et al., 2002). Consequently, different habitats shape specific microbial communities (Uroz et al., 2012; Bardgett et al., 2014). The divergence in environmental and nutritional requirements leads to niche differentiation and selection of organisms with different lifestyle strategies (Deveau et al., 2018).
Recent studies focusing on microbial functions, i.e., considering different bacterial metabolic pathways and fungal guilds, considered either the impact of forest composition or different soil depths. For instance, the use of predicted functional gene profiles based on 16S rRNA gene sequences facilitated the understanding of the vertical distribution of functional microbial groups (Uroz et al., 2013; Mushinski et al., 2018; Zhang et al., 2019) in relation to nutrient use across soil profiles in different ecosystems (Tripathi et al., 2019; Luan et al., 2020). Concerning fungal guilds, mycorrhizal fungal composition differed significantly in spruce compared to beech-dominated plots, and spruce-dominated plots had a higher abundance of saprotrophic relative to mycorrhizal fungi (Asplund et al., 2019). Likewise, analyzing fungal guilds in different soil layers indicated that competition for resources between saprotrophic and ectomycorrhizal (EcM) fungi led to restrictions in saprotrophic fungal dominance only in the topsoil in forests dominated by EcM trees (Carteron et al., 2020).
Considering the mentioned constraints, our study investigated how bacterial and fungal soil communities vary in relation to forest composition and three soil depths, while deciphering the potential microbial functional roles and contribution to forest ecosystem processes. We aimed to identify soil microbial taxa associated with deciduous forests, evergreen forests, or their mixtures vertically along the soil profiles until 30 cm depth. We expected that forest composition and soil depth exert different effects in shaping soil bacterial and fungal diversity and community structure. Precisely that microbial diversity decreases with deeper soil depth and an increasing proportion of evergreen trees. Due to their differences in the ability to process organic substrates, connection to plant roots, and cellular structures, we further anticipated that bacterial communities vary more than fungal communities along the soil profiles. In addition, we hypothesize that fungal communities are stronger shaped by forest composition than soil bacteria. At the functional level, we expect that bacterial pathway and fungal guild abundances will reduce with the deeper soil depth and with a high proportion of evergreen trees.
Materials and Methods
Study Sites
The study sites were composed of different “forest types” located in different European countries, which included boreal forests (Finland), hemi-boreal forests (Poland), and mountainous beech forests (Romania), all studied in the frame of the project SoilForEurope (http://websie.cefe.cnrs.fr/soilforeurope/). The plot design and sampling procedure have been described previously (Prada-Salcedo et al., 2021b). Our study used samples retrieved from 44 mature forest plots of 30 × 30 m2 from the FunDivEUROPE platform (Baeten et al., 2013). In spring 2017, the soil was sampled at subplot level (10 × 10 m2), considering five repetitions per plot using a split-tube sampler (diameter 5.3 cm). In addition, three different “soil depths” (0–10, 10–20, and 20–30 cm) were sampled from intact soil cores. Our plots contained monospecific and multispecies forest systems, and the “forest composition” categories consisted of either deciduous forests on eleven plots that comprise the tree species Acer pseudoplatanus, Betula pendula, Carpinus betulus, Fagus sylvatica, and Quercus robur; or evergreen forests on eleven plots including the tree species Abies alba, Picea abies, and Pinus sylvestris; or 22 mixed forests plots with evergreen and deciduous composition comprising eleven tree combinations (Supplementary Table 1).
Molecular Procedures and Bioinformatics
Soil samples for molecular analyses were transported at 4°C and stored at −20°C until processing. The molecular procedure was described previously (Prada-Salcedo et al., 2021a,b). Briefly, total genomic DNA was extracted for each subplot and layer sample using the Power Soil™ DNA Isolation Kit (QIAGEN Laboratories Inc., Solana Beach, USA). After pooling DNA at the plot level, bacterial and fungal amplicon libraries were prepared by performing amplification of 16S rRNA gene V4 region (primers: P5_8N_515F and P5_7N_515F together with P7_2N_806R and P7_1N_806R) (Caporaso et al., 2011; Moll et al., 2018) and ITS2 (primers: P5-5 N-ITS4 and P5-6 N-ITS4 together with P7-3 N-fITS7 and P7-4 N-fITS7) (Gardes and Bruns, 1993; Ihrmark et al., 2012; Leonhardt et al., 2019), respectively. We performed a paired-end sequencing of 2 × 300 bp, using the MiSeq Reagent kit v3 on an Illumina MiSeq platform at the Department Soil Ecology of the Helmholtz-Center for Environmental Research—UFZ in Halle (Saale), Germany.
Raw sequences were extracted based on the presence of primers and then subsampled to 140,000 reads to compare samples in equal conditions. We used the standardized workflow of the dadasnake pipeline (Weißbecker et al., 2020; https://github.com/a-h-b/dadasnake), which is based on the DADA2 algorithm (Callahan et al., 2016). DADA2 produces amplicon sequence variants (ASVs), which bypass the previous common clustering at 97% sequence similarity. The ASV approach is more robust and realistic, and leads to a better resolution of microbial communities, especially in very rich environments such as soil (Joos et al., 2020). We used mostly default conditions of dadasnake pipeline, adjusting only the following parameters: minimum read length was set to 170 bp for bacterial read sequences and to 70 bp for fungal forward and reverse read sequences; maximum expected error after truncation was 0.5 for bacterial and 3 for each fungal sequence, respectively. Forward and reverse sequence variants were merged, and the chimeric sequences were removed by the consensus method. The bacterial ASVs were classified against the bacterial sequences in the SILVA database (Quast et al., 2013) (release 138). The prediction of functional profiles based on the bacterial ASVs was done using the Tax4Fun2 default database and the “group pathways” output (Wemheuer et al., 2020). Tax4Fun2 predicts gene abundances from 16S rRNA gene sequences and returns gene groups from the KEGG database that correspond to pathways at different levels, e.g., pathway IDs corresponding to Level 1: ko00190, ko00195, ko00680, ko00710, ko00720, ko00910, and ko00920; this ID pathways are corresponding to Level 2: oxidative phosphorylation, photosynthesis, methane metabolism, carbon fixation, photosynthetic organisms, nitrogen metabolism, and sulfur metabolism; and these pathways are grouped to the corresponding level 3: energy metabolism (see Supplementary Table 2 for all retrieved and used pathways). Fungal ASVs were checked and trimmed using ITSx (Bengtsson-Palme et al., 2013), and taxonomy was assigned using the UNITE (2017—Version 01.12.2017) database (Nilsson et al., 2019). Finally, we classified fungal guilds according to the assigned ASV taxonomy using FunGuild, using the main guilds of pathogens/parasites, saprotrophic, symbiotic (i.e., endophyte and epiphyte), mycorrhizal (i.e., ectomycorrhizal, arbuscular mycorrhizal, and orchid mycorrhizal), multi-lifestyle (taxa that shift guilds during their lifecycle), and undefined taxa (Nguyen N. H. et al., 2016).
Statistical Methods
Abundances, Microbial Diversity, Microbial Community Composition, and Functional Abundances
From a total of 16,358 bacterial ASVs and 4,635 fungal ASVs, we evaluated the fungal and bacterial communities by using the R packages “phyloseq” (McMurdie and Holmes, 2013) and “vegan” (Oksanen et al., 2017) in the R software (R 2021-R version 4.1.2). For evaluation of abundances and microbial diversity, sequencing data was rarefied at the cutoff of 25,000 reads per sample for bacteria and 8,000 reads per sample for fungal communities, maintaining sample saturation and the balance between different soil depths (Supplementary Figure 1). Kruskal–Wallis tests followed by Dunn post-hoc tests were used to compare relative abundances at phyla level between forest composition and soil depth variables (Mendiburu, 2010). To evaluate relationships between microbial diversity, forest composition, soil depth, and possible interactions, we performed linear mixed models (LMMs) using forest type (country) as a random effect, since we consider this as the main source of variation, with the glmmTMB package (Brooks et al., 2017). Distances of microbial communities between different samples were calculated using Bray-Curtis dissimilarity, visualized by non-metric multidimensional scaling (NMDS), and statistically tested by permutational multivariate analysis of variance (PERMANOVA) using “vegan.” To see how forest composition and soil depth impacted microbial diversity at the functionality level, we used bacterial pathways and fungal guild abundances. The abundances of bacterial pathways were normalized by the “normalize” function with the “scale” method using the “BBmisc” package (Bischl et al., 2017). Then, the abundances of bacterial pathways and abundances of fungal guild were evaluated in relation to the proportion of evergreen trees in the plot and soil depth to fit generalized linear models (GLMs) with the glm function.
Random Forest Model Classifier for Forest Composition and Soil Depth Analysis
Random forest (RF) classifiers (Breiman, 2001) were trained on bacterial and fungal Hellinger-transformed abundance tables (Legendre and Gallagher, 2001), which could predict different classes (referred to as “labels”) of forest composition and soil depths from various microbial features (phyla, families, genera, and ASVs). In addition to the taxonomy levels, the Tax4Fun2-predicted functional profiles were used as features of specific bacterial pathways. For each feature-class pair, 100 random forests of 500 decision trees were trained using the default settings of the “randomForest” function implemented in the randomForest package (Liaw and Wiener, 2002), performed in the R and using 130 samples for bacteria and 131 samples for fungi. We kept the best random forest out of the 100 computed ones according to the out of bag (OOB) error. Furthermore, a group of 100 random forest classifiers was trained and tested by randomly created groups of 102 training samples and 29 testing samples for each feature-class pair. The test accuracy of each feature-class pair was averaged to validate the accuracy of the random forest models.
We obtained feature importance measures from the best RF according to the OOB error. These measurements were the mean decrease of the Gini index and the mean decrease in accuracy. The mean decrease of the Gini index gives an insight into how often a certain feature was successfully used to split the samples with respect to the used labels. Thereby, this importance measure is a proxy for an association which can be understood as an ecological indicator. Likewise, the mean decrease accuracy is obtained by permuting the value of the feature and evaluating the effect (decrease in accuracy) of the permutation (Breiman, 2001; Thompson et al., 2019). Neither importance measure necessarily indicates ecological importance, but signifies suitability as an indicator for a condition.
Results
Bacterial and Fungal Abundances, Diversity, and Community Composition in Relation to Forest Composition and Soil Depth
To gain a general overview of how forest composition and soil depth impact the community structure of soil bacterial and fungal communities, we first evaluated relative abundances at the phylum level. At this level, forest composition appeared not to be a driver of bacterial structure. However, slight differences were observed in Acidobacteriota, Actinobacteriota, Nitrospirota, Planctomycetota, Proteobacteria, RCP2-54, and Verrucomicrobiota. In contrast, soil depth resulted in certain patterns. Proteobacteria showed the greatest relative sequence abundances, which decreased with soil depth (Supplementary Table 3). Similar trends were found for Actinobacteriota and Bacteroidota. Conversely, Chloroflexi, Gemmatimonadota, and Methylomirabilota increased in abundance along with deeper soil layers. The relative abundances of Acidobacteriota, Verrucomicrobiota, Nitrospirota, and Planctomycetota slowly differed across soil depths, and Firmicutes did not show differences across depth (Figure 1A; Supplementary Table 3). For fungi, the greatest relative abundances of reads were affiliated to Ascomycota and Basidiomycota, followed by Mortierellomycota, Mucoromycota, and unclassified fungi. Other fungal phyla were scarcely detected. For purely deciduous forest plots, we detected an increasing relative sequence abundance of Basidiomycota with soil depth and decreasing relative sequence abundances of Ascomycota and Mortierellomycota (Supplementary Table 4). In mixed forests, the Mucoromycota increased in relative sequence abundance with soil depth. Furthermore, Ascomycota and Mortierellomycota showed higher relative sequence abundances in mixed forest plots compared with either deciduous or evergreen ones, particularly at 20–30 cm (Figure 1B; Supplementary Table 4).
Figure 1. Relative abundances of (A) bacterial and (B) fungal phyla in relation to forest composition, split between deciduous, mixture, and evergreen forests (vertical panels), and soil depths. The data represent values of taxa with a relative abundance higher than 0.1 (for statistical details see Supplementary Tables 3, 4).
In a second step, we analyzed microbial Shannon diversity and ASV richness in relation to forest composition and sampling depth. Our results revealed general patterns of Shannon diversity in bacteria and fungi in relation to forest composition. With more evergreen trees present, the microbial diversity decreased (bacterial P = 0.0041; ANOVA, P = 0.005624; PseudoR2 = 0.54; fungal P = 0.0292; ANOVA, P = 7.567e-06; PseudoR2 = 0.52; Figures 2A,B). Moreover, bacterial Shannon diversity and ASV richness differed along with soil depth, whereby the topsoil harbored the highest diversities (Figure 2A; Supplementary Table 5). In the upper 20 cm, the bacterial evenness was higher than at lower depths (Supplementary Table 5). Fungi also displayed the highest Shannon diversity, ASV richness, and evenness in the topsoil, but the values were similar in the 20–30 cm soil layer (Figure 2B; Supplementary Table 5).
Figure 2. Diversity and community composition associated with forest compositions and soil depth: (A) Bacterial and (B) Fungal Shannon diversity in relation to evergreen proportions, lines represents linear model estimates and shaded areas represent 95% confidence; (C) Bacterial and (D) Fungal community composition depicted as NMDS scaling based on Bray-Curtis dissimilarity, and different colors represent the strongest impacting variables according to PERMANOVA (for statistical details see Supplementary Tables 5, 6).
The forest type was the greatest driver of microbial community composition (Supplementary Table 6). Additionally, soil depth impacted bacterial communities more than forest composition (Figures 2C,D). Conversely, fungal communities were shaped by forest composition, which was verified by PERMANOVA performed individually for each forest type (Figures 2C,D; Supplementary Table 6).
Bacterial Pathways and Fungal Guilds Highly Associated With Forest Composition and Soil Depth
We used GLM models to estimate microbial roles at the functional level using bacterial pathways and fungal guilds. In general, of the 96 bacterial pathways, 22 responded to forest composition and 77 to soil depth (Supplementary Tables 7, 8). The models with bacterial pathway abundances using evergreen tree proportion as an explanatory variable did not show a positive association with functional core pathways linked to energy metabolisms, such as carbon fixation, degradation, photosynthesis, or nitrogen. However, significant responses were detected in more specific bacterial pathways related to metabolisms of xenobiotics, lipids, amino acids, biosynthesis of secondary metabolites, and the pathways of carbohydrate degradation like starch, sucrose, and propanoate metabolism. For most of these positive response pathways, with the exception of quorum sensing, the results suggested an increasing bacterial pathway abundance with a higher proportion of evergreen trees (Figure 3A; Supplementary Table 7). Most of the core pathways, such as carbon fixation, nitrogen metabolism, or lipid metabolism, responded to soil depth, and most of the pathways were significantly more represented in the topsoil. In contrast, quorum sensing and atrazine degradation were positively related to deeper soil layers (Figure 3B; Supplementary Table 8).
Figure 3. GLM of bacterial pathways related to evergreen tree proportion (A) and soil depth (B). Estimates are the intercept of the model and the color represents the slope. Diamonds represent specific pathways that respond with significant differences according to ANOVA (p < 0.05) of the GLM. Pathways with significant effects in relation to evergreen tree proportion and soil depth were labeled with specific pathways (Level 2), and model details are provided in Supplementary Tables 7, 8.
Regarding the fungal guilds associated with evergreen tree proportion, the results showed that the saprotrophs, arbuscular mycorrhizal fungi, and plant pathogens had significant responses and showed higher abundances in deciduous forests and mixtures compared to evergreen forests (Figure 4A; Supplementary Table 9). Our results indicated a response of soil saprotrophs, orchid mycorrhizal fungi, and different pathogens/parasites particularly in the topsoil (Figure 4B; Supplementary Table 10).
Figure 4. GLM of fungal guilds related to evergreen tree proportion (A) and soil depth (B). Estimates are the intercept of the model and the color represents the slope. Diamonds represent specific guilds that respond significantly according to ANOVA (p < 0.05) of the GLM. Fungal guilds with significant effects in relation to evergreen tree proportion and soil depth were labeled with specific guilds (Guild 2), and model details are given in Supplementary Tables 9, 10.
Microbial Taxa Closely Associated With Forest Composition and Soil Depth
We used the RF approach to determine microbial taxa indicative of forest composition and soil depth. According to our expectations, the RF model detected microbial indicator features, i.e., microbial taxa with a potential relationship to one or both considered variables. Throughout, the estimated prediction accuracy of our RF models improved from broad to fine taxonomical levels, with ASVs being the best predicting feature for bacteria and fungi. At the bacterial ASV level, the RF validation with the training/test (the mean of the validation performance to classified 29 tested samples) models reached an estimated accuracy of 75% and 84% for forest composition and soil depth, respectively. The fungal ASV-based RF models delivered an estimated accuracy of 81% and 69% for forest composition and soil depth, respectively (Supplementary Figure 2). Moreover, the RF models were used to predict bacterial metabolic pathways and fungal functional guilds. The RF validation estimated accuracy at these functional levels was lower than the ASV taxonomy level. The bacterial pathways displayed an estimated accuracy of 57% and 62% for forest compositions and soil depth, respectively. Concerning fungal guilds, RF reached an estimated accuracy of 58% and 51% for forest composition and soil depth, respectively (Supplementary Figure 2).
We identified the top 10 bacterial and fungal genera associated with forest composition and soil depths by considering the individual mean importance in the best RF model (Tables 1, 2). Most of the microbial genera indicative of soil depth differed from those related to forest composition. For example, the bacterial genus Devosia was highly linked to soil depth, but not to forest composition. Contrastingly, the genus GOUTA6 (Nitrosomonadacea) was highly indicative of forest composition, but less indicative of soil depth. However, the bacterial genera Puia and Caulobacter were comparably associated with both forest composition and soil depth (Tables 1, 2).
Table 1. Top 10 bacterial and fungal genera associated with forest composition based on our RF model.
Table 2. Top 10 bacterial and fungal genera associated with soil depth based on our RF model.
In relation to forest composition, deciduous, mixed, and evergreen forests diverged in their associated bacterial communities. For example, the genus GOUTA6 was an important feature in the RF indicating deciduous and evergreen forests of bacterial communities. Sediminibacterium was highly associated with deciduous and mixed forests, while Inquilinus was linked to evergreen forests (Table 1). Considering the RF importance for the genera associated with soil depth, for some bacteria, importance decreased with depth, e.g., Devosia or Granulicella, whereas other bacterial genera were more indicative of a particular depth, such as Chthoniobacter and Acidipila for 0–10 cm, Cytophaga and Conexibacter for 10–20 cm, and 1921–3 (Ktedonobacteraceae) and Rhodoplanes for 20–30 cm (Table 2).
Fungal genera linked to forest composition and soil depths were Meliniomyces, Cladophialophora, and Hygrophorus, besides unclassified taxa (Tables 1, 2). Fungi associated with forest composition were Trichophaea with the highest RF importance in deciduous and mixed forests, and Sistotrema, Inocybe, and Solicoccozyma with importance in mixed and evergreen forests (Table 2). The fungal genera linked to soil depth appeared to be rather related to the uppermost or lowest depths considered, e.g., Cortinarius and, in particular, Tomentella showed a high RF importance in 0–10 and 20–30 cm. In contrast, Mycena was highly associated with the first 20 cm, while Ilyonectria appeared to be more linked to 20–30 cm (Table 2).
Discussion
General Bacterial and Fungal Niche Preferences
Our study aimed to investigate bacterial and fungal diversity in relation to forest composition, in terms of the proportion of both evergreen and deciduous trees, and soil depth. In accordance with previous studies, our results showed that an increasing proportion of evergreen trees negatively affected microbial Shannon diversity, which also decreased along with soil depth (Jobbágy and Jackson, 2001; Goebes et al., 2019). In line with Osburn et al. (2019), bacterial communities were shaped by soil depth physicochemical variables, while fungal assemblies were rather shaped by forest composition. The strong link between soil fungi, above- and belowground litter, and also root exudates might be responsible for these patterns (Eisenhauer et al., 2017). The vertical gradient from copiotrophic conditions in the topsoil, i.e., rich in organic compounds, to oligotrophic conditions, i.e., simple minerals, with increasing depth conditions (Jobbágy and Jackson, 2001; Goebes et al., 2019) was apparently a more dominant driver than forest composition for bacterial communities. The bacterial presence in oligotrophic environments is associated with their capability to solubilize key nutrients, such as phosphorus and potassium, from primary minerals (Uroz et al., 2009, 2012; Prada-Salcedo et al., 2014). In contrast, the higher association of fungi with forest composition is perhaps due to their greater prevalence to interact with plants and their wider enzymatic potential to decompose complex organic substrates and recalcitrant compounds from plant and fungal necromass (Fabian et al., 2017; Algora Gallardo et al., 2021). Though we found no interactions between the proportion of evergreen and soil depth affecting microbial diversity, there are reports of such effects. For instance, the different litter types and decay rates between broadleaf and coniferous (needleleaf) trees differentially modify pH and base cation content of soils, and consequently, particular taxa are indirectly selected at surface soil layers (Augusto et al., 2015). Moreover, at deeper soil horizons, the rooting growth system of deep vs. shallow roots, and different rates between the lifespan of short and long roots of broadleaf and coniferous trees affect humus form and soil structure, causing another indirect effect on microbial diversity (Berger and Berger, 2012). These effects can be more pronounced in rhizospheres because the amount and nature of root exudates also change depending on the tree species (Prescott and Grayston, 2013). Therefore, complementary studies on particular microbial groups, e.g., within the rhizosphere or symbionts, could extend the knowledge on how tree–soil interactions affect the microbial communities.
Potential Bacterial Predicted Pathways and Fungal Guilds
The pathway predictions showed the impact of the proportion of evergreen trees on bacterial pathways, suggesting that the forest composition and mixed forest might affect the microbial processes. By using the same soil samples as this study, Gillespie et al. (2020) experimentally tested microbial responses to drought. Thereby, they found a less decrease in microbial respiration and denitrification in soils obtained from mixed compared to monospecific forests. Likewise, soil microbial stress levels were lower in mixed than in monospecific forests. This enhancement of microbial functional resistance in forest mixtures could be attributed to a higher number of microbial interactions, which primarily carry positive plant–microbe feedbacks (Prada-Salcedo et al., 2021a). Our results revealed that forest composition has a greater relationship with specific pathways, like xenobiotic metabolism, suggesting that more diverse forests could take up xenobiotics at enhanced rates and through different paths (Salem et al., 2017; Desai et al., 2019). Similarly, the higher response of secondary metabolism biosynthesis in top soils with higher evergreen proportions could induce more heterogeneous microenvironments with greater opportunities for interactions and competition that trigger the production of secondary metabolites (Sharrar et al., 2020).
Responses of the fungal guilds, particularly with regard to saprotrophic and root-associated fungi, have been noticed before (Bödeker et al., 2016; Peršoh et al., 2018). Additionally, fungal taxa with multiple lifestyles differ across depths. This group supposedly combines different lifestyles (i.e., biotrophic, necrotrophic, and saprotrophic) and often cope with limited nutritional dependencies, thus several taxa with undefined functions can adapt depending on the niche conditions (Lewis, 1973; Suzuki and Sasaki, 2019). Our results suggest higher responses and abundances (Supplementary Figure 3) of potentially pathogenic fungi in deciduous and mixtures forest, due to the higher number of potential niches and hosts in more diverse forests. Our results are consistent with those of (Nguyen D. et al., 2016), which used the same platform and plots to evaluate the incidence of foliar fungal disease in relation to tree species diversity. The authors found a decrease in the disease incidence with tree species richness in conifers, but not in broadleaved trees for which the incidence tended to increase with tree diversity. These results could suggest that deciduous and mixture forests favor conditions to maintain higher overall pathogen loads. But this does not necessarily mean higher disease levels, as forests dominated by one tree species have shown to be more susceptible to pests and diseases because the huge nutrient source and homogeneous habitat of monocultures normally rise to rapid infection rates (Liu et al., 2018). Fungal saprotrophs were predominantly linked to deciduous forests. In accordance, Chen et al. (2019) evaluated differences between fungi associated with subtropical evergreen and deciduous forests, finding higher abundances of saprophytic fungi in deciduous forests and deducing that these forests mainly harbor fast-growing copiotrophic fungi. Our results also show that within the subsoil, symbiotrophs, mainly mycorrhizal fungi were found (Supplementary Figure 3). Fungi of this guild receive photo-assimilated carbon from plants in exchange for the supply of poorly soluble nutrients that mycorrhizal fungi mobilize from rock weathering. Therefore, it is evident that this guild also responds in the greatest depth (Thorley et al., 2015). The structural and functional diversity of mycorrhizal fungi concerning short- to long-distance extension of their extra-radical mycelium and rhizomorphs enable them an optimal soil exploration (Agerer, 2006) and their dominance in the subsoil.
Random Forest Modeling Related to Microbial Communities in Forest Ecosystem
The RF approach contributes to identifying indicator species for a particular class of samples or conditions. Thus, we could show the need to study microbial abundances with different methods that may reveal the role of less abundant or unknown microbes, which also participate in the biogeochemical cycles (Jousset et al., 2017). The results of the RF using fine taxonomic ranks, such as genera or ASV, increased classifier accuracy. Functional predictions either from bacterial pathways or fungal guilds did not offer such high levels of accuracy compared to the taxonomical output. The lack of high accuracy values by using functional features, such as RF features, could be explained by the high functional redundancy in soils (Louca et al., 2018). Hence, the model accuracy is probably affected, especially, at broad functional levels like bacterial pathways and fungal guilds because these can be potentially performed by a wide range of extant taxa that cover “core” functions (e.g., respiration, nitrogen, and phosphorus cycling) (Jia and Whalen, 2020).
Indicator Taxa and Community Structure According to Forest Composition and Soil Depth
Based on the literature, we characterized some of the top indicator taxa from this study to evaluate whether specific conditions of the forest composition or soil depth favor particular microbes and lifestyle strategies. We found bacterial genera, such as Sediminibacterium, Rhodomicrobium, and Jatrophihabitans, that were indicator taxa within deciduous forests and are reported to display microbial traits common to r-strategists, i.e., higher performances in resource-rich environments (Duchow and Douglas, 1949; Qu and Yuan, 2008; Kim et al., 2015; Papp et al., 2020). These taxa are involved in the fast turnover of litter and consequently fast degradation of cellulose and hemicellulose compounds (Vesterdal et al., 2008; Adams et al., 2019). In contrast, linked to evergreen plots, we found bacterial indicator genera like Inquilinus, Acidocella, or Alicyclobacillus that could be adapted to typically acidic soils (Jung et al., 2011; Okamoto et al., 2017), and Longimycelium and Actinomadura (Zakharova et al., 2003; Xia et al., 2013), characterized by slow growth and efficient use of recalcitrant carbon resources and secondary metabolite production, which are traits rather associated with K-strategists and environments where resource supply rates are lower (Ho et al., 2017). Concerning fungi, interestingly, our RF revealed that free-living unicellular fungi were indicators of evergreen forest plots. Yeasts like Solicoccozyma, Apiotrichum, and unclassified Saccharomycetales are commonly isolated from forest soils (Liu et al., 2015; Yurkov et al., 2015; James et al., 2016; Mašínová et al., 2016). However, our results showed their clear niche preference for the evergreen forests. According to Birkhofer et al. (2012) and Yarwood et al. (2010), these yeasts can tolerate low pH values, and sandy and drained soils, which are common in evergreen forests (Retallack, 2005; Adams et al., 2019).
The link of the bacterial phyla Bacteroidota and Proteobacteria with the first soil centimeters corresponds to previous studies (Feng et al., 2019). This link is likely due to the copiotrophic properties and associations to aerobic niches in genera like Devosia, Ferruginibacter, and Granulicella (Nakagawa et al., 1996; Lim et al., 2009; Pankratov and Dedysh, 2010). In general, there were fewer bacterial and fungal taxa of importance for 10–20 cm compared to that in the 0–10 and 20–30 cm, which is probably related to the changing conditions within this intermediate soil compartment, with more microaerophilic and mixing of copiotrophic/oligotrophic conditions (Fierer et al., 2003; Moll et al., 2015; Neira Román et al., 2015) acting as filters for soil microbes, and hence it is complicated to find a clear microbial-niche association or clear indicator taxa. Accordingly, the intermediate depth (10–20 cm) could be seen as a transition zone for soil microbial communities (Eilers et al., 2012; Feng et al., 2019).
Conclusion
Our analyses provide a comprehensive understanding of bacteria and fungi related to forest composition in terms of evergreen and deciduous trees, and to soil depth. Microbial diversity and functionality of bacterial pathways and fungal guilds were reduced across the vertical gradient. We demonstrate that bacterial community composition varied primarily with soil depth, whereas fungal communities were strongly influenced by tree composition. The microbes within the intermediate soil depth were less predictable compared to topsoil and subsoil. A higher proportion of evergreen trees decrease taxonomic microbial diversity, while increasing the prediction of bacterial pathways abundances. The fungal guilds, especially saprotrophs, arbuscular mycorrhizal fungi, and plant pathogens, were highly represented in deciduous and mixed forests. Likewise, we emphasized the potential impact of forest composition by identifying several microbial functional pathways and guilds, which can be critical in providing higher biosynthetic capacities that support resistance against disturbances, and maintaining forest soil functionality.
Finally, the RF classifier performed with higher accuracy at finer taxonomic resolution compared to predicted bacterial metabolic pathways, and fungal guilds revealed microbial taxa linked to forest composition and soil depth. The identification of microbial indicator taxa by random forest approach can be linked to their functionality, and consequently contributes to explaining plant–microbial interactions under different forest compositions and provides support for forest management decisions.
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 below: https://www.ebi.ac.uk/ena, PRJEB33611 and experiment accession ERX5307385–ERX5307648.
Author Contributions
FB conceptualized the study and acquired funding. LP-S together with SoilforEurope members did the sampling work, performed molecular data set analysis, analyzed the data, and wrote the manuscript. JP-S contributed to the machine learning method. KG and AH-B provided substantial contributions to ecological significance. FB, KG, and AH-B provided revisions for the data analyses. All authors contributed to the completion of the manuscript.
Funding
This research was part of the SoilForEUROPE project funded through the 2015–2016 BiodivERsA COFUND call for research proposals, with the national funders French National Research Agency (ANR, France), Belgian Science Policy Office (Belpo, Belgium), German Research Foundation (DFG, Germany), Research Foundation Flanders (FWO, Belgium), and the Swedish Research Council (Formas, Sweden). LP-S and FB are grateful to the DFG Grant BU 941 28-1 that supported this project.
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.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank the site managers Leena Finér from the Finnish Forest Research Institute (Metla), Bogdan Jaroszewicz from the University of Warsaw (Poland), Olivier Bouriaud from the Forest Research and Management Institute (ICAS, Romania), and Filippo Bussotti from the University of Florence (Italy) and all associates, as well as the SoilForEUROPE consortium for their assistance while soil sampling and design contributions. We kindly thank Beatrix Schnabel, Luise Kaiser, and Melanie Günther for support in the molecular laboratory and Illumina sequencing. The community composition data have been computed at the high-performance computing cluster EVE, a joint effort of both the Helmholtz-Center for Environmental Research-UFZ and the German Center for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig. We kindly thank the server admins.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.920618/full#supplementary-material
References
Adams, M. B., Kelly, C., Kabrick, J., and Schuler, J. (2019). “Chapter 6—temperate forests and soils,” in Developments in Soil Science. Busse, M., Giardina, C.P., Morris, D.M., and Page-Dumroese, D.S. eds (London: Elsevier), pp. 83–108. doi: 10.1016/B978-0-444-63998-1.00006-9
Agerer, R. (2006). Fungal relationships and structural identity of their ectomycorrhizae. Mycologic. Progr. 5, 67–107. doi: 10.1007/s11557-006-0505-x
Algora Gallardo, C., Baldrian, P., and López-Mondéjar, R. (2021). Litter-inhabiting fungi show high level of specialization towards biopolymers composing plant and fungal biomass. Biol. Fertil. Soils 57, 77–88. doi: 10.1007/s00374-020-01507-3
Ana, R., Blanca, S. P., G, R., S, Aurore, C., Benoit, M., et al. (2015). Compartmentalized and contrasted response of ectomycorrhizal and soil fungal communities of Scots pine forests along elevation gradients in France and Spain. Environ. Microbiol. 17, 3009–3024. doi: 10.1111/1462-2920.12894
Asplund, J., Kauserud, H., Ohlson, M., and Nybakken, L. (2019). Spruce and beech as local determinants of forest fungal community structure in litter, humus and mineral soil. FEMS Microbiol. Ecol. 95, 232. doi: 10.1093/femsec/fiy232
Augusto, L., De Schrijver, A., Vesterdal, L., Smolander, A., Prescott, C., and Ranger, J. (2015). Influences of evergreen gymnosperm and deciduous angiosperm tree species on the functioning of temperate and boreal forests. Biol. Rev. 90, 444–466. doi: 10.1111/brv.12119
Baeten, L., Verheyen, K., Wirth, C., Bruelheide, H., Bussotti, F., Finér, L., et al. (2013). A novel comparative research platform designed to determine the functional significance of tree species diversity in European forests. Perspect. Plant Ecol. Evol. Systemat. 15, 281–291. doi: 10.1016/j.ppees.2013.07.002
Bardgett, R. D., Mommer, L., and De Vries, F. T. (2014). Going underground: root traits as drivers of ecosystem processes. Trends Ecol. Evol. 29, 692–699. doi: 10.1016/j.tree.2014.10.006
Bengtsson-Palme, J., Ryberg, M., Hartmann, M., Branco, S., Wang, Z., Godhe, A., et al. (2013). Improved software detection and extraction of ITS1 and ITS2 from ribosomal ITS sequences of fungi and other eukaryotes for analysis of environmental sequencing data. Methods Ecol. Evol. 4, 914–919. doi: 10.1111/2041-210X.12073
Berger, T. W., and Berger, P. (2012). Greater accumulation of litter in spruce (Picea abies). compared to beech (Fagus sylvatica). stands is not a consequence of the inherent recalcitrance of needles. Plant Soil 358, 349–369. doi: 10.1007/s11104-012-1165-z
Binkley, D., and Giardina, C. (1998). Why do tree species affect soils? the warp and woof of tree-soil interactions. Biogeochemistry 42, 89–106. doi: 10.1007/978-94-017-2691-7_5
Birkhofer, K., Schöning, I., Alt, F., Herold, N., Klarner, B., Maraun, M., et al. (2012). General relationships between abiotic soil properties and soil biota across spatial scales and different land-use types. PLoS One 7, e43292. doi: 10.1371/journal.pone.0043292
Bischl, B., Lang, M., Bossek, J., Horn, D., Richter, J., and Surmann, D. (2017). BBmisc: Miscellaneous Helper Functions for B. Bischl. R package version 1.11. Available online at: https://CRAN.R-project.org/package=BBmisc (accessed June 27, 2022).
Blume, E., Bischoff, M., Reichert, J. M., Moorman, T., Konopka, A., and Turco, R. F. (2002). Surface and subsurface microbial biomass, community structure and metabolic activity as a function of soil depth and season. Appl. Soil Ecol. 20, 171–181. doi: 10.1016/S0929-1393(02)00025-2
Bödeker, I. T. M., Lindahl, B. D., Olson, Å., and Clemmensen, K. E. (2016). Mycorrhizal and saprotrophic fungal guilds compete for the same organic substrates but affect decomposition differently. Funct. Ecol. 30, 1967–1978. doi: 10.1111/1365-2435.12677
Bonkowski, M. (2004). Protozoa and plant growth: the microbial loop in soil revisited. New Phytologist 162, 617–631. doi: 10.1111/j.1469-8137.2004.01066.x
Brooks, M. E., Kristensen, K., van Benthem, K. J., Magnusson, A., Berg, C. W., Nielsen, A., et al. (2017). Modeling zero-inflated count data with glmmTMB. bioRxiv 17, 132753. doi: 10.1101/132753
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: High-resolution sample inference from Illumina amplicon data. Nature Methods 13, 581–583. doi: 10.1038/nmeth.3869
Caporaso, J. G., Lauber, C. L., Walters, W. A., Berg-Lyons, D., Lozupone, C. A., Turnbaugh, P. J., et al. (2011). Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proceed. Nat. Acad. Sci. 108, 4516–4522. doi: 10.1073/pnas.1000080107
Carteron, A., Beigas, M., Joly, S., Turner, B. L., and Laliberté, E. (2020). Temperate forests dominated by arbuscular or ectomycorrhizal fungi are characterized by strong shifts from saprotrophic to mycorrhizal fungi with increasing soil depth. Microb. Ecol. 82, 377–390. doi: 10.1007/s00248-020-01540-7
Certini, G., Campbell, C. D., and Edwards, A. C. (2004). Rock fragments in soil support a different microbial community from the fine earth. Soil Biol. Biochemistr. 36, 1119–1128. doi: 10.1016/j.soilbio.2004.02.022
Chen, L., Xiang, W., Wu, H., Ouyang, S., Lei, P., Hu, Y., et al. (2019). Contrasting patterns and drivers of soil fungal communities in subtropical deciduous and evergreen broadleaved forests. Appl. Microbiol. Biotechnol. 103, 5421–5433. doi: 10.1007/s00253-019-09867-z
Desai, M., Haigh, M., and Walkington, H. (2019). Phytoremediation: metal decontamination of soils after the sequential forestation of former opencast coal land. Sci. Total Environ. 656, 670–680. doi: 10.1016/j.scitotenv.2018.11.327
Deveau, A., Bonito, G., Uehling, J., Paoletti, M., Becker, M., Bindschedler, S., et al. (2018). Bacterial–fungal interactions: ecology, mechanisms and challenges. FEMS Microbiol. Rev. 42, 335–352. doi: 10.1093/femsre/fuy008
Du, C., Geng, Z., Wang, Q., Zhang, T., He, W., Hou, L., et al. (2017). Variations in bacterial and fungal communities through soil depth profiles in a Betula albosinensis forest. J. Microbiol. 55, 684–693. doi: 10.1007/s12275-017-6466-8
Duchow, E., and Douglas, H. C. (1949). Rhodomicrobium vannielii, a new photoheterotrophic bacterium. J. Bacteriol. 58, 409–416. doi: 10.1128/jb.58.4.409-416.1949
Dukunde, A., Schneider, D., Schmidt, M., Veldkamp, E., and Daniel, R. (2019). Tree species shape soil bacterial community structure and function in temperate deciduous forests. Front. Microbiol. 10, 19. doi: 10.3389/fmicb.2019.01519
Eilers, K. G., Debenport, S., Anderson, S., and Fierer, N. (2012). Digging deeper to find unique microbial communities: the strong effect of depth on the structure of bacterial and archaeal communities in soil. Soil Biol. Biochemistr. 50, 58–65. doi: 10.1016/j.soilbio.2012.03.011
Eisenhauer, N., Lanoue, A., Strecker, T., Scheu, S., Steinauer, K., Thakur, M. P., et al. (2017). Root biomass and exudates link plant diversity with soil bacterial and fungal biomass. Scientific Rep.s 7, 44641. doi: 10.1038/srep44641
Fabian, J., Zlatanovic, S., Mutz, M., and Premke, K. (2017). Fungal-bacterial dynamics and their contribution to terrigenous carbon turnover in relation to organic matter quality. ISME J 11, 415–425. doi: 10.1038/ismej.2016.131
Feng, H., Guo, J., Wang, W., Song, X., and Yu, S. (2019). Soil depth determines the composition and diversity of bacterial and archaeal communities in a poplar plantation. Forests 10, 550. doi: 10.3390/f10070550
Fierer, N., Allen, A. S., Schimel, J. P., and Holden, P. A. (2003). Controls on microbial CO2 production: a comparison of surface and subsurface soil horizons. Global Change Biol. 9, 1322–1332. doi: 10.1046/j.1365-2486.2003.00663.x
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
Geisen, S., Mitchell, E. A. D., Adl, S., Bonkowski, M., Dunthorn, M., Ekelund, F., et al. (2018). Soil protists: a fertile frontier in soil biology research. FEMS Microbiol. Rev. 42, 293–323. doi: 10.1093/femsre/fuy006
Gillespie, L. M., Fromin, N., Milcu, A., Buatois, B., Pontoizeau, C., and Hättenschwiler, S. (2020). Higher tree diversity increases soil microbial resistance to drought. Commun. Biol. 3, 377. doi: 10.1038/s42003-020-1112-0
Goebes, P., Schmidt, K., Seitz, S., Both, S., Bruelheide, H., Erfmeier, A., et al. (2019). The strength of soil-plant interactions under forest is related to a Critical Soil Depth. Scientific Rep. 9, 8635. doi: 10.1038/s41598-019-45156-5
Goldmann, K., Schöning, I., Buscot, F., and Wubet, T. (2015). Forest management type influences diversity and community composition of soil fungi across temperate forest ecosystems. Front. Microbiol. 6, 1300. doi: 10.3389/fmicb.2015.01300
He, L., Mazza Rodrigues, J. L., Soudzilovskaia, N. A., Barceló, M., Olsson, P. A., Song, C., et al. (2020). Global biogeography of fungal and bacterial biomass carbon in topsoil. Soil Biol. Biochemistr. 151, 108024. doi: 10.1016/j.soilbio.2020.108024
Ho, A., Di Lonardo, D. P., and Bodelier, P. L. (2017). Revisiting life strategy concepts in environmental microbial ecology. FEMS Microbiol. Ecol. 93, 6. doi: 10.1093/femsec/fix006
Ihrmark, K., Bödeker, I. T. M., Cruz-Martinez, K., Friberg, H., Kubartova, A., Schenck, J., et al. (2012). New primers to amplify the fungal ITS2 region—evaluation by 454-sequencing of artificial and natural communities. FEMS Microbiol. Ecol. 82, 666–677. doi: 10.1111/j.1574-6941.2012.01437.x
James, S. A., Bond, C. J., Stanley, R., Ravella, S. R., Péter, G., Dlauchy, D., et al. (2016). Apiotrichum terrigenum sp. nov., a soil-associated yeast found in both the UK and mainland Europe. Int. J. Systemat. Evol. Microbiol. 66, 5046–5050. doi: 10.1099/ijsem.0.001467
Jia, Y., and Whalen, J. K. (2020). A new perspective on functional redundancy and phylogenetic niche conservatism in soil microbial communities. Pedosphere 30, 18–24. doi: 10.1016/S1002-0160(19)60826-X
Jobbágy, E., and Jackson, R. (2001). The distribution of soil nutrients with depth: global patterns and the imprint of plants. Biogeochemistry 53, 51–77. doi: 10.1023/A:1010760720215
Joos, L., Beirinckx, S., Haegeman, A., et al. (2020). Daring to be differential: metabarcoding analysis of soil and plant-related microbial communities using amplicon sequence variants and operational taxonomical units. BMC Genom. 21, 733. doi: 10.1186/s12864-020-07126-4
Jousset, A., Bienhold, C., Chatzinotas, A., Gallien, L., Gobet, A., Kurm, V., et al. (2017). Where less may be more: how the rare biosphere pulls ecosystems strings. The ISME J. 11, 853–862. doi: 10.1038/ismej.2016.174
Jung, H.-M., Lee, J.-S., Bae, H.-M., Yi, T.-H., Kim, S.-Y., Lee, S.-T., et al. (2011). Inquilinus ginsengisoli sp. nov., isolated from soil of a ginseng field. Int. J. Systemat. Evol. Microbiol. 61, 201–204. doi: 10.1099/ijs.0.018689-0
Kim, S.-J., Moon, J.-Y., Lim, J.-M., Hamada, M., Ahn, J.-H., Weon, H.-Y., et al. (2015). Jatrophihabitans soli sp. nov., isolated from soil. Int. J. Systemat. Evol. Microbiol. 65, 1759–1763. doi: 10.1099/ijs.0.000173
Legendre, P., and Gallagher, E. D. (2001). Ecologically meaningful transformations for ordination of species data. Oecologia 129,: 271–280. doi: 10.1007/s004420100716
Leonhardt, S., Hoppe, B., Stengel, E., Noll, L., Moll, J., Bässler, C., et al. (2019). Molecular fungal community and its decomposition activity in sapwood and heartwood of 13 temperate European tree species. PLoS ONE 14, e0212120. doi: 10.1371/journal.pone.0212120
Lewis, D. H. (1973). Concepts in fungal nutrition and the origin of biotrophy. Biologic. Rev. 48, 261–277. doi: 10.1111/j.1469-185X.1973.tb00982.x
Lim, J. H., Baek, S. H., and Lee, S. T. (2009). Ferruginibacter alkalilentus gen. nov., sp. nov. and Ferruginibacter lapsinanis sp. nov., novel members of the family 'Chitinophagaceae' in the phylum Bacteroidetes, isolated from freshwater sediment. Int. J. Syst. Evol. Microbiol. 59, 2394–2399. doi: 10.1099/ijs.0.009480-0
Liu, C. L. C., Kuchma, O., and Krutovsky, K. V. (2018). Mixed-species vs. monocultures in plantation forestry: Development, benefits, ecosystem services and perspectives for the future. Glob. Ecol. Conserv. 15, e00419. doi: 10.1016/j.gecco.2018.e00419
Liu, X. Z., Wang, Q. M., Göker, M., Groenewald, M., Kachalkin, A. V., Lumbsch, H. T., et al. (2015). Towards an integrated phylogenetic classification of the Tremellomycetes. Stud. Mycol. 81, 85–147. doi: 10.1016/j.simyco.2015.12.001
Llado, S., Lopez-Mondejar, R., and Baldrian, P. (2018). Drivers of microbial community structure in forest soils. Appl. Microbiol. Biotechnol. 102, 4331–4338. doi: 10.1007/s00253-018-8950-4
Loeppmann, S., Blagodatskaya, E., Pausch, J., and Kuzyakov, Y. (2016). Enzyme properties down the soil profile—a matter of substrate quality in rhizosphere and detritusphere. Soil Biol. Biochemistr. 103, 23. doi: 10.1016/j.soilbio.2016.08.023
Lopez, S., van der Ent, A., Sumail, S., Sugau, J. B., Buang, M. M., Amin, Z., et al. (2020). Bacterial community diversity in the rhizosphere of nickel hyperaccumulator plant species from Borneo Island (Malaysia). Environ. Microbiol. 22, 1649–1665. doi: 10.1111/1462-2920.14970
Louca, S., Polz, M. F., Mazel, F., Albright, M. B. N., Huber, J. A., O'Connor, M. I., et al. (2018). Function and functional redundancy in microbial systems. Nat. Ecol. Evol. 2, 936–943. doi: 10.1038/s41559-018-0519-1
Luan, L., Liang, C., Chen, L., Wang, H., Xu, Q., Jiang, Y., et al. (2020). Coupling bacterial community assembly to microbial metabolism across soil profiles. mSystems 5, e00298–e00220. doi: 10.1128/mSystems.00298-20
Mašínová, T., Bahnmann, B. D., Větrovský, T., Tomšovský, M., Merunková, K., and Baldrian, P. (2016). Drivers of yeast community composition in the litter and soil of a temperate forest. FEMS Microbiol. Ecol. 93, 223. doi: 10.1093/femsec/fiw223
McMurdie, P. J., and Holmes, S. (2013). phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLOS ONE 8, e61217. doi: 10.1371/journal.pone.0061217
Mendiburu, F. (2010). Agricolae: statistical procedures for agricultural research. R package version 1, 1–8.
Moll, J., Goldmann, K., Kramer, S., Hempel, S., Kandeler, E., Marhan, S., et al. (2015). Resource type and availability regulate fungal communities along arable soil profiles. Microb. Ecol. 70, 390–399. doi: 10.1007/s00248-015-0569-8
Moll, J., Kellner, H., Leonhardt, S., Stengel, E., Dahl, A., Bässler, C., et al. (2018). Bacteria inhabiting deadwood of 13 tree species are heterogeneously distributed between sapwood and heartwood. Environ. Microbiol. 20, 3744–3756. doi: 10.1111/1462-2920.14376
Moore, J., McCann, K., Setälä, H., and Ruiter, P. (2002). Top-down is bottom-up: Does predation in the rhizosphere regulate above ground production. Ecology 84, 2. doi: 10.1890/0012-9658(2003)084[0846:TIBDPI]2.0
Mori, A. S., Lertzman, K. P., and Gustafsson, L. (2017). Biodiversity and ecosystem services in forest ecosystems: a research agenda for applied forest ecology. J. Appl. Ecol. 54, 12–27. doi: 10.1111/1365-2664.12669
Mushinski, R. M., Zhou, Y., Gentry, T. J., and Boutton, T. W. (2018). Bacterial metataxonomic profile and putative functional behavior associated with C and N cycle processes remain altered for decades after forest harvest. Soil Biol. Biochemistr. 119, 184–193. doi: 10.1016/j.soilbio.2018.01.008
Nacke, H., Goldmann, K., Schöning, I., Pfeiffer, B., Kaiser, K., Castillo-Villamizar, G. A., et al. (2016). Fine spatial scale variation of soil microbial communities under european beech and norway spruce. Front. Microbiol. 7, 2067. doi: 10.3389/fmicb.2016.02067
Nakagawa, Y., Sakane, T., and Yokota, A. (1996). Transfer of “Pseudomonas riboflavina” (Foster 1944), a gram-negative, motile rod with long-chain 3-hydroxy fatty acids, to Devosia riboflavina gen. nov., sp. nov., nom. rev. Int. J. Syst. Bacteriol. 46, 16–22. doi: 10.1099/00207713-46-1-16
Neira Román, J., Ortiz, M., Morales-Salinas, L., and Acevedo, E. (2015). Oxygen diffusion in soils: Understanding the factors and processes needed for modeling. Chilean J. Agricult. Res. 75, 35–44. doi: 10.4067/S0718-58392015000300005
Nguyen, D., Castagneyrol, B., Bruelheide, H., Bussotti, F., Guyot, V., Jactel, H., et al. (2016). Fungal disease incidence along tree diversity gradients depends on latitude in European forests. Ecol. Evol. 6, 2426–2438. doi: 10.1002/ece3.2056
Nguyen, N. H., Song, Z., Bates, S. T., Branco, S., Tedersoo, L., Menke, J., et al. (2016). FUNGuild: An open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 20, 241–248. doi: 10.1016/j.funeco.2015.06.006
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
Okamoto, R., Kojima, H., and Fukui, M. (2017). Acidocella aquatica sp. nov., a novel acidophilic heterotrophic bacterium isolated from a freshwater lake. Int. J. Systemat. Evol> Microbiol. 67, 4773–4776. doi: 10.1099/ijsem.0.002376
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2017). vegan: Community Ecology Package. Ordination methods, diversity analysis and other functions for community and vegetation ecologists. Version 2.4-4. URL Available online at: https://cran.r-project.org/package=vegan (accessed June 27, 2022).
Osburn, E. D., McBride, S. G., Aylward, F. O., Badgley, B. D., Strahm, B. D., Knoepp, J. D., et al. (2019). Soil bacterial and fungal communities exhibit distinct long-term responses to disturbance in temperate forests. Front. Microbiol. 10, 872. doi: 10.3389/fmicb.2019.02872
Pankratov, T. A., and Dedysh, S. N. (2010). Granulicella paludicola gen. nov., sp. nov., Granulicella pectinivorans sp. nov., Granulicella aggregans sp. nov. and Granulicella rosea sp. nov., acidophilic, polymer-degrading acidobacteria from Sphagnum peat bogs. Int. J. Syst. Evol. Microbiol. 60, 2951–2959. doi: 10.1099/ijs.0.021824-0
Papp, K., Hungate, B. A., and Schwartz, E. (2020). Glucose triggers strong taxon-specific responses in microbial growth and activity: insights from DNA and RNA qSIP. Ecology 101, e02887. doi: 10.1002/ecy.2887
Peršoh, D., Stolle, N., Brachmann, A., Begerow, D., and Rambold, G. (2018). Fungal guilds are evenly distributed along a vertical spruce forest soil profile while individual fungi show pronounced niche partitioning. Mycologic. Progr. 17, 925–939. doi: 10.1007/s11557-018-1405-6
Ponge, J.-F. (2003). Humus forms in terrestrial ecosystems: a framework to biodiversity. Soil Biol. Biochemistr. 35, 935–945. doi: 10.1016/S0038-0717(03)00149-4
Prada-Salcedo, L. D., Goldmann, K., Heintz-Buschart, A., Reitz, T., Wambsganss, J., Bauhus, J., et al. (2021b). Fungal guilds and soil functionality respond to tree community traits rather than to tree diversity in European forests. Mol. Ecol. 30, 572–591. doi: 10.1111/mec.15749
Prada-Salcedo, L. D., Prieto, C., and Correa, M. F. (2014). Screening phosphate solubilizing actinobacteria isolated from the rhizosphere of wild plants from the Eastern Cordillera of the Colombian Andes. Afric. J. Microbiol. Res. 8, 734–742. doi: 10.5897/AJMR2013.5940
Prada-Salcedo, L. D., Wambsganss, J., Bauhus, J., Buscot, F., and Goldmann, K. (2021a). Low root functional dispersion enhances functionality of plant growth by influencing bacterial activities in European forest soils. Environ. Microbiol. 23, 1889–1906. doi: 10.1111/1462-2920.15244
Prescott, C. E., and Grayston, S. J. (2013). Tree species influence on microbial communities in litter and soil: current knowledge and research needs. For. Ecol. Manage. 309, 19–27. doi: 10.1016/j.foreco.2013.02.034
Qu, J.-H., and Yuan, H.-L. (2008). Sediminibacterium salmoneum gen. nov., sp. nov., a member of the phylum Bacteroidetes isolated from sediment of a eutrophic reservoir. Int. J. Systemat. Evolut. Microbiol. 58, 2191–2194. doi: 10.1099/ijs.0.65514-0
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Retallack, G. J. (2005). “SOILS | Modern,” in Encyclopedia of Geology, eds Selley, R.C., Cocks, L.R.M., and Plimer, I.R. (Oxford: Elsevier), pp. 194–202.
Richter, D. D., and Markewitz, D. (1995). How Deep Is Soil?: Soil, the zone of the earth's crust that is biologically active, is much deeper than has been thought by many ecologists. BioScience 45, 600–609. doi: 10.2307/1312764
Rinnan, R., and Bååth, E. (2009). Differential utilization of carbon substrates by bacteria and fungi in tundra soil. Appl. Environ. Microbiol. 75, 3611–3620. doi: 10.1128/AEM.02865-08
Salem, H. M., Abdel-Salam, A., Abdel-Salam, M. A., and Seleiman, M. F. (2017). “Soil xenobiotics and their phyto-chemical remediation,” in Xenobiotics in the Soil Environment: Monitoring, Toxicity and Management, eds M.Z. Hashmi, V. Kumar and A. Varma (Cham: Springer International Publishing), pp. 267–280.
Sharrar, A. M., Crits-Christoph, A., Méheust, R., Diamond, S., Starr, E. P., and Banfield, J. F. (2020). Bacterial secondary metabolite biosynthetic potential in soil varies with phylum, depth, and vegetation type. mBio 11, e00416–00420. doi: 10.1128/mBio.00416-20
Sun, S., Li, S., Avera, B. N., Strahm, B. D., and Badgley, B. D. (2017). Soil Bacterial and Fungal communities show distinct recovery patterns during forest ecosystem restoration. Appl. Environ. Microbiol. 83, e00966–e00917. doi: 10.1128/AEM.00966-17
Suzuki, S. U., and Sasaki, A. (2019). Ecological and evolutionary stabilities of biotrophism, necrotrophism, and saprotrophism. Am. Nat. 194, 90–103. doi: 10.1086/703485
Swift, M. J., Heal, O. W., and Anderson, J. M. (1979). Decomposition in Terrestrial Ecosystems. Oxford: Blackwell.
Thompson, J., Johansen, R., Dunbar, J., and Munsky, B. (2019). Machine learning to predict microbial community functions: an analysis of dissolved organic carbon from litter decomposition. PLoS One 14, e0215502. doi: 10.1371/journal.pone.0215502
Thoms, C., Gattinger, A., Jacob, M., Thomas, F. M., and Gleixner, G. (2010). Direct and indirect effects of tree diversity drive soil microbial diversity in temperate deciduous forest. Soil Biol. Biochemistr. 42, 1558–1565. doi: 10.1016/j.soilbio.2010.05.030
Thorley, R. M., Taylor, L. L., Banwart, S. A., Leake, J. R., and Beerling, D. J. (2015). The role of forest trees and their mycorrhizal fungi in carbonate rock weathering and its significance for global carbon cycling. Plant Cell Environ. 38, 1947–1961. doi: 10.1111/pce.12444
Tripathi, B. M., Kim,1, H. M., Jung, J. Y., Nam, S., Ju, H. T., Kim, M., et al. (2019). Distinct taxonomic and functional profiles of the microbiome associated with different soil horizons of a moist tussock Tundra in Alaska. Front. Microbiol. 10, 1442. doi: 10.3389/fmicb.2019.01442
Uri, V., Varik, M., Aosaar, J., Kanal, A., Kukumägi, M., and Lõhmus, K. (2012). Biomass production and carbon sequestration in a fertile silver birch (Betula pendula Roth). forest chronosequence. For. Ecol. Manage. 267, 117–126. doi: 10.1016/j.foreco.2011.11.033
Uroz, S., Calvaruso, C., Turpault, M.-P., and Frey-Klett, P. (2009). Mineral weathering by bacteria: ecology, actors and mechanisms. Trends Microbiol. 17, 378–387. doi: 10.1016/j.tim.2009.05.004
Uroz, S., Ioannidis, P., Lengelle, J., Cébron, A., Morin, E., Buée, M., et al. (2013). Functional assays and metagenomic analyses reveals differences between the microbial communities inhabiting the soil horizons of a Norway spruce plantation. PloS one 8, e55929–e55929. doi: 10.1371/journal.pone.0055929
Uroz, S., Turpault, M., Delaruelle, C., Mareschal, L., Pierrat, J., and Frey-Klett, P. (2012). Minerals affect the specific diversity of forest soil bacterial communities. Geomicrobiol. J. 29, 88–98. doi: 10.1080/01490451.2010.523764
Van Der Heijden, M. G. A., Bardgett, R. D., and Van Straalen, N. M. (2008). The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol. Lett. 11, 296–310. doi: 10.1111/j.1461-0248.2007.01139.x
Vesterdal, L., Schmidt, I. K., Callesen, I., Nilsson, L. O., and Gundersen, P. (2008). Carbon and nitrogen in forest floor and mineral soil under six common European tree species. For. Ecol. Manage. 255, 35–48. doi: 10.1016/j.foreco.2007.08.015
Weißbecker, C., Schnabel, B., and Heintz-Buschart, A. (2020). Dadasnake, a Snakemake implementation of DADA2 to process amplicon sequencing data for microbial ecology. GigaScience 9, 135. doi: 10.1093/gigascience/giaa135
Wemheuer, F., Taylor, J. A., Daniel, R., Johnston, E., Meinicke, P., Thomas, T., et al. (2020). Tax4Fun2: prediction of habitat-specific functional profiles and functional redundancy based on 16S rRNA gene sequences. Environ. Microbiome 15, 11. doi: 10.1186/s40793-020-00358-7
Xia, Z.-F., Guan, T.-W., Ruan, J.-S., Huang, Y., and Zhang, L.-L. (2013). Longimycelium tulufanense gen. nov., sp. nov., a filamentous actinomycete of the family Pseudonocardiaceae. Int. J. Systemat. Evol. Microbiol. 63, 2813–2818. doi: 10.1099/ijs.0.044222-0
Yarwood, S. A., Bottomley, P. J., and Myrold, D. D. (2010). Soil microbial communities associated with douglas-fir and red alder stands at high- and low-productivity forest Sites in Oregon, USA. Microb. Ecol. 60, 606–617. doi: 10.1007/s00248-010-9675-9
Yost, J. L., and Hartemink, A. E. (2020). How deep is the soil studied—an analysis of four soil science journals. Plant and Soil 452, 5–18. doi: 10.1007/s11104-020-04550-z
Yurkov, A., Inácio, J., Chernov, I. Y., and Fonseca, Á. (2015). Yeast biogeography and the effects of species recognition approaches: the case study of widespread basidiomycetous species from birch forests in Russia. Curr. Microbiol. 70, 587–601. doi: 10.1007/s00284-014-0755-9
Zakharova, O. S., Zenova, G. M., and Zvyagintsev, D. G. (2003). Some approaches to the selective isolation of actinomycetes of the genus actinomadura from soil. Microbiology 72, 110–113. doi: 10.1023/A:1022294526830
Keywords: bacterial pathways, deciduous/evergreen, fungal guilds, microbial indicator taxa, random forest, r/K-strategists
Citation: Prada-Salcedo LD, Prada-Salcedo JP, Heintz-Buschart A, Buscot F and Goldmann K (2022) Effects of Tree Composition and Soil Depth on Structure and Functionality of Belowground Microbial Communities in Temperate European Forests. Front. Microbiol. 13:920618. doi: 10.3389/fmicb.2022.920618
Received: 14 April 2022; Accepted: 14 June 2022;
Published: 11 July 2022.
Edited by:
Jonathan M. Adams, Nanjing University, ChinaReviewed by:
Dennis Goss-Souza, Federal Institute of Education, BrazilLin Zhang, China Agricultural University, China
Copyright © 2022 Prada-Salcedo, Prada-Salcedo, Heintz-Buschart, Buscot and Goldmann. 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: Luis Daniel Prada-Salcedo, bHVpcy5zYWxjZWRvJiN4MDAwNDA7dWZ6LmRl