Soil organic matter and water content affect the community characteristics of arbuscular mycorrhizal fungi in Helan mountain, an arid desert grassland area in China

Introduction: Arbuscular mycorrhizal fungi (AMF) are vital in terrestrial ecosystems. However, the community structure characteristics and influencing factors of AMF in the forest ecosystems of arid desert grassland areas require further investigation.

Methods: Therefore, we employed high-throughput sequencing technology to analyze the soil AMF community characteristics at different elevations in the Helan mountains.

Results: The results revealed that significant differences (P < 0.05) were observed in the soil physicochemical properties among different elevations, and these properties exhibited distinct trends with increasing elevation. Through high-throughput sequencing, we identified 986 operational taxonomic units (OTUs) belonging to 1 phylum, 4 classes, 6 orders, 12 families, 14 genera, and 114 species. The dominant genus was Glomus. Furthermore, significant differences (P < 0.05) were observed in the α-diversity of the soil AMF community across different elevations. Person correlation analysis, redundancy analysis (RDA), and Monte Carlo tests demonstrated significant correlations between the diversity and abundance of AMF communities with soil organic matter (OM) (P < 0.01) and soil water content (WC) (P < 0.05).

Discussion: This study provides insights into the structural characteristics of soil AMF communities at various altitudes on the eastern slope of Helan mountain and their relationships with soil physicochemical properties. The findings contribute to our understanding of the distribution pattern of soil AMF and its associations with environmental factors in the Helan mountains, as well as the stability of forest ecosystems in arid desert grassland areas.

1 Introduction

Mycorrhiza is a symbiotic association formed between the mycelium of fungi and the root system of higher plants. Symbiotic fungi expand the absorptive capacity of plant roots by obtaining nutrients for their growth and development from the plant. Simultaneously, they enhance the plant’s ability to acquire water and nutrients and facilitate interactions between the plant and other soil microorganisms (Tedersoo et al., 2020). Mycorrhizal fungi are important components of natural ecosystems and have a significant impact on plant community composition and ecosystem function (Koide and Mosse, 2004). About 320,000–340,000 plants can form mycorrhizae, with the most widespread distribution of arbuscular mycorrhizal fungi (AMF), which account for about 72% of mycorrhizal fungi (Genre et al., 2020). Moreover, in addition to nutrient uptake, AMF contributes to enhancing plant resistance to environmental and biological stresses, thereby improving plant disease and stress tolerance (Hashem et al., 2021; Dowarah et al., 2022). Several studies comparing different ecosystems have demonstrated a greater reliance on AMF for plant survival in semi-arid habitats (Li and Zhao, 2005). Other studies have shown that the diversity and abundance of AMF communities change along environmental gradients, mainly due to the influence of different external environmental factors such as altitude, soil physicochemical properties, and climate (Bonfim et al., 2016; Prudent et al., 2020). Particularly, factors such as soil moisture, organic matter (OM), and the content of nitrogen and phosphorus elements have a significant impact on the diversity and abundance of AMF communities (Gomoryova et al., 2013; Egan et al., 2017).

Mountain ecosystems serve as natural experimental areas for studying the biodiversity of organisms along environmental gradients, including continuous temperature changes (Barcenas-Moreno et al., 2009), soil physicochemical properties, and vegetation composition (Wang et al., 2022), providing valuable research opportunities. Previous research results have shown that the abundance and diversity of microorganisms exhibit different trends with changes in elevation gradient. Fungi in Andean Yungas forests (Wang et al., 2015) and bacteria near Gunnison County, Colorado (Bryant et al., 2008) exhibited a decrease in both diversity and abundance with increasing elevation. The relative abundance of Glomus in the Guangxi hills of China increased with elevation (Yu et al., 2023). In the inter-root soils of the Taibai Mountains in China, microbial α-diversity initially increased and then decreased with elevation (Ren et al., 2021). In the arid valley of Panzhihua, there was a trend of decreasing and then increasing microbial community richness with elevation (Zhang R. J. et al., 2022), and it has also been demonstrated that the diversity of the AMF community is not correlated with elevation (Shen et al., 2014). These observations can be attributed to the fact that environmental conditions vary across regions and elevations, creating complex circumstances for microbial survival and distribution. Additionally, slope orientation is an important influencing factor, as there are significant differences in vegetation types, soil physicochemical properties, light exposure, and temperature between shady and sunny slopes, which may have a significant impact on microbial abundance and diversity (Guo et al., 2022).

As an important ecological barrier in northwest China, the Helan mountains exhibit significant vertical vegetation zonation and a complete vegetation spectrum. Along the altitude gradient, vegetation types transition from desert grasslands to shrublands, sparse forests, subalpine coniferous forests, and finally to alpine shrub meadows (Jiang et al., 2000). The diversity and abundance of AMF are closely related to these different vegetation types (Uhlmann et al., 2004; Dassen et al., 2017), and vegetation diversity influences the community structure of AMF (Prommer et al., 2020). At the same time, AMF plays an important role in maintaining the stability of plant communities. There are relatively limited studies on the diversity and distribution of AMF communities and the main factors affecting them in the Helan mountains. In this study, we focused on the eastern slope of the Helan mountains as the study area and used high-throughput sequencing technology to investigate the characteristics of AMF communities in the rhizosphere soil of plants at different altitudes on the sunny slope of the eastern slope. We also analyzed the correlation between various soil physicochemical indicators and AMF community characteristics. This study aimed to reveal the distribution patterns of AMF communities and their influencing factors, offering valuable reference for the maintenance of mountainous ecosystem stability.

