Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Microbiol., 25 July 2025

Sec. Microbiological Chemistry and Geomicrobiology

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1627840

Ecosystem carbon storage and carbon metabolizing microorganisms in three types of grasslands on the Qinghai-Tibet Plateau

Qian LiuQian Liu1Wenquan YangWenquan Yang2Jiancun Kou
Jiancun Kou1*Qinyao LiQinyao Li1Yangcan ZhangYangcan Zhang1Xilai LiXilai Li3Jing ZhangJing Zhang3Zhiting HaoZhiting Hao1Lu ChiLu Chi1Yuze NingYuze Ning1
  • 1College of Grassland Agriculture, Northwest A&F University, Yangling, China
  • 2College of Life Sciences, Northwest A&F University, Yangling, China
  • 3College of Agriculture and Animal Husbandry, Qinghai University, Xining, China

Background: The response of soil microorganisms to environmental changes can affect the storage and stability of carbon pools in ecosystems. However, the intrinsic link between the structure of soil carbon-metabolizing microbial communities and their roles in different types of alpine grasslands remains unclear.

Methods: This study explores how carbon storage varies among alpine meadow (AM), alpine wetland (AW), and alpine desert (AD) on the Qinghai-Tibet Plateau and assesses the influence of a wide range of soil microbial and vegetation factors, so as to identify microbial predictors of ecosystem carbon storage. The study revealed four types of carbon metabolizing microbial communities responded to changes in vegetation types and their impact on the storage and stability of carbon pools.

Results: The carbon storage of three grassland types followed the relativity of AW > AM > AD. Soil water content (SWC) was identified as the major factor affecting the carbon storage of grassland ecosystems by increasing vegetation belowground biomass and soil total carbon content, directly or indirectly influencing the diversity of four types of soil microorganisms through its effects on soil physicochemical properties. The community structure of these four types of carbon metabolizing microorganisms in AW significantly differed from that of AM and AD. The diversity of carbon-fixing microorganisms significantly reduced ecosystem carbon storage to a great extent. The relative abundance of carbon-fixing microorganisms Thiobacillus, Mesorhizobium, Azospirillum, and Methylibium significantly increased grassland carbon storage, while the relative abundance of chitinase-producing microorganisms Cellulomonas and Stenotrophomonas significantly decreased it.

1 Introduction

The Qinghai-Tibet Plateau, with an average altitude of over 4,000 m, is home to the world’s highest ecosystem, possessing an enormous potential for carbon sequestration and playing a crucial role in the global carbon balance (Chen H. et al., 2022; Canedoli et al., 2020). The alpine grassland on the Plateau covers more than 60% of its total area (Gilmanov et al., 2007). Although the soil is generally poor, the region is rich in vegetation and is one of the global biodiversity hotspots (Liu et al., 2018; Berdugo et al., 2019). Therefore, it is highly sensitive to changes in environmental conditions such as climate, vegetation, and soil physicochemical properties (Bai et al., 2012; Cleland et al., 2013).

Major grasslands atop the Qinghai-Tibet Plateau fall into three types of alpine meadow, alpine wetland, and alpine desert. Their distribution pattern is jointly shaped by water and heat. In the alpine environment of the Plateau, differential soil water content (SWC) is generally regarded as the key environmental factor driving the formation of different types of alpine grasslands (Du and Gao, 2020). Environmental changes (such as alterations in climate, vegetation, and soil properties) profoundly influence the carbon cycle processes in plateau grasslands (Ganjurjav et al., 2016; Yang et al., 2021). Among them, water availability (rather than mere precipitation) acts on carbon storage through multiple pathways (Du and Gao, 2020; Sun et al., 2022). It not only directly modulates plant growth (Deng et al., 2017; Sun et al., 2020) and vegetation carbon accumulation, but also, more critically, drives the soil carbon cycle by altering SWC, and thus affecting soil microbial communities (Li Y. et al., 2023; Li Z. et al., 2023). Soil microorganisms dominate the carbon transformation processes through catabolic and anabolic metabolism (Philippot et al., 2013; Liang et al., 2017). For instance, soil carbon-fixing microorganisms mediate the Calvin cycle (Ge et al., 2013), while methanogens and methanotrophs control methane production and oxidation, respectively (Juottonen et al., 2006), and chitinase-producing microorganisms are involved in the mineralization and decomposition of organic matter (such as chitin) (Brankatschk et al., 2011). Therefore, SWC regulates the core internal mechanisms of soil carbon pools (carbon fixation, transformation, and loss) in different grasslands by influencing the community structure, diversity, and activity of such key carbon-fixing microorganisms as methanogens, methanotrophs, and chitinase-producing microorganisms.

Previous studies have confirmed that soil microbial communities are susceptible to the nature of grassland and the associated environmental variables (such as vegetation characteristics and soil physicochemical properties) (Li Y. et al., 2023; Li Z. et al., 2023; Prober et al., 2015; Jones et al., 2018; Wang et al., 2024). Current research on microorganisms in alpine grasslands has focused primarily on the overall bacterial and fungal communities (Li Y. et al., 2023). Extremely scarce research has been carried out to specifically target functional microbial communities (carbon-fixing bacteria, methanogens, methanotrophs, and chitinase-producing bacteria) that drive core carbon cycling processes in different grasslands. As a result, it remains unknown how variations in SWC influence these key carbon-functional microbial communities and ultimately lead to differences in soil carbon pools among alpine meadows, wetlands, and deserts. Therefore, it is necessary to explore the community variation characteristics of carbon-metabolizing microorganisms in different types of grassland.

This study analyzes the differences in ecosystem carbon storage among three types of alpine grasslands, alpine meadow (AM), alpine wetland (AW), and alpine desert (AD) in the Qinghai-Tibet Plateau using the illumina sequencing technology. It focuses on the impact mechanisms of SWC on ecosystem carbon storage and soil carbon-metabolizing microbial communities in different grasslands. The objectives of this study are to: (1) investigate the differences in grassland carbon storage affected by SWC; (2) assess the influence of SWC on soil properties and vegetation characteristics among three types of alpine grasslands; (3) explore the influence of SWC on the diversity of carbon metabolizing microorganisms; and (4) ascertain the impact of soil carbon metabolizing microorganisms on grassland carbon storage. The results not only provide basic data for the study of carbon cycle and transformation in different types of alpine grasslands, but also provide reference and guidance for the study of carbon source and carbon storage control mechanism in different water conditions of alpine grassland.

2 Materials and methods

2.1 Study site

The experiments were conducted in Menyuan County, Haibei Prefecture, Qinghai Province, China (37°19′-37°47’N, 101°14′-102°1′E) (Figure 1a). The area has an average altitude of 3,681 m and a typical plateau alpine climate. The annual average temperature is 0.8°C, with a daily temperature range of 11.6–17.5°C, and the annual average precipitation is 520 mm. Alpine meadow is the main grassland type, and the soil is plateau meadow soil (Nachtergaele et al., 2012). The main dominant plants are Kobresia myosuroides, Carex spp., and Elymus dahuricus. Additionally, there are also pockets of alpine wetland and alpine desert. The former is mainly distributed in low-lying flat areas between mountains, with the dominant plants being Carex spp., Kobresia myosuroides, and Polygonum aviculare. The latter is primarily located near seasonal snow-covered zones on mountain tops, dominated by Saussurea pulchra and Leontopodium leontopodioides. The appearance of the three types of alpine grasslands is shown in Figure 1b. To account for regional heterogeneity, each type was replicated by selecting three spatially distributed sites at a distance of >5 km from each other. Each site-type combination was treated as an independent replicate to validate vegetation-specific differences in carbon storage and microbial communities.

Figure 1
Map and terrain visuals highlighting a region with varying altitudes and vegetation types. The map shows alpine meadows, wetlands, and deserts with a color gradient indicating altitude. Insets (b) feature images of terrain types: AM (Alpine Meadow) with lush green grass, AW (Alpine Wetland) with sparse vegetation, and AD (Alpine Desert) with rocky, barren ground. A compass and scale are included for orientation.

Figure 1. Sampling locations (a) and grassland types (b). AM, alpine meadow; AW, alpine wetland; AD, alpine desert.

2.2 Soil and plant sample collection

