Variation in soil bacterial community characteristics inside and outside the West Ordos National Nature Reserve, northern China

Nature reserves are crucial for protecting biological habitats and maintaining biodiversity. Soil bacterial community plays an irreplaceable role in the structure and function of ecosystem. However, the impact of nature reserves on soil bacterial communities is still unclear. To explore the effects of desert grassland nature reserve management on soil microbial communities, we compared the differences in soil bacterial community composition, α-diversity and community structure inside and outside a desert grassland nature reserve, and explored the correlation between soil bacterial communities and plant biomass and soil chemical index. We found that (1) the relative abundance of Acidobacteriota is highest in the soil both inside and outside the nature reserve in shrub grassland; (2) the Chao1 index of soil bacterial communities in the core protected zone and general control zone of the reserve was significantly higher than that outside the reserve (p < 0.05) in the shrub grassland. Similarly, in the herbaceous grassland, the Shannon index of soil bacterial communities was significantly higher in the core protected zone of the reserve than that outside the reserve (p < 0.05). (3) While we found no significant difference in soil bacterial community structure between inside and outside the reserve in the shrub grassland, we found that the soil bacterial community structure in the core protected zone was significantly different from that outside the reserve in the herbaceous grassland (p < 0.05); (4) we also found that higher plant productivity and soil nutrients promoted most soil dominant bacterial phyla, while higher soil pH and salinity inhibited most soil dominant bacterial phyla. Our findings thus help better understand the influencing factors of and the mechanisms behind variation in soil bacterial communities inside and outside desert grassland nature reserves.


Introduction
Desert grassland is the driest type of grassland vegetation.The desert grassland in Inner Mongolia Autonomous Region serves as a crucial base for animal husbandry and an important ecological security barrier in northern China (Guo et al., 2022).A nature reserve that focuses on protecting the natural ecosystem of desert biotic and abiotic Guo et al. 10.3389/fmicb.2024.1404848Frontiers in Microbiology 02 frontiersin.orgenvironments is known as a desert ecosystem type nature reserve (Ma et al., 2023).These reserves are situated in fragile environments that suffer from issues such as overgrazing, drought, and land degradation (Sankey et al., 2015;Xiong et al., 2016;Zuo et al., 2022).Therefore, it is crucial to enhance the management and protection of these reserves.The West Ordos National Nature Reserve is home to a variety of rare plants, such as Tetraena mongolica, Helianthemum songaricum, Potaninia mongolica, Reaumuria songarica, which hold significant ecological value (Li E. G. et al., 2018;Dong et al., 2023).Soil microorganisms play a crucial role in soil ecosystems (Griffiths et al., 2008;Xia et al., 2020).They reflect the health of soil ecosystems (Agbenin and Adeniyi, 2014) and are essential for soil material cycling, energy flow, and promoting plant growth and development (Mlambo et al., 2007;Zhang et al., 2023).Soil bacteria are a crucial component of soil microorganisms.Numerous studies have been conducted on the characteristics of soil bacterial communities, with a primary focus on the impact of natural and anthropogenic factors on their structure and function.These factors include climate change, plant community succession, changes in land use types, alterations in the physical and chemical properties of the soil (Li G. et al., 2018;Ma S. et al., 2018;Chai et al., 2019;Sünnemann et al., 2021).For instance, in alpine meadows, soil bacterial biomass typically first decreases and then increases with grazing intensity (Yang et al., 2018).In semi-arid meadows, heavy grazing may significantly reduce soil bacterial biomass (Liu et al., 2016).These findings suggest that the impact of human activities on soil bacterial communities could vary with vegetation type.Studying the role and mechanism of nature reserves in protecting the structure and function of soil bacteria can guide the improvement and development of conservation measures.Some studies have focused on the impact of nature reserves on plants and animals, while others have studied changes in land use and ecosystem services within and outside nature reserves, and some studies have found that nature reserves can enhance local plant diversity and ecosystem services (Tranquilli et al., 2012;Carranza et al., 2013;Gray et al., 2016;Yu et al., 2024).However, few studies have compared the variation in characteristics of soil bacterial communities inside and outside nature reserves.Even a comparative study of soil bacterial communities in core protected areas and general control areas in nature reserves has not been reported.Based on previous studies, nature reserves enhance local plant diversity and ecosystem services, which may have a positive impact on soil bacterial communities.However, the impact and mechanism of management measures in nature reserves on soil microorganisms are still unclear.Therefore, studying the structure and function of soil bacterial communities inside and outside nature reserves is significant for the sustainable management and optimization of soil ecosystems (Atwell et al., 2018).
In this study, we systematically and comprehensively compare the differences in soil bacterial community characteristics between inside and outside the West Ordos National Nature Reserve using ground survey data.We also explore the relationship between soil bacteria and soil physicochemical properties and plant productivity.Our aim is to provide a reference basis for scientifically and objectively evaluating the effectiveness of soil ecosystem protection in desert nature reserves in arid zones and upgrading the management level of the nature reserves.