2 Materials and methods

2.1 Overview of the sampling area

Helan mountain is situated between the Alashan Plateau and the Yinchuan Plain, spanning coordinates 38°27′-39°30′N and 105°41′-106°41′E. The highest altitude of the sampling point reaches 2,580 m. The Helan mountains have a typical continental climate, with an average annual temperature of −1 to 8°C, an average annual precipitation of 200–400 mm, and a large average annual evaporation of more than 2,000 mm. The sampling area is situated on the sunny slope of the eastern side of the Helan mountains. Common plant species in this region include Roegneria alashanica, Aristida adscensionis, Stipa breviflora, Prunus mongolica, Ulmus glaucescens, Picea crassifolia, and so on. The soil is mostly sandy loam.

2.2 Experimental design

The sample plots were established at altitudes ranging from 1,349 to 2,580 m, covering typical plant communities on the eastern slope of the Helan mountains, including grassland (GR), shrubs (SH), U. glaucescens forests (UL), Pinus tabuliformis forests (PF), coniferous and broad-leaved mixed forests (CB), P. crassifolia forests (PI), and subalpine meadow (SU), as shown in Table 1. Within each sampling area, three 20 m × 20 m quadrats were established. The five-point sampling method was used to collect plant inter-root soil from 0 to 20 cm depth with a 4 cm diameter soil sampler, which was mixed and sealed in a sterilized bag, transferred to the laboratory in a refrigerated box, and a part of the soil was air-dried and stored at room temperature, while the other part was stored in a refrigerator at −80°C to be preserved and used for soil microbial sequencing (Schleuss and Muller, 2001).

TABLE 1

Table 1. Basic information of sample plots.

2.3 Experimental measurement methods

2.3.1 Determination of soil physical and chemical properties

Soil pH was determined using the potentiometric method, with a soil-to-liquid ratio of 1:2.5. Soil water content (WC) was determined using the drying method. Soil-available phosphorus (AP) was measured using a 0.5 mol⋅ L^–1 NaHCO₃ solution. Soil total phosphorus (TP) was analyzed using the perchloric acid-sulfuric acid solubilization method, followed by the molybdenum antimony antimicrobial colorimetric method. Soil OM was quantified using the potassium dichromate oxidation-external heating method. Total nitrogen (TN) in the soil was determined using the Kjeldahl method. Soil available nitrogen (AN) was measured using the alkaline dissolution-diffusion method. Soil-available calcium (AK) was analyzed using the ammonium acetate-flame photometer method (Lv and Li, 2010).

2.3.2 Determination of soil AMF community diversity

Arbuscular mycorrhizal fungi community diversity was determined by Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China). Soil AMF community DNA was extracted by MP kit using AMV4-5NF/AMDGR primers with primer sequences of 5′-AAGCTCGTAGTTGAATTTCG-3′, 5′-CCCAACTATCCCTATTAATCAT-3′, and PCR for 18S rRNA gene amplification. The TransGen AP221-02 kit containing TransStart Fastpfu DNA polymerase was used for the PCR reaction. (1) PCR products from the same sample were mixed and detected by 2% agarose gel electrophoresis, (2) PCR products were recovered by cutting the gel using the AxyPrepDNA Gel Recovery Kit, (3) Tris_HCl elution; 2% agarose electrophoresis was detected. (4) Sodium hydroxide denaturation to generate single-stranded DNA fragments (Kowalchuk et al., 2002). DNA bipartite sequencing was performed on the Illumina MiSeq PE300 platform (Illumina, San Diego, CA, USA). Using UPARSE (Edgar, 2013) software (version 11),¹ OTUs were classified based on a 97% sequence similarity threshold. The representative sequences of OTUs with 97% similarity were taxonomically analyzed using the RDP classifier Bayesian algorithm. Finally, the taxonomic information was obtained from the maarjam20220506/AM species classification database.

2.4 Data analysis

The data were organized using Excel 2021, and SPSS 26 was utilized for data processing. To analyze the differences in soil physicochemical properties among vegetation zones (P < 0.05), one-way ANOVA and Tukey’s honest Significant Difference test were employed. Pearson’s correlation analysis (PCA) with a two-tailed significance test was used to explore the correlations between different variables. Non-metric multidimensional scaling analysis (NMDS) and analysis of similarity test (ANOSIM) were selected to examine the differences among AMF communities. Bar scale plots were generated using Origin 2022. The computation platform provided by Shanghai Meiji Biomedical Technology Co. was used for the analysis and mapping of soil AMF community diversity and its correlation with environmental factors.

Soil AMF community α-diversity analysis was conducted using the Shannon-Wiener diversity index (H), Simpson dominance index (D), Chao1 index (S_obs), and ACE index (S_ACE).

① Shannon-Wiener index (H):

$H=-\sum _{i=1}^{S_{obs}}\frac{n_i}{N}ln\frac{n_i}{N}$

② Simpson dominance index (D):

$\text{D}=-\frac{{\sum }_{\text{i}=1}^{{\text{S}}_{\text{obs}}}{\text{n}}_{\text{i}}({\text{n}}_{\text{i}}-1)}{\mathrm{N}(\mathrm{N}-1)}$

Where: S_obs is the number of OTUs actually observed; n_i is the number of sequences contained in the i-th OTU and N is the number of all sequences.

③ Chao1 index (S_obs):

