Effect of elevation on composition and diversity of fungi in the rhizosphere of a population of Deyeuxia angustifolia on Changbai Mountain, northeastern China

Soil fungi are a key component of terrestrial ecosystems and play a major role in soil biogeochemical cycling. Although the diversity and composition of fungal communities are regulated by many abiotic and biotic factors, the effect of elevation on soil fungal community diversity and composition remains largely unknown. In this study, the soil fungal composition and diversity in Deyeuxia angustifolia populations along an elevational gradient (1,690 m to 2020 m a.s.l.) were assessed, using Illumina MiSeq sequencing, on the north-facing slope of the Changbai Mountain, northeastern China. Our results showed that soil physicochemical parameters changed significantly along with the elevational gradients. The Ascomycota and Basidiomycota were the most dominant phyla along with the gradient. Alpha diversity of soil fungi decreased significantly with elevation. Soil nitrate nitrogen (NO₃⁻-N) was positively correlated with fungal richness and phylogenetic diversity (PD), indicating that soil nitrate nitrogen (NO₃⁻-N) is a key soil property determining fungal community diversity. In addition to soil nitrate content, soil pH and soil moisture were the most important environmental properties determining the soil fungal diversity. Our results suggest that the elevational changes in soil physicochemical properties play a key role in shaping the community composition and diversity of soil fungi. This study will allow us to better understand the biodiversity distribution patterns of soil microorganisms in mountain ecosystems.

1. Introduction

Mountains have attracted increasing curiosity of ecologists because of their high sensitivity to global climate change (Frey et al., 2016). Elevational gradients in the high mountains are often characterized by dramatic changes in the abiotic and biotic factors within short geographical distances (Rime et al., 2016; Donhauser and Frey, 2018; Adamczyk et al., 2019; Zhang et al., 2021), including vegetation, soil physicochemical properties, precipitation, temperature and illumination (Shen et al., 2013; Lin et al., 2017; Zhang et al., 2022), which may lead to significant changes in soil microorganisms within short distance. For example, soil bacterial community diversity varied significantly with elevation, while soil fungal community diversity did not change across an elevational gradient of 400 m (Yu et al., 2019), but both the soil fungal and bacterial community diversities changed significantly with elevation across an elevational gradient of 1,500 m in the Tibetan Plateau (Han et al., 2022).

In recent years, many microorganisms have been discovered in a variety of cold environments (Sahay et al., 2013; Shivaji et al., 2013; Prasad et al., 2014; Rime et al., 2016; Donhauser et al., 2020; Zong and Fu, 2021; Zhang et al., 2022). On the Changbai Mountain, northeastern China, many studies have investigated the structure and function of microbial community in forest soils (Zhou, 2006; Wang et al., 2013; Li et al., 2017; Ping et al., 2017), but no studies have investigated the elevational patterns of fungi in grassy soils on the Changbai Mountain. Deyeuxia angustifolia is a typical grass in the mountain ecosystems that plays an important role in biogeochemistry (Sui et al., 2021; Weng et al., 2022). This species is distributed on the Chnagbai Mountain from 1,690 m to 2020 m a.s.l., which provides an ideal field platform to study the diversity of soil fungi with an elevational gradient but within the same herbaceous plant population.

Soil fungi are key components, playing an important role in biogeochemical cycling and litter decomposition of terrestrial ecosystem, and are closely related to soil properties and aboveground vegetation community characteristics (Wang et al., 2013; Ping et al., 2017; Ni et al., 2018; Hanif et al., 2019; Ren et al., 2019; Sui et al., 2021; Yang et al., 2021; Zhou et al., 2021). Soil fungal community diversity showed a significant relationship with soil pH (Fouts et al., 2012; Zhou et al., 2021), C/N ratio, soil temperature and soil organic carbon (Ping et al., 2017; Sui et al., 2021; Zhou et al., 2021; Deng et al., 2023). Different ecosystems with different characteristics determine the characteristics of soil fungal communities (Shi et al., 2014). The traditional explanation for this phenomenon is that aboveground vegetation affects soil fungal communities by altering the physicochemical properties of the soils. Soil, vegetation, and climatic factors change gradually with increasing elevation, suggesting that the fungal community may vary along with elevational gradients on mountains (Tedersoo et al., 2016).

The Changbai Mountain is an important gene pool of biodiversity in Northeast China (Xue and Tisdell, 2001; Tang et al., 2011). Its rich species diversity has made it a research hotspot. As one of the main vegetation species on Changbai Mountain, D. angustifolia is distributed from 1,690 to 2020 m and is indispensable for the protection of ecosystem functions. According to previous reports, D. angustifolia in the alpine tundra are invaded from lower elevations due to climate change (Zong et al., 2013, 2014), indicating that the soil fungal community in the alpine tundra ecosystem could also correspondingly change. Because soil fungi play an important role in litter degradation and nutrient cycling (Frey et al., 2016; Deng et al., 2023), the changes in fungal community composition and diversity caused by D. angustifolia invasion may directly affect the ecosystem structure and function. Therefore, understanding the changes in soil fungal community composition and diversity under D. angustifolia at different elevations can help us predict changes in ecosystem structure and function following D. angustifolia invasion. Unfortunately, there are no comprehensive studies on soil fungi in D. angustifolia population along elevational gradients.