Study site
The research area is situated in the West Ordos National Nature Reserve (106°40′-107°44′ E, 39°14′-40°11′ N, 1,100-2,000 m asl), Otog Banner, Inner Mongolia Autonomous Region (Figure 1).The area has an arid and semi-arid continental climate in the mesotemperate zone.The winter is long and cold, the summer is short and warm, the spring is windy, and the fall has stable weather.The frost-free period in this area is 129 days.The average annual temperature is 7°C, with the lowest temperature recorded at −36.8°C and the highest at 37°C.Annual precipitation averages at 272.3 mm, with uneven distribution throughout the year, mostly concentrated in June-September.Annual evaporation is 2470.4mm, which is 9.1 times the amount of precipitation.The average annual wind speed is 3 m/s.The average wind speed in the area is 3.4 m/s, with an average of 36 days of winds of grade 6 or above in a calendar year.Of these days, 19, or 52.5%, occur in spring.The soil is mainly chestnut-calcium soil.

Sample selection and setting
Nature reserves are divided into core protected zones and general controlled zones.The core protected zones are primarily responsible for limiting human activities to the greatest extent possible in order to fulfill their protection function.The general controlled zones, on the other hand, are responsible for fulfilling their protection function while also considering public service functions such as scientific research, education, and recreation.In this research area, the core protected zone is closed and human activities are strictly prohibited.Necessary human activities, such as grazing, may be limited to the general control zone.Production activities outside the reserve are not restricted.The primary human activities in this study area are cattle and sheep grazing and crop planting.The plant biomass and diversity of the core protected area in the nature reserve are higher and more abundant, followed by the general control area, while the plant biomass and diversity outside the nature reserve are lower than those inside the protected area.The degree of soil desertification outside the protected area is more serious than that inside the protected area, and there is obvious soil salinization.The survey was conducted in August 2021, with sample plots selected from the core protected zone, general controlled zone, and outside the reserve.And the three areas are 30 km apart from each other.Two vegetation types, shrub grassland and herbaceous grassland, were surveyed in each sample site.Soil surveys were conducted in three randomly selected squares of each planting type, separated by 500 meters from each other, for a total of 18 squares.The "S"-shaped five-point mixed sampling method was used to collect soil samples from each plot at depths of 0 to 20 cm, and the samples were then mixed to create a homogeneous sample of the site (Wei et al., 2024).The soil samples were divided into two portions.One portion was used to determine soil physicochemical properties, while the other was stored in a refrigerator at −80°C after impurities were removed.This portion was used for high-throughput sequencing of soil bacteria (Teng et al., 2018).

