Wen Zhao
Kaichuan Huang
Reyila Mumin
Yifei Sun
Baokai Cui- State Key Laboratory of Efficient Production of Forest Resources, School of Ecology and Nature Conservation, Beijing Forestry University, Beijing, China
Soil fungi play a critical role in the biogeochemical cycles of forest ecosystems. Larix gmelinii is a strong and important timber tree species, which forms close associations with a wide range of soil fungi. However, the temporal-spatial disparity effects on the assembly of soil fungal communities in L. gmelinii forests are poorly understood. To address these questions, a total of 120 samples, including 60 bulk soil and 60 root samples, were collected from Aershan and Genhe in July (summer) and October (autumn)2021. We obtained 7,788 operational taxonomic units (OTUs) after merging, filtering, and rarefying using high-throughput sequencing. The dominant phyla are Basidiomycota, Ascomycota, Mortierellomycota, and Mucoromycota. There were 13 dominant families, among which the families with average relative abundance more than 5% included Thelephoraceae, Mortierellaceae, Archaeorhizomycoaceae, and Inocybaceae. In the functional guilds, symbiotrophic fungi had a relative advantage in the identified functions, and the relative abundances of pathotrophic and saprotrophic fungi varied significantly between sites. There were 12 families differentially expressed across compartments, 10 families differentially expressed between seasons, and 69 families were differentially expressed between sites. The variation in alpha diversity in the bulk soil was greater than that in the rhizosphere soil. Among the three parts (compartment, season, and site), the site had a crucial effect on the beta diversity of the fungal community. Deterministic processes dominated fungal community assembly in Genhe, whereas stochastic processes dominated in Aershan. Soil physicochemical properties and climatic factors significantly affected fungal community structure, among which soil total nitrogen and pH had the greatest effect. This study highlights that spatial variations play a vital role in the structure and assembly of soil fungal communities in L. gmelinii forests, which is of great significance for us in maintaining the health of the forests.
1 Introduction
As an important component of soil microorganisms, soil fungi play a vital role in the cycling of carbon, nitrogen, and phosphorus in the ecosystem (Lladó et al., 2018) and promote the biogeochemical cycle of the earth (Crowther et al., 2019; Jiang et al., 2021). Among these soil fungi, mycorrhizal fungi established mutualistic interaction with their hosts (Zhong et al., 2021), saprophytic fungi carried out material and energy cycling in soil (Setälä and McLean, 2004), and pathogenic fungi were partly transferring in soil and causing disease (Maron et al., 2011). The closest connection between soil fungi and plants occurs in the rhizosphere, which is due to the root exudates of the plant and nutrient availability of the soil (Li et al., 2016; Praeg and Illmer, 2020). At the same time, plants also choose relevant fungi, so vegetation types have a significant impact on the fungal community in rhizosphere soil (Berendsen et al., 2012; Liu et al., 2023).
Larix gmelinii is a deciduous tree of the Pinaceae family with rich timber accumulation and extreme cold tolerance (Kajimoto et al., 1999; Sun et al., 2007). Northeast China is mainly located in the Greater Khingan Mountains, which are important forest resources. Larix gmelinii is an ectomycorrhizal tree species (Ouimette et al., 2013; Miyamoto et al., 2021; Meng et al., 2023). Numerous fungi attach to their host for nutrient uptake (Ishida et al., 2007; Orwin et al., 2011). It can also help plants absorb nutrients (Kuang et al., 2021) and maintain health (Berendsen et al., 2012). Complex plant-associated fungal communities are directly influenced by organic compounds secreted from plant roots (Broeckling et al., 2008).
Fungal communities in the forest soil exhibited strong temporal and spatial differences. Previous studies showed that ectomycorrhizal fungi play a significant role in plant growth season, whereas saprophytic fungi evidently increase during the cold season (Voříšková et al., 2014; Santalahti et al., 2016). Seasonal changes in temperature and precipitation are also responsible for variations in soil fungal communities (Tedersoo et al., 2014; Zhang et al., 2018). In our previous study, temperature and precipitation were found to be the most important factors affecting fungal diversity with seasonal variations (Zhao et al., 2023). As for spatial differences, soil fungal communities have different structures between geographic sites, and their physicochemical properties directly constitute different fungal living environments, such as nitrogen (Zhou et al., 2016), pH (Zhang et al., 2016), and moisture (Zumsteg et al., 2013). Our previous study showed that differences in temperature and precipitation are also the main driving factors for differences in fungal communities at large regional scales (Wang et al., 2022). Even at a small local scale, there are significant differences in soil fungal communities among different vegetation types (Huang et al., 2023).
Community assembly, which is the balance between fungal community activity and environmental filtration, is categorized into deterministic and stochastic processes. Deterministic processes involve biotic and abiotic factors, including interspecific interactions such as competition, predation, symbiosis, and environmental filtering. Stochastic processes are understood as the drift, migration, speciation, and species extinction of microorganisms (Nemergut et al., 2013). When the environment limited the life of fungi, deterministic processes were the main factors affecting communities. Therefore, the communities dominated by deterministic processes were more unstable. In contrast, the communities dominated by stochastic processes were resistant to environmental change (Powell et al., 2015). It is difficult to understand the key factors that affect soil fungal community assembly processes, because we cannot directly measure an indicator to gauge the process. Ecologists have obtained a deeper and more scientific understanding of the community assembly process and shown that the stochastic and deterministic processes can be strongly influenced by selection, dispersal, diversification, and drift (Vellend, 2010; Stegen et al., 2012), and compute a value through a framework to delineate the processes (Stegen et al., 2013). Zhou and Ning (2017) suggested that one of the main goals of microbial ecology is to project future scenarios of microbial community structure and functions in a changing environment.
In this study, we aimed to determine the differences in rhizosphere and bulk soil fungal communities with different temporal and spatial variations and to determine the relationship between the fungal community and environmental variables. Considering the long winters and frozen soil in Frigid-Temperate Zones, and combining with the annual variation in temperature and precipitation, we finally chose July and October to represent the seasons (summer and autumn). We selected two representative sites in the main distribution of L. gmelinii in the Greater Khingan Mountains and collected 120 rhizosphere and bulk soil samples for analysis and exploration in July and October of 2021. Furthermore, we identified the main taxa and their functional guilds, summarized their variations, and calculated the community assembly processes and the main factors that caused these results. This was a multifaceted systematic analysis to explore the temporal-spatial patterns and driving factors of the soil fungal community in L. gmelinii forest.
2 Materials and methods
2.1 Study sites
The present study was conducted in two typical areas of L. gmelinii forests in the Greater Khingan Mountains, China (Figure 1, Table 1). One study area was located in the Aershan National Forest Park, Aershan (47°15′17″N, 120°17′34″E, shortened to Aershan), the southernmost area of L. gmelinii growing in China. The other was located in the Jinlin Forest Farm, Genhe (51°21′51″N, 121°20′31″E, shortened to Genhe), which is the core area for the distribution of L. gmelinii. The study area is located in a high-latitude permafrost distribution area in China (Ding et al., 2019). The winters in the study areas were long and cold, the summers were short and rainy, and spring and autumn were cool with less precipitation (Liu et al., 2020). The dominant tree species in these areas were L. gmelinii and Betula platyphylla, with a forest coverage rate of 80%. The zonal soil at both sites is Retisol (IUSS Working Group, 2015). The soil in Aershan is a typical black soil that is soft, moist, and nutritious, whereas the soil in Genhe is rich in gravel. A shovel was used to turn over the stones and obtain the soil in the stone cracks.
Figure 1 Distribution of the two forests.
Table 1 Localities, climates and forest conditions of the sampling sites.
2.2 Collection of bulk and rhizosphere soil
The rhizosphere and bulk soil samples were collected in July and October 2021, representing summer and autumn, respectively. At each sampling site, we established three 20 m × 20 m plots at least 100 m apart and collected five soil samples from five points (four corners and the center) as independent samples. Roots were extracted by digging the trees near the points to collect rhizosphere soil. In total, 120 samples (5 independent replicates× 3 plots × 2 compartments × 2 sites × 2 seasons) were used for DNA extraction. Additionally, 60 bulk soil samples were used to determine the soil physicochemical properties. Specifically, bulk soil from 0 to –20 cm was obtained, sifted through a 2 mm mesh sieve, mixed, and then packed into centrifuge tubes and sealed bags. The roots were stored in sealed bags. The samples in centrifuge tubes were stored with dry ice during the sampling process and transported to the laboratory stored at −80°C. The roots were stored with ice in the field and then transported to the laboratory stored at 4°C for the rhizosphere soil. The bulk soil was stored in sealed bags and transported at room temperature.
Rhizosphere soil was then separated from roots by shaking root and soil in 25 mL 0.1 M sterile phosphate buffer (7.1 g Na2HPO4 and 4.4 g NaH2PO4·H2O were added to 820 mL deionized water; pH 7.0) in a 50 mL plastic conical centrifuge tube at 67 g for 5 min. The soil suspension was centrifuged at 4000 g for 15 min, and the obtained soil pellet was regarded as the rhizosphere soil and stored at −80°C until DNA extraction (Liu et al., 2021).
2.3 Environmental variables
To determine soil properties, we followed the method described by Bao (2000). The soil in the sealed bags was dried naturally to determine the following physicochemical properties: soil organic carbon (SOC) was determined by the potassium dichromate volumetric method (Li et al., 2020), available phosphorus (AP) was determined by molybdenum antimony resistance colorimetry (Cai et al., 2018), pH was determined by the potential method (Li et al., 2020), total nitrogen (TN) was determined by Kjeldahl method (Shi et al., 2014), and the cation exchange capacity (CEC) was determined by the ammonium acetate exchange method (Liang et al., 2015).
Historical monthly weather data were obtained from WorldClim (https://www.worldclim.org/data/monthlywth.html#; accessed July 20, 2023), including average minimum temperature (°C), average maximum temperature (°C), and total precipitation (mm) in July and October 2021. The mean annual temperature and precipitation for 1970–2000 were obtained from WorldClim.
2.4 Soil DNA extraction and PCR amplification
DNA was extracted using the DNeasy PowerSoil Pro Kit (Qiagen, Frankfurt, Germany). A NanoDrop NC2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) was used to quantify the DNA, and the quality of the DNA was detected using 1.2% agarose gel electrophoresis. The ITS1 region of fungi was amplified, and the primers were ITS5F (5-GGAAGTAAAGTAAAAGTCGTAAAGG-3) and ITS2R (5-GCTGCGTTCTTCATCGATGC-3) (White et al., 1990). PCR system (20 μL): 2 μL (2.5 mM) of dNTP. Of note, 1 μL (10 µM) forward primers and 1 μL (10 µM) reverse primers, 2 μL template DNA, 10 μL ddH2O, and 4 μL Fast pfu DNA polymerases. Circulating system: 95°C for 2 min; 30 cycles at 95°C for 30 seconds, 55°C for 30 seconds, and 72°C for 30 seconds; 72°C for 5 min. PCR amplification was performed using an Applied Biosystems 2720 Thermal Cycler (Thermo Fisher Scientific, Waltham, MA, USA). The Quant-iT PicoGreen dsDNA assay was performed using a Microplate Reader FLx800 (BioTek, Burlington, Vermont, USA). PCR amplicons were purified and recovered by adding Vazyme VAHTSTM DNA Clean Beads (Vazyme, Nanjing, China), quantified with the fluorescent reagent of Quant-iT PicoGreen dsDNA assay kit (Invitrogen, Carlsbad, CA, USA), and then mixed in proportion to the sequencing amount. Pair-end 2 × 250 bp sequencing was performed using the Illumina MiSeq platform with a NovaSeq 6000 SP reagent kit at Shanghai Personal Biotechnology Co., Ltd. (Shanghai, China). All raw sequencing data from this study have been deposited in the NCBI database with the Short Read Archive (SRA) accession number PRJNA1037725.
2.5 Sequence analysis
The early sequence processing in this experiment was based on QIIME 2–2021.2 (Bolyen et al., 2019). A Demux plug-in was used to split the samples. The DADA2 plug-in was used to perform quality control, such as filtering and noise removal; the positions at which the forward and reverse read sequences should be truncated due to a decrease in quality were 222 and 232, respectively, and the minimum and maximum sequence lengths were 230 and 438, respectively. The Vsearch plug-in was used to cluster the sequences into operational taxonomic units (OTUs) according to the 97% similarity principle. The Phylogeny plug-in was used to generate the phylogenetic trees. According to the fungus UNITE v.8.2 database (Nilsson et al., 2019), we used the Feature Classifier plug-in to annotate the species. Functional guilds and trophic modes were obtained from the FUNGuild database (Nguyen et al., 2016).
ITS is not a conserved region, making the constructed phylogenetic trees unstable. In this case, a hybrid tree based on the Silva 138 18S database and Unite 8.2 fungal dynamic database was constructed using q2-ghost-tree (Fouquier et al., 2016). The results of species composition and diversity analysis based on the OTU selected according to the ghost tree were basically consistent with those of the ITS tree, and the community assembly processes of the two phylogenetic trees were consistent. However, compared to phylogenetic trees based on ITS primers, the minimum frequency of samples decreased from 39,602 to 9,499, resulting in the loss of 76% of OTUs. Therefore, we decided to construct the ITS phylogenetic tree in this study.
2.6 Statistical analysis
The data processing part was mainly completed with R 4.1.2 (R Core Team, 2022) and SPSS 26.0 (SPSS Inc., Chicago, IL, USA). The vegan package provides tools for describing the community ecology. First, we used the rarefy command of the vegan package to standardize the OTU table according to the minimum sequence number of samples and constructed the rarefaction curve using the rare-curve command. Counts of different taxonomic categories were obtained using the Phyloseq package. Then we grouped taxa and functions which was not identified or whose average relative abundance was less than 1% into “others,” and defined dominant taxa as which a relative abundance more than 1%. The MicrobiotaProcess package was used to obtain alpha diversity indices. For alpha diversity analysis, the Shannon diversity index, Chao1 richness index, and Pielou evenness index were selected for one-way analysis of variance (ANOVA) and Tukey’s honest significant difference (HSD) test at P < 0.05 to compare species diversity. Beta diversity was tested by Principal coordinate analysis (PCoA) based on the OTU level, and the differences between fungal communities were tested by Permutational multivariate analysis of variance (PERMANOVA) of 999 permutations. Differential expression analysis was performed using the DESeq2 package based on the family level except for unassigned families. Comparisons were performed in four groups for each part to avoid the influence of other parts, and differentially expressed families (P < 0.05) were categorized as up regulated (log2FoldChange > 2) or down regulated (log2FoldChange < −2). A family that was up regulated or down regulated in two or more groups in each part was defined as differentially expressed in this part. We selected the null model described by Stegen et al. (2013) to classify community assembly processes. We calculated the β-nearest taxon index (βNTI) and modified Raup-Crick index (RCBray) using iCAMP package. βNTI values < −2 (homogeneous selection) or > +2 (heterogeneous selection) indicated deterministic processes; −2 < βNTI values <+2 indicated stochastic processes including homogenizing dispersal (RCBray < −0.95), dispersal limitation (RCBray > 0.95) and undominated (−0.95 < RCBray < 0.95). The Mantel test was performed on community structure and environmental variables. A heatmap was used to show Spearman’s correlation between the relative abundance of dominant phyla and environmental variables. The data were visualized using the ggplot2 package. All analyses of variance (ANOVA) were calculated using SPSS, with P < 0.05 considered statistically significant.
3 Results
3.1 Soil properties and climatic factors
We only analyzed the soil properties in bulk soil, as we can only extract a small amount of rhizosphere soil for DNA extraction. There were significant differences in the physicochemical properties of the soil among the different variables (Table 2). Soil total nitrogen (TN), SOC, pH, and CEC in Aershan were significantly higher than those in Genhe. Available phosphorus (AP) in Genhe was significantly higher than that in Aershan. Between the two seasons, there were no significant differences in TN, SOC, and CEC in Aershan, and no significant difference in pH in Genhe. In summary, compared with the season, the differences in soil physicochemical properties between the two sites were greater.
Table 2 Soil properties compared by ANOVA (mean values ± S.E.s) and monthly average climate factors in 2021.
3.2 Fungal community composition and functional guilds
In total, 7,593,884 quality-filtered sequences were obtained from the 120 samples. After rarefying 39,602 sequences per sample (Figure 2A), 4,752,240 unique sequences were clustered into 7,788 operational taxonomic units (OTUs) with 97% identity. These OTUs belonged to 17 phyla, 60 classes, 138 orders, 362 families, and 864 genera.
Figure 2 Fungal community structure in different compartments, seasons, and sites. (A) The rarefaction curve of the total 120 samples. Fungal community composition at phylum (B) and family (C) levels, and the functional guilds (D) of the community. In Figure (D), red represents pathotroph, a series of blue represents saprotroph, a series of green represents symbiotroph. A, Aershan; G, Genhe; J, July; O: October; B, Bulk; R, Rhizosphere. BAJ, bulk soil in Aershan, July; BAO, bulk soil in Aershan, October; BGJ, bulk soil in Genhe, July; BGO, bulk soil in Genhe, October; RAJ, rhizosphere soil in Aershan, July; RAO, rhizosphere soil in Aershan, October; RGJ, rhizosphere soil in Genhe, July; RGO, rhizosphere soil in Genhe, October.
We found that the differences in fungal composition were significant (Figures 2B, C), and there were also significant differences in functional guilds between the groups (Figure 2D). Basidiomycota, Ascomycota, Mortierellomycota, and Mucoromycota were the dominant phyla (average relative abundance > 1%). There were 13 dominant families, among which the families with a relative abundance more than 5% included Thelephoraceae, Mortierellaceae, Archaeorhizomycoaceae, and Inocybaceae. In terms of functional guilds, symbiotrophic fungi had an absolute advantage in the identified functions and the relative abundances of pathotrophic and saprotrophic fungi in Aershan were significantly higher than those in Genhe.
More families were differentially expressed between sites (Figure 3C), while there were fewer families with differential expression in different compartments and seasons (Figures 3A, B). Clavulinaceae, Hoehnelomycetaceae, Hydnangiaceae, Hygrophoraceae, Pyronemataceae, Rhizopogonaceae, Sporidiobolaceae, and Tritirachiaceae were significantly express families in bulk soil. Mytilinidiaceae, Russulaceae, Trichosporonaceae, and Tulasnellaceae were significantly expressed in rhizosphere soil. In July, only Archaeorhizomycetaceae and Trimorphomycetaceae were significantly expressed. As for October, Clavulinaceae, Cordycipitaceae, Leucosporidiaceae, Mortierellaceae, Mycosphaerellaceae, Pezizaceae, Schizoporaceae, and Vibrisseaceae were significantly expressed. A total of 69 families showed differential expression between these two sites, with 40 families dominating at Aershan and 29 families dominating at Genhe (Table 3).
Figure 3 Differential expression analyses at family level in different compartments, seasons, and sites. (A) The ratio of rhizosphere to bulk. (B) The ratio of October to July. (C) The ratio of Genhe to Aershan. BAJ, bulk soil in Aershan, July; BAO, bulk soil in Aershan, October; BGJ, bulk soil in Genhe, July; BGO, bulk soil in Genhe, October; RAJ, rhizosphere soil in Aershan, July; RAO, rhizosphere soil in Aershan, October; RGJ, rhizosphere soil in Genhe, July; RGO, rhizosphere soil in Genhe, October.
Table 3 The families of differential expression analyses in different compartments, seasons, and sites.
These results indicated that there were significant differences in community composition and functional guilds between different compartments, seasons, and sites.
3.3 Fungal diversity and community assembly processes among different compartments, seasons, and sites
The differences between the Shannon diversity, Chao1 richness, and Pielou evenness indices were significant in different groups, which consist of compartments, seasons, and sites. Shannon diversity and Pielou evenness indices of RGO were highest among these groups (Figures 4A, C). The group with the highest Chao1 richness index is BAJ (Figures 4B). Three-way ANOVA was used to analyze the effects of compartment, season, and site on alpha diversity. Among the single factors by Three-way ANOVA (Table 4), the compartment had a significant effect on the Shannon (F = 10.179, P = 0.002) and Pielou (F = 12.58, P = 0.001) indices, whereas both season (F = 24.934, P < 0.0001) and site (F = 85.252, P < 0.0001) had a significant effect on the Chao1 index.
Figure 4 Differences in soil fungal communities among different compartments, seasons, and sites. The Shannon diversity index (A), Chao1 richness index (B) and Pielou evenness index (C) of the fungal community. (D) Principal coordinate analysis (PCoA) based on the Bray-Curtis distance of fungal community. (E) The β-Nearest Taxon Index (βNTI) of the communities. (F) The proportion of deterministic and stochastic processes in fungal community assembly process in the two sites. Different lower case letters indicate significant difference among groups (P < 0.05). BAJ, bulk soil in Aershan, July; BAO, bulk soil in Aershan, October; BGJ, bulk soil in Genhe, July; BGO, bulk soil in Genhe, October; RAJ, rhizosphere soil in Aershan, July; RAO, rhizosphere soil in Aershan, October; RGJ, rhizosphere soil in Genhe, July; RGO, rhizosphere soil in Genhe, October.
Table 4 Effect of compartment, season, site and their interactions on fungal alpha diversity index by Three-way ANOVA.
Principal coordinate analysis (PCoA) based on the Bray–Curtis distance (Figure 4D) indicated significant differences among multiple groups, especially between different sites (R² = 0.202, P = 0.001). The difference in fungal community structure between the rhizosphere and bulk soil in Genhe was relatively small, whereas the difference between the two seasons was significant. According to the results of the permutational multivariate analysis of variance (PERMANOVA), the grouping was credible (F = 9.088, P = 0.001).
According to the result of the β-nearest taxon index (βNTI, Figure 4E), the values were dispersed around −2. Homogeneous selection was more important than other processes in Genhe fungal community assembly, whereas dispersal limitation was more important than other processes in Aershan (Figure 4F). Soil fungal communities in Aershan were dominated by stochastic processes with average relative importance of 70.48−84.76%, whereas those in Genhe were dominated by deterministic processes with average relative importance of 59.05−68.57%. The proportions of deterministic and stochastic processes in the fungal community assembly process significantly differed between the two sites.
According to the results, there were significant differences in alpha diversity and beta diversity of fungal communities in the soil of L. gmelinii. Differently dominant processes in the community assembly may cause the above results.
3.4 The relationship between fungal community and environmental factors in bulk soil
We performed mantle tests on the relationship between individual environmental factors and the fungal community structure of bulk soil and found that all soil properties, temperature, and precipitation significantly affected the fungal community structure (Table 5). Among the environmental factors, TN (r = 0.734, P = 0.001) and pH (r = 0.174, P = 0.001) had the greatest effect on fungal community structure. Altogether, the correlation between the total environmental factors and community structure was not high (r = 0.174, P = 0.001).
Table 5 The Spearman’s correlations (r) between single environmental factor (Euclidean distance) and fungal community structure (Bray-Curtis distance) of the bulk soil determined by the Mantel test.
In the correlation relationships between alpha diversity indices and environmental factors (Figure 5A), climatic factors had a significant effect on all three alpha diversity indices, and alpha diversity increased as temperature and precipitation increased. Additionally, the Chao1 evenness index was significantly influenced by the CEC, TN, pH, SOC, and AP. The effects of the CEC and pH on the Shannon diversity index were also significant.
Figure 5 The impact of environmental factors on fungal communities. The Spearman’s correlation relationships between alpha diversity indices (A), the relative abundance of dominate families (B) and environmental factors. *, P < 0.05; **, P < 0.01; ***, P < 0.001. TN, total nitrogen; SOC, soil organic carbon; AP, available phosphorus; CEC: cation exchange capacity; Tmax, the monthly minimum temperature; Tmin, the monthly minimum temperature; Prec: mean monthly precipitation.
To determine which taxa dominate the variation in fungal community structure in L. gmelinii forests under environmental changes, we evaluated the relationship between fungal families and the environment. In total, 16 of 17 dominant families were significantly correlated with environmental factors (Figure 5B), whereas the correlation between Hygrophoraceae and environmental factors was low and not significant. Russulaceae was significantly and negatively correlated with climatic factors. Inocybaceae were only significantly correlated with soil physicochemical properties. Mortierellaceae, Atheliaceae, Pyronemataceae, Pseudeurotiaceae, Hyaloscyphaceae, Tricholomataceae, and Thelephoraceae significantly correlated with AP, CEC, TN, pH, Tmax, and Tmin. A few of them were significantly correlated with Prec or SOC. Umbelopsidaceae was significantly and positively correlated with CEC, TN, and climatic factors. Archaeorhizomycetaceae and Serendipi,taceae were significantly correlated with environmental factors except Prec. Myxotrichaceae, Cortinariaceae, and Clavulinaceae were correlated with some soil physicochemical properties and Prec. The 17 dominant families were divided into two categories based on their correlations as a result of clustering.
Our results demonstrated that soil physicochemical properties and climatic factors played key roles in the diversity, composition and community structure of soil fungi.
4 Discussion
4.1 The fungal community of rhizosphere soil
Hosts produce exudates to regulate the rhizosphere microenvironment and provide carbon to microorganisms (Zhuang et al., 2020). Fungi combined with host roots to assist plants in absorbing carbon, nitrogen, phosphorus, and nutrients and play an important role in plant growth, development (Lu et al., 2018), and health (Berendsen et al., 2012). Many factors affect rhizosphere fungal communities, including vegetation type (Philippot et al., 2013), soil properties (Li et al., 2016), and biotic interactions (Paungfoo-Lonhienne et al., 2015). In our study, the fluctuations in Shannon and Chao1 indices were smaller in rhizosphere soil than in bulk soil. We supposed that the composition of fungal community in the rhizosphere was more stable. The alpha diversity of the fungal community in rhizosphere of L. gmelinii forest in Genhe was higher than that in bulk soil. However, the result was completely opposite in the soybean soil (Zhang et al., 2018). The selective effects of the root on specific taxa caused all the differences (Philippot et al., 2013). There were four families enriched in the rhizosphere soil of L. gmelinii. Russulaceae and Tulasnellaceae formed obligate mutualistic root symbionts with L. gmelinii, and they were also widely distributed in the mycorrhizae of orchid species (Li et al., 2021). Mytilinidiaceae and Trichosporonaceae were mainly involved in soil saprophytic processes and animal pathogenic transmission (Lourenço et al., 2020; Eagar et al., 2023). The families enriched in rhizosphere suggesting that those fungi were active in assisting nutrient uptake and transforming organic matter for L. gmelinii. Our results further confirmed that the host root system played a decisive role in soil fungi, and selected the specific fungal species. Exploring the status of rhizosphere fungal communities helps us understand the root health status of plants.
4.2 The effect of seasonal changes on soil fungal communities
Seasonal changes in temperature and precipitation have important effects on aboveground and underground organisms. A previous study has shown that fungal communities have better adaptability and resistance to cold environments compared with bacterial communities (Buckeridge et al., 2013). Regardless of seasonal changes, Basidiomycetes and Ascomycetes are the main components of fungal communities (Xie and Yin, 2022). However, Zeng et al. (2020) reported that the relative abundance of Ascomycota decreases as the annual average precipitation increases. Previous studies have shown that saprophytic fungi dominate the soil during the cold season, whereas ectomycorrhizal fungi have a greater advantage during the growing season (Voříšková et al., 2014; Santalahti et al., 2016; Zeng et al., 2022). At our sampling sites, temperature and precipitation were concentrated in July; the temperature in October was below zero and the soil was relatively dry. High soil moisture promotes the decomposition of dead leaves and increases soil nutrients, thereby affecting the functional groups of the soil fungal community (Brockett et al., 2012). During the drier season, soil fungal abundance decreases compared to the wet season, but fungal diversity increases (Kivlin and Hawkes, 2016). In our previous study, variations in temperature and precipitation were important factors influencing seasonal changes in the fungal communities (Zhao et al., 2023). In fact, seasonal changes were masked by sites in this study.
4.3 The effect of community assembly on fungal community
Deterministic and stochastic processes were used to determine soil fungal community structure. Selection, drift, speciation, and dispersal are usually responsible for the deterministic fitness differences, species abundance, creates new species, and movement by the four processes, among which selection is a deterministic process and diffusion is a stochastic process (Vellend, 2010). The fungal community assembly in Aershan was dominated by stochastic processes, whereas it was dominated by deterministic processes in Genhe. The body size is crucial for dispersal (Luan et al., 2020), it is also an important stochastic factor that leads to the low similarity of the two communities in our study. The assembly process of different communities can cause differences in the community structure. Nevertheless, the role of stochastic processes in microbial community assembly remains unclear (Nemergut et al., 2013; Gao et al., 2020; Jiao et al., 2020). Combining niche breadth, deterministic processes have a strong impact on specialist species with narrow niche breadth, while the generalist species with wide niche breadth has stronger competitiveness (Graham et al., 2016). In general, the community dominated by stochastic process has more resources around it and less competitive pressure (Jiao et al., 2020). Based on this understanding of community assembly, we consider that the soil condition of Genhe was more limited, resulting in contemporary selection dominating the fungal communities. The soil conditions in Aershan were better, providing more opportunities for fungal communities to drift, diversify, and disperse.
4.4 The effect of contemporary environment heterogeneity and historical contingencies on the structure of fungal communities
Notably, contemporary environmental heterogeneity and historical events (geographic distance) were the two major factors affecting the biographical distribution of soil fungal communities (Ramette and Tiedje, 2007). The two sampling sites for L. gmelinii were 460 km apart, and significant differences were found in the soil physicochemical properties between the two sites, so we considered the effect of geographical distance. In previous studies, the geographical scale played an important role in fungal community structure (Tedersoo et al., 2014; Zhang et al., 2018; Adamo et al., 2021). He et al. (2017) even suggested that spatial heterogeneity has a more significant effect than temporal consistency.
In the biogeographical distribution of soil fungi, the contemporary environment played a more important role (Zeng et al., 2019). As important soil nutrient factors, C, N, and P not only have a profound effect on aboveground vegetation but also provide an environment for the living dynamics of soil fungi. Nutritional factors such as TN, AP, and SOC are important in the regulation of fungal communities (Lauber et al., 2008; He et al., 2016; Dang et al., 2017; Li et al., 2020; Luo et al., 2023). Long-term nitrogen enrichment decreases the relative abundance and diversity of ectomycorrhizal taxa; however, the saprotrophs respond positively to nitrogen enrichment. (Zhou et al., 2016; Morrison et al., 2016). Nitrogen deposition affects the soil fungal community in L. gmelinii forests annually and is a major cause of global environmental change (Galloway et al., 2008). SOC significantly affects the structure of soil fungal communities (Sul et al., 2013; Li et al., 2017). A few families in our study were significantly related to SOC; however, all related families belonged to saprophytic and symbiotic fungi, such as ericoid mycorrhizal fungi, orchid mycorrhizal fungi, and ectomycorrhizal fungi. They degrade cellulose (Li et al., 2017), and assist hosts in obtaining carbon (Suetsugu et al., 2021). Nevertheless, when the other four physicochemical properties were significantly correlated with alpha diversity indices and some dominant families, AP showed a significant opposite correlation (Figure 5). Similar to our results for paddy fields in subtropical China, fungal richness and diversity indices declined significantly with increasing phosphorus rates (Liu et al., 2018). Bao et al. (2013) also found that the relationship between available phosphorus and fungal biodiversity showed a hump trend, slowly increasing before peaking and rapidly decreasing thereafter. As a non-nutritional factor, pH has a significant effect on soil fungal communities (Tedersoo et al., 2020). Compared with bacteria, fungi have a wider pH adaptation range (Rousk et al., 2010), but the dominant fungal taxa in different soil pH environments differ (Zhang et al., 2016). In our study, the correlation between the most dominant families and pH was high and extremely significant. CEC also affects the structure of fungal communities (Docherty et al., 2015). Some symbiotic taxa and the family Archaeorhiziomyceae were negatively correlated with pH, CEC, TN, and SOC. The relationship between these fungi and their hosts allows them to gain advantages in poor and acidic soils (Rosling et al., 2011; Vohník et al., 2016; Zhang and Sun, 2021). In addition to the factors studied, soil moisture content (Docherty et al., 2015) and salinity (Cao et al., 2021) also play important roles.
There were close relationships between TN, SOC, AP, pH, CEC and community structure, with TN and pH having the greatest effect (Table 5). The soil total nitrogen of our samples ranged from 0.81 g/kg to 6.88 g/kg, and the pH values ranged from 3.71 to 6.23, which was consistent with the values in the study of Miyamoto et al. (2021). We did not conduct further experiments on its optimal gradient, but pH around 4.89 (4.89 ± 0.08) may be its superior pH by ANOVA. Combined with the results that soil total nitrogen was positively correlated with alpha diversity indices (including Shannon diversity index, Chao1 richness index and Pielou evenness index), we supposed that high soil total nitrogen is more advantageous for soil fungal communities. So appropriate nitrogen supplementation and soil pH regulation can maintain and regulate the health of L. gmelinii forest.
5 Conclusion
Our findings revealed that compared to compartment and season, the impact of site on the soil fungal community of L. gmelinii was more important. Firstly, geographical distance determined low community similarity in terms of the stochastic process of community assembly. Secondly, there were significant differences in soil physicochemical properties between the two sites, which belonged to deterministic process. Environmental factors significantly affected fungal community structure, with soil TN and pH being the most important influencing factors. The geographical distance and environmental differences between the two sites result in significant differences in the soil fungal community assembly process. Meanwhile, seasonal changes appeared insignificant. This further deepens our understanding of the soil fungi of L. gmelinii forest and grasps the important soil environment in which L. gmelinii lives, which helps us to maintain the health of L. gmelinii.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.
Author contributions
WZ: Data curation, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing. KH: Data curation, Investigation, Methodology, Software, Writing – review & editing. RM: Data curation, Investigation, Methodology, Writing – review & editing. JL: Data curation, Investigation, Methodology, Writing – review & editing. YS: Conceptualization, Supervision, Writing – original draft, Writing – review & editing. BC: Conceptualization, Funding acquisition, Investigation, Project administration, Supervision, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Natural Science Foundation of China (Nos. 32325001, 32270010), the Fundamental Research Funds for the Central Universities (No. QNTD202307).
Acknowledgments
We would like to thank Dan-Dan Wang, Shun Liu, and Kun Qian for their generous help with the field sampling. We thank the reviewers and the editor for their constructive comments through the revision process of this manuscript.
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.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Adamo, I., Castaño, C., Bonet, J. A., Colinas, C., Martínez de Aragón, J., Alday, J. G. (2021). Soil physico-chemical properties have a greater effect on soil fungi than host species in Mediterranean pure and mixed pine forests. Soil Biol. Biochem. 160, 108320. doi: 10.1016/j.soilbio.2021.108320
Bao, Z., Matsushita, Y., Morimoto, S., Hoshino, Y. T., Suzuki, C., Nagaoka, K., et al. (2013). Decrease in fungal biodiversity along an available phosphorous gradient in arable Andosol soils in Japan. Can. J. Microbiol. 59, 368–373. doi: 10.1139/cjm-2012-0612
Berendsen, R. L., Pieterse, C. M. J., Bakker, P. A. H. M. (2012). The rhizosphere microbiome and plant health. Trends Plant Sci. 17, 478–486. doi: 10.1016/j.tplants.2012.04.001
Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C. C., Al-Ghalith, G. A., et al. (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37, 852–857. doi: 10.1038/s41587-019-0209-9
Brockett, B. F. T., Prescott, C. E., Grayston, S. J. (2012). Soil moisture is the major factor influencing microbial community structure and enzyme activities across seven biogeoclimatic zones in western Canada. Soil Biol. Biochem. 44, 9–20. doi: 10.1016/j.soilbio.2011.09.003
Broeckling, C. D., Broz, A. K., Bergelson, J., Manter, D. K., Vivanco, J. M. (2008). Root exudates regulate soil fungal community composition and diversity. Appl. Environ. Microbiol. 74, 738–744. doi: 10.1128/AEM.02188-07
Buckeridge, K. M., Banerjee, S., Siciliano, S. D., Grogan, P. (2013). The seasonal pattern of soil microbial community structure in mesic low arctic tundra. Soil Biol. Biochem. 65, 338–347. doi: 10.1016/j.soilbio.2013.06.012
Cai, Z. Q., Zhang, Y. H., Yang, C., Wang, S. (2018). Land-use type strongly shapes community composition, but not always diversity of soil microbes in tropical China. Catena 165, 369–380. doi: 10.1016/j.catena.2018.02.018
Cao, C., Tao, S., Cui, Z., Zhang, Y. (2021). Response of soil properties and microbial communities to increasing salinization in the meadow grassland of northeast China. Microb. Ecol. 82, 722–735. doi: 10.1007/s00248-021-01695-x
Crowther, T. W., van den Hoogen, J., Wan, J., Mayes, M. A., Keiser, A. D., Mo, L., et al. (2019). The global soil community and its influence on biogeochemistry. Science 365, eaav0550. doi: 10.1126/science.aav0550
Dang, P., Yu, X., Le, H., Liu, J., Shen, Z., Zhao, Z. (2017). Effects of stand age and soil properties on soil bacterial and fungal community composition in Chinese pine plantations on the Loess Plateau. PloS One 12, e0186501. doi: 10.1371/journal.pone.0186501
Ding, L. Z., Man, X. L., Xiao, R. H. (2019). Characteristics of nitrogen content in rhizosphere soil of main tree species in northern part of Daxinganling during growing seasons. J. Cent. South Univ. For. Technol. 39, 65–71 + 92. doi: 10.14067/j.cnki.1673–923x.2019.02.011
Docherty, K. M., Borton, H. M., Espinosa, N., Gebhardt, M., Gil-Loaiza, J., Gutknecht, J. L. M., et al. (2015). Key edaphic properties largely explain temporal and geographic variation in soil microbial communities across four biomes. PloS One 10, e0135352. doi: 10.1371/journal.pone.0135352
Eagar, A. C., Smemo, K. A., Phillips, R. P., Blackwood, C. B. (2023). Context-dependence of fungal community responses to dominant tree mycorrhizal types in Northern hardwood forests. Soil Biol. Biochem. 178, 108971. doi: 10.1016/j.soilbio.2023.108971
Fouquier, J., Rideout, J. R., Bolyen, E., Chase, J., Shiffer, A., McDonald, D., et al. (2016). Ghost-tree: creating hybrid-gene phylogenetic trees for diversity analyses. Microbiome 4, 1–10. doi: 10.1186/s40168-016-0153-6
Galloway, J. N., Townsend, A. R., Erisman, J. W., Bekunda, M., Cai, Z., Freney, J. R., et al. (2008). Transformation of the nitrogen cycle: recent trends, questions, and potential solutions. Science 320, 889–892. doi: 10.1126/science.1136674
Gao, C., Montoya, L., Xu, L., Madera, M., Hollingsworth, J., Purdom, E., et al. (2020). Fungal community assembly in drought-stressed sorghum shows stochasticity, selection, and universal ecological dynamics. Nat. Commun. 11, 34. doi: 10.1038/s41467-019-13913-9
Graham, E. B., Crump, A. R., Resch, C. T., Fansler, S., Arntzen, E., Kennedy, D. W., et al. (2016). Coupling spatiotemporal community assembly processes to changes in microbial metabolism. Front. Microbiol. 7. doi: 10.3389/fmicb.2016.01949
He, D., Xiang, X., He, J. S., Wang, C., Cao, G., Adams, J., et al. (2016). Composition of the soil fungal community is more sensitive to phosphorus than nitrogen addition in the alpine meadow on the Qinghai-Tibetan Plateau. Biol. Fertil. Soils 52, 1059–1072. doi: 10.1007/s00374-016-1142-4
He, J., Tedersoo, L., Hu, A., Han, C., He, D., Wei, H., et al. (2017). Greater diversity of soil fungal communities and distinguishable seasonal variation in temperate deciduous forests compared with subtropical evergreen forests of eastern China. FEMS Microbiol. Ecol. 93, fix069. doi: 10.1093/femsec/fix069
Huang, K., Guo, Z., Zhao, W., Song, C., Wang, H., Li, J., et al. (2023). Response of fungal communities to afforestation and its indication for forest restoration. For. Ecosyst. 10, 100125. doi: 10.1016/j.fecs.2023.100125
Ishida, T. A., Nara, K., 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
IUSS Working Group WRB (2015). “World Reference Base for soil resources 2014, update 2015 international soil classification system for naming soils and creating legends for soil maps,” in World Soil Resources Reports No. 106 (FAO, Rome), 106.
Jiang, S., Xing, Y., Liu, G., Hu, C., Wang, X., Yan, G., et al. (2021). Changes in soil bacterial and fungal community composition and functional groups during the succession of boreal forests. Soil Biol. Biochem. 161, 108393. doi: 10.1016/j.soilbio.2021.108393
Jiao, S., Yang, Y., Xu, Y., Zhang, J., Lu, Y. (2020). Balance between community assembly processes mediates species coexistence in agricultural soil microbiomes across eastern China. ISME J. 14, 202–216. doi: 10.1038/s41396-019-0522-9
Kajimoto, T., Matsuura, Y., Sofronov, M. A., Volokitina, A. V., Mori, S., Osawa, A., et al. (1999). Above-and belowground biomass and net primary productivity of a Larix gmelinii stand near Tura, central Siberia. Tree Physiol. 19, 815–822. doi: 10.1093/treephys/19.12.815
Kivlin, S. N., Hawkes, C. V. (2016). Tree species, spatial heterogeneity, and seasonality drive soil fungal abundance, richness, and composition in Neotropical rainforests. Environ. Microbiol. 18, 4662–4673. doi: 10.1111/1462-2920.13342
Kuang, J., Han, S., Chen, Y., Bates, C. T., Wang, P., Shu, W. (2021). Root-associated fungal community reflects host spatial co-occurrence patterns in a subtropical forest. ISME Commun. 1, 65. doi: 10.1038/s43705-021-00072-6
Lauber, C. L., Strickland, M. S., Bradford, M. A., Fierer, N. (2008). The influence of soil properties on the structure of bacterial and fungal communities across land-use types. Soil Biol. Biochem. 40, 2407–2415. doi: 10.1016/j.soilbio.2008.05.021
Li, S., Huang, X., Shen, J., Xu, F., Su, J. (2020). Effects of plant diversity and soil properties on soil fungal community structure with secondary succession in the Pinus yunnanensis forest. Geoderma 379, 114646. doi: 10.1016/j.geoderma.2020.114646
Li, Y., Li, Y., Chang, S. X., Liang, X., Qin, H., Chen, J., et al. (2017). Linking soil fungal community structure and function to soil organic carbon chemical composition in intensively managed subtropical bamboo forests. Soil Biol. Biochem. 107, 19–31. doi: 10.1016/j.soilbio.2016.12.024
Li, J. Y., Zhao, Z. Y., Gao, Y., Xing, X. K. (2021). Mycorrhizal fungal community composition of medicinal plant Spiranthes sinensis (Orchidaceae). Mycosystema 40, 1317–1327. doi: 10.13346/j.mycosystema.210033
Li, Z., Zu, C., Wang, C., Yang, J., Yu, H., Wu, H. (2016). Different responses of rhizosphere and non-rhizosphere soil microbial communities to consecutive Piper nigrum L. monoculture. Sci. Rep. 6, 35825. doi: 10.1038/srep35825
Liang, Y., Jiang, Y., Wang, F., Wen, C., Deng, Y., Xue, K., et al. (2015). Long-term soil transplant simulating climate change with latitude significantly alters microbial temporal turnover. ISME J. 9, 2561–2572. doi: 10.1038/ismej.2015.78
Liu, H., Li, J., Carvalhais, L. C., Percy, C. D., Prakash Verma, J., Schenk, P. M., et al. (2021). Evidence for the plant recruitment of beneficial microbes to suppress soil-borne pathogens. New Phytol. 229, 2873–2885. doi: 10.1111/nph.17057
Liu, M., Liu, J., Chen, X., Jiang, C., Wu, M., Li, Z. (2018). Shifts in bacterial and fungal diversity in a paddy soil faced with phosphorus surplus. Biol. Fertil. Soils 54, 259–267. doi: 10.1007/s00374-017-1258-1
Liu, L., Ma, L., Zhu, M., Liu, B., Liu, X., Shi, Y. (2023). Rhizosphere microbial community assembly and association networks strongly differ based on vegetation type at a local environment scale. Front. Microbiol. 14. doi: 10.3389/fmicb.2023.1129471
Liu, Y., Trancoso, R., Ma, Q., Yue, C., Wei, X., Blanco, J. A. (2020). Incorporating climate effects in Larix gmelinii improves stem taper models in the Greater Khingan Mountains of Inner Mongolia, northeast China. For. Ecol. Manage. 464, 118065. doi: 10.1016/j.foreco.2020.118065
Lladó, S., López-Mondéjar, R., Baldrian, P. (2018). Drivers of microbial community structure in forest soils. Appl. Microbiol. Biotechnol. 102, 4331–4338. doi: 10.1007/s00253-018-8950-4
Lourenço, K. S., Suleiman, A. K. A., Pijl, A., Cantarella, H., Kuramae, E. E. (2020). Dynamics and resilience of soil mycobiome under multiple organic and inorganic pulse disturbances. Sci. Total Environ. 733, 139173. doi: 10.1016/j.scitotenv.2020.139173
Lu, T., Ke, M., Lavoie, M., Jin, Y., Fan, X., Zhang, Z., et al. (2018). Rhizosphere microorganisms can influence the timing of plant flowering. Microbiome 6, 231. doi: 10.1186/s40168-018-0615-0
Luan, L., Jiang, Y., Cheng, M., Dini-Andreote, F., Sui, Y., Xu, Q., et al. (2020). Organism body size structures the soil microbial and nematode community assembly at a continental and global scale. Nat. Commun. 11, 6406. doi: 10.1038/s41467-020-20271-4
Luo, H., Wang, C., Zhang, K., Ming, L., Chu, H., Wang, H. (2023). Elevational changes in soil properties shaping fungal community assemblages in terrestrial forest. Sci. Total Environ. 900, 165840. doi: 10.1016/j.scitotenv.2023.165840
Maron, J. L., Marler, M., Klironomos, J. N., Cleveland, C. C. (2011). Soil fungal pathogens and the relationship between plant diversity and productivity. Ecol. Lett. 14, 36–41. doi: 10.1111/j.1461-0248.2010.01547.x
Meng, Z. Y., Li, M., Yang, X. J., Lv, G. F., Niu, Y. F., Zhang, Y. P. (2023). Diversity, community composition and influencing factors of ectomycorrhizal fungi in typical forest types in cold temperate regions. Acta Ecol. Sin. 43, 38–47. doi: 10.5846/stxb202206071619
Miyamoto, Y., Danilov, A. V., Bryanin, S. V. (2021). The dominance of Suillus species in ectomycorrhizal fungal communities on Larix gmelinii in a post-fire forest in the Russian Far East. Mycorrhiza 31, 55–66. doi: 10.1007/s00572-020-00995-3
Morrison, E. W., Frey, S. D., Sadowsky, J. J., van Diepen, L. T. A., Thomas, W. K., Pringle, A. (2016). Chronic nitrogen additions fundamentally restructure the soil fungal community in a temperate forest. Fungal Ecol. 23, 48–57. doi: 10.1016/j.funeco.2016.05.011
Nemergut, D. R., Schmidt, S. K., Fukami, T., O’Neill, S. P., Bilinski, T. M., Stanish, L. F., et al. (2013). Patterns and processes of microbial community assembly. Microbiol. Mol. Biol. Rev. 77, 342–356. doi: 10.1128/MMBR.00051-12
Nguyen, N. H., Song, Z., Bates, S. T., Branco, S., Tedersoo, L., Menke, J., et al. (2016). FUNGuild: An open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 20, 241–248. doi: 10.1016/j.funeco.2015.06.006
Nilsson, R. H., Larsson, K. H., Taylor, A. F. S., Bengtsson-Palme, J., Jeppesen, T. S., Schigel, D., et al. (2019). The UNITE database for molecular identification of fungi: Handling dark taxa and parallel taxonomic classifications. Nucleic Acids Res. 47, D259–D264. doi: 10.1093/nar/gky1022
Orwin, K. H., Kirschbaum, M. U. F., St John, M. G., Dickie, I. A. (2011). Organic nutrient uptake by mycorrhizal fungi enhances ecosystem carbon storage: A model-based assessment. Ecol. Lett. 14, 493–502. doi: 10.1111/ele.2011.14.issue-5
Ouimette, A., Guo, D., Hobbie, E., Gu, J. (2013). Insights into root growth, function, and mycorrhizal abundance from chemical and isotopic data across root orders. Plant Soil 367, 313–326. doi: 10.1007/s11104-012-1464-4
Paungfoo-Lonhienne, C., Yeoh, Y. K., Kasinadhuni, N. R. P., Lonhienne, T. G. A., Robinson, N., Hugenholtz, P., et al. (2015). Nitrogen fertilizer dose alters fungal communities in sugarcane soil and rhizosphere. Sci. Rep. 5, 8678. doi: 10.1038/srep08678
Philippot, L., Raaijmakers, J. M., Lemanceau, P., van der Putten, W. H. (2013). Going back to the roots: The microbial ecology of the rhizosphere. Nat. Rev. Microbiol. 11, 789–799. doi: 10.1038/nrmicro3109
Powell, J. R., Karunaratne, S., Campbell, C. D., Yao, H., Robinson, L., Singh, B. K. (2015). Deterministic processes vary during community assembly for ecologically dissimilar taxa. Nat. Commun. 6, 8444. doi: 10.1038/ncomms9444
Praeg, N., Illmer, P. (2020). Microbial community composition in the rhizosphere of Larix decidua under different light regimes with additional focus on methane cycling microorganisms. Sci. Rep. 10, 22324. doi: 10.1038/s41598–020-79143-y
Ramette, A., Tiedje, J. M. (2007). Multiscale responses of microbial life to spatial distance and environmental heterogeneity in a patchy ecosystem. Proc. Natl. Acad. Sci. 104, 2761–2766. doi: 10.1073/pnas.0610671104
R Core Team (2022). R: A language and environment for statistical computing (Vienna, Austria: R foundation for statistical computing).
Rosling, A., Cox, F., Cruz-Martinez, K., Ihrmark, K., Grelet, G. A., Lindahl, B. D., et al. (2011). Archaeorhizomycetes: unearthing an ancient class of ubiquitous soil fungi. Science 333, 876–879. doi: 10.1126/science.1206958
Rousk, J., Bååth, E., Brookes, P. C., Lauber, C. L., Lozupone, C., Caporaso, J. G., et al. (2010). Soil bacterial and fungal communities across a pH gradient in an arable soil. ISME J. 4, 1340–1351. doi: 10.1038/ismej.2010.58
Santalahti, M., Sun, H., Jumpponen, A., Pennanen, T., Heinonsalo, J. (2016). Vertical and seasonal dynamics of fungal communities in boreal Scots pine forest soil. FEMS Microbiol. Ecol. 92, fiw170. doi: 10.1093/femsec/fiw170
Setälä, H., McLean, M. A. (2004). Decomposition rate of organic substrates in relation to the species diversity of soil saprophytic fungi. Oecologia 139, 98–107. doi: 10.1007/s00442-003-1478-y
Shi, L. L., Mortimer, P. E., Ferry Slik, J. W., Zou, X. M., Xu, J., Feng, W. T., et al. (2014). Variation in forest soil fungal diversity along a latitudinal gradient. Fungal Divers. 64, 305–315. doi: 10.1007/s13225-013-0270-5
Stegen, J. C., Lin, X., Fredrickson, J. K., Chen, X., Kennedy, D. W., Murray, C. J., et al. (2013). Quantifying community assembly processes and identifying features that impose them. ISME J. 7, 2069–2079. doi: 10.1038/ismej.2013.93
Stegen, J. C., Lin, X., Konopka, A. E., Fredrickson, J. K. (2012). Stochastic and deterministic assembly processes in subsurface microbial communities. ISME J. 6, 1653–1664. doi: 10.1038/ismej.2012.22
Suetsugu, K., Haraguchi, T. F., Okada, H., Tayasu, I. (2021). Stigmatodactylus sikokianus (Orchidaceae) mainly acquires carbon from decaying litter through association with a specific clade of Serendipitaceae. New Phytol. 231, 1670–1675. doi: 10.1111/nph.17523
Sul, W. J., Asuming-Brempong, S., Wang, Q., Tourlousse, D. M., Penton, C. R., Deng, Y., et al. (2013). Tropical agricultural land management influences on soil microbial communities through its effect on soil organic carbon. Soil Biol. Biochem. 65, 33–38. doi: 10.1016/j.soilbio.2013.05.007
Sun, Y. J., Zhang, J., Han, A. H., Wang, X. J., Wang, X. J. (2007). Biomass and carbon pool of Larix gmelini young and middle age forest in Xing’an mountains Inner Mongolia. Acta Ecol. Sin. 27, 1756–1762. doi: 10.3321/j.issn:1000–0933.2007.05.011
Tedersoo, L., Anslan, S., Bahram, M., Drenkhan, R., Pritsch, K., Buegger, F., et al. (2020). Regional-scale in-depth analysis of soil fungal diversity reveals strong pH and plant species effects in northern Europe. Front. Microbiol. 11. doi: 10.3389/fmicb.2020.01953
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
Vellend, M. (2010). Conceptual synthesis in community ecology. Q Rev. Biol. 85, 183–206. doi: 10.1086/652373
Vohník, M., Pánek, M., Fehrer, J., Selosse, M. A. (2016). Experimental evidence of ericoid mycorrhizal potential within Serendipitaceae (Sebacinales). Mycorrhiza 26, 831–846. doi: 10.1007/s00572-016-0717-0
Voříšková, J., Brabcová, V., Cajthaml, T., Baldrian, P. (2014). Seasonal dynamics of fungal communities in a temperate oak forest soil. New Phytol. 201, 269–278. doi: 10.1111/nph.12481
Wang, D. D., Zhao, W., Reyila, M., Huang, K. C., Liu, S., Cui, B. K. (2022). Diversity of Microbial Communities of Pinus sylvestris var. mongolica at Spatial Scale. Microorganisms 10, 371. doi: 10.3390/microorganisms10020371
White, T. J., Bruns, T. D., Lee, S. B., Taylor, J. W. (1990). Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. PCR Protocols: A Guide to Methods and Applications (New York: Academic Press, Inc.), 315–322. doi: 10.1016/B978-0-12-372180-8.50042-1
Xie, L., Yin, C. (2022). Seasonal variations of soil fungal diversity and communities in subalpine coniferous and broadleaved forests. Sci. Total Environ. 846, 157409. doi: 10.1016/j.scitotenv.2022.157409
Zeng, Q., Jia, P., Wang, Y., Wang, H., Li, C., An, S. (2019). The local environment regulates biogeographic patterns of soil fungal communities on the Loess Plateau. Catena 183, 104220. doi: 10.1016/j.catena.2019.104220
Zeng, Q., Lebreton, A., Man, X., Jia, L., Wang, G., Gong, S., et al. (2022). Ecological drivers of the soil microbial diversity and composition in primary old-growth forest and secondary woodland in a subtropical evergreen broad-leaved forest biome in the Ailao Mountains, China. Front. Microbiol. 13. doi: 10.3389/fmicb.2022.908257
Zeng, Q., Liu, Y., Xiao, L., An, S. (2020). Climate and soil properties regulate soil fungal communities on the Loess Plateau. Pedobiologia. 81–82, 150668. doi: 10.1016/j.pedobi.2020.150668
Zhang, Y. H., Sun, L. F. (2021). Research advances on the mycorrhizas of Ericaceae plants. Mycosystema 40, 1299–1316. doi: 10.13346/j.mycosystema.210077
Zhang, T., Wang, N. F., Liu, H. Y., Zhang, Y. Q., Yu, L. Y. (2016). Soil pH is a key determinant of soil fungal community composition in the Ny-Ålesund Region, Svalbard (High Arctic). Front. Microbiol. 7. doi: 10.3389/fmicb.2016.00227
Zhang, J., Zhang, B., Liu, Y., Guo, Y., Shi, P., Wei, G. (2018). Distinct large-scale biogeographic patterns of fungal communities in bulk soil and soybean rhizosphere in China. Sci. Total Environ. 644, 791–800. doi: 10.1016/j.scitotenv.2018.07.016
Zhao, W., Wang, D. D., Huang, K. C., Liu, S., Reyila, M., Sun, Y. F., et al. (2023). Seasonal variation in the soil fungal community structure of Larix gmelinii forests in Northeast China. Front. Microbiol. 14. doi: 10.3389/fmicb.2023.1106888
Zhong, Y., Chu, C., Myers, J. A., Gilbert, G. S., Lutz, J. A., Stillhard, J., et al. (2021). Arbuscular mycorrhizal trees influence the latitudinal beta-diversity gradient of tree communities in forests worldwide. Nat. Commun. 12, 3137. doi: 10.1038/s4146
Zhou, J., Jiang, X., Zhou, B., Zhao, B., Ma, M., Guan, D., et al. (2016). Thirty four years of nitrogen fertilization decreases fungal diversity and alters fungal community composition in black soil in northeast China. Soil Biol. Biochem. 95, 135–143. doi: 10.1016/j.soilbio.2015.12.012
Zhou, J., Ning, D. (2017). Stochastic community assembly: does it matter in microbial ecology? Microbiol. Mol. Biol. Rev. 81, e00002–e00017. doi: 10.1128/MMBR.00002-17
Zhuang, W., Yu, X., Hu, R., Luo, Z., Liu, X., Zheng, X., et al. (2020). Diversity, function and assembly of mangrove root-associated microbial communities at a continuous fine-scale. NPJ Biofilms Microbiomes 6, 52. doi: 10.1038/s41522-020-00164-6
Keywords: climatic factors, fungal community assembly, Larix gmelinii forest, physicochemical properties, soil fungi
Citation: Zhao W, Huang K, Mumin R, Li J, Sun Y and Cui B (2024) Spatial variations impact the soil fungal communities of Larix gmelinii forests in Northeast China. Front. Plant Sci. 15:1408272. doi: 10.3389/fpls.2024.1408272
Received: 28 March 2024; Accepted: 29 April 2024;
Published: 24 May 2024.
Edited by:
Kai Sun, Nanjing Normal University, ChinaReviewed by:
Yu Shi, Henan University, ChinaXu Liu, Nanjing Agricultural University, China
Bo Zhu, Zhejiang Chinese Medical University, China
Copyright © 2024 Zhao, Huang, Mumin, Li, Sun and Cui. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yifei Sun, yifeisun2016@163.com; Baokai Cui, cuibaokai@bjfu.edu.cn