Community Assembly of Endophytic Fungi in Ectomycorrhizae of Betulaceae Plants at a Regional Scale

The interaction between aboveground and belowground biotic communities drives community assembly of plants and soil microbiota. As an important component of belowground microorganisms, root-associated fungi play pivotal roles in biodiversity maintenance and community assembly of host plants. The Betulaceae plants form ectomycorrhizae with soil fungi and widely distribute in various ecosystems. However, the community assembly of endophytic fungi in ectomycorrhizae is less investigated at a large spatial scale. Here, we examined the endophytic fungal communities in ectomycorrhizae of 22 species in four genera belonging to Betulaceae in Chinese forest ecosystems, using Illumina Miseq sequencing of internal transcribed spacer 2 amplicons. The relative contribution of host phylogeny, climate and soil (environmental filtering) and geographic distance (dispersal limitation) on endophytic fungal community was disentangled. In total, 2,106 endophytic fungal operational taxonomic units (OTUs) were obtained at a 97% sequence similarity level, dominated by Leotiomycetes, Agaricomycetes, Eurotiomycetes, and Sordariomycetes. The endophytic fungal OTU richness was significantly related with host phylogeny, geographic distance, soil and climate. The endophytic fungal community composition was significantly affected by host phylogeny (19.5% of variation explained in fungal community), geographic distance (11.2%), soil (6.1%), and climate (1.4%). This finding suggests that environmental filtering by plant and abiotic variables coupled with dispersal limitation linked to geographic distance determines endophytic fungal community assembly in ectomycorrhizae of Betulaceae plants, with host phylogeny being a stronger determinant than other predictor variables at the regional scale.

Introduction

The interaction between aboveground and belowground biotic communities drives community dynamics and ecosystem functioning (Hooper et al., 2000; Wardle et al., 2004; Van Dam and Heil, 2011). As an important component of belowground microorganisms, root-associated fungi exhibit a wide range of strategies including mycorrhizal, endophytic, and pathogenic species (Peay et al., 2008), and these fungi can influence plant growth and fitness, plant community dynamics and biogeochemical cycling (Landeweert et al., 2001; van der Heijden et al., 2008; Liang et al., 2016; Almario et al., 2017). Therefore, revealing mechanisms underlying root-associated fungal community assembly will greatly improve our understanding of biodiversity maintenance and community assembly of plants.

Environmental filtering (niche process) and dispersal limitation (neutral process) are two of the fundamental processes in structuring biotic community assembly (Cottenie, 2005). The effect of environmental filtering by host plants (e.g., identity and phylogeny) on root-associated fungal communities has been widely documented in previous studies (e.g., Ishida et al., 2007; Tedersoo et al., 2013; van der Linde et al., 2018; Wu et al., 2018; Schroeder et al., 2019; Sepp et al., 2019; Wang et al., 2019a, b). For example, plant identity significantly shaped root-associated ectomycorrhizal (EM) fungal communities (Ishida et al., 2007; Tedersoo et al., 2012; van der Linde et al., 2018) and arbuscular mycorrhizal (AM) fungal communities (Martinez-Garcia et al., 2015; Ciccolini et al., 2016; Zheng et al., 2016; Sepp et al., 2019) through host specificity, generating diverse substrates and changing microhabitats (Wardle, 2006; Dickie, 2007; Aponte et al., 2010). Likewise, an increasing amount of studies have indicated that host phylogeny strongly predicts root-associated fungal communities, such as total (Wehner et al., 2014; Schroeder et al., 2019), EM (Põlme et al., 2013; Tedersoo et al., 2013; Wu et al., 2018; Wang et al., 2019a), AM (Wang et al., 2019b), pathogenic (Schroeder et al., 2019; Wang et al., 2019b), and saprotrophic (Schroeder et al., 2019) fungi. According to phylogenetic niche conservatism, closely related plants are more similar in morphological and functional traits than distantly related ones (Losos, 2008), and maybe tend to associate with more similar fungal partners.

Environmental filtering by abiotic variables (e.g., soil and climate) has also been found to impact root-associated fungal communities (e.g., Jarvis et al., 2013; Gao et al., 2015; Miyamoto et al., 2015; Mundra et al., 2016; Barnes et al., 2018; Li et al., 2018; Van Geel et al., 2018; Wu et al., 2018; Wang et al., 2019a). For instance, climate differences explained the largest variation of EM fungal community composition in two Japanese mountains (Miyamoto et al., 2015), while soil variables were the most important drivers of AM fungal community at the Europe (Van Geel et al., 2018). In addition, with increasing nitrogen (N) addition, the abundance of fungal species changed (Cox et al., 2010). This could be ascribable to the fact that fungal taxa commonly prefer to certain niche properties and differ in response to environmental changes (Cox et al., 2010; Mcguire et al., 2010; Mundra et al., 2016; Barnes et al., 2018; Miyamoto et al., 2018).

Generally, geographic distance is widely hypothesized to predict biotic community distribution as the stochastic dispersal of individuals from one site to another could be hindered by the spatial distance between sites, suggesting the effect of dispersal limitation on community assembly (Hubbell, 2001). Many previous studies have demonstrated that dispersal limitation strongly determines community assembly of root-associated fungi, such as total fungi (Wehner et al., 2014; Beck et al., 2015), AM fungi (Dumbrell et al., 2010; Davison et al., 2015; Van Geel et al., 2018), and EM fungi (Tedersoo et al., 2012; Põlme et al., 2013; Matsuoka et al., 2016; van der Linde et al., 2018; Wu et al., 2018; Wang et al., 2019a) at various spatial scales and in different ecosystems. Besides the geographic distance, fungal taxa commonly show varying functional traits, such as dispersal, germination, and colonization abilities, which also contribute to dispersal limitation (Nara, 2009; Peay et al., 2012; Kivlin et al., 2014). Additionally, priority effect could also give rise to the variation in fungal community, that is, early-arriving fungal species could colonize much more habitats and resources, and being competitive dominant over the later ones (Kennedy and Bruns, 2005; Kennedy et al., 2009). However, most studies on root-associated fungal community assembly have mainly focused on mycorrhizal fungi (e.g., Põlme et al., 2013; van der Linde et al., 2018; Van Geel et al., 2018; Wu et al., 2018; Sepp et al., 2019; Wang et al., 2019a). In contrast, the community assembly of root endophytic fungi has less been investigated at a large spatial scale.

The Betulaceae plants are typical EM hosts, widely distribute in Chinese forest ecosystems, and have importantly economic and ecological values (Li and Skvortsov, 1999). In addition to EM fungi in ectomycorrhizae, there are also endophytic fungi which can live in plant tissues without causing apparent harm to the host plants (Petrini, 1991). Recently, we have revealed that both environmental filtering (i.e., host phylogeny, soil and climate) and dispersal limitation (i.e., geographic distance) shaped EM fungal community assembly of Betulaceae plants in Chinese forest ecosystems, and environmental filtering by host plant phylogeny was the strongest determinant (Wang et al., 2019a). As previous studies showed that there was high endophytic fungal diversity in ectomycorrhizae of host plants (Gao et al., 2013, 2015; Matsuoka et al., 2016; Wu et al., 2018; Wang et al., 2019a), we analyzed the community assembly of endophytic fungi in ectomycorrhizae associated with Betulaceae plants in this study. We hypothesize that (1) environmental filtering and dispersal limitation affect endophytic fungal community and (2) environmental filtering by host plant phylogeny has a stronger effect on endophytic fungal community than other variables at a large scale.