Sample collection and determination method
The biomass survey for herbaceous plants involved cutting the aboveground portion of all grasses in the sampling plots by species.The cut grasses were then brought back to the laboratory, dried in an oven at 65°C to a constant weight, and weighed.For shrubs, standard plants of representative sizes and heights were selected from the surveyed sampling plots.The aboveground portion of the standard plant was cut near the ground using scissors.The samples were then placed in paper bags and taken to the laboratory for biomass dry weight determination using the drying method.Sample biomass was calculated using the formula: biomass of shrubs = biomass of standard plants × number of shrub plants.Following plant sampling, a single root sample measuring 0-40 cm was taken from each sample plot using a root drill.The root samples were rinsed with water using a 2 mm sieve.After rinsing, the root samples were dried in an oven at 65°C until a constant weight was achieved.The weight was then used to calculate the belowground biomass of the plants.Soil pH, salinity, organic carbon, total nitrogen, and effective phosphorus were measured using a pH meter, conductivity analysis, oxidation of potassium dichromate with external heating, automatic nitrogen fixation, and ICP-AES analysis, respectively (Crepin and Johnson, 1993).

Extraction, sequencing and bioinformatics analysis of soil bacterial DNA
In the laboratory, soil samples were subjected to DNA extraction and PCR amplification.The PowerSoil DNA extraction kit (MoBio Inc., United States) was used to extract total soil DNA from 0.3-0.5 g of soil samples, following the experimental steps specified in the kit.The gene fragments of the V3-V4 region were amplified using primers 338F (5′-ACTCCTACGGGGAGGCAGCA-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′).The PCR reaction system consisted of 5 × FastPfu Buffer (4.0 μL), dNTPs (2.5 mmol/L) (2.0 μL), forward and reverse primers (5 μmol/L) (0.8 μL each), FastPfu Polymerase (2.5 U/μL) (0.4 μL), BSA (0.2 μL), Template DNA (10 ng), and ddH 2 O (supplemented to 20 μL).PCR reaction conditions were as follows: 95°C for 5 min; 25 cycles of 95°C for 30 s, 55°C for 30 s, and 72°C for 40 s; 72°C for 7 min, and 4°C for infinite conservation.The amplification products were loaded onto a 3% agarose gel and visualized using a gel imaging system.No changes in content were Location of the study area.Guo et al. 10.3389/fmicb.2024.1404848Frontiers in Microbiology 04 frontiersin.orgmade.The PCR products were then quantified using the Qubit fluorescence quantification system and mixed in the appropriate proportion for sequencing.MiSeq library construction was performed followed by online sequencing using Illumina HiSeq2500 (Liu et al., 2014).The data for each sample was split based on the barcode sequence.The barcode and primer sequences were then spliced using FLASH (v1.2.7) pairs (Magoc and Salzberg, 2011).The highthroughput sequence data were processed using QIIME (v1.7.0).Valid sequences were clustered using UPARSE (v7.0) (Caporaso et al., 2013).The sequences were clustered to obtain operational taxonomic units (OTUs) using a default 97% sequence similarity.From each OTU, one representative sequence was selected, and the species annotation was performed using the RDP Classifier with the Greengenes database as a reference at a confidence level of 80% (Caporaso et al., 2010).Finally, the sample data were standardized to the lowest number in the sample.