$S_{chao1}=S_{obs}+\frac{F_1(F_1-1)}{2(F_2+1)}$

Where: S_chao1 is the estimated number of OTUs; F₁ is the number of OTUs containing only one sequence; and F₂ is the number of OTUs containing only two sequences.

④ ACE index (S_ACE):

$S_{ACE}=S_{abund}+\frac{S_{rare}}{C_{ACE}}+\frac{F_1}{C_{ACE}}{\mathrm{\gamma }}_{ACE}^2$

Where, ${\mathrm{\gamma }}_{ACE}^2=\max [\frac{S_{rare}{\sum }_{i=1}^{10}i\left(i-1\right)F_i}{C_{ACE}{N}_{rare}\left(N_{rare}-1\right)}-1,0]$

$C_{ACE}=1-\frac{F_1}{N_{rare}}$

$N_{rare}={\sum }_{i=1}^{10}i{F}_i$

S_rare is the number of OTUs with no more than “abund”” sequences; S_abund is the number of OTUs with more than “abund” sequences; abund is the threshold for “dominance” OTUs. abund is the threshold of “dominance” OTUs, and the default is 10.

3 Results

3.1 Physical and chemical properties of soils at different altitudes

The physicochemical properties of soils at different altitudes on the eastern slope of Helan mountain exhibited significant variations (P < 0.05), as presented in Table 2. WC and TP displayed an increasing trend with altitude, reaching their highest values in the high-altitude area (13.07% and 161.84 mg/kg, respectively). On the other hand, OM, TN, AN, AP, and AK demonstrated a pattern of initially increasing and then decreasing along the altitude gradient. The middle altitude area had the highest levels of OM, AN, and TN (9.25%, 137.26 mg/kg, and 165.77 mg/kg, respectively), while AP reached its peak in the high-altitude area (12.43 mg/kg). AK, on the other hand, exhibited its highest concentration in the low-altitude region (49.83 mg/kg). Soil pH decreased with altitude, ranging from 7.95 to 8.13, with minor fluctuations, indicating a slightly alkaline nature of the soil on the eastern slope of Helan mountain. In summary, the nutrient content of soils on the eastern slope of Helan mountain increased with altitude, in agreement with the findings of previous studies.

TABLE 2

Table 2. Physicochemical properties of soils of different vegetation types.

3.2 α-Diversity of soil AMF communities at different altitude gradients of Helan mountain

As shown in Figure 1, when the end of the dilution curve of the Shannon index tends to flatten toward flat, it indicates that the amount of sequencing data in the sample is reasonable enough. That is, the biodiversity in the sample is high, and there exists a relative abundance of different species or gene expression levels that can be reasonably analyzed for subsequent data.

FIGURE 1

Figure 1. Dilution curve of soil AMF community.

Table 3 displays significant differences (P < 0.05) in the α-diversity of AMF communities among different vegetation types and altitude gradients. Specifically, the soil AMF community in GR exhibited the highest Shannon index and the lowest Simpson index, while the soil AMF community in CB showed the lowest Shannon index and the highest Simpson index. This indicates that the α-diversity of soil AMF communities was highest in desert grasslands and lowest in mixed coniferous forests. Regarding the ACE index and Chao1 index, UL had the highest values, whereas CB had the lowest. No significant differences were observed among vegetation types at different altitudes.

TABLE 3

Table 3. Soil AMF diversity index for different vegetation types.

In summary, the α-diversity of soil AMF communities on the eastern slope of the Helan mountains exhibited significant variations among different vegetation types. The overall trend indicated a decrease followed by an increase in altitude. Overall, gray elm forests showed the highest alpha diversity among soil AMF communities, while mixed coniferous and broadleaf forests displayed the lowest diversity.

3.3 Taxonomic composition and distribution of soil AMF at different altitudes of Helan mountain

A total of 986 OTUs of AMF were identified through high-throughput sequencing in this study. These OTUs belonged to 1 phylum, 4 classes, 6 orders, 12 families, 14 genera, and 114 species.

Figure 2 illustrates significant variations in the relative abundance of soil AMF at the genus level among samples at different altitudes. Among all the samples, Glomus exhibited the highest relative abundance, with an average percentage of 85.13%. The highest relative abundance (98.27%) was observed in sample SH at a low altitude, indicating that Glomus was the dominant species in the soil AMF community on the eastern slope of the Helan mountains. In addition, the 2nd to 5th relative abundances were unclassified__p__Glomeromycota, Claroideoglomus, Archaeospora, and unclassified__o_Archaeosporales in the order of relative abundance. As shown in Figure 3, further analysis of the species differences at the genus level across different samples revealed that Glomus, unclassified_p_Glomeromycota, Claroideoglomus, and Diversispora were present in all altitudinal gradient samples, indicating their wide distribution on the eastern slopes of Mount Helan. Sclerocystis was found exclusively in sample UL, while unclassified_p_Glomeromycota, Archaeospora, and Paraglomus exhibited significant differences (P < 0.05) among samples.

FIGURE 2

Figure 2. Community composition of soil AMF.

FIGURE 3

Figure 3. Significance test for difference in AMF species at different altitude. The symbol “*” heatmap between soil factors represents significant correlations at the 0.05 level.