Materials and Methods

Study Site, Sampling, and Molecular Analysis

The detailed description of study site, sampling and molecular analysis can be found in Wang et al. (2019a). Briefly, 300 root samples of 23 plant species belonging to four genera (Alnus, Betula, Carpinus, and Corylus) of Betulaceae and 33 rhizosphere soil samples were collected from 26 forest sites in China, extending from 24°30′ N to 52°56′ N, and 94°12′ E to 128°52′ E (Supplementary Figure 1). The geographic coordinates (latitude and longitude) and altitude of each site were recorded using a M-241 GPS instrument (Holux Technology, Inc., Taiwan, China). Mean annual temperature (MAT) and mean annual precipitation (MAP) of each site were compiled from the WorldClim global data set at a resolution of 2.5 min, based on 1970–2000 MAT and MAP (Hijmans et al., 2005). Information on the selected plant species, soil properties, geographical location and climate is available in Wang et al. (2019a).

In total, 44,400 EM root tips from 296 samples (four samples with less than 150 EM root tips were discarded) were used for Illumina MiSeq sequencing of the fungal internal transcribed spacer 2 (ITS2) region of the rDNA (see Wang et al., 2019a). A total of 6,482,896 high-quality ITS2 sequences were obtained from 8,321,293 raw sequences after quality control, and were grouped into 5,763 non-singleton operational taxonomic units (OTUs) at a 97% similarity level. Of these 5,763 OTUs, 4,770 (6,286,352 reads) were identified as fungi based on search against the international nucleotide sequence databases collaboration (INSDC) and the unified system for the DNA-based fungal species linked to the classification (UNITE) databases (Kõljalg et al., 2013) by using the basic local alignment search tool (BLAST) (Altschul et al., 1990). Among the fungal data, 2,667 OTUs (846,012) were endophytic fungi using for subsequent analysis. To eliminate the effect of different sequence numbers on the endophytic fungal community analysis, the number of sequences per sample was normalized to the smallest sample size using the sub. sample command in MOTHUR 1.31.2 (Schloss et al., 2009). The representative endophytic fungal OTU sequences have been submitted to the European Nucleotide Archive (ENA) under study accession nos. LR586058-LR588163. Detailed information about endophytic fungi is summarized in Supplementary Table 1.

Statistical Analysis

All statistical analyses were conducted in R v.3.5.2 (R Development Core Team, 2018), except for the construction of a plant phylogeny tree. ITS sequences of 22 Betulaceae species (one species was removed as samples of which were discarded after normalization) and Quercus acutissima (used as an outgroup) were downloaded from GenBank or sequenced by ourselves (see Wang et al., 2019a). The plant sequences were firstly aligned using the L-INS-i algorithm in MAFFT v. 7. 215 (Katoh and Standley, 2013), and then used to construct a maximum likelihood phylogenetic tree by adopting a general time reversible model with 1000 bootstrap replicates in MEGA ver. 7.0 (Kumar et al., 2016) (Supplementary Figure 2). Pairwise patristic distances (pairwise sums of the branch lengths connecting two terminal taxa) were calculated based on the phylogenetic tree by using the cophenetic. phylo command in the ape package (Paradis et al., 2004). Principal coordinates analysis (PCoA) was adopted to transform the pairwise patristic distance matrix to phylogenetic eigenvectors by applying the cmdscale command in the vegan package (Oksanen et al., 2013). Spatial eigenvectors of principal coordinate of neighbor matrices (PCNMs) were constructed based on the transformation of geographic distance (latitude and longitude) and only those with positive eigenvalues were retained by using the pcnm command in the PCNM package (Dray et al., 2006). Significant host phylogenetic PCoA and spatial PCNM vectors were forward-selected prior to subsequent statistical analyses using the forward. sel command in the packfor package (Dray et al., 2009).

Rarefaction curves depicting endophytic fungal OTU numbers in each plant genus were obtained by using the specaccum command in the vegan package (Oksanen et al., 2013). The endophytic fungal OTU richness and relative abundances of four most abundant endophytic fungal classes (each > 9% of total endophytic fungal reads) did not satisfy the normality of distribution and homogeneity of variance among the four plant genera before and after logarithmic and square root transformation. Thus, Kruskal–Wallis tests were adopted to examine the effect of plant genus on the endophytic fungal OTU richness and relative abundances of the four abundant endophytic fungal classes. Then Dunn’s tests with Bonferroni adjustment were carried out to detect the significant differences in these respects among the four plant genera using the post hoc.kruskal.dunn. test function in the PMCMR package (Pohlert, 2014). Generalized linear models (GLMs) with Poisson distribution and log link function were constructed to determine the effect of spatial, host, soil and climatic variables on endophytic fungal OTU richness. After fitting the global model, the dredge command in the MuMIn package (Bartón, 2018) was firstly used to build all possible candidate models those derived from the different combinations of variables in the global model. Because the response variable was over-dispersed in the models, QAICc (quasi-likelihood adjusted Akaike Information Criterion) was used for the model selection. In order to account for the uncertainty in model selection, model averaging method was applied based on candidate models within 2 ΔQAICc units from the model with the lowest QAICc value by using the model. avg command in the MuMIn package (Bartón, 2018). The relative importance of each variable was generated by summing the QAICc weights from candidate models where they appeared. The variables with significant relative importance in the averaged model were regarded as the best predictors for endophytic fungal OTU richness.

A distance matrix of endophytic fungal community (Hellinger transformed OTU read data) was constructed by calculating dissimilarities using the Bray–Curtis method (Clarke et al., 2006). Endophytic fungal community dissimilarities were visualized by PCoA ordination using the pcoa command in the ape package (Paradis et al., 2004), and the ellipses reflecting 95% confidence intervals were draw around clusters using the ordiellipse function in the vegan package (Oksanen et al., 2013). Meanwhile, analysis of similarity (ANOSIM) was applied to assess the significance of difference in fungal community composition among four plant genera by using the anosim command in the vegan package (Oksanen et al., 2013). Permutational multivariate analysis of variance (PerMANOVA) was implemented to determine the relative importance of host, spatial, soil and climatic variables on endophytic fungal community composition using the adonis command in the vegan package (Oksanen et al., 2013). We also adopted variation partitioning in redundancy analysis to quantify the contribution of host, spatial, soil and climatic variables to the endophytic fungal community richness and composition by using the varpart command in the vegan package (Oksanen et al., 2013).