To comprehensively understand these changes, we investigated the composition and diversity of soil fungal community in D. angustifolia population at 1,690, 1,800, 1,910, and 2,020 m above sea level (a.s.l.) along an elevational gradient, using Illumina Miseq sequencing, on the Chnagbai Mountain. We hypothesize that the soil fungal composition and diversity in D. angustifolia population change significantly along the elevational gradient, because the soil environmental characteristics that shape the soil fungal composition and diversity change with increasing elevation. Therefore, the objectives of this study were (1) to compare the fungal diversity and community composition in D. angustifolia population in response to elevation, and (2) to evaluate the relationships between soil fungal communities and soil physicochemical properties across the elevational gradient.

2. Materials and methods

2.1. Research site

This study was performed on the Changbai Mountain (126^°55′-129^°00′E, 41^°23′-42^°36’N) in northeastern China. The local climate is a typical continental temperate monsoon climate with a daily average temperature of 5.9°C during the growing season (June to September). The average annual precipitation during the growing season can reach 958 mm. The mean annual precipitation and temperature is approximately 600 mm and 4°C, respectively.

To study the elevational pattern of soil fungal composition and diversity, we selected pure D. angustifolia population at 1,690, 1,810, 1,910, and to 2020 m a.s.l. along an elevational gradient on the north-facing slope of the Changbai Mountain. During October 1 to 7, 2018, three independent plots (10 m × 10 m) were set up in D. angustifolia population at each elevation. Ten to fifteen soil samples (0 ~ 20 cm organic layer) were sampled and pooled for each plot, using a sterile soil drill (5 cm in diameter, 20 cm deep). After removing the surface litter and humus layer, approximately 1 kg soils for each plot were collected. The soil samples were sieved (2 mm mesh) to remove stones, visible roots and residues and other debris, and divided into two sub-samples: one stored at −80°C for sequencing, and the other one stored at 4°C for soil physicochemical properties.

2.2. Measurements of soil chemical properties

A soil-water (deionized water) (1:2.5 w/v) suspension was shaken for 30 min prior to measuring the pH with a pH meter (Thermo Scientific Orion 3-Star Benchtop, Cambridge, United Kingdom). Soil moisture content (SMC) was measured by comparing the fresh wet weight with the dry weight after drying at 120°C for 24 h. Soil organic carbon (SOC) and the total nitrogen (TN) content were measured using an elemental analyzer (Elementar, Langenselbold, Germany). Ammonium (NH₄⁺-N) and nitrate (NO₃⁻-N) nitrogen content were measured using a continuous flow analysis system (SKALAR SAN++, Breda, the Netherlands). The total phosphorus (TP) content was measured using a spectrophotometer, and available phosphorus (AP) content was measured using the colorimetric method upon extraction with 0.5 M NaHCO₃. The total potassium (TK) content was measured by digesting the soil with concentrated hydrofluoric acid, and available potassium (AK) content was extracted by acetic acid and ammonium leaching method. The extracted TK and AK content were determined using inductively coupled plasma atomic emission spectrometry (ICP-AES-7500, Shimadzu, Japan). Soil microbial biomass C (MBC) and biomass N (MBN) content were measured with a TOC analyzer (TOC-LCPH, Shimadzu, Japan). Soil mechanical compositions (Sand, Silt, Clay) were determined according to the method of Zhang et al. (2021). Three independent replicates per sample were performed for all the soil physicochemical properties.

2.3. Soil DNA extraction and ITS rRNA sequencing

Using the MOBIO Power Soil Extraction Kit (Mo Bio Laboratories, Carlsbad, CA, United States), soil total DNA was extract from 1 g of fresh soil according to the manufacturer’s instructions. The DNA was diluted in TE buffer (DNA Elution Solution-Ultra Pure Water). The DNA quantity and quality were detected using a NanoDrop ND-1000 spectrophotometer (Thermo Scientific, United States).

