Effects of precipitation change and nitrogen addition on the composition, diversity, and molecular ecological network of soil bacterial communities in a desert steppe

Currently, the impact of changes in precipitation and increased nitrogen(N) deposition on ecosystems has become a global problem. In this study, we conducted a 8-year field experiment to evaluate the effects of interaction between N deposition and precipitation change on soil bacterial communities in a desert steppe using high-throughput sequencing technology. The results revealed that soil bacterial communities were sensitive to precipitation addition but were highly tolerant to precipitation reduction. Reduced precipitation enhanced the competitive interactions of soil bacteria and made the ecological network more stable. Nitrogen addition weakened the effect of water addition in terms of soil bacterial diversity and community stability, and did not have an interactive influence. Moreover, decreased precipitation and increased N deposition did not have a superimposed effect on soil bacterial communities in the desert steppe. Soil pH, moisture content, and NH4+-N and total carbon were significantly related to the structure of bacterial communities in the desert steppe. Based on network analysis and relative abundance, we identified Actinobacteria, Proteobacteria, Acidobacteria and Cyanobacteria members as the most important keystone bacteria that responded to precipitation changes and N deposition in the soil of the desert steppe. In summary, we comprehensively analyzed the responses of the soil bacterial community to precipitation changes and N deposition in a desert steppe, which provides a model for studying the effects of ecological factors on bacterial communities worldwide.


Introduction
Change in precipitation and increase in N deposition are current global issues [1,2] that greatly affect the composition of terrestrial ecosystems. The global deposition of reactive N was 100 Tg N year -1 in 1995, which is predicted to increase to 200 Tg N year -1 by 2050 due to industrial pollution and agricultural practices [3,4]. Further, the effects of climate change are becoming obvious, especially the change in precipitation [5]. It is predicted that precipitation will increase in high latitudes and decrease in most subtropical regions [5]. In fact, compared with other grasslands, changes in precipitation have a greater impact on desert grassland ecosystems [6]. Compared with other steppe types, desert steppe is less stable and highly sensitive to disturbances due to climate change [7,8]. Water and N are key limiting factors of the survival of organisms inhabiting a desert steppe and are also two coupling factors in steppe ecosystems [9,10]. The dissolution of inorganic N is highly dependent on water [11][12][13]. Further, soil N addition enhances the photosynthetic capability of plants [14][15][16]. Therefore, understanding the responses of desert steppe ecosystems to changes in water and N content will play important roles in elucidating complex ecosystems and predicting the responses of ecosystems to global change.
Soil microbes are key components of below-ground ecosystems, as their diversity, composition, and activities are the major drivers of terrestrial ecosystem productivity and diversity [17][18][19]. Some studies have reported that precipitation and nitrogen deposition indirectly affect microbial communities by affecting plant communities [20,21]. At the same time, [22] pointed out that precipitation and nitrogen deposition directly affect the composition of microbial communities. Both soil chemistry and soil microbes were affected by N input and increased precipitation, especially when they were applied simultaneously [23,24]. For example, some studies show that changes in N addition and precipitation alter microbial communities through changes in soil pH [25]. The addition of nitrogen mainly affects the soil microbial community by lowering the pH of the soil, while increasing water promotes the increase of soil pH [26]. Similarly, precipitation affects microbial communities by increasing plant diversity, while the addition of nitrogen does the opposite [24]. In addition, N addition significantly decreases the relative abundance of soil fungi and increases the proportion of bacteria, whereas water addition has the opposite effects [22,27,28]. These results indicate that the response of the soil microbial community to N availability is highly dependent on changes in precipitation [29]. There are many studies on global change factors affecting the structure and function of microbial communities, some of them report on short term environmental changes from a few months to two years [30,31], and others report on the impact of singular climate factors on microbial communities [32]. However, there are few studies on microbial community responses to long-term precipitation changes and N deposition, especially in desert steppes [29]. Therefore, it is highly valuable to study the effects of long-term changes in water and N deposition on microbial communities in desert steppe. Moreover, the interaction between N deposition and precipitation change has practical significance in the responses of microbial ecosystems to complex global changes in multiple environmental factors.
Soil microbes do not typically live in isolation, but instead form complex inter-species networks that substantially regulate ecological community structure [33] and ecosystem function [34]. Ecological network analysis, a system analysis method based on random matrix theory, analyzes the interactions between different entities in a system [35]. It has recently been used to study complex microbial systems, such as ecosystem food web [36] and microbial community structure [34,[37][38][39]. Analyzing the ecological network structure can reveal the relationship between ecosystem complexity and stability [40]. Complex interrelationships between species (such as predation, symbiosis, and competition) are important to the stability of the communities. Therefore, the complexity and stability of the ecological network established through the interrelationships of species can reflect the ability of the ecosystem to respond to changes in the external environment.
In this study, we conducted an 8-year field experiment to study the effects of interaction between N deposition and precipitation change on soil bacterial communities in a desert steppe using high-throughput sequencing technology. network analysis was an important method to examine the response of soil bacteria to changes in water and N. The objectives of this study were the followings: 1) Investigate whether soil bacterial communities are sensitive to long-term precipitation changes and N deposition, 2) explore the key species in the bacterial community affected by environmental factors, and 3) infer the impact of changes in precipitation and increased N deposition on desert steppe ecosystem.