In order to eliminate the influence of sites on assessing the plant/fungus preferences, two sites (i.e., Changbaishan and Taiyueshan) with the number of plant species ≥ 4 were selected for the preference analysis according to Toju et al. (2016). Additionally, we also calculated index $H_2^{\prime }$ characterizing interaction specialization at community level and checkerboard scores reflecting co-occurrences of fungal species within community (see Wang et al., 2019a).

Results

Identification of Endophytic Fungi

A total of 2,667 (846,012 reads) endophytic fungal OTUs was obtained from 4,770 (6,286,352 reads) fungal OTUs at 97% sequence similarity level. As the number of endophytic fungal sequences varied from 30 to 12,791 across all samples, a random re-sampling procedure was implemented to construct a subset to a depth of 500 sequences per sample (46 samples with reads less than 500 were excluded), resulting a normalized dataset containing 2,106 endophytic fungal OTUs across 250 samples. The 100 most abundant OTUs accounted for 66.7% of the total endophytic fungal reads (Supplementary Figure 3A). The frequency distribution of endophytic fungal OTUs had a long tail, with 1,842 OTUs occurring in no more than 10 samples (Supplementary Figure 3B). However, the rarefaction curves for each of four plant genera did not reach plateau (Supplementary Figure 4).

Endophytic Fungal OTU Richness

The endophytic fungal OTU richness was 44.3 ± 2.2 in Alnus, 48.1 ± 1.6 in Betula, 47.9 ± 2.1 in Carpinus, and 64.8 ± 4.1 (mean ± SE) in Corylus, respectively. Kruskal–Wallis test indicated that plant genus significantly affected the OTU richness of endophytic fungi (χ² = 21.048, P < 0.001). For example, endophytic fungal OTU richness was significantly higher in Corylus than in other three plant genera, but no significant difference was observed in the three plant genera (Figure 1). The final averaged model based on 11 candidate models, suggested that host phylogeny (PCoA1, PCoA3, PCoA8, and PCoA16), spatial vectors (PCNM1, PCNM2, PCNM3, PCNM7, and PCNM10), MAT and soil total N were the best predictor variables for endophytic fungal OTU richness (Table 1 and Supplementary Table 2). In addition, variation partitioning showed that 35.2% of variation in endophytic fungal OTU richness was explained by host phylogeny (19.0%), spatial distance (18.1%), soil (16.5%) and climate (4.4%), with pure effects of 11.3, 9.6, 0.9, and 1.6%, respectively (Figure 2).

FIGURE 1

Figure 1. Richness of endophytic fungal operational taxonomic units (OTUs) in four plant genera. Bars without shared letters indicate significant differences in richness of fungal OTUs among the four host genera according to Dunn’s tests with Bonferroni adjustment at P < 0.05.

TABLE 1

Table 1. Plant and abiotic predictor variables for the richness of endophytic fungal operational taxonomic units following the model-averaging procedure.

FIGURE 2

Figure 2. Result of variation partitioning analysis showing pure and shared effects of plant, soil, climatic and spatial variables on endophytic fungal richness. Numbers indicate proportions of explained variation (adjusted R² values). PCNM, spatial principal coordinates of neighbor matrices; PCoA, principal coordinates analysis of host phylogeny; N, soil total nitrogen; P, soil total phosphorus; MAT, mean annual temperature; MAP, mean annual precipitation.

Endophytic Fungal Community Composition

The endophytic fungal community was dominated by Leotiomycetes, Agaricomycetes, Eurotiomycetes, and Sordariomycetes (Figure 3A). Kruskal–Wallis tests showed that plant genus had a significant effect on the relative abundance of Leotiomycetes (χ² = 15.829, P = 0.001), Eurotiomycetes (χ² = 23.32, P < 0.001) and Sordariomycetes (χ² = 28.773, P < 0.001), but not Agaricomycetes (χ² = 3.428, P = 0.33) (Figures 3B–E). For example, the relative abundance of Leotiomycetes was significantly higher in Betula than in Alnus (Figure 3B). The relative abundances of Eurotiomycetes and Sordariomycetes were significantly lower in Alnus and Betula than in Carpinus and Corylus (Figures 3D,E).

FIGURE 3

Figure 3. Relative abundance of endophytic fungal classes in four plant genera. (A) Abundant classes (>1%) and others (unclassified and rare classes), (B) Leotiomycetes, (C) Agaricomycetes, (D) Eurotiomycetes, and (E) Sordariomycetes. Bars without shared letters indicate significant differences according to Dunn’s tests with Bonferroni adjustment at P < 0.05.

The PCoA ordination indicated that the endophytic fungal community composition was significantly different across the four plant genera (ANOSIM, R = 0.11, P = 0.001), with the first two principle coordinates explaining 7.1 and 6.8% of variation in fungal community, respectively (Figure 4). PerMANOVA revealed that 39.2% of variation in endophytic fungal community composition was explained, in which host phylogeny (PCoA1-16) explained the most (19.5%), followed by spatial vectors (PCNM1-10; 11.2%), soil (6.1%), climate (1.4%) and altitude (0.9%) (Table 2). Similarly, variation partitioning demonstrated that 28.9% of variation in endophytic fungal community composition was explained by host phylogeny (14.0%), spatial vectors (13.1%), soil (10.8%) and climate (8.1%), with pure effects of 6.7, 5.5, 4.4, and 1.7%, respectively (Figure 5). Taken together, environmental filtering coupled with dispersal limitation drives endophytic fungal community assembly, with a stronger effect of host phylogeny than that of other variables.

FIGURE 4

Figure 4. Principal coordinates analysis (PCoA) of endophytic fungal community composition. Percentage of variance explained by each PCoA axis is shown in parenthesis.

TABLE 2

Table 2. Relative importance of plant and abiotic variables on endophytic fungal community composition revealed by Permutational multivariate analysis of variance.

FIGURE 5

Figure 5. Result of variation partitioning analysis showing pure and shared effects of plant, soil, climatic and spatial variables on endophytic fungal community composition. Numbers indicate proportions of explained variation (adjusted R² values). PCNM, spatial principal coordinates of neighbor matrices; PCoA, principal coordinates analysis of host phylogeny; C, soil total carbon; N, soil total nitrogen; P, soil total phosphorus; Ca, soil total calcium; PSD, soil particle size distribution; MAT, mean annual temperature; MAP, mean annual precipitation.

Host/Fungus Preference and Checkerboard Score

Host/fungus preference analysis performed in Changbaishan site indicated that six of seven plant species showed significant preferences for abundant fungi (>0.5% of total sequences), 25 of 52 abundant fungal OTUs showed significant preferences for host plants, and 42 pairs of plants and abundant fungi exhibited remarkably strong preferences (Figure 6A). In Taiyueshan site, one of four plant species and 10 of 37 abundant fungi showed significant preferences toward partners, and 13 pairs of plants and abundant fungi indicated remarkably strong preferences (Figure 6B). In addition, the plant–fungus associations in both sites were significantly specialized at community level and there was less co-occurrence within endophytic fungal communities, as the observed values were all significantly higher than that generated by randomized null models (Figures 7A,B).

