- 1Qinghai Academy of Animal and Veterinary Science, Academy of Animal and Veterinary Science, Qinghai University, Xinin, Qinghai, China
- 2Qinghai University Qinghai Provincial Key Laboratory of Adaptive Management on Alpine Grassland, Xining, Qinghai, China
Introduction: Soil properties and extracellular enzyme activities are the key factors affecting microbial communities at different altitudes. It is very important to understand their distribution patterns along the elevation gradient for predicting the function of alpine ecosystems in response to climate change, especially in alpine shrub ecosystems, which are rarely studied.
Methods: In eastern Qinghai-Tibet Plateau Set up four sample plots with different altitudes (3400,3700,4000,4300 and 4300 m), Illumina I6S gene sequencing was used to analyze the changing law of soil bacterial community, and its diversity and influencing factors were discussed by combining correlation and redundancy analysis.
Results: The results showed that the contents of soil total carbon, total nitrogen, water content, microbial biomass carbon, nitrogen, phosphorus and alkaline phosphatase were the highest, but the bacterial diversity was the lowest (Shannon and Pielou index); Dominant bacteria are different at different altitude gradients. Acidobacteriota is the most abundant at 3700, 4000 and 4300 m above sea level, while Verrucomicrobiota is the most abundant at 3700m above sea level; In addition, PCoA analysis showed that the community structure of soil bacteria changed significantly, with similar structures at 3700 , 4000 and 4300 m above sea level; Redundancy analysis (RDA) showed that soil properties (TC and BD) and enzyme activities (AG, BG and CBH) were the main factors affecting the bacterial community, while soil properties (TN), microbial biomass (MBC and MBN) and enzyme activities (AG and BG) all had significant effects on the functional groups of soil bacteria.
Conclusion: To sum up, these results show that soil physical and chemical properties, microbial biomass, enzyme activities and bacterial communities have different responses to different altitude gradients.These studies provide a new perspective for us to deeply understand the driving factors of soil bacterial community variation along the elevation gradient. It is helpful to strengthen our awareness of the protection of Potentilla fruticosa shrub, and also provide basic information for the study of alpine ecosystem in this area.
1 Introduction
Soil microorganisms and enzymes are important components of soil and play a significant role in the ecosystem (Shu et al., 2023). They are not only regulators of ecosystem functions but also the main drivers of nutrient transformation and material cycling in the ecosystem (Wang et al., 2025). The main groups of soil microorganisms include bacteria, fungi, actinomycetes, and protozoa. These microorganisms can secrete enzymes such as alkaline phosphatase, β-1,4-glucosidase, cellulase, amidase, and nitrogenase (Arunrat et al., 2024). Through decomposition, they can mineralize soil organic matter, including remains of animals and plants, and release nutrients into the surrounding soil, playing an important role in the normal growth of plants and the stability of ecosystem structure and function (Wyszkowska et al., 2019). As the main group of soil microorganisms, soil bacteria have complex ecological functions and promote the decomposition of soil organic matter by participating in the cycle and metabolism of soil nutrients (Yan et al., 2023). Soil moisture, organic carbon, and total nitrogen are generally important environmental factors that influence the composition of soil bacterial communities (Pang et al., 2023). In addition, changes in the environment of soil ecological processes can lead to significant changes in soil bacterial communities at the regional scale (Crowther et al., 2019). Especially in alpine shrub ecosystems, the distribution pattern of soil bacterial communities along altitudinal gradients has always been a hot topic in ecological research (Zhu et al., 2025). Therefore, analyzing the soil bacterial communities in the P.fruticosa shrubs of the Yellow River source area on the Qinghai-Tibet Plateau is of great significance for the biodiversity conservation and ecological management of alpine ecosystems.
In the alpine shrub ecosystem, there are significant differences in climate, plant and soil characteristics in a short spatial distance (Bayranvand et al., 2021; Zhao et al., 2025). Altitude gradient is considered as a powerful “natural laboratory” to study the response and feedback of soil function to climate change (Feng et al., 2021). Especially the altitude model to detect soil microbial community and extracellular enzyme activity, because the climate and biological properties change dramatically in a small spatial scale (Klimek et al., 2020). The variation patterns of soil bacterial community composition, biomass and soil extracellular enzyme activity along altitude gradient were studied extensively (Ren et al., 2018). For example, some studies show that microbial biomass increases, decreases or has no obvious trend with the elevation (Chang et al., 2016; Klimek et al., 2020). In Nyainqentanglha Mountain (Si et al., 2014) and Italian Alps (Margesin et al., 2014), it was found that soil extracellular enzyme activity decreased along the elevation gradient, while in south-central Chile (Reyes et al., 2011) and Wuyishan (Jin et al., 2011), it was found that soil enzyme activity increased with altitude. A study of Kohala volcano in Hawaii by Peay et al. (2017) found that the bacterial diversity showed a hump-shaped upward trend. In a study in Mount Fuji, Singh et al. (2012) found that the diversity of soil bacteria and the abundance of Acinetobacter decreased monotonously with the increase of altitude gradient. Cui et al. (2019) observed the different changes of α-diversity of soil bacteria in Qinghai-Tibet Plateau at high altitude gradient, and the results showed that Shannon-Wiener diversity of bacteria first increased and then decreased. The above research found that the research results of bacterial community and extracellular enzyme activity were inconsistent, and the consensus was still difficult to reach. This shows that the influencing factors are not only the temperature fluctuation caused by altitude rise, but these contradictory observations may be attributed to the interaction of various abiotic and biological factors on bacterial communities and extracellular enzyme activities (Chang et al., 2016; Siles et al., 2017; Zhao et al., 2025). Because this area is very sensitive to climate change, there is still a lack of synchronous understanding of the changes of soil bacterial community and extracellular enzyme activity along the elevation gradient, which needs to be solved urgently (Fan et al., 2021).
At present, the research on the response of soil bacteria to altitude gradient mainly focuses on large-scale areas and different vegetation types (Yuan et al., 2014). However, the mechanism behind the small-scale changes of soil bacterial community structure in the same vegetation type, especially in the P. fruticosa shrub, is limited. In view of that fact that the response of soil microbial community structure to altitude gradient change has become an important part of predict the response of ecosystem to environmental change (Cui et al., 2019), It is necessary to study the response mode and driving factors of soil microorganisms in shrub communities with the same dominant species to altitude gradient changes, so as to improve people’s awareness. P. fruticosa shrub is widely distributed on shady slopes, semi-shady slopes and diluvial fans in mountainous areas at an altitude of 3 200~4 500 m, and is an important shrub species in alpine ecosystem, which plays an important role in maintaining the function of alpine meadow ecosystem in Qinghai-Tibet Plateau (Li et al., 2010).
Therefore, in this study,P. fruticosa shrub investigated along different altitudes in alpine shrub ecosystem was taken as the research object, and Illumina MiSeq technique was used to analyze the composition of soil bacterial community and extracellular enzyme activities related to C, N and P acquisition in the meadow of P.fruticosa shrub at altitude. By studying the important correlation between soil bacterial community and extracellular enzyme activities and soil variables, Characterizing changes in soil bacterial communities and extracellular enzyme activities along altitude, and identifying possible drivers of changes in soil bacterial communities and extracellular enzyme activities along altitude are of great significance in guiding ecological restoration and regeneration of Source region of the Yellow River in the eastern part of the Qinghai-Tibetan Plateau and realize the sustainable management of shrub meadow. We assume that: (i) at high altitude (3700 m), due to the combined influence of soil variables, soil microbial biomass and extracellular enzyme activities related to N and P cycles increased significantly, while extracellular enzyme activities related to carbon cycle decreased significantly; (ii) α diversity and composition of soil bacterial communities will change with altitude. (iii) Soil factor (TC) will be an important factor affecting the bacterial community, because it can affect the soil microbial community by driving C and nutrients into the soil microorganisms during the nutrient flow.
2 Materials and methods
2.1 Study area
Golog Tibetan Autonomous Prefecture in Qinghai Province is located in the hinterland of the Qinghai-Tibet Plateau (96°54′∼101°51′E,32°31′∼35°37′N), at the Source region of the Yellow River, with an average elevation of over 4,200 m, with an annual average temperature of-0.4°C∼3.7°C, annual precipitation of 400-760 mm and annual sunshine hours of 1,988. The vegetation type is alpine meadow and the soil is mainly alpine shrub meadow soil (Xie L. L. et al., 2024). Common herbs are Cyperaceae, such asCarex capillifolia, Carex alatauensis and Carex atrofusca Schkuhr; Gramineae includes Elymus nutans Griseb, Poa pratensis L. and Deschampsia cespitosa; Miscellaneous grasses include Bistorta vivipara, Pedicularis kansuensis Maxim, and Anaphalis lactea Maxim.
2.2 Experimental design and sample collection
In the study area, the experimental sites were selected according to the following principles: (1) The slope and aspect were consistent, and P. fruticosa was the dominant shrub species in the sample plot. (2) In all the experimental sites, the growth status of P. fruticosa shrub is similar; (3) The environmental conditions of each experimental site are relatively consistent, that is, the topography and soil types are basically the same; (4) Based on the above principles, this study selected four experimental sites (3400, 3700, 4000, and 4300 m) with a span of 300 m. The detailed distribution is shown in Figure 1.