Statistical analysis
We use relative abundance to characterize the dominance of soil bacterial phyla and orders.The higher the relative abundance of soil bacterial phyla and orders, the higher the degree of dominance.In this study, the top 10 bacterial phyla and order in relative abundance were selected as the dominant bacteria.The soil bacterial community's α-diversity was measured using the Chao1 index and Shannon-Wiener index (Jones, 2010).The Chao1 index was calculated using the formula: Chao1 index: F "S obs " denotes the number of OTUs observed in the sample."F 1 " and "F 2 " denote the number of OTUs when there is only one sequence and when there are two sequences, respectively.
The Shannon-Wiener index was calculated using the formula: Shannon-Wiener index: "S obs " denotes the number of OTUs observed in the sample."n i " is the number of OTUs containing "i" sequences."N" is the number of OTUs for all sequences.
The Chao1 index is used to estimate the number of OTUs that are present in a sample.The Shannon-Wiener index is used to measure the diversity of a community.In the formula, n is the total number of bacteria of a particular species and N is the total number of bacteria of all species.The α-diversity was calculated using R and single variance analysis of soil bacterial alpha diversity index was performed in different regions, and multiple comparisons were made with Tukey HSD.Principal component analysis was used to analyze the differences in soil bacterial community structure between core protected areas, general control areas and outside protected areas, and permutational multivariate analysis of variance was used for statistical testing.Spearman's method was used to perform correlation analysis, including positive and negative correlation, based on the abundance of each species and the changes.Statistical tests were then conducted to filter out data groups with absolute values of correlation coefficients greater than 0.1 and p ≤ 0. 05.In this study, correlation network analysis graphs were created using the Python language to demonstrate the top twenty phyla with the highest correlation.The analysis was performed in Cytoscape 2.6 (Otasek et al., 2019).We define the bacterial phylum with the largest number of relationships with other bacterial phyla as the core bacterial phyla.If some bacterial phyla have the same number of correlations with other bacterial phyla, the one with the highest relative abundance is selected as the core bacterial phyla.A correlation heatmap was created using Spearman's correlation coefficient to analyze the correlation between the top 10 dominant bacterial phyla and orders in soil and both plant productivity and soil factors.In the shrub grassland, the top three dominant bacterial phyla of the soil in the core protected zone, general controlled zone, and outside the protected areas were Acidobacteriota (29.28 ± 2.03, 27.22% ± 1.22, 29.21% ± 2.42%), Actinobacteriota (21.47% ± 1.85, 26.36% ± 2.61, 27.11% ± 2.61%) and Proteobacteriota (15.23% ± 2.00, 14.80 ± 4.35, 14.07% ± 4.14%), respectively.In herbaceous grassland, the top 10 dominant bacterial phyla in the soil of the core protected zone were Acidobacteriota (29.21% ± 2.42%), Actinobacteriota (21.99% ± 1.93%), and Proteobacteria (15.61% ± 1.26%); and the dominant bacterial phyla in the soil of the general controlled zones and the soils outside the protected areas were Actinobacteriota (38.58% ± 3.84 and 32.36% ± 4.76%), Acidobacteriota (23.35% ± 1.45 and 26.06% ± 3.65%) and Proteobacteria (10.58% ± 2.09 and 11.48% ± 1.39%) in general control areas and outside protected areas, respectively (Figure 2A).

Alpha diversity of the soil bacterial community
Figure 3 displays the α-diversity indices of soil bacterial communities inside and outside the protected areas.In the shrub grassland, the Chao1 indices of soil bacterial communities in both core protected zones and general controlled zones were significantly higher than those outside the protected areas (p < 0.05) (Figure 3A).The soil bacterial communities' Shannon index was higher in both core protected zones and general controlled zones than outside the protected areas (Figure 3B).In herbaceous grassland, the Chao1 index of soil bacterial communities was higher in core protected zones than in general controlled zones and outside protected zones (Figure 3C).The Shannon index indicated that the soil bacterial community in core protected zones had a significantly higher diversity than that in outside protected zones (p < 0.05), with general controlled zones falling in between (Figure 3D).

Principal component analysis of soil bacterial community structure
The soil bacterial community structure in the shrub grassland was not significantly different inside and outside the protected area, although the sum of the contribution of the first and second principal components to the difference in soil bacterial community structure was 73.91% (Figure 4A).The first and second principal components contributed to 59.59% of the difference in soil bacterial community structure in the herbaceous grassland (Figure 4B).There was a significant difference (p < 0.05) between the soil bacterial community structure inside and outside the protected area, indicating that the establishment of the protected area had a significant impact on the soil bacterial community structure of the herbaceous grassland.