Soil and plant samples were collected from July 28 to July 31, 2023 on sunny days for the three typical types of alpine grasslands (AM, AW and AD) at three experimental sites each (Figure 1). Within each site, three random plots (10 m × 10 m) serving as three replicates were established, with a distance of >10 m between adjoining plots. Three quadrats (1 m × 1 m) were set up in each plot (Supplementary Figure S1). After vegetation survey in each quadrat, above-ground parts were collected and all above-ground parts of plants in the same quadrat were clipped to determine vegetation biomass by types of plants (grasses, sedges and forbs), the total biomass of the quadrat and subsequent carbon content analysis. The three-replicate average was regarded as the final results for each analyzed parameter. After clearing the surface cover, soils were sampled to a depth of 0–20 cm using a soil auger with a diameter of 5 cm. Five soil cores were taken from each plot, resulting in a total of 15 soil cores from the same plot (10 m × 10 m) (total sampled area = 15 × π × (0.05/2)2 ≈ 0.029 m2). They were then mixed to form one composite soil sample for the subsequent collection of vegetation roots. A total of 9 soil samples were collected from 3 kinds of grasslands for calculating the root biomass of vegetation (3 types of grassland × 3 sites). According to the above method, five soil cores were additionally collected from each plot and mixed for the analysis of microorganisms and soil physical and chemical properties. Immediately, a portion of it was transferred into 1.5 mL cryogenic vials, rapidly placed into a low-temperature storage container, and dispatched immediately to Shanghai Personalbio Biotechnology Co., Ltd. (Shanghai, China) for microbial analysis. In total, 27 soil samples were collected for the three types of grassland (3 types of grassland × 3 sites × 3 replicates).

2.3 Soil physicochemical properties and vegetation characteristics analysis

In the laboratory, the collected soil samples were treated by differential processing. The composite sample mixed with 15 soil cores was continuously rinsed with flowing water to separate plant roots from soil particles, and then dried in an oven at 65°C until the weight was stable. The dried roots and aboveground biomass were homogenized through 0.25 mm mesh grinding for plant carbon quantification using potassium dichromate oxidation (Li et al., 2017). The soil collected for the analysis of soil physical and chemical properties was divided into two parts. The first portion was air-dried and sieved through a 0.149 mm mesh to determine soil total carbon content and soil physicochemical properties, including soil total nitrogen (TN), ammonium nitrogen (NH4+-N) and nitrate nitrogen (NO3-N) using the AA3 continuous flow analyzer (Auto Analyzer 3-AA3). Soil organic matter (SOM) was determined by the potassium dichromate-sulfuric acid oxidation method (Fu et al., 2004). Total phosphorus (TP) and available phosphorus (AP) of the soil were analyzed using ultraviolet spectrophotometer (UV-2450). Soil total carbon was measured with a soil total carbon analyzer (Primacs ATC 100 analyzer, Skalar, Netherlands). Soil total potassium (TK) was determined using an atomic absorption spectrometer (Shimadzu A6300, Japan). The second fraction was preserved for characterizing soil physical properties, including soil bulk density (BD), SWC, and pH. The first two were measured gravimetrically after 2 mm sieving (Wang et al., 2016), BD was determined by drying the sampled soil at 105°C to a constant weight. SWC was determined by oven-drying fresh samples in pre-weighed aluminum tins at 105°C to constant weight, with values calculated from mass loss. Soil pH was determined using a Palintest SKW500 sensor at a soil-water ratio of 1:2.5.

2.4 Carbon storage calculation

Soil carbon storage was calculated using the following formula (Kumar et al., 2012) (Equation 1):

Soil TCstock = i = 1 n ( TC × BD × di )     (1)

TC: total carbon; BD: soil bulk density (g cm−3); di: soil layer thickness (cm) (its value was taken as the maximum of the 0–20 cm depth).

Plant carbon stock, including above-ground and below-ground parts of plants, was calculated according to the following formula (Li et al., 2017) (Equations 2, 3):

X stock = X × C     (2)
Plant C stock = shoot C stock + root C stock     (3)

Where X represents the measured plant biomass (g m−2) (above-ground or root biomass), while C stands for the carbon content of above-ground parts and roots (g kg−1) of plants.

The formula for calculating ecosystem carbon stock was as follows (Wang et al., 2023) (Equation 4):

Ecosystem C stock = Plant C stock + Soil TC stock     (4)

2.5 Sequencing and raw data analysis of pmoA, cbbL, mcrA, and chiA genes

Soil genomic DNA was extracted with the OMEGA Soil DNA Kit (Omega Bio-Tek, United States). Microbial community profiling involved PCR amplification of target regions was carried out using primers listed in Table 1, with the thermal cycling parameters set as: initial denaturation at 95°C (3 min), 27 cycles of 95°C (30 s), 55°C (30 s), and 72°C (45 s), concluding with final extension at 72°C (10 min). Purified amplicons were sequenced on the Illumina MiSeq PE300 platform (Majorbio Bio-Pharm, China).

Table 1
www.frontiersin.org

Table 1. Specific primers, fragment lengths and annealing temperature of functional genes.

Bioinformatic processing included primer removal via cutadapt and sequence deduplication using the derep_fulllength module. Chimeric sequences were eliminated through uchime_denovo analysis. High-quality sequences were clustered into operational taxonomic units (OTUs) at the 97% similarity using the cluster_size module. Sequence normalization was achieved using rarefaction analysis (Kemp and Aller, 2004) to standardize sequencing depth across samples, enabling comparative assessment of OTU distributions and relative abundances. Representative sequences and OTU tables were generated for subsequent ecological analyzes.

2.6 Statistical analysis

SPSS 26.0 software (IBM Corporation, Armonk, NY, United States) was utilized to conduct a one-way analysis of variance (ANOVA), followed by a Tukey test to evaluate the characteristics of different types of alpine grassland plants, soil physicochemical properties, α-diversity of soil carbon metabolizing microorganisms, and the relative abundance of soil carbon metabolizing microbial phyla and genera. Origin 2018 was used to create bar charts of carbon components across different grassland types. Based on the OTU analysis results, the vegan package in R (v.3.5.2) was employed to calculate α diversity indices (Chao1 and Shannon), perform principal coordinate analysis (PCoA), and construct random forest models for species. The vegan and ggplot2 modules in R (v.3.5.2) were utilized to create redundancy analysis (RDA) plots of environmental factors and carbon-metabolizing microbial communities. Differences in the composition of carbon-metabolizing microorganisms among different types of alpine grasslands were tested using a matrix based on Bray–Curtis dissimilarity, and visualized through non-metric multidimensional scaling (NMDS) and PCoA based on Bray–Curtis dissimilarity. Non-parametric tests consistent with multivariate analysis of variance (ADONIS, ANOSIM, MRPP) were used to assess the significance of these differences. Venn analysis was conducted to determine the number of shared and unique OTUs among the three types of alpine grasslands. Potential microbial biomarkers were obtained using the Linear Discriminant Analysis Effect Size (LEfSe) method. Correlation analysis and Mantel tests were conducted to investigate the relationships between carbon storage, vegetation characteristics, soil physicochemical properties, and β-diversity and species richness of carbon-metabolizing microbial communities. Random forest modeling was used to predict the relative importance of α and β diversity indices of soil carbon-metabolizing microorganisms, and the relative abundance of key genera of carbon metabolizing microorganisms and their impact on carbon storage. Predictors with negative IncMSE values were excluded prior to model fitting. The final random forest model, optimized at 1000 trees via cross-validation, quantified variable importance through permutation-based significance testing.

Partial Least Squares Structural Equation Modeling (PLS-SEM) was carried out using SmartPLS 4.0 to quantify the direct and indirect effects of SWC on carbon storage and carbon metabolizing microorganisms across the grassland types. PLS-SEM was selected for its robustness with a limited sample size and abnormally distributed data (Jöreskog and Goldberger, 1975; Hair et al., 2019). Within this framework, we analyzed the influence of SWC on soil physical and chemical properties (pH, BD, TC, TN, TP, TK, SOM, NO3-N, NH4+-N and AP), vegetation characteristics (AGB and BGB), and the diversity of carbon-metabolizing microorganisms (Shannon and Simpson indices), as well as the pathways and extent of their impact on the carbon storage of grassland ecosystems. Model validity was assessed using three metrics: the coefficient of determination (R2), predictive relevance (Q2), and goodness of fit (GoF). An R2 value of approximately 0.67 is considered substantial, a value around 0.33 indicates a moderate explanatory power, and a value around 0.19 suggests a weak explanatory power. When Q2 is greater than 0, the SEM model has a higher predictive relevance. A GoF value less than 0.1 indicates a weak explanatory power, around 0.25 a moderate explanatory power, and around 0.36 a strong explanatory power (Lu et al., 2024). The predictive performance of PLS-SEM was compared with that of the data-driven random forest method to illustrate the applicability of PLS-SEM for predictive purposes.

3 Results

3.1 Characteristics of plants and soil in different types of alpine grassland

There were significant differences in plant characteristics among the three types of alpine grassland (p < 0.05). The aboveground biomass of grasses, sedges and forbs in AM was significantly higher than in the other two types of alpine grassland, and the Shannon-Wiener index of biomass and vegetation on miscellaneous grassland in AW was the lowest (p < 0.05). In contrast, vegetation coverage and belowground biomass were significantly higher in AW than in AM and AD (p < 0.05). Soil pH ranged from 6.52 to 7.13, and the soil pH of AD was significantly higher than that of AM and AW (p < 0.05). The contents of SWC, TC, TN, SOM and NH4+-N in AW were significantly higher than those in AM and AD, while the contents of BD and TK in AW were significantly lower than those in AM and AD (p < 0.05). Except for soil pH, there was no significant difference in the content of other soil physical and chemical properties between AM and AD. There was no significant difference in the contents of TP and NO3-N among the three types of alpine grassland (Table 2).