Figure 1. Distribution (A) and sampling method (B) of four experimental sites in Guoluo Tibetan Autonomous Prefecture.
The sampling time of this study is August, 2023, the peak season of forage growth, in the altitude range of 3400 -4300 m, taking 300 m as an altitude span to illustrate the difference of different altitude gradients. In order to reduce the spatial heterogeneity, we set up 6 repeated plots in each experimental site, and the interval between each plot is 10 m. Therefore, we have a total of 24 plots (4 experimental sites and 6 repeated plots). Field quadrat survey was used to record the species composition and coverage of each plot, and the latitude, longitude and altitude of each experimental point were recorded by GPS. A total of 24 soil samples were collected in this study (Bai et al., 2020).
At each elevation gradient, six 10m×10m shrub quadrats were established. Within each quadrat, five sampling points were set up using the five-point sampling method. After removing surface litter and stones from each sampling point, soil samples were collected from the 0-20 cm depth using a soil auger with a diameter of 3.5 cm. The soil samples from the five sampling points within each quadrat were thoroughly mixed to form a composite soil sample. Therefore, a total of 24 soil samples were collected in this study (4 elevation gradients × 6 replicates). Pass the soil sample through a 2 mm sieve to remove roots and stones, and then the soil samples are immediately sent back to the cooler in the laboratory. The soil sample is divided into two parts, which are used to determine the chemical properties and microorganisms existing in the soil. A portion soil samples are naturally air-dried indoors to determine soil physical and chemical properties and other indicators. The other part was stored at -80°C for Qualcomm gene sequencing (Zhou et al., 2024).
2.3 Analysis of soil nutrients and enzyme activities
The contents of total nitrogen (TC) and total nitrogen (TN) in soil were determined by C and N analyzer (Elementar, Langenselbold, Germany);Total phosphorus (TP) in soil was determined by molybdenum antimony colorimetric method and flow analyzer (Proxima, AMS Alliance, Paris, France); pH measurement: Measure the soil pH from the sample containing the ratio of soil (air-dried) to water of 1:5; Soil moisture content (SM): measure the soil moisture content after drying in an oven at 105°C for 48 h to constant weight; BD: The density (BD) of soil samples was measured by circular knife method and dried at 105°C for 24 hours. Soil microbial biomass carbon (MBC) and microbial biomass nitrogen (MBN) were determined by chloroform steaming-K2SO4 leaching method, and soil microbial biomass phosphorus (MBP) was determined by chloroform steaming-NaHCO3 leaching method (Zhao et al., 2023).
Six important enzymes in soil carbon, nitrogen and phosphorus cycling were selected to determine their activities, among which β-D-glucosidase (BG), cellulose-degrading cellobiohydrolase (CBH) and α-1,4-glucosidase (AG) were C-harvesting enzymes (Cenini et al., 2015). β-1,4-N-acetylglucosamine glucosidase (NAG) and Leucine aminopeptidase (LAP) are common acquisition enzymes (Salazar et al., 2020). Alkaline phosphatase (ALP) is p-acquiring enzyme (Delgado-Baquerizo et al., 2017). According to the method of Saiya-Cork et al. (2002) the activity of extracellular enzymes in soil was determined: firstly, 1g of soil was weighed and put into a jar filled with 125ml acetic acid buffer with pH 5.05 mmol/L, and stirred for 5min with a magnetic stirrer to homogenize the soil and solution; Then use a pipette to take 200μL of suspension. Add it to a 96-well plate, put it in the dark at 20°C for 4 hours, and then add 10μL 1 mol/L NaOH solution; Finally, the soil enzyme activity was determined by a multifunctional enzyme marker (Synergy H, USA) (Saiya-Cork et al., 2002).
2.4 16S amplicon sequencing method and analysis
DNA was extracted from 0.5g of soil sample using the Hipure Soil DNA Kit (Magen, Guangzhou, China). 2µL of the DNA sample was taken and the OD value of the nucleic acid was measured using a NanoDrop 2000 micro-spectrophotometer to detect the purity of the nucleic acid. Agarose (1% agarose) gel electrophoresis was used to detect the integrity of the nucleic acid sample and the degree of protein contamination. 2µL of the DNA sample was taken and the DNA concentration of each sample was detected using Qubit fluorescence quantification (Qubit 3.0, Thermo Fisher, USA), and the total DNA amount of the sample was calculated based on the concentration quantification results. Specific primers with barcodes were used to amplify the V4 region of the bacterial 16S rRNA, and the primers used were 515 F: (GTGYCAGCMGCCGCGGTAA) and 806 R (GGACTACNVGGGTWTCTAAT) (Zhao et al., 2023).
After obtaining the raw sequence (Raw reads), low-quality sequences (reads) are filtered, followed by assembly and re-filtering to obtain representative sequences of operational taxonomic units (OTUs). Uparse software is used to cluster all effective sequences of all samples, and sequences are clustered into OTUs with 97% consistency, and the absolute abundance and relative information of each OTU in each sample are calculated (Li et al., 2016). Select all OTUs with an average abundance greater than 1 in the comparative groups (i.e., high-abundance union OTUs) for Venn analysis. During the construction of OTUs, Uparse will select representative sequences (the Tag sequence with the highest abundance in the OTUs). These representative sequences are then used with the Naïve Bayesian assignment algorithm of the RDP Classifier to annotate species against the Silva database (Edgar, 2013), with a confidence threshold set between 0.8 and 1.0 (Xie L. L. et al., 2024). Based on the species annotation information of OTUs, count the number of Tags sequences for each sample at the phylum level. Based on the counting statistics of OTUs, analyze the differences in soil bacterial communities using the Bray-Curtis distance coefficient. Use Faprotax to predict the ecological functions of soil bacteria.
2.5 Statistical analysis
Excel 2021 is used to sort out the data, and then Kolmogorov-Smirnov is used to test the normality of the data to determine whether it meets the normality hypothesis of the subsequent statistical test. SPSS 24 (IBM Corp., Armonk, NY, USA) was used to analyze the diversity index (Simpson, Shannon, Pielou and good_coverage index) by one-way analysis of variance (one-way ANOVA; α=0.05) to determine the significance between different altitude gradients, and the least significant difference method (LSD) is used for multiple comparisons. All data in the table are average standard error. Principal Coordinate Analysis (PCoA) is a dimensionality reduction analysis of microbial community based on Unweight unifrac distance, which evaluates the degree of explanation of bacterial community structure on each coordinate axis by percentage. The relationship between bacterial colony structure and soil physical and chemical properties was obtained by redundancy analysis,RDA) of Canoco 5.0 software package, and was plotted by Origin 2021. Linear discriminant analysis(LDA) effect size (Lefse) method is used to determine the microbial groups with significant differences among soil samples, and the unique microbial groups are detected with LDA threshold of 2.5 and significance less than 0.05, and all samples are compared (Zhang and Laanbroek, 2020). One-way similarity analysis (ANOSIM) based on unweighted UniFrac distance was used to test whether the difference of β diversity between treatments was significantly higher than that within groups. The relationship between α and β-diversity and bacterial community was discussed by Mantel test. The correlation network diagram between microbial community at the door level and soil physical and chemical properties was constructed by Spearman correlation coefficient using vegan software package. Then, a molecular ecological network including all samples was created, and the abundance of OTU was normalized by using the Hmisc package of R, and the Spearman correlation among parameters was estimated. Then, the molecular network was constructed by using the RMT model. Subsequently, Gephi 9.2 was used to visualize the network and calculate the parameters, including the number of nodes and edges, the average degree (avgK), the average path length (GD) and the average clustering coefficient (avgCC), to determine the interaction of bacteria between different altitude gradients (Barberán et al., 2012). According to OTU classification information and FAPROTAX database, the function of soil bacterial community was predicted.
3 Results
3.1 Soil properties, microbial biomass and enzyme activities at different altitude gradients
The results show that altitude has significant effects on soil total carbon, total nitrogen, soil bulk density, microbial biomass carbon, microbial biomass nitrogen and microbial biomass phosphorus (Table 1). Soil total carbon and total nitrogen first increased and then decreased with the elevation, all of which were significantly higher than 3400, 4000 and 4300 m. With the elevation, the soil bulk density increased first, then decreased, and then increased, and the soil bulk density content reached the maximum at 3400 m. The biomass carbon, nitrogen, and phosphorus of microorganisms first increased, then decreased, and then increased again along the altitude gradient. Specifically, within the range of 3400-4300m, the content of all three reached their peak at 3700m (MBP: 21.50 ± 4.91mg/kg; MBC: 2334.49 ± 268.49 mg/kg; MBN: 162.16 ± 29.19 mg/kg), then significantly decreased at 4000m, and showed an increasing trend again at 4300m. Compared to the microbial biomass at 4000m, phosphorus, carbon, and nitrogen increased by 10.27 mg/kg, 478.05 mg/kg, and 67.49 mg/kg, respectively, at 4300m.
As shown in Figure 2, the single factor analysis of variance shows that altitude significantly affects the enzyme activities related to soil carbon, nitrogen and phosphorus cycling (Figure 2). With the increase of altitude, the activities of enzymes related to soil phosphorus cycle first increased and then decreased, reaching the highest at 3700 m (Figure 2A). Altitude significantly affected the enzyme activities related to soil carbon cycle (β-1,4-glucosidase, α-1,4-glucosidase and β-D-cellobiohydrolase). Among them, β-1,4-glucosidase and α-1,4-glucosidase decreased with the elevation (Figures 2B, C), while β-D-cellobiohydrolase increased first and then decreased with the elevation (Figure 2D). L-leucine aminopeptidase related to soil nitrogen cycle decreased with the elevation, reaching the lowest at 4300 m (Figure 2E), while β-1,4-N-acetylglucosaminidase first increased and then decreased with the elevation, reaching the highest at 4000 m (Figure 2F).

Figure 2. Changes of soil enzyme activities at different altitudes. alkaline phosphatase (A);BG, β-1,4-glucosidase (B); AG,α-1,4-glucosidase (C);CBH, β-D-cellobiosidase (D); LAP, L-leucine aminopeptidase (E); NAG, β-1,4-N-acetylglucosaminidase (F). Different lowercase letters indicate significant differences at the 0.05 confidence level.
3.2 α diversity of bacterial communities under different altitude gradients
As shown in Figure 3, there are significant differences in bacterial α diversity among different altitude gradients (Figure 3). With the increase of altitude, the Shannon index first increased and then decreased, and the Shannon index reached its peak at 4300 m (9.52) (Figure 3A). However, the difference of good_coverage index between different altitudes is not significant (Figure 3D). Simpson index and Pielou index first decreased and then increased with the elevation, and both reached the peak at 4300 m, which were 0.99 and 0.76 respectively (Figures 3B, C).

