Wanqing Dong
Tingting Ding
Tingyu Duan
94% of researchers rate our articles as excellent or good
Learn more about the work of our research integrity team to safeguard the quality of each article we publish.
Find out more
- 1State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystems, Lanzhou University, Lanzhou, China
- 2Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Engineering Research Center of Grassland Industry, Ministry of Education, Gansu Tech Innovation Center of Western China Grassland Industry, Lanzhou, China
- 3College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, China
Introduction: Arbuscular mycorrhizal (AM) fungi play a crucial role in maintaining diversity and ensuring the proper functioning of grassland ecosystems. A comprehensive understanding of the diversity, distribution patterns, and drivers of AM fungi in different habitats is essential for exploring the ecological roles in grassland ecosystems.
Methods: In this study, we utilized high-throughput sequencing technology to explore the diversity of AM fungi and their distribution at an altitude of approximately 3,500 m in the alpine grassland of the eastern Tibetan Plateau. Additionally, we investigated the impacts of seasonal variation on AM fungal communities.
Results: A total of 97 species of AM fungi, comprising 937 operational taxonomic units (OTUs) belong to 9 families and 10 genera, were identified from the soil samples. Notably, the genera Glomus and Paraglomus were the most abundant and dominant within the identified communities. The composition, distribution, and diversity of the AM fungal communities in the alpine grassland of the eastern Tibetan Plateau were significantly affected by seasonal variation (p < 0.05), with geographic distance being a determining factor. Total nitrogen (TN), soil organic matter (SOM), and pH were identified as the key soil factors driving changes in AM fungal communities.
Discussion: The results demonstrated that the characteristics of AM fungal communities in the alpine grasslands of the eastern Tibetan Plateau were affected by seasonal variations and geographic location, and these findings are significant for the application of AM fungi in the restoration of grassland in similar ecosystems.
1 Introduction
Grasslands are vital ecosystems on Earth, providing a diverse range of ecological, economic, and social functions (Zhao et al., 2020). They are significant for environmental protection, biodiversity conservation, livestock economy development, and addressing global climate change (Bengtsson et al., 2019). Arbuscular mycorrhizal (AM) fungi, an important type of soil fungi, establish symbiotic relationships with the roots of numerous grassland plants (Miller et al., 2012). This symbiotic relationship not only enhances nutrient absorption for plants but also improves soil health, thereby maintaining the stability and development of ecosystems (Miller et al., 2012). AM fungi expand the contact area between plant roots and soil through the extension of hyphae. This expansion increases the absorption and transformation of essential nutrients, such as nitrogen and phosphorus, as well as water, ultimately promoting plant growth and resilience (Allen, 2007). Additionally, AM fungi can enhance the physical structure of the soil by secreting substances like glomalin, which promotes the formation of soil aggregates, improving soil moisture retention and aeration (Miller and Jastrow, 2000).
AM fungi can also significantly influence vegetation communities and ecosystems in grasslands (Zhang et al., 2012), and accelerate the establishment of plant communities in degraded grasslands (Torrez et al., 2016). During the process of grassland succession, AM fungi increase the richness and diversity of plants, particularly among late-successional species (Koziol and Bever, 2017). Moreover, AM fungi suppress the occurrence of poisonous plants by regulating the interspecific competition among grassland plants (Jin et al., 2011; Neuenkamp et al., 2019), and alleviate the negative impacts of environmental issues such as nitrogen deposition, climate change, and soil pollution (Lenoir et al., 2016; Kang et al., 2019; Cui et al., 2021). Additionally, in the restoration of degraded grasslands, the recovery of AM fungal communities is faster than plant communities, facilitating the early establishment of mycorrhizal plants and thus accelerating the grassland recovery process (Mao et al., 2019). Therefore, conducting in-depth research and providing protection for these valuable fungal resources is essential for enhancing the stability and productivity of grassland ecosystems.
The composition and diversity of microbial communities undergo significant changes with variations in environmental conditions, serving as “indicators” in ecosystems (Bissett et al., 2013). Among these communities, AM fungi are highly sensitive to environmental changes, and their community characteristics reflect the health status of grassland ecosystems. For example, the dominant species and distribution of AM fungi effectively indicate soil health and the extent of grassland degradation (Vasconcellos et al., 2016; Wu et al., 2021). This sensitivity has increased the value of AM fungi in environmental monitoring and ecological restoration efforts (Asmelash et al., 2016). The composition, diversity, and functional characteristics of AM fungal communities are influenced by various factors, including vegetation communities, soil physicochemical properties, climate change, and human activities (Mohammad et al., 2003; Jamiołkowska et al., 2018; Jerbi et al., 2021). Among these factors, the type of host plant is a key biological determinant affecting AM fungal communities (López-García et al., 2017). Additionally, other factors, such as soil pH, directly influence the growth, reproduction, and interactions of AM fungi with plants (Van Aarle et al., 2002; Li et al., 2021). Moreover, variations in soil nutrient content, particularly nitrogen and phosphorus, due to management practices or climate change significantly impact the structure and function of AM fungal communities (Treseder and Allen, 2002). Research indicates that soil available nitrogen can directly or indirectly alter the composition of AM fungi (Xiao et al., 2021; Shi et al., 2022), while soil available phosphorus shows a significant negative correlation with the richness of AM fungi (Ceulemans et al., 2019).
Recognizing the interrelated impacts of various environmental and anthropogenic factors on AM fungal communities is crucial, as these factors are often interconnected rather than isolated. For example, significant differences in the composition and diversity of AM fungal communities have been observed across altitudinal gradients. As altitude increases, changes in climate and soil conditions typically occur, directly impacting the distribution and activity of AM fungi. Additionally, these changes can produce indirect effects through alterations in vegetation communities (Gai et al., 2012; Geml, 2017). External conditions such as mowing, grazing, fertilization, and climate change have complex and patterned influences on AM fungal communities (Chen et al., 2014; Faghihinia et al., 2022; Zubek et al., 2022). Therefore, ongoing in-depth research on the responses of AM fungal communities to environmental changes will enhance our understanding of the adaptive mechanisms and future trends of grassland ecosystems.
The Tibetan Plateau, recognized as the highest geomorphological unit on Earth, is characterized by its unique geography, alpine climate, and natural geographic distribution pattern, which play a vital role in maintaining ecological stability (Dong et al., 2022). Alpine grasslands encompass over three-fifths of the plateau and mainly include alpine meadow and alpine steppe, which are the two most representative types of alpine grassland ecosystems (Duan et al., 2021). Over the past decades, the Tibetan Plateau has experienced the impacts of climate change (Latif et al., 2019), rodent foraging (Zhang et al., 2022) and overgrazing (Niu et al., 2019). These factors have led to a rapid decline in the diversity, structure, and function of grassland vegetation communities (Zhang et al., 2019). Additionally, there has been a loss of soil nutrients and degradation of soil structure, which hinders plant growth and reduces above-ground productivity (Zhang et al., 2019). The highly fragile nature of the alpine grassland ecosystem on the Tibetan Plateau makes restoration difficult once degradation has occurred (Dong, 2023). Therefore, research on restoration of degraded grasslands has gradually become a research hotspot.
At present, most studies on the diversity of AM fungi have been conducted in regions such as Tibet and Qinghai in the southwest and northwest. In contrast, relatively few investigations have been conducted in eastern areas such as Sichuan and Gansu (Zhang et al., 2017). Given the vast extent of the Tibetan Plateau and the significant habitat differences between different regions, further research is needed to expand the scope of research. Moreover, a multitude of studies have revealed that seasonal variation can affect spore density, root colonization, community composition, and diversity of AM fungi (Li et al., 2005; Lara-Pérez et al., 2020; Liu et al., 2020). However, little attention has been given to how seasonal variation affects the diversity and composition of AM fungal communities specifically on the Tibetan Plateau. In light of this gap, the current study used high-throughput sequencing technology to explore the composition of AM fungal communities and their seasonal dynamics in alpine grasslands of the eastern Tibetan Plateau. Additionally, the relationships between soil properties, seasonal variation, and AM fungal communities were analyzed. We hypothesized that seasonal variations modulate soil properties, thereby indirectly affecting the structure and composition of AM fungal communities.
2 Materials and methods
2.1 Study region
The research sites are located in the Aba Tibetan and Qiang Autonomous Prefecture, situated on the southeastern edge of the Tibetan Plateau, and constitutes one of the principal distribution regions of natural grasslands in Sichuan Province. The main forms of grassland degradation at the sites include desertification and the emergence of black soil patches. The present study selected 2 distinct types of degraded grasslands in Aba Prefecture in the eastern Tibetan Plateau. Site 1 is located in Hongyuan county (HY, 32°53′7″ - 32°53′25″ N, 102°35′58″ - 102°36′2″ E) and is mainly afflicted by desertification. The dominant plants are Kobresia setchwanensis and K. pygmaea, with accompanying species such as Potentilla anserina, Aster alpina, Anemone trullifolia, A. rivularis, and Thalictrum alpinum. Site 2 is located in Ruoergai county (R, 33°55′40″-33°55′43″ N, 102°23′30″-102°23′33″ E) and is mainly characterized by bared black soil patches. The dominant plants are Carex moorcroftii and K. pygrwaea, with subsidiary species such as Elymus nutans, Leymus paboanus and Artemisia wellbyi. Both sites exhibit a moderate degree of degradation and are situated at an elevation of approximately 3,500 m. The mean annual temperature ranges from 1.1 to 2.9°C, while the average annual precipitation varies from 648 to 860 mm, indicative of a continental plateau cold temperate monsoon climate.
2.2 Sample collection
To evaluate the AM fungal community and its seasonal response in the study area, soil samples were collected in July (Summer, S) and September (Autumn, A) in the year 2020 and 2021 from two sites: Hongyuan (HY) and Ruoergai (RA) in the alpine grassland on the Tibetan Plateau. The eight sampling events are designated as follows: HYS20 for Hongyuan Summer 2020, RS20 for Ruoergai Summer 2020, HYA20 for Hongyuan Autumn 2020, RA20 for Ruoergai Autumn 2020, HYS21 for Hongyuan Summer 2021, RS21 for Ruoergai Summer 2021, HYA21 for Hongyuan Autumn 2021, and RA21 for Ruoergai Autumn 2021. Within each site, 12 sampling points were randomly selected, each at least 100 m apart, with each sampling point considered a replicate. At each sampling point, 3 soil samples (0–15 cm) were randomly collected using a shovel and thoroughly mixed to form a composite sample. A total of 96 soil samples were collected (2 sites ×2 seasons ×2 years ×12 replicates). Fresh soil samples had visible plant residues, roots, and stones removed and were then passed through a sieve (< 2 mm). The sieved soil samples were split into two parts. One subsample was used for sequencing, while the other was air-dried to evaluate the basic properties of the soil.
2.3 Soil property measurements
Refer to Lu (2000) for the determination of soil properties. Soil pH was measured using a pH meter (1∶2.5, Vsoil ∶ Vwater) (PHS-3C, Shanghai Yidian, China). Total nitrogen (TN) content in soil was determined through acid digestion and the Kjeldahl method using a decoction oven (SH220F, Shandong Haineng, China) and an Automated Discrete Analyzer (SMARTCHEM 450, Beijing Lijia, China). Total phosphorus (TP) content in soil was dissolved in a mixture of concentrated H2SO4 and HClO4 using a decoction oven (SH220F, Shandong Haineng, China) and assessed according to the molybdate-blue colorimetric method in an Automated Discrete Analyzer (SMARTCHEM 450, Beijing Lijia, China). Available potassium (AK) content in soil was extracted with NH4OAc and measured using a Flame Photometer (Model 410 Flame Photometer, Sherwood Scientific, UK). Soil organic matter (SOM) content was determined by the potassium dichromate oxidation method using an oil bath (HH-S, Changzhou Putian, China).
2.4 Illumina HiSeq sequencing and bioinformatics analyses
Illumina HiSeq sequencing was employed to investigate the community diversity and distribution of AM fungi. Firstly, soil DNA was extracted from 0.5 g of fresh soil using the E.Z.N.A.® soil DNA kit (Omega Biotek, United States) following the manufacturer’s instructions. Purity and concentration were assessed using a NanoDrop 2000 ultra-micro spectrophotometer (Thermo Scientific, Inc., United States). The V4-V5 hypervariable regions of the fungal gene (18S rRNA) were amplified with primers AML1 and AML2 in a thermal cycle PCR system, followed by a second amplification using the primers AMV4.5NF and AMDGR (Lee et al., 2008; Suzuki et al., 2020). In the PCRs, a mixture consisting (20 μL) of 5 × FastPfu Buffer (4 μL), 2.5 mM dNTPs (2 μL), 0.5 μM each primer (0.8 μL), FastPfu Polymerase (0.4 μL), BSA (0.2 μL) and Template DNA (10 ng), was utilized for the reactions. The standard procedure for database building operations was conducted according to the NEBNext® Ultra™ DNA Library Prep Kit for Illumina® (New England Biolabs, United States). Amplicon sequencing was performed using Illumina Hiseq2500 platform (Meige Company, Shenzhen, China).
Finally, the sequencing data were processed. The quality of the Raw Reads data was filtered at both ends independently using the Trimmomatic software (V0.33). Leveraging the primer and barcode information at the start and end of the sequences, the sequences were allocated to the corresponding samples with the help of Mothur software (V1.35.1). Subsequently, the barcodes and primers were removed, thereby obtaining paired-end clean reads after quality control. Next, the PE reads were spliced by means of FLASH (V1.2.11) software and filtered to acquire the raw spliced sequences (raw tags). The spliced sequences were subject to quality control and filtration using Mothur software to obtain valid spliced fragments (Clean Tags). OTUs clustering was then performed and 18S reads were purified from the MaarjAM database (Opik et al., 2010). Subsequently, the OTU was assigned to different fungal taxa (BLAST), with a 97% percentage similarity (Altschul et al., 1990). The sequence data has been uploaded to the NCBI database (BioProject ID: PRJNA1173197).
2.5 Statistical analyses
Based on the Richness-rarefaction curve, the sequencing depth in this study is sufficient (Supplementary Figure S1). We compared the number of OTUs and community composition at different sites and their seasonal dynamics using the original OTU data. Subsequently, based on the minimum sequence, we rarefied the original OTU table, and the resulting data were utilized to analyze the Alpha diversity, Beta diversity, differential microbial communities, and correlation analysis of AM fungi. R software (V2.15.3) was utilized to conduct statistics on the relative abundances of OTUs, followed by the creation of a bar chart depicting community composition. Subsequently, QIIME software (V1.9.1) was employed to calculate four diversity indices (Chao1, Richness, Shannon, Simpson). To evaluate effects of season, site, their interactions, and year on the composition of the AM fungal community, we utilized “strata” in the “adonis” function to account for difference in site and season. Furthermore, to illustrate the difference in AM fungal community composition for site × season, Principal Coordinate Analysis (PCoA) was employed to order the Bray–Curtis dissimilarity matrix.
For the difference analysis of species between samples, we used LEfSe. First, the taxa with significant abundance difference between samples were compared and screened. Next, the paired Wilcoxon rank-sum test was employed to assess the differences between each pair of samples. Finally, Linear Discriminant Analysis (LDA) was utilized for dimensionality reduction and to evaluate the impact magnitude of taxa with significant differences (LDA > 2), followed by plotting. Redundancy Analysis (RDA) was conducted to explore the relationship between AM fungal community characteristics and soil properties. Pearson’s correlations analysis was employed to investigate the relationships between AM fungi and soil properties. The Co-occurrence network diagrams were calculated and plotted using R (version 4.2.3) and Gephi (version 0.9.2) (threshold = 0.0001). Bar graphs were plotted in GraphPad Prism 8.0, while Venn diagrams, PCoA, RDA, and Heatmaps were generated using Wekemo Bioincloud1 (Gao et al., 2024). One-way Analysis of Variance (ANOVA) was used to assess soil properties and AM fungal diversity (SPSS, version 22.0). At the 5% significance level, the Tukey (HSD) test was employed for comparative analysis of the significance of different samples.
3 Results
3.1 Soil properties
There were significant differences in soil properties between sites and seasons as measured over 2 years (Table 1). In 2020, the contents of TN, TP and SOM in the soil ranged from 5.01–14.10, 0.82–1.00, and 8.56–18.09 g/kg, respectively, while the contents of AK in the soil ranged from 81.35–121.55 mg/kg. Seasons (summer vs. autumn) did not affect TP or AK, however, TN and SOM exhibited significant differences. The TN and SOM in summer soils were significantly higher those in autumn soils (p < 0.05). Moreover, site location (Hongyuan vs. Ruoergai) had no effect on TP, but TN, SOM and AK showed significant differences. The TN in Hongyuan was significantly lower than in Ruoergai in summer soils (p < 0.05), while no significant difference was observed in autumn soils. Additionally, SOM and AK in Hongyuan were significantly lower than those in Ruoergai across all soil samples (p < 0.05). The soil pH ranged from 4.48 to 5.41 in Hongyuan and from 6.48 to 7.49 in Ruoergai. Consistent with TN, the pH of summer soil samples was significantly higher than that of autumn soil samples (p < 0.05).
Table 1
Table 1. AM fungi alpha-diversity indices in all soil samples.
In 2021, the contents of TN, TP, and SOM in soil samples ranged from 2.54–6.21, 1.01–1.21, and 4.62–13.82 g/kg, respectively, while AK contents ranged from 69.49–111.53 mg/kg. In Hongyuan, the season significantly affected TN, TP, and SOM content, however AK showed no significant differences. In Ruoergai, the season had no effect on TN, TP, SOM, and AK. Additionally, TN and SOM were significantly higher in Hongyuan than in Ruoergai for both summer and autumn soils (p < 0.05). For TP and AK, levels were significantly higher in HYS than in RS (p < 0.05), with no significant difference between HYA and RA. The soil pH ranged from 5.71–6.04 in Hongyuan and from 7.42–7.61 in Ruoergai. In Hongyuan, the pH of summer samples was significantly higher than that of autumn samples (p < 0.05), while no significant difference was observed in Ruoergai. These results showed that seasonal variation significantly affected the physicochemical properties of the soils in Hongyuan and Ruoergai, particularly TN, TP, SOM and pH.
3.2 AM fungal community composition and structure
We characterized 15,933 OTUs (operational taxonomic units) of AM fungi using a 97% sequence similarity cutoff from 449,063 to 66,657 high-quality reads of the 96 soil samples. Most OTUs of AM fungi were classified as Glomus (15.73%), Paraglomus (15.60%), Claroideoglomus (5.06%), Diversispora (2.49%), followed by Entrophospora (2.15%), Ambispora (2.71%), and Others (1.52%) at the genus level (Supplementary Table S1). A total of 937 OTUs were detected, belonging to 4 order, 9 families, 10 genera and 97 species (Supplementary Table S2). Of the 97 species, 50 belonged to the genus Glomus, accounting for 51.55% of the total, making it the dominant genus. The remaining genera, in order from highest to lowest, are Acaulospora, Claroideoglomus, Archaeospora, Diversispora, Scutellospora, Paraglomus, Ambispora, Entrophospora, and Pacispora, which account for 11.34, 11.34, 5.15, 5.15, 5.15, 5.15, 3.09, 1.03, and 1.03% of the total, respectively (Supplementary Figure S2).
The number of unique and shared OTUs varied across different seasons and sites in 2020 and 2021 (Figure 1). In 2020, we identified 1,453 shared OTUs (Figure 1A), while in 2021, this number decreased to 992 (Figure 1B). The trend for unique OTUs was consistent in both years: RA > HYA > HYS > RS. Subsequently, we analyzed the number of AM fungal species in different soil samples (Supplementary Figure S3). The number of AM fungal species showed a gradual increase over time in Hongyuan, whereas no similar variation was observed in Ruoergai. In addition, in 2020, the number of AM fungal species in autumn soils from Hongyuan was significantly higher than in summer samples (HYA20 > HYS20) (p < 0.05), but no significant differences were found in 2021. In Ruoergai, there was no significant variation in the number of AM fungal species across different years and seasons. Thus, the effect of seasonal variation on the soil AM fungal community appears to differ among sites.
Figure 1
Figure 1. Venn. Number of unique and shared OTUs across different sites and seasons in 2020 (A) and 2021 (B). HYS20, HYA20, RS20, RA20, HYS21, HYA21, RS21, and RA21 represent different groups of soil samples.
To illustrate difference in the AM fungal community for site × season, we applied a PCoA using the Bray-Curtis algorithm (Figure 2). The PCoA assigned the AM fungal community to two groups: (HYS and HYA) and (RS and RA) in 2020 (Figure 2A). The HYS20 and HYA20 were clustered together and were distant from RS20 and RA20. In 2021, the PCoA assigned the AM fungal community to four groups: HYS, HYA, RS and RA (Figure 2B), which were distant from one another. Thus, the AM fungal community structure significantly changed with the season in 2020 and significantly changed with both season and site in 2021. PERMANOVA revealed that the interaction of season and site significantly affected the AM fungal community in both 2020 and 2021 (p = 0.001).
Figure 2
Figure 2. The principal coordinates analysis (PCoA) of AM fungal communities from soil samples collected in 2020 (A) and 2021 (B) is based on the Bray-Curtis distances. HYS20, HYA20, RS20, RA20, HYS21, HYA21, RS21, and RA21 represent different groups of soil samples.
Additionally, we analyzed community composition at the family and genus levels based on the relative abundance of AM fungi in soil samples (Figure 3). In 2020, the relative abundance of Paraglomeraceae (Paraglomus) was higher, whereas that of Ambisporaceae (Ambispora) and Claroideoglomeraceae (Claroideoglomus) were lower in HYA20 compared to HYS20. The relative abundance of Glomeraceae (Glomus) and Paraglomeraceae (Paraglomus) were higher, while that of Ambisporaceae (Ambispora), Claroideoglomeraceae (Claroideoglomus), and Diversisporaceae (Diversispora) were lower in RA20 than in RS20. In 2021, the relative abundance of Ambisporaceae (Ambispora), Claroideoglomeraceae (Claroideoglomus), and Glomeraceae (Glomus) was higher, whereas that of Paraglomeraceae (Paraglomus) were lower in HYA21 than HYS21. The relative abundance of Ambisporaceae (Ambispora), Claroideoglomeraceae (Claroideoglomus), Diversisporaceae (Diversispora), and Glomeraceae (Glomus) was higher in RA21 than in RS21. This result was contrary to that of 2020. These findings indicate that the community composition of AM fungi in the soils of Hongyuan and Ruoergai varies by years and is significantly influenced by seasonal variation.
Figure 3
Figure 3. Community composition of AM fungi. The relative abundance (%) of AM fungi (based on the number of OTUs) at the family and genus levels in 2020 (A, B) and 2021 (C, D). HYS20, HYA20, RS20, RA20, HYS21, HYA21, RS21, and RA21 represent different groups of soil samples. Others, Taxa with an abundance of <1% are summarized as “others,” including Archaeosporaceae (Archaeospora), Acaulosporaceae (Acaulospora), Gigasporaceae (Scutellospora), and Pacisporaceae (Pacispora).
To identify the key differential species of AM fungi in different soil samples, we conducted a LEfSe analysis (Figure 4; Supplementary Figure S4). The linear discriminant analysis (LDA) of LEfSe across all soil samples indicated that more AM fungal taxa (4 families and 4 genus) were detected in the HYS20 compared to the other three subgroups in 2020 (Figure 4A). Additionally, the summer soil samples exhibited a greater number of differential AM fungal taxa than those in autumn (HYS20 > HYA20; RS20 > RA20). In 2021, the LEfSe analysis revealed that more AM fungal taxa (3 families and 3 genus) were detected in the HYA21 than in the other three subgroups (Figure 4B). Furthermore, the autumn soil samples in Hongyuan contained more differential AM fungal taxa than the summer samples, while the reverse was true in Ruoergai. We further demonstrated that seasonal variation significantly affected the community composition of AM fungi in both Hongyuan and Ruoergai, particularly in the latter, with effects varying between years.
Figure 4
Figure 4. LEfSe analysis of AM fungal communities in 2020 (A) and 2021 (B) shows LDA scores exceeding 2.0 across all soil samples. (C–F) Represent the AM fungal taxa that exhibited significant changes due to seasonal variations across different years and sites. The LDA bars represent AM fungal taxa (at the family and genus levels) with LDA scores greater than 2.0 in soil samples from different sites and seasons (p < 0.05). HYS20, HYA20, RS20, RA20, HYS21, HYA21, RS21, and RA21 represent different groups of soil samples.
Subsequently, we analyzed the AM fungal taxa in Hongyuan and Ruoergai that exhibited significantly changes due to seasonal influences across different years (Figures 4C–F). In Hongyuan, there were 2 and 13 AM fungal taxa with significant changes in 2020 and 2021, respectively. Specifically, the AM fungal taxa that changed significantly under seasonal influence included 1 family and 1 genus: Claroideoglomeraceae (Claroideoglomus) in 2020 (Figure 4C). In 2021, the differential AM fungal taxa included 6 families and 7 genera, namely Archaeosporaceae (Archaeospora), Ambisporaceae (Ambispora), Diversisporaceae (Diversispora), Entrophospora, Gigasporaceae (Scutellospora), Glomeraceae (Glomus), and Paraglomeraceae (Paraglomus) (Figure 4D). In Ruoergai, there were 13 AM fungal taxa with significant changes in 2020 and 2021, respectively. The differential AM fungal taxa included 6 families and 7 genera, namely Acaulosporaceae (Acaulospora), Ambisporaceae (Ambispora), Claroideoglomeraceae (Claroideoglomus), Diversisporaceae (Diversispora), Entrophospora, Glomeraceae (Glomus), and Pacisporaceae (Pacispora) in 2020 (Figure 4E). The differential AM fungal taxa consisted of 6 families and 7 genera: Acaulosporaceae (Acaulospora), Claroideoglomeraceae (Claroideoglomus), Diversisporaceae (Diversispora), Entrophospora, Glomeraceae (Glomus), Paraglomeraceae (Paraglomus), and Pacisporaceae (Pacispora) in 2021 (Figure 4F). This result serves as an important data support for subsequent analysis and research.
3.3 AM fungal richness and diversity
We calculated the Chao1, Richness, Shannon and Simpson indices of AM fungi to evaluate the α-diversity of the AM fungal community in 2020 and 2021 across different sites and seasons (Table 2). The results showed that the Chao1, Richness, Shannon and Simpson indices of AM fungal communities in soil samples varied significantly between sites and seasons (p < 0.01), although the trends differed between years.
Table 2
Table 2. Comparison of soil properties under different sites and seasons.
In 2020, the Chao1 and Richness indices of AM fungal communities in autumn soils were similar to those in summer soils. The Shannon index was higher in autumn soils compared to summer soils, with significant differences observed in Ruoergai (p < 0.05). Conversely, the Simpson index in autumn soil samples was significantly lower than in summer samples (p < 0.05). Additionally, the Chao1, Richness and Shannon were significantly higher in RA20 than in HYA20 (p < 0.05), while Shannon index was significantly lower in RS20 than in HYS20 (p < 0.05).
In 2021, the Chao1 and Richness indices of AM fungal communities in autumn soils were significantly higher than in summer soils (HYA21 > HYS21; RA21 > RS21) (p < 0.05). The Shannon index in autumn soils was similar to those in summer soils in Hongyuan. The Simpson index was higher in HYS21 than in HYA21. The Chao1, Richness, and Shannon were significantly higher in RA21 than in HYA21 (p < 0.05), Shannon was significantly higher in RS21 than in HYS21 (p < 0.05), and Simpson was significantly lower in RS21 than in HYS21 (p < 0.05). These findings indicate potential site-specific environmental factors influencing AM fungal communities. Furthermore, the significant differences in diversity indices within the same season across different sites suggest that local conditions, such as soil properties and microclimates, play a crucial role in shaping AM fungal diversity.
3.4 AM fungal network structure
The network structure based on the OTU data of AM fungal revealed differences between sites and seasons, with the AM fungal in Ruoergai exhibiting a significantly more complex network than those in Hongyuan (Figure 5). We further analyzed and compared the number of nodes, edges, weighted average degree, modularity, average clustering coefficient, and average path length across eight network structures (Figures 5I–N). The results indicated that the number of nodes, edges, weighted average degree, average clustering coefficient, and average path length were higher in HYS20 than in HYA20, while these properties were lower in HYS21 than in HYA21. The network modularization was lower in HYS20 than in HYA20, but higher in HYS21 than in HYA21. In Ruoergai, the number of nodes, edges, weighted average degree, and average path length were lower in summer soils than in autumn soils (RS20 < RA20; RS21 < RA21). Conversely, the network modularization was higher in summer soils than in autumn soils (RS20 > RA20; RS21 > RA21). The average clustering coefficient was higher in RS20 compared to RA20, while it was lower in RS21 than in RA21. These findings suggest that seasonal variation influences the complexity of AM fungal communities in the soils of Hongyuan and Ruoergai.
Figure 5
Figure 5. Network structure analysis of AM fungal community. (A–H) Represent the co-occurrence pattern of HYS20, HYA20, HYS21, HYA21, RS20, RA20, RS21, and RA21, respectively. Nodes are colored according to AM fungal taxa at the genus level, while the width and color of edges indicate the correlation coefficient. Different colors represent different AM fungal genera. (I–N) Represent different network topology properties, including the number of nodes, number of edges, weighted average degree, modularity, average clustering coefficient, and average path length.
3.5 Correlations of AM fungal community with soil properties
The relationship between AM fungal communities and soil properties was analyzed at the genus level (Figure 6). A redundancy analysis (RDA) was conducted to assess the influences of soil factors on the AM fungal community across all soil samples. The first two ordinal axes of the RDA accounted for 72.02 and 20.42% of the total variation in the soil AM fungal community from 2020, respectively (Figure 6A). The first axis was negatively correlated with soil TN, AK, SOM, and pH, while showing a positive correlation with soil TP. The second axis exhibited a positive correlation with soil TN, TP, and SOM, and a negative correlation with soil AK and pH. The most significant soil factors were pH (R2 = 0.77, p < 0.001) and TN (R2 = 0.48, p < 0.001) (Figure 6A; Supplementary Table S3).
Figure 6
Figure 6. Redundancy analysis (RDA) at the genus level and Spearman’s rank correlation heatmap (p < 0.05 *, p < 0.01 **, p < 0.001 ***), were employed to examine the correlation between AM fungal communities and soil properties in 2020 (A,B) and 2021 (C,D). The genera with abundance <1% are summarized as “others.” HYS20, HYA20, RS20, RA20, HYS21, HYA21, RS21, and RA21 represent different groups of soil samples. TN, total nitrogen; TP, total phosphorus; AK, available potassium; SOM, soil organic matter; pH, soil pH.
We also employed Spearman’s rank correlation to evaluate the relationships among the different AM fungal genera and soil properties in 2020 (Figure 6B; Supplementary Table S4). Acaulospora was significantly negatively correlated with pH (R2 = −0.36, p = 0.012). Ambispora was significantly negatively correlated with SOM (R2 = −0.35, p = 0.014), AK (R2 = −0.38, p = 0.008), and pH (R2 = −0.46, p = 0.001). Claroideoglomus, Diversispora, and Entrophospora were significantly positively correlated with TN, SOM, and pH (p < 0.05). In contrast, Glomus and Paralomus were significantly negatively correlated with TN, SOM, and pH (p < 0.05). Additionally, Claroideoglomus and Entrophospora were significantly positively correlated with AK, while Paralomus was significantly negatively correlated (p < 0.05). Others were significantly positively correlated with soil TN (R2 = 0.34, p = 0.020), SOM (R2 = 0.39, p = 0.007), and pH (R2 = 0.55, p < 0.001), while negatively correlated with TP (R2 = −0.33, p = 0.023).
The first two ordinal axes of RDA accounted for 59.80 and 34.14% of the total variation in the soil AM fungal community from 2021, respectively (Figure 6C). The first axis was positively correlated with soil TP and AK, but negatively correlated with soil TN, SOM, and pH. The second axis was positively correlated with soil TP, AK, SOM, and TN, while negatively correlated with soil pH. The most significant soil factors were TN (R2 = 0.79, p < 0.001), pH (R2 = 0.77, p < 0.001), and SOM (R2 = 0.73, p < 0.001) (Figure 6C; Supplementary Table S3). Spearman’s rank correlation was employed to assess the relationship among the abundant AM fungal genera and soil properties in 2021 (Figure 6D; Supplementary Table S4). Overall, the correlation between AM fungi (Acaulospora, Ambispora, Claroideoglomus, Diversispora, and Entrophospora) and TN, TP, SOM, and pH in 2021 was similar to that in 2020, but with stronger correlations observed in 2021. Glomus was significantly negatively correlated with TP (p < 0.05), and Paraglomus was positively correlated with SOM and AK, which was exactly the opposite of the results in 2020.
Furthermore, we developed a correlation coefficient matrix to examine the relationship between soil factors and the diversity and richness of AM fungal communities (Figure 7). Soil TN was significantly positively correlated with the Simpson index of the AM fungal community, whereas it was significantly negatively correlated with the Shannon index in 2020 (p < 0.05). Soil TN, TP, SOM, and pH were significantly positively correlated with the Chao1 and Richness indices, while TP and AK were significantly negatively correlated with the Chao1 and Richness indices (p < 0.05). Additionally, soil TP was significantly positively correlated with the Simpson index, while it was significantly negatively correlated with the Shannon index (p < 0.05).
Figure 7
Figure 7. Correlation coefficient matrix. Factors influencing the AM fungal community were analyzed using Pearson correlation (p < 0.05 *, p < 0.01 **, p < 0.001 ***). Red colors indicate negative correlations blue colors indicate positive connections. Darker color represents a greater correlation coefficient. (A) 2020, (B) 2021. TN, total nitrogen; TP, total phosphorus; AK, available potassium; SOM, soil organic matter; pH, soil pH.
4 Discussion
Extreme habitats often host groups of organisms with unique physiological and ecological functions. The Tibetan plateau alpine grassland ecosystem is a typical extreme habitat, making the study of AM fungal diversity particularly significant. In this research, we analyzed the characteristics and seasonal dynamics of soil AM fungal communities in Hongyuan and Ruoergai. We determined the diversity of AM fungal communities and their response mechanisms to seasonal variation in different degraded grasslands on the eastern Tibetan Plateau. The findings suggest that both spatial (location-specific) and temporal (interannual and seasonal) factors significantly influence the diversity and composition of AM fungal communities.
4.1 The AM fungal species
A total of 97 species of AM fungi were identified in this study, representing approximately 66.90% of the known species of AM fungi in China, indicating that the soil of Hongyuan and Ruoergai are abundant in AM fungal resources. The 97 AM fungi identified belong to 10 genera: Acaulospora, Archaeospora, Claroideoglomus, Diversispora, Entrophospora, Gigaspora, Glomus, Pacispora, Paraglomus, and Scutellospora, all of which have been previously reported in the Tibetan Plateau (Gai et al., 2009; Mao et al., 2024; Zhang et al., 2024). Among these, Glomus and Paraglomus were the dominant genera, detected in all soil samples, with relative abundances of 15.73 and 15.60%, respectively. The spatial distribution of AM fungal communities in 15 alpine grasslands on the Tibetan Plateau showed that Glomus was the dominant genus in alpine meadow and alpine steppe, followed by Paraglomus (Zhang et al., 2024). Glomus has repeatedly been shown to be widely distributed in different ecosystems worldwide and to be highly adaptable to extreme environments (Coutinho et al., 2015; Bonfim et al., 2016). Additionally, Acaulospora occurred at different sites and seasons, although its relative abundance was below 1% in this study. It has been found to have a higher occurrence frequency in alpine grasslands (Gai et al., 2009), and its relative abundance gradually increases with elevation (3105–4,556 m) (Li et al., 2018). This result further confirms that Acaulospora is more suited to survive in high-altitude alpine grasslands.
We also found that the AM fungal taxa exhibiting significant changes mainly belonged to Paraglomeraceae (Paraglomus), Gigasporaceae (Scutellospora), Glomeraceae (Glomus), Acaulosporaceae (Acaulospora), and Ambisporaceae (Ambispora) in Hongyuan, as well as Diversisporaceae (Diversispora), Claroideoglomeraceae (Claroideoglomus), Pacisporaceae (Pacispora), Glomeraceae (Glomus), Entrophospora, and Archaeosporaceae (Archaeospora) in Ruoergai (Figures 4A,B). This indicates that these AM fungal taxa are highly sensitive to changes in the external environment. Future research should focus on these AM fungal taxa, in particularly Glomus and Paraglomus.
4.2 Changes in AM fungal community characteristics
Seasonal variation is an important factor influencing AM fungal communities, but its role in the alpine grassland ecosystem on the Tibetan Plateau is still lacking. In this study, we found seasonal variation affected the species, community structure, composition, diversity, and network complexity of AM fungi in Hongyuan and Ruoergai. The results are consistent with previous findings in other grassland ecosystems. For example, Dumbrell et al. (2011) found that the composition and structure of AM fungal communities in temperate grassland significantly changed with seasonal variations. Su et al. (2011) explored the seasonal variation of AM fungi in the rhizosphere of five plant species in Inner Mongolia grassland, and found that seasons had significant effects on the spore density, Richness, and Shannon index of AM fungi. Our hypothesis that seasonal variations are a crucial factor affecting the AM fungal community in the alpine grasslands of the eastern Tibetan Plateau, and this viewpoint is supported. Previous studies have shown that the response of AM fungal communities to seasonal variation differs and is affected by plant species and habitat (Helgason et al., 2014). In this study, we found that the response of AM fungal communities to seasonal variation is also affected by site-specific environmental factors, as evidenced by the contrasting results of Shannon index and the network complexity of AM fungal communities in Hongyuan soils and Ruoergai soils in 2021. Therefore, we speculate that, in addition to seasonal variations, geographical distance is also an important factor influencing the AM fungal communities in the eastern Tibetan Plateau. Previous studies have indicated that geographic distance significantly affects the community dynamics of AM fungi in the Tibetan Plateau (Zhang et al., 2024), and has a more significant impact on fungal communities compared to other environmental factors (Kivlin et al., 2011; Beck et al., 2015). This phenomenon is also known as the distance decay rule (Wang et al., 2022).
However, some studies have found the community composition, richness, and spore density of AM fungi do not change significantly with seasonal variations (Pereira et al., 2018; Álvarez-Lopeztello et al., 2019). Thus, the effect of seasonal variation on the characteristics of AM fungal communities is complex and influenced by a various of environmental factors. Therefore, further research is necessary to clarify the impact of climate and environmental factors on the characterization of AM fungal communities in these ecosystems.
4.3 Correlations of changes in AM fungal community with soil properties
Soil, as an important medium for the interaction between plants and AM fungi, is usually regarded as a key environmental factor (Davison et al., 2015; Wang et al., 2021). In this study, seasonal variations significantly impacted the TN, TP, SOM, and pH of the soils in Hongyuan and Ruoergai. The changes in these indicators were significantly correlated with the variations in AM fungal communities, and similar patterns of change were observed across different years. A comprehensive analysis of the results indicated that soil TN, SOM, and pH were the key soil factors influencing the characteristics of AM fungal communities. In previous studies, TN and pH have also been identified as crucial soil factors influencing the structure of AM fungal communities (Guo and Gong, 2014; Ma et al., 2021). Additionally, the pH can affect the root colonization of AM fungi and determine their niche space (Van Aarle et al., 2002; Davison et al., 2021). The compounds produced by the decomposition of SOM contribute to the colonization and hyphal growth of AM fungi (Gryndler et al., 2009). Moreover, AM fungi can, in turn, inhibit the decomposition of SOM by competing with saprophytic fungi (Frey, 2019), or increase soil TN content by reducing nitrogen leaching (Fang et al., 2020). Therefore, we propose that the response of soil AM fungal community characteristics to seasonal variation in the eastern Tibetan Plateau is linked to soil properties.
5 Conclusion
Through high-throughput sequencing, the present study reveals the diversity, seasonal dynamics, and soil driving factors of AM fungal communities in the eastern Tibetan Plateau. A total of 937 OTUs were detected, belonging to 4 orders, 9 families, 10 genera, and 97 species, suggesting that the research sites possess a high abundance and rich diversity of AM fungi. Seasonal variation exerts a significantly impact on the community structure, composition, diversity, and network complexity of AM fungi. Moreover, the influence of this seasonal variation is modulated by geographic distance and year. Additionally, soil properties also play a vital role as factors influencing AM fungal communities. Among these properties, total nitrogen, soil organic matter and pH are the key factors. The dominant AM fungi genera, such as Glomus and Paraglomus, are highly sensitive to seasonal variation and thus deserve particular attention in ecological monitoring as well as grassland restoration efforts.
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 at: https://www.ncbi.nlm.nih.gov/, PRJNA1173197.
Author contributions
WD: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. TDi: Writing – original draft, Writing – review & editing. TDu: Conceptualization, Investigation, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was supported by the National Key Research and Development Program of China (2019YFC0507703).
Acknowledgments
We thank Jie Deng, Yingde Li, and Faxi Li for their help in field sampling and lab analysis.
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.
Generative AI statement
The authors declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1511979/full#supplementary-material
Footnotes
References
Allen, M. F. (2007). Mycorrhizal fungi: highways for water and nutrients in arid soils. Vadose Zone J. 6, 291–297. doi: 10.2136/vzj2006.0068
Altschul, S. F., Gish, W., Miller, W., Myers, E. W., and Lipman, D. J. (1990). Basic local alignment search tool. J. Mol. Biol. 215, 403–410. doi: 10.1016/S0022-2836(05)80360-2
Álvarez-Lopeztello, J., del Castillo, R. F., Robles, C., and Hernández-Cuevas, L. V. (2019). Spore diversity of arbuscular mycorrhizal fungi in human-modified neotropical ecosystems. Ecol. Res. 34, 394–405. doi: 10.1111/1440-1703.12004
Asmelash, F., Bekele, T., and Birhane, E. (2016). The potential role of arbuscular mycorrhizal fungi in the restoration of degraded lands. Front. Microbiol. 7:1095. doi: 10.3389/fmicb.2016.01095
Beck, S., Powell, J. R., Drigo, B., Cairney, J. W., and Anderson, I. C. (2015). The role of stochasticity differs in the assembly of soil-and root-associated fungal communities. Soil Biol. Biochem. 80, 18–25. doi: 10.1016/j.soilbio.2014.09.010
Bengtsson, J., Bullock, J. M., Egoh, B., Everson, C., Everson, T., O'connor, T., et al. (2019). Grasslands—more important for ecosystem services than you might think. Ecosphere 10:e02582. doi: 10.1002/ecs2.2582
Bissett, A., Brown, M. V., Siciliano, S. D., and Thrall, P. H. (2013). Microbial community responses to anthropogenically induced environmental change: towards a systems approach. Ecol. Lett. 16, 128–139. doi: 10.1111/ele.12109
Bonfim, J. A., Vasconcellos, R. L. F., Gumiere, T., de Lourdes Colombo Mescolotti, D., Oehl, F., and Nogueira 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
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
Chen, Y. L., Zhang, X., Ye, J. S., Han, H. Y., Wan, S. Q., and Chen, B. D. (2014). Six-year fertilization modifies the biodiversity of arbuscular mycorrhizal fungi in a temperate steppe in Inner Mongolia. Soil Biol. Biochem. 69, 371–381. doi: 10.1016/j.soilbio.2013.11.020
Coutinho, E. S., Fernandes, G. W., Berbara, R. L. L., Valério, 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
Cui, N., Shi, L., Guo, J., and Zhang, T. (2021). Arbuscular mycorrhizal fungi alleviate elevated temperature and nitrogen deposition-induced warming potential by reducing soil N2O emissions in a temperate meadow. Ecol. Indic. 131:108193. doi: 10.1016/j.ecolind.2021.108193
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
Davison, J., Moora, M., Semchenko, M., Adenan, S. B., Ahmed, T., Akhmetzhanova, A. A., et al. (2021). Temperature and pH define the realised niche space of arbuscular mycorrhizal fungi. New Phytol. 231, 763–776. doi: 10.1111/nph.17240
Dong, S. (2023). Revitalizing the grassland on the Qinghai–Tibetan plateau. Grassland Res. 2, 241–250. doi: 10.1002/glr2.12055
Dong, S., Shang, Z., Gao, J., and Boone, R. (2022). Enhancing the ecological services of the Qinghai-Tibetan Plateau’s grasslands through sustainable restoration and management in era of global change. Agric. Ecosyst. Environ. 326:107756. doi: 10.1016/j.agee.2021.107756
Duan, H., Xue, X., Wang, T., Kang, W., Liao, J., and Liu, S. (2021). Spatial and temporal differences in alpine meadow, alpine steppe and all vegetation of the Qinghai-Tibetan plateau and their responses to climate change. Remote Sens. 13:669. doi: 10.3390/rs13040669
Dumbrell, A. J., Ashton, P. D., Aziz, N., Feng, G., Nelson, M., Dytham, C., et al. (2011). Distinct seasonal assemblages of arbuscular mycorrhizal fungi revealed by massively parallel pyrosequencing. New Phytol. 190, 794–804. doi: 10.1111/j.1469-8137.2010.03636.x
Faghihinia, M., Zou, Y., Bai, Y., Dudáš, M., Marrs, R., and Staddon, P. L. (2022). Grazing intensity rather than host plant’s palatability shapes the community of arbuscular mycorrhizal fungi in a steppe grassland. Microb. Ecol. 84, 1062–1071. doi: 10.1007/s00248-021-01920-7
Fang, F., Wang, C., Wu, F., Tang, M., and Doughty, R. (2020). Arbuscular mycorrhizal fungi mitigate nitrogen leaching under poplar seedlings. Forests 11:325. doi: 10.3390/f11030325
Frey, S. D. (2019). Mycorrhizal fungi as mediators of soil organic matter dynamics. Annu. Rev. Ecol. Evol. Syst. 50, 237–259. doi: 10.1146/annurev-ecolsys-110617-062331
Gai, J. P., Christie, P., Cai, X. B., Fan, J. Q., Zhang, J. L., Feng, G., et al. (2009). Occurrence and distribution of arbuscular mycorrhizal fungal species in three types of grassland community of the Tibetan plateau. Ecol. Res. 24, 1345–1350. doi: 10.1007/s11284-009-0618-1
Gai, J. P., Tian, H., Yang, F. Y., Christie, P., Li, X. L., and Klironomos, J. N. (2012). Arbuscular mycorrhizal fungal diversity along a Tibetan elevation gradient. Pedobiologia 55, 145–151. doi: 10.1016/j.pedobi.2011.12.004
Gao, Y. Y., Zhang, G. X., Jiang, S. Y., and Liu, L. X. (2024). Wekemo Bioincloud: a user-friendly platform for meta-omics data analyses. iMeta 3:e175. doi: 10.1002/imt2.175
Geml, J. (2017). “Altitudinal gradients in mycorrhizal symbioses: the current state of knowledge on how richness and community structure change with elevation” in Biogeography of mycorrhizal symbiosis, 230.
Gryndler, M., Hršelová, H., Cajthaml, T., Havránková, M., Řezáčová, V., Gryndlerová, H., et al. (2009). Influence of soil organic matter decomposition on arbuscular mycorrhizal fungi in terms of asymbiotic hyphal growth and root colonization. Mycorrhiza 19, 255–266. doi: 10.1007/s00572-008-0217-y
Guo, X., 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
Helgason, T., Feng, H., Sherlock, D. J., Young, J. P. W., and Fitter, A. H. (2014). Arbuscular mycorrhizal communities associated with maples (Acer spp.) in a common garden are influenced by season and host plant. Botany 92, 321–326. doi: 10.1139/cjb-2013-0263
Jamiołkowska, A., Księżniak, A., Gałązka, A., Hetman, B., Kopacki, M., and Skwaryło-Bednarz, B. (2018). Impact of abiotic factors on development of the community of arbuscular mycorrhizal fungi in the soil: a review. Int. Agrophys. 32, 133–140. doi: 10.1515/intag-2016-0090
Jerbi, M., Labidi, S., Bahri, B. A., Laruelle, F., Tisserant, B., Jeddi, F. B., et al. (2021). Soil properties and climate affect arbuscular mycorrhizal fungi and soil microbial communities in Mediterranean rainfed cereal cropping systems. Pedobiologia 87-88:150748. doi: 10.1016/j.pedobi.2021.150748
Jin, L., Zhang, G., Wang, X., Dou, C., Chen, M., Lin, S., et al. (2011). Arbuscular mycorrhiza regulate inter-specific competition between a poisonous plant, Ligularia virgaurea, and a co-existing grazing grass, Elymus nutans, in Tibetan plateau alpine meadow ecosystem. Symbiosis 55, 29–38. doi: 10.1007/s13199-011-0141-3
Kang, F., Yang, B., Wujisiguleng,, Yang, X., Wang, L., Guo, J., et al. (2019). Arbuscular mycorrhizal fungi alleviate the negative effect of nitrogen deposition on ecosystem functions in meadow grassland. Land Degrad. Dev. 31, 748–759. doi: 10.1002/ldr.3491
Kivlin, S. N., Hawkes, C. V., and Treseder, K. K. (2011). Global diversity and distribution of arbuscular mycorrhizal fungi. Soil Biol. Biochem. 43, 2294–2303. doi: 10.1016/j.soilbio.2011.07.012
Koziol, L., and Bever, J. D. (2017). The missing link in grassland restoration: arbuscular mycorrhizal fungi inoculation increases plant diversity and accelerates succession. J. Appl. Ecol. 54, 1301–1309. doi: 10.1111/1365-2664.12843
Lara-Pérez, L. A., Oros-Ortega, I., Córdova-Lara, I., Estrada-Medina, H., O’Connor-Sánchez, A., Góngora-Castillo, E., et al. (2020). Seasonal shifts of arbuscular mycorrhizal fungi in Cocos nucifera roots in Yucatan, Mexico. Mycorrhiza 30, 269–283. doi: 10.1007/s00572-020-00944-0
Latif, A., Ilyas, S., Zhang, Y., Xin, Y., Zhou, L., and Zhou, Q. (2019). Review on global change status and its impacts on the Tibetan plateau environment. J. Plant Ecol. 12, 917–930. doi: 10.1093/jpe/rtz038
Lee, J., Lee, S., and Young, J. P. (2008). Improved PCR primers for the detection and identification of arbuscular mycorrhizal fungi. FEMS Microbiol. Ecol. 65, 339–349. doi: 10.1111/j.1574-6941.2008.00531.x
Lenoir, I., Lounes-Hadj Sahraoui, A., and Fontaine, J. (2016). Arbuscular mycorrhizal fungal-assisted phytoremediation of soil contaminated with persistent organic pollutants: a review. Eur. J. Soil Sci. 67, 624–640. doi: 10.1111/ejss.12375
Li, X., Qi, Z., Yu, X., Xu, M., Liu, Z., Du, G., et al. (2021). Soil pH drives the phylogenetic clustering of the arbuscular mycorrhizal fungal community across subtropical and tropical pepper fields of China. Appl. Soil Ecol. 165:103978. doi: 10.1016/j.apsoil.2021.103978
Li, X., Xu, M., Christie, P., Li, X., and Zhang, J. (2018). Large elevation and small host plant differences in the arbuscular mycorrhizal communities of montane and alpine grasslands on the Tibetan plateau. Mycorrhiza 28, 605–619. doi: 10.1007/s00572-018-0850-z
Li, L. F., Yang, A. N., and Zhao, Z. W. (2005). Seasonality of arbuscular mycorrhizal symbiosis and dark septate endophytes in a grassland site in Southwest China. FEMS Microbiol. Ecol. 54, 367–373. doi: 10.1016/j.femsec.2005.04.011
Liu, M., Yue, Y., Wang, Z., Li, L., Duan, G., Bai, S., 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
López-García, Á., Varela-Cervero, S., Vasar, M., Öpik, M., Barea, J. M., and Azcón-Aguilar, C. (2017). Plant traits determine the phylogenetic structure of arbuscular mycorrhizal fungal communities. Mol. Ecol. 26, 6948–6959. doi: 10.1111/mec.14403
Lu, R. (2000). Methods for analysis of soil agrochemistry. Beijing: China Agricultural Science and Technology Press.
Ma, X., Geng, Q., Zhang, H., Bian, C., Chen, H. Y., Jiang, D., et al. (2021). Global negative effects of nutrient enrichment on arbuscular mycorrhizal fungi, plant diversity and ecosystem multifunctionality. New Phytol. 229, 2957–2969. doi: 10.1111/nph.17077
Mao, L., Pan, J., Jiang, S., Shi, G., Qin, M., Zhao, Z., et al. (2019). Arbuscular mycorrhizal fungal community recovers faster than plant community in historically disturbed Tibetan grasslands. Soil Biol. Biochem. 134, 131–141. doi: 10.1016/j.soilbio.2019.03.026
Mao, H., Whalen, J. K., Zhang, Z., Sheng, X., Hu, G., Chen, B., et al. (2024). Arbuscular mycorrhizal fungi associated with alpine meadow multifunctionality in a warmer climate with variable precipitation. Soil Biol. Biochem. 198:109555. doi: 10.1016/j.soilbio.2024.109555
Miller, R. M., and Jastrow, J. D. (2000). “Mycorrhizal fungi influence soil structure” in Arbuscular mycorrhizas: physiology and function, 3–18.
Miller, R. M., Wilson, G. W., and Johnson, N. C. (2012). “Arbuscular mycorrhizae and grassland ecosystems” in Biocomplexity of plant–fungal interactions, 59–84.
Mohammad, M. J., Hamad, S. R., and Malkawi, H. I. (2003). Population of arbuscular mycorrhizal fungi in semi-arid environment of Jordan as influenced by biotic and abiotic factors. J. Arid Environ. 53, 409–417. doi: 10.1006/jare.2002.1046
Neuenkamp, L., Zobel, M., Lind, E., Gerz, M., and Moora, M. (2019). Arbuscular mycorrhizal fungal community composition determines the competitive response of two grassland forbs. PLoS One 14:e0219527. doi: 10.1371/journal.pone.0219527
Niu, Y. J., Zhu, H. M., Yang, S. W., Ma, S. J., Zhou, J. W., Chu, B., et al. (2019). Overgrazing leads to soil cracking that later triggers the severe degradation of alpine meadows on the Tibetan plateau. Land Degrad. Dev. 30, 1243–1257. doi: 10.1002/ldr.3312
Opik, M., Vanatoa, A., Vanatoa, E., Moora, M., Davison, J., Kalwij, J. M., et al. (2010). The online database MaarjAM reveals global and ecosystemic distribution patterns in arbuscular mycorrhizal fungi (Glomeromycota). New Phytol. 188, 223–241. doi: 10.1111/j.1469-8137.2010.03334.x
Pereira, C. M. R., da Silva, D. K. A., Goto, B. T., Rosendahl, S., and Maia, L. C. (2018). Management practices may lead to loss of arbuscular mycorrhizal fungal diversity in protected areas of the Brazilian Atlantic Forest. Fungal Ecol. 34, 50–58. doi: 10.1016/j.funeco.2018.05.001
Shi, G., Yang, Y., Liu, Y., Uwamungu, J. Y., Liu, Y., Wang, Y., et al. (2022). Effect of Elymus nutans on the assemblage of arbuscular mycorrhizal fungal communities enhanced by soil available nitrogen in the restoration succession of revegetated grassland on the Qinghai-Tibetan plateau. Land Degrad. Dev. 33, 931–944. doi: 10.1002/ldr.4201
Su, Y. Y., Sun, X., and Guo, L. D. (2011). Seasonality and host preference of arbuscular mycorrhizal fungi of five plant species in the Inner Mongolia steppe, China. Braz. J. Microbiol. 42, 57–65. doi: 10.1590/S1517-83822011000100008
Suzuki, K., Takahashi, K., and Harada, N. (2020). Evaluation of primer pairs for studying arbuscular mycorrhizal fungal community compositions using a MiSeq platform. Biol. Fertil. Soils 56, 853–858. doi: 10.1007/s00374-020-01431-6
Torrez, V., Ceulemans, T., Mergeay, J., de Meester, L., and Honnay, O. (2016). Effects of adding an arbuscular mycorrhizal fungi inoculum and of distance to donor sites on plant species recolonization following topsoil removal. Appl. Veg. Sci. 19, 7–19. doi: 10.1111/avsc.12193
Treseder, K. K., and Allen, M. F. (2002). Direct nitrogen and phosphorus limitation of arbuscular mycorrhizal fungi: a model and field test. New Phytol. 155, 507–515. doi: 10.1046/j.1469-8137.2002.00470.x
Van Aarle, I. M., Olsson, P. A., and Söderström, B. (2002). Arbuscular mycorrhizal fungi respond to the substrate pH of their extraradical mycelium by altered growth and root colonization. New Phytol. 155, 173–182. doi: 10.1046/j.1469-8137.2002.00439.x
Vasconcellos, R. L., Bonfim, J. A., Baretta, D., and Cardoso, E. J. (2016). Arbuscular mycorrhizal fungi and glomalin-related soil protein as potential indicators of soil quality in a recuperation gradient of the Atlantic Forest in Brazil. Land Degrad. Dev. 27, 325–334. doi: 10.1002/ldr.2228
Wang, J., Hu, A., Meng, F., Zhao, W., Yang, Y., Soininen, J., et al. (2022). Embracing mountain microbiome and ecosystem functions under global change. New Phytol. 234, 1987–2002. doi: 10.1111/nph.18051
Wang, J., Yang, C. X., Sun, X. M., Zhang, H. O., Guo, Z., Cao, T. T., et al. (2021). Growth responses of Sorghum bicolor (L.) Moench to arbuscular mycorrhizal fungi under simulated nitrogen deposition. Bangladesh J. Botany 50, 933–938. doi: 10.3329/bjb.v50i0.56447
Wu, X., Yang, J., Ruan, H., Wang, S., Yang, Y., Naeem, I., et al. (2021). The diversity and co-occurrence network of soil bacterial and fungal communities and their implications for a new indicator of grassland degradation. Ecol. Indic. 129:107989. doi: 10.1016/j.ecolind.2021.107989
Xiao, D., Chen, Y., He, X., Xu, Z., Bai, S. H., Zhang, W., et al. (2021). Temperature and precipitation significantly influence the interactions between arbuscular mycorrhizal fungi and diazotrophs in karst ecosystems. For. Ecol. Manag. 497:119464. doi: 10.1016/j.foreco.2021.119464
Zhang, C. J., Kang, Y. K., Yao, B. H., An, K., Pu, Q. S., Wang, Z. C., et al. (2022). Increased availability of preferred food and decreased foraging costs from degraded grasslands lead to rodent pests in the Qinghai-Tibet plateau. Front. Ecol. Evol. 10:971429. doi: 10.3389/fevo.2022.971429
Zhang, F., Li, Y., Ji, B., and Dong, S. (2024). Spatial distribution and drivers of arbuscular mycorrhizal fungi on the Tibetan plateau. Front. Plant Sci. 15:1427850. doi: 10.3389/fpls.2024.1427850
Zhang, T., Sun, Y., Shi, Z., and Feng, G. (2012). Arbuscular mycorrhizal fungi can accelerate the restoration of degraded spring grassland in Central Asia. Rangel. Ecol. Manag. 65, 426–432. doi: 10.2111/REM-D-11-00016.1
Zhang, L., Wang, X. J., Wang, Q., and Jin, L. (2017). Advances in the study of arbuscular mycorrhizal fungi in high altitude and cold habitats on Tibetan plateau. J. Fungal Res. 15:12. doi: 10.13341/j.jfr.2014.1710
Zhang, W., Xue, X., Peng, F., You, Q., and Hao, A. (2019). Meta-analysis of the effects of grassland degradation on plant and soil properties in the alpine meadows of the Qinghai-Tibetan plateau. Glob. Ecol. Conserv. 20:e00774. doi: 10.1016/j.gecco.2019.e00774
Zhao, Y., Liu, Z., and Wu, J. (2020). Grassland ecosystem services: a systematic review of research advances and future directions. Landsc. Ecol. 35, 793–814. doi: 10.1007/s10980-020-00980-3
Keywords: AM fungi, seasonal changes, Illumina HiSeq sequencing, soil factors, community composition
Citation: Dong W, Ding T and Duan T (2025) Diversity of arbuscular mycorrhizal fungi and its response to seasonal variation in alpine grassland of the eastern Tibetan Plateau. Front. Microbiol. 16:1511979. doi: 10.3389/fmicb.2025.1511979
Edited by:
Ming Liu, Chinese Academy of Sciences (CAS), ChinaReviewed by:
Jiangbing Xu, Nanjing University of Information Science and Technology, ChinaRuibo Sun, Anhui Agricultural University, China
Copyright © 2025 Dong, Ding and Duan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tingyu Duan, ZHVhbnR5QGx6dS5lZHUuY24=