Study site
The experiment was conducted from June 24, 2006 to October 11, 2014 in a natural desert steppe ecosystem in Siziwang Banner (41 46 0 43.6"N, 111 53 0 41.7"E), Inner Mongolia, northern China. From 2006 to 2014, the mean soil moisture content was 5.42%, with the highest monthly soil moisture from June to September. The mean soil temperature was 6.44˚C, with the highest monthly mean temperature in July (24.91˚C). The average annual rainfall was 205.29 mm, and more than 75% of measured rainfall occurred from June to September (S1 Fig). The soil in the study site is Kastanozem (Haplic Calcisols according to the FAO classification) with a sandy loam texture., the basic soil properties at 0-30 cm depth, according to our background survey, are shown in S1 Table. The dominant plant species are Stipa breviflora, Artemisia frigida and Cleistogenes Keng.

Experimental design
Thirty-six 6 m × 15 m experimental plots were established in June 2006 (Fig 1). The following six treatments were randomly assigned: control (CK), water addition (WA), water reduction (WR), N addition (NA), simultaneous N addition and water addition (NAWA) and simultaneous N addition and water reduction treatments (NAWR). The distance between plots was approximately 2 m. A plastic baffle with a depth of one meter was inserted between each group of plots to prevent the exchange of moisture and nutrients between differently treated adjacent plots. An artificial drought shelter was used for water translocation. Twelve 6 m × 7.5 m artificial drought sheds were randomly constructed in the main plots, with a rain shelter of 100% transmittance on top of the shed. Each artificial drought shed eliminated 30% of natural precipitation in the WR treatment. The rainwater from the WR subplot is evenly sprayed to the WA subplot to supplement 30% of natural precipitation after each rainfall. In the NAWA and NAWR subplots, besides water regulation, N was applied once a year using granular NH 4 NO 3 (10 g N m -2 ) [41] before rain in late April or early May. Each experimental treatment had six replicate plots. The experiment lasted 8 years, from 2006 to 2014. Rainfall and air temperature were obtained from meteorological stations established in the experimental site. A diagrammatic representation of the experimental device is shown in Fig 1. All necessary permits were obtained for the described study, which complied with all relevant regulations.

Soil sampling
Soil samples were collected in the middle of August in 2014, because the standing crop of the steppe communities reaches its annual peak from the middle to end of August. To obtain soil samples, 10 soil cores (30 cm depth and 5 cm diameter) were obtained from random locations in each plot, and then mixed to form one composite sample. Each soil core was divided into the following five layers: 0-2, 2-5, 5-10, 10-20 and 20-30 cm soil depths. After removing the roots, gravel and coarse fragments (2mm sterilized nylon mesh), soil samples were subdivided into two subsamples and stored at 4˚C and -80˚C for soil chemical analysis and DNA extraction, respectively.