FIGURE 6

Figure 6. Preferences observed in plant and endophytic fungus associations in Changbaishan (A) and Taiyueshan (B). Standardized d′ estimates of preferences for fungal OTUs for indicated plant species (columns). Likewise, the standardized d′ estimate of preferences for plant species is indicated for each of the observed fungal OTUs (row). A cell in the matrix indicates a two-dimensional preference (2DP) estimate, indicating the extent an association of a focal plant–fungus pair was observed more/less frequently than expected by chance. The cell with asterisk inside represents significant preferences in plant–fungus pair. The relationship between 2DP and FDR-adjusted P-values shows that 2DP values larger than 2.3 (Changbaishan) and 1.8 (Taiyueshan) and those smaller than –2.6 (Changbaishan) and –2.5 (Taiyueshan) represented strong preferences and avoidance, respectively. Red line indicates that the 2DP value is not available as the standard deviation of the number of samples for the focal plant–fungal OTU pair across randomized matrices is zero. Because multiple species/OTUs were tested, the P values are shown as false discovery rates (FDRs) in the plant/fungus analysis. ^∗P < 0.05, ^∗∗P < 0.01, ^∗∗∗P < 0.001.

FIGURE 7

Figure 7. Community-level interaction specialization and checkerboard score of endophytic fungal community. (A) $H_2^{\prime }$ metric and (B) Checkerboard score reflecting co-occurrence within fungal community. The significant difference was estimated by comparing the observed value (red) with randomized values (yellow; mean ± SD) based on two-tailed t-test. ^∗P < 0.05, ^∗∗P < 0.01, ^∗∗∗P < 0.001.

Discussion

We found that the endophytic fungal OTU richness was significantly higher in Corylus than in other three plant genera. Furthermore, the endophytic fungal OTU richness was affected by host phylogeny, geographic distance, MAT and soil total N. Similar results were reported in previous studies on root-associated fungal communities, such as EM fungi (Tedersoo et al., 2013; Wu et al., 2018; Wang et al., 2019a), AM fungi (Davison et al., 2015; Van Geel et al., 2018) and pathogenic fungi (Wang et al., 2019b).

We found that host phylogeny significantly shaped endophytic fungal community composition. This finding suggests the important role of environmental filtering derived from host plants in determining fungal communities. All host phylogenetic variables (PCoA1-16) were retained as drivers for endophytic fungal community assembly in present study, suggesting that the splits among host plants at various taxonomic levels could affect fungal community (Tedersoo et al., 2013; Wu et al., 2018; Wang et al., 2019a). The host effect can be explained by phylogenetic niche conservatism, namely, closely related host species are more morphologically and functionally similar than distantly related species (Losos, 2008), and tend to share more similar fungal community. Indeed, in present study the distribution pattern of endophytic fungal community in the ordination corresponded well with phylogenetic position of the four plant genera in the phylogenetic tree. Namely, plant phylogenetic tree demonstrated that Alnus and Betula were relatively closely related genera, whilst Carpinus and Corylus were another two closely related genera (Supplementary Figure 2). Meanwhile, the PCoA ordination showed a relatively close cluster of endophytic fungal communities in Alnus and Betula, whilst Carpinus and Corylus exhibited another one (Figure 4). Additionally, the specificity between plants and fungi that derived from their long history of co-evolution has been widely documented to influence root-associated fungal communities including EM fungi (Tedersoo et al., 2010; Wu et al., 2018; Wang et al., 2019a), AM fungi (Sepp et al., 2019), endophytic fungi (Maciá-Vicente et al., 2008; Luo et al., 2017) and total fungi (Toju et al., 2016; Schroeder et al., 2019). Indeed, our plant/fungus preference analysis conducted in two sites showed that some abundant endophytic fungi and plant species significantly preferred to some certain partners. Moreover, community-level analyses in both sites also suggested that significantly specificity existed in plant–fungus associations in present study.

In addition to environmental filtering by host plants, abiotic environmental variables including soil (pH, nutrients, and particle size distribution) and climate (MAT and MAP) have also been observed to affect endophytic fungal community in this study. Similar findings were also reported in some studies on root-associated fungal communities, such as endophytic fungi (Maciá-Vicente et al., 2012; Glynou et al., 2016), EM fungi (Jarvis et al., 2013; Miyamoto et al., 2015; Wang et al., 2019a), AM fungi (Dumbrell et al., 2010; Zheng et al., 2014; Van Geel et al., 2018), pathogenic fungi (Wang et al., 2019b), and total fungi (Barnes et al., 2018; Li et al., 2018). Fungal taxa are known to prefer to different niche properties including nutrients and climate (Cox et al., 2010; Mcguire et al., 2010; Mundra et al., 2016; Barnes et al., 2018; Miyamoto et al., 2018). For example, some fungal taxa preferred to nutrient-rich habitats, while some taxa could grow well in nutrient-poor habitats (Cox et al., 2010; Mundra et al., 2016). In addition, some fungi only distributed in habitats with lower MAT and exhibited narrower temperature breadths, while some fungi were detected across wider temperature ranges and in warmer habitats (Miyamoto et al., 2018).

The endophytic fungal community composition was also affected by geographic distance, which is in accordance with the community assembly of other root-associated fungi, such as EM fungi (Tedersoo et al., 2011; Matsuoka et al., 2016; Wang et al., 2019a), AM fungi (Dumbrell et al., 2010; Wang et al., 2019b), pathogenetic fungi (Schroeder et al., 2019; Wang et al., 2019b), and total fungi (Wehner et al., 2014; Beck et al., 2015; Schroeder et al., 2019). All spatial vectors (PCNM1-10) were retained as predictors for endophytic fungal community composition in this study, reflecting fungal community assembly at various spatial variation ranging from relatively small to large scales (Tedersoo et al., 2011; Wu et al., 2018; Wang et al., 2019a). This pattern suggests that dispersal limitation is an important driver of endophytic fungal community turnover. The dispersal limitation effect could be attributable to the geographic isolation generating dispersal barriers that hamper the movement of fungal propagules (e.g., spore and hyphae) from the original to new habitats (Peay et al., 2010, 2012). In addition, fungal taxa per se commonly differ in their dispersal abilities, and the order of fungal propagules arriving new habitats could influence the outcome of fungal competition, and then cause fungal community variation (Kennedy and Bruns, 2005; Kennedy et al., 2009). Indeed, our checkerboard score analysis indicated that endophytic fungal community was competitively structured. This is in agreement with previous findings in root-associated fungal communities, such as EM fungi (Pickles et al., 2012; Wu et al., 2018; Wang et al., 2019a) and total fungi (Wehner et al., 2014; Toju et al., 2016).