Table 2
www.frontiersin.org

Table 2. Plant characteristics and soil physicochemical properties of three types of alpine grassland.

3.2 Carbon component differences across alpine grassland types

Among the three types of alpine grassland, the highest carbon storage was found in the soil, accounting for over 90% of the ecosystem’s total carbon, followed by plant roots, with the lowest carbon storage in the belowground parts of plants (Figure 2a). Significant differences in carbon components were observed among AM, AW, and AD (p < 0.05). Each carbon component was lower in AD than in AM and AW (p < 0.05). The carbon storage in plant roots, vegetation, soil, and the entire ecosystem was the highest in AW (p < 0.05). However, the belowground vegetation carbon storage was significantly higher in AM than in AW and AD (p < 0.05) (Figure 2b).

Figure 2
Bar graphs showing carbon storage (g per cubic centimeter) in various components: (a) Shoot, Root, and Plant carbon across groups AM, AW, and AD. (b) Soil and Ecosystem carbon for the same groups. Group AW consistently shows higher carbon storage than AM and AD. Error bars indicate variability.

Figure 2. Differences in carbon storage in plant (a) and in soil and ecosystems (b) among three types of alpine grassland. AM, alpine meadow; AW, alpine wetland; AD, alpine desert. Different letters indicate significant differences at the 0.05 level.

3.3 Microbial community diversity and community composition of soil carbon metabolism

The α diversity indices of soil microbial communities encompass richness indices and diversity indices, characterized by the Chao1, Shannon, and Simpson indices, respectively. Results from the Spearman analysis, as shown in Figure 3, indicate differences in the diversity and richness of soil carbon-metabolizing microbial communities among the three types of alpine grasslands. The Chao1, Shannon, and Simpson indices of the four types of carbon-metabolizing microorganisms were the lowest in AW, while the Chao1 index and Shannon index of methanotrophs were significantly higher in AM than in AW (p < 0.05) (Figure 3a). The trends in the Chao1 index, Shannon index, and Simpson diversity of soil carbon-fixing microorganisms were consistent across the three types of alpine grasslands, with the lowest values found in AW and the highest in AD (Figure 3b). The Chao1 value of the methanogens community showed little variation, whereas the Shannon and Simpson indices were significantly higher in AM than in AW (p < 0.05) (Figure 3c). The Chao1 value of the chitinase-producing microorganism community was significantly higher in AD than in AM, while no significant difference in the Shannon and Simpson indices existed between AM and AD (Figure 3d).

Figure 3
Grid of boxplots with four rows labeled (a), (b), (c), and (d). Each row has three boxplots showing Chao1, Shannon, and Simpson indices for groups AM, AW, and AD. The boxplots show data distribution with median lines and letter annotations indicating statistical differences.

Figure 3. Box plots of α diversity index of soil carbon metabolizing microorganisms in three types of alpine grassland. They show Chao1, Shannon and Simpson indices of methanotrophs (a), carbon-fixing microorganism (b), methanogens (c), and chitinase-producing microorganism (d) communities. AM, alpine meadow; AW, alpine wetland; AD, alpine desert. Different letters indicate significant differences at the 0.05 level.

The Venn diagram showed that the number of unique OTUs for methanotrophs and methanogenic microorganisms was the highest in AW (Figures 4a,c), while the number of unique OTUs for carbon-fixing microorganisms and chitinase-producing microorganisms was the highest in AD (Figures 4b,d). The three types of alpine grasslands shared 582 (56.0%) OTUs for methanotrophs, 1,066 (37.4%) OTUs for carbon-fixing microorganisms, 263 (65.7%) OTUs for methanogens, and 369 (52.0%) OTUs for chitinase-producing microorganisms (Figures 4ad).

NMDS and PCoA results revealed discernible differences in the composition of the soil carbon-metabolizing microbial community among the alpine grasslands. Specifically, the communities of methanotrophs, carbon-fixing microorganisms, methanogens and chitinase-producing microorganisms in AW were significantly more distinctive from each other than they were in AM and AD (Figures 4eh), indicating significant differences in the structure of the soil carbon metabolizing microbial community between AW and AM/AD. The carbon-fixing microorganism communities were significantly separated among AM, AW, and AD (Figure 4f), suggesting significant differences in the soil carbon-fixing microbial community structure among the alpine grasslands. PCoA revealed that, at the OTU level, there were significant differences (p < 0.05) in the soil carbon metabolizing microbial communities among different types of alpine grasslands, with the first two principal coordinate axes explaining 42.2, 32.9, 35.2, and 31.8% of the total variation for methanotrophs, carbon-fixing microorganisms, methanogens, and chitinase-producing microorganisms, respectively (Figures 4c,f,i,l).

Figure 4
The image contains twelve sub-panels. Panels (a) to (d) are Venn diagrams showing overlapping data sets labeled AM, AW, and AD with various intersections, percentages, and numbers. Panels (e) to (h) depict NMDS plots with analyses of group separations, marked by stress values. Panels (i) to (l) are PCoA plots showing data point distributions in different colored ellipses, labeled by percentage contributions of PC1 and PC2 axes. Each panel uses colors red for AM, green for AW, and blue for AD.

Figure 4. Soil carbon metabolizing microbial β diversity in three types of alpine grassland. The Venn diagram shows the number of unique and common OTUs of methanotrophs (a), carbon-fixing microorganism (d), methanogens (g), and chitinase-producing (j) microorganisms communities; NMDS analysis of methanotrophs (b), carbon-fixing microorganism (e), methanogenic (h) and chitinase-producing microorganisms (k); Principal coordinate analysis (PCoA) of OTU level methanotrophs (c), carbon-fixing microorganism (f), methanogens (i), and chitinase-producing microorganisms (l) communities based on Bray-Curtis distances. AM, alpine meadow; AW, alpine wetland; AD, alpine desert.

In terms of classification and composition, at the phylum level, the dominant phylum of all types of methanotrophs and carbon-fixing microorganisms in AM was Pseudomonadota, and the dominant phylum of methanogens was Euryarchaeota. The dominant phyla of chitinase-producing microorganisms were Pseudomonadota and Actinomycetota (Figure 5). At the genus level, the three types of alpine grasslands differed in their taxonomic composition of methanotrophs. Specifically, AM and AD were dominated by the Methylococcus (33.08–38.56%) species in the methanotrophs community, AW was dominated by the Methylocystis (10.25–92.84%) genus (Figure 6a). The top-10 genera components (in relative abundance) of carbon-fixing microorganism community were identified as Mesorhizobium (4.64–25.08%), Azospirillum (7.30–10.67%), Variovorax (1.41–14.69%), Methylibium (4.85–10.37%), Rubrivivax (6.61–9.21%), Thiobacillus (0.05–17.34%), Bradyrhizobium (3.67–6.14%), Cupriavidus (2.43–5.31%), Dokdonella (3.06–3.62%), and Rhodopseudomonas (1.07–4.61%) (Figure 6b). Among them, the relative abundance of Mesorhizobium, Azospirillum, and Methylibium was the highest in AM and the lowest in AW. In contrast, the relative abundance of Thiobacillus and Rhodopseudomonas was significantly higher in AW than in AM and AD (p < 0.05). The dominant components of methanogens community were Methanobrevibacter (27.52–48.57%) and Methanoregula (10.06–22.46%), and there were no significant differences in the composition of methanogens community among the three types of alpine grasslands (Figure 6c). In the composition of chitinase-producing microorganism community, the genera with an average relative abundance greater than 1% were Streptomyces (19.49–25.56%), Lysobacter (8.57–17.68%), Stenotrophomonas (4.34–17.03%), and Janthinobacterium (5.04–16.42%). Among them, the relative abundance of Streptomyces and Stenotrophomonas in AD was significantly higher than that in AM (p < 0.05) (Figure 6d).

Figure 5
Four stacked bar charts show the relative abundance of different microbial taxa across three sites labeled AM, AW, and AD. Chart (a) includes Pseudomonadota and others. Chart (b) includes Pseudomonadota, Actinomycetota, among others. Chart (c) focuses on Euryarchaeota, Candidatus Thermoplasmatota, and others. Chart (d) features Pseudomonadota, Actinomycetota, and additional groups. Each chart is differentiated by color-coded legends indicating specific taxa.

Figure 5. Soil community composition of (a) methanotrophs, (b) carbon-fixing microorganism, (c) methanogens and (d) chitinase-producing microorganism at the phylum level. AM, alpine meadow; AW, alpine wetland; AD, alpine desert.

Figure 6
Four bar graphs (a, b, c, d) showing the relative abundance percentages of different microbial taxa across three sites: AM, AW, and AD. Each graph highlights the proportion of various genera, with a legend indicating corresponding color-coded taxa. Graph (a) focuses on methanotrophic bacteria, (b) on nitrogen-fixing bacteria, (c) on methanogenic archaea, and (d) on actinobacteria. The graphs illustrate variations in microbial community structures at these sites.