Analysis of soil chemical parameters
The soil pH was measured using a pH electrode at a soil:water ratio of 1:2.5 (w/v) (Delta 320; Mettler-Toledo Instruments Co., Shanghai, China). Soil samples were dried in an oven at 105˚C for 12 h and weighed to determine the soil moisture content. Soil nitrate (NO 3 -) and ammonium (NH 4 + ) were extracted by shaking 5.0 g of fresh soil with 2 M KCl for 1 h, and their content was determined using a continuous flowing analyzer (SAN++, Skalar, Holland). The total C and total N concentrations of the soil were determined by dry combustion of duplicate subsamples using the LECO 2000 CHN Analyzer (LECO, Chicago, USA).

DNA extraction and high throughput sequencing
Total metagenomic DNA was isolated from soil samples using the Fast DNA SPIN Kit for Soil (MP Biotechnology, USA) according to manufacturer's instructions. The quality of DNA was evaluated by agarose gel electrophoresis. The V4 region of 16S rRNA was amplified using specific primers (515F, 5'-GTGCCAGCMGCCGCG GTAA-3' and 806R, 5'-GGACTACHVGGGTWTC TAAT-3') with designated barcodes [42]. PCR reactions were carried out in 30 μL reactions with 15 μL of Phusion High-Fidelity PCR Master Mix (New England Biolabs); 1μl forward primer (10μM), 1μl reverse primer (10μM), and 10 ng template DNA. Thermal cycling consisted of initial denaturation at 98˚C for 1 min, followed by 30 cycles of denaturation at 98˚C for 10 s, annealing at 50˚C for 30 s, elongation at 72˚C for 30 s, and final 72˚C for 5 min. Use TruSeq DNA PCR-Free Sample Preparation Kit for library construction. Libraries were assessed by the Qubit 2.0 Fluorometer (Thermo Scientific) and deep sequenced on the Miseq platform with PE300 strategy. The date has been uploaded to the National Microbiology Date Centers (NMDC10017737).

Bioinformatic analysis of 16S rRNA sequences
The raw readings are demultiplexed and assigned to the samples based on their unique barcodes. Trimmomatic (V0.35) was used for quality filtering with default parameters [43]. Then, high-quality paired-end reads were merged by using FLASH [44]. After removal of chimeric sequences based on the Silva database [45] and UChime algorithm [46], clean reads were subject to further analysis. Operational taxonomical units (OTUs) were clustered at a similarity threshold of 97% by using Uparse (v7.0.1001) [47]. All samples had sequencing depth of over 10 000 reads. Samples were normalized to 4269 sequences per sample. The Silva SSUrRNA (SSU128) Database was used based on the Mothur algorithm (thresholds: 0.8) to annotate taxonomic information for each representative sequence.