Figure 3. Diversity and richness index of bacteria in different altitude gradients. Shannon index (A), Simpson index (B), Pielou index (C), and good coverage index (D). Different lowercase letters indicate significant differences at the 0.05 confidence level.
3.3 Composition and structure of bacterial community under different altitude gradients
At the phylum level (Figure 4A), the dominant groups of soil bacteria in different altitude gradients are mainly Acidobacteriota, Verrucomicrobiota, Proteobacteria and Planctomycetota, with abundance of 25.69%, 19.24%, 12.36% and 6.85% respectively. The relative abundance of bacteria at different phylum levels varied significantly with altitude gradient (P<0.05). With the elevation, Acidobacteriota first increased and then decreased, reaching the maximum relative abundance (30.27%) at 4000 m, while Verrucomicrobiota showed a downward trend with the elevation gradient, and the relative abundance was the lowest (16.77%) at 4000 m. Proteobacteria first increased and then decreased with the elevation, with the highest relative abundance (13.57%) at 4000 m, while Planctomycetota showed an upward trend with the elevation, with the highest relative abundance (7.49%) at 4300m (Figure 4A). Venn shows that there are 35 phylum between different altitude gradients, and the total number of bacteria OTUs under different altitude gradients is 8144 (Figure 4B).

Figure 4. Composition of soil bacterial community at the door level (A), Venn analysis of bacteria OTU at different altitude gradients (B).
A total of 2281 operational taxonomic units (OTUs) were identified. The number of OTUs at different altitudes was 4820 (3400 m), 4345 (3700 m), 4612 (4000 m), and 4782 (4300 m), respectively. Correspondingly, the number of unique bacterial OTUs at these altitudes was 1147 (3400 m), 665 (3700 m), 677 (4000 m), and 871 (4300 m), respectively.
UPGMA method and principal coordinate analysis (PCoA) were used to cluster soil bacterial communities (Figure 5A). UPGMA classification can study the similarity between different altitude gradients, and UPGMA identified four groups, namely 3700, 4000, 4300 and 3400 m, in which the community structure at 3700 and 4000 m is similar, while the community structure at 3400 m is different from other altitudes (Figure 5B). In addition, PcoA, ANOSIM and Adonis analysis showed that there were significant differences in bacterial communities among different altitude gradients (Adonis P=0.001; Anosim P=0.001), the first and second principal components explain 29.89% and 16.19% of the variation, respectively.

Figure 5. Arithmetic average unweighted paired grouping (UPGMA) clustering (A), principal coordinate analysis of bacterial community diversity at different altitude gradients (B).
3.4 Indicator microbial groups of bacterial communities under different altitude gradients
LefSe analysis showed that altitude significantly changed the bacterial community. Specifically, at the level of phylum, class and genus, some bacterial groups are significantly enriched at different altitude gradients (Figure 6). The number of indicator bacteria at each altitude is as follows: 3400 (22), 3700 (9), 4000 (5) and 4300 (12) for 22 kinds of indicator bacteria determined by LDA>2.0. Most bacterial groups are mainly concentrated at an altitude of 3400 m. Specifically, at the phylum level, Crenarchaeota, Patescibacteria, Desulfobacterota and Gemmatimonadota are indicator bacteria at an altitude of 3400 m. Verrucomicrobiota and Methylomirabilota are indicator bacteria at an altitude of 3700 m. Acidobacteriota is an indicator bacterium at an altitude of 4000 m. Chloroflexi is an indicator bacterium at an altitude of 4300 m.

Figure 6. The phylogenetic tree diagram shows significant differences between bacterial enrichment groups, the LDA score diagram shows biomarkers at different altitudes, the groups with significant differences in abundance at different altitudes are represented by colored dots, and the phylogenetic groups from phylum to genus are represented by phylogenetic tree circles, only the LDA > 2.
3.5 Changes of bacterial functional groups under different altitude gradients
Based on FAPROTAX database, the soil bacterial community was functionally labeled, and 56 functional groups were obtained. The results show that the altitude gradient significantly affects the functional groups related to carbon and nitrogen cycle (Figure 7). With the increase of altitude, functional groups related to nitrogen cycle, such as fumarate_respire, sulfate_respiratic and xylanolysis, showed a downward trend (Figures 7C, E, F), while celluloiysis first decreased and then increased (Figure 7A), and chitinolysis first increased and then decreased (Figure 7B). However, with the increase of altitude, mthylotrophy first decreases, then increases and then decreases, and the lowest is at 3700m (Figure 7D). Functional groups related to nitrogen cycle, tend to decrease with the increase of altitude (Figures 7G, J, K). However, the aerobic_chemoheterotrophy and Chemoheterotopy showed an upward trend with the elevation ladder (Figures 7H, I). Ureolysis first increases and then decreases with altitude, reaching the highest at 3700 m (Figure 7L).

Figure 7. Changes of functional groups of soil bacteria at different altitudes. Different letters indicate significant differences between treatments (P< 0.05). (A) cellulolysis, (B) chitinolysis, (C) fumarate_respire, (D) methylotrophy, (E) sulfate_respiratic, (F) xylanolysis, (G) aerobic_ammonia_oxidation, (H) aerobic_chemoheterotrophy, (I) chemoheterotrophy, (J) nitrate reduction, (K) nitrification, (L) ureolysis.
3.6 The relationship between bacteria, soil, microbial biomass and enzyme activity
Spearman correlation analysis showed that there was a significant relationship between dominant bacteria and soil factors (Figure 8A). NAG, ALP, TC and TN were negatively correlated with Chloroflexi. TC and TN were negatively correlated with Actinobacteriota and Gemmatimonadota. ALP,TC,TN and Methylomirabilota is a very significant positive correlation. There is a significant positive correlation between AG and Crenarchaeota (Figure 8A). Mantel test was used to evaluate the α diversity and β diversity of bacteria related to soil physical and chemical properties, microbial biomass and enzyme activity (Figure 8B). Mantel test revealed that the observed α diversity of bacteria was positively correlated with TC, TN, SM, BD and ALP. β diversity was positively correlated with BD, BG, NAG, LAP, CBH and AG (Figure 8B).

Figure 8. Thermal map of Spearman correlation coefficient between the top 10 dominant bacteria and soil factors (A). Mental test of the relationship between α and β diversity and environmental factors (B). *P<0.05,**P<0.01,***P<0.001.
RDA was used to analyze the relationship between soil bacterial communities and functional groups and soil nutrients and enzyme activities (Figure 9). In the first two axes, 32.79% and 15.61% were explained respectively (Figure 9A). The results showed that the main driving factors of soil bacterial community were TC, CBH, AG, BD and BG. The arrow of TC is the longest, which indicates that TC has the greatest influence on the overall distribution of bacterial communities in different altitude gradients, followed by CBH, AG, BD and BG. As shown in Figure 8A, the arrow of TC is the longest, which indicates that TC has the greatest influence on the overall distribution of bacterial communities in different altitude gradients, followed by CBH, AG, BD and BG. The results of Monte Carlo experiments show that the important values of soil physical and chemical properties decrease in the order of TC>CBH>Ag>BD>BG>pH>SM>LAP>MBC>MBP>NAG>TP>TN>MBN>ALP (Table 2). In addition, TC, CBH and AG are significantly different at the level of 0.01, while BD and BG are significantly different at the level of 0.05, and the contribution rates of these five factors are 22.3%, 11.3%, 11.4%, 6.6% and 5.8% respectively (Table 2). These results show that TC, CBH, AG, BD and BG are the main driving factors affecting the bacterial community structure at different altitude gradients. RDA was used to analyze the relationship between soil bacterial functional groups and soil nutrients and enzyme activities (Figure 8). In the first two axes, 22.33% and 22.16% were explained respectively (Figure 9B). The main factors affecting bacterial functional groups are AG, TN, BG, MBC and MBN. The difference of AG is significant at the level of 0.01, while that of TN, BG, MBC and MBN is significant at the level of 0.05. The contribution rates of these five factors are 18.60%, 8.8%, 7.0%, 7.0% and 6.6% respectively (Table 3).