Figure 6. Soil community composition of (a) methanotrophs, (b) carbon-fixing microorganism, (c) methanogens and (d) chitinase-producing microorganism at the genus level. AM, alpine meadow; AW, alpine wetland; AD, alpine desert.

LEfSe analysis revealed that there are 65 taxonomic groups of soil carbon metabolism microorganisms, including 2 kingdoms, 4 phyla, 10 classes, 13 orders, 12 families, and 14 genera. They showed statistically significant differences among the three types of grassland. Specifically, Methylosinus and Methylacidiphilum are the dominant methanotrophs in AW (Supplementary Figure S2a). Mesorhizobium, Variovorax, and Methylibium dominate the carbon-fixing microorganisms in AM, while Thiobacillus, Rhodopseudomonas, and Sulfuricaulis exhibit an absolute dominance in AW (Supplementary Figure S2b). Methanobrevibacter and Methanosphaera show a clear dominance in AW and AD, respectively (Supplementary Figure S2c). Stenotrophomonas, Chitiniphilus, and Amycolatopsis are the dominant chitinase-producing microorganisms in AD, while Janthinobacterium exhibit a clear dominance in AM (Supplementary Figure S2d).

3.4 Relationships between soil carbon-metabolizing microbial communities and environmental factors

Based on the Mantel test results, soil characteristics (SOM, NO3-N, NH4+-N, AP, TN, TC, SWC) and BGB were found to be significantly positively correlated with carbon storage (p < 0.01) (Figure 7a). Carbon storage was significantly correlated with the community structure of methanotrophs and carbon-fixing microorganisms. The community structure of methanotrophs was significantly correlated with NH4+-N (r = 0.424, p = 0.001; r = 0.495, p = 0.001) (Supplementary Table S1). Except for the methanogens community, soil physicochemical properties and vegetation characteristics were significantly correlated with the community structure of methanotrophs, carbon-fixing microorganisms, and chitinase-producing microorganisms. Figure 7b showed that carbon storage was significantly negatively correlated with the diversity indices (Chao1, Shannon, and Simpson) of carbon-fixing microorganisms and the diversity indices (Chao1 and Shannon) of chitinase-producing microorganisms (p < 0.05).

Figure 7
Four panels display species and environmental factor ordinations alongside bar graphs. Panels (a), (c), (e), (g) present factor ordination plots with species vectors. Panels (b), (d), (f), (h) show bar graphs of contribution rates for vegetation, soil, and carbon storage. Vegetation and soil contribute predominantly, with minimal contributions from carbon storage in each respective graph.

Figure 7. RDA results showing the relationship between environmental factors and soil (a) methanotrophs, (c) carbon-fixing microbial, (e) methanogens, and (g) chitinase-producing microbial communities; contribution of influencing factors to the diversity of (b) methanotrophs, (d) carbon-fixing microbial, (f) methanogens, and (h) chitinase-producing microbial communities. The contribution rate of soil physicochemical properties is the sum of the contribution rates of pH, BD, SWC, TC, TN, T, TK, SOM, NO3-N, NH4+-N, and AP. The contribution rate of vegetation is the sum of the contribution rates of ABG, BGB, Sedge−biom, Grass−biom, Forb−biom, Shannon, Simpson, and coverage.

RDA analysis of environmental factors and dominant phyla and genera of methanotrophs revealed a strong positive correlation between NO3-N and Pseudomonadota, a strong positive correlation between pH and Methylomagnum, a strong positive correlation between carbon storage and Methylocystis, and a positive correlation between TP and Methylococcus, but a strong negative correlation between TP and Methylomonas (Figure 7a). RDA analysis of environmental factors and dominant phyla and genera of carbon-fixing microorganisms revealed a strong positive correlation between SWC and the dominant phylum Pseudomonadota of carbon-fixing microorganisms, a strong positive correlation between AP and Cupriavidus, and a strong positive correlation between TP and Methylibium. Additionally, sedge biomass, grass biomass, and total belowground biomass of vegetation were strongly positively correlated with Mesorhizobium but strongly negatively correlated with Pseudomonadota (Figure 7c). Soil TN, TC, and SOM were strongly positively correlated with the dominant phylum Euryarchaeota of methanogens, while vegetation Shannon and Simpson diversity indices were strongly positively correlated with the dominant genera Methanobrevibacter and Methanosarcina of methanogens, respectively, but strongly negatively correlated with the dominant phylum Euryarchaeota of methanogens (Figure 7e). Among plant characteristics, aboveground biomass and vegetation coverage were the main factors influencing the dominant genera Janthinobacterium and Chromobacterium of chitinase-producing microorganisms, and the vegetation Shannon and Simpson indices were strongly negatively correlated with the dominant phylum Pseudomonadota (Figure 7g). Furthermore, among soil physicochemical properties, TK was strongly positively correlated with the dominant phylum Actinomycetota of chitinase-producing microorganisms, TP was strongly positively correlated with the dominant genus Actinoplanes, and soil NH4+-N, TC, TN, and SOM were strongly positively correlated with the dominant genus Chromobacterium (Figure 7g).

Based on the RDA analysis results, the influential factors contributing to methanotrophs diversity were identified as soil physicochemical properties (a contribution of 65.6%), also the main factors affecting methanotrophs diversity (Figure 7b). Vegetation characteristics were the main factors influencing the diversity of carbon-fixing microorganisms, with a contribution rate of 65.6% (Figure 7d). The main factors influencing methanogens diversity were vegetation and soil characteristics, with a contribution rate of 47.7 and 50.4%, respectively (Figure 7f). Plant characteristics were the main factor affecting chitinase-producing microorganisms, with a contribution of 73.0% (Figure 7h).

Random forest modeling reveals microbial predictors of ecosystem carbon storage (Figure 9), and indicates that the diversity of microorganisms involved in the four types of soil carbon metabolism significantly affected soil carbon storage (p < 0.05). Among them, the α diversity indices (Shannon and Simpson indices) and β diversity index of carbon-fixing microorganisms were important predictors of carbon storage of grassland ecosystems. In terms of relative abundance at the genus level, Methylococcus, Methylomonas, Methylocapsa, and Methylocystis were the important taxa among methanotrophs for predicting carbon storage, while Azospirillum, Methylibium, Variovorax, Bradyrhizobium, Cupriavidus, Thiobacillus, and Mesorhizobium were the important taxa among carbon-fixing microorganisms. Additionally, among methanogens, Methanosarcina and Methanoregula were the important groups influencing carbon storage, and Amycolatopsis, Actinoplanes, Chitiniphilus, Cellulomonas, and Stenotrophomonas were the important taxa among chitinase-producing microorganisms for predicting carbon storage.

Figure 9
Structural equation model diagram showing relationships between

Figure 9. Partial Least Squares Path Modeling (PLS-SEM) of the impact of SWC on plants, soil characteristics, ecosystem carbon storage, and carbon metabolism microorganisms among three types of grassland. Solid and dashed lines indicate significant effects (p < 0.05) and non-significant effects (p > 0.05), respectively. The width of the paths indicates the level of significance, and R2 represents the proportion of the explained variance, Q2 represents the predictive relevance (model’s goodness of fit: 0.529).

Figure 9 showed that SWC explained 89.3% of the variation in carbon storage in the alpine grassland. The R2 values for ecosystem carbon storage, methanotrophs diversity, carbon-fixing microorganism diversity, methanogens diversity, and chitinase-producing microorganism diversity were 0.893, 0.571, 0.775, 0.212, and 0.545, respectively, indicating that soil and vegetation characteristics of the grasslands significantly influenced their carbon storage and carbon metabolism microorganism diversity. SWC affected ecosystem carbon storage (path coefficient = −0.923), mainly by changing vegetation characteristics (path coefficient = 0.884) and soil physicochemical properties (path coefficient = 0.867). SWC enhanced the diversity of methanotrophs (path coefficient = 0.721) and reduced the diversity of carbon-fixing microorganisms (path coefficient = −0.471) by increasing the biomass of vegetation. It also indirectly affected the diversity of methanotrophs by influencing soil physicochemical properties (path coefficient = −1.028).

4 Discussion

4.1 Effects of SWC on carbon storage in different alpine grassland ecosystems

The study found significant differences in carbon storage among the three types of alpine grassland ecosystems, with the relativity following the order of alpine wetland (AW) > alpine meadow (AM) > alpine desert (AD) (Figure 2b), an outcome consistent with previous findings (Li Z. et al., 2023; Wang et al., 2022). This study further revealed that surface SWC in the alpine grasslands is a key factor contributing to these differences in carbon storage (Figure 9). SWC influences vegetation carbon content by affecting vegetation indices (vegetation diversity and biomass) and also alters soil carbon storage in alpine grasslands by influencing soil physicochemical properties (Figure 9).

