Abstract
The assembly of biological communities depends on deterministic and stochastic processes whose influence varies across spatial and temporal scales. Although ectomycorrhizal (ECM) fungi play a key role in forest ecosystems, our knowledge on ECM community assembly processes and their dependency on spatial scales is still scarce. We analysed the assembly processes operating on ECM fungal communities associated with Cistus albidus L. and Quercus spp. in Mediterranean mixed forests (Southern Spain), for which root tip ECM fungi were characterized by high-throughput sequencing. The relative contribution of deterministic and stochastic processes that govern the ECM fungal community assembly was inferred by using phylogenetic and compositional turnover descriptors across spatial scales. Our results revealed that stochastic processes had a significantly higher contribution than selection on root tip ECM fungal community assembly. The strength of selection decreased at the smallest scale and it was linked to the plant host identity and the environment. Dispersal limitation increased at finer scales, whilst drift showed the opposite pattern likely suggesting a main influence of priority effects on ECM fungal community assembly. This study highlights the potential of phylogeny to infer ECM fungal community responses and brings new insights into the ecological processes affecting the structure and dynamics of Mediterranean forests.
Similar content being viewed by others
Introduction
Understanding the ecological processes behind community organization is key to predict the establishment and coexistence of local species pools (Pearson et al. 2018). They have been classified into selection, dispersal, drift and speciation (Vellend 2010). Speciation operates at a wide temporal scale and it is expected to have little effect in organizing communities able to exchange individuals, i.e. belonging to the same metacommunity. The other processes can be classified into deterministic that may conduct selection from the regional species pool (e.g. species interactions and environmental filters), and stochastic, based on probabilistic events associated with population demography: drift (i.e. populations fluctuating by chance) and dispersal (i.e. the ability of individuals to move across communities) (Stegen et al. 2012; Ning et al. 2020). Dispersal may have contrasting effects: when it is limited, it enhances stochastic drift, whereas high dispersal may cause mass effects that homogenize communities (Stegen et al. 2013; Evans et al. 2017). Deterministic processes fit into the niche theory that assumes that species environmental requirements determine community assembly (Chesson 2000; Goberna et al. 2019; Viana and Chase 2019).
This theoretical scheme is dynamic and the relative importance of deterministic and stochastic processes varies across spatial and temporal scales (Chase and Myers 2011; Kivlin et al. 2014; Davison et al. 2016; Viana and Chase 2019). In fact, the effect of deterministic mechanisms such as environmental filtering is usually observed across broad abiotic gradients. Meanwhile, the importance of biotic interactions, together with stochastic processes, increases as the environment becomes more homogeneous at smaller scales (Götzenberger et al. 2012; Vályi et al. 2016).
Soil microbial communities show complex structural and functional responses to the environment that might hinder our understanding of community assembly processes (Pérez-Izquierdo et al. 2019). Microbial community features, such as phylogenetic relatedness (Miller et al. 2017) or functional traits (Martiny et al. 2013; López-García et al. 2018), are increasingly relevant in studying the assembly of microbial communities. Thus, to infer ecological processes, Stegen et al. (2013) proposed an analytical framework combining both the phylogenetic turnover (i.e. the averaged evolutionary distance among taxa in different communities) and compositional turnover of microbial communities. In fact, both abundance-based and phylogenetic information are currently steadily explored in microbial ecology by using null model approaches (Tripathi et al. 2018; Zhao et al. 2019; Pereira et al. 2020). Null modelling allows estimating the extent to which the structural patterns of the observed community differ from random expectations derived from stochastic processes (Chase and Myers, 2011; Münkemüller et al. 2020). This analytical approach may also link the dependency of assembly processes to spatial or temporal scales (Kivlin et al. 2014; Viana and Chase 2019; Zhao et al. 2019). Thus, Davison et al. (2016) demonstrated the shift from environmental filtering towards biotic interactions from broader to finer spatial scales by studying the phylogenetic structure of arbuscular mycorrhizal fungal communities collected worldwide. However, despite this potential, the null modelling approach has been barely implemented to analyse the scale dependency of ectomycorrhizal (ECM) fungal community assembly rules (but see Pickles et al. 2012).
The ECM symbiosis has key ecological and biogeochemical implications in forests, such as affecting net primary production and promoting nutrient mobilization through belowground mycelial networks (Lu and Hedin 2019; Rog et al. 2020; Tedersoo et al. 2020). The species composition of ECM fungal communities is influenced by geographical and abiotic environmental conditions (e.g. climate, seasonality or soil properties, Rincón et al. 2015; Pérez-Izquierdo et al. 2017; van der Linde et al. 2018), biotic factors such as the host partner (Pérez-Izquierdo et al. 2017; Wang et al. 2019), dispersal ability and fungal life-history differences across taxa (Peay et al. 2012) and the interaction of both biotic and abiotic factors (Buscardo et al. 2010; Põlme et al. 2013). However, the spatial hierarchy of assembly processes governing ECM fungal communities remains unclear (Zhao et al. 2019) and particularly in the Mediterranean area where such kinds of studies are scarce. The Mediterranean basin is a biodiversity hotspot where species have diversified by adapting to particular environmental conditions, such as large soil heterogeneity, severe seasonal droughts or fire recurrence (Rundel et al. 2016; Pérez-Valera et al. 2018; Pérez-Izquierdo et al. 2020). These adaptations have led to highly diverse biological communities with complex assembly processes (Rincón et al. 2014; Alcántara et al. 2018).
Under this environmental context, we aimed to identify the main potential factors and underlying processes explaining ECM fungal community structure in Mediterranean mixed forests. We also wanted to study the habitat scale dependency of these assembly processes. We expected that ECM fungal communities would be deterministically assembled due to the relatively high dispersal potential of this fungal group (Peay et al. 2012) (hypothesis 1). That notwithstanding, and as formerly stated (Stegen et al. 2013; Vellend et al. 2014; Peay 2018), the relative importance of deterministic (environmental filtering) vs. stochastic (dispersal limitation and drift) assembly processes should change with the spatial scale, with stochastic forces increasing in importance at smaller scales (hypothesis 2). To test these hypotheses, we analysed, at different spatial scales, the structure of root tip ECM fungal communities of three dominant plant species of Mediterranean mixed forests in Southern Spain.
Materials and methods
Study site and sampling design
This study was conducted on mixed forests located in two protected areas in Southern Spain (hereafter sites) (Fig. 1): Natural Park of Sierras de Cazorla, Segura y Las Villas (38.29°N, − 2.57°W, hereafter called Segura) and Sierra Sur de Jaén Park (37.64°N, − 3.73°W, hereafter called Jaén). Mixed forests of Pinus halepensis Mill., Quercus ilex L. and Quercus faginea Lam. dominated in Jaén, and mixed forests of Pinus nigra subsp. salzamanii J.F. Arnold, Pinus pinaster Ait., Q. ilex, Q. faginea and Quercus pyrenaica Willd. in Segura (Pulgar et al. 2017; Alcántara et al. 2018). The understory plant species often found was highly diverse, belonging to the genera Cistus, Crataegus, Juniperus, Rosmarinus and Thymus, among others. Both sites show similar geological context (Siles et al. 2010). In the region, the climate is continental Mediterranean, with an annual mean rainfall of 800–1000 mm and a mean annual temperature of 10–12 °C.
Experimental design and sample processing and analyses are detailed in Suppl. Info. Appendix 2. Briefly, three plant species were sampled at both sites (Cistus albidus L., Q. faginea and Q. ilex). We collected ECM root tips from four individuals per plant species in four plots per site that were molecularly analysed (final n = 92) by Illumina Miseq (ITS-1 rDNA region). Gravimetric moisture (GM), pH and organic matter (OM) content of soils were also determined. Raw sequences were processed with DADA2 pipeline v1.16 (Callahan et al. 2016; R Core Team 2021) and the LULU algorithm (Frøslev et al. 2017). Operational taxonomic units were obtained by clustering amplicon sequence variants at a 97% similarity and their taxonomy was checked against UNITE database v7.2 (Abarenkov et al. 2010). OTUs were classified by fungal guilds with the FUNGuild database v1.0 (Nguyen et al. 2016). The final output yielded 449 OTUs and 6,582,941 reads related to the ECM lifestyle. Homogeneous sequencing depth across samples was confirmed by rarefaction (vegan R package, Oksanen et al. 2019) (Fig. S1).
Plant and ECM fungal phylogenies
The phylogeny of plant species used in this study has been previously published by Alcántara et al. (2018). The ECM fungal community phylogeny was estimated with the program Phylomatic as implemented in Phylocom 4.2 (Webb et al. 2008) and the BEAST software v.1.10.4 (Suchard et al. 2018). The phylogenetic tree was generated with a reference fungal mega-tree, whose topology and age for major nodes were inferred from the phylogenetic information available in the literature (Pérez-Izquierdo et al. 2019). The input for Phylomatic was the list of fungal OTUs that was matched to the most resolved position in the mega-tree previously constructed, so that if any genus was missing from the mega-tree, a polytomy of genera within that family was returned (Moles et al. 2005). For the remaining undated nodes, ages were estimated with the BLADJ algorithm in Phylocom (Webb et al. 2008) that distributes undated nodes evenly between nodes of known ages. To check for the robustness of results to the topological and chronological uncertainty introduced by the Phylomatic and BLADJ procedures, an additional phylogenetic tree based on a branch length adjustment procedure that follows a birth–death evolutionary model while randomly resolves the polytomies in the BEAST program was constructed (Drummond and Rambaut 2007; Kuhn et al. 2011). Markov Chain Monte Carlo (MCMC) analyses for 5 × 106 iterations were run, sampling trees every 103 iterations, discarded a 25% burn-in and recovered the maximum clade credibility tree using the TreeAnnotator v1.5.4 software (Drummond and Rambaut 2007).
Data analyses
Drivers of ECM fungal community composition
The distribution of measured soil variables (pH, OM and gravimetric soil moisture previously log transformed) was tested via linear mixed-effect models (LMMs) with host plant species identity and season as fixed factors, and site and plot nested in site as random factors. The interaction host plant × season was not significant and it was discarded from models. The same model was used to analyse the nearest taxon index (NTI) (Webb et al. 2002) as a measure of phylogenetic distance between OTUs in a single community. NTI allows the inference of potential assembly processes, i.e. phylogenetic similarity or clustering of fungi would indicate the action of an environmental filter, while overdispersion that of competition. NTI was quantified with the ses.mntd function and taxa.labels as null model (999 null communities in picante R package) and the relative abundance and the phylogenetic distance matrices of the ECM fungal communities. The significance of fixed and random factors was determined with the Anova (car R package) and rand functions (lmerTest R package) respectively, and their relative contribution by the coefficient of determination (pseudo-R2) calculated with the r.squaredGLMM function (MuMIn R package). Significant differences of predicted means were separated by the Tukey test (HSD.test function, agricolae R package).
Phylogenetic relatedness is a good proxy for shared traits of species, allowing to identify environmental filters driving community assembly (Martiny et al. 2015; Goberna and Verdú 2016; Rog et al. 2020). Thus, we used the β-nearest taxon index (βNTI), which is the between-assemblage analogue of NTI (Fine and Kembel 2011). βNTI was calculated with qpen function and taxa shuffle as a null model (999 null communities in iCAMP R package).
To disentangle the weight of biotic (plant species), abiotic (environmental properties) factors, and stochastic events in driving phylogenetic turnover (βNTI) of the fungal community, we subjected to principal component analysis (PCA) the soil abiotic variables, the spatial variables—previously decomposed via principal coordinates of neighbourhood matrix (PCNM, Borcard and Legendre 2002)—and the two principal coordinates axes of phylogenetic distances among plant species (Stegen et al. 2013). The significant influence of PCA axes on βNTI was tested by forward selection based on distance-based redundancy analysis (dbRDA) (ordistep and capscale functions in vegan R package).
To determine if any environmental or spatial variable could be driving species composition independently of their phylogeny, we applied the approach described for βNTI to the Raup Crick metric based on Bray–Curtis dissimilarity (RCBray) that measures the turnover of OTUs composition against that expected by chance (Stegen et al. 2013). This metric was calculated with the qpen function and taxa shuffle as the null model (999 null communities in the iCAMP R package).
Finally, to check the fungal taxonomic groups related to the community phylogenetic turnover, we calculated the principal coordinate axes (PCoA) of βNTI (function pcoa in the ape R package). PCoA axes were correlated with the significant PCA axes, derived from the forward selection, by applying the Spearman method with Bonferroni adjustment, using the corr.test function in the psych R package. PCoAs significantly related to PCA axes were tested against the relative abundance of main fungal phyla (Ascomycetes, Basidiomycetes) and families, following the same correlation analysis procedure as described above.
The main trends of ECM fungal community composition were represented by using a non-metric multidimensional scaling (nMDS) (Bray–Curtis dissimilarity, metaMDS function, vegan R package), where host plant species and site were plotted together with the abundance of ECM fungal taxa (envfit function 999 permutations, vegan R package).
Ecological assembly processes acting on ECM fungal communities across spatial scales
To estimate the relative importance of the different assembly processes operating in the ECM fungal communities (hypothesis 2), we followed the analytical approach described by Stegen et al. (2012, 2013) (Fig. 2). We first discriminate among deterministic (i.e. selection) and stochastic processes by the βNTI values calculated for each pair of samples. When |βNTI|> 2, i.e. when the phylogenetic turnover is significantly smaller (βNTI < − 2) or greater (βNTI > 2) than expected by chance, it indicates that community assembly is affected by habitat filtering processes (selection). When |βNTI|< 2 (i.e. − 2 < βNTI < + 2), βNTI values are non-different from random expectations showing communities mainly driven by dispersal and/or drift processes (stochasticity) (sensu Stegen et al. 2013). To more finely delimit the stochastic processes, the matrix of RCBray values was calculated with the qpen function in iCAMP R package. RCBray > 0.95 that did not show |βNTI|> 2 assigned dispersal limitation processes operating in community assembly, RCBray < − 0.95 that did not show |βNTI|> 2 indicated homogenizing dispersal (i.e. high dispersal ability of ECM fungal species) and |RCBray|< 0.95 assigned drift as the determinant of community assembly (Fig. 2, Stegen et al. 2013). Both metrics (βNTI and RCBray) were calculated across spatial/habitat scales (regional, site, plot nested in site and host plant species nested in plot and site) by restricting the input abundance matrices, used for calculating the null communities, to regional, site, plot and host plant species, respectively. The relative contribution of each assembly process at each scale was determined as the number of community-community turnover assigned to each assembly process. A χ2 test was used to detect significant differences in the proportion of assembly processes across spatial scales.
Finally, to seek for the presence of species interactions driving the assembly of ECM fungal communities across scales, NTI values were calculated for each spatial scale, similar to what was done with the previous metrics. The significance of averaged NTI values at each scale respect to the null expectation was assessed by means of t test (p < 0.05) and differences between spatial scales were evaluated by the Tukey test with the function HSD.test in the agricolae R package.
All analyses conducted in this study were carried out with functions and packages of the R free software v.4.1.1 (R Core Team 2021).
Results
The studied sites (Fig. 1) differed in soil properties, with higher pH and OM found in Jaén (pH 7.7 ± 0.1 and OM 16.5 ± 2.1%) than Segura (pH 6.6 ± 0.2 and OM 8.9 ± 1.0%) (Tables 1 and S1). Soil pH varied significantly between host plant species, but differences were mainly driven by samples from Segura, where C. albidus (7.1 ± 0.2) showed higher values than Quercus spp. (6.4 ± 0.1; 6.3 ± 0.2).
The ECM fungal community dissimilarity did also vary across sampling sites and host plant species (Tables 1 and S1; Fig. 3). The nMDS revealed that the abundance of Pezizales significantly correlated with C. albidus in Jaén and that abundance of Russulales and Helotiales did correlate with Quercus spp. in Segura (Fig. 3). NTI values were positive at all spatial scales analysed, and different from null expectations (Tables S1 and S2), pointing to phylogenetic clustering of ECM fungal communities.
Drivers of ECM fungal community composition
Two PCA axes, composed by environmental and spatial variables, were found to drive the phylogenetic composition of ECM fungal communities (βNTI) (Fig. 4; Table 2). Depending of the components of significant PCA axes, they were interpreted as follows: those correlated to soil or plant phylogenetic variables were considered as ‘measured’ environmental filters, while PCA axes correlating with PCNM axes (spatial variables), but not with measured variables, were considered ‘unmeasured’ environmental filters (Table S3). Our results showed that PCA 4 mainly correlated with plant phylogeny (first PCoA axis, differentiating C. albidus from Quercus spp.; see Table S3), hence representing a ‘measured’ environmental filter. PCA 9 did also significantly correlate with βNTI, mainly due to PCNM 6, which represents the narrowest spatial scale and a spatially structured ‘unmeasured’ environmental filter.
βNTI was decomposed into PCoA axes, and those axes significantly correlating with PCA 4 and 9 (Table S4) revealed the variation of certain fungal families. In fact, PCoA 5 and PCoA 36 (related to PCA 4, r = − 0.34, p < 0.001 and r = − 0.24, p = 0.02), correlated with the relative abundance of Gomphidiaceae plus Russulaceae (negatively), and Tuberaceae (positively), respectively (Table S4). By contrast, PCoA 9 (related to PCA 9, r = − 0.21, p = 0.04) did not correlate with any fungal family. In the case of PCoA axes 44 and 57 (related to PCA 9, r = 0.21, p = 0.05, and r = − 0.23, p = 0.03, respectively), only PCoA 57 significantly correlated with the abundance of Helvellaceae (positively) and Pezizaceae (negatively) (Table S4).
Focusing on compositional ECM fungal community turnover (RCBray), whether a PCA axis did significantly drive it, but not that of βNTI, it was interpreted as an indication of dispersal limitation (i.e., limiting species distribution) or, alternatively, as an environmental driver acting on traits non-phylogenetically conserved. The RCBray outcomes were only marginally influenced by PCA 3, which was related to the second PCoA axis of plant phylogeny, i.e. differentiating equally between the three plant species (Table 2).
Scale dependency of ECM fungal community assembly processes
ECM fungal community turnover was explained to a different extent by different assembly processes, depending on the habitat scale (Fig. 5). Homogenizing dispersal was the process explaining the least number of community turnovers and it did not change with scale (χ2 = 0.94, p = 0.82). By contrast, selection (χ2 = 20.79; p < 0.001), dispersal limitation (χ2 = 163.00; p < 0.001) and drift (χ2 = 124.45; p < 0.001) did show significant variations across habitat scales. Selection explained a small percentage of community turnover and its importance decreased between the regional and the rest of the scales (Fig. 5). Drift decreased in importance when reducing the scale, while dispersal limitation followed the opposite pattern, increasing in importance when reducing the scale (Fig. 4).
Despite the general significant phylogenetic clustering revealed by the NTI analysis at every analysed habitat scale, it showed no variation across habitat scales (Table S2).
Discussion
Drivers of ECM fungal community composition
The fact that phylogenetic clustering was predominant across habitat scales indicates that environmental filtering is a widely distributed assembly force of root tip ECM fungal communities of the Mediterranean mixed forests, as previously reported in other forest ecosystems (Koide et al. 2011; Pena et al. 2017). Contrary to our expectations, we did not observe strong relationships between ECM fungal communities and the studied soil variables. If any, they were probably masked by the great effect of host plant identity (e.g. pH was significantly higher under C. albidus than Quercus spp.). The pH gradient, mediated by the plant species identity, could be a determinant to filter the regional ECM fungal species pool (i.e. Ascomycetes are usually tied to higher soil pH), as other studies have evidenced in soil microbial communities (Vályi et al. 2016; Glassman et al. 2017; Tripathi et al. 2018), and would explain the correlation of this fungal phylum with C. albidus (particularly Pezizales and the family Tuberaceae within this fungal order). Contrary, Gomphidiaceae and Russulaceae (Basidiomycetes) showed the opposite trend. Although the three soil variables analysed are commonly reported as the main environmental predictors of ECM fungal community assembly (Glassman et al. 2017; van der Linde et al. 2018), the phylogenetic turnover observed in our study system revealed that other non-measured spatially structured environmental variables (e.g. temperature, soil texture, quality of organic matter or soil key nutrients) might also explain ECM fungal community outcomes (Pena et al. 2017; van der Linde et al. 2018). This phylogenetic preference points towards the conservation of ecological traits in the evolution of ECM fungi (e.g. production of fruit bodies and spores, dispersal strategies, host preference or mycelium exploration type) (Treseder and Lennon 2015; van der Heijden et al. 2015; van der Linde et al. 2018). These results partially confirmed our first hypothesis: the ECM fungal community significantly depended on the host plant species identity and soil properties.
Despite the significant effect of environmental predictors, our results further revealed that all types of assembly processes described in Vellend (2010) had a role in configuring the ECM fungal communities. The magnitude of the influence of each process did largely vary. In fact, dispersal limitation and drift had a much greater role than deterministic (i.e. selection) processes contrary to our expectations (first hypothesis). This result suggests that stochastic mechanisms, such as priority effects (i.e. arrival timing of species into the community), could be key to understanding ECM community assembly outcomes (Peay 2018). Indeed, it is known that the early arrival of ECM fungal species to colonize root tips can limit the establishment of later arriving taxa, what would explain the role of dispersal limitation in the metacommunity (Kennedy et al. 2009; Napoli et al. 2010; Thoen et al. 2019). However, the lower extent of selection structuring ECM fungal communities could be an artefact of the used methodological approach. This approach relies on the assumption of phylogenetic conservatism of functional traits that, in the case of ECM fungi, could be challenged by their paraphyletic nature (Tedersoo et al. 2010). In any case, the fact that community compositional turnover (measured through RCBRAY) did not reveal any contribution of selection would also confirm the greater role of stochasticity over the selection on the ECM community assembly in our study system.
ECM fungal community assembly processes were habitat scale dependent
The influence of assembly processes revealed differences with the analysed scale confirming our second hypothesis. Drift and dispersal limitation did explain most of this variation and both processes showed opposite patterns, being the former more influent at lower scale (host plant) and the last at higher scale (regional). Our results were consistent with the framework proposed by Zobel (1997) and previously tested by Davison et al. (2016) on arbuscular mycorrhizal fungal communities: the higher importance of dispersal limitation at finer scales, which might be related to the potential main role of priority effects also driven by competition, i.e. first taxa arriving within the system determine the establishment of the later taxa (Pickles et al. 2012). On the other hand, homogenizing dispersal showed a general low contribution to assembly and did not vary across habitat scales, suggesting that communities were not assembled by few highly abundant or highly dispersed taxa (i.e. mass effect) (Evans et al. 2017). As expected, selection generally diminished with habitat scales (Chase 2014; Zhao et al. 2019), likely influenced by those dominant ECM fungal taxa whose relatively high abundance could be reinforced by priority effects (Moeller and Peay 2016). The lack of differences in phylogenetic clustering across habitat scales seems to point out towards a lack of hierarchical effect of spatial scale on biotic interactions governing ECM fungal community assembly that has been previously described for other microbial groups (Götzenberger et al. 2012; Davison et al. 2016; Vályi et al. 2016; Goberna et al. 2019). In any case, the eventual absence of spatial hierarchical assembly might be supported by the functional redundancy of distantly related ECM fungal taxa (i.e. lack of phylogenetic conservatism on functional traits) (Pena et al. 2017). This is particularly important in environments such as the Mediterranean, where inter- and intra-annual environmental variability can lead to increasing stochastic events that would mask hierarchical assembly outcomes within communities. This idea would argue the need to complement the phylogenetic information obtained from the ECM fungal communities studied with other aspects of functional diversity (e.g. fungal traits, Põlme et al. 2021) in order to reveal the mechanisms underlying ECM fungal community assembly in Mediterranean forest ecosystems.
Conclusions
In this study, we have found that (1) spatial factors and host plant species are determinants of ECM fungal community assembly in Mediterranean mixed forests; (2) phylogenetic and compositional turnover are good community structure proxies to evaluate the contribution of deterministic and stochastic assembly processes in ECM fungal communities and (3) the contribution of assembly processes is habitat scale dependent. Thus, ECM community assembly in Mediterranean forests is dependent on both deterministic and stochastic processes, particularly dispersal limitation and drift, increasing the first and diminishing the second at smaller scales. These patterns confirm the habitat scale dependency of assembly processes, as previous studies described in other biological guilds but not in ECM fungal communities. Further studies are needed to disentangle the role of fine-tuned biotic interactions in forests, as well as how ECM phylogeny complements with functional diversity and environmental drivers, particularly in Mediterranean ecosystems. Additional information may be found in the online version of this article at the publisher’s website.
Data accessibility
Data are deposited in the Sequence Read Archive https://www.ncbi.nlm.nih.gov/bioproject/PRJNA787911.
References
Abarenkov K, Henrik Nilsson R, Larsson KH et al (2010) The UNITE database for molecular identification of fungi–recent updates and future perspectives. New Phytol 186:281–285. https://doi.org/10.1111/j.1469-8137.2009.03160.x
Alcántara JM, Pulgar M, Trøjelsgaard K et al (2018) Stochastic and deterministic effects on interactions between canopy and recruiting species in forest communities. Funct Ecol 32:2264–2274. https://doi.org/10.1111/1365-2435.13140
Borcard D, Legendre P (2002) All-scale spatial analysis of ecological data by means of principal coordinates of neighbour matrices. Ecol Modell 153:51–68. https://doi.org/10.1016/S0304-3800(01)00501-4
Buscardo E, Rodríguez-Echeverría S, Martín MP et al (2010) Impact of wildfire return interval on the ectomycorrhizal resistant propagules communities of a Mediterranean open forest. Fungal Biol 114:628–636. https://doi.org/10.1016/j.funbio.2010.05.004
Callahan BJ, McMurdie PJ, Rosen MJ et al (2016) DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods 13:581–583. https://doi.org/10.1038/nmeth.3869
Chase JM (2014) Spatial scale resolves the niche versus neutral theory debate. J Veg Sci 25:319–322. https://doi.org/10.1111/jvs.12159
Chase JM, Myers JA (2011) Disentangling the importance of ecological niches from stochastic processes across scales. Philos Trans R Soc B Biol Sci 366:2351–2363. https://doi.org/10.1098/rstb.2011.0063
Chesson P (2000) Mechanisms of maintenance of species diversity. Annu Rev Ecol Syst 31:343–366. https://doi.org/10.1146/annurev.ecolsys.31.1.343
Davison J, Moora M, Jairus T et al (2016) Hierarchical assembly rules in arbuscular mycorrhizal (AM) fungal communities. Soil Biol Biochem 97:63–70. https://doi.org/10.1016/j.soilbio.2016.03.003
Drummond AJ, Rambaut A (2007) BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol 7:214. https://doi.org/10.1186/1471-2148-7-214
Evans S, Martiny JBH, Allison SD (2017) Effects of dispersal and selection on stochastic assembly in microbial communities. ISME J 11:176–185. https://doi.org/10.1038/ismej.2016.96
Fine PV, Kembel SW (2011) Phylogenetic community structure and phylogenetic turnover across space and edaphic gradients in western Amazonian tree communities. Ecography 34:552–565. https://doi.org/10.1111/j.1600-0587.2010.06548.x
Frøslev TG, Kjøller R, Bruun HH et al (2017) Algorithm for post-clustering curation of DNA amplicon data yields reliable biodiversity estimates. Nature Commun 8:1–11. https://doi.org/10.1038/s41467-017-01312-x
Glassman SI, Wang IJ, Bruns TD (2017) Environmental filtering by pH and soil nutrients drives community assembly in fungi at fine spatial scales. Mol Ecol 26:6960–6973. https://doi.org/10.1111/mec.14414
Goberna M, Montesinos-Navarro A, Valiente-Banuet A et al (2019) Incorporating phylogenetic metrics to microbial co-occurrence networks based on amplicon sequences to discern community assembly processes. Mol Ecol Resour 19:1552–1564. https://doi.org/10.1111/1755-0998.13079
Goberna M, Verdú M (2016) Predicting microbial traits with phylogenies. ISME J 10:959–967. https://doi.org/10.1038/ismej.2015.171
Götzenberger L, de Bello F, Bråthen KA et al (2012) Ecological assembly rules in plant communities-approaches, patterns and prospects. Biol Rev 87:111–127. https://doi.org/10.1111/j.1469-185X.2011.00187.x
Kennedy PG, Peay KG, Bruns TD (2009) Root tip competition among ectomycorrhizal fungi: are priority effects a rule or an exception? Ecology 90:2098–2107. https://doi.org/10.1890/08-1291.1
Kivlin SN, Winston GC, Goulden ML, Treseder KK (2014) Environmental filtering affects soil fungal community composition more than dispersal limitation at regional scales. Fungal Ecol 12:14–25. https://doi.org/10.1016/j.funeco.2014.04.004
Koide RT, Fernandez C, Petprakob K (2011) General principles in the community ecology of ectomycorrhizal fungi. Ann for Sci 68:45–55. https://doi.org/10.1007/s13595-010-0006-6
Kuhn TS, Mooers A, Thomas GH (2011) A simple polytomy resolver for dated phylogenies. Methods Ecol Evol 2:427–436. https://doi.org/10.1111/j.2041-210X.2011.00103.x
López-García Á, Gil-Martínez M, Navarro-Fernández CM et al (2018) Functional diversity of ectomycorrhizal fungal communities is reduced by trace element contamination. Soil Biol Biochem 121:202–211. https://doi.org/10.1016/j.soilbio.2018.03.021
Lu M, Hedin LO (2019) Global plant–symbiont organization and emergence of biogeochemical cycles resolved by evolution-based trait modelling. Nat Ecol Evol 3:239–250. https://doi.org/10.1038/s41559-018-0759-0
Martiny AC, Treseder K, Pusch G (2013) Phylogenetic conservatism of functional traits in microorganisms. ISME J 7:830–838. https://doi.org/10.1038/ismej.2012.160
Martiny JBH, Jones SE, Lennon JT, Martiny AC (2015) Microbiomes in light of traits: a phylogenetic perspective. Science 350. https://doi.org/10.1126/science.aac9323
Miller ET, Farine DR, Trisos CH (2017) Phylogenetic community structure metrics and null models: a review with new methods and software. Ecography 40:461–477. https://doi.org/10.1111/ecog.02070
Moeller HV, Peay KG (2016) Competition-Function Tradeoffs in Ectomycorrhizal Fungi Peerj 4:e2270. https://doi.org/10.7717/peerj.2270
Moles AT, Ackerly DD, Webb CO et al (2005) Factors that shape seed mass evolution. Proc Natl Acad Sci 102:10540–10544. https://doi.org/10.1073/pnas.0501473102
Münkemüller T, Gallien L, Pollock LJ et al (2020) Dos and don’ts when inferring assembly rules from diversity patterns. Glob Ecol Biogeogr 29:1212–1229. https://doi.org/10.1111/geb.13098
Napoli C, Mello A, Borra A et al (2010) Tuber melanosporum, when dominant, affects fungal dynamics in truffle grounds. New Phytol 185:237–247. https://doi.org/10.1111/j.1469-8137.2009.03053.x
Ning D, Yuan M, Wu L et al (2020) A quantitative framework reveals ecological drivers of grassland microbial community assembly in response to warming. Nature Commun 11:1–12. https://doi.org/10.1038/s41467-020-18560-z
Nguyen NH, Song Z, Bates ST et al (2016) FUNGuild: an open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol 20:241–248. https://doi.org/10.1016/J.FUNECO.2015.06.006
Oksanen J, Blanchet FG, Friendly M et al (2019) vegan: community ecology package. R package version 2.5–4. https://CRAN.R-project.org/package=vegan
Pearson DE, Ortega YK, Eren Ö, Hierro JL (2018) Community assembly theory as a framework for biological invasions. Trends Ecol Evol 33:313–325. https://doi.org/10.1016/j.tree.2018.03.002
Peay KG (2018) Timing of mutualist arrival has a greater effect on Pinus muricata seedling growth than interspecific competition. J Ecol 106:514–523. https://doi.org/10.1111/1365-2745.12915
Peay KG, Schubert MG, Nguyen NH, Bruns TD (2012) Measuring ectomycorrhizal fungal dispersal: macroecological patterns driven by microscopic propagules. Mol Ecol 21:4122–4136. https://doi.org/10.1111/j.1365-294X.2012.05666.x
Pena R, Lang C, Lohaus G et al (2017) Phylogenetic and functional traits of ectomycorrhizal assemblages in top soil from different biogeographic regions and forest types. Mycorrhiza 27:233–245. https://doi.org/10.1007/s00572-016-0742-z
Pereira CMR, López-García Á, Silva DKA et al (2020) Tropical forest type influences community assembly processes in arbuscular mycorrhizal fungi. J Biogeogr 47:434–444. https://doi.org/10.1111/jbi.13739
Pérez-Izquierdo L, Zabal-Aguirre M, Flores-Rentería D et al (2017) Functional outcomes of fungal community shifts driven by tree genotype and spatial-temporal factors in Mediterranean pine forests. Environ Microbiol 19:1639–1652. https://doi.org/10.1111/1462-2920.13690
Pérez-Izquierdo L, Zabal-Aguirre M, González-Martínez SC et al (2019) Plant intraspecific variation modulates nutrient cycling through its below ground rhizospheric microbiome. J Ecol 107:1594–1605. https://doi.org/10.1111/1365-2745.13202
Pérez-Izquierdo L, Zabal-Aguirre M, Verdú M et al (2020) Ectomycorrhizal fungal diversity decreases in Mediterranean pine forests adapted to recurrent fires. Mol Ecol 29:1–14. https://doi.org/10.1111/mec.15493
Pérez-Valera E, Verdú M, Navarro-Cano JA, Goberna M (2018) Resilience to fire of phylogenetic diversity across biological domains. Mol Ecol 27:2896–2908. https://doi.org/10.1111/mec.14729
Pickles BJ, Genney DR, Anderson IC, Alexander IJ (2012) Spatial analysis of ectomycorrhizal fungi reveals that root tip communities are structured by competitive interactions. Mol Ecol 21:5110–5123. https://doi.org/10.1111/j.1365-294X.2012.05739.x
Põlme S, Bahram M, Yamanaka T et al (2013) Biogeography of ectomycorrhizal fungi associated with alders (Alnus spp.) in relation to biotic and abiotic variables at the global scale. New Phytol 198:1239–1249. https://doi.org/10.1111/nph.12170
Põlme S, Abarenkov K, Henrik Nilsson R et al (2021) Correction to: FungalTraits: a user friendly traits database of fungi and fungus-like stramenopiles. Fungal Diver 107:129–132. https://doi.org/10.1007/s13225-021-00470-0
Pulgar M, Alcántara JM, Rey PJ (2017) Effects of sampling effort on estimates of the structure of replacement networks. J Veg Sci 28:445–457. https://doi.org/10.1111/jvs.12492
R Core Team (2021) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/
Rincón A, Santamaría BP, Ocaña L, Verdú M (2014) Structure and phylogenetic diversity of post-fire ectomycorrhizal communities of maritime pine. Mycorrhiza 24:131–141. https://doi.org/10.1007/s00572-013-0520-0
Rincón A, Santamaría-Pérez B, Rabasa SG 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. https://doi.org/10.1111/1462-2920.12894
Rog I, Rosenstock NP, Körner C, Klein T (2020) Share the wealth: trees with greater ectomycorrhizal species overlap share more carbon. Mol Ecol 29:2321–2333. https://doi.org/10.1111/mec.15351
Rundel PW, Arroyo MT, Cowling RM et al (2016) Mediterranean biomes: evolution of their vegetation, floras, and climate. Annu Rev Ecol Evol Syst 47:383–407. https://doi.org/10.1146/annurev-ecolsys-121415-032330
Siles G, Rey PJ, Alcántara JM et al (2010) Effects of soil enrichment, watering and seedling age on establishment of Mediterranean woody species. Acta Oecologica 36:357–364
Stegen JC, Lin X, Fredrickson JK et al (2013) Quantifying community assembly processes and identifying features that impose them. ISME J 7:2069–2079. https://doi.org/10.1038/ismej.2013.93
Stegen JC, Lin X, Konopka AE, Fredrickson JK (2012) Stochastic and deterministic assembly processes in subsurface microbial communities. ISME J 6:1653–1664. https://doi.org/10.1038/ismej.2012.22
Suchard MA, Lemey P, Baele G et al (2018) Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol 4:1–5. https://doi.org/10.1093/ve/vey016
Tedersoo L, Bahram M, Zobel M (2020) How mycorrhizal associations drive plant population and community biology. Science 367. https://orcid.org/0000-0001-7957-6704
Tedersoo L, May TW, Smith ME (2010) Ectomycorrhizal lifestyle in fungi: global diversity, distribution, and evolution of phylogenetic lineages. Mycorrhiza 20:217–263. https://doi.org/10.1126/science.aba1223
Thoen E, Aas AB, Vik U et al (2019) A single ectomycorrhizal plant root system includes a diverse and spatially structured fungal community. Mycorrhiza 29:167–180. https://doi.org/10.1007/s00572-019-00889-z
Treseder KK, Lennon JT (2015) Fungal traits that drive ecosystem dynamics on land. Microbiol Mol Biol Rev 79:243–262. https://doi.org/10.1128/MMBR.00001-15
Tripathi BM, Stegen JC, Kim M et al (2018) Soil pH mediates the balance between stochastic and deterministic assembly of bacteria. ISME J 12:1072–1083. https://doi.org/10.1038/s41396-018-0082-4
Vályi K, Mardhiah U, Rillig MC, Hempel S (2016) Community assembly and coexistence in communities of arbuscular mycorrhizal fungi. ISME J 10:2341–2351. https://doi.org/10.1038/ismej.2016.46
Van der Heijden MG, Martin FM, Selosse MA, Sanders IR (2015) Mycorrhizal ecology and evolution: the past, the present, and the future. New Phytol 205:1406–1423. https://doi.org/10.1111/nph.13288
Van der Linde S, Suz LM, Orme CDL et al (2018) Environment and host as large-scale controls of ectomycorrhizal fungi. Nature 558:243–248. https://doi.org/10.1038/s41586-018-0189-9
Vellend M (2010) Conceptual synthesis in community ecology. Q Rev Biol 85:183–206. https://doi.org/10.1086/652373
Vellend M, Srivastava DS, Anderson KM et al (2014) Assessing the relative importance of neutral stochasticity in ecological communities. Oikos 123:1420–1430. https://doi.org/10.1111/oik.01493
Viana DS, Chase JM (2019) Spatial scale modulates the inference of metacommunity assembly processes. Ecology 100:1–9. https://doi.org/10.1002/ecy.2576
Wang YL, Gao C, Chen L et al (2019) Host plant phylogeny and geographic distance strongly structure Betulaceae-associated ectomycorrhizal fungal communities in Chinese secondary forest ecosystems. FEMS Microbiol Ecol 95:1–15. https://doi.org/10.1093/femsec/fiz037
Webb CO, Ackerly DD, Kembel SW (2008) Phylocom: software for the analysis of phylogenetic community structure and trait evolution. Bioinformatics 24:2098–2100. https://doi.org/10.1093/bioinformatics/btn358
Webb CO, Ackerly DD, McPeek MA, Donoghue MJ (2002) Phylogenies and community ecology. Annu Rev Ecol Syst 33:475–505. https://doi.org/10.1146/annurev.ecolsys.33.010802.150448
Zhao J, Gao Q, Zhou J et al (2019) The scale dependence of fungal community distribution in paddy soil driven by stochastic and deterministic processes. Fungal Ecol 42:100856. https://doi.org/10.1016/j.funeco.2019.07.010
Zobel M (1997) The relative role of species pools in determining plant species richness: an alternative explanation of species coexistence? Trends Ecol Evol 12:266–269. https://doi.org/10.1016/S0169-5347(97)01096-3
Acknowledgements
We gratefully acknowledge L. Pomarède and R. Gómez for their help in fieldwork, L. López for her support in laboratory and L. Auer and M. Buée for bioinformatics support.
Funding
Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature. This work was supported by the project COEXMED-II (grant CGL2015-69118-C2-2P) founded by the Spanish Ministry for Economy and Competitiveness (MINECO). J. P. R. held a pre-doctoral fellowship awarded by MINECO (grant BES-2016–078055). A. L. G was supported by the University of Jaén under the Plan 6-UJA post-doctoral fellowship.
Author information
Authors and Affiliations
Contributions
J. M. A., C. A. and A. R. designed the research; J. P. R., L. P. I., J. L. G., J. M. A, C. A. and A.R. performed sampling; J. P. R. did laboratory work; J. P. R, A. L. G., L. P. I. and A. R. analysed data and J. P. R. wrote the manuscript. All authors revised the manuscript.
Corresponding author
Ethics declarations
Conflict of interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Below is the link to the electronic supplementary material.
Additional information may be found in the online version of this article at the publisher’s website:
Appendix 1
. Supporting information with figures and tables (DOCX 118 KB)
Appendix 2
. Supporting information related to Experimental Procedure (DOCX 36 KB)
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
J., PR., J. L., G., L., PI. et al. Scale dependency of ectomycorrhizal fungal community assembly processes in Mediterranean mixed forests. Mycorrhiza 32, 315–325 (2022). https://doi.org/10.1007/s00572-022-01083-4
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00572-022-01083-4