Additionally, we found that host phylogeny showed a stronger effect on endophytic fungal community than other predictor variables tested in present study. This finding coincides with previous studies concerning root-associated fungi involved in a wide phylogenetic spectra of plants (17 to 51 species of two to 39 genera belonging to one to 22 families) (Tedersoo et al., 2013; Wehner et al., 2014; Wu et al., 2018; Schroeder et al., 2019; Wang et al., 2019a, b), yet conflicts with studies only considered a set of closely related plant species (fewer than seven species in one genus) (Tedersoo et al., 2009; Glassman et al., 2015, 2017; Erlandson et al., 2016). For example, host phylogeny explained more variation of fungal community than spatial and soil variables in studies involving in 25 species of 22 genera in Asteraceae (Wehner et al., 2014) and 45 plant species in 39 genera from 22 families (Wang et al., 2019b). By contrast, in studies those only considering less phylogenetically distant host plants (<7 species in one genus) (Tedersoo et al., 2009; Glassman et al., 2015, 2017; Erlandson et al., 2016), geographic distance and/or abiotic environmental variables showed stronger effects than host plants in determining fungal community composition. In our study, 22 plant species in four genera belonging to Betulaceae represented a relatively wide phylogenetic range and resulted in a stronger host effect relative to other variables on fungal community assembly. These studies suggest that host effect on fungal community assembly is phylogenetic scale-dependent.

We found that some endophytic fungal taxa obtained in present study could be classified as saprophytic and pathogenic groups in previous studies. For example, Cladophialophora, Mycena, Mortierella, and Meliniomyces could be assigned to saprophytic fungi, and Fusarium and Neonectria have been identified as plant pathogenic fungi (Rossman and Palm-Hernández, 2008; Tedersoo et al., 2014; Voříšková et al., 2014). Additionally, members of Phialocephala, as globally distributed root endophytic fungal taxa, have been proven to harbor a broad gene inventory that links pathogenic and saprophytic lifestyles (Schlegel et al., 2016). Although some endophytic fungal OTUs in our study showed latent saprophytic or pathogenic functions as reported in previous studies, the ectomycorrhizae selected for analysis in present study were all apparently healthy root tips without any disease. This finding suggests that endophytic fungi could include some potential saprophytic and pathogenic fungi, and their function need to be explored in future research. In conclusion, this study firstly characterizes endophytic fungal community assembly in the ectomycorrhizae of Betulaceae plants at a large scale (ca. 2.3 million km²). The endophytic fungal OTU richness was significantly related with host phylogeny, geographic distance, soil, and climate. The endophytic fungal community composition was significantly influenced by host phylogeny, geographic distance, soil and climatic variables, in which host phylogeny explained a larger variation than other predictor variables. This finding suggests that both environmental filtering (i.e., host plant phylogeny, soil and climate) and dispersal limitation (i.e., geographic distance) significantly drive community turnover of endophytic fungi, with host phylogeny being a better predictor than any other variables. Our study contributes to a growing literature showing the effect of environmental filtering and dispersal limitation on fungal community.

Data Availability Statement

The datasets generated for this study can be found in the European Nucleotide Archive, LR586058–LR588163.

Author Contributions

Y-LW and L-DG conceived and designed the study. Y-LW performed the experiments, analyzed the data, and drafted the manuscript. Y-LW, CG, LC, N-NJ, B-WW, P-PL, X-CL, and YZ performed the field sampling work. XQ, PM, and BB gave some suggestion for English improvement. L-DG revised and approved the final manuscript.

Funding

This work was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant Number XDB31000000), National Natural Science Foundation of China (Grant Numbers 91751113, 31210103910, and 31470545), and Beijing Natural Science Foundation (Grant Number 5182022).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We are grateful to Prof. Ye Chen (Jiujiang University), Prof. Xian-Kun Li (Guangxi Institute of Botany of Chinese Academy of Sciences), Prof. Shui-Sheng You (Fujian Agriculture and Forestry University), Dr. Ji Ye (Institute of Applied Ecology of Chinese Academy of Sciences), and Dr. Hong-Jin Dong (Kunming Institute of Botany of Chinese Academy of Sciences) for their help in identifying plant species. We also thank Dr. Hirokazu Toju (Center for Ecological Research of Kyoto University) for his help regarding the analysis of host/fungus preferences.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.03105/full#supplementary-material

References

Almario, J., Jeena, G., Wunder, J., Langen, G., Zuccaro, A., Coupland, G., et al. (2017). Root-associated fungal microbiota of nonmycorrhizal Arabis alpina and its contribution to plant phosphorus nutrition. Proc. Natl. Acad. Sci. U.S.A. 114, E9403–E9412. doi: 10.1073/pnas.1710455114

PubMed Abstract | CrossRef Full Text | Google Scholar

Altschul, S. F., Gish, W., Miller, W., Myers, E. W., and Lipman, D. J. (1990). Basic local alignment search tool. J. Mol. Biol. 215, 403–410. doi: 10.1016/S0022-2836(05)80360-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Aponte, C., García, L. V., Marañón, T., and Gardes, M. (2010). Indirect host effect on ectomycorrhizal fungi: leaf fall and litter quality explain changes in fungal communities on the roots of co-occurring Mediterranean oaks. Soil Biol. Biochem. 42, 788–796. doi: 10.1016/j.soilbio.2010.01.014

CrossRef Full Text | Google Scholar

Barnes, C. J., Van Der Gast, C. J., Mcnamara, N. P., Rowe, R., and Bending, G. D. (2018). Extreme rainfall affects assembly of the root-associated fungal community. New Phytol. 220, 1172–1184. doi: 10.1111/nph.14990

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartón, K. (2018). MuMIn: Multi-Model Inference. R Package Version 1. 42. 1. Available at: http://cran.r-project.org/package=MuMIn (accessed March 11, 2019).

Google Scholar

Beck, S., Powell, J. R., Drigo, B., Cairney, J. W. G., and Anderson, I. C. (2015). The role of stochasticity differs in the assembly of soil- and root-associated fungal communities. Soil Biol. Biochem. 80, 18–25. doi: 10.1016/j.soilbio.2014.09.010

CrossRef Full Text | Google Scholar

Ciccolini, V., Ercoli, L., Davison, J., Vasar, M., Opik, M., and Pellegrino, E. (2016). Land-use intensity and host plant simultaneously shape the composition of arbuscular mycorrhizal fungal communities in a Mediterranean drained peatland. FEMS Microbiol. Ecol. 92:fiw186. doi: 10.1093/femsec/fiw186

PubMed Abstract | CrossRef Full Text | Google Scholar

Clarke, K. R., Somerfield, P. J., and Chapman, M. G. (2006). On resemblance measures for ecological studies, including taxonomic dissimilarities and a zero-adjusted Bray-Curtis coefficient for denuded assemblages. J. Exp. Mar. Biol. Ecol. 330, 55–80. doi: 10.1016/j.jembe.2005.12.017

CrossRef Full Text | Google Scholar