SWC directly affects vegetation carbon storage by altering plant community structure and biomass (Figure 9). A high SWC promotes the growth of sedges with a high carbon content and increases the availability of soil nutrients and moisture, thereby enhancing belowground vegetation biomass and resulting in the highest vegetation carbon storage in AW (Yuan et al., 2019; Zhang et al., 2011). In contrast, the alpine desert with the lowest SWC has a low vegetation biomass and a high proportion of forbs with a low carbon content (Table 2), leading to the lowest vegetation carbon storage (Yuan et al., 2019). Additionally, SWC has a direct and significantly positive effect on soil carbon storage (Figure 9), with high SWC environments being more conducive to soil carbon accumulation (He et al., 2012; Balasubramanian et al., 2020). Therefore, AW with the highest SWC has the highest vegetation and soil carbon storage, while AD with the lowest SWC has the lowest carbon storage.

4.2 Differences in soil physicochemical properties among different grassland types

SOM, NH4+-N, TN and TC in AW were significantly higher than those in AM and AD, while pH was significantly lower (p < 0.05) (Table 2). These findings are consistent with those in previous studies (Liu et al., 2019). This is mainly because increased SWC promotes plant growth and litter input, thereby increasing soil SOM and TC content (Jobbágy and Jackson, 2000). In addition, a high SWC also promotes the nitrogen mineralization process, thereby increasing soil NH4+-N content. It can also reduce nitrogen volatilization and leaching losses, leading to an increase in soil total nitrogen (TN) (Delgado-Baquerizo et al., 2013). Conversely, a higher soil moisture content enhances soil microbial activity, promoting the dissolution and migration of potassium ions, resulting in a decrease in TK content (Chen Y. S. et al., 2022). Additionally, increased SWC may lead to an increase in acidic substances in water (such as carbonic acid formed by the dissolution of carbon dioxide), thereby lowering soil pH (Qian et al., 2023).

4.3 Effects of SWC on the diversity and composition of soil carbon-metabolizing microorganisms

The PLS-SEM model indicated that SWC in different grasslands indirectly affects the diversity of methanotrophs and soil carbon-fixing microbes by influencing soil physicochemical properties (SOM, TK, NH4+-N, TN, TC, pH) and vegetation characteristics (BGB) (Figure 9).

There were significant differences in the α-diversity (richness, Shannon, Simpson indices) of carbon-metabolizing microorganisms among the three types of alpine grassland, with the lowest α-diversity for all four types of microorganisms (methanotrophs, carbon-fixing microbes, methanogens, and chitinase-producing microorganism) observed in AW (Figure 3). Additionally, the β-diversity of these microorganisms in AW also differed significantly from that in AM and AD (Figures 4eh). SWC was a key factor contributing to these differences. The high SWC in AW led to poor soil aeration, significantly inhibiting the activity and diversity of aerobic or facultative anaerobic microorganisms (such as soil carbon-fixing microbes, methanotrophs, and chitinase-producing microorganism) (Delgado-Baquerizo et al., 2016; Tang et al., 2023; Li et al., 2024). Although SWC was positively correlated with the mcrA gene (Wei et al., 2018), the hypoxic environmental conditions in AW may have selectively favored them to become the dominant groups, thereby reducing the overall diversity of methanogens (Narrowe et al., 2017). Previous studies have also found that plant biomass can highly influence soil microbial activity (Hassan et al., 2019; Schultz et al., 2011). SWC significantly increased the diversity of methanotroph communities by increasing vegetation BGB but decreased the diversity of carbon-fixing microbial communities (Figure 9).

This study also found that the four carbon-metabolizing microorganisms had a similar community composition in all grasslands, with Pseudomonadota being the dominant phylum for soil methanotrophs and carbon-fixing microorganisms, Euryarchaeota being the dominant phylum for soil methanogens, and Pseudomonadota and Actinomycetota being the dominant phyla for chitinase-producing microorganisms (Figure 5). However, the relative abundances of key dominant phyla varied among the grasslands, due to the influence of SWC. Specifically, the relative abundance of Pseudomonadota, the dominant phylum for carbon-fixing microorganisms, followed the order of AW > AM > AD (Figure 5b), which was strongly positively correlated with SWC (LEfSe results in Figure 7c). The lowest SWC in AD not only directly limited water supply but also reduced NH4+-N content, belowground biomass, and vegetation cover (Figure 7c). These environmental conditions, together with poor soil nutrients and low vegetation biomass due to the low SWC, jointly inhibited the growth of Pseudomonadota (Rambia et al., 2024). Conversely, the extremely high SWC in AW directly created a hypoxic environment that was highly favorable for the growth and methane production of the dominant phylum Euryarchaeota for methanogens, resulting in its highest relative abundance in AW (Reitschuler et al., 2016). The low SWC in AD led to a low vegetation biomass (Table 2), reduced soil organic input, and poor nutrient availability, while Actinomycetota gained a competitive advantage in this resource-limited environment due to its ability to decompose recalcitrant carbon (Tang et al., 2023). Therefore, SWC not only directly regulates the anaerobic environment required by methanogens (Euryarchaeota) but also affects the relative abundance of dominant phyla in the four carbon-metabolizing microbial communities by influencing water availability, nutrient cycling, and vegetation conditions.

LDA analysis further revealed the effects of SWC on key functional microorganisms (Figure 10). The dominant methanotrophic genera Methylosinus and Methylacidiphilum in AW can only grow using methane as the sole carbon and energy source (Oshkin et al., 2020; Di et al., 2024), indicating that methane concentration produced by the hypoxic environment with a high SWC in AW was higher than that in AM and AD. The strict anaerobic carbon-fixing genus Sulfuricaulis and the methanogenic genus Methanobrevibacter dominated AW due to the strong selection for their growth and functions by the anaerobic conditions created by a high SWC (Spain et al., 2009; Kuo et al., 2024). The distribution patterns of these marker microorganisms demonstrated that SWC influences specific carbon-metabolizing functional microbial communities by regulating soil environmental conditions.

Figure 10
Diagram consisting of two panels labeled A and B. Panel A features a correlation matrix with colored lines connecting variables like pmoA, cbbL, mcrA, and chiA community diversity to soil properties and biomass factors. Panel B is a heatmap showing Pearson correlation coefficients among microbial diversity indices, C storage, and environmental factors. Red and blue indicate positive and negative correlations respectively. Significance levels are indicated by asterisks.

Figure 10. (a) Results of Mantel test showing the relationship between environmental factors and the community structure of soil methanotrophs, carbon-fixing microorganisms, methanogens, and chitinase-producing microorganisms; (b) Correlation coefficient between the diversity of methanotrophs, carbon-fixing microorganisms, methanogens, and chitinase-producing microorganisms and ecosystem carbon storage. * p < 0.05, ** p < 0.01, *** p < 0.001.

4.4 Relationships between soil carbon-metabolizing microorganisms and carbon storage

The diversity and composition of soil carbon-metabolizing microorganisms affect carbon storage in alpine grasslands (Figure 8a). The β-diversity of methanotrophs significantly influenced carbon storage (p < 0.05) (Figure 8a), with key genera such as Methylococcus (Type I), Methylomonas (Type I), Methylocapsa (Type I), and Methylocystis (Type II) directly participating in methane oxidation (Liu et al., 2016; Li et al., 2022). The dominant Methylocystis (Type II) in AW had a higher methane oxidation efficiency, significantly enhancing the carbon fixation capacity of the system (Singh and Tate, 2007). The diversity of carbon-fixing microorganisms had an extremely significant effect on carbon storage (p < 0.01), but its α-diversity was negatively correlated with carbon storage (Figure 10b), possibly reflecting the inhibition of aerobic carbon-fixing microbial activity by the high SWC of AW. Key carbon-fixing microbial genera such as Thiobacillus (directly fixing CO₂), Mesorhizobium and Azospirillum (promoting plant carbon fixation), and Methylibium (oxidizing methane) influenced carbon accumulation through diverse pathways (Komaniecka et al., 2010; Estendorfer et al., 2020). The higher abundances of Mesorhizobium, Azospirillum, and Methylibium in AM than in AD explained why AM had a higher carbon storage than AD (Figure 6b). Methanogenic genera such as Methanosarcina and Methanoregula influenced carbon decomposition through methane production. The Shannon diversity of chitinase-producing microorganism significantly influenced carbon storage (Figure 8a), with key genera such as Amycolatopsis, Actinoplanes, Chitiniphilus, Cellulomonas, and Stenotrophomonas being important taxonomic groups for predicting carbon storage (Figure 8e). Among them, Cellulomonas and Stenotrophomonas had extremely significant effects on carbon storage (p < 0.01). These microorganisms can participate in the carbon cycle by degrading chitin, and an increase in their relative abundance may lead to an increase in carbon decomposition capacity in grasslands, thereby reducing carbon storage (Cai et al., 2023; Zhang et al., 2022). The dominance of Stenotrophomonas, Chitiniphilus, and Amycolatopsis in AD explained why it had the lowest carbon storage. Therefore, the abundance of specific functional microbial groups and their metabolic activities are the core microbiological mechanisms driving differences in carbon storage among different types of alpine grassland.