Correlation network analysis of soil bacterial phyla
To investigate the differences in soil bacterial community structure inside and outside the nature reserve, we compared the correlation networks of the top 20 bacterial phyla in relative abundance.In shrub grassland, the core phylum of the soil bacterial network in the core protected zones, general controlled zones and outside the protected areas were Acidobacteriota, Proteobacteria and Gemmatimonadota, respectively (Figures 5A,C,E).There were 29, 31, and 31 sets of positive correlations and 37, 30, and 32 sets of negative correlations between soil bacterial phyla in core protected areas, general control areas, and outside protected areas, respectively.There was a significant positive correlation between soil Acidobacteriota and Nitrospirota in core protected areas and general control areas within nature reserves (Figures 5A,C), while there was no significant correlation between the two in soils outside protected areas (Figure 5E).
In herbaceous grassland, the core phyla of the soil bacterial network in the core protected zones, general controlled zones and outside the protected areas were Acidobacteriota, Acidobacteriota and Actinobacteriota, respectively (Figures 5B,D,F).There were 31, 33, and 30 sets of positive correlations and 30, 28, and 39 sets of negative correlations between soil bacterial phyla in the core protected area, general control area, and outside the protected area, respectively.There was a significant negative correlation between soil Actinobacteriota and Acidobacteriota outside the protected area (Figures 5B,D), but there was no significant correlation between the two in the soil of the core protected area and the general control area within the nature reserve (Figure 5F).

Correlation between soil bacterial communities and plant biomass and soil chemical indices
Figure 6 shows the relationship between the top 10 bacterial phyla in abundance and plant biomass, as well as soil factors.There was a significant positive correlation between Bacteroidota in soil and both aboveground and belowground plant biomass.In contrast, Chloroflexi showed a negative correlation with both plants aboveground and belowground biomass.Actinobacteriota was negatively correlated only with aboveground plant biomass.The aboveground and belowground biomass of plants showed a negative correlation with Chloroflexi.Meanwhile, Actinobacteriota and Chloroflexi exhibited a positive correlation with soil salt content.On the other hand, Bacteroidota, Nitrospirota, Proteobacteria, and Acidobacteriota showed a negative correlation with soil salt content.Additionally, there was a significant positive correlation between Bacteroidota and soil organic carbon.Lastly, Chloroflexi exhibited a negative correlation with soil organic carbon.The abundance of Bacteroidota and Nitrospirota showed a positive correlation with soil total nitrogen, while Chloroflexi exhibited a negative correlation.No significant correlation was observed between the dominant bacteria and soil pH or available phosphorus (Figure 6).The relationship between the 10 most abundant bacterial orders and plant biomass, as well as soil factors, is illustrated in Figure 7.Both Solirubrobacterales and Gaiellales in soil were significantly positively correlated with soil salt content; while both Rhizobiales and Burkholderiales were significantly negatively correlated with soil salt content.Gaiellales was significantly negatively correlated with soil organic carbon.Both Solirubrobacterales and Gaiellales were significantly negatively correlated with soil total nitrogen.