Fungal ITS rRNA region was amplified using primers ITS1 (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and ITS2 (5′- GCTGCGTTCATCGATGC -3′) (Fouts et al., 2012). A 6-bp barcode sequence unique to each sample was added to the primers for distinguishing multiple samples. The PCR reaction was performed in triplicate in a 25 μL mixture containing 2.5 μL of TransStart Buffer, 2 μL of dNTPs, 1 μL of each primer (10 ng/μL), and 30 ng of template DNA. The PCR conditions were as follows: pre-denaturation at 94°C for 5 min, 30 cycles of dunaturation at 94°C for 30 s, 55°C for 30 s, 72°C for 45 s, and a final extension at 72°C for 10 min. The PCR products were inspected by 2% agarose electrophoresis, and were purified using the AxyPrep DNA purification kit. Three independent PCR replicates per sample and then three PCR samples were pool at equal amount and PE300 paired-end sequenced on the Illumina Miseq v3 platfrom (2 × 300 bp). The raw sequences were uploaded to the Sequence Read Archive (SRA) database and accession number was SUB10527794.

2.4. Bioinformatics and statistical analysis

Sequences were analyzed using QIIME (version 1.8¹) software on the Allwegene bioinformation cloud platform.² The original PE reads were quality filtering following criteria: if the mean score < 20 or the length < 200 bp, and the ambiguities sequence were removed. The forward and reverse reads merged using PEAR software (version 0.9.8). The chimeras removed using Usearch (version 7.1).³Operational taxonomic units (OTU) were generated at a similarity level of 97% using the UPARSE algorithm (Edgar, 2013). Taxonomic analysis was performed on the representative sequences of OTU, with a confidence threshold of 0.7, and the UNITE 8.2 fungi database was used for comparison.⁴

Before further analysis, the sequences were normalized according to the lowest number of sequences for a single sample. The sequences occurring <1% were classified into “others.” Alpha diversity indexes (Richness, Chao1, Shannon index, PD=Phylogenetic diversity) were calculated in QIIME1 based on OTU table. Principal coordinate analysis (PCoA) was finished using R software (v.3.2.5, R Development Core Team, 2016) “Vegan” package based on Bray–Curtis dissimilarity at OTU level. Rarefaction curve was finished using R software (v.3.2.5, R Development Core Team, 2016) “microeco” package. Redundancy analysis (RDA) was performed using R software (v.3.2.5, R Development Core Team, 2016) “microeco” package based on OTU table and soil physicochemical parameters. Fungal functional guilds were assigned by using FUNGuild v1.0 and the differences of guilds among elevations were performed by one-way analysis of variance (ANOVA), Duncan test (p < 0.05). One-way analysis of variance (ANOVA) was used to detect the difference of soil physicochemical parameters among elevations using SPSS software (version 22.0). Pearson correlation analysis among fungal community composition (both phyla and genera) and soil physicochemical properties was also performed using SPSS software (version 22.0). Permutational multivariate analysis of variance (PERMANOVA) was used to test the differences in soil fungal composition among the four elevations, using Bray-curtis distance matrices (i.e., the adonis2 function of the vegan package) (Zhong and Fu, 2022).

3. Results

3.1. Soil physicochemical properties

All soil physicochemical parameters, except SOC and sand content, were significantly different (p < 0.05) among the four elevations (Table 1). Soil pH, and the content of SMC, nitrate nitrogen and available potassium declined with increasing elevation, while other soil physicochemical properties did not show a clear tendency (Table 1). Silt ranged from 0.5% (1,910 m a.s.l.) to 1.6% (2,020 m a.s.l.), and the clay ranged from 4.3% (1,690 m a.s.l.) to 5.1% (1,910 m a.s.l.).

TABLE 1

Table 1. Soil physicochemical characteristics along an elevational gradient on the Changbai Mountains, northeastern China.

3.2. Rarefaction curve and fungal diversity

The rarefaction curve (Figure 1) tended to flatten, indicating that the sequencing number was sufficient and reasonable to cover the fungal communities. The Chao1, richness, phylogenetic diversity (PD), and Shannon index of the soil fungi were significant different (one way ANOVA, p < 0.01) among elevations along the elevational gradient (Table 2). The alpha diversity indices including Chao1, richness, PD_whole_tree, and Shannon-Wiener decreased from 1,690 m to 1800 m a.s.l., then increased to 1910 m a.s.l., followed by decreasing tendency up to 2,020 m a.s.l.

FIGURE 1

Figure 1. Rarefaction curve of fungal sequences in soils along an elevational gradient on the Changbai Mountain, northeastern China. Note A (1–3) is 1,690 m a.s.l.; B (1–3) is 1,800 m a.s.l.; C (1–3) is 1910 m a.s.l.; D (1–3) is 2020 m a.s.l.

TABLE 2

Table 2. Fungal alpha diversity along an elevational gradient on the Changbai Mountain, northeastern China¹.

Pearson correlation analysis showed that the soil fungal alpha diversity (richness and PD index) were significantly positively correlated with the content of soil NO₃⁻-N, but significantly negatively correlated with the content of soil NH₄⁺-N and TN (Table 3). PD index of soil fungi was significantly positively correlated with soil pH, while the Shannon-Wiener index of soil fungi was significantly negatively correlated with the content of soil TN (Table 3).

TABLE 3

Table 3. Pearson’s rank correlation coefficients between fungal alpha-diversity and soil physicochemical characteristics.

3.3. Composition of the soil fungal community

The ordination showed that the fungal communities were clearly separated by elevation (PCo1, 44.63% and PCo2, 25.36%) (Figure 2). The soil fungal beta diversity differed significantly among elevations (Figure 2, Anosim R = 0.97, p < 0.01). Moreover, there were significant differences in species composition between any two out of the four elevations (Supplementary Table S2, p < 0.05).

FIGURE 2

Figure 2. Principal coordinate analysis (PCoA) of fungal communities along an elevational gradient on Changbai Mountain, northeastern China. The dominant plant species at the sites were D. angustifolia. Beta-Diversity index was calculated at the OTU level (97%).

The soil fungal community composition (both at the phyla and genera level) was significantly different among the four elevations (Figure 3 and Supplementary Tables S1, S3). All the obtained sequences belonged to 8 phyla. The prevailing phyla was Basidiomycota (53% relative abundance), Ascomycota (26%), and Mortierellomycota (19%) across all the soil samples (Figure 3A). The relative abundance of these dominant fungal phyla changed remarkably with elevation (Supplementary Table S1), showing that the relative abundance of Basidiomycota had an increased tendency but Ascomycota showed a decreased tendency with increasing elevation (Figure 3A and Supplementary Table S1).

FIGURE 3

Figure 3. Relative abundance of the dominant fungal phyla (A) and genera (B) in soils along an elevational gradient on the Changbai Mountain, northeastern China.

The dominant genera with a relative abundance of >1% were Mortierella (18.7% relative abundance), Russula (13.3%), Inocybe (7.7%), Archaeorhizomyces (4.7%), Laccaria (4.2%), Tricholoma (3.2%), Leotia (2.5%), Cortinarius (2.5%), Clavulina (2.4%), Entoloma (2.1%), Elaphomyces (2.0%), Piloderma (1.6%), Tomentella (1.3%), Solicoccozyma (1.1%), Gymnomyces (1.1%), Sebacina (1.0%) (Figure 3B). Russula was most abundant at the highest site (2,020 m a.s.l.), while Inocybe was most abundant at the lowest elevational site (Figure 3B). Similar to the phyla, the relative abundance of these most dominant fungal genera also changed remarkably with increasing elevation (Supplementary Table S3).

3.4. Relationships between fungal community and soil properties

RDA revealed that soil properties (i.e., NH₄⁺-N, NO₃⁻-N, TN, SOC, MC, pH, AK, TK, and TP) were the key environmental factors that shaped the soil fungal community (Figure 4). The first two axes of the RDA accounted for 37.57% of the total variance. As shown in Figure 4, MC (p < 0.05), pH (p < 0.05), NH₄⁺-N (p < 0.05), NO₃⁻-N (p < 0.05), TN (p < 0.05) significantly influenced the fungal community. The soil fungal community structure at 2,020 m a.s.l. was significantly positively correlated with NH₄-N and TN, while that at 1,690 m a.s.l. and 1,800 m a.s.l. was significantly positively correlated with soil moisture content, pH, AK and NO₃^—N (Figure 4).

FIGURE 4

Figure 4. Redundancy analysis (RDA) of soil fungal community structures and environmental characteristics (arrows) at different elevations.

At the phyla level, the abundance of phyla was closely correlated with some certain, phyla-specific soil physicochemical factors (Table 4). For the two dominant phyla, for example, the abundance of Ascomycota was significantly negatively correlated with the content of soil TN, NH₄⁺-N, and AP (Table 4), while the abundance of Basidiomycota was significantly negatively correlated with soil TK (Table 4).

TABLE 4

Table 4. Pearson’s rank correlations between the relative abundances of dominant fungal taxa and soil physicochemical variables.

At the genus level, the abundance of genus was also closely correlated with some certain, genus-specific soil physicochemical factors (Table 4). For instance, the abundance of Ilyonectria, Inocybe and Archaeorhizomyces was significantly positively correlated with SMC, while the abundance of Russula was significantly negatively correlated with the content of SMC (Table 4).

3.5. Fungal functional guilds

The dominated functional groups of fungi were ectomycorrhizal, Endophyte, Undefined Saprotroph, Plant Pathogen, Animal Pathogen, Wood Saprotrophc, Lichenized, Ericoid Mycorrhizal and Arbuscular Mycorrhizal (Supplementary Table S4). Except Ericoid Mycorrhizal, Lichenized and Plant Pathogen, other soil fungal functional groups were significantly different among elevations along the elevational gradient (Supplementary Table S4). The absolute abundance of Arbuscular Mycorrhizal and Lichenized was highest at 1,690 m a.s.l., and the absolute abundance of Animal Pathogen, Ericoid Mycorrhizal and Wood Saprotroph was highest at 1,810 m a.s.l. while the absolute abundance of Ectomycorrhizal was highest at 2020 m a.s.l. (Supplementary Table S4).

4. Discussion

4.1. Changes in soil physicochemical properties with elevation

In this study, we found significant differences in soil physicochemical properties among elevations along the elevational gradient from 1,690 m to 2020 m a.s.l. on the Changbai Mountain (Table 1), this is consistent with the results of Zong et al. (2014) who investigated the soil physicochemical properties in D. angustifolia population along elevational gradients on the same mountain. The climate changes dramatically with increasing elevation on mountains, which causes markedly changes in soil environment such as biogeochemical cycling and soil nutrients (Zong et al., 2014). In our study, for example, the soil water content and pH value were significantly higher at 1,690 m than at 2,020 m a.s.l., which indicates that soil moisture holding capacity decreases with the invasion of D. angustifolia. We found a significant lower value of soil organic carbon at higher elevation (2,020 m) and at lower elevation (1,800 m), which is consistent with results of previous studies (Meng et al., 2018; Luo et al., 2020).

4.2. Changes in soil fungal alpha diversity with elevation

Soil fungi play an important role in biogeochemical cycling and ecological process (Buee et al., 2009; Russo et al., 2012), but previous studies have focused on bacterial diversity and composition, with only a few focusing on fungi in mountain ecosystems. Therefore, we have only limited information about changes in fungi with elevation (Djukic et al., 2010; Adamczyk et al., 2019; Zhou et al., 2021). Our results showed that the alpha diversity (Shannon index, Chao1 index, Richness and PD index) of soil fungi changed significantly with elevation, showing that the alpha diversity decreased with increasing elevation (Table 3). However, this pattern does not seem to be widespread. The distribution pattern of soil fungal diversity in the literature showed declining, humped, U-shaped, or no change with increasing elevation. Yang et al. (2017) reported that soil fungal diversity decreased monotonically from 700 m to 2,600 m a.s.l. across various ecosystems on the Changbai Mountain. Ping et al. (2017) reported that soil fungal diversity in Pinus koraiensis forest showed a hollow curve’s pattern along an elevational gradient from 699 m to 1,044 m a.s.l. on the Changbai Mountain.

Our study indicated that the alpha diversity of soil fungi showed a decreasing trend with increasing elevation. Shen et al. (2014) showed that the Chao1 index of soil fungi on the Changbai Mountain was not correlated with elevation, but it had a strong correlation with soil pH. Ni et al. (2018) found that soil fungal Chao1 index increased with increasing elevation, and soil C/N was the most important environmental factor determining the soil fungal Chao1 index along an elevational gradient from 2,000 m to 2,500 m a.s.l. on the Changbai Mountain. However, Zhang et al. (2022) found that both the Chao1 and Shannon index of soil fungi decreased significantly from 2,785 m a.s.l. to 4,578 m, and the climate and soil properties had opposite effects on them in Tibetan Plateau.

In our present study, soil TN, NH₄⁺-N, and NO₃⁻-N were significantly correlated with fungal richness (r = −0.53, p < 0.05; r = 0.57, p < 0.05; r = −0.59, p < 0.01) along the elevational gradient. As the soil N content decreased with increasing elevation (Zong et al., 2014), resulting thus in a reduction in soil fungal diversity found in the present study. Similar studies have also demonstrated that soil nutrients are an important environmental factor affecting the distribution pattern of soil fungal diversity (Newsham et al., 2016; Yang et al., 2017).

4.3. Fungal compositions

In our present study, the dominant phyla are Basidiomycota and Ascomycota, which is consistent with findings of previous studies on fungal composition along elevational gradients on Changbai Mountain. However, the elevational patterns of Basidiomycota and Ascomycota differed among studies. Our study found that Basidiomycota showed an increasing trend, while Ascomycota showed a decreasing trend with elevation from 1,690 m to 2,020 m a.s.l. However, Shen et al. (2014) and Ni et al. (2018) found that Basidiomycota showed an overall increasing trend, while Ascomycota’s did not change or tended to decrease with elevation on Changbai Montain. Ping et al. (2017) found that the abundance of Basidiomycota first decreased from 699 m to 937 m a.s.l., then increased from 937 m to 1,177 m a.s.l., while the abundance of Ascomycota increased from 699 m to 937 m a.s.l., and then decreased from 937 m to 1,177 m a.s.l. Yu et al. (2019) and Zhang et al. (2021) found that Ascomycota and Mortierellomycota are the predominant fungi in Tibetan grassland communities. This difference in dominant phyla may be closely related to the fact that the vegetation (or plant ecosystem) and soil physicochemical properties differ significantly between the two regions (Jin et al., 2018; Zong and Fu, 2021).

Soil fungal composition is affected by multiple factors such as vegetation composition, soil physicochemical properties, and microclimate along elevational gradients (Ni et al., 2018). A previous study carried ou on the Changbai Mountain found that the aboveground vegetation composition was closely related to the fungal composition, and dominant plant species significantly affected the fungal composition (Ni et al., 2018). In our present study, the dominant Basidiomycota and Ascomycota were mainly affected by soil physicochemical properties because the vegetation population (i.e., D. angustifolia population) did not change along with the elevational gradient. The fungal Basidiomycota and Ascomycota are involved in the soil organic metabolism (Luo et al., 2021) and thus their abundances were significantly determined by soil organic matter content as a result of decomposition of plant residues (Li et al., 2019).

Our results also showed that the functional guilds of soil fungi at high elevation were mainly ectomycorrhizal fungi and plant pathogens. Ni et al. (2018) also reported that the dominant functional fungi were ectomycorrhizal fungi in the alpine tundra on the Changbai Mountain. Similarly, Timling et al. (2014) reported that the dominant functional guilds in the arctic were the ectomycorrhizal fungi. However, in this study, the alpine tundra population of D. angustifolia is an invasive plant, and the effect of the original tundra species on the soil has not completely replaced by the invasive plant, so that the soil fungal functions were still ectomycorrhizal fungi and plant pathogens.

5. Conclusion

The present study revealed that a small elevational difference on mountains may lead to marked difference in soil physicochemical properties, fungal diversity, and community composition. Our results showed that fungal alpha diversity was higher at lower elevations, which may be a result of higher soil nutrient levels at lower elevations. The phyla of Ascomycota and Basidiomycota, as well as the genera of Mortierella and Russula dominated the soil fungal communities across the entire elevational gradient. Generally, Basidiomycota increased but Ascomycota decreased with increasing elevation. The changes in soil pH and nutrients were the most important soil environment factors leading to changes in soil fungal beta diversity. Our results highlight the different patterns of fungal communities across elevational gradients, and further elucidate the variation in fungal community composition and ecological functions in temperate mountain ecosystems.

Data availability statement

The data presented in the study are deposited in the NCBI Sequence Read Archive repository, accession number SUB10527794.

Author contributions

XS and ML designed and performed the experiment and prepared this manuscript. BF and M-HL revised this manuscript and language editing. GD and LY helped to do the experiment and finish the bioinformatic analysis. All coauthors contributed to manuscript editing, read, and agreed to the published version of the manuscript.

Funding

This work was funded by the Natural Sciences Foundation of Heilongjiang Province (LH2020C088), Heilongjiang Province Postdoctoral Research Start-up Fund Project (LBH-Q21167), Outstanding Youth Foundation of Heilongjiang University (JCL202006), the China Scholarship Council Visiting Scholar Program (201908230401), and the Basic Scientific Research of Provincial Higher Education Institutions in Heilongjiang Province of 2022. Heilongjiang Provincial Ecological Environmental Protection Research Project (HST2022ST008) and the central government guides local science and technology development special projects (ZY20B15). Open access funding by WSL - Swiss Federal Institute For Forest, Snow And Landscape Research.

Conflict of interest

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

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

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

Footnotes

1. ^http://qiime.org

2. ^http://218.2.224.234

3. ^http://drive5.com/aparse

4. ^https://unite.ut.ee/

References

Adamczyk, M., Hagedorn, F., Wipf, S., Donhauser, J., and Frey, B. (2019). The soil microbiome of gloria mountain summits in the swiss alps. Front. Microbiol. 10:1080. doi: 10.3389/fmicb.2019.01080

PubMed Abstract | CrossRef Full Text | Google Scholar

Buee, M., Reich, M., Murat, C., Morin, E., Nilsson, R. H., Uroz, S., et al. (2009). 454 Pyrosequencing analyses of forest soils reveal an unexpectedly high fungal diversity. New Phytol. 184:449–456.

Google Scholar

Deng, J. J., Zhou, W. M., Dai, L. M., Yuan, Q., Zhou, L., Qi, L., et al. (2023). The effects of shrub removal on soil microbial communities in primary forest, secondary forest and plantation forest on Changbai Mountain. Microb. Ecol. 85:642–658. doi: 10.1007/s00248-021-01943-0

CrossRef Full Text | Google Scholar

Djukic, I., Zehetner, F., Mentler, A., and Gerzabek, M. H. (2010). Microbial community composition and activity in different alpine vegetation zones. Soil Biol. Biochem. 42, 155–161. doi: 10.1016/j.soilbio.2009.10.006

CrossRef Full Text | Google Scholar

Donhauser, J., and Frey, B. (2018). Alpine soil microbial ecology in a changing world. FEMS Microbiol. Ecol. 94:fiy099. doi: 10.1093/femsec/fiy099

CrossRef Full Text | Google Scholar

Donhauser, J., Niklaus, P. A., Rousk, J., Larose, C., and Frey, B. (2020). Temperatures beyond the community optimum promote the dominance of heat-adapted, fast growing and stress resistant bacteria in alpine soils. Soil Biol. Biochem. 148:107873. doi: 10.1016/j.soilbio.2020.107873

CrossRef Full Text | Google Scholar

Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604

PubMed Abstract | CrossRef Full Text | Google Scholar

Fouts, D. E., Sebastian, S., Janaki, P., Manolito, T., Waterman, R. C., and Macneil, M. D. (2012). Next generation sequencing to define prokaryotic and fungal diversity in the bovine rumen. PLoS One 7:e48289. doi: 10.1371/journal.pone.0048289

PubMed Abstract | CrossRef Full Text | Google Scholar

Frey, B., Rime, T., Phillips, M., Stierli, B., Hajdas, I., Widmer, F., et al. (2016). Microbial diversity in European alpine permafrost and active layers. FEMS Microbiol. Ecol. 92, 1–17. doi: 10.1093/femsec/fiw018

CrossRef Full Text | Google Scholar

Han, F. S., Yu, C. Q., and Fu, G. (2022). Warming alters elevation distributions of soil bacterial and fungal communities in alpine grasslands. Glob. Ecol. Conserv. 39:e02306. doi: 10.1016/j.gecco.2022.e02306

CrossRef Full Text | Google Scholar

Hanif, M. A., Guo, Z., Moniruzzaman, M., He, D., Yu, Q., Rao, X., et al. (2019). Plant taxonomic diversity better explains soil fungal and bacterial diversity than functional diversity in restored forest ecosystems. Plan. Theory 8:479. doi: 10.3390/plants8110479

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, Y. H., Zhang, Y. J., Xu, J. W., Tao, Y., He, H. S., Guo, M., et al. (2018). Comparative assessment of tundra vegetation changes between north and southwest slopes of Changbai Mountains, China, in response to global warming. Chin. Geogr. Sci. 28, 665–679. doi: 10.1007/s11769-018-0978-y

CrossRef Full Text | Google Scholar

Li, J., Delgado-Baquerizo, M., Wang, J. T., Hu, H. W., and Singh, B. K. (2019). Fungal richness contributes to multifunctionality in boreal forest soil. Soil Biol. Biochem. 136:107526. doi: 10.1016/j.soilbio.2019.107526

CrossRef Full Text | Google Scholar

Li, L., Xing, M., Lv, J., Wang, X., and Chen, X. (2017). Response of rhizosphere soil microbial to Deyeuxia angustifolia encroaching in two different vegetation communities in alpine tundra. Sci. Rep. 7, 1–13. doi: 10.1038/srep43150

CrossRef Full Text | Google Scholar

Lin, Y. T., Jia, Z., Wang, D., and Chiu, C. Y. (2017). Effects of temperature on the composition and diversity of bacterial communities in bamboo soils at different elevations. Biogeosciences 14, 4879–4889. doi: 10.5194/bg-14-4879-2017

CrossRef Full Text | Google Scholar

Luo, Z. M., Liu, J. X., Bao, J. B., Wang, X., Wu, Q., and Chai, B. F. (2020). Elevational distribution patterns and assembly mechanisms of soil fungal community in Mount Wutai, Shanxi, China. Acta Ecol. Sin. 40, 7009–7017. doi: 10.5846/stxb201905171023

CrossRef Full Text | Google Scholar

Luo, X., Liu, K., Shen, Y., Yao, G., Yang, W., Mortimer, P. E., et al. (2021). Fungal community composition and diversity vary with soil horizons in a subtropical forest. Front. Microbio. 12:650440.

Google Scholar

Meng, M. J., Guo, X. P., and Zhang, J. C. (2018). Effects of altitude on soil microbial community in Fengyang Mountain coniferous and broad-leaved forest. Acta Ecol. Sin. 38, 7057–7065. doi: 10.5846/stxb201708211503

CrossRef Full Text | Google Scholar

Newsham, K. K., Hopkins, D. W., Carvalhais, L. C., Fretwell, P. T., Rushton, S. P., O’Donnell, A. G., et al. (2016). Relationship between soil fungal diversity and temperature in the maritime Antarctic. Nat. Clim. Change. 6, 182–186.

Google Scholar

Ni, Y., Yang, T., Zhang, K., Shen, C., and Chu, H. (2018). Fungal communities along a small-scale elevational gradient in an alpine tundra are determined by soil carbon nitrogen ratios. Front. Microbiol. 9:1815. doi: 10.3389/fmicb.2018.01815

PubMed Abstract | CrossRef Full Text | Google Scholar

Ping, Y., Han, D. X., Wang, N. N., Hu, Y., Mu, L. Q., and Feng, F. J. (2017). Vertical zonation of soil fungal community structure in a Korean pine forest on Changbai Mountain, China. World J. Microbiol. Biotechnol. 33, 1–10. doi: 10.1007/s11274-016-2133-1

CrossRef Full Text | Google Scholar

Prasad, S., Manasa, P., Buddhi, S., Tirunagari, P., Begum, Z., Rajan, S., et al. (2014). Diversity and bioprospective potential (cold-active enzymes) of cultivable marine bacteria from the subarctic glacial fjord, Kongsfjorden. Curr. Microbiol. 68, 233–238. doi: 10.1007/s00284-013-0467-6

PubMed Abstract | CrossRef Full Text | Google Scholar

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

Google Scholar

Ren, C. J., Liu, W. C., Zhao, F. Z., Zhong, Z. K., Deng, J., Han, X. H., et al. (2019). Soil bacterial and fungal diversity and compositions respond differently to forest development. Catena 181:104071. doi: 10.1016/j.catena.2019.104071

CrossRef Full Text | Google Scholar

Rime, T., Hartmann, M., and Frey, B. (2016). Potential sources of microbial colonizers in an initial soil ecosystem after retreat of an alpine glacier. ISME J. 10, 1625–1641. doi: 10.1038/ismej.2015.238

PubMed Abstract | CrossRef Full Text | Google Scholar

Russo, S. E., Legge, R., Weber, K. A., Brodie, E. L., Goldfarb, K. C., Benson, A. K., et al. (2012). Bacterial community structure of contrasting soils underlying Bornean rain forests: Inferences from microarray and next-generation sequencing methods. Soil Biol. Biochem. 55, 48–59.

Google Scholar

Sahay, H., Babu, B. K., Singh, S., Kaushik, R., Saxena, A. K., and Arora, D. K. (2013). Cold-active hydrolases producing bacteria from two different sub-glacial Himalayan lakes. J. Basic Microbiol. 53, 703–714. doi: 10.1002/jobm.201200126

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, C. C., Liang, W., Shi, Y., Lin, X., Zhang, H., Wu, X., et al. (2014). Contrasting elevational diversity patterns between eukaryotic soil microbes and plants. Ecology 95, 3190–3202. doi: 10.1890/14-0310.1

CrossRef Full Text | Google Scholar

Shen, C. C., Xiong, J., Zhang, H., Feng, Y., Lin, X., Li, X., et al. (2013). Soil pH drives the spatial distribution of bacterial communities along elevation on Changbai Mountain. Soil Biol. Biochem. 57, 204–211. doi: 10.1016/j.soilbio.2012.07.013

CrossRef Full Text | Google Scholar

Shi, L. L., Mortimer, P. E., Ferry Slik, J. W., Zou, X. M., Xu, J., Feng, W. T., et al. (2014). Variationin forest soil fungal diversity along a latitudinal gradient. Fungal Divers. 64, 305–315. doi: 10.1007/s13225-013-0270-5

CrossRef Full Text | Google Scholar

Shivaji, S., Begum, Z., Rao, S. S. S. N., Reddy, P. V. V. V., Manasa, P., Sailaja, B., et al. (2013). Antarctic ice core samples: culturable bacterial diversity. Res. Microbiol. 164, 70–82. doi: 10.1016/j.resmic.2012.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Sui, X., Zhang, R. T., Frey, B., Yang, L. B., Liu, Y. N., Ni, H. W., et al. (2021). Soil physicochemical properties drive the variation in soil microbial communities along a forest successional series in a degraded wetland in northeastern China. Ecol. Evol. 11, 2194–2208. doi: 10.1002/ece3.7184

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, L. N., Li, A. X., and Shao, G. F. (2011). Landscape-level forest ecosystem conservation on Changbai Mountain, China and North Korea (DPRK). Mt. Res. Dev. 31, 169–175. doi: 10.1659/MRD-JOURNAL-D-10-00120.1

CrossRef Full Text | Google Scholar

Tedersoo, L., Bahram, M., Cajthaml, T., Polme, S., Hiiesalu, I., Anslan, S., et al. (2016). Tree diversity and species identity effects on soil fungi, protists and animals are context dependent. ISME J. 10, 346–362. doi: 10.1038/ismej.2015.116

PubMed Abstract | CrossRef Full Text | Google Scholar

Timling, I., Walker, D. A., Nusbaum, C., Lennon, N. J., and Taylor, D. L. (2014). Rich and cold: diversity, distribution and drivers of fungal communities in patterned‐ground ecosystems of the North American Arctic. Mol. Ecol. 23, 3258–3272.

Google Scholar

Wang, N. N., Yang, X., Li, S. L., Sui, X., Han, S. J., and Feng, F. J. (2013). Effects of precipitation variation on the distribution pattern of soil fungal diversity in broad-leaved Korean pine mixed forest. Chin. J. Appl. Ecol. 24, 1985–1990.

Google Scholar

Weng, X. H., Sui, X., Liu, Y. N., Yang, L. B., and Zhang, R. T. (2022). Effect of nitrogen addition on the carbon metabolism of soil microorganisms in a Calamagrostis angustifolia wetland of the Sanjiang Plain, northeastern China. Ann. Microbiol. 72, 1–14. doi: 10.1186/s13213-022-01674-8

CrossRef Full Text | Google Scholar

Xue, D., and Tisdell, C. (2001). Valuing ecological functions of biodiversity in Changbaishan Mountain Biosphere reserve in Northeast China. Biodivers. Conserv. 10, 467–481. doi: 10.1023/A:1016630825913

CrossRef Full Text | Google Scholar

Yang, W., Li, S., Wang, X., Liu, F., Li, X., and Zhu, X. (2021). Soil properties and geography shape arbuscular mycorrhizal fungal communities in black land of China. Appl. Soil Ecol. 167:104109. doi: 10.1016/j.apsoil.2021.104109

CrossRef Full Text | Google Scholar

Yang, H., Lü, G., Jiang, H., Shi, D. N., and Liu, Z. (2017). Diversity and distribution of soil micro-fungi along an elevation gradient on the north slope of Changbai Mountain. J. For. Res. 28, 831–839. doi: 10.1007/s11676-016-0344-9

CrossRef Full Text | Google Scholar

Yu, C. Q., Han, F. S., and Fu, G. (2019). Effects of 7 years experimental warming on soil bacterial and fungal community structure in the Northern Tibet alpine meadow at three elevations. Sci. Total Environ. 655, 814–822. doi: 10.1016/j.scitotenv.2018.11.309

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, T., Liu, Y. F., Sui, X., and Song, F. Q. (2021). Land use patterns: effects on soil physical and chemical properties and enzyme activities in the Western Heilongjiang. J. Agric. 11, 33–41. doi: 10.11923/j.issn.2095-4050.cjas20190700110

CrossRef Full Text | Google Scholar

Zhang, G. Y., Shen, Z. X., and Fu, G. (2022). Geo-distribution patterns of soil fungal Community of Pennisetum flaccidum in Tibet. J. Fungi 8:1230. doi: 10.3390/jof8111230

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhong, Z., and Fu, G. (2022). Response of soil fungal species, phylogenetic and functional diversity to diurnal asymmetric warming in an alpine agricultural ecosystem. Agric. Ecosyst. Environ. 335:107993. doi: 10.1016/j.agee.2022.107993

CrossRef Full Text | Google Scholar

Zhou, Y. (2006). Investigation of the alpine tundra seed plants in Changbai Mountain, Guihaia. 26, 152–156.

Google Scholar

Zhou, Y., Jia, X., Han, L., Liu, Z., Kang, S., and Zhao, Y. (2021). Fungal community diversity in soils along an elevation gradient in a Quercus aliena var. acuteserrata forest in Qinling Mountains, China. Appl. Soil Ecol. 167:104104. doi: 10.1016/j.apsoil.2021.104104

CrossRef Full Text | Google Scholar

Zong, N., and Fu, G. (2021). Variations in species and function diversity of soil fungal community along a desertification gradient in an alpine steppe. Ecol. Indic. 131:108197. doi: 10.1016/j.ecolind.2021.108197

CrossRef Full Text | Google Scholar

Zong, S. W., Xu, J. W., and Wu, Z. F. (2013). Investigation and mechanism analysis on the invasion of Deyeuxia. Angustifolia to tundra zone in western slope of Changbai Mountain. J. Mount. Sci. 31, 448–455.

Google Scholar

Zong, S. W., Xu, J. W., Wu, Z. F., Qiao, L. L., Wang, D. D., Meng, X. J., et al. (2014). Analysis of the process and impacts of Deyeuxia angustifolia invasion on the Alpine tundra, Changbai Mountain. Acta Ecol. Sin. 34, 6837–6846. doi: 10.5846/stxb201302270319

CrossRef Full Text | Google Scholar