Cottenie, K. (2005). Integrating environmental and spatial processes in ecological community dynamics. Ecol. Lett. 8, 1175–1182. doi: 10.1111/j.1461-0248.2005.00820.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cox, F., Barsoum, N., Lilleskov, E. A., and Bidartondo, M. I. (2010). Nitrogen availability is a primary determinant of conifer mycorrhizas across complex environmental gradients. Ecol. Lett. 13, 1103–1113. doi: 10.1111/j.1461-0248.2010.01494.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Davison, J., Moora, M., Opik, M., Adholeya, A., Ainsaar, L., Bâ, A., et al. (2015). Global assessment of arbuscular mycorrhizal fungus diversity reveals very low endemism. Science 349, 970–973. doi: 10.1126/science.aab1161

PubMed Abstract | CrossRef Full Text | Google Scholar

Dickie, I. A. (2007). Host preference, niches and fungal diversity. New Phytol. 174, 230–233. doi: 10.1111/j.1469-8137.2007.02055.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dray, S., Legendre, P., and Blanchet, G. (2009). Packfor: Forward Selection with Permutation. R package version 0.0-7. Available at: http://r-forge.rproject.org/R/?group_id=195 (accessed July 25, 2016)Google Scholar

Dray, S., Legendre, P., and Peres-Neto, P. R. (2006). Spatial modelling: a comprehensive framework for principal coordinate analysis of neighbour matrices (PCNM). Ecol. Model. 196, 483–493. doi: 10.1016/j.ecolmodel.2006.02.015

CrossRef Full Text | Google Scholar

Dumbrell, A. J., Nelson, M., Helgason, T., Dytham, C., and Fitter, A. H. (2010). Relative roles of niche and neutral processes in structuring a soil microbial community. ISME J. 4, 337–345. doi: 10.1038/ismej.2009.122

PubMed Abstract | CrossRef Full Text | Google Scholar

Erlandson, S. R., Savage, J. A., Cavender-Bares, J. M., and Peay, K. G. (2016). Soil moisture and chemistry influence diversity of ectomycorrhizal fungal communities associating with willow along an hydrologic gradient. FEMS Microbiol. Ecol. 92:fiv148. doi: 10.1093/femsec/fiv148

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, C., Shi, N. N., Liu, Y. X., Peay, K. G., Zheng, Y., Ding, Q., et al. (2013). Host plant genus−level diversity is the best predictor of ectomycorrhizal fungal diversity in a Chinese subtropical forest. Mol. Ecol. 22, 3403–3414. doi: 10.1111/mec.12297

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, C., Zhang, Y., Shi, N. N., Zheng, Y., Chen, L., Wubet, T., et al. (2015). Community assembly of ectomycorrhizal fungi along a subtropical secondary forest succession. New Phytol. 205, 771–785. doi: 10.1111/nph.13068

PubMed Abstract | CrossRef Full Text | Google Scholar

Glassman, S. I., Peay, K. G., Talbot, J. M., Smith, D. P., Chung, J. A., Taylor, J. W., et al. (2015). A continental view of pine-associated ectomycorrhizal fungal spore banks: a quiescent functional guild with a strong biogeographic pattern. New Phytol. 205, 1619–1631. doi: 10.1111/nph.13240

PubMed Abstract | CrossRef Full Text | Google Scholar

Glassman, S. I., Wang, I. J., and Bruns, T. D. (2017). Environmental filtering by pH and soil nutrients drives community assembly in fungi at fine spatial scales. Mol. Ecol. 26, 6960–6973. doi: 10.1111/mec.14414

PubMed Abstract | CrossRef Full Text | Google Scholar

Glynou, K., Ali, T., Buch, A. K., Haghi Kia, S., Ploch, S., Xia, X., et al. (2016). The local environment determines the assembly of root endophytic fungi at a continental scale. Environ. Microbiol. 18, 2418–2434. doi: 10.1111/1462-2920.13112

PubMed Abstract | CrossRef Full Text | Google Scholar

Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G., and Jarvis, A. (2005). Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25, 1965–1978. doi: 10.1002/joc.1276

CrossRef Full Text | Google Scholar

Hooper, D. U., Bignell, D. E., Brown, V. K., Brussaard, L., Dangerfield, J. M., Wall, D. H., et al. (2000). Interactions between aboveground and belowground biodiversity in terrestrial ecosystems: patterns, mechanisms, and feedbacks. Bioscience 50, 1049–1061.

Google Scholar

Hubbell, S. P. (2001). The Unified Neutral Theory of Biodiversity and Biogeography. Princeton: Princeton University Press.

Google Scholar

Ishida, T. A., Nara, K., and Hogetsu, T. (2007). Host effects on ectomycorrhizal fungal communities: insight from eight host species in mixed conifer-broadleaf forests. New Phytol. 174, 430–440. doi: 10.1111/j.1469-8137.2007.02016.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Jarvis, S., Woodward, S., Alexander, I. J., and Taylor, A. F. S. (2013). Regional scale gradients of climate and nitrogen deposition drive variation in ectomycorrhizal fungal communities associated with native Scots pine. Glob. Change Biol. 19, 1688–1696. doi: 10.1111/gcb.12178

PubMed Abstract | CrossRef Full Text | Google Scholar

Katoh, K., and Standley, D. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780. doi: 10.1093/molbev/mst010

PubMed Abstract | CrossRef Full Text | Google Scholar

Kennedy, P. G., and Bruns, T. D. (2005). Priority effects determine the outcome of ectomycorrhizal competition between two Rhizopogon species colonizing Pinus muricata seedlings. New Phytol. 166, 631–638. doi: 10.1111/j.1469-8137.2005.01355.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kennedy, P. G., Peay, K. G., and Bruns, T. D. (2009). Root tip competition among ectomycorrhizal fungi: are priority effects a rule or an exception? Ecology 90, 2098–2107. doi: 10.1890/08-1291.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kivlin, S. N., Winston, G. C., Goulden, M. L., and Treseder, K. K. (2014). Environmental filtering affects soil fungal community composition more than dispersal limitation at regional scales. Fungal Ecol. 12, 14–25. doi: 10.1016/j.funeco.2014.04.004

CrossRef Full Text | Google Scholar

Kõljalg, U., Nilsson, R. H., Abarenkov, K., Tedersoo, L., Taylor, A. F., Bahram, M., et al. (2013). Towards a unified paradigm for sequence−based identification of fungi. Mol. Ecol. 22, 5271–5277. doi: 10.1111/mec.12481

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Stecher, G., and Tamura, K. (2016). MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054

PubMed Abstract | CrossRef Full Text | Google Scholar

Landeweert, R., Hoffland, E., Finlay, R. D., Kuyper, T. W., and van Breemen, N. (2001). Linking plants to rocks: ectomycorrhizal fungi mobilize nutrients from minerals. Trends Ecol. Evol. 16, 248–254. doi: 10.1016/S0169-5347(01)02122-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, P. Q., and Skvortsov, A. K. (1999). Betulaceae. Flora of China Cycadaceae through Fagaceae, Vol. 4. Beijing: Science Press, 286–313.

Google Scholar