Figure 9. Ranking chart of redundant analysis (RDA) results to determine the relationship between bacterial communities and soil factors (A) and the relationship between bacterial functional groups and soil factors (B).
3.7 Co-occurrence network structure of different altitude gradients
Cluster analysis divides different altitude gradients into three groups (Figure 3), and the co-occurrence network reveals their significant differences (Figure 6). The complex relationship between nodes is described by calculating the topological properties of the network. Generally speaking, the number of 4300 m network nodes (191) is lower than 3400 (195), 3700 (198) and 4000 m (199) (Figure 6A), while the number of links is higher than 3700 (920) and 4000 m (1132) (Figure 6B). In addition, the altitude of 4300 m has a higher average degree and modularity (Figures 6G, H), indicating the cluster topology and modular structure. The network complexity at 4300 m is higher, indicating that the relationship between species is more complicated and the network stability is higher. On the contrary, the 3400 m network is more susceptible to environmental interference, with a modularity index of 0.53 (Figure 6H).
4 Discussion
4.1 Changes of soil physical and chemical properties, microbial biomass and extracellular enzymes in different altitude gradients
Altitude is one of the important topographical factors. Due to different altitudes, changes in climate characteristics and soil types all lead to differences in soil nutrients (Gilgen and Buchmann, 2009).Therefore, soil is affected by various ecological factors, and complex biological and abiotic factors closely interact to form a soil ecosystem (Denef et al., 2009). In this study, total carbon and total nitrogen first increased and then decreased with the elevation, which is consistent with the research conclusion of (Zhou, 2019). Altitude, as a natural geographical variation, is one of the important factors affecting the distribution of soil carbon and nitrogen (Xie L. L. et al., 2024). Due to altitude, continuous changes in heat, temperature, precipitation, and sunlight occur within the region (Bin et al., 2022). This is consistent with the findings of Gao et al. (2021). The total carbon (TC) and total nitrogen (TN) in the study area appeared at 3700 meters rather than 4300 meters. The main reason is that as the altitude rises to 3700m, human interference (grazing activities) gradually decreases compared to lower altitudes (3400m), soil moisture content gradually increases, and the content of plant litter and organic matter increases, which can provide more nutrients for the growth of soil microorganisms, promote the reproduction and activity of soil microorganisms, accelerate the utilization and mineralization rate of organic matter by soil microorganisms, thereby improving the physical and chemical properties of the soil (Pingree and Deluca, 2018). However, as the altitude continues to rise, wind erosion gradually worsens, the soil layer becomes thinner, the hydrothermal conditions deteriorate, and the shelter conditions decrease, causing significant changes in soil nutrients (Deng et al., 2019). Therefore, the total carbon and total nitrogen content of soil at different altitude gradients shows significant heterogeneity. Soil total phosphorus increased with the elevation (Table 1). A main driving factor seems to be that the degree of soil weathering decreases with the decrease of altitude, because soil phosphorus mainly comes from weathered soil (Helfenstein et al., 2018). Soil microbial biomass first increased and then decreased and increased with the elevation, and the microbial biomass was the highest at 3700 m above sea level (Table 1). The reason may be that the surface soil temperature at 3700 m above sea level is high, which is beneficial to plant growth and there are more exudates from vegetation roots, Higher soil temperature promoted microbial proliferation and increased microbial biomass (Uselman et al., 2000). It is worth noting that at 4300 m, microbial biomass phosphorus, microbial biomass carbon and microbial biomass nitrogen are higher than 4000 m, indicating that microbial biomass is affected by many factors (Ren et al., 2021). With the elevation, its climatic conditions will change (Xu et al., 2013). At the same time, soil microorganisms have a certain physiological adaptability to climate conditions. Although altitude increases will cause temperature drops, microorganisms can maintain balance by reducing metabolic activity (Li et al., 2014). Secondly, the low temperatures in high-altitude areas may increase physiological stress by altering osmotic pressure and reducing the ability to obtain effective nutrients Sun et al., 2025).
In this low-temperature environment, the decomposition of soil organic matter is slow, while the transformation of microbial biomass remains high (Wang B. R. et al., 2021; Wang R. et al., 2021). Soil extracellular enzymes can reveal the cycling law of nutrients such as C, N and P in soil, and reflect the intensity and direction of various biochemical processes (Kotroczo et al., 2014). At present, the most widely studied soil extracellular enzymes are BG,CBH,NAG,ALP and LAP, because they can be used as indicators of C demand, N demand and P demand respectively (Schimel and Weintraub, 2003). Specifically, BG hydrolyzes cellobiose into glucose during the carbon cycle, CBH hydrolyzes cellulose to produce sucrose during the carbon cycle, NAG hydrolyzes chitin during the carbon-nitrogen cycle, LAP hydrolyzes protein peptide during the nitrogen cycle, releasing amino acids from the N-terminal, and ALP hydrolyzes phosphate from phosphate sugars and phospholipids during the phosphorus cycle (Sinsabaugh et al., 2009). Previous studies have shown that altitude difference causes the change of micro-habitat, which has a gradient effect on soil physical and chemical properties, resulting in different soil enzyme activities (He et al., 2009). Studies such as Zhen et al. (2019) showed that soil enzyme activities showed significant differences between different altitudes. Our study found that the activities of ALP, BG, AG and LAP decreased significantly with the elevation, and our results were inconsistent with those recently reported, which reported that the soil enzyme activities in southern Patagonia did not respond significantly to the elevation gradient (Truong et al., 2019). This difference may be attributed to the different sampling points used in the two studies. In our study, the sampling points covered 3400-4300 m, while in the other study, the sampling points only covered 130-640 m. Different vegetation types in different regions respond differently to soil microbial communities and enzyme activities. In addition, In our research, we found that ALP is related to soil properties, especially to total carbon and total nitrogen content (S2). Our results are consistent with previous research reports that soil enzyme activities are strongly influenced by soil properties, especially total carbon (Truong et al., 2019). These results further confirmed that the extracellular enzyme activities associated with N and P cycling increased significantly with altitude, while those associated with C cycle decreased significantly with altitude, which supported our first hypothesis. To sum up, there is no unified conclusion on the influence of altitude gradient on soil enzyme activity, which is mainly due to the similarities and differences of regional and soil substrate conditions in the response of soil enzyme activity to microclimate change. When analyzing soil microbial activity, soil extracellular enzyme activity is also the most important analysis item. Therefore, the study on the dynamic relationship between soil extracellular enzyme activity and soil characteristics is helpful to understand the influence of microorganisms on ecosystem processes such as litter decomposition, soil carbon cycle and nutrient cycle.
4.2 Changes of bacterial community composition and structure at different altitude gradients
Soil microorganisms are very sensitive to environmental changes, and macro variables such as altitude gradient cause changes in micro factors such as soil properties, temperature and moisture, thus affecting soil microbial community structure to varying degrees (Shen et al., 2013). In this study, UPGMA and principal coordinate (PCA) analysis found that bacterial communities at different altitudes can be well distinguished, and the characteristics of soil bacterial communities changed significantly with the change of altitude (Figure 5). The bacterial communities at 3400 m above sea level are obviously different from those at other altitudes. In addition, the bacterial communities at 3700 m, 4000 m and 4300 m above sea level are similar in composition and structure (Wang B. R. et al., 2021; Wang R. et al., 2021). These results show that different altitude gradients can lead to soil bacteria with different evolutionary characteristics. Specifically, in this study, Acidobacteria is the dominant bacteria, followed by Verrucomicrobiota and Proteobacteria. This is similar to Cui et al. (2019) finding that the dominant bacteria phylum is Actinomycetes, Bacteroides and Blastomycetes in the alpine ecosystem at an altitude of 2800-3500 m on the Qinghai-Tibet Plateau. Just as we found in our study that the relative abundance of various groups of microbial communities changed significantly with the elevation (Figure 4A), the co-occurrence model also supported our findings (Figure 10). Acidobacteria, Verrucomicrobiota and Proteobacteria are always the dominant populations of soil bacterial community at different altitudes, which have made great contributions to the variation of community composition, but the highest relative abundance of the three species appears at different altitudes. This result is similar to the bacterial diversity of Xie L. L. et al. (2024) in the desert habitat with an average elevation of 2800 m on the Qinghai-Tibet Plateau, in which Proteobacteria, Actinomycetes, Bacteroides and Chlorocurvata are the dominant groups, indicating that there are similar bacterial compositions in soils in different habitats, and at the same time, their respective dominant groups are similar.This is mainly related to the wide niche of Acidobacteria, Verrucomicrobiota and Proteobacteria, and their strong adaptability to different environments (Gao et al., 2019). In addition, the difference in adaptability of flora to microenvironment is the reason why the maximum relative abundance of the three major fungi appears at different altitudes. Studies have shown that Acidobacteria can improve soil nutrient content by degrading complex lignin and cellulose (Pankratov et al., 2011). Soil microbial communities respond to abiotic and biotic conditions. In mountain plant communities with vertical gradients, altitude can create unique microenvironments, which have unique microbial responses (Hernández-Cáceres et al., 2022). This shows that the distribution of soil bacteria is not only affected by altitude, but also by soil heterogeneity, which shows different research results. On the one hand, soil bacterial species have high spatial heterogeneity and are very sensitive to climate, soil and vegetation conditions, which leads to significant differences in dominant bacterial species in different research areas (Cao et al., 2022). This also shows the relative abundance of soil bacterial species in the high altitude area of Qinghai-Tibet Plateau, and the study of soil microorganisms still has great exploration space in this area. On the other hand, due to the existence of altitude gradient, the natural, soil and biological environment suddenly changed in a short geographical distance, and these gradients created a variety of bacterial ecological habitats, which led to great differences in soil bacteria species in the Pfruticosa shrub meadow under different altitude gradients on the Qinghai-Tibet Plateau (Hernández-Cáceres et al., 2022).

