Diversity characteristics of arbuscular mycorrhizal fungi communities in the soil along successional altitudes of Helan Mountain, arid, and semi-arid regions of China

Introduction: Arbuscular mycorrhizal fungi (AMF) perform a vital role in terrestrial ecosystems.

Methods: To investigate the diversity of AMF communities on the western slope of Helan Mountain at different altitudes and their influence factors, high-throughput sequencing was used to study the structure and diversity of soil AMF communities under different environments and their interrelationships between AMF and environmental factors.

Results: The results revealed that there were significant differences (p < 0.05) in the physical and chemical properties of the soil along the different altitudes. A total of 1,145 OTUs were obtained by high-throughput sequencing, belonging to 1 phylum, 4 class, 6 orders, 13 families, 18 genera and 135 species, with the dominant genus being Glomus, which accounted for 75.27% of the relative abundance of the community. Soil AMF community structure was shown to be variable at the generic level according to NMDS analysis. Correlation analysis showed that soil pH, water content (WC), organic matter (OM), available K, available P and N were significantly correlated with AMF community diversity and species abundance (p < 0.05, p < 0.01). Based on redundancy analysis (RDA) and Monte Carlo test results, soil pH, WC and OM had highly significant effects (p < 0.01) on AMF community diversity and species abundance.

Discussion: This study investigates the relationship between AMF community structure and diversity and soil physicochemical properties at different elevations on the western slope of Helan Mountain, which is of great significance to the study of the Helan Mountain ecosystem.

1. Introduction

Arbuscular mycorrhizal fungi (AMF) are mycorrhizal fungi that exist in the inter-rooted soil and can form symbiotic relationships with 80% of terrestrial vascular plants (Davison et al., 2015; Shi et al., 2020) and are widely found in various terrestrial ecosystems, playing an irreplaceable role in maintaining and improving vegetation communities, soil conditions and stabilizing ecosystem functions, and are an important component of terrestrial ecosystems (Elhindi et al., 2018; Qiang et al., 2019). It has been found that AMF can regulate the composition of soil inter-root microbial communities, soil structure and nutrient cycling (Cheng et al., 2012; Veresoglou et al., 2012); and after symbiosis with plants, it can not only improve the uptake of N, P and other nutrients by plants, but also enhance the resistance of plants to drought, salinity and pests and diseases, and promote plant growth (Auge, 2001; Garg and Chandel, 2011; Hodge and Storer, 2015). AMF communities are also influenced by a variety of environmental factors, such as altitude, plant communities, soil physicochemical properties and climatic factors. Especially, soil C, N, and P content and host plant species have a significant impact on the diversity and abundance of AMF communities (Kruger et al., 2017; Ezeokoli et al., 2019). Therefore, research on the relationship between AMF communities and environmental factors contributes to a better understanding of terrestrial ecosystems.

In mountain ecosystems, environmental factors can change considerably at smaller scales with altitude, causing differences in plant and animal communities between different altitudes (Deng et al., 2020; Zanzottera et al., 2020). Both changes in plant communities and soil characteristics affect the community structure of microorganisms in the soil (Li et al., 2018). Vieira et al. (2019) showed that AMF community diversity decreased with elevation increased in the zone of temperate climates. However, some studies have shown that AMF community diversity is independent or negatively correlated with elevations (Coutinho et al., 2015; Shen et al., 2015). In addition, slope orientation, an important topographic factor, and soil physicochemical properties and vegetation type differences between shady and sunny slopes (Sun et al., 2019; Guo et al., 2022) may affect AMF community structure. In conclusion, the different vegetation types, temperature, light and soil characteristics that different altitudes have can affect the AMF community. In recent years, there have been many achievements in the study of AMF. Among them, Dong et al. (2023) showed that drought conditions can have an impact on AMF growth and symbiosis with plants. Zhang et al. (2021) showed the effect of altitude on AMF community in Qinling Mountain. Zhao et al. (2019) showed that plant community succession has an important influence on AMF community. Liu et al. (2020) showed the influence of season and soil on AMF community in mountainous areas. The results of these studies show that soil and vegetation are the main environmental factors affecting AMF community, and that climate, seasonality and precipitation all indirectly affect AMF community through soil, but the types, magnitude and trends of the effects of soil factors on AMF community vary from different studies and need to be explored in depth.

Helan Mountain, located in the northwest of Ningxia, is an important natural geographical boundary and the dividing line between grassland and desert in northwest China, and plays the role of an important ecological barrier with a unique ecosystem. The vegetation types vary significantly along the successional altitudes, and the vegetation types are in order of desertification grassland, mountain scrub, sparse forest grassland, mountain coniferous forest and subalpine scrub meadow, etc., along the altitude rises (Wu et al., 2021). Vegetation diversity is often closely related to AMF community diversity and there is a positive interaction between them. The rich AMF resources contained in different elevations of Helan Mountain may play an important role in maintaining and protecting the stability of different vegetation ecosystems in Helan Mountain, but there is still a gap in research on the AMF communities in Helan Mountain.

Therefore, this study uses high-throughput sequencing technology to study the diversity of AMF communities in rhizosphere soil at different altitudes on the western slope of Helan Mountain, and to analyze the correlation between the physical and chemical properties of rhizosphere soil and the structural characteristics of AMF communities on the western slope of Helan Mountain, to provide a theoretical basis for maintaining the stability and development of Helan Mountain ecosystems and exploring the role of AMF communities in terrestrial ecosystems.

2. Materials and methods

2.1. Study area