Comparative analysis of AMF community species diversity among different samples was conducted at the OTU level. The results of NMDS analysis, utilizing the Bray–Curtis distance algorithm, and the ANOSIM test revealed significant differences in the community structure of AMF across the samples (R = 0.5556, P = 0.001). As depicted in Figure 4, the samples exhibited distinct clustering patterns in a two-dimensional space. The PF sample points at the middle altitude appeared more scattered, indicating lower species similarity among those samples. Conversely, the SH and GR sample points at the low altitude were more closely grouped, suggesting higher species similarity. Notably, the NMDS1 axis clearly distinguished the AMF communities between high and low altitudes. The greatest species differences were observed between the SH and GR communities at low altitudes and the PF and CB communities at medium altitudes.

FIGURE 4

Figure 4. Non-metric multidimensional scaling analysis of β-diversity of soil AMF communities.

To further investigate the differences between the PF and GR samples, LEfse analysis (Figure 5) was performed. The analysis revealed that the abundance of Glomus was significantly higher in GR compared to PF. Conversely, the abundance of unclassified__p__Glomeromycota, Archaeospora, and Paraglomus was significantly higher in PF compared to GR.

FIGURE 5

Figure 5. LEfSe bar analysis of GR vs. PF.

3.4 Correlation analysis between α-diversity of soil AMF communities and soil physicochemical properties at different altitude gradients of Helan mountain

The results of Pearson correlation analysis (Table 4) revealed significant associations between soil AMF community α-diversity indices and various soil factors. The Shannon index exhibited a highly significant positive correlation with soil OM and a significant positive correlation with WC. Conversely, the Simpson index showed a highly significant negative correlation with OM and a significant negative correlation with WC. The ACE index was significantly positively correlated with ALT and OM, while the Chao1 index demonstrated a highly significant positive correlation with WC and OM, as well as a significant positive correlation with ALT.

TABLE 4

Table 4. Person correlation analysis between soil physical and chemical properties and soil AMF community α-diversity.

To explore the relationship between soil AMF community α-diversity and environmental factors, redundancy analysis (RDA) analysis was conducted (Figure 6). The first and second ordination axes explained 86.60% and 3.63% of the variance in AMF community composition, respectively. The Shannon, ACE, and Chao indices exhibited positive correlations with AP, OM, ALT, WC, total TN, AN, and TP, while showing negative correlations with AK and pH. Simpson’s index displayed the opposite pattern. The results of the Monte Carlo significance test indicated that OM (r² = 0.3926, P = 0.003) and WC (r² = 0.2715, P = 0.049) were the primary environmental factors influencing the α-diversity of AMF communities.

FIGURE 6

Figure 6. Redundancy analysis of soil environmental factors and diversity coefficients of AMF communities.

The results of the heatmap plot of environmental factors and soil AMF genus level showed (Figure 7) that different environmental factors had different degrees of influence on soil AMF genera. Among them, soil OM showed a highly significant positive correlation with the genera unclassified_ p__Glomeromycota and Archaeospora, a significant positive correlation with the genera Paraglomus and unclassified_o_Glomerales, and a highly significant negative correlation with Glomus. WC was positively correlated with Ambispora and Claroideoglomus, AN was positively correlated with Paraglomus, Scutellospora was most affected by ALT and AK, and Claroideoglomus was most affected by WC and pH. In addition, the effects of environmental factors on unclassified__o_Archaeosporales, Acaulospora, Sclerocystis, unclassified__c__Glomeromycetes, unclassified__f__Glomeraceae, and Diversispora had no significant effect.

FIGURE 7

Figure 7. Person correlation heatmap between soil factors and AMF community. The symbols * and ** represent significant correlations at 0.05 and 0.01 levels, respectively.

4 Discussion

4.1 Characterization of AMF community structure in soils at different altitudes

Consistent with the results of many studies (Zhang M. G. et al., 2022), Glomus, as the most common genus of AMF, is widely distributed in the natural environment (Schwarzott et al., 2001). In our study, Glomus also emerged as the dominant genus in the AMF community of the soil on the east slope of Helan mountain. This can be attributed to Glomus’ exceptional reproductive capacity and its ability to establish mycorrhizal associations rapidly with host plants. These symbiotic relationships enhance the nutrient and water absorption capabilities of host plants, promote plant growth efficiency, and improve stress resistance (Beimforde et al., 2014; Zhao et al., 2017). Consequently, Glomus contributes to the stability of plant communities and indirectly enriches soil biodiversity. Precipitation enhances vegetation photosynthesis efficiency, leading to increased biomass in host plants and providing more nutrients for AMF (Manzoni et al., 2012). However, as altitude rises, precipitation and soil WC increase and plant roots become less reliant on AMF for water uptake. Furthermore, the decrease in temperature at higher altitudes hinders the growth and development of AMF (Yang et al., 2022). Since Glomus encompasses a greater number of species, the adverse effects of environmental changes associated with increased altitude are more pronounced for this genus compared to others (Zhang M. G. et al., 2022).