Network analysis
Network analysis is a method to determine the role of nodes in complex networks based on topological properties, and it offers new insights to determine key species and significant module members in microbe communities [48]. Network analysis was used to explore the cooccurrence patterns among the bacterial taxa. Analysis of ecological networks was performed using the Network Analyses Pipeline (http://ieg4.rccc.ou.edu/mena). More information on the pipeline, including its properties and relevant theories, can be found in [49][50][51]. Six co-occurrence networks of bacteria for the six treatments were built. For each network analysis, we used the data of samples from the five soil layers. The analysis was performed as follows. First, individual OTUs assigned into phylum level, was submitted to construct a network using the default settings. (OTUs appeared in more than half of all samples to be retained for subsequent analysis.) A similarity matrix was constructed based on Spearman's Rho between pairwise OTUs. A cutoff value (similarity threshold, St = 0.6) for the similarity matrix was generated automatically using the default settings. A link was assigned between a pair of OTUs if the correlation between their abundance was larger than the St. Second, a set of measures (the number of nodes and links, average path length, number of positive and negative links, average connectivity, average clustering coefficient and modularity) was calculated to describe the topology of the resulting networks. The different roles of each node in the network were identified using the Zi plots (of values that measure within-module connectivity) and Pi plots (of values that measure among-module connectivity). These parameters can be used to develop Zi-Pi diagrams, which can be used to study key species in a community [49][50][51]. Third, Use the "igraph" and "psych" packages in the R (3.5.2) environment to calculate and draw network plots. Finally, the "randomize the network structure and then calculate network" was run. Random networks were generated for comparison with the topology of the empirical network; each link had the same probability of being assigned to any node [52].

Statistical analyses
The Shannon-Wiener diversity index is used to evaluate the diversity of different treatment methods [53]. Using "vegan" packaging in the R environment, redundant analysis is used to quantify the relative contribution of soil chemical factors to the composition of the bacterial community. The Monte Carlo displacement test was used to test the significance of soil chemical factors related to changes in bacterial communities. One-way analysis of variance was used to test the significance of the chemical properties of each soil layer under different treatment conditions.
The result of bacteria is expressed as the value of six replicate samples mixed into one sample (no replicate) by high-throughput sequencing. The results of soil physical and chemical properties are standardized average values of six replicate plots. Moreover, the statistical significance was accepted at P < 0.05. The software Sigmaplot 12.5 and SPSS were used for the statistical analyses.

Ethics statements
All experimental procedures in our research have been approved by the College of Grassland, Resources and Environment, Inner Mongolia Agricultural University. The sampling site was the temperature and nitrogen increase experiment plot in Siziwang Banner, College of Grassland, Resources and Environment, Inner Mongolia Agricultural University.

Relationships between soil chemical properties and bacteria communities
Changes in the soil chemical properties among all the treatments in the five layers are listed in Table 1. Precipitation enrichment enhanced soil moisture in all the layers, especially in the 5-10, 10-20 and 20-30 cm layers, which reached a significant level (P < 0.05). Precipitation reduction decreased soil moisture in all the layers, but not significantly. However, precipitation enrichment and reduction had no significant effect on other chemical properties, including the soil pH, content of total C, TN, NH 4 + -N and NO 3 -N. At the same time, the changes of the total soil C and N concentrations were not significant under different treatments.
Redundancy analysis was performed to identify the relationship between environmental variables using soil physical and chemical properties and the soil bacterial community structures (Fig 2). The first two dimensions explained 50.1% and 28.35% of the total variances in the bacterial communities, respectively.

Diversity and composition of soil bacterial communities
We obtained 287 286 high-quality soil bacterial 16S rRNA sequences. Bacterial sequences were affiliated with 30 phyla. Because the top 10 most abundant phyla contain the most community information (accounting for more than 95% of the community), this was used by us to describe the main components of soil bacterial communities. Shannon index and genera richness were used to evaluate the diversity of the soil bacterial community in desert steppe (Fig 3). Overall, compared with CK treatment, the soil bacterial community diversity of WA treatment increased, and the soil bacterial community of NA treatment decreased. The soil bacterial community diversity of WR treatment and CK treatment was the most similar, and there was no difference. Under the NAWA treatment, the diversity of the soil bacterial community increased compared with that under the CK. However, compared with that under the WA treatment, the bacterial community biodiversity under the NAWA treatment decreased. Moreover, the diversity of soil bacterial communities under the NAWR treatment did not decrease compared with that under the CK. That is, there was no superimposed effect of N deposition and drought on the diversity of soil bacterial communities in desert steppe.
As shown in Fig 4, the differences in the relative abundance of all bacterial phyla of each treatment at the 0-2 cm soil layer was greater than that in other soil layers. In the 2-5 cm soil layer, the relative abundance of Acidobacteria was lower in all treatments compared to control. Especially, compared to other treatments, WR treatment caused highest relative decrease (26.62%). Compared with that under the CK, the relative abundance of Proteobacteria increased by 10.98% (WA), 14.96% (WR), 9.91% (NA), 14.33% (NAWA) and 11.15% (NAWR). The relative abundance of Cyanobacteria increased by 16.82% under the WR treatment compared with that under CK. We detected no differences in the relative abundance of any other bacterial phyla between the treatments.
In the 2-5 cm soil layer, the change in soil bacterial phyla was similar to that in the 0-2 cm soil layer. For example, the relative abundance of Proteobacteria increased by 14.78% under the WA treatment, but that of Acidobacteria decreased by 7.20% under the WR treatment, compared with that under the CK. In the 5-30 cm soil layer, the relative abundance of soil bacterial phyla was relatively stable except for that of Actinobacteria. Cyanobacteria appeared mainly in the 0-2 cm soil layer, and its relative abundance decreased gradually with increasing soil depth. In non-surface soil, the relative abundance of Cyanobacteria was higher than other treatments under NA and NAWA. Overall, Actinobacteria, Proteobacteria and Chloroflexi were the dominant bacterial phyla in the 0-30 cm soil layer, accounting for 47.31%, 12.87% and 6.97% of all bacterial phyla, respectively. The relative abundance of Actinobacteria was stable among different treatments. Compared with that under the CK, the relative abundance of Proteobacteria increased by 8.37% (WA) and 7.31% (NAWA).

Network analysis of soil bacterial communities
The visualized networks of soil bacterial communities showed the roles of different phyla and their interactions with the other phyla under different treatments (Fig 5). There were obvious differences among the six empirical networks of the bacterial communities under the CK, NA, WA, WR, NAWA and NAWR treatments and their corresponding random networks in terms of the average path distance, average clustering coefficient, and average modularity ( Table 2). Compared with the CK treatment, the topological structure of the soil bacterial network under different treatments showed greater differences, corresponding to differences in the major characteristics of bacterial communities (Fig 5, Table 2). The total number of links in the WA network increased by 11.76% ( Table 2) compared with that in the CK network, whereas there were no obvious changes in the average connectivity. The total number of links and average connectivity in the NAWA treatment network increased by 7.19% and 1.33 compared with those in the CK network, respectively. In the WR treatment network, these values decreased by 9.15% and 3.76, respectively, and in the NAWR treatment network, these values decreased by 14.05% and 3.95, respectively. The percentage of negative links in the WA, NAWA and NAWR networks decreased, whereas that in the WR network increased. The average path distance of WA and NAWA networks decreased, while the average path distance of WR and NAWR networks increased. Compared with the CK network, the modularity of these five networks have been improved.
The Zi-Pi plot in Fig 6 showed the distribution of bacteria based on their topological roles in soil bacterial networks. As shown in Fig 6, the majority of the nodes (> 95%) under different treatment networks were categorized as peripherals (specialists, Zi < 2.5 and Pi < 0.62, indicating nodes with only a few links that are predominant to other nodes within their modules). 21 nodes were classified as connectors (generalists, Zi < 2.5 and Pi � 0.62, which are predominantly linked to several modules). Fifteen of these 21 nodes belonged to Planctomycetes (3 OTU), Proteobacteria (2 OTU), Bacteroidetes (2 OTU), Chloroflexi (4 OTU),

PLOS ONE
The response of desert steppe to environment Acidobacteria (2 OTU), Actinobacteria (1 OTU) and Cyanobacteria (1 OTU) from the WR treatment network; two belonged to Proteobacteria (1 OTU) and Verrucomicrobia (1 OTU) from the WA treatment network; three belonged to Chlorobi (1 OTU), Chloroflexi (1 OTU) and Cyanobacteria (1 OTU) from the NA treatment network and one belonged to Bacteroidetes (1 OTU) from the CK treatment network. Our results suggest that there was no node in the network hub (super generalists, Zi � 2.5 and Pi � 0.62) category that existed as both a module hub and connector. As expected, more connectors (generalists) existed in the water regulation networks than in the CK network.

Discussion
In the 10-year experiment, we studied the effects of precipitation changes and N deposition on the soil bacterial community in the desert steppe environment. The results showed that precipitation changes and N deposition not only potentially affect the composition and diversity of soil bacterial communities, but also change their ecological network, and thus affect the

Effects of precipitation changes on the soil bacterial community
Our study area is located in arid and semi-arid steppes where the water content is a major limiting factor of steppe productivity. We expected that drought simulation for 8 consecutive years would significantly reduce the diversity of soil microbial communities, but our results were contrary to our expectations. Our results indicated that the diversity of the soil bacterial community was sensitive to increased precipitation, but they showed a strong tolerance to drought. This might be because soil bacteria inhabiting the arid environment long term have a great potential for drought tolerance. The distribution of precipitation throughout the year in our experimental area was uneven. More than 75% of the measured rainfall occurred from June to September (S1 Fig). Long-term and frequent drying-wetting cycle stress might be one of the reasons that the microbial community has great potential for drought tolerance [32]. , both N addition and water addition (NAWA), and both N addition and water reduction (NAWR). Network plots were drawn using using the "igraph" and "psych" package in the R environment. The node colors represent bacterial phyla of the operational taxonomic units (OTUs). The node size is proportional to the square root of the abundance of the OTU. The red edges indicate positive relationships and green edges indicate negative relationships. The thickness of the connection between two phyla represents the size of the correlation coefficient. (All depths of the same treatment are included). https://doi.org/10.1371/journal.pone.0248194.g005 Compared to that in other soils, [32] pointed out that microbes had a greater ability to maintain their function in response to repeated water stress in grassland soil. In addition, some research shows that a change in precipitation has no significant effect on the composition of the bacterial community of the top soil or rhizosphere soil, but does affect the abundance of Note: CK represents the control, WA represents the water addition treatment, WR represents the water reduction treatment, NA represents the N addition, NAWA represents the N addition and water addition treatment, and NAWR represents the N addition and water reduction treatment.
https://doi.org/10.1371/journal.pone.0248194.t002 bacteria [23,54]. This is consistent with our results. Precipitation changes did not significantly affect the composition of the soil bacterial community at the phylum level, but only affected a few species sensitive to soil moisture. Precipitation change also affected the stability of the soil bacterial community in the desert steppe. The co-occurrence of bacteria within the community leads to the complexity of the community structure and maintains the stability of the community [30,37,55]. In terms of community stability and sensitivity, increased modularity reduced the influence of environmental factors on the entire network and enhanced the stability of the soil bacterial community. The shorter average path distance indicated higher sensitivity and more rapid response to water addition. Increased modularity suggested that water reduction reduced the influence of environmental factors on the entire network and enhanced the stability of the soil bacterial community. The longer average path distance indicated a lower sensitivity and slower response to water reduction. This might be related to the adaptation of bacterial communities to the drought conditions of desert steppe.

Effects of water on the soil bacterial communities with N enrichment
Water and N addition increased the influence of environmental factors on the entire network and weakened the stability of the soil bacterial community, whereas N addition and water reduction had opposite effects on the stability and sensitivity of the soil bacterial community. Actually, the effect of water and N on soil microbial communities in desert steppe has been reported previously. Currently, there are two opinions about the effect of water and N addition on the soil bacterial communities in desert steppe. First, water and N addition might interactively affect soil microbial communities in semiarid steppe [56]. Second, water and N addition can have opposite effects on soil bacterial communities [13,27,57]. Our study indicated that N addition weakened the effect of water addition in terms of bacterial diversity and community stability, and they did not generate interactive influence. Moreover, N addition and water reduction together did not affect the diversity and stability of soil bacterial communities; that is, N deposition and drought did not have a cumulative effect on soil bacterial communities in desert steppe. This complexity might be related to the efficiency of water and N use in soil. Although precipitation enrichment could compensate the N requirement for microbe growth in arid soil by promoting N mineralization, increasing water availability induces substrate loss through nitrate leaching [58]. It has been demonstrated that precipitation conditions could lead to nutrients loss, thereby inducing soil nutrient deficiency [59]. These results indicated that precipitation plays an essential role in nutrient budget, which in turn influences microbial functional groups. In addition, precipitation regimes could alleviate the negative effects of soil pH induced by N addition on soil microbial communities [60]. Thus, the absence of interactive effects confirmed that the response of the soil microbial community to N addition could be mediated by precipitation management in desert steppe.
Moreover, N addition has a negligible effect on soil bacterial communities, which can also be explained by the limitation in soil C. As reported previously, N addition can also lead to a decrease in soil microbial community biomass and activity in systems that are C limited [61]. Total soil N and C concentrations were not significantly influenced by N addition in our field experiment site. A previous study has shown that N addition caused an increase in aboveground net primary productivity, but root biomass did not proportionally increase with aboveground biomass in desert steppe [62]. This can explain why the total soil C concentration did not increase consistently in plots with N addition. Heterotrophic microbes are thought to be C limited when the C:N ratio is below 30, and N limited when the C:N ratio is above 30 [63]. The C:N ratio in our study ranged from 7.35 to 9.98. At the same time, the results of redundancy analysis (RDA) (Fig 2B) also showed that total organic carbon was significantly correlated with the bacterial community structure (P < 0.05). These results suggest that soil microbes in our study were C limited, which might greatly influence the soil microbe community in desert steppe.

Responses of the keystone bacteria to precipitation changes and N deposition
Generalists are considered to be key organisms in a community [50]. In this study, 21 nodes from the WR, WA and CK treatments were categorized as connectors (generalists), which were highly connected to several modules. The 21 generalists represented nine phyla, viz. Planctomycetes, Proteobacteria, Bacteroidetes, Chloroflexi, Acidobacteria, Actinobacteria, Chlorobi, Cyanobacteria and Verrucomicrobia, which can be the key phyla affecting the stability of soil microbial communities in desert steppe.
At the same time, the difference in the relative abundance of bacteria in the community composition is also an indicator to determine the keystone bacteria. We found that there are seven bacterial phyla that account for more than 90% of the total community, namely Actinobacteria, Proteobacteria, Cyanobacteria, Chloroflexi, Acidobacteria, Planctomycetes and Gemmatimonadetes. Among them, the relative abundance of Actinobacteria, Proteobacteria, Acidobacteria and Cyanobacteria greatly differed with precipitation changes and N deposition treatments. In addition, the intersection of the key phyla that affect the ecological network and community composition may be the most important keystone bacteria responding to changes in precipitation and N deposition in desert desert steppe. In this study, we found that these four phyla are the keystone phyla of the ecological network, and the relative abundance of their bacterial communities will change with precipitation and N deposition treatments. Therefore, Actinobacteria, Proteobacteria, Acidobacteria and Cyanobacteria are the key phyla affecting the composition of soil microbial communities in desert steppe.
In fact, there are still many deficiencies in our research on network analysis. The effects of water-nitrogen on the stability and sensitivity of the ecological network can help us understand the overall response of soil bacterial communities to environmental changes. At the same time, keystone species that respond significantly to environmental changes can be identified. Although we can identify keystone phyla by the method of network analysis, we still cannot recognize sensitive species and conservative species in keystone phyla. We still think it is a practical approach to narrow down the range of key species. The key species we identified need to be further validated by systematic experiments.

Conclusions
The response of the composition, diversity, and network of soil bacterial communities to precipitation change and N deposition exhibited a certain degree of uniqueness in desert steppe. Decreasing precipitation and increasing N deposition did not exhibit a cumulative effect on soil bacterial communities in the desert steppe. This means that even if the diversity and functions of the ecosystem are expected to be affected by N deposition in the future, combination with precipitation may help maintain the stability of the semi-arid steppe ecosystem under climate change. Based on network analysis and relative abundance, we identified Actinobacteria, Proteobacteria, Acidobacteria and Cyanobacteria as the most important keystone bacteria responding to precipitation changes and N deposition in the soil of desert steppe. In future research, we should focus on the underlying mechanisms of these key phyla to climate change and find the key species which are sensitive and conservative to global change in the keystone phyla.