Li, S., Shakoor, A., Wubet, T., Zhang, N., Liang, Y., and Ma, K. (2018). Fine-scale variations of fungal community in a heterogeneous grassland in Inner Mongolia: effects of the plant community and edaphic parameters. Soil Biol. Biochem. 122, 104–110. doi: 10.1016/j.soilbio.2018.04.007

CrossRef Full Text | Google Scholar

Liang, M., Liu, X., Gilbert, G. S., Zheng, Y., Luo, S., Huang, F., et al. (2016). Adult trees cause density−dependent mortality in conspecific seedlings by regulating the frequency of pathogenic soil fungi. Ecol. Lett. 19, 1448–1456. doi: 10.1111/ele.12694

PubMed Abstract | CrossRef Full Text | Google Scholar

Losos, J. B. (2008). Phylogenetic niche conservatism, phylogenetic signal and the relationship between phylogenetic relatedness and ecological similarity among species. Ecol. Lett. 11, 995–1003. doi: 10.1111/j.1461-0248.2008.01229.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, J., Walsh, E., Miller, S., Blystone, D., Dighton, J., and Ning, Z. (2017). Root endophytic fungal communities associated with pitch pine, switchgrass, and rosette grass in the pine barrens ecosystem. Fungal Bio. 121, 478–487. doi: 10.1016/j.funbio.2017.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Maciá-Vicente, J. G., Ferraro, V., Burruano, S., and Lopez-Llorca, L. V. (2012). Fungal assemblages associated with roots of halophytic and non-halophytic plant species vary differentially along a salinity gradient. Microb. Ecol. 64, 668–679. doi: 10.1007/s00248-012-0066-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Maciá-Vicente, J. G., Jansson, H. B., and Abdullah, S. K. (2008). Fungal root endophytes from natural vegetation in Mediterranean environments with special reference to Fusarium spp. FEMS Microbiol. Ecol. 64, 90–105. doi: 10.1111/j.1574-6941.2007.00443.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Martinez-Garcia, L. B., Richardson, S. J., Tylianakis, J. M., Peltzer, D. A., and Dickie, I. A. (2015). Host identity is a dominant driver of mycorrhizal fungal community composition during ecosystem development. New Phytol. 205, 1565–1576. doi: 10.1111/nph.13226

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsuoka, S., Mori, A. S., Kawaguchi, E., Hobara, S., and Osono, T. (2016). Disentangling the relative importance of host tree community, abiotic environment and spatial factors on ectomycorrhizal fungal assemblages along an elevation gradient. FEMS Microbiol. Ecol. 92:fiw044. doi: 10.1093/femsec/fiw044

PubMed Abstract | CrossRef Full Text | Google Scholar

Mcguire, K. L., Bent, E., Borneman, J., Majumder, A., Allison, S. D., and Treseder, K. K. (2010). Functional diversity in resource use by fungi. Ecology 91, 2324–2332. doi: 10.1890/09-0654.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Miyamoto, Y., Sakai, A., Hattori, M., and Nara, K. (2015). Strong effect of climate on ectomycorrhizal fungal composition: evidence from range overlap between two mountains. ISME J. 9, 1870–1879. doi: 10.1038/ismej.2015.8

PubMed Abstract | CrossRef Full Text | Google Scholar

Miyamoto, Y., Terashima, Y., and Nara, K. (2018). Temperature niche position and breadth of ectomycorrhizal fungi: reduced diversity under warming predicted by a nested community structure. Glob. Change Biol. 24, 5724–5737. doi: 10.1111/gcb.14446

PubMed Abstract | CrossRef Full Text | Google Scholar

Mundra, S., Bahram, M., and Eidesen, P. B. (2016). Alpine bistort (Bistorta vivipara) in edge habitat associates with fewer but distinct ectomycorrhizal fungal species: a comparative study of three contrasting soil environments in Svalbard. Mycorrhiza 26, 809–818. doi: 10.1007/s00572-016-0716-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Nara, K. (2009). Spores of ectomycorrhizal fungi: ecological strategies for germination and dormancy. New Phytol. 181, 245–248. doi: 10.1111/j.1469-8137.2008.02691.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Oksanen, J., Kindt, R., Legendre, P., O’hara, B., Stevens, M. H. H., Oksanen, M. J., et al. (2013). Vegan: Community Ecology Package. R Package Version 2.0-10. Available at: http://CRAN.R-project.org/package=vegan(accessed September 21, 2013)

Google Scholar

Paradis, E., Claude, J., and Strimmer, K. (2004). APE: analyses of phylogenetics and evolution in R language. Bioinformatics 20, 289–290. doi: 10.1093/bioinformatics/btg412

PubMed Abstract | CrossRef Full Text | Google Scholar

Peay, K. G., Garbelotto, M., and Bruns, T. D. (2010). Evidence of dispersal limitation in soil microorganisms: isolation reduces species richness on mycorrhizal tree islands. Ecology 91, 3631–3640. doi: 10.1890/09-2237.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Peay, K. G., Kennedy, P. G., and Bruns, T. D. (2008). Fungal community ecology: a hybrid beast with a molecular master. Bioscience 58, 799–810. doi: 10.1641/B580907

CrossRef Full Text | Google Scholar

Peay, K. G., Schubert, M. G., Nguyen, N. H., and Bruns, T. D. (2012). Measuring ectomycorrhizal fungal dispersal: macroecological patterns driven by microscopic propagules. Mol. Ecol. 21, 4122–4136. doi: 10.1111/j.1365-294X.2012.05666.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Petrini, O. (1991). “Fungal endophytes of tree leaves,” in Microbial Ecology of Leaves, eds J. H. Andrews and S. S. Hirano, (Berlin: Springer Verlag), 179–197. doi: 10.1007/978-1-4612-3168-4_9

CrossRef Full Text | Google Scholar

Pickles, B. J., Genney, D. R., Anderson, I. C., and Alexander, I. J. (2012). Spatial analysis of ectomycorrhizal fungi reveals that root tip communities are structured by competitive interactions. Mol. Ecol. 21, 5110–5123. doi: 10.1111/j.1365-294X.2012.05739.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pohlert, T. (2014). The Pairwise Multiple Comparison of Mean Ranks Package (PMCMR). R Package Version 4.2. 2016. Available at: http://CRAN.R-project.org/package=PMCMR (accessed January 25, 2017)

Google Scholar

Põlme, S., Bahram, M., Yamanaka, T., Nara, K., Dai, Y. C., Grebenc, 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. doi: 10.1111/nph.12170

PubMed Abstract | CrossRef Full Text | Google Scholar

R Development Core Team, (2018). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.

Google Scholar

Rossman, A. Y., and Palm-Hernández, M. E. (2008). Systematics of plant pathogenic fungi: why it matters. Plant Dis. 92, 1376–1386. doi: 10.1094/pdis-92-10-1376

PubMed Abstract | CrossRef Full Text | Google Scholar