Significant differences in soil AMF community α-diversity were observed among different altitudes. Consistent with previous studies (Guo et al., 2020), the α-diversity of soil AMF communities on the eastern slopes of the Helan mountains exhibited a positive correlation with ALT. This suggests that altitude and topography influence variations in soil properties and vegetation types, which, in turn, impact the AMF and mycorrhizal environments. The influence of host plants on soil AMF community α-diversity is substantial (Li et al., 2014; Davison et al., 2015). The Helan mountains exhibit distinct vertical distribution patterns in vegetation zones, giving rise to significant differences in soil AMF across these zones (Li et al., 2015). Research has established a direct correlation between soil AMF community diversity and plant community diversity and stability (Stover et al., 2018; Zhao et al., 2022). The region with the highest soil AMF diversity on the eastern slopes of the Helan mountains is situated in the lower altitude GR. Although the vegetation diversity and richness in the GR are not particularly high, the presence of dense shrubs and grasses, coupled with lower altitudes and favorable temperatures, fosters increased wildlife activity and appropriate disturbances. These factors accelerate material cycling, and energy flow, and increase the diversity of soil AMF. UL located at mid-altitude showed the highest AMF abundance among the seven experimental sample sites. This may be because the diversity and richness of vegetation in this region is significantly higher than that at lower altitudes, providing a wider range of hosting options for AMF. The PI also showed relatively higher abundance and diversity at higher altitudes. This may be because the PI is a relatively mature forest area with a more stable plant community. It may also be that in arid and semi-arid regions, mycorrhizal hyphae are more important to plants due to low soil WC (Bahadur et al., 2019). The increase in soil WC at high altitudes promotes the reproduction and growth of plant-AMF symbionts.

4.2 Main factors influencing soil AMF community structure

The composition and distribution of AMF communities on the eastern slopes of the Helan mountains are significantly influenced by differences in soil physical and chemical properties, as indicated by the findings of this study. Correlation analysis revealed that soil OM and WC exerted a stronger impact on the α-diversity of AMF communities compared to other soil physical and chemical properties. This can be attributed to WC promoting the reproductive development of AMF and influencing the selection and colonization of AMF by host plants (Shukla et al., 2013). OM serves as a carbon source necessary for microbial reproduction and positively affects AMF spore density and colonization, thereby enhancing AMF reproductive capacity to a certain extent (Bending et al., 2002; Shao et al., 2019; Sadeghi et al., 2023). Interestingly, the relative abundance of Glomus, the dominant genus in the study area, exhibited a pattern of decreasing and then increasing along different altitudinal gradients, which contradicted the trend observed for OM. This suggests that factors such as temperature at lower altitudes and soil moisture content at higher altitudes positively promote the growth and reproduction of Glomus as the dominant genus. It indicates that Glomus is not solely reliant on soil OM and possesses adaptability to nutrient-poor soil environments. Alternatively, the growth of Glomus may be limited by chemical inhibition or nutrient competition from plants (Singh et al., 2007).

This study conducted Person correlation analysis and found a highly significant positive correlation between altitude and soil OM, demonstrating the significant influence of altitude on the α-diversity of AMF communities. As a bridge for nutrient conversion between plant roots and the soil environment, the diversity of AMF is inevitably affected by inter-root soil OM (Zinger et al., 2011). Excessively high or low nutrient levels can inhibit AMF colonization and growth, with soil nitrogen and OM content playing a key role in determining the abundance and diversity of AMF (Xiao et al., 2018). Although soil nitrogen content was not the primary factor influencing AMF diversity and richness in this study, TN showed a significant positive correlation with altitude in Person’s correlation analysis, indicating that altitude indirectly influenced the α-diversity of AMF communities by affecting soil nitrogen content. Contrary to the results of previous studies by Shen et al. (2013), Guo and Gong (2014), and Li et al. (2022), this study found that pH was not the main factor influencing the AMF community structure on the east slope of Helan mountain. This discrepancy may be attributed to the limited range of soil pH variation in the samples of the present experiment.

Among the seven sample sites in this study, the middle altitude area, represented by the UL, exhibited medium levels of soil physicochemical indices, and the richness and diversity of the AMF community were relatively high across all sample sites. In the high-altitude region, represented by the PI, all soil physicochemical properties were at the upper range among the experimental sample plots, but the AMF community diversity was not high. In the low-altitude regions (SH and GR), the levels of soil physical and chemical properties were relatively lower, yet the AMF community diversity was high. This may be attributed to the fact that low and middle altitudes have more suitable temperatures that positively impact AMF compared to soil nutrients and moisture. Additionally, the timing of sampling may have influenced the results (Liu et al., 2021).

5 Conclusion

The present study investigated the diversity and abundance of soil AMF along altitudinal gradients on the eastern slopes of the Helan mountains. The findings are as follows:

1. Significant differences in soil physicochemical properties were observed at different altitudes. Soil moisture increased with altitude, while nutrient levels showed an initial increase followed by a decrease.

2. High-throughput sequencing detected a total of 986 OTUs belonging to 1 phylum, 4 classes, 6 orders, 12 families, 14 genera, and 114 species. The dominant genus was Glomus, accounting for an average proportion of 85.13%. The α-diversity of soil AMF communities exhibited significant differences across different altitudes, generally showing an increasing trend followed by a decrease.

3. Person correlation analysis and RDA analysis revealed that soil OM and WC were the primary factors influencing the diversity and abundance of soil AMF communities on the eastern slopes of the Helan mountains. Nitrogen, phosphorus, and quick-acting potassium were identified as secondary factors affecting the diversity and abundance of soil AMF communities.

As the study progresses, a deeper understanding of soil microbial diversity, distribution patterns, and the underlying mechanisms in the Helan mountains will be revealed. This knowledge will enable better utilization of the ecological functions of soil microbes and ultimately contribute to environmental improvement and ecological restoration.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI – https://www.ncbi.nlm.nih.gov/sra/PRJNA1109613.

Author contributions