Figure 8
Grouped bar charts display the increase in mean squared error (MSE) in percentage across various gene markers and bacterial taxa. Chart (a) categorizes data by gene markers: pmoA, cbbL, mcrA, and chiA, with different colors for each. Charts (b), (c), (d), and (e) show MSE increases for different taxa, each identified by name, with percentage increases shown on the x-axis. Stars denote significance levels.

Figure 8. (a) Random forest modeled results of soil carbon storage based on the diversity of carbon metabolism microorganisms and the relative abundance of dominant genera among (b) methanotrophs, (c) carbon-fixing microorganisms, (d) methanogens, and (e) chitinase-producing microorganisms. The percent increase in mean squared error (MSE) of the variables was used to estimate the importance of the predictors, with higher MSE (%) values indicating greater importance. * p < 0.05, ** p < 0.01, *** p < 0.001.

5 Conclusion

This study found that the carbon storage capacity of the three types of alpine grassland ecosystems follows the ranking order of AW > AM > AD. SWC can increase vegetation carbon storage and soil carbon storage by enhancing belowground biomass and total soil carbon content, ultimately affecting the carbon storage of grassland ecosystems. SWC enhances methanogen diversity and reduces soil carbon-fixing microbes diversity by promoting vegetation biomass, and it also influences soil carbon-fixing microbes diversity by affecting soil physicochemical properties. Random forest model revealed that the diversity of the four types of carbon metabolizing microorganisms can all influence grassland carbon storage, among which there is a significant negative correlation between the diversity of carbon-fixing microorganisms and grassland carbon storage. The carbon-fixing microorganisms Thiobacillus, Mesorhizobium, Azospirillum, and Methylibium, as well as the chitinase-producing microorganisms Cellulomonas and Stenotrophomonas, have a highly significant impact on carbon storage.

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

QiaL: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft. WY: Formal analysis, Funding acquisition, Project administration, Resources, Supervision, Validation, Writing – review & editing. JK: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Validation, Writing – review & editing. QinL: Investigation, Visualization, Writing – review & editing. YZ: Investigation, Writing – review & editing. XL: Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing. JZ: Funding acquisition, Project administration, Resources, Supervision, Validation, Writing – review & editing. ZH: Investigation, Writing – review & editing. LC: Investigation, Writing – review & editing. YN: Investigation, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by the National Natural Science Foundation of China [grant numbers U23A20159, U21A20191] and the Qinghai Science and Technology Department [grant number 2023-QY-210].

Conflict of interest

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

The reviewer HL declared a shared affiliation with authors JZ and XL to the handling editor at the time of review.

Generative AI statement

The authors declare that no Gen 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.1627840/full#supplementary-material

References

Bai, Y. F., Wu, J. G., Clark, C. M., Pan, Q. M., Zhang, L. X., Chen, S. P., et al. (2012). Grazing alters ecosystem functioning and C:N:P stoichiometry of grasslands along a regional precipitation gradient. J. Appl. Ecol. 49, 1204–1215. doi: 10.1111/j.1365-2664.2012.02205.x

Crossref Full Text | Google Scholar

Balasubramanian, D., Zhou, W. J., Ji, H. L., Grace, J., Bai, X. L., Song, Q. H., et al. (2020). Environmental and management controls of soil carbon storage in grasslands of southwestern China. J. Environ. Manag. 254:109810. doi: 10.1016/j.jenvman.2019.109810

PubMed Abstract | Crossref Full Text | Google Scholar

Berdugo, M., Maestre, F. T., Kéfi, S., Gross, N., Le Bagousse-Pinguet, Y., and Soliveres, S. (2019). Aridity preferences alter the relative importance of abiotic and biotic drivers on plant species abundance in global drylands. J. Ecol. 107, 190–202. doi: 10.1111/1365-2745.13006

Crossref Full Text | Google Scholar

Bourne, D. G., McDonald, I. R., and Murrell, J. C. (2001). Comparison of pmoA PCR primer sets as tools for investigating methanotroph diversity in three Danish soils. Appl. Environ. Microbiol. 67, 3802–3809. doi: 10.1128/aem.67.9.3802-3809

Crossref Full Text | Google Scholar

Brankatschk, R., Töwe, S., Kleineidam, K., Schloter, M., and Zeyer, J. (2011). Abundances and potential activities of nitrogen cycling microbial communities along a chronosequence of a glacier forefield. ISME J. 5, 1025–1037. doi: 10.1038/ismej.2010.184

PubMed Abstract | Crossref Full Text | Google Scholar

Cai, R., Zuo, S., Cao, X., and Xu, C. (2023). Carbon and nitrogen transformation mechanism and functional microorganisms of sheep manure composting on the Qinghai-Tibet plateau under different moisture content. J. Environ. Chem. Eng. 11:111341. doi: 10.1016/j.jece.2023.111341

Crossref Full Text | Google Scholar

Canedoli, C., Ferrè, C., Abu El Khair, D., Comolli, R., Liga, C., Mazzucchelli, F., et al. (2020). Evaluation of ecosystem services in a protected mountain area: soil organic carbon stock and biodiversity in alpine forests and grasslands. Ecosyst. Serv. 44:101135. doi: 10.1016/j.ecoser.2020.101135

Crossref Full Text | Google Scholar

Chen, H., Ju, P. J., Zhu, Q., Xu, X. L., Wu, N., Gao, Y. H., et al. (2022). Carbon and nitrogen cycling on the Qinghai-Tibetan plateau. Nat. Rev. Earth Environ. 3, 701–716. doi: 10.1038/s43017-022-00344-2

Crossref Full Text | Google Scholar

Chen, Y. S., Zhang, S. H., Li, H. L., and Wang, Y. D. (2022). Drivers of nutrient content and spatial variability of soil multifunctionality in the topsoil of Kyrgyzstan. Front. Environ. Sci. 10:1001984. doi: 10.3389/fenvs.2022.1001984

Crossref Full Text | Google Scholar

Cleland, E. E., Collins, S. L., Dickson, T. L., Farrer, E. C., Gross, K. L., Gherardi, L. A., et al. (2013). Sensitivity of grassland plant community composition to spatial vs. temporal variation in precipitation. Ecology 94, 1687–1696. doi: 10.1890/12-1006.1

PubMed Abstract | Crossref Full Text | Google Scholar

Delgado-Baquerizo, M., Maestre, F. T., Gallardol, A., Bowker, M. A., Wallenstein, M. D., Quero, J. L., et al. (2013). Decoupling of soil nutrient cycles as a function of aridity in global drylands. Nature 502, 672–676. doi: 10.1038/nature12670

Crossref Full Text | Google Scholar

Delgado-Baquerizo, M., Maestre, F. T., Reich, P. B., Trivedi, P., Osanai, Y., Liu, Y. R., et al. (2016). Carbon content and climate variability drive global soil bacterial diversity patterns. Ecol. Monogr. 86, 373–390. doi: 10.1002/ecm.1216

Crossref Full Text | Google Scholar

Deng, L., Shangguan, Z. P., Wu, G. L., and Chang, X. F. (2017). Effects of grazing exclusion on carbon sequestration in China's grassland. Earth-Sci. Rev. 173, 84–95. doi: 10.1016/j.earscirev.2017.08.008

Crossref Full Text | Google Scholar

Di, X., Li, P., Wang, J. X., Nowak, V., Zhi, S., Jin, M. L., et al. (2024). Genome mining analysis uncovers the previously unknown biosynthetic capacity for secondary metabolites in Verrucomicrobia. Mar. Biotechnol. 26, 1324–1335. doi: 10.1007/s10126-024-10374-0

PubMed Abstract | Crossref Full Text | Google Scholar

Du, C. J., and Gao, Y. H. (2020). Opposite patterns of soil organic and inorganic carbon along a climate gradient in the alpine steppe of northern Tibetan plateau. Catena 186:104366. doi: 10.1016/j.catena.2019.104366

Crossref Full Text | Google Scholar

Estendorfer, J., Stempfhuber, B., Vestergaard, G., Schulz, S., Rillig, M. C., Joshi, J., et al. (2020). Definition of core bacterial taxa in different root compartments of Dactylis glomerata, grown in soil under different levels of land use intensity. Diversity 12:392. doi: 10.3390/d12100392

Crossref Full Text | Google Scholar

Fu, B. J., Liu, S. L., Chen, L. D., Lü, Y. H., and Qiu, J. (2004). Soil quality regime in relation to land cover and slope position across a highly modified slope landscape. Ecol. Res. 19, 111–118. doi: 10.1111/j.1440-1703.2003.00614.x

Crossref Full Text | Google Scholar

Ganjurjav, H., Gao, Q. Z., Gornish, E. S., Schwartz, M. W., Liang, Y., Cao, X. J., et al. (2016). Differential response of alpine steppe and alpine meadow to climate warming in the Central Qinghai-Tibetan plateau. Agric. For. Meteorol. 223, 233–240. doi: 10.1016/j.agrformet.2016.03.017