Schlegel, M., Münsterkötter, M., Güldener, U., Bruggmann, R., Duò, A., Hainaut, M., et al. (2016). Globally distributed root endophyte Phialocephala subalpina links pathogenic and saprophytic lifestyles. BMC Genomics 17:1015. doi: 10.1186/s12864-016-3369-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Schroeder, J. W., Martin, J. T., Angulo, D. F., Razo, I. A. D., Barbosa, J. M., Perea, R., et al. (2019). Host plant phylogeny and abundance predict root-associated fungal community composition and diversity of mutualists and pathogens. J. Ecol. 107, 1557–1566. doi: 10.1111/1365-2745.13166

CrossRef Full Text | Google Scholar

Sepp, S. K., Davison, J., Jairus, T., Vasar, M., Moora, M., Zobel, M., et al. (2019). Non-random association patterns in a plant-mycorrhizal fungal network reveal host-symbiont specificity. Mol. Ecol. 28, 365–378. doi: 10.1111/mec.14924

PubMed Abstract | CrossRef Full Text | Google Scholar

Tedersoo, L., Bahram, M., Jairus, T., Bechem, E., Chinoya, S., Mpumba, R., et al. (2011). Spatial structure and the effects of host and soil environments on communities of ectomycorrhizal fungi in wooded savannas and rain forests of Continental Africa and Madagascar. Mol. Ecol. 20:3071–3080. doi: 10.1111/j.1365-294X.2011.05145.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Tedersoo, L., Bahram, M., Põlme, S., Kõljalg, U., Yorou, N. S., Wijesundera, R., et al. (2014). Global diversity and geography of soil fungi. Science 346:1256688. doi: 10.1126/science.1256688

PubMed Abstract | CrossRef Full Text | Google Scholar

Tedersoo, L., Bahram, M., Toots, M., Diedhiou, A. G., Henkel, T. W., Kjoller, R., et al. (2012). Towards global patterns in the diversity and community structure of ectomycorrhizal fungi. Mol. Ecol. 21, 4160–4170. doi: 10.1111/j.1365-294X.2012.05602.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Tedersoo, L., Mett, M., Ishida, T. A., and Bahram, M. (2013). Phylogenetic relationships among host plants explain differences in fungal species richness and community composition in ectomycorrhizal symbiosis. New Phytol. 199, 822–831. doi: 10.1111/nph.12328

PubMed Abstract | CrossRef Full Text | Google Scholar

Tedersoo, L., Sadam, A., Zambrano, M., Valencia, R., and Bahram, M. (2010). Low diversity and high host preference of ectomycorrhizal fungi in Western Amazonia, a neotropical biodiversity hotspot. ISME J. 4, 465–471. doi: 10.1038/ismej.2009.131

PubMed Abstract | CrossRef Full Text | Google Scholar

Tedersoo, L., Suvi, T., Jairus, T., Ostonen, I., and Põlme, S. (2009). Revisiting ectomycorrhizal fungi of the genus Alnus: differential host specificity, diversity and determinants of the fungal community. New Phytol. 182, 727–735. doi: 10.1111/j.1469-8137.2009.02792.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Toju, H., Tanabe, A., and Ishii, H. (2016). Ericaceous plant–fungus network in a harsh alpine–subalpine environment. Mol. Ecol. 25, 3242–3257. doi: 10.1111/mec.13680

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Dam, N. M., and Heil, M. (2011). Multitrophic interactions below and above ground: en route to the next level. J. Ecol. 99, 77–88. doi: 10.1111/j.1365-2745.2010.01761.x

CrossRef Full Text | Google Scholar

van der Heijden, M. G., 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

PubMed Abstract | CrossRef Full Text | Google Scholar

van der Linde, S., Suz, L. M., Orme, C. D. L., Cox, F., Andreae, H., Asi, E., et al. (2018). Environment and host as large-scale controls of ectomycorrhizal fungi. Nature 558, 243–258. doi: 10.1038/s41586-018-0189-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Geel, M., Jacquemyn, H., Plue, J., Saar, L., Kasari, L., Peeters, G., et al. (2018). Abiotic rather than biotic filtering shapes the arbuscular mycorrhizal fungal communities of European seminatural grasslands. New Phytol. 220, 1262–1272. doi: 10.1111/nph.14947

PubMed Abstract | CrossRef Full Text | Google Scholar

Voříšková, J., Brabcová, V., Cajthaml, T., and Baldrian, P. (2014). Seasonal dynamics of fungal communities in a temperate oak forest soil. New Phytol. 201, 269–278. doi: 10.1111/nph.12481

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y. L., Gao, C., Chen, L., Ji, N. N., Wu, B. W., Li, X. C., et al. (2019a). Host plant phylogeny and geographic distance strongly structure Betulaceae-associated ectomycorrhizal fungal communities in Chinese secondary forest ecosystems. FEMS Microbiol. Ecol. 95:fiz037. doi: 10.1093/femsec/fiz037

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z. H., Jiang, Y., Deane, D. C., He, F. L., Shu, W. S., and Liu, Y. (2019b). Effects of host phylogeny, habitat and spatial proximity on host−specificity and diversity of pathogenic and mycorrhizal fungi in a subtropical forest. New Phytol. 223, 462–474. doi: 10.1111/nph.15786

PubMed Abstract | CrossRef Full Text | Google Scholar

Wardle, D. A. (2006). The influence of biotic interactions on soil biodiversity. Ecol. Lett. 9, 870–886. doi: 10.1111/j.1461-0248.2006.00931.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wardle, D. A., Bardgett, R. D., Klironomos, J. N., Setala, H., Van Der Putten, W. H., and Wall, D. H. (2004). Ecological linkages between aboveground and belowground biota. Science 304, 1629–1633. doi: 10.1126/science.1094875

PubMed Abstract | CrossRef Full Text | Google Scholar

Wehner, J., Powell, J. R., Muller, L. A. H., Caruso, T., Veresoglou, S. D., Hempel, S., et al. (2014). Determinants of root-associated fungal communities within Asteraceae in a semi-arid grassland. J. Ecol. 102, 425–436. doi: 10.1111/1365-2745.12197

CrossRef Full Text | Google Scholar

Wu, B. W., Gao, C., Chen, L., Buscot, F., Goldmann, K., Purahong, W., et al. (2018). Host phylogeny is a major determinant of Fagaceae-associated ectomycorrhizal fungal community assembly at a regional scale. Front. Microbiol. 9:2409. doi: 10.3389/fmicb.2018.02409

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, Y., Chen, L., Luo, C. Y., Zhang, Z. H., Wang, S. P., and Guo, L. D. (2016). Plant identity exerts stronger effect than fertilization on soil arbuscular mycorrhizal fungi in a sown pasture. Microb. Ecol. 72, 647–658. doi: 10.1007/s00248-016-0817-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, Y., Kim, Y. C., Tian, X. F., Chen, L., Yang, W., Gao, C., et al. (2014). Differential responses of arbuscular mycorrhizal fungi to nitrogen addition in a near pristine Tibetan alpine meadow. FEMS Microbiol. Ecol. 89, 594–605. doi: 10.1111/1574-6941.12361

PubMed Abstract | CrossRef Full Text | Google Scholar