Surveyed sites are located on the western slope of Helan Mountain in Alxa Left Banner, Inner Mongolia Autonomous Region. It is located in the eastern edge of the Alxa Plateau, the west side of the Yinchuan Plain (between 38°27′-39°30′N and 105°41′-106°41′E), with a typical continental climate and mountainous climatic characteristics. The average annual temperature is 8.6°C, the annual precipitation is 200–400 mm, mainly concentrated in May–September, and the annual evaporation is high, reaching over 2,000 mm. The main plant species in the study area are Agropyron mongolicum (Keng), Stipa breviflora (Griseb), Prunus mongolica (Maxim), Caragana stenophylla (Pojark), Ulmus glaucescens (Franch), Juniperus rigida (Siebold & Zucc), Pinus tabuliformis (Carriere), and Picea crassifolia (Kom), etc.

2.2. Experimental design

Based on the field survey, eight altitudes (ALT) were selected along the altitude gradient according to different community types on the western slope of Helan Mountain from 1800 to 2,750 m above sea level in the summer of 2021, and three sample plots were set up within each vegetation zone with a sample area of 20 m × 20 m. The soil samples were collected by digging up 5 plant rhizosphere soil samples using a 4-cm diameter soil sampler according to the five-point sampling method, mixing them well, placing them in sterilized bags and storing them in a refrigerated incubator. The soil was taken back to the laboratory and stored in a refrigerator at −80°C for the determination of soil AMF community diversity and naturally dried and sieved for the determination of soil physical and chemical properties (Table 1).

TABLE 1

Table 1. Basic information about the sample site.

2.3. Measurement items and methods

Soil pH was determined using pH meter with a 5:1 ratio of water to soil; water content (WC) was determined using the drying method; organic matter (OM) was determined using the potassium dichromate dilution calorimetric method; total nitrogen (TN) was determined using the semi-micro Kjeldahl method; alkali-hydrolyzable nitrogen (AN) was determined by alkaline diffusion; total phosphorus (TP) was determined by HClO₄-H₂SO₄ digestion with molybdenum-antimony colorimetry; available phosphorus (AP) was determined by NaHCO₃ leaching with molybdenum-antimony colorimetry; available potassium (AK) was determined by NH₄OAc extraction and flame photometric method.

Soil AMF community diversity was sent to Majorbio Bio-Pharm Technology Co. Ltd (Shanghai, China) for determination, and soil AMF community DNA was extracted using the MP kit, using AMV4-5NF/AMDGR primers (Sato et al., 2005) with primer sequences of 5′- AAGCTCGTAGTTGAATTTCG-3′, 5’-CCCAACTATCCCTATTAATCAT-3′, and PCR amplification of the 18S rRNA gene was performed. PCR amplification cycling conditions were as follows: initial denaturation at 95°C for 3 min, followed by 27 cycles of denaturing at 95°C for 30 s, annealing at 55°C for 30 s and extension at 72°C for 45 s, and single extension at 72°C for 10 min, and end at 4°C.

Library construction of PCR products using TruSeq^TM DNA Sample Prep Kit: (1) Addition of official Illumina splice sequences to the outer end of the target region by PCR; (2) Recovery of PCR products by gel cutting using a gel recovery kit; (3) Elution in Tris–HCl buffer and detection by 2% agarose electrophoresis; (4) Denaturation by sodium hydroxide to produce single-stranded DNA fragments. Paired-end sequencing was performed on an Illumina MiSeq PE300 platform (Illumina, San Diego, United States). The raw sequencing reads were deposited into the NCBI Sequence Read Archive (SRA) database (Accession Number: PRJMNA932928).

Using UPARSE (Stackebrandt and Goebel, 1994; Edgar, 2013) software¹ (version 7.1), OTUs were classified at 97% sequence similarity and the RDP classifier Bayesian algorithm was used to taxonomically analyze representative sequences of OTUs at 97% similarity level by maarjam20220506/AM species database for comparison and to obtain taxonomic information.

2.4. Data analysis

Data were collated using EXCEL 2019 and one-way ANOVA (One-way ANOVA) and Duncan’s test (Duncan’s test) were used to analyze the differences in soil physicochemical properties between vegetation types using SPSS22 (p < 0.05), and one-tailed significance test was used to test Person between different variables. Correlation analysis was performed, and differences between AM fungal communities were analyzed using non-metric multidimensional scaling analysis (NMDS) and analysis of similarity test (ANOSIM), with r and p values associated with NMDS calculated based on 999 permutations of ANOSIM. Origin 2019 was used for mapping. Soil AMF community diversity, β-diversity and structure were analyzed and mapped using the computing platform of Shanghai Meiji Biomedical Technology Company Limited and the R language “Vegan” package.

The AMF community diversity was analyzed using the Shannon index, the Simpson index, the Chao1 index and the ACE index.

2.4.1. Shannon-Wiener index

The formula is calculated as:

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

2.4.2. Simpson index

The formula is calculated as:

$D_{Simpson}=\frac{{{\displaystyle \sum }}_{i=1}^{S_{obs}}{n}_i\left(n_i-1\right)}{N\left(N-1\right)}$

Where, Sobs is the number of OTUs observed in practice; n_i is the number of sequences contained in the i-th OTU; and N is the number of all sequences.

2.4.3. Chao1 index

The formula is calculated as:

$S_{Chao1}=S_{obs}+\frac{n_1\left(n_1-1\right)}{2\left(n_2-1\right)}$

Where, S_Chao1 is the estimated number of OTUs; n₁ is the number of OTUs containing only one sequence; n₂ is the number of OTUs containing only two sequences.

2.4.4. ACE index

The formula is calculated as:

$f\left(x\right)=\left\{\begin{array}{c}{S}_{abund}+\frac{S_{rare}}{C_{ACE}}+\frac{n_1}{C_{ACE}}{\hat{\gamma }}_{ACE}{}^2,for{\gamma }_{ACE}<0.80\\ S_{abund}+\frac{S_{rare}}{C_{ACE}}+\frac{n_1}{C_{ACE}}{\tilde{\gamma }}_{ACE}{}^2,for{\hat{\gamma }}_{ACE}\ge 0.80\end{array}\right.$

Among these:

$N_{rare}={\displaystyle \sum }_{i=1}^{abund}i{n}_i,C_{ACE}=1-\frac{n_1}{N_{rare}}$

${\hat{\gamma }}^2{}_{ACE}=\max \left[\frac{S_{rare}}{C_{ACE}}\frac{{{\displaystyle \sum }}_{i=1}^{abund}i\left(i-1\right)n_i}{N_{rare}\left(N_{rare}-1\right)}-1\mathrm{,}0\right]$

${\tilde{\gamma }}^2{}_{ACE}=\max \left[{\tilde{\gamma }}^2{}_{ACE}\left\{1+\frac{N_{rare}\left(1-C_{ACE}\right){{\displaystyle \sum }}_{i=1}^{abund}i\left(i-1\right)n_i}{N_{rare}\left(N_{rare}-C_{ACE}\right)}\right\}\mathrm{,}0\right]$

Where, S_rare is the number of OTUs containing “abund” sequences or less than “abund”; S_abund is the number of OTUs with more than “abund” sequences; abund is the threshold of “dominant” OTUs, default is 10.

3. Results

3.1. Physico-chemical properties of soil at different altitudes

As can be seen from Table 2, there were significant differences in soil physicochemical properties between different altitudes (p < 0.05). WC, OM, AK and AN increased with the increasing altitude, and the maximum values occurring in the high-altitude region (they are 14.95, 11.70%, 54.27 mg/kg, 96.73 mg/kg respectively). Both pH and TP increased and then decreased with increasing altitude, with the maximum values occurring in the mid-altitude region (they are 8.20, 213.81 mg/kg respectively), where the pH varied from 7.94 to 8.20, indicating that the soil on the western slope of Helan Mountain is alkaline. TN and AP decreased and then increased with increasing altitude, with the maximum TN value occurring at high-altitude areas (202.61 mg/kg) and the maximum AP value at low-altitude areas (55.37 mg/kg). Overall, the nutrient content of the soil increases with altitude.

TABLE 2

Table 2. Physico-chemical properties of soil at different altitudes.

3.2. Soil arbuscular mycorrhizal fungi community diversity at different altitudes

As can be seen from Figure 1, the dilution curve flattens out as the number of sequenced bars increases, indicating that the amount of sample sequencing data is large enough to reflect the true situation of the AMF community in the soil and that subsequent data analysis can be carried out.

FIGURE 1

Figure 1. AMF community dilution curves for soil at different elevations.

A total of 2,158,539 valid sequences with a total base number of 466,973,432 bp and an average sequence length of 216 bp were obtained from the AMF assay for each sample. The coverage of each sample was above 99%, indicating that the sequencing results are representative of the AMF diversity in the samples. As can be seen from Table 3, there were significant differences in soil AMF community diversity between elevations (p < 0.05). The soil AMF community of SU had the largest Shannon index and the smallest Simpson index, while the soil AMF community of PI had the smallest Shannon index and the largest Simpson index. The maximum values of the ACE and Chao1 indices are at SU and the minimum values are at GR and PI, respectively. Overall, the diversity of the soil AMF community increases with altitude.

TABLE 3

Table 3. AMF diversity indices for soils at different altitudes.

3.3. Taxonomic composition and distribution of arbuscular mycorrhizal fungi in soil at different altitudes

The results of the high-throughput sequencing data analysis showed that the comparison with the maarjam20220506/AM species database revealed a total of one phylum, four class, six orders, 13 families, 18 genera and 135 species at different altitudes.

According to Figure 2, the relative abundance of soil AMF at the genus level varied between samples. The relative abundance of Glomus was the highest in all samples, averaging 75.27%, reaching a maximum relative abundance in GR (94.04%), and exceeding 90% in GR, UL and SH, indicating that Glomus is the dominant species in the soil AMF community on the western slope of Helan Mountain. Analysis of species differences at different sample subordinate levels showed (Figure 3) that Acaulospora, Sclerocystis, and Entrophospora were only found in the SU. Glomus, unclassified_p__Glomeromycota, Diversispora, Archaeospora, Acaulospora, and Sclerocystis were significantly different in different samples (p < 0.05).

FIGURE 2

Figure 2. Genus level composition of AMF communities in soil at different altitudes.

FIGURE 3

Figure 3. Significance test for differences in AMF species at different altitudes. *represent significant correlations at the 0.05, **represent significant correlations at the 0.01, and ***represent significant correlations at the 0.001.

NMDS of soil AMF community structure based on genus level and bray_curtis distance calculations showed (Figure 3) that the PI group was relatively more dispersed within its group, indicating poor intra-group repeatability; the remaining samples were relatively more concentrated within their group, indicating good inter-replicate similarity within the remaining sample groups (Figure 4). The larger area of overlap between the GR, UL, and SF groups at low-altitude indicates that there is less difference in community structure between the three samples, whereas there is no area of overlap between PI and SC at high-altitude, indicating that there is a greater difference in community structure between the two samples and a greater difference in community structure between the low and high altitudes. LEfse analysis was further carried out for the two groups of samples with the greatest differences (GR and SU; Figure 5), of which the abundance of Glomus was significantly higher in GR than in SU, and the abundance of unclassified_p__Glomeromycota, unclassified_o__Glomerales, Diversispora, Archaeospora, Paraglomus, and Ambispora was significantly higher in SU than in GR.

FIGURE 4

Figure 4. NMDS analysis of the β-diversity of soil AMF communities at different altitudes.

FIGURE 5

Figure 5. LEfSe analysis of GR and SU.

3.4. Correlation analysis between the physico-chemical properties of soil and the diversity of arbuscular mycorrhizal fungi communities at different altitudes

The correlation analysis showed (Table 4) that the ALT, WC, OM, and AK were highly significantly positively correlated with the Shannon Index, ACE Index and Chao1 Index, and highly significantly negatively correlated with the Simpson Index. The pH was highly significantly negatively correlated with the Shannon index, significantly positively correlated with the ACE and Chao1 indices and highly significantly positively correlated with the Simpson index. The AN was significantly negatively correlated with the Simpson Index and significantly positively correlated with the ACE Index. The TP was significantly and positively correlated with the Shannon index, the ACE index and the Chao1 index. The AP was highly significantly negatively correlated with the Shannon and Chao1 indices, significantly negatively correlated with the AVCE index and highly significantly positively correlated with the Simpson index.

TABLE 4

Table 4. Person correlation analysis between soil AMF community diversity and soil physicochemical properties.

The RDA showed (Figure 6) that the Shannon index, ACE index and Chao1 index were positively correlated with ALT, AK, WC, OM, TP, AN and TN, and negatively correlated with pH and AP, with ALT having the greatest effect on the structure of the AMF community. According to Monte Carlo tests, ALT (r² = 0.6623, p = 0.001), WC (r² = 0.5036, p = 0.002), OM (r² = 0.4820, p = 0.004) and pH (r² = 0.3888, p = 0.006) had significant effects on AMF community diversity and species abundance. The first ordination axis explained 72.57% of the variation in AMF community composition and the second ordination axis explained 11.79% of the variation in AMF community composition. The first ordination axis explains most of the variation in soil AMF community diversity and richness at different elevations on the western slope of Helan Mountain.

FIGURE 6

Figure 6. RDA between AMF diversity and physico-chemical properties of soil.

According to the results of Heatmap plots of environmental factors and soil AMF genus levels (Figure 7), soil physicochemical properties and altitude have different effects on each genus. Glomus was most affected by physico-chemical properties of soil and altitude, with Glomus being significantly negatively correlated with ALT, WC, OM, and AN, highly significantly positively correlated with pH, significantly negatively correlated with AK and TN, and significantly positively correlated with AP. The pH was highly significantly negatively correlated with unclassified_o__Glomerales and significantly positively correlated with unclassified_p__Glomeromycota. Acaulospora, Diversispora, and Pacisporawere highly significantly positively correlated with OM and AK, and highly significantly negatively correlated with AP. In addition, soil physicochemical properties and elevation had no significant effect on Ambispora, unclassified_f__Diversisporaceae and Scutellospora.

FIGURE 7

Figure 7. Heatmap of Person correlation between AMF community genus level and physico-chemical properties of soil. *represent significant correlations at the 0.05, **represent significant correlations at the 0.01, and ***represent significant correlations at the 0.001.

4. Discussion

4.1. Characteristics of changes in soil arbuscular mycorrhizal fungi community structure at different altitudes

A total of 1,418 AMF OTUs were obtained in this study, indicating that AMF species are abundant and widely distributed on the western slope of Helan Mountain, and are an important part of the ecosystem. The results of this study show that Glomus is the dominant genus in the soil AMF community on the western slope of Helan Mountain, which is consistent with the results of most studies on AMF communities in temperate mountainous areas (Zhang et al., 2022a,b), because Glomus contains a lots of species and can produce a large number of mycelial fragments and spores, allowing it to spread and colonize plant roots more widely and with greater environmental adaptability (Vandenkoornhuyse et al., 2002; Zhao et al., 2017). In the present study, the relative abundance of Glomus was highly significantly negatively correlated with altitude, which is consistent with the findings of Haug et al. (2019), while most other AMF genera were positively correlated with altitude. Several studies have shown that Glomus can enhance plant stress tolerance and biomass under drought conditions (Nacoon et al., 2021); appropriate increases in temperature and precipitation can enhance the photosynthetic efficiency of the host plant, allowing the plant to provide more photosynthetic products to AMF and promote mycelial and spore development (Zhao et al., 2022). As altitude increases, on the one hand, increased precipitation may make host plants less dependent on Glomus; on the other hand, lower temperatures limit the development of AMF mycelium and spores, which may have a greater impact on Glomus than other AMF genera. It also suggests that there is an effect of altitude on the community structure and intergeneric relationships of soil AMF communities.

In addition to differences in distribution, AMF community diversity varies during different altitudes. Some studies have shown (Egan et al., 2017; Shi et al., 2019) that AMF community diversity increases and then decreases with increasing altitude, with different correlations at different altitude regions. Then, our results showed that AMF community diversity increases with altitude, unlike the results of other studies (Lugo et al., 2008; Guo et al., 2022). This may be due to the small study area, the low altitude and limited soil nutrients in this study. It has been shown that both increases in soil nutrients and decreases in anthropogenic disturbance positively affect AMF community diversity (Opik et al., 2008; Helgason and Fitter, 2009). In this study, soil nutrients increased significantly with altitude and were less disturbed by anthropogenic disturbances, which together with a variety of environmental factors influenced the AMF community, resulting in significant differences in community diversity during different altitudes. Therefore, the mechanisms by which different altitudes and vegetation in Helan Mountain affect the AMF community need to be further explored in future studies.

4.2. Influencing factors of soil arbuscular mycorrhizal fungi community structure

The soil physico-chemical properties in Helan mountain have an important influence on AMF community diversity, in agreement with the findings of many studies (Ma et al., 2021; Maitra et al., 2021). Besides, pH, WC, AP, OM, and AK were the main influencing factors to affect the diversity of the AMF community. Furthermore, pH and AP were highly significantly and negatively correlated with community diversity. It is possible that pH can affect not only the formation and development of AMF spores (Zhang et al., 2022a), but also indirectly by affecting plant community composition and plant uptake and utilization of soil nutrients (Xu et al., 2017). The effect of AP on AMF community diversity is consistent with the findings of Ji and Bever (2016) and Ceulemans et al. (2019), possibly because elevated AP levels inhibit AMF spore germination and mycelial growth and cause plants to reduce root symbiosis with AMF and take up P directly from the soil (Johnson et al., 2015; Shi et al., 2021). In addition, WC, OM, and AK were highly significantly and positively correlated with community diversity, most likely because that the increased WC not only promotes the reproductive development of AMF (Zhao et al., 2017), but also changes the selection and colonization of AMF by host plants, which has a significant effect on AMF community structure and diversity (Symanczik et al., 2015). OM could improve AMF reproduction to some extent (Zhao et al., 2022), and AK can promote AMF infestation of plants and affect AMF communities. Soil N was a minor factor affecting the AMF community, unlike the results of some other study (Ji and Bever, 2016), probably because the western slope of the Helan Mountain is gently sloping, making the variation of soil N between different altitudes insignificant.

The soil characteristics has an important influence on the AMF community, and the analysis of Heatmap (Figure 7) indicated that soil physicochemical properties had different effects on different genera of AMF. Glomus, as the dominant genus, was most affected by soil physicochemical properties and most genera were affected by soil factors, in agreement with some studies (Xu et al., 2017; Kim et al., 2022), explaining to some extent the differences in AMF communities between different altitudes. It further demonstrates that soil factors have an important influence on AMF community structure and diversity.

In addition, it has been shown that host plants are one of the main factors affecting AMF communities in montane ecosystems (Martinez-Garcia et al., 2015), which may be related to the fact that different vegetation has different nutrient requirements (Ezeokoli et al., 2020), and that plant species diversity has been found to be positively correlated with fungal community diversity in some studies (Johnson et al., 2003). In the middle and low-altitude regions (GR-CO), AMF communities differed less structurally and more in diversity, with higher plant diversity in UL (Ulmus glaucescens Forests) and SH (Shrub) plants and lower plant diversity in GR (Grass; Yang et al., 2022), the same variation as in AMF communities diversity, indicating the influence of plant diversity on AMF community diversity. The relative abundance of Claroideoglomus was higher in the mid and high altitude regions (PF-SU) than in the low-altitude regions, and Claroideoglomus was found to promote nutrient uptake by plants (Jamiolkowska et al., 2020), which may allow rising soil nutrients in the mid and high altitude regions to promote symbiosis between plants and Claroideoglomus, allowing plants to take up more nutrients and increase their biomass. It has been found that under more severe climatic factors at high-altitude, plants will select more for symbiosis with AMFs that can enhance plant cold resistance (Bauer et al., 2017; Dong et al., 2021). This may account for the marked differences in AMF community structure in the high-altitude region (PI, SU) compared to the low-altitude region. In terms of community genus level composition, the relative abundance of unclassified_p_Glomeromycota, unclassified_o_Glomerales and unclassified_c_Glomeromycetes was significantly higher in the high-altitude region and may be key genera for improving plant cold resistance. However, further studies are needed to determine the specific species and functions. In addition, there were also significant differences between the AMF community structures of PI and SU, and the community diversity of SU was much higher than that of PI, probably due to the significant differences in vegetation types between SU and PI, with SU having significantly higher plant diversity and soil nutrients than PI, and the fact that Picea crassifolia in PI was mainly symbiotic with ectomycorrhizal fungi and the AMF community was not a major fungal phylum (Wang and Han, 2022). It further demonstrates the influence of vegetation type and plant diversity on AMF communities.

In future research, the interrelationship between soil AMF communities and plant community characteristics and soil nutrients at different altitudes on the western slope of Helan Mountain can be investigated to further demonstrate the ameliorating effect of soil AMF communities on plants and soil, and to provide a theoretical basis for the protection of Helan Mountain ecology, vegetation restoration and soil management.

5. Conclusion

In the present study, soil AMF on the western slope of Helan Mountain belonged to one phylum, four class, six orders, 13 families, 18 genera and 135 species, with the dominant genus being Glomus. Soil AMF community structure and diversity differed significantly between different altitudes, with AMF community diversity increasing with increasing altitude. The physico-chemical properties of the soil and vegetation types have significant effects on the AMF community structure and diversity in soil, with soil pH, WC, AP, OM and AK being the main factors affecting the AMF community.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

PY: conceptualization, data curation, methodology, software, and writing-original draft preparation. HH: data curation, methodology, software, and writing-original draft preparation. YL: investigation. HZ: investigation and experimental help. JL: investigation. LS: investigation. QX: investigation. YL: investigation. JL: investigation. XN: investigation, supervision, funding acquisition, manuscript revising, and editing. All authors contributed to the article and approved the submitted version.

Funding

This work was financially supported by the Key Project of Research and Development of Ningxia, China (Nos. 2020BFG03006 and 2021BEG02005), Project of Natural Science Foundation of Ningxia, China (No. 2020AAC03107), and Central Government Guides Local Science and Technology Development Project (No. 2022FRD05001).

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

Auge, R. M. (2001). Water relations, drought and vesicular-arbuscular mycorrhizal symbiosis. Mycorrhiza 11, 3–42. doi: 10.1007/s005720100097

CrossRef Full Text | Google Scholar

Bauer, J. T., Blumenthal, N., Miller, A. J., Ferguson, J. K., and Reynolds, H. L. (2017). Effects of between-site variation in soil microbial communities and plant-soil feedbacks on the productivity and composition of plant communities. J. Appl. Ecol. 54, 1028–1039. doi: 10.1111/1365-2664.12937

CrossRef Full Text | Google Scholar

Ceulemans, T., van Geel, M., Jacquemyn, H., Boeraeve, M., Plue, J., Saar, L., et al. (2019). Arbuscular mycorrhizal fungi in European grasslands under nutrient pollution. Glob. Ecol. Biogeogr. 28, 1796–1805. doi: 10.1111/geb.12994

CrossRef Full Text | Google Scholar

Cheng, L., Booker, F. L., Tu, C., Burkey, K. O., Zhou, L. S., Shew, H. D., et al. (2012). Arbuscular mycorrhizal fungi increase organic carbon decomposition under elevated CO₂. Science 337, 1084–1087. doi: 10.1126/science.1224304

PubMed Abstract | CrossRef Full Text | Google Scholar

Coutinho, E. S., Fernandes, G. W., Berbara, R. L. L., Valerio, H. M., and Goto, B. T. (2015). Variation of arbuscular mycorrhizal fungal communities along an altitudinal gradient in rupestrian grasslands in Brazil. Mycorrhiza 25, 627–638. doi: 10.1007/s00572-015-0636-5

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, W., Wang, J. L., Scott, M. B., Fang, Y. H., Liu, S. R., Yang, X. Y., et al. (2020). Sampling methods affect nematode-trapping fungi biodiversity patterns across an elevational gradient. BMC Microbiol. 20:15. doi: 10.1186/s12866-020-1696-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, Q., Guo, X., Chen, K. Y., Ren, S. J., Muneer, M. A., Zhang, J., et al. (2021). Phylogenetic correlation and symbiotic network explain the interdependence between plants and arbuscular mycorrhizal fungi in a Tibetan alpine meadow. Front. Plant Sci. 12:4861. doi: 10.3389/fpls.2021.804861

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, J., Jiang, Y. M., Lyu, M., Cao, C., Li, X. J., Xiong, X. L., et al. (2023). Drought changes the trade-off strategy of root and arbuscular mycorrhizal fungi growth in a subtropical Chinese fir plantation. Forests 14:15. doi: 10.3390/f14010114

CrossRef Full Text | Google Scholar

Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10:2604. 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

Elhindi, K. M., Al-Mana, F. A., El-Hendawy, S., Al-Selwey, W. A., and Elgorban, A. M. (2018). Arbuscular mycorrhizal fungi mitigates heavy metal toxicity adverse effects in sewage water contaminated soil on Tagetes erecta L. Soil Sci. Plant Nutr. 64, 662–668. doi: 10.1080/00380768.2018.1490631

CrossRef Full Text | Google Scholar

Ezeokoli, O. T., Mashigo, S. K., Maboeta, M. S., Bezuidenhout, C. C., Khasa, D. P., and Adeleke, R. A. (2020). Arbuscular mycorrhizal fungal community differentiation along a post-coal mining reclamation chronosequence in South Africa: a potential indicator of ecosystem recovery. Appl. Soil Ecol. 147:103429. doi: 10.1016/j.apsoil.2019.103429

CrossRef Full Text | Google Scholar

Ezeokoli, O. T., Nwangburuka, C. C., Adeleke, R. A., Roopnarain, A., Paterson, D. G., Maboeta, M. S., et al. (2019). Assessment of arbuscular mycorrhizal fungal spore density and viability in soil stockpiles of south African opencast coal mines. South African J. Plant Soil 36, 91–99. doi: 10.1080/02571862.2018.1537011

CrossRef Full Text | Google Scholar

Garg, N., and Chandel, S. (2011). Effect of mycorrhizal inoculation on growth, nitrogen fixation, and nutrient uptake in Cicer arietinum (L.) under salt stress. Turk. J. Agric. For. 35, 205–214. doi: 10.3906/tar-0908-12

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:115193. doi: 10.1016/j.jenvman.2022.115193

PubMed Abstract | CrossRef Full Text | Google Scholar

Haug, I., Setaro, S., and Suarez, J. P. (2019). Species composition of arbuscular mycorrhizal communities changes with elevation in the Andes of South Ecuador. PLoS One 14:e0221091. doi: 10.1371/journal.pone.0221091

PubMed Abstract | CrossRef Full Text | Google Scholar

Helgason, T., and Fitter, A. H. (2009). Natural selection and the evolutionary ecology of the arbuscular mycorrhizal fungi (phylum Glomeromycota). J. Exp. Bot. 60, 2465–2480. doi: 10.1093/jxb/erp144

PubMed Abstract | CrossRef Full Text | Google Scholar

Hodge, A., and Storer, K. (2015). Arbuscular mycorrhiza and nitrogen: implications for individual plants through to ecosystems. Plant Soil 386, 1–19. doi: 10.1007/s11104-014-2162-1

CrossRef Full Text | Google Scholar

Jamiolkowska, A., Thanoon, A. H., Skwarylo-Bednarz, B., Patkowska, E., and Mielniczuk, E. (2020). Mycorrhizal inoculation as an alternative in the ecological production of tomato (Lycopersicon esculentum mill.). Int. Agrophys. 34, 253–264. doi: 10.31545/intagr/118196

CrossRef Full Text | Google Scholar

Ji, B. M., and Bever, J. D. (2016). Plant preferential allocation and fungal reward decline with soil phosphorus: implications for mycorrhizal mutualism. Ecosphere 7:1256. doi: 10.1002/ecs2.1256

CrossRef Full Text | Google Scholar

Johnson, D., Booth, R. E., Whiteley, A. S., Bailey, M. J., Read, D. J., Grime, J. P., et al. (2003). Plant community composition affects the biomass, activity and diversity of microorganisms in limestone grassland soil. Eur. J. Soil Sci. 54, 671–678. doi: 10.1046/j.1351-0754.2003.0562.x

CrossRef Full Text | Google Scholar

Johnson, N. C., Wilson, G. W. T., Wilson, J. A., Miller, R. M., and Bowker, M. A. (2015). Mycorrhizal phenotypes and the law of the minimum. New Phytol. 205, 1473–1484. doi: 10.1111/nph.13172

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, K., Neuberger, P., Daly, E. J., Gorzelak, M., and Hernandez-Ramirez, G. (2022). Arbuscular mycorrhizal fungi community linkages to soil nutrient availability across contrasting agroecosystems. Appl. Soil Ecol. 176:104464. doi: 10.1016/j.apsoil.2022.104464

CrossRef Full Text | Google Scholar

Kruger, C., Kohout, P., Janouskova, M., Puschel, D., Frouz, J., and Rydlova, J. (2017). Plant communities rather than soil properties structure arbuscular mycorrhizal fungal communities along primary succession on a mine spoil. Front. Microbiol. 8:719. doi: 10.3389/fmicb.2017.00719

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J. B., Shen, Z. H., Li, C. N., Kou, Y. P., Wang, Y. S., Tu, B., et al. (2018). Stair-step pattern of soil bacterial diversity mainly driven by pH and vegetation types along the elevational gradients of Gongga Mountain China. Front. Microbiol. 9:569. doi: 10.3389/fmicb.2018.00569

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, M., Yue, Y. J., Wang, Z. H., Li, L., Duan, G. Z., Bai, S. L., et al. (2020). Composition of the arbuscular mycorrhizal fungal community and changes in diversity of the rhizosphere of Clematis fruticosa over three seasons across different elevations. Eur. J. Soil Sci. 71, 511–523. doi: 10.1111/ejss.12884

CrossRef Full Text | Google Scholar

Lugo, M. A., Ferrero, M., Menoyo, E., Estevez, M. C., Sineriz, F., and Anton, A. (2008). Arbuscular mycorrhizal fungi and rhizospheric bacteria diversity along an altitudinal gradient in south American puna grassland. Microb. Ecol. 55, 705–713. doi: 10.1007/s00248-007-9313-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, X. C., Geng, Q. H., Zhang, H. G., Bian, C. Y., Chen, H. Y. H., Jiang, D. L., et al. (2021). Global negative effects of nutrient enrichment on arbuscular mycorrhizal fungi, plant diversity and ecosystem multi-functionality. New Phytol. 229, 2957–2969. doi: 10.1111/nph.17077

PubMed Abstract | CrossRef Full Text | Google Scholar

Maitra, P., Zheng, Y., Wang, Y. L., Mandal, D., Lu, P. P., Gao, C., et al. (2021). Phosphorus fertilization rather than nitrogen fertilization, growing season and plant successional stage structures arbuscular mycorrhizal fungal community in a subtropical forest. Biol. Fertil. Soils 57, 685–697. doi: 10.1007/s00374-021-01554-4

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Nacoon, S., Ekprasert, J., Riddech, N., Mongkolthanaruk, W., Jogloy, S., Vorasoot, N., et al. (2021). Growth enhancement of sunchoke by arbuscular mycorrhizal fungi under drought condition. Rhizosphere 17:100308. doi: 10.1016/j.rhisph.2021.100308

CrossRef Full Text | Google Scholar

Opik, M., Moora, M., Zobel, M., Saks, U., Wheatley, R., Wright, F., et al. (2008). High diversity of arbuscular mycorrhizal fungi in a boreal herb-rich coniferous forest. New Phytol. 179, 867–876. doi: 10.1111/j.1469-8137.2008.02515.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiang, W., He, X. L., Wang, J. J., and Zhao, L. L. (2019). Temporal and spatial variation of arbuscular mycorrhizal fungi under the canopy of Hedysarum scoparium in the northern desert, China. Appl. Soil Ecol. 136, 139–147. doi: 10.1016/j.apsoil.2019.01.003

CrossRef Full Text | Google Scholar

Sato, K., Suyama, Y., Saito, M., and Sugawara, K. (2005). A new primer for discrimination of arbuscular mycorrhizal fungi with polymerase chain reaction-denature gradient gel electrophoresis. Grassl. Sci. 51, 179–181. doi: 10.1111/j.1744-697X.2005.00023.x

CrossRef Full Text | Google Scholar

Shen, C. C., Ni, Y. Y., Liang, W. J., Wang, J. J., and Chu, H. Y. (2015). Distinct soil bacterial communities along a small-scale elevational gradient in alpine tundra. Front. Microbiol. 6:582. doi: 10.3389/fmicb.2015.00582

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, Z. Y., Yin, K. J., Wang, F. Y., Mickan, B. S., Wang, X. G., Zhou, W. L., et al. (2019). Alterations of arbuscular mycorrhizal fungal diversity in soil with elevation in tropical forests of China. Div. Basel 11:10. doi: 10.3390/d11100181

CrossRef Full Text | Google Scholar

Shi, Z. Y., Zhang, J. C., Lu, S. C., Li, Y., and Wang, F. Y. (2020). Arbuscular mycorrhizal fungi improve the performance of sweet sorghum grown in a Mo-contaminated soil. J. Fungi 6:14. doi: 10.3390/jof6020044

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, J. C., Zhao, B. Y., Zheng, S., Zhang, X. W., Wang, X. L., Dong, W. T., et al. (2021). A phosphate starvation response-centered network regulates mycorrhizal symbiosis. Cells 184:30. doi: 10.1016/j.cell.2021.09.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Stackebrandt, E., and Goebel, B. M. (1994). Taxonomic note: a place for DNA-DNA re-association and 16S rRNA sequence analysis in the present species definition in bacteriology. Int. J. Syst. Bacteriol. 44, 846–849. doi: 10.1099/00207713-44-4-846

CrossRef Full Text | Google Scholar

Sun, Y., Wang, Y. B., Yang, W. J., Sun, Z., and Zhao, J. P. (2019). Variation in soil hydrological properties on shady and sunny slopes in the permafrost region Qinghai-Tibetan Plateau. Environ. Earth Sci. 78, 1–11. doi: 10.1007/s12665-019-8067-9

CrossRef Full Text | Google Scholar

Symanczik, S., Courty, P. E., Boller, T., Wiemken, A., and Al-Yahya'ei, M. N. (2015). Impact of water regimes on an experimental community of four desert arbuscular mycorrhizal fungal (AMF) species, as affected by the introduction of a non-native AMF species. Mycorrhiza 25, 639–647. doi: 10.1007/s00572-015-0638-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Vandenkoornhuyse, P., Husband, R., Daniell, T. J., Watson, I. J., Duck, J. M., Fitter, A. H., et al. (2002). Arbuscular mycorrhizal community composition associated with two plant species in a grassland ecosystem. Mol. Ecol. 11, 1555–1564. doi: 10.1046/j.1365-294X.2002.01538.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Veresoglou, S. D., Shaw, L. J., Hooker, J. E., and Sen, R. (2012). Arbuscular mycorrhizal modulation of diazotrophic and denitrifying microbial communities in the (mycor)rhizosphere of Plantago lanceolata. Soil Biol. Biochem. 53, 78–81. doi: 10.1016/j.soilbio.2012.05.007

CrossRef Full Text | Google Scholar

Vieira, L. C., da Silva, D. K. A., de Melo, M. A. C., Escobar, I. E. C., Oehl, F., and da Silva, G. A. (2019). Edaphic factors influence the distribution of arbuscular mycorrhizal fungi along an altitudinal gradient of a Tropical Mountain. Microb. Ecol. 78, 904–913. doi: 10.1007/s00248-019-01354-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X. B., and Han, Q. S. (2022). A closer examination of the “Abundant-Center” for ectomycorrhizal fungal community associated with Picea crassifolia in China. Front. Plant Sci. 13:9801. doi: 10.3389/fpls.2022.759801

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, M. Y., Pang, D. B., Chen, L., Li, X. B., Liu, L. Z., Liu, B., et al. (2021). Chemical composition of soil organic carbon and aggregate stability along an elevation gradient in Helan Mountains, Northwest China. Ecol. Indic. 131:108228. doi: 10.1016/j.ecolind.2021.108228

CrossRef Full Text | Google Scholar

Xu, X. H., Chen, C., Zhang, Z., Sun, Z. H., Chen, Y. H., Jiang, J. D., et al. (2017). The influence of environmental factors on communities of arbuscular mycorrhizal fungi associated with Chenopodium ambrosioides revealed by MiSeq sequencing investigation. Sci. Rep. 7:5134. doi: 10.1038/srep45134

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J. L., Yu, J., Wang, J. F., and Li, X. W. (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

Zanzottera, M., Fratte, M. D., Caccianiga, M., Pierce, S., and Cerabolini, B. E. L. (2020). Community-level variation in plant functional traits and ecological strategies shapes habitat structure along succession gradients in alpine environment. Community Ecol. 21, 55–65. doi: 10.1007/s42974-020-00012-9

CrossRef Full Text | Google Scholar

Zhang, M. G., Shi, Z. Y., Xu, X. F., and Wang, X. G. (2022a). Arbuscular mycorrhizal fungi associated with roots reveal high diversity levels at different elevations in tropical montane rainforests. Div. Basel 14:587. doi: 10.3390/d14080587

CrossRef Full Text | Google Scholar

Zhang, M. G., Shi, Z. Y., Yang, M., Lu, S. C., Cao, L. B., and Wang, X. G. (2021). Molecular diversity and distribution of arbuscular mycorrhizal fungi at different elevations in Mt. Taibai of Qinling Mountain. Front. Microbiol. 12:9386. doi: 10.3389/fmicb.2021.609386

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, M. G., Yang, M., Shi, Z. Y., Gao, J. K., and Wang, X. G. (2022b). Biodiversity and variations of arbuscular mycorrhizal fungi associated with roots along elevations in Mt Taibai of China. Div. Basel 14:626. 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:e4155. doi: 10.7717/peerj.4155

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, C., Long, J., Liao, H. K., Zheng, C. L., Li, J., Liu, L. F., et al. (2019). Dynamics of soil microbial communities following vegetation succession in a karst mountain ecosystem Southwest China. Sci. Rep. 9:10. doi: 10.1038/s41598-018-36886-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, L. L., Zhang, K. X., Sun, X., and He, X. L. (2022). Dynamics of arbuscular mycorrhizal fungi and glomalin in the rhizosphere of Gymnocarpos przewalskii in Northwest Desert China. Appl. Soil Ecol. 170:104251. doi: 10.1016/j.apsoil.2021.104251

CrossRef Full Text | Google Scholar