Crossref Full Text | Google Scholar

Ge, T. D., Wu, X. H., Chen, X. J., Yuan, H. Z., Zou, Z. Y., Li, B. Z., et al. (2013). Microbial phototrophic fixation of atmospheric CO2 in China subtropical upland and paddy soils. Geochim. Cosmochim. Acta 113, 70–78. doi: 10.1016/j.gca.2013.03.020

Crossref Full Text | Google Scholar

Gilmanov, T. G., Soussana, J. F., Aires, L., Allard, V., Ammann, C., Balzarolo, M., et al. (2007). Partitioning European grassland net ecosystem CO2 exchange into gross primary productivity and ecosystem respiration using light response function analysis. Agric. Ecosyst. Environ. 121, 93–120. doi: 10.1016/j.agee.2006.12.008

Crossref Full Text | Google Scholar

Hair, J. F., Risher, J. J., Sarstedt, M., and Ringle, C. M. (2019). When to use and how to report the results of PLS-SEM. Eur. Bus. Rev. 31, 2–24. doi: 10.1108/EBR-11-2018-0203

Crossref Full Text | Google Scholar

Hassan, M. K., McInroy, J. A., and Kloepper, J. W. (2019). The interactions of rhizodeposits with plant growth-promoting rhizobacteria in the rhizosphere: a review. Agriculture 9:142. doi: 10.3390/agriculture9070142

Crossref Full Text | Google Scholar

He, N. P., Chen, Q. S., Han, X. G., Yu, G. R., and Li, L. H. (2012). Warming and increased precipitation individually influence soil carbon sequestration of inner Mongolian grasslands, China. Agric. Ecosyst. Environ. 158, 184–191. doi: 10.1016/j.agee.2012.06.010

Crossref Full Text | Google Scholar

Jobbágy, E. G., and Jackson, R. B. (2000). The vertical distribution of soil organic carbon and its relation to climate and vegetation. Ecol. Appl. 10, 423–436. doi: 10.2307/2641104

Crossref Full Text | Google Scholar

Jones, D. L., Olivera-Ardid, S., Klumpp, E., Knief, C., Huil, P. W., Lehndorff, E., et al. (2018). Moisture activation and carbon use efficiency of soil microbial communities along an aridity gradient in the Atacama Desert. Soil Biol. Biochem. 117, 68–71. doi: 10.1016/j.soilbio.2017.10.026

Crossref Full Text | Google Scholar

Jöreskog, K. G., and Goldberger, A. S. (1975). Estimation of a model with multiple indicators and multiple causes of a single latent variable. J. Am. Stat. Assoc. 70, 631–639. doi: 10.2307/2285946

Crossref Full Text | Google Scholar

Juottonen, H., Galand, P. E., and Yrjälä, K. (2006). Detection of methanogenic Archaea in peat: comparison of PCR primers targeting the mcr a gene. Res. Microbiol. 157, 914–921. doi: 10.1016/j.resmic.2006.08.006

PubMed Abstract | Crossref Full Text | Google Scholar

Kemp, P. F., and Aller, J. Y. (2004). Bacterial diversity in aquatic and other environments: what 16S rDNA libraries can tell us. FEMS Microbiol. Ecol. 47, 161–177. doi: 10.1016/s0168-6496(03)00257-5

PubMed Abstract | Crossref Full Text | Google Scholar

Komaniecka, I., Zdzisinska, B., Kandefer-Szerszen, M., and Choma, A. (2010). Low endotoxic activity of lipopolysaccharides isolated from Bradyrhizobium, Mesorhizobium, and Azospirillum strains. Microbiol. Immunol. 54, 717–725. doi: 10.1111/j.1348-0421.2010.00269.x

PubMed Abstract | Crossref Full Text | Google Scholar

Kumar, S., Lal, R., Liu, D. S., Kumar, S., Lal, R., and Liu, D. A. (2012). A geographically weighted regression kriging approach for mapping soil organic carbon stock. Geoderma. 8. doi: 10.1016/j.geoderma.2012.05.022

Crossref Full Text | Google Scholar

Kuo, J., Liu, D., Wen, W. H., Chiu, C. Y., Chen, W. Y., Wu, Y. W., et al. (2024). Different microbial communities in paddy soils under organic and nonorganic farming. Braz. J. Microbiol. 55, 777–788. doi: 10.1007/s42770-023-01218-5

PubMed Abstract | Crossref Full Text | Google Scholar

Liang, C., Schimel, J. P., and Jastrow, J. D. (2017). The importance of anabolism in microbial control over soil carbon storage. Nat. Microbiol. 2:17105. doi: 10.1038/nmicrobiol.2017.105

PubMed Abstract | Crossref Full Text | Google Scholar

Li, C. Y., Li, X. L., Shi, Y., Yang, Y. W., and Li, H. L. (2022). Effects of nitrogen addition on soil carbon-fixing microbial diversity on different slopes in a degraded alpine meadow. Front. Plant Sci. 13:921278. doi: 10.3389/fpls.2022.921278

PubMed Abstract | Crossref Full Text | Google Scholar

Li, N., Wang, B. R., Zhou, Y., Li, H. J., Zhu, Z. L., Dou, Y. X., et al. (2024). Response of the C-fixing bacteria community to precipitation changes and its impact on bacterial necromass accumulation in semiarid grassland. J. Environ. Manag. 354:120289. doi: 10.1016/j.jenvman.2024.120289

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, G. M., Zhang, X. L., Wu, T. H., Wu, X. L., Smoak, J. M., Li, X. X., et al. (2019). Seasonal changes in labile organic matter as a function of environmental factors in a relict permafrost region on the Qinghai-Tibetan plateau. Catena 180, 194–202. doi: 10.1016/j.catena.2019.04.026

Crossref Full Text | Google Scholar

Liu, S. B., Zamanian, K., Schleuss, P. M., Zarebanadkouki, M., and Kuzyakov, Y. (2018). Degradation of Tibetan grasslands: consequences for carbon and nutrient cycles. Agric. Ecosyst. Environ. 252, 93–104. doi: 10.1016/j.agee.2017.10.011

Crossref Full Text | Google Scholar

Liu, S. L., Du, Y. G., Zhang, F. W., Lin, L., Li, Y. K., Guo, X. W., et al. (2016). Distribution of soil carbon in different grassland types of the Qinghai-Tibetan plateau. J. Mt. Sci. 13, 1806–1817. doi: 10.1007/s11629-015-3764-7

Crossref Full Text | Google Scholar

Li, X. B., Bai, Y. X., Wen, W. Y., Wang, H., Li, R. H., Li, G. Q., et al. (2017). Effects of grassland degradation and precipitation on carbon storage distributions in a semi-arid temperate grassland of Inner Mongolia, China. Acta Oecol. 85, 44–52. doi: 10.1016/j.actao.2017.09.008

Crossref Full Text | Google Scholar

Li, Y., Xiong, L., Zeng, K., Wei, Y., Li, H., and Ji, X. (2023). Microbial-driven carbon fixation in natural wetland. J. Basic Microbiol. 63, 1115–1127. doi: 10.1002/jobm.202300273

PubMed Abstract | Crossref Full Text | Google Scholar

Li, Z., Yang, Y., Zheng, H., Hu, B., Dai, X., Meng, N., et al. (2023). Environmental changes drive soil microbial community assembly across arid alpine grasslands on the Qinghai-Tibetan plateau, China. Catena 228:107175. doi: 10.1016/j.catena.2023.107175

Crossref Full Text | Google Scholar

Lu, S. S., Zhou, N. Q., Cai, Y., Guo, M. S., Sheng, D., and Liu, Y. (2024). Ecological stoichiometry of carbon, nitrogen and phosphorus in deep sediments within the critical zone of south dongting lake wetland, China. J. Soil Sci. Plant Nutr. 24, 6389–6403. doi: 10.1007/s42729-024-01976-7

Crossref Full Text | Google Scholar

Nachtergaele, F., Velthuizen, H., Verelst, L., Wiberg, D., Batjes, N., Dijkshoorn, K., et al. (2012). Harmonized world soil database (version 1.2). Laxenburg, Austria: FAO, Rome, Italy and IIASA.

Google Scholar

Narrowe, A. B., Angle, J. C., Daly, R. A., Stefanik, K. C., Wrighton, K. C., and Miller, C. S. (2017). High-resolution sequencing reveals unexplored archaeal diversity in freshwater wetland soils. Environ. Microbiol. 19, 2192–2209. doi: 10.1111/1462-2920.13703

PubMed Abstract | Crossref Full Text | Google Scholar

Oshkin, I. Y., Miroshnikov, K. K., Grouzdev, D. S., and Dedysh, S. N. (2020). Pan-genome-based analysis as a framework for demarcating two closely related methanotroph menera Methylocystis and Methylosinus. Microorganisms 8:768. doi: 10.3390/microorganisms8050768