LS: Data curation, Formal analysis, Methodology, Software, Writing – original draft, Writing – review & editing. PY: Data curation, Methodology, Software, Writing – review & editing. SY: Investigation, Writing – review & editing. HB: Investigation, Writing – review & editing. RZ: Investigation, Writing – review & editing. GS: Investigation, Writing – review & editing. YH: Investigation, Writing – review & editing. DP: Investigation, Writing – review & editing. XLN: Funding acquisition, Project administration, Supervision, Writing – review & editing. XN: Funding acquisition, Project administration, Supervision, 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 Ningxia Science Foundation of Ningxia Province (Award Number: 2023AAC02020), Ningxia Key Research and Development Plan project (Award Number: 2020bfg03006), the central government guides local construction projects (Award Number: 2022FRG05001).

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.

Footnotes

1. ^ http://drive5.com/uparse/

References

Bahadur, A., Batool, A., Nasir, F., Jiang, S. J., Qin, M. S., Zhang, Q., et al. (2019). Mechanistic insights into arbuscular mycorrhizal fungi-mediated drought stress tolerance in plants. Int. J. Mol. Sci. 20:18. doi: 10.3390/ijms20174199

PubMed Abstract | Crossref Full Text | Google Scholar

Barcenas-Moreno, G., Gomez-Brandon, M., Rousk, J., and Baath, E. (2009). Adaptation of soil microbial communities to temperature: Comparison of fungi and bacteria in a laboratory experiment. Glob. Change Biol. 15, 2950–2957. doi: 10.1111/j.1365-2486.2009.01882.x

Crossref Full Text | Google Scholar

Beimforde, C., Feldberg, K., Nylinder, S., Rikkinen, J., Tuovila, H., Dorfelt, H., et al. (2014). Estimating the Phanerozoic history of the Ascomycota lineages: Combining fossil and molecular data. Mol. Phylogenet. Evol. 78, 386–398. doi: 10.1016/j.ympev.2014.04.024

PubMed Abstract | Crossref Full Text | Google Scholar

Bending, G. D., Turner, M. K., and Jones, J. E. (2002). Interactions between crop residue and soil organic matter quality and the functional diversity of soil microbial communities. Soil Biol. Biochem. 34, 1073–1082. doi: 10.1016/s0038-0717(02)00040-8

Crossref Full Text | Google Scholar

Bonfim, J. A., Vasconcellos, R. L. F., Gumiere, T., Mescolotti, D. D. C., Oehl, F., and Cardoso, E. J. B. (2016). Diversity of arbuscular mycorrhizal fungi in a brazilian atlantic forest Toposequence. Microb. Ecol. 71, 164–177. doi: 10.1007/s00248-015-0661-0

PubMed Abstract | Crossref Full Text | Google Scholar

Bryant, J. A., Lamanna, C., Morlon, H., Kerkhoff, A. J., Enquist, B. J., and Green, J. L. (2008). Microbes on mountainsides: Contrasting elevational patterns of bacterial and plant diversity. Proc. Natl. Acad. Sci. U.S.A. 105, 11505–11511. doi: 10.1073/pnas.0801920105

PubMed Abstract | Crossref Full Text | Google Scholar

Dassen, S., Cortois, R., Martens, H., de Hollander, M., Kowalchuk, G. A., van der Putten, W. H., et al. (2017). Differential responses of soil bacteria, fungi, archaea and protists to plant species richness and plant functional group identity. Mol. Ecol. 26, 4085–4098. doi: 10.1111/mec.14175

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

Dowarah, B., Gill, S. S., and Agarwala, N. (2022). Arbuscular Mycorrhizal fungi in conferring tolerance to biotic stresses in plants. J. Plant Growth Regul. 41, 1429–1444. doi: 10.1007/s00344-021-10392-5

Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

Egan, C. P., Callaway, R. M., Hart, M. M., Pither, J., and Klironomos, J. (2017). Phylogenetic structure of arbuscular mycorrhizal fungal communities along an elevation gradient. Mycorrhiza 27, 273–282. doi: 10.1007/s00572-016-0752-x

PubMed Abstract | Crossref Full Text | Google Scholar

Genre, A., Lanfranco, L., Perotto, S., and Bonfante, P. (2020). Unique and common traits in mycorrhizal symbioses. Nat. Rev. Microbiol. 18, 649–660. doi: 10.1038/s41579-020-0402-3

PubMed Abstract | Crossref Full Text | Google Scholar

Gomoryova, E., Strelcova, K., Skvarenina, J., and Gomory, D. (2013). Responses of soil microorganisms and water content in forest floor horizons to environmental factors. Eur. J. Soil Biol. 55, 71–76. doi: 10.1016/j.ejsobi.2012.12.001

Crossref Full Text | Google Scholar

Guo, X. H., and Gong, J. (2014). Differential effects of abiotic factors and host plant traits on diversity and community composition of root-colonizing arbuscular mycorrhizal fungi in a salt-stressed ecosystem. Mycorrhiza 24, 79–94. doi: 10.1007/s00572-013-0516-9

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, Y. N., Zhang, H. D., Bao, Y. Y., Tan, H. Z., Liu, X. H., and Rahman, Z. U. (2022). Distribution characteristics of soil AM fungi community in soft sandstone area. J. Environ. Manag. 316:8. doi: 10.1016/j.jenvman.2022.115193

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, Y. X., Ren, C. J., Yi, J. J., Doughty, R., and Zhao, F. Z. (2020). Contrasting responses of rhizosphere bacteria, fungi and arbuscular mycorrhizal fungi along an elevational gradient in a temperate Montane forest of China. Front. Microbiol. 11:2042. doi: 10.3389/fmicb.2020.02042