Differences in dominant bacterial phyla, genus and α-diversity of soil bacterial communities in and outside nature reserves
In this study, the soil dominant bacterial phyla in the nature reserve were more similar to that outside the nature reserve, which were Acidobacteriota, Actinobacteriota, Bacteroidota, Chloroflexi, Gemmatimonadota, Methylomirabilota, Myxococcota, Nitrospirota, Proteobacteria, and unclassified_Bacteria; the relative abundance of each dominant phylum varied slightly.In herbaceous grassland, the highest abundance of soil Acidobacteriota was found in the core protected zone within the protected zone (29.21% ± 2.42%), while the Comparison of chao1 index and Shannon index of shrub grassland (A,B) and herbaceous grassland (C,D) inside and outside the nature reserves.SC, shrub grassland in core protected zone; SG, shrub grassland in general controlled zone; SO, shrub grassland outside the protected area; HC, herbaceous grassland in core protected zone; HG, herbaceous grassland in general controlled zone; HO, herbaceous grassland outside the protected area. 10.3389/fmicb.2024.1404848 Frontiers in Microbiology 07 frontiersin.orghighest abundance of soil Actinobacteriota was found in the general controlled zones and the soil outside the protected areas (38.58% ± 3.84 and 32.36% ± 4.76%, respectively), which suggests that the strict protective measures altered relative abundance of soil dominant phyla in herbaceous grassland.Acidobacteriota plays an important role in degrading plant residues and participating in ecosystem iron cycling (Eichorst et al., 2011;Stursová et al., 2012).The increase in the relative abundance of soil Acidobacteriota under strict protection measures favors the decomposition and transformation of organic matter in plant communities and promotes soil nutrient cycling.Species diversity is an indispensable factor in maintaining ecosystem function, and for soil ecosystems, soil bacterial community diversity plays an irreplaceable role in promoting soil material cycling and energy flow (Lladó et al., 2017;Yu et al., 2021).The Chao1 index is a crucial metric for characterizing the species richness of soil bacterial communities (Wang et al., 2023).Our study revealed that in shrub grasslands, the Chao1 index of soil bacterial communities in core protected zones and general controlled zones was significantly higher than that outside the protected areas (p < 0.05).This increase can be attributed to the protective measures implemented by nature reserves.Human activities, such as grazing, were effectively restricted in the protected areas of desert grassland.This restriction led to the effective restoration of vegetation and soil (Huang et al., 2021).As a result, soil bacterial communities were able to thrive, leading to an increase in species richness.The Shannon index of soil bacterial communities in the herbaceous grassland was higher in the core protected zone than outside of it (p < 0.05).This suggests that management measures in protected zones have a positive effect on improving the diversity of soil bacterial communities.However, the changes in the patterns of the Chao1 index and Shannon index of soil bacterial communities inside and outside the protected areas under different vegetation types (shrub grassland and herbaceous grassland) were not entirely consistent.This suggests that there are differences in the species richness and diversity of soil bacterial communities among different vegetation types under the same management measures.These differences may be related to the protective effect of scrub, the "fertility islands effect, " and micro environmental improvement (Gregory et al., 2001;Chen et al., 2015).Further investigation is needed to understand the mechanisms behind these differences.
4.2 Comparison of differences in soil bacterial community structure inside and outside the nature reserves The analysis of principal components reveals that the soil bacterial community structure in the shrub grassland within the nature reserve did not differ significantly from that outside the reserve.The distance between each sample point was relatively far, which may be attributed to the spatial heterogeneity of the shrub grassland (Cheng et al., 2024).In herbaceous grassland, the soil bacterial community structure inside the protected area differed significantly from that outside, indicating that the management measures implemented in the protected area had a significant impact on the soil bacterial community structure of herbaceous grassland.This finding is consistent with the results of the study conducted by Esposito (Esposito et al., 2024).In desert grasslands, the primary human activity is grazing on herbaceous grasses.Overgrazing can damage both the vegetation and soil structure and function (Milazzo et al., 2023), which can then impact the structure of the soil bacterial community.Taxa with high connectivity are typically considered core taxa in the network and play a crucial role in the bacterial network (Ma B. et al., 2018).The correlation network analysis among soil dominant bacterial phyla in this study reveals that in shrub grassland, the core phyla of soil bacterial networks in core protected zones, general controlled zones, and outside protected zones were Acidobacteriota, Proteobacteria, and Gemmatimonadota, respectively.In herbaceous grassland, the soil bacterial networks in core protected zones, general controlled zones, and outside protected zones of nature reserves were Acidobacteria, Proteobacteria, and Gemmatimonadota, respectively.In herbaceous Analysis of the principal components of shrub grassland (A) and herbaceous grassland (B) inside and outside nature reserves.SC, shrub grassland in core protected zone; SG, shrub grassland in general controlled zone; SO, shrub grassland outside the protected area; HC, herbaceous grassland in core protected zone; HG, herbaceous grassland in general controlled zone; HO, herbaceous grassland outside the protected area.The soil bacterial network in the core protected zone (A), general controlled zone (C) and outside the protected areas (E) in shrub grassland and the soil bacterial network in the core protected zone (B), general controlled zone (D) and outside the protected areas (F) in herbaceous grassland. 10.3389/fmicb.2024.1404848 Frontiers in Microbiology 09 frontiersin.orggrasslands, the Acidobacteriota phylum dominated the core protected zones, while the general control zones were dominated by Acidobacteriota as well.Actinobacteriota was the dominant phylum in the outside protected zones of nature reserves.The core phyla of the soil bacterial network structure changed as the protection measures became stricter (Qiu et al., 2022).The core phylum of the soil bacterial network structure of different vegetation types in the core protected zone was Acidobacteriota, while the core phylum of the soil bacterial network structure of different vegetation types in the general controlled zone and outside the protected areas was different, which may be due to the different intensities of human activities in the two zones (Hua et al., 2024).The correlation between various types of soil bacteria also changes with the environment, and it has been found that the correlation between various types of soil bacteria in different stages of soil development showed significant differences due to different environments (Dini-Andreote et al., 2014).In this study, there were 29, 31, and 31 sets of positive correlations and 37, 30, and 32 sets of negative correlations between soil bacterial phyla in core protected areas, general control areas, and outside protected areas, respectively.However, in herbaceous grassland, the correlations between soil bacterial phyla were different of that in shrub grassland, which was related to the habitat differences between shrub and herbaceous grassland.We need to collect more samples to improve further research.And we need to conduct in-depth analysis of correlation networks in the future, including conducting research at the OTU level.