PubMed Abstract | Crossref Full Text | Google Scholar

Philippot, L., Spor, A., Hénault, C., Bru, D., Bizouard, F., Jones, C. M., et al. (2013). Loss in microbial diversity affects nitrogen cycling in soil. ISME J. 7, 1609–1619. doi: 10.1038/ismej.2013.34

PubMed Abstract | Crossref Full Text | Google Scholar

Prober, S. M., Leff, J. W., Bates, S. T., Borer, E. T., Firn, J., Harpole, W. S., et al. (2015). Plant diversity predicts beta but not alpha diversity of soil microbes across grasslands worldwide. Ecol. Lett. 18, 85–95. doi: 10.1111/ele.12381

PubMed Abstract | Crossref Full Text | Google Scholar

Qian, F. K., Yu, Y. J., Dong, X. R., and Gu, H. L. (2023). Soil quality evaluation based on a minimum data set (MDS)-a case study of Tieling county, Northeast China. Land 12:1263. doi: 10.3390/land12061263

Crossref Full Text | Google Scholar

Rambia, A., Veluchamy, C., Rawat, J. M., Jamdhade, M. D., Purohit, S., Pawar, K. D., et al. (2024). Revealing bacterial and fungal communities of the untapped forest and alpine grassland zones of the Western-Himalayan region. Int. Microbiol. 27, 781–795. doi: 10.1007/s10123-023-00430-5

PubMed Abstract | Crossref Full Text | Google Scholar

Reitschuler, C., Hofmann, K., and Illmer, P. (2016). Abundances, diversity and seasonality of (non-extremophilic) archaea in alpine freshwaters. Antonie Van Leeuwenhoek 109, 855–868. doi: 10.1007/s10482-016-0685-6

PubMed Abstract | Crossref Full Text | Google Scholar

Schultz, R., Andrews, S., O’Reilly, L., Bouchard, V., and Frey, S. (2011). Plant community composition more predictive than diversity of carbon cycling in freshwater wetlands. Wetlands 31, 965–977. doi: 10.1007/s13157-011-0211-6

Crossref Full Text | Google Scholar

Singh, B. K., and Tate, K. (2007). Biochemical and molecular characterization of methanotrophs in soil from a pristine New Zealand beech forest. FEMS Microbiol. Lett. 275, 89–97. doi: 10.1111/j.1574-6968.2007.00885.x

PubMed Abstract | Crossref Full Text | Google Scholar

Spain, A. M., Krumholz, L. R., and Elshahed, M. S. (2009). Abundance, composition, diversity and novelty of soil Proteobacteria. ISME J. 3, 992–1000. doi: 10.1038/ismej.2009.43

PubMed Abstract | Crossref Full Text | Google Scholar

Sun, J., Ye, C. C., Liu, M., Wang, Y., Chen, J., Wang, S., et al. (2022). Response of net reduction rate in vegetation carbon uptake to climate change across a unique gradient zone on the Tibetan plateau. Environ. Res. 203:111894. doi: 10.1016/j.envres.2021.111894

Crossref Full Text | Google Scholar

Sun, J., Zhou, T. C., Liu, M., Chen, Y. C., Liu, G. H., Xu, M., et al. (2020). Water and heat availability are drivers of the aboveground plant carbon accumulation rate in alpine grasslands on the Tibetan plateau. Glob. Ecol. Biogeogr. 29, 50–64. doi: 10.1111/geb.13006

Crossref Full Text | Google Scholar

Tang, X., Zhang, M., Fang, Z., Yang, Q., Zhang, W., Zhou, J., et al. (2023). Changing microbiome community structure and functional potential during permafrost thawing on the Tibetan plateau. FEMS Microbiol. Ecol. 99:fiad117. doi: 10.1093/femsec/fiad117

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, C., Wang, G., Wang, Y., Rafique, R., Ma, L., Hu, L., et al. (2016). Fire alters vegetation and soil microbial community in alpine meadow. Land Degrad. Dev. 27, 1379–1390. doi: 10.1002/ldr.2367

Crossref Full Text | Google Scholar

Wang, Q. Q., Huang, Q., Liu, Y., Khan, M. A., Guo, G. M., Lu, W. K., et al. (2022). Temporal variation of the coupling relationship between methanogens and biogeochemical process in soil-microbes-rice system. Chemosphere 303:135099. doi: 10.1016/j.chemosphere.2022.135099

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, S. L., Chen, X. D., Li, W., Gong, W. L., Wang, Z. W., and Cao, W. X. (2023). Grazing exclusion alters soil methane flux and methanotrophic and methanogenic communities in alpine meadows on the Qinghai-Tibet plateau. Front. Microbiol. 14:1293720. doi: 10.3389/fmicb.2023.1293720

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, X. T., Zhao, N. N., Li, W. C., Pu, X., Xu, P., and Wang, P. C. (2024). Core bacterial taxa determine formation of forage yield in fertilized soil. Microorganisms 12:1679. doi: 10.3390/microorganisms12081679

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Y., Xiao, J., Ma, Y., Luo, Y., Hu, Z., Li, F., et al. (2021). Carbon fluxes and environmental controls across different alpine grassland types on the Tibetan plateau. Agric. For. Meteorol. 311:108694. doi: 10.1016/j.agrformet.2021.108694

Crossref Full Text | Google Scholar

Wei, S. P., Cui, H. P., Zhu, Y. H., Lu, Z. Q., Pang, S. J., Zhang, S., et al. (2018). Shifts of methanogenic communities in response to permafrost thaw results in rising methane emissions and soil property changes. Extremophiles 22, 447–459. doi: 10.1007/s00792-018-1007-x

PubMed Abstract | Crossref Full Text | Google Scholar

Xiao, X., Yin, X. B., Lin, H., Sun, L. G., You, Z. Y., Wang, P., et al. (2005). Chitinase genes in lake sediments of Ardley Island, Antarctica. Appl. Environ. Microbiol. 71, 7904–7909. doi: 10.1128/aem.71.12.7904-7909.2005

Crossref Full Text | Google Scholar

Yang, W. S., Liu, Y., Zhao, J. X., Chang, X. F., Wiesmeier, M., Sun, J., et al. (2021). SOC changes were more sensitive in alpine grasslands than in temperate grasslands during grassland transformation in China: a meta-analysis. J. Clean. Prod. 308:127430. doi: 10.1016/j.jclepro.2021.127430

Crossref Full Text | Google Scholar

Yuan, Z. Q., Jiang, X. J., Liu, G. J., Jin, H. J., Chen, J., and Wu, Q. B. (2019). Responses of soil organic carbon and nutrient stocks to human-induced grassland degradation in a Tibetan alpine meadow. Catena 178, 40–48. doi: 10.1016/j.catena.2019.03.001

Crossref Full Text | Google Scholar

Zhang, G. G., Kang, Y. M., Han, G. D., Mei, H., and Sakurai, K. (2011). Grassland degradation reduces the carbon sequestration capacity of the vegetation and enhances the soil carbon and nitrogen loss. Acta Agric. Scand. Sect. B - Soil Plant Sci. 61, 356–364. doi: 10.1080/09064710.2010.495079

Crossref Full Text | Google Scholar

Zhang, Y. L., Sun, C. X., Wang, S. Q., Xie, H. T., Jiang, N., Chen, Z. H., et al. (2022). Stover and biochar can improve soil microbial necromass carbon, and enzymatic transformation at the genetic level. GCB Bioenergy 14, 1082–1096. doi: 10.1111/gcbb.12984

Crossref Full Text | Google Scholar

Zhou, Z. C., Han, P., and Gu, J. D. (2014). New PCR primers based on mcrA gene for retrieving more anaerobic methanotrophic archaea from coastal reedbed sediments. Appl. Microbiol. Biotechnol. 98, 4663–4670. doi: 10.1007/s00253-014-5599-5

Crossref Full Text | Google Scholar

Keywords: soil carbon functional microbial, carbon storage, alpine grassland, Tibetan Plateau, soil water content

Citation: Liu Q, Yang W, Kou J, Li Q, Zhang Y, Li X, Zhang J, Hao Z, Chi L and Ning Y (2025) Ecosystem carbon storage and carbon metabolizing microorganisms in three types of grasslands on the Qinghai-Tibet Plateau. Front. Microbiol. 16:1627840. doi: 10.3389/fmicb.2025.1627840

Received: 13 May 2025; Accepted: 14 July 2025;
Published: 25 July 2025.

Edited by:

Reeta Goel, G. B. Pant University of Agriculture and Technology, India

Reviewed by:

Honglin Li, Qinghai University, China
Xinwei Li, Qinghai University, China
Zhang Minghui, Beihua University, China

Copyright © 2025 Liu, Yang, Kou, Li, Zhang, Li, Zhang, Hao, Chi and Ning. 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: Jiancun Kou, amlhbmN1bjAyQDE2My5jb20=

Disclaimer: 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.