Figure 10. Based on correlation analysis, 3400 (A), 3700 (B), 4000 (C) and 4300 m (D) bacterial community coexistence networks. (E–H) represent the main characteristics of the co-occurrence network. The color and size of nodes represent the relative abundance of OTU groups and specific OTUs at the gate level.
The α diversity of soil bacterial community directly reflects the change of soil bacterial community quantity and richness (Ding et al., 2023). Many studies show that α diversity decreases with the elevation (Yao et al., 2017). On the contrary, this study found that the diversity of soil bacteria in the study area was different at different altitude gradients (Figure 3), and the differences between different altitude gradients were obvious. At an altitude of 3400 m, there is the largest number of species unique to soil bacteria (Figure 4B). Studies have shown that specific soil microbial communities exist at specific altitudes (Liu et al., 2013). Shannon index and Pielou index of soil bacteria decreased first and then increased along the elevation gradient, and the maximum and minimum values appeared at the altitude of 3400 m and 3700 m respectively (Figures 3A,C). The reason may be related to the temperature and precipitation on the altitude gradient in Qinghai-Tibet Plateau (Gong et al., 2022). Therefore, the number and species of soil microorganisms decrease at the 3700 m. At the 3400 m, the bacterial diversity may increase slightly due to the increase of soil water content. In addition, the changes of soil Simpson index and good_coverage index along the elevation gradient are not obvious in this study (Figures 3B, D), which is consistent with the research results of Li et al. (2022) indicating that bacterial diversity may be difficult to form a unified pattern along the elevation gradient. These results further confirm our second hypothesis.
LEfSe analysis provides a new perspective for studying the response of soil bacteria to soil niche change.In the current research, most bacterial groups are obviously enriched at an altitude of 3400 m, indicating that these bacteria have a stable niche at an altitude of 3400 m. However, at an altitude of 3700 m, only Verrucomicrobiae is abundant. It is generally believed that in a long-term stable community (that is, a community with little interference), Some highly competitive species will reproduce, making a few species dominant in the community (Doležal et al., 2013).
The change of soil microbial community function is similar to the change of community structure, and different altitude gradients form different microbial functional groups. Different functional groups are often filtered into different environments, and these environments are the characteristics of biochemical cycles in these areas (Crowther et al., 2019). This study showed that at different altitudes, the bacterial functional group related to Nitrogen cycle changed obviously. Aerobics_ammonia_oxidation, nitrate_reduction and nitrification decrease with altitude, while aerobic_chemoheteotrophy, chemotetropy and ureolysis increase obviously with altitude. Aerobic_ammonia_oxidation and nitrification (Figures 7G, K) have the most obvious changes, which are much higher than other functions. This phenomenon and characteristics prove the protection strategy of limiting element nitrogen in different altitude gradients (Xi et al., 2017). It shows that the effect of nitrate_reduction is more intense at the altitude of 3400 m, and these different nitrogen-reducing bacteria groups may be related to the mechanism of soil nitrogen stability. For example, certain microbial communities (such as Bacillus, Pseudomonas and Bacillus subtilis) limit nitrogen loss by reducing nitrate to a more stable form (Kuypers et al., 2018). It is attributed to the fact that shrubs have nitrogen fixation characteristics (Hao et al., 2024). Furthermore, because Proteobacteria has a variety of species involved in nitrogen fixation, it can form a symbiotic relationship and fix nitrogen as a nitrogen compound that plants can use (Creamer et al., 2016). Similar to the change of microbial community structure, microbial functional groups at different altitudes are also different (Figure 7). The Carbon cycle of bacteria related to carbon cycle, such as cellulolysis, fumarate_respire, methylotrophy, sulfate_respiratic and xylanolysis, varies significantly with altitude, and their contents are all the highest at 3400 m. The relative abundance of methylotrophy is the largest, which may be due to the high vegetation coverage at 3400 m above sea level, and the great change of litter input will promote aerobic environment, thus accelerating the oxidation of CH4 (Qu et al., 2020).
4.3 Correlation between bacterial community and soil properties
The changes of soil microorganisms and functional groups in different altitude gradients of P. fruticosa shrub meadow are closely related to soil environmental factors (Figure 9). Spearman correlation shows that TC and TN are positively correlated with Methylomirabilota and Verrucomicrobilota (Figure 9A). TC is the most important factor affecting bacterial community and diversity (Figure 8B; Table 1). This is consistent with the research results of Jiang et al. (2021), which showed that microbial communities were significantly correlated with soil nutrients. The research shows that soil carbon accumulation affects the structure of soil microbial community under vegetation by affecting the diversity, productivity and stability of plant community, which explains why the change of TC will significantly affect the structure of bacterial community (Chen et al., 2021). This supports our second hypothesis. Secondly, the soil bulk density is the second environmental factor that affects the bacterial community, which may be because the rhizosphere soil of P. fruticosa shrub will produce the “ Fertility islands effect”, which makes the soil nutrient content enriched around its roots (Michalet et al., 2016). However, most microorganisms, animals and plants in the soil will draw enough nutrients for their own growth, so they concentrate on grabbing nutrients near the rhizosphere of the P. fruticosa shrub, which makes the soil nutrient content in short supply. However, due to the increase in altitude, the number of plants and animals decreases, and the exposed areas are relatively extensive, which makes the number of animals and plants in this area decrease obviously, and the number of soil microorganisms also decreases. Therefore, the dead bodies of litter and dead animals and plants cannot be effectively utilized, and the soil organic matter is low. It further aggravates the degree of soil hardening, resulting in high soil bulk density, poor water stability and serious soil water loss. The local micro-habitat situation is not conducive to the reproduction and metabolism of soil microorganisms, so soil bulk density has gradually become the main environmental limiting factor to limit bacteria (Yang et al., 2013). Soil extracellular enzymes are closely related to a series of biological and abiotic factors (Ivashchenko et al., 2021). In our study, TC and TN were positively correlated with ALP (S2), and CBH and AG were positively correlated with Chloroflexi, Planctomycetota and Crenarchaeota (Figure 8A). This confirms the view that soil total nutrients play a key role in the altitude model of regulating extracellular enzyme activity (Cao et al., 2021).
The main environmental limiting factor in bacterial functional groups is total nitrogen (Figures 9B, Table 2). This may be because different microbial functional groups perform different functions at different stages of the material cycle, and they work together to promote the rapid and orderly flow of matter and energy in the ecosystem and ensure the normal succession of grassland ecosystem (Gordon et al., 2008). The change of environmental factors has a direct impact on microbial functional groups. Nitrogen-fixing bacteria transform N2 into NH4+ for plants to synthesize organic nitrogen and supplement nitrogen in grassland soil; Ammoniated bacteria nitrated NH4+ into NO3-; Nitrifying bacteria convert NO3- into N2 to complete the nitrogen cycle, which is related to the enzyme activity, so soil enzyme activity is also the main factor affecting microbial functional groups (Wang B. R. et al., 2021; Wang R. et al., 2021). Nevertheless, due to the large reduction of plants and microorganisms and the intensification of soil erosion caused by bare plots, the only soil nutrients are still not enough to maintain the normal growth of plants and microorganisms, and the competition among species is intensified (Ahmed et al., 2019). Soil nitrogen content is an indispensable environmental factor to maintain the normal flow of microbial functional groups, so soil total nitrogen is the main environmental factor to limit microbial functional groups (Yang et al., 2008). Secondly, the influencing factors are microbial biomass carbon, nitrogen and extracellular enzymes (AG and BG) related to carbon cycle. Microbial biomass carbon is an important indicator reflecting the size of microbial community, and extracellular enzymes are mainly produced by microorganisms, and the size of microbial biomass determines the potential of enzyme production (De Vries et al., 2018). Soil microbial biomass carbon is an important part of active organic carbon pool. Because of its fast turnover rate and easy decomposition, it is also a great supply source of soil effective nutrient pool (Wang B. R. et al., 2021; Wang R. et al., 2021 ). These fast-moving carbon sources can continuously provide energy for soil microorganisms and keep them relatively high in activity, so as to drive the biogeochemical cycle of nutrients such as carbon and nitrogen, and then affect microbial functional groups (Zhou et al., 2022).
4.4 Co-occurrence patterns of bacteria under different altitude gradients
Species form a complex network system by secreting metabolites and interacting with other microorganisms (Peura et al., 2015). The co-occurrence network among microorganisms reflects the function and stability of biological communities (Ma et al., 2020). It provides a preliminary means to explore the organization and dynamics of microbial interaction and niche (Banerjee et al., 2019). It is widely used to estimate co-occurrence patterns in complex microbial relationships (Deng et al., 2012). It has been pointed out that the positive correlation connection in the co-occurrence network represents the mutual synergy between microorganisms, while the negative correlation connection represents the potential antagonism between microorganisms (Blanchet et al., 2020). In this study, the connection of soil bacterial co-occurrence network is mainly positive correlation (bacterial network: 97.26%-99.3% positive correlation link), and the complexity of soil microbial network is the highest at 4300 m, followed by 3400 and 4000 m, and the complexity of 3700 m is the lowest (Figure 10). These differences are reflected in factors such as average degree and modularity index. Our research is consistent with Duan et al. (2021) research on microbial co-occurrence network in Qinghai-Tibet Plateau. This may indicate that when microbial communities are stressed by cold and harsh environment (such as strong solar radiation, low soil temperature and low soil oxygen content), soil microbial communities can maintain their ability to resist external interference through mutual cooperation (Hernández-Cáceres et al., 2022). The complexity of the co-occurrence network tends to increase with the altitude, which indicates that the scale of the microbial co-occurrence network is large in high altitude areas, and the interaction between microorganisms has been enhanced. Similar conclusions have been drawn in the previous research on the bacterial co-occurrence network in Gongga Mountain (Zhu et al., 2020). The possible reason is that the microbial community has been stressed by the cold and harsh environment (Hernández-Cáceres et al., 2022). This shows that the soil microbial community in high altitude areas is under greater environmental stress than that in low altitude areas, and the co-occurrence network structure of soil microorganisms in high altitude areas will become more complicated. Therefore, more diverse and richer microbial groups create greater possibilities, that is, there may be more potential taxa in the environment to participate in potential interactions, thus leading to more complex association networks.
5 Conclusion
This study explored the distribution patterns and driving factors of soil bacterial community structure, diversity, and enzyme activity along an altitude gradient in P. fruticosa shrub meadows. It was found that the soil bacteria in P. fruticosa shrubs have different adaptabilities, which are closely related to soil nutrients and enzyme activity. This study provides a new perspective for further understanding the distribution pattern and driving forces of soil bacterial communities in P. fruticosa shrubs along the altitude gradient, enhancing our understanding of microbial ecology in alpine shrub environments. We observed significant differences in soil bacterial α-diversity and community composition at different altitudinal gradients. The altitudinal gradient influences soil bacterial communities by correlating with soil properties, microbial biomass, and enzyme activity. Soil properties (TC and BD) and enzyme activities (CBH, AG, and BG) are the main factors affecting changes in soil bacterial communities, while AG, TN, BG, MBC, and MBN have significant impacts on bacterial functional groups.
We also found that TC, TN, SM, MBC, MBN, and MBP were highest at 3700 m altitude, while ALP was the most active. Moreover, the dominant bacterial phylum shifted from Acidobacteria (at 3400 m) to Verrucomicrobia (at 3700 m). This shows that this area has special ecological adaptability and should be protected first to avoid interference. It is recommended that protective land use measures (such as restricting overgrazing and development) be implemented in this area, and a monitoring system for P. fruticosa shrub vegetation be established to keep track of its growth and water conservation status in real time. Based on the monitoring data, protection and management measures should be adjusted in a timely manner to maintain its unique microbial community and ecological functions, thereby ensuring the long-term stability and sustainable development of the ecosystem.
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 authors.
Author contributions
LX: Conceptualization, Investigation, Methodology, Visualization, Writing – review & editing. YM (2nd author): Conceptualization, Funding acquisition, Methodology, Resources, Supervision, Writing – review & editing. YW: Investigation, Writing – original draft, Writing – review & editing. YM (4th author): Validation, Writing – review & editing. YL: Formal Analysis, Software, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work has been supported by the Chief Scientist Program of Qinghai Province (2024-SF-101).
Acknowledgments
We sincerely thank this fund for its support in carrying out the experiment. At the same time, we thank Zhou Xuanbo, Wang Xinyou and Zhang Lijun for their help in field work and laboratory work.
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 author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
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/fpls.2025.1539945/full#supplementary-material
References
Ahmed, I. U., Mengistie, H. K., Godbold, D. L., and Sandén, H. (2019). Soil moisture integrates the influence of land-use and season on soil microbial community composition in the Ethiopian highlands. Appl. Soil Ecology. 135, 85–90. doi: 10.1016/j.apsoil.2018.11.010
Arunrat, N., Sansupa, C., Sereenonchai, S., Hatano, R., and Lal, R. (2024). Fire-induced changes in soil properties and bacterial communities in rotational shifting cultivation fields in Northern Thailand. Biology 13, 383. doi: 10.3390/biology13060383
Bai, Y. X., She, W. W., Miao, L., Qin, S. G., and Zhang, Y. Q. (2020). Soil microbial interactions modulate the effect of Artemisia ordosica on herbaceous species in a desert ecosystem, northern China. Soil Biol. Biochem. 150, 108013. doi: 10.1016/j.soilbio.2020.108013
Banerjee, S., Walder, F., Büchi, L., Meyer, M., Held, A. Y., Gattinger, A., et al. (2019). Agricultural intensification reduces microbial network complexity and the abundance of keystone taxa in roots. ISME J. 13, 1722–1736. doi: 10.1038/s41396-019-0383-2
Barberán, A., Bates, S. T., Casamayor, E. O., and Fierer, N. (2012). Using network analysis to explore co-occurrence patterns in soil microbial communities. ISME J. 6, 343–351. doi: 10.1038/ismej.2011.119
Bayranvand, M., Akbarinia, M., Jouzani, G. S., Gharechahi, J., Kooch, Y., and Petr, B. (2021). Composition of soil bacterial and fungal communities in relation to vegetation composition and soil characteristics along an altitudinal gradient. FEMS Microbiol. ecology. 97, fiaa201. doi: 10.1093/femsec/fiaa201
Bin, H., Li, Q., Zhang, P., Li, W. J., Xue, X. H., Zou, S., et al. (2022). Effects of Elevation on ecological stoichiometry of plant leaves, litter, and soils in pseudotsuga sinensis forest in the Karst Mountain region, Southwest China. J. Soil Sci. Plant Nutr. 22, 3582–3597. doi: 10.1007/s42729-022-00911-y
Blanchet, F. G., Cazelles, K., and Gravel, D. (2020). Co-occurrence is not evidence of ecological interactions. Ecol. letters. 23, 1050–1063. doi: 10.1111/ele.13525
Cao, L. H., Liu, H. M., Yang, H., and Lian, Y. Z. (2022). Soil microbial distribution and fungal community composition at different altitudes on sejila mountain, southeastern tibet. J. Soil Water Conserv. 36, 371–378. doi: 10.13870/j.cnki.stbcxb.2022.06.045
Cao, R., Yang, W. Q., Chang, C. H., Wang, Z., Wang, Q., Li, H., et al. (2021). Differential seasonal changes in soil enzyme activity along an altitudinal gradient in an alpine-gorge region. Appl. Soil Ecol. 166, 104078. doi: 10.1016/j.apsoil.2021.104078
Cenini, V. L., Fornara, D. A., McMullan, G., Ternan, N., Lajtha, K., and Crawley, M. J. (2015). Chronic nitrogen fertilization and carbon sequestration in grassland soils: evidence of a microbial enzyme link. Biogeochemistry. 126, 301–313. doi: 10.1007/s10533-015-0157-5
Chang, E. H., Chen, T. H., Tian, G. L., and C.Yu, C. (2016). The effect of altitudinal gradient on soil microbial community activity and structure in moso bamboo plantations. Appl. Soil Ecol. 98, 213–220. doi: 10.1016/j.apsoil.2015.10.018
Chen, L. L., Xu, H. B., Wu, S. Y., and Baoyin, T. (2021). Plant and soil properties mediate the response of soil microbial communities to moderate grazing in a semiarid grassland of northern China. J. Environ. Management. 284, 112005. doi: 10.1016/j.jenvman.2021.112005
Creamer, C. A., Filley, T. R., Boutton, T. W., and Rowe, H. I. (2016). Grassland to woodland transitions: Dynamic response of microbial community structure and carbon use patterns. J. Geophysical Research: Biogeosciences. 121, 1675–1688. doi: 10.1002/2016JG003347
Crowther, T. W., Van den Hoogen, J., Wan, J., Mayes, M. A., Keiser, A. D., Mo, L., et al. (2019). The global soil community and its influence on biogeochemistry. Science 365, eaav0550. doi: 10.1126/science.aav0550
Cui, Y. X., Bing, H. J., Fang, L. C., Wu, Y. H., Yu, J. L., Shen, G. T., et al. (2019). Diversity patterns of the rhizosphere and bulk soil microbial communities along an altitudinal gradient in an alpine ecosystem of the eastern Tibetan Plateau. Geoderma 338, 118–127. doi: 10.1016/j.geoderma.2018.11.047
Delgado-Baquerizo, M., Eldridge, D. J., Ochoa, V., Gozalo, B., Singh, B. K., and Maestre, F. T. (2017). Soil microbial communities drive the resistance of ecosystem multifunctionality to global change in drylands across the globe. Ecol. Lett. 20, 1295–1305. doi: 10.1111/ele.12826
Denef, K., Roobroeck, D., Wadu, M., Mihiri, C. W., Lootens, P., and Boeckx, P. (2009). Microbial community composition and rhizodeposit-carbon assimilation in differently managed temperate grassland soils. Soil Biol. Biochem. 41, 144–153. doi: 10.1016/j.soilbio.2008.10.008
Deng, Y., Jiang, Y. H., Yang, Y. F., He, Z. L., Luo, F., and Zhou, J. Z. (2012). Molecular ecological network analyses. BMC Bioinf. 13, 1–20. doi: 10.1186/1471-2105-13-113
Deng, Q. Y., Zhang, X. L., Niu, J. J., and Qin, H. (2019). Species diversity of plant communities along an altitude gradient in yinmachi mountain, northwestern shanxi, China. Ecol. Environ. Sci. 28, 865–872. doi: 10.16258/j.cnki.1674-5906.2019.05.002
De Vries, F. T., Griffiths, R. I., Bailey, M., Craig, H., Girlanda, M., Gweon, H. S., et al. (2018). Soil bacterial networks are less stable under drought than fungal networks. Nat. Commun. 9, 3033. doi: 10.1038/s41467-018-05516-7
Ding, B. J., Feng, L., Ba, S., Jiang, X. L., Liu, G. H., and Liu, W. Z. (2023). Temperature drives elevational diversity patterns of different types of organisms in Qinghai-Tibetan Plateau wetlands. Iscience 26, 6–7. doi: 10.1016/j.isci.2023.107252
Doležal, J., Yakubov, V., and Hara, T. (2013). Plant diversity changes and succession along resource availability and disturbance gradients in Kamchatka. Plant ecology. 214, 477–488. doi: 10.1007/s11258-013-0184-z
Duan, Y. L., Lian, J., Wang, L. L., Wang, X. Y., Luo, Y. Q., Wang, W. F., et al. (2021). Variation in soil microbial communities along an elevational gradient in alpine meadows of the Qilian Mountains, China. Front. Microbiol. 12. doi: 10.3389/fmicb.2021.684386
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/NMETH.2604
Fan, S. Y., Sun, H., Yang, J. Y., Qin, J. H., Shen, D. J., and Chen, Y. X. (2021). Variations in soil enzyme activities and microbial communities along an altitudinal gradient on the eastern Qinghai–Tibetan plateau. Forests. 12, 681. doi: 10.3390/f12060681
Feng, J., Zeng, X. M., Zhang, Q. G., Zhou, X. Q., Liu, Y. R., and Huang, Q. Y. (2021). Soil microbial trait-based strategies drive metabolic efficiency along an altitude gradient. ISME Commun. 1, 71. doi: 10.1038/s43705-021-00076-2
Gao, X. H., Li, M., Lu, P., Lu, G. F., and Niu, Y. F. (2019). Bacterial community in the rhizosphere soil of Betula platyphylla in the Daqing Mountains, Hohhot. Acta ECOLOGICA Sin. 39, 3586–3596. doi: 10.5846/stxb201805221123
Gao, H. N., Li, C. X., Sun, X. M., Zhang, Y., and Chen, N. L. (2021). Stoichiometry characteristics of soil at different altitudes in the Qilian Mountains. J. desert Res. 41, 219–227. doi: 10.7522/j.issn.1000-694X.2020.00125
Gilgen, A. K. and Buchmann, N. (2009). Response of temperate grasslands at different altitudes to simulated summer drought differed but scaled with annual precipitation. Biogeosciences 6, 2525–2539. doi: 10.5194/bg-6-2525-2009
Gong, X. W., Wang, X. Y., Li, Y. Q., Ma, L., Li, M. Y., and Si, H. T. (2022). Observed changes in extreme temperature and precipitation indices on the Qinghai-Tibet Plateau 1960–2016. Front. Environ. Sci. 10. doi: 10.3389/fenvs.2022.888937
Gordon, H., Haygarth, P. M., and Bardgett, R. D. (2008). Drying and rewetting effects on soil microbial community composition and nutrient leaching. Soil Biol. Biochem. 40, 302–311. doi: 10.1016/j.soilbio.2007.08.008
Hao, J. H., Zhang, X. J., and Li, J. J. (2024). Characterization of the community structure of rhizosphere soil and root-endophytic bacteria in different shrubs. Environ. Sci. 45, 6756–6765. doi: 10.13227/j.hjkx.202311212
He, R., Wang, J. S., Zheng, Z., Fang, Y. H., Xu, Z. K., Quan, W., et al. (2009). Variations of soil microbial biomass across four different plant communities along an elevation gradient in Wuyi Mountains,China. Acta ECOLOGICA Sin. 29, 5138–5144. doi: 10.3321/j.issn:1000-0933.2009.09.065
Helfenstein, J., Tamburini, F., von Sperber, C., Massey, M. S., Pistocchi, C., Chadwick, O. A., et al. (2018). Combining spectroscopic and isotopic techniques gives a dynamic view of phosphorus cycling in soil. Nat. Commun. 9, 3226. doi: 10.1038/s41467-018-05731-2
Hernández-Cáceres, D., Stokes, A., Angeles-Alvarez, G., Abadie, J., Anthelme, F., Bounous, M., et al. (2022). Vegetation creates microenvironments that influence soil microbial activity and functional diversity along an elevation gradient. Soil Biol. Biochem. 165, 108485. doi: 10.1016/j.soilbio.2021.108485
Ivashchenko, K., Sushko, S., Selezneva, A., Ananyeva, N., Zhuravleva, A., Kudeyarov, V., et al. (2021). Soil microbial activity along an altitudinal gradient: Vegetation as a main driver beyond topographic and edaphic factors. Appl. Soil Ecol. 168, 104197. doi: 10.1016/j.apsoil.2021.104197
Jiang, S., Xing, Y. J., Liu, G. C., Hu, C. Y., Wang, X. C., Yan, G. Y., et al. (2021). Changes in soil bacterial and fungal community composition and functional groups during the succession of boreal forests. Soil Biol. Biochem. 161, 108393. doi: 10.1016/j.soilbio.2021.108393
Jin, Y. H., Wang, J. S., Li, L. G., Ruan, H. H., Xu, Z. K., and Han, L. Y. (2011). Soil enzyme activities in typical vegetation zones along an altitude gradient in Wuyi Mountains. Chin. J. Ecol. 30, 1955–1961. doi: 10.13292/j.1000-4890.2011.0327
Klimek, B., Jaźwa, M., Choczyński, M., Stolarczyk, M., and Niklińska, M. (2020). The drivers of soil microbial communities structure on forest stands along the altitudinal gradient in western carpathians. Acta Oecologica 108, 103643. doi: 10.1016/j.actao.2020.103643
Kotroczo, Z., Veres, Z., Fekete, I., Krakomperger, Z., Tóth, J. A., Lajtha, K., et al. (2014). Soil enzyme activity in response to long-term organic matter manipulation. Soil Biol. Biochem. 70, 237–243. doi: 10.1016/j.soilbio.2013.12.028
Kuypers, M. M. M., Marchant, H. K., and Kartal, B. (2018). The microbial nitrogen-cycling network. Nat. Rev. Microbiol. 16, 263–276. doi: 10.1038/nrmicro.2018.9
Li, M. S., Dai, G. H., and Mu, L. Q. (2022). Composition and diversity of soil bacterial communities under identical vegetation along an elevational gradient in Changbai Mountains, China. Front. Microbiol. 13. doi: 10.3389/fmicb.2022.1065412
Li, C., Shimono, A., Shen, H. H., and Tang, Y. H. (2010). Phylogeography of Potentilla fruticosa, an alpine shrub on the Qinghai-Tibetan Plateau. J. Plant Ecology. 3, 9–15. doi: 10.1093/jpe/rtp022
Li, Y. M., Wang, S. P., and Jiang, L. L. (2016). Changes of soil microbial community under different degraded gradients of alpine meadow. Agriculture Ecosyst. Environment. 222, 213–222. doi: 10.1016/j.agee.2016.02.020
Li, P., Yang, Y. H., Han, W. X., and Fang, J. Y. (2014). Global patterns of soil microbial nitrogen and phosphorus stoichiometry in forest ecosystems. Global Ecol. Biogeography 23, 979–987. doi: 10.1111/geb.12190
Liu, B. R., Zhang, X. Z., Hu, T. H., and Li, W. J. (2013). Soil microbial diversity under typical vegetation zones along an elevation gradient in Helan Mountains. Acta ECOLOGICA Sin. 33, 7211 7220. doi: 10.5846/stxb201208061110
Ma, B., Wang, Y. L., Ye, S. D., Liu, S., Stirling, E., Gilbert, J. A., et al. (2020). Earth microbial co-occurrence network reveals interconnection pattern across microbiomes. Microbiome. 8, 1–12. doi: 10.1186/s40168-020-00857-2
Margesin, R., Minerbi, S., and Schinner, F. (2014). Long-term monitoring of soil microbiological activities in two forest sites in South Tyrol in the Italian Alps. Microbes environments. 29, 277–285. doi: 10.1264/jsme2.ME14050
Michalet, R., Schöb, C., Xiao, S., Zhao, L., Chen, T., An, L. Z., et al. (2016). Beneficiary feedback effects on alpine cushion benefactors become more negative with increasing cover of graminoids and in dry conditions. Funct. Ecol. 30, 79–87. doi: 10.1111/1365-2435.12507
Pang, D. B., Wu, M. Y., Zhao, Y. R., Yang, J., Dong, L. G., Wu, X. D., et al. (2023). Soil microbial community characteristics and the influencing factors at different elevations on the eastern slope of Helan Mountain,Northwest China. Chin. J. Appl. Ecology. 34, 1 957–1 967. doi: 10.13287/j.1001-9332.202307.031
Pankratov, T. A., Ivanova, A. O., Dedysh, S. N., and Liesack, W. (2011). Bacterial populations and environmental factors controlling cellulose degradation in an acidic Sphagnum peat. Environ. Microbiol. 13, 1800–1814. doi: 10.1111/j.1462-2920.2011.02491.x
Peay, K. G., von Sperber, C., Cardarelli, E., Toju, H., Francis, C. A., Chadwick, O. A., et al. (2017). Convergence and contrast in the community structure of Bacteria, Fungi and Archaea along a tropical elevation–climate gradient. FEMS Microbiol. Ecol. 93, fix045. doi: 10.1093/femsec/fix045
Peura, S., Bertilsson, S., Jones, R. I., and Eiler, A. (2015). Resistant microbial cooccurrence patterns inferred by network topology. Appl. Environ. Microbiol. 81, 2090–2097. doi: 10.1128/AEM.03660-14
Pingree, M. R. A. and Deluca, T. H. (2018). The influence of fire history on soil nutrients and vegetation cover in mixed-severity fire regime forests of the eastern Olympic Peninsula, Washington,USA. For. Ecol. Manage. 422, 95–107. doi: 10.1016/j.foreco.2018.03.037
Qu, Z. L., Liu, B., Ma, Y., Xu, J., and Sun, H. (2020). The response of the soil bacterial community and function to forest succession caused by forest disease. Funct. Ecology. 34, 2548–2559. doi: 10.1111/1365-2435.13665
Ren, C. J., Zhang, W., Zhong, Z. K., Han, X. H., Yang, G. H., Feng, Y. Z., et al. (2018). Differential responses of soil microbial biomass, diversity, and compositions to altitudinal gradients depend on plant and soil characteristics. Sci. Total Environ. 610, 750–758. doi: 10.1016/j.scitotenv.2017.08.110
Ren, C. J., Zhou, Z. H., Guo, Y. X., Yang, G. H., Zhao, F. Z., Wei, G. H., et al. (2021). Contrasting patterns of microbial community and enzyme activity between rhizosphere and bulk soil along an elevation gradient. Catena. 196, 104921. doi: 10.1016/j.catena.2020.104921
Reyes, F., Lillo, A., Ojeda, N., Reyes, M., and Alvear, M. (2011). Influence of slope and orientation on biological activities of a relict temperate forest from South-Central Chile. Rev. Bosque 32, 255–265. doi: 10.4067/S0717-92002011000300007
Saiya-Cork, K. R., Sinsabaugh, R. L., and Zak, D. R. (2002). The effects of long term nitrogen deposition on extracellular enzyme activity in an Acer saccharum forest soil. Soil Biol. Biochem. 34, 1309–1315. doi: 10.1016/S0038-0717(02)00074-3
Salazar, A., Rousk, K., Jónsdóttir, I. S., Bellenger, J. P., and Andrésson, Ó.S. (2020). Faster nitrogen cycling and more fungal and root biomass in cold ecosystems under experimental warming: a meta-analysis. Ecology 101, e02938. doi: 10.1002/ecy.2938
Schimel, J. P. and Weintraub, M. N. (2003). The implications of exoenzyme activity on microbial carbon and nitrogen limitation in soil: a theoretical model. Soil Biol. Biochem. 35, 549–563. doi: 10.1016/S0038-0717(03)00015-4
Shen, C. C., Xiong, J. B., Zhang, H. Y., Feng, Y. Z., Lin, X. G., Li, X. Y., et al. (2013). Soil pH drives the spatial distribution of bacterial communities along elevation on Changbai Mountain. Soil Biol. Biochem. 57, 204–211. doi: 10.1016/j.soilbio.2012.07.013
Shu, X. Y., Hu, Y. F., Liu, W. J., Xia, L. L., Zhang, Y. Y., Zhou, W., et al. (2023). Linking between soil properties, bacterial communities, enzyme activities, and soil organic carbon mineralization under ecological restoration in an alpine degraded grassland. Front. Microbiol. 14. doi: 10.3389/fmicb.2023.1131836
Si, G. C., Wang, J., Xia, Y. Q., Yuan, Y. L., Zhang, G., and Lei, T. (2014). Change characteristics of microbial communities and enzyme activities in soils of marshes in Nyaiqentanglha Mountains with heights above sea level. Wetland Science. 12, 340–348. doi: 10.13248/j.cnki.wetlandsci.2014.03.011
Siles, J. A., Cajthaml, T., Filipova, A., Minerbi, S., and Margesin, R. (2017). Altitudinal, seasonal and interannual shifts in microbial communities and chemical composition of soil organic matter in Alpine forest soils. Soil Biol. Biochem. 112, 1–13. doi: 10.1016/j.soilbio.2017.04.014
Singh, D., Takahashi, K., Kim, M., Chun, J., and Adams, J. M. (2012). A hump-backed trend in bacterial diversity with elevation on Mount Fuji, Japan. Microbial ecology. 63, 429–437. doi: 10.1007/s00248-011-9900-1
Sinsabaugh, R. L., Hill, B. H., and Shah Jennifer, J. F. (2009). Ecoenzymatic stoichiometry of microbial organic nutrient acquisition in soil and sediment. Nature. 462, 795–798. doi: 10.1038/nature08632
Sun, P. C., Hao, R. R., Fan, F. J., Wang, Y., and Zhu, F. Y. (2025). Adaptation of high-altitude plants to plateau abiotic stresses: A case study of the qinghai-tibet plateau. Int. J. Mol. Sci. 26, 2292. doi: 10.3390/ijms26052292
Truong, C., Gabbarini, L. A., Corrales, A., Mujic, A. B., Escobar, J. M., Moretto, A., et al. (2019). Ectomycorrhizal fungi and soil enzymes exhibit contrasting patterns along elevation gradients in southern Patagonia. New Phytologist. 222, 1936–1950. doi: 10.1111/nph.15714
Uselman, S. M., Qualls, R. G., and Thomas, R. B. (2000). Effects of increased atmospheric CO2, temperature, and soil N availability on root exudation of dissolved organic carbon by a N-fixing tree (Robinia pseudoacacia L.). Plant Soil. 222, 191–202. doi: 10.1023/A:1004705416108
Wang, B. R., An, S. S., Liang, C., Liu, Y., and Kuzyakov, Y. (2021). Microbial necromass as the source of soil organic carbon in global ecosystems. Soil Biol. Biochem. 162, 108422. doi: 10.1016/j.soilbio.2021.108422
Wang, L., Li, Z. M., Gao, C. J., Li, J., and Cui, K. (2025). The altitude shapes the heterogeneity of tree growth by modulating the soil characteristics and microorganisms in the Calocedrus macrolepis plantation. For. Ecol. Manage. 578, 122459. doi: 10.1016/j.foreco.2024.122459
Wang, R., Wang, M., Wang, J., Yao, J. H., Li, X. W., Lin, Y. H., et al. (2021). Soil bacterial characteristics under four habitats with different vegetation communities on the Qinghai-Tibetan Plateau. Wetlands 41, 58. doi: 10.1007/s13157-021-01455-0
Wyszkowska, J., Borowik, A., Olszewski, J., and Kucharski, J. (2019). Soil bacterial community and soil enzyme activity depending on the cultivation of Triticum aestivum, Brassica napus, and Pisum sativum ssp. arvense. Diversity 11, 246. doi: 10.3390/d11120246
Xi, N. X., Zhu, B. R., and Zhang, D. Y. (2017). Contrasting grass nitrogen strategies reflect interspecific trade-offs between nitrogen acquisition and use in a semi-arid temperate grassland. Plant Soil. 418, 267–276. doi: 10.1007/s11104-017-3296-8
Xie, L. L., Ma, Y. S., Wang, Y. L., Ma, Y., and Wang, X. L. (2024). Changes in soil bacterial and fungal community composition and functional groups during the artificial restoration of degraded grassland of “Black-soil mountain. Ecol. Evol. 14, e70361. doi: 10.1002/ece3.70361
Xu, X. F., Thornton, P. E., and Post, W. M. (2013). A global analysis of soil microbial biomass carbon, nitrogen and phosphorus in terrestrial ecosystems. Global Ecol. Biogeography 22, 737–749. doi: 10.1111/geb.12029
Yan, R. Y., Zhao, X. M., Li, P. H., Si, Z. Y., Gao, Y., and Li, J. F. (2023). Composition and diversity of soil microbial communities in walnut orchards at different altitudes in southeastern Tibet. Land 12, 1419. doi: 10.3390/land12071419
Yang, C. D., Long, R. J., Chen, X. R., and Xue, L. (2008). Advanced research in soil microbial functional groups. Chin. J. Soil Sci. 02), 421–425. doi: 10.19336/j.cnki.trtb.2008.02.043
Yang, N., Zou, D. S., Yang, M. Y., Zhao, L. F., Song, G. T., and Lin, Z. G. (2013). The change of soil microbial biomass carbon and the relationship between it and soil physio-chemical factors in different restoration stages on sloping-land with purple soils in Hengyang. Ecol. Environ. Sci. 22, 25–30. doi: 10.16258/j.cnki.1674-5906.2013.01.014
Yao, F., Yang, S., Wang, Z. R., Wang, X., Ye, J., Wang, X. G., et al. (2017). Microbial taxa distribution is associated with ecological trophic cascades along an elevation gradient. Front. Microbiol. 8. doi: 10.3389/fmicb.2017.02071
Yuan, Y. L., Si, G. C., Wang, J., Luo, T. X., and Zhang, G. X. (2014). Bacterial community in alpine grasslands along an altitudinal gradient on the Tibetan Plateau. FEMS Microbiol. ecology. 87, 121–132. doi: 10.1111/1574-6941.12197
Zhang, Q. F. and Laanbroek, H. J. (2020). Tannins from senescent Rhizophora mangle mangrove leaves have a distinctive effect on prokaryotic and eukaryotic communities in a Distichlis spicata salt marsh soil. FEMS Microbiol. ecology. 96, fiaa148. doi: 10.1093/femsec/fiaa148
Zhao, L. R., Li, K. X., Fang, N. F., Jia, S. Y., Cui, Z., Liu, Y. F., et al. (2025). Shrub encroachment and degradation impact water balance by altering water fluxes in alpine meadows. Catena 258, 109231. doi: 10.1016/j.catena.2025.109231
Zhao, W., Yin, Y. L., Li, S. X., Dong, Y. L., and Su, S. F. (2023). Changes in soil fungal community composition and functional groups during the succession of Alpine grassland. Plant Soil. 484, 201–216. doi: 10.1007/s11104-022-05781-y
Zhen, Z., Wang, S. B., Luo, S. W., Ren, L., Liang, Y. Q., Yang, R. C., et al. (2019). Significant impacts of both total amount and availability of heavy metals on the functions and assembly of soil microbial communities in different land use patterns. Front. Microbiol. 10. doi: 10.3389/fmicb.2019.02293
Zhou, H. (2019). Characteristics and adaptation mechanism of soil microbial communities of alpine steppe under different altitudes on Qinghai-Tibetan Plateau. Gansu Agric. University, 29–31. doi: 10.27025/d.cnki.ggsnu.2019.000033
Zhou, S., Ma, H. F., Wang, J., Ren, C. J., Guo, Y. X., Wang, J., et al. (2022). Latitudinal distribution of forest soil microbial biomass carbon and its affecting factors in China. SCIENTIA Silvae SINICAE 58, 49–57. doi: 10.11707/j.1001-7488.20220206
Zhou, X. B., Wang, X. L., Wang, Y. L., Ma, Y., Liu, Y., and Ma, Y. S. (2024). Vegetation restoration has an implication for fungal diversity and composition in a degraded temperate desert type rangeland of China. Ecol. Eng. 207, 107348. doi: 10.1016/j.ecoleng.2024.107348
Zhu, N., Li, K. X., Gao, J. M., Wei, B. M., Zhao, L. R., Liu, L., et al. (2025). Vegetation degradation dominates over elevation in structuring fungal communities in alpine meadows. Front. Microbiol. 16. doi: 10.3389/fmicb.2025.1596407
Keywords: alpine ecosystem, bacterial community, different altitudes, co-occurrence network, soil properties
Citation: Xie L, Ma Y, Wang Y, Ma Y and Liu Y (2025) Changes of soil bacterial community composition and functional groups in different altitude gradients of Potentilla fruticosa shrub in eastern Qinghai-Tibet Plateau. Front. Plant Sci. 16:1539945. doi: 10.3389/fpls.2025.1539945
Received: 05 December 2024; Accepted: 12 August 2025;
Published: 20 October 2025.
Edited by:
Yongjun Tan, Hunan Hybrid Rice Research Center, ChinaReviewed by:
Yu Liu, Northwestern Polytechnical University, ChinaLyudmyla Symochko, University of Coimbra, Portugal
Copyright © 2025 Xie, Ma, Wang, Ma and Liu. 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: Lele Xie, MjY3ODYzMzA5OUBxcS5jb20=; Yushou Ma, bWF5dXNob3VAc2luYS5jb20=; Yanglong Wang, d2FuZ3lsNTA2QDE2My5jb20=