Correlation between soil bacterial communities and plant biomass and soil chemical indices
Research has demonstrated that an increase in litter amount of plant community and belowground biomass within plant communities can lead to an increase in soil nutrient content and promote soil microbial diversity (Mao et al., 2024).Microorganisms require nutrients mainly derived from soil.Studies have shown that when soil nutrients are abundant, the number and biomass of microbial communities are positively correlated with soil organic carbon and nitrogen, as well as total soil nitrogen (He et al., 2023).Additionally, these studies have found a correlation between soil bacterial communities and plant productivity.Our study revealed that Bacteroidota, Nitrospirota, Proteobacteria, and Acidobacteriota are significantly negatively correlated with soil salinity.This suggests that high salinity soil environments have a negative effect on these bacterial phyla, which is unfavorable to their growth and development, and negatively affects the structure of the entire bacterial community.Bacteroidota exhibited significant positive correlations with aboveground biomass, belowground biomass, soil organic carbon, and soil total nitrogen.The vegetation and soil restoration measures were found to be favorable for the survival and development of Bacteroidota.In this study, Chloroflexi exhibited significant negative correlations with plant aboveground biomass, belowground biomass, and soil organic carbon.This finding contrasts with the results of (Zhang et al., 2024) in Jilin Momog, which may be due to differences in the response of soil Chloroflexi to plant productivity and soil factors in arid versus humid environments.We found that higher plant productivity and soil nutrients had a positive effect on most soil dominant phyla, while higher soil pH and salinity had a negative effect on most soil dominant phyla.

Conclusion
The species richness and diversity of soil bacterial communities are higher inside than outside the West Ordos National Nature Reserve.Higher plant productivity and soil nutrients have positive effects on most soil dominant bacterial phyla, while higher soil pH and salinity have negative effects on most soil dominant bacterial phyla.Therefore, improving the management of nature reserves and implementing measures to restore saline and alkaline lands in protected areas can

FIGURE 1
FIGURE 1 of dominant bacterial phyla, orders and α-diversity of soil bacterial communities in and outside nature reserves 3.1.1Dominant phyla and orders in the soil bacterial community

FIGURE 2
FIGURE 2Distribution of soil dominant bacterial phyla (A) and orders (B) inside and outside protected areas.SC, shrub grassland in core protected zone; SG, shrub grassland in general controlled zone; SO, shrub grassland outside the protected area; HC, herbaceous grassland in core protected zone; HG, herbaceous grassland in general controlled zone; HO, herbaceous grassland outside the protected area.
The main plant species in the area include Zygophyllum xanthoxylum, Tetraena mongolica, Helianthemum songaricum, Reaumuria songarica, and Stipa tianschanica.Established in 1995, West Ordos National Nature Reserve primarily protects endangered plants and desert ecosystems, such as T. mongolica and H. songaricum.It is classified as a desert ecosystem nature reserve.