PubMed Abstract | Crossref Full Text | Google Scholar

Hashem, A., Akhter, A., Alqarawi, A. A., Singh, G., Almutairi, K. F., and Abd Allah, E. F. (2021). Mycorrhizal fungi induced activation of tomato defense system mitigates Fusarium wilt stress. Saudi J. Biol. Sci. 28, 5442–5450. doi: 10.1016/j.sjbs.2021.07.025

PubMed Abstract | Crossref Full Text | Google Scholar

Jiang, Y., Kang, M. Y., Liu, S., Tian, L. S., and Lei, M. D. (2000). A study on the vegetation in the east side of Helan Mountain. Plant Ecol. 149, 119–130. doi: 10.1023/A:1026571016286

Crossref Full Text | Google Scholar

Koide, R. T., and Mosse, B. (2004). A history of research on arbuscular mycorrhiza. Mycorrhiza 14, 145–163. doi: 10.1007/s00572-004-0307-4

PubMed Abstract | Crossref Full Text | Google Scholar

Kowalchuk, G. A., De Souza, F. A., and Van Veen, J. A. (2002). Community analysis of arbuscular mycorrhizal fungi associated with Ammophila arenaria in Dutch coastal sand dunes. Mol. Ecol. 11, 571–581. doi: 10.1046/j.0962-1083.2001.01457.x

PubMed Abstract | Crossref Full Text | Google Scholar

Li, H., Wang, X. G., Liang, C., Hao, Z. Q., Zhou, L. S., Ma, S., et al. (2015). Aboveground-belowground biodiversity linkages differ in early and late successional temperate forests. Sci. Rep. 5:11. doi: 10.1038/srep12234

PubMed Abstract | Crossref Full Text | Google Scholar

Li, P., Yin, R. B., Zhou, H. M., Yuan, X. Y., and Feng, Z. Z. (2022). Soil pH drives poplar rhizosphere soil microbial community responses to ozone pollution and nitrogen addition. Eur. J. Soil Sci. 73:14. doi: 10.1111/ejss.13186

Crossref Full Text | Google Scholar

Li, T., and Zhao, Z. W. (2005). Arbuscular mycorrhizas in a hot and arid ecosystem in southwest China. Appl. Soil Ecol. 29, 135–141. doi: 10.1016/j.apsoil.2004.11.005

Crossref Full Text | Google Scholar

Li, X. L., Gai, J. P., Cai, X. B., Li, X. L., Christie, P., Zhang, F. S., et al. (2014). Molecular diversity of arbuscular mycorrhizal fungi associated with two co-occurring perennial plant species on a Tibetan altitudinal gradient. Mycorrhiza 24, 95–107. doi: 10.1007/s00572-013-0518-7

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, R. C., Xiao, Z. Y., Hashem, A., Abd Allah, E. F., and Wu, Q. S. (2021). Mycorrhizal fungal diversity and its relationship with soil properties in Camellia oleifera. Agric. Basel 11:10. doi: 10.3390/agriculture11060470

Crossref Full Text | Google Scholar

Lv, Y. Z., and Li, B. G. (2010). Soil science experiment. Beijing: China Agricultural Press.

Google Scholar

Manzoni, S., Schimel, J. P., and Porporato, A. (2012). Responses of soil microbial communities to water stress: Results from a meta-analysis. Ecology 93, 930–938. doi: 10.1890/11-0026.1

PubMed Abstract | Crossref Full Text | Google Scholar

Prommer, J., Walker, T. W. N., Wanek, W., Braun, J., Zezula, D., Hu, Y. T., et al. (2020). Increased microbial growth, biomass, and turnover drive soil organic carbon accumulation at higher plant diversity. Glob. Change Biol. 26, 669–681. doi: 10.1111/gcb.14777

PubMed Abstract | Crossref Full Text | Google Scholar

Prudent, M., Dequiedt, S., Sorin, C., Girodet, S., Nowak, V., Duc, G., et al. (2020). The diversity of soil microbial communities matters when legumes face drought. Plant Cell Environ. 43, 1023–1035. doi: 10.1111/pce.13712

PubMed Abstract | Crossref Full Text | Google Scholar

Ren, C. J., Zhou, Z. H., Guo, Y. X., Yang, G. H., Zhao, F. Z., Wei, G. H., et al. (2021). Contrasting patterns of microbial community and enzyme activity between rhizosphere and bulk soil along an elevation gradient. Catena 196:11. doi: 10.1016/j.catena.2020.104921

Crossref Full Text | Google Scholar

Sadeghi, S., Petermann, B. J., Steffan, J. J., Brevik, E. C., and Gedeon, C. (2023). Predicting microbial responses to changes in soil physical and chemical properties under different land management. Appl. Soil Ecol. 188:12. doi: 10.1016/j.apsoil.2023.104878

Crossref Full Text | Google Scholar

Schleuss, U., and Muller, F. (2001). Requirements for soil sampling in the context of ecosystem research. Sci. Total Environ. 264, 193–197. doi: 10.1016/s0048-9697(00)00624-0

PubMed Abstract | Crossref Full Text | Google Scholar

Schwarzott, D., Walker, C., and Schussler, A. (2001). Glomus, the largest genus of the arbuscular mycorrhizal fungi (Glomales), is nonmonophyletic. Mol. Phylogenet. Evol. 21, 190–197. doi: 10.1006/mpev.2001.1007

PubMed Abstract | Crossref Full Text | Google Scholar

Shao, P. S., Liang, C., Rubert-Nason, K., Li, X. Z., Xie, H. T., and Bao, X. L. (2019). Secondary successional forests undergo tightly-coupled changes in soil microbial community structure and soil organic matter. Soil Biol. Biochem. 128, 56–65. doi: 10.1016/j.soilbio.2018.10.004

Crossref Full Text | Google Scholar

Shen, C. C., Liang, W. J., Shi, Y., Lin, X. G., Zhang, H. Y., 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. B., Zhang, H. Y., Feng, Y. Z., Lin, X. G., Li, X. Y., 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

Shukla, A., Vyas, D., and Jha, A. (2013). Soil depth: An overriding factor for distribution of arbuscular mycorrhizal fungi. J. Soil Sci. Plant Nutr. 13, 23–33. doi: 10.4067/s0718-95162013005000054

PubMed Abstract | Crossref Full Text | Google Scholar

Singh, S., Ghoshal, N., and Singh, K. P. (2007). Variations in soil microbial biomass and crop roots due to differing resource quality inputs in a tropical dryland agroecosystem. Soil Biol. Biochem. 39, 76–86. doi: 10.1016/j.soilbio.2006.06.013

Crossref Full Text | Google Scholar

Stover, H. J., Naeth, M. A., and Boldt-Burisch, K. (2018). Soil disturbance changes arbuscular mycorrhizal fungi richness and composition in a fescue grassland in Alberta Canada. Appl. Soil Ecol. 131, 29–37. doi: 10.1016/j.apsoil.2018.07.008

Crossref Full Text | Google Scholar

Tedersoo, L., Bahram, M., and Zobel, M. (2020). How mycorrhizal associations drive plant population and community biology. Science 367:aba1223. doi: 10.1126/science.aba1223

PubMed Abstract | Crossref Full Text | Google Scholar

Uhlmann, E., Gorke, C., Petersen, A., and Oberwinkler, F. (2004). Arbuscular mycorrhizae from semiarid regions of Namibia. Can. J. Bot. 82, 645–653. doi: 10.1139/b04-039

Crossref Full Text | Google Scholar

Wang, J. J., Hu, A., Meng, F. F., Zhao, W. Q., Yang, Y. F., Soininen, J., et al. (2022). Embracing mountain microbiome and ecosystem functions under global change. New Phytol 234, 1987–2002. doi: 10.1111/nph.18051

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, J. T., Cao, P., Hu, H. W., Li, J., Han, L. L., Zhang, L. M., et al. (2015). Altitudinal distribution patterns of soil bacterial and archaeal communities along mt. Shegyla on the tibetan plateau. Microb. Ecol. 69, 135–145. doi: 10.1007/s00248-014-0465-7

PubMed Abstract | Crossref Full Text | Google Scholar

Xiao, S. H., Zhang, Z. R., You, W. B., Liu, J. H., Wu, J. Q., Cai, C. T., et al. (2018). Soil microbial community composition in four Nothotsuga longibracteata forests in Southern China. Pol. J. Environ. Stud. 27, 917–925. doi: 10.15244/pjoes/76116

Crossref Full Text | Google Scholar

Yang, J., Yu, J., Wang, J., and Li, X. (2022). Soil Microbial community responds to elevation gradient in an arid Montane ecosystem in Northwest China. Forests 13:1360. doi: 10.3390/f13091360

Crossref Full Text | Google Scholar

Yu, L. M., Zhang, Z. F., Zhou, L. W., and Huang, K. C. (2023). Effects of altitude and continuous cropping on arbuscular mycorrhizal fungi community in Siraitia grosvenorii rhizosphere. Agric. Basel 13:18. doi: 10.3390/agriculture13081548

Crossref Full Text | Google Scholar

Zhang, R. J., Tian, X. R., Xiang, Q. J., Penttinen, P., and Gu, Y. F. (2022). Response of soil microbial community structure and function to different altitudes in arid valley in Panzhihua, China. BMC Microbiol. 22:11. doi: 10.1186/s12866-022-02500-6

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, M. G., Yang, M., Shi, Z. Y., Gao, J. K., and Wang, X. G. (2022). Biodiversity and variations of arbuscular mycorrhizal fungi associated with roots along elevations in Mt. Taibai of China. Divers. Basel 14:13. doi: 10.3390/d14080626

Crossref Full Text | Google Scholar

Zhao, H., Li, X. Z., Zhang, Z. M., Zhao, Y., Yang, J. T., and Zhu, Y. W. (2017). Species diversity and drivers of arbuscular mycorrhizal fungal communities in a semi-arid mountain in China. PeerJ 5:18. doi: 10.7717/peerj.4155

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, W., Yin, Y. L., Li, S. X., Liu, J. J., Dong, Y. L., and Su, S. F. (2022). Soil microbial community varied with vegetation types on a small regional scale of the Qilian mountains. Sustainability 14:15. doi: 10.3390/su14137910

Crossref Full Text | Google Scholar

Zinger, L., Lejon, D. P. H., Baptist, F., Bouasria, A., Aubert, S., Geremia, R. A., et al. (2011). Contrasting Diversity patterns of crenarchaeal, bacterial and fungal soil communities in an alpine landscape. PLoS One 6:e0019950. doi: 10.1371/journal.pone.0019950

PubMed Abstract | Crossref Full Text | Google Scholar