Next Article in Journal
Interaction of Naturally Occurring Phytoplankton with the Biogeochemical Cycling of Mercury in Aquatic Environments and Its Effects on Global Hg Pollution and Public Health
Previous Article in Journal
Diagnosis, Characterization and Treatment of Emerging Pathogens
Previous Article in Special Issue
Escherichia coli Strains Display Varying Susceptibility to Grazing by the Soil Amoeba Dictyostelium discoideum
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High Ammonium Addition Changes the Diversity and Structure of Bacterial Communities in Temperate Wetland Soils of Northeastern China

1
Engineering Research Center of Agricultural Microbiology Technology, Ministry of Education & Heilongjiang Provincial Key Laboratory of Ecological Restoration and Resource Utilization for Cold Region & Key Laboratory of Microbiology, College of Heilongjiang Province & School of Life Sciences, Heilongjiang University, Harbin 150080, China
2
Swiss Federal Institute for Forest, Snow and Landscape Research WSL, 8903 Birmensdorf, Switzerland
3
Institute of Nature and Ecology, Heilongjiang Academy of Sciences, Harbin 150001, China
4
Heilongjiang Academy of Forestry, Harbin 150022, China
5
Key Laboratory of Geographical Processes and Ecological Security in Changbai Mountains, Ministry of Education, School of Geographical Sciences, Northeast Normal University, Changchun 130024, China
6
School of Life Science, Hebei University, Baoding 071002, China
*
Authors to whom correspondence should be addressed.
Microorganisms 2023, 11(8), 2033; https://doi.org/10.3390/microorganisms11082033
Submission received: 4 July 2023 / Revised: 5 August 2023 / Accepted: 6 August 2023 / Published: 8 August 2023
(This article belongs to the Special Issue Microbial Interactions in Soil 2.0)

Abstract

:
The soil microbiome is an important component of wetland ecosystems and plays a pivotal role in nutrient cycling and climate regulation. Nitrogen (N) addition influences the soil’s microbial diversity, composition, and function by affecting the soil’s nutrient status. The change in soil bacterial diversity and composition in temperate wetland ecosystems in response to high ammonium nitrogen additions remains unclear. In this study, we used high-throughput sequencing technology to study the changes of soil bacterial diversity and community structure with increasing ammonium concentrations [CK (control, 0 kg ha−1 a−1), LN (low nitrogen addition, 40 kg ha−1 a−1), and HN (high nitrogen addition, 80 kg ha−1 a−1)] at a field experimental site in the Sanjiang Plain wetland, China. Our results showed that except for soil organic carbon (SOC), other soil physicochemical parameters, i.e., soil moisture content (SMC), dissolved organic nitrogen (DON), total nitrogen (TN), pH, ammonium nitrogen (NH4+), and dissolved organic carbon (DOC), changed significantly among three ammonium nitrogen addition concentrations (p < 0.05). Compared to CK, LN did not change soil bacterial α-diversity (p > 0.05), and HN only decreased the Shannon (p < 0.05) and did not change the Chao (p > 0.05) indices of soil bacterial community. Ammonium nitrogen addition did not significantly affect the soil’s bacterial community structure based on non-metric multidimensional scaling (NMDS) and PERMANOVA (ADONIS) analyses. Acidobacteriota (24.96–31.11%), Proteobacteria (16.82–26.78%), Chloroflexi (10.34–18.09%), Verrucomicrobiota (5.23–11.56%), and Actinobacteriota (5.63–8.75%) were the most abundant bacterial phyla in the soils. Nitrogen addition changed the complexity and stability of the bacterial network. SMC, NO3, and pH were the main drivers of the bacterial community structure. These findings indicate that enhanced atmospheric nitrogen addition may have an impact on bacterial communities in soil, and this study will allow us to better understand the response of the soil microbiome in wetland ecosystems in the framework of increasing nitrogen deposition.

1. Introduction

In recent decades, agricultural fertilization, combustion of fossil fuels, and other human activities have increased nitrogen inputs to wetland ecosystems [1]. Therefore, external nitrogen addition and atmospheric nitrogen deposition can directly or indirectly affect soil’s aboveground biological processes and underground biochemistry [2,3]. Nitrogen addition can influence soil’s microbial communities by alleviating resource limitations [4] and altering its physiochemical properties [5]. Conversely, as an important part of regulating the underground ecological process, the dynamic changes in soil microorganisms can also provide feedback to ecosystem functions [6] by affecting litter decomposition [7,8], nutrient cycling, and greenhouse gas emissions [9]. Nevertheless, despite the critical contributions of soil microorganisms to the underground ecological process and ecosystem stability, the studies that characterized the effects of nitrogen addition on soil microbial community structure and diversity were still limited.
Soil bacteria are a major component of soil microorganisms and play crucial functions in the ecosystem [10,11]. Moreover, their diversity, composition, complexity, and stability are particularly sensitive to changes in nitrogen content [12]. Several studies have shown that nitrogen addition can lead to major shifts in various edaphic parameters, which affect soil bacterial communities. For example, N addition indirectly affected the bacterial community composition by decreasing soil pH [3], and Song et al. [8] found that increased N availability enhanced bacterial activity. Li et al. [13] and Xie et al. [14] found that increasing N input significantly affected the soil’s microbial species interactions and stability, which in turn affects the stability of the microbial community structure [13,15]. Moreover, soil N availability also played an important role in shaping microbial networks [16], such as ammonium nitrogen (NH4+-N) concentration enhanced inter-microbial interactions [17]. Some studies indicated that nitrogen addition significantly changed the relative abundance of bacteria [18,19]. For example, nitrogen addition enhanced the relative abundance of Chloroflexi, and reduced the relative abundance of Chlorobi and Verrucomicrobia in a coastal wetland [20]. Zhang et al. [21] also found that N addition increased the relative abundance of Chloroflexi but decreased the relative abundance of Firmicutes and Bacteroidetes in Deyeuxia angustifolia wetland. In contrast, Hu et al. [2] reported that high nitrogen treatment significantly reduced the relative abundance of Chloroflexi in alpine wetland. The varied results may be attributed to the different forms of nitrogen application because the above studies usually used urea or ammonium nitrate fertilizer (NH4NO3) to simulate the nitrogen deposition. However, the responses of bacterial communities to ammonium nitrogen (NH4+) inputs remain lacking in studies.
Ammonium nitrogen is the dominant form of atmospheric nitrogen deposition [22]. Ammonium nitrogen may increase aboveground biomass [23], stimulate soil bacterial activities [24], increase the ratio of NH4+ to NO3, and affect the rate of soil N transformation [25]. In addition, ammonium nitrogen can produce more protons than nitrate nitrogen in the soil [26], and bacteria preferentially assimilate NH4+ over NO3 due to the lower energy expenditure for NH4+ assimilation [27]. At present, many studies have been conducted on the effects of ammonium nitrogen addition to soil bacteria. For example, Zhang et al. [28] evaluated the biomass of soil microbial communities responding to ammonium (NH4Cl) application in a slash pine (Pinus elliottii) plantation in a subtropical forest ecosystem, and found that ammonium treatments changed the bacterial content of phospholipid fatty acids (PLFAs). Liu et al. [29] examined soil bacteria in Chinese fir plantations using plate counts and PLFA analysis and found that the addition of ammonium (NH4+) increased bacterial abundance and biomass. Geng et al. [30] found that ammonium addition changed bacterial community structure and increased the ratio of fungi/bacteria (F/B) in a boreal forest. However, these studies were mainly carried out in forest ecosystems, and the response of bacteria to ammonium nitrogen in wetland ecosystems remains unclear.
Sanjiang Plain is a typical wetland, representative of the global temperate wetland ecosystem [31]. This freshwater marsh wetland is the largest concentrated distribution in the northeast of China and is one of the most sensitive areas to global changes [32]. D. angutifolia wetland is the main wetland type in the Sanjiang Plain, accounting for 34.45% of Sanjiang Plain wetland [33]. The nitrogen deposition concentrations of D. angutifolia wetland have increased severely, and this is attributed to agricultural and industrial activities [34,35]. These activities’ ecological impacts on soil microbial communities have been of great concern in recent years. For example, Song et al. [8] found that N addition increased litter decomposition rates through enhanced microbial activity. Zhang et al. [21] showed that N addition changed the structure and abundance of bacterial communities. Sui et al. [36] indicated that N addition reduced the abundance of some major bacterial groups, but that there was no significant impact on the bacterial community structure. The above studies mainly used NH4NO3 as a nitrogen source, however, the NH4+ form was one of the primary forms of nitrogen deposition in the Sanjiang Plain [34] and, to our knowledge, there are no studies on the response of the bacterial community to NH4+ addition in Sanjiang Plain. Therefore, our study investigated the changes in the soil bacterial community, diversity, and stability using high-throughput sequencing technology and the correlation with the soil’s physicochemical parameters in response to NH4+ addition on the experimental platform of the Sanjiang Plain. We propose the following hypotheses: (1) NH4+ addition would change the soil’s bacterial structure, changing the complexity and stability of the network structure of the bacterial community. (2) High nitrogen would reduce bacterial diversity more strongly than low nitrogen. This research will be of great significance for comprehensively evaluating wetland nitrogen cycling and improving sustainable wetland management in Northeastern China.

2. Materials and Methods

2.1. Study Area

This study area is located in the Honghe National Nature Reserve of Sanjiang Plain (47°42′18″–47°52′07″ N, 133°34′38″–133°46′29″ E) (Supplementary Figure S1). The annual average temperature and evaporation are 1.9 °C and 1166 mm, respectively. The annual precipitation is 585 mm, mostly concentrated from July to September. Northwest winds prevail throughout the year, and the ≥10 °C effective accumulated temperature is 2165~2624 °C, with an annual frost-free period of 131 day. The soil belongs to the categories of bleached stagnant soil and fibrous organic soil [37]. The main wetland types are meadow wetland, meadow swamp wetland, and swamp wetland. The dominant vegetations are D. angustifolia, Glyceria spiculose, Carex lasiocarpa, and Carex pseudocuraica.

2.2. Soil Sampling

The field experiment was established in May 2016 with three concentration levels: control (CK: 0 kg NH4+ ha−1 year), low nitrogen (LN: 40 kg NH4+ ha−1 year), and high nitrogen (HN: 80 kg NH4+ ha−1 year). The low nitrogen addition concentrations (40 kg N ha−1 year) were set up based on the local fertilizer application and the current status of atmospheric N deposition in China [21]. The high nitrogen addition concentrations (80 kg N ha−1 year) were set up to predict the impact of continued N input on soil microorganisms in the future. Each treatment had three replicates (9 plots in total) and each plot was 10 m × 10 m in area. NH4Cl was used as the N source, dissolved in water, and sprayed uniformly in May each year, and CK was sprayed with an equal amount of water. In October 2021 (five and half years after the beginning of the experiment), 15 to 20 soil samples (0~20 cm) were randomly obtained using a soil auger (5 cm in diameter, 20 cm deep) after removing the litter layer from each plot. The soil samples were mixed into one soil sample to represent the soil sample of the plot, after removing impurities such as vegetation fallout from the sample, and then placed in a carrying incubator at 4 °C and quickly transported to the laboratory. The soil samples were divided into two parts. One was stored at −80 °C for the determination of the soil’s bacterial community, and another was air-dried, thoroughly ground, and sieved by 2 mm for the determination of its physicochemical properties.

2.3. Determination of Soil Physicochemical Properties

Soil moisture content (SMC) was determined via the drying method [33]; soil pH was measured by a water–soil ratio of 1:2.5; soil total nitrogen (TN) was measured by the Autoanalyzer (AutoAnalyzer 3, Hamburg, Germany); soil organic carbon (SOC) was determined by an elemental analyzer (VarioEL III, Frankfurt, Germany); nitrate nitrogen (NO3) and ammonium nitrogen (NH4+) were measured by the Autoanalyzer (AutoAnalyzer 3, Hamburg, Germany); dissolved organic carbon (DOC) and dissolved organic nitrogen (DON) were determined by a K2SO4 leaching method at a ratio of 1:5 soil to liquid, extracted for 30 min, centrifuged and filtered, and then measured by a TOC analyzer [38].

2.4. DNA Extraction

Soil bacterial DNA was extracted using the E.Z.N.A.® Soil DNA Kit (Omega Bio-tek, Norcross, GA, USA). PCR amplification was performed with the primer pairs 515F (5′-barcode-GTGCCAGCMGCCGCGG-3′) and 907R (5′-CCGTCAATTCMTTTRAGTTT-3′) spanning the variable region V4-V5 of the 16S rRNA gene [39], where the barcode is an eight-base sequence unique to each sample. PCR reactions were made in a total volume of 20 μL that contained 4 μL of 5 × FastPfu Buffer, 2 μL of 2.5 mM dNTPs, 0.8 μL of each primer (5 μM), 0.4 μL of FastPfu Polymerase, 10 ng of template DNA, and finally ddH2O up to 20 μL. PCR amplification conditions were as follows: 95 °C for 2 min, followed by 25 cycles of 95 °C for 30 s, 55 °C for 30 s, 72 °C for 30 s, and a final extension at 72 °C for 5 min. Three replicates of each sample were made and PCR products from the same sample were pooled and subjected to agarose electrophoresis and purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA). Purified PCR products were quantified using Qubit® 3.0 (Life Invitrogen, Waltham, USA) and mixed in the appropriate proportions as required for the sequencing volume of each sample. Sequencing was performed on an Illumina MiSeq platform (Shanghai BIOZERON Co., Ltd., Shanghai, China). All raw sequences have been deposited into an NCBI Sequence Read Archive with the BioProject accession: PRJNA880534.

2.5. Bioinformatics

Raw bacterial fastq files were de-multiplexed, quality filtered, and assessed through QIIME 2 (version 2022.2). Forward and reverse reads were merged using PEAR (Paired-End reAd mergeR, version 0.9.6) software. Low quality sequences less than 200 bp in length and sequences with an average quality score less than 20 were removed before further analysis. Exact barcode matching was implemented, allowing for a two-nucleotide mismatch in the primer matching. Reads including unclear characters were also included for removal. The trimmed sequences were checked for chimera, and those with chimeras were removed using the Uchime algorithm. Only sequences with an overlap of more than 10 bp were assembled based on their overlapping sequences; reads that were not able to be assembled were discarded. Sequences were clustered into operational taxonomic units (OTUs) at 97% similarity using UPARSE. With respect to taxonomy, the obtained OTUs were annotated using the SILVA database (version. 138). We rarefied the data set to the lowest number of sequences for all the samples and used this data set for further community analysis.

2.6. Statistical Analyses

One-way ANOVA and a Duncan test were used to analyze the effect of different levels of N addition on the soil’s physicochemical properties, and p < 0.05 was considered statistically significant via using the SPSS 20.0 (SPSS, Inc., Chicago, IL, USA). The rarefaction curve was created using R (version 4.4.3) software. Stamp analysis [40] was used to produce box plots of diversity indices and to analyze differences in bacterial community composition (phylum and genus level) between the different nitrogen addition treatments. A correlation matrix was constructed to look for relationships between soil properties and bacterial community by using the “Hmisc” package of R. The α diversity indices (Chao1, Shannon, and Simpson) were performed by Mothur (version 1.30.1). Redundancy analysis (RDA) and non-metric multidimensional scaling (NMDS) (using Bray–Curtis dissimilarity) and a Venn diagram were achieved using the vegan package in the R program (version 4.1.0). Indicator species analysis of soil bacterial composition under different N treatments were analyzed using the LEfSe component of the software Galaxy 21.09 (LDA = 3.5). For the network analysis, only OTUs with >1% abundance were selected and imported into the microplatform (http://www.cloud.biomicroclass.com/CloudPlatform/SoftPage/IGC) for intra-group correlation analysis (Pearson correlation) (accessed on 26 March 2023). Data files for edges and nodes were downloaded, and edge data with p < 0.05 were selected to import into Cytoscape for network plotting.

3. Result

3.1. Changes of Soil Physicochemical Characteristics under Different NH4+ Addition

Different nitrogen addition treatments had different effects on the soil’s physicochemical parameters (Table 1). The SMC, DON, and TN significantly increased (p < 0.05) and, inversely, pH significantly decreased (p < 0.05) with increasing concentrations of N addition. Compared with CK, HN treatment increased the contents of NO3 and DOC, whereas it decreased the NH4+, while LN treatment decreased the NH4+ and DOC (p < 0.05).

3.2. Sequencing Data Analysis and Distribution of Soil Bacterial Community

After high-throughput sequencing and optimization, a total of 429,316 sequences were obtained from the amplicons. A total of 178,818,450 bp was obtained, with a mean read length of 416.59 bp. The reads of 401–450 bp comprised 99.56% of the total sequences obtained and rarefaction curves were close to saturation, indicating that the sequencing data were representative of the real situation of the wetland soil bacterial community (Supplementary Figure S2).
In total, 7727 OTUs were detected, of which 3247 (42.02%) were shared between the CK, LN, and HN treatments (Supplementary Figure S3). CK samples produced the highest number of specific OTUs (912, accounting for 11.80%), while 743 OTUs were specific to LN and 675 OTUs were specific to HN, accounting for 9.62% and 8.74%, respectively. A large proportion of OTUs were shared in CK and LN (16.08%), while 5.99% of OTUs were shared between LN and HN, and 5.75% were isolated from both CK and HN. Therefore, the distribution of soil bacteria OTUs in the ecological network was significantly influenced by the nitrogen addition in the D. angustifolia wetlands of the Sanjiang Plain.

3.3. Increasing Ammonium Addition Influenced Bacterial α-Diversity

The α-diversities of the three nitrogen addition treatments are shown in Figure 1. Compared with CK, LN did not significantly change the Shannon index and the Simpson index (p > 0.05), but HN significantly decreased the Shannon index and increased the Simpson index (p < 0.05). Thus, the HN treatment significantly reduced the heterogeneity of the soil bacterial community. The Chao1 index and OTU numbers did not significantly differ (p > 0.05) between the CK, LN, and HN treatment, indicating that N addition did not change the richness of the soil bacterial community.
The effect of ammonium addition on bacterial community composition was demonstrated by non-metric multidimensional scaling (NMDS) based on Bray–Curtis dissimilarity (Figure 2). The soil bacterial community structure of LN, HN, and CK were significantly separated (Table S1, R2 = 0.48 and p < 0.05), indicating that the higher ammonium addition significantly affected the bacterial community structure. Interestingly, the PERMANOVA analysis showed that soil bacterial community structure of LN and HN did not significantly change compared to CK (Table S1).

3.4. Soil Bacterial Community Composition Was Changed by Ammonium Addition

A total of 7727 bacterial OTUs belonged to 44 phyla. The Acidobacteriota (24.96–31.11%), Proteobacteria (16.82–26.78%), Chloroflexi (10.34–18.09%), Verrucomicrobiota (5.23–11.56%), and Actinobacteriota (5.63–8.75%) were the most abundant bacterial phyla across all N treatments, and they represented approximately 80.85% of the reads from bacterial phyla (Figure 3). LN did not significantly change bacterial phyla (p > 0.05). HN treatment significantly increased the relative abundance of Acidobacteriota, but significantly decreased the relative abundance of Proteobacteria and Myxococcota (Supplementary Figure S5, p < 0.05).
At the genus level, the dominated bacterial genera were Candidatus Udaeobacter (11.66–39.45%), Candidatus Solibacter (14.32–16.51%), Bryobacter (5.66–6.92%), Bradyrhizobium (4.48–7.54%), Acidothermus (2.99–5.55%), Haliangium (2.13–4.91%), Candidatus Koribacter (2.55–4.44%), and Anaeromyxobacter (1.91–2.10%), and they represented approximately 64.28% of the reads from bacterial genera (Figure 3). Furthermore, compared with CK, the LN treatment significantly increased the relative abundances of Bradyrhizobium, Acidothermus, and Conexibacter and decreased the relative abundances of Haliangium (Supplementary Figure S6, p < 0.05). HN treatment significantly decreased the relative abundances of Gemmatimonas, Sphingomonas, Haliangium, Puia, Reyanella, and Citrifermentans (Supplementary Figure S7, p < 0.05).

3.5. Indicator Analysis of Soil Bacterial Community Responding to Ammonium Addition

The cladogram of the soil bacterial community under different N addition treatments with an LDA score of 3.5 and 23 taxa exhibited significant differences among the three N treatments (Figure 4). Each node in the graph represents a bacterial taxon, with non-significant bacteria labelled using yellow nodes and significant bacteria given other colors (red, blue, and green). The colors represent the subgroups in which they were abundant, i.e., they indicate that these branches had a higher abundance in that subgroup. As seen in the Figure 4, Minicystis of Myxococcota, Rhodococcus of Actinobacterita, Sphingomonas, Ralstonia, Phenylobacterium, and Rhizobacter of Proteobacteria were indicator genera in CK; Edaphobaculum was an indicator genus in LN; and Candidatus_Udaeobacter was an indicator genus in HN.

3.6. Ammonium Addition Changed the Bacterial Network Structure and Complexity

A network can be developed to represent the interaction and connectivity among different bacterial populations carrying the OTUs (Figure 5). The nodes of all the networks were assigned to 15 bacterial phyla. The OTU numbers, average geodesic distance, and number of modules under LN and HN were smaller than those of CK, but the average clustering coefficient increased from CK to HN (Table 2), indicating that LN and HN reduced the stability of the bacterial network. Total links (positive and negative links), average degree, and network density all showed a tendency as HN > CK > LN, indicating that HN increased the complexity of the bacterial network, whereas LN decreased the complexity compared to CK. These results indicated that N addition changed the complexity and stability of the bacterial network.

3.7. Relationships between Microbial Communities and Soil Properties

The effects of soil properties on bacterial communities were analyzed using RDA (Figure 6). SMC, pH, and NO3 were the three most influential factors contributing to the changes in bacterial communities (Table S3), with the first two RDA axes accounting for 46.05% of the variation in bacterial communities. The soil bacterial communities under both CK and LN treatments were positively correlated with NH4+ and pH, whereas they were inhibited by NO3, SMC, DOC, TN, DON, and SOC. The bacterial community under HN treatment showed positive correlations with NO3, SMC, DOC, TN, DON, and SOC, and negative correlations with the soil pH and NH4+.
The relative abundance of the top 10 phyla and top 20 genera of soil bacteria were correlated with environmental factors, and the results showed that almost all phyla and 14 genera were significantly correlated with environmental factors (Figure 7). At the phylum level (Figure 7a), the abundances of Acidobacteriota, Proteobacteriota, Myxococcota, Desulfobacterota, and Bacteroidota presented a positive relationship with pH and NH4+, while they showed a negative relationship with DOC, DON, SOC, NO3, SMC, and TN. Verrucomicrobiota, Methylomirabilota, Acidobacteriota, and Chloroflexi showed the opposite results. At the genus level (Figure 7b), 13 genera positively correlated with pH and NH4+, while they showed a negative relationship with DOC, DON, SOC, and NO3. In contrast, Candidatus_Udaeobacter showed a positive relationship with DOC, DON, SOC, and NO3.

4. Discussion

4.1. Changes in Soil Bacterial Community Diversity and Structure

In accordance with our second hypothesis, our results showed that HN decreased the Shannon diversity and community richness, but LN had no significant impact on the bacterial α diversity (Figure 1). Nitrogen input to soil may change the nitrogen’s ecological niche occupied by different species in the ecosystem [41], which in turn may cause changes in bacterial α diversity [42]. The difference between the Shannon index and the Simpson index in HN was related to species-relative abundance and species richness. The Simpson index was particularly sensitive to the abundance of dominant species, and HN significantly reduced the relative abundance of Gemmatimonas, Sphingomonas, Haliangium, Puia, Reyranella, and Citrifermentans (Figure S7), which may have caused an increase in the Simpson index. As for the Shannon index, which is closely related to the complexity of the communities, the reduction in species richness under HN would reduce the Shannon index (Figure 1 and Figure S3). N addition may increase N availability, changing the substrate utilization patterns in bacterial communities [43], reducing the utilization of carbon sources by the bacterial community, and decreasing the metabolic activity and diversity of the bacteria [29]. Our results were consistent with the findings of previous studies, which showed that increasing N addition concentrations reduced bacterial community diversity [18,44]. High nitrogen addition reduced soil pH, causing acidification of the soil, and therefore changing the diversity of the soil bacterial community [45]. However, previous studies are inconsistent with our findings [2,46]. For example, Hu et al. [2] conducted research on soil bacterial community in a swamp meadow wetland and found that short-term nitrogen application has no significant effect on the microbial diversity. This may be due to the duration of nitrogen addition. In our study, the 5-year N addition may have alleviated the limitation of nitrogen use by the soil bacterial community and caused a saturated state of N in the soil, therefore affecting the soil bacterial community. In addition, the differences in plant types on the ground may have caused differences in plant root exudates and ground litter, resulting in different soil environments, which in turn have different effects on the soil bacterial community [47,48].
According to the NMDS analysis, NH4+ addition changed the bacterial community structure (Figure 2, Supplementary Table S1), which is in accordance with our first hypothesis. The CK and LN were located in the negative half-axis of NMDS1, and HN was in the positive half-axis of NMDS1, indicating that HN had a greater effect on the soil bacterial community structure than LN. The reason may be that the continuous high NH4+ input altered the soil environment more than the low NH4+ concentration, and RDA analysis showed that SMC, NO3, and pH were the main factors affecting the bacterial community structure, which was sensitive to NH4+ concentration. Therefore, the continuous input of N changed the soil’s environmental factors and affected the bacterial community structure [2,15]. However, Sha et al. [49] investigated six years of nitrogen addition on soil microorganisms in a desert ecosystem (NH4NO3 as nitrogen source) and found that enhanced atmospheric nitrogen addition did not change the bacterial community structure. The differences with our results may be due to the inconsistency of the nitrogen sources applied in the study. Different forms of N input may have different effects on soil bacteria owing to preferential uptake by plants and competition between plant and microbial communities [23,50]. Furthermore, bacteria have different ways of using NH4+ and NO3; NH4+ can be directly used by bacteria, while bacteria need to reduce NO3 to NH4+ before assimilating it [50,51,52,53]. In addition, the wetland is formed by land–water interaction, and the wetland ecosystem is more vulnerable and sensitive to environmental changes when compared with the desert ecosystem [2]. Therefore, the response of wetland ecosystems to nitrogen inputs would be more significant.

4.2. Effects of N Addition on the Soil Bacterial Community Composition

Nitrogen addition significantly changed the composition of the soil’s bacterial phyla communities (p < 0.05) (Figure 3, Supplementary Figures S4–S7). In our study, we found that Acidobacteriota, Proteobacteria, and Chloroflexi were essential components of the soil bacteria occurring in the temperate wetland, which is consistent with those found in other wetland ecosystems [2,54]. According to oligotrophic–copiotrophic theory, fast-growing copiotrophic taxa are more likely to grow under nutrient-rich conditions. Conversely, oligotrophic taxa can survive under lower nutrient conditions, since they have a low growth rate [55]. Acidobacteria and Chloroflexi are regarded as oligotrophic microorganisms, while Proteobacteria are considered to be copiotrophic microorganisms [55,56]. However, in our study, we observed that nitrogen addition significantly increased the relative abundance of Acidobacteriota and Chloroflexi, whereas it decreased the Proteobacteria (p < 0.05) (Supplementary Figures S4 and S5). However, these results were contrary to the oligotrophic–copiotrophic theory results [57,58]. The reason for this result could possibly be explained by the changes in the soil’s physicochemical properties resulting from the continued ammonium nitrogen addition. The Acidobacteriota is generally acidophilic bacteria [56], and Acidobacteria is usually negatively correlated with soil pH (Figure 7; [59,60]). Therefore, ammonium nitrogen addition may be increasing the relative abundance of Acidobacteriota by reducing soil pH [45,61]. Soil bacterial communities have a narrow pH range for optimal growth, and acidification makes the environment more restrictive, which may also be a reason for the change in the relative abundance of the Proteobacteria [28,62]. Low pH would lead to a decrease in the relative abundance of Proteobacteria, which positively correlated with soil pH [63,64]. According to previous studies, the Chloroflexi were highly tolerant to eutrophication, and they contributed to nutrient removal [65,66]; nitrogen addition increased the nitrogen availability, which may significantly increase the relative abundance of Chloroflexi. Furthermore, the RDA showed that SMC was the most influential factor, resulting in changes in the bacterial community (Figure 6, Supplementary Table S3), and bacteria are sensitive to the moisture content of the soil (Figure 6, [67,68]). The changes in SMC caused by ammonium nitrogen addition may have an impact on the growth and metabolic activities of bacteria, causing changes in the composition of the soil bacterial phyla [69,70,71].
At the genus level, Candidatus Udaeobacter, Candidatus Solibacter, Bryobacter, Bradyrhizobium, and Acidothermus were the main bacterial genera in all soil samples. Nitrogen addition significantly altered the relative abundance of Bradyrhizobium and Acidothermus (Figure 3, Supplementary Figure S5). These are important bacterial communities participating in the carbon and nitrogen cycles in soils [72,73,74]. In our study, LN treatment significantly increased the relative abundance of Bradyrhizobium and Acidothermus. These changes can be explained by the co-trophic hypothesis. Bradyrhizobium and Acidothermus are classified as copiotrophs [56], and nitrogen addition enriched the level of available N in the soil and therefore stimulated their growth [75]. Xu et al. [76] found that nitrogen addition increased the relative abundance of Bradyrhizobium, which is consistent with our results. Bradyrhizobium is a parthenogenic nitrogen-fixing microorganism involved in the nitrate and nitrite reduction pathways [77,78]. The process of biological nitrogen fixation requires large amounts of energy and is mainly dependent on soil nutrients, which are increased by the addition of nitrogen, thus increasing the abundance of the nitrogen-fixing bacteria Bradyrhizobium. However, Wan et al. (2021) found that nitrogen addition decreased the relative abundance of Bradyrhizobium. The differences may be explained by differences in the concentration of nitrogen applied. In our study, LN significantly increased the relative abundance of Bradyrhizobium, but HN showed inhibition, although not significantly (Supplementary Figure S7). The implication is that there is a degree of concentration between LN and HN, which when exceeded will inhibit the growth of Bradyrhizobium [79]. Furthermore, the pH in the soil may have a strong influence on Acidothermus, which are acidophilic genus [73]. The reduction in soil pH provided a suitable environment for the growth of Acidothermus (Table 1), increasing the relative abundance of Acidothermus [80]. However, Wan et al. [81] found that N addition reduced the relative abundance of Acidothermus in orchard ecosystems, and previous studies found that Acidothermus was positively correlated with soil organic matter content [73], which is contrary to our findings. This may be explained by differences in the ecosystem types and the geographical regions and vegetations of different studies. This study was in a temperate natural wetland and the aboveground vegetations were dominated by D. angustifolia, so the difference may be attributed to differences in land use patterns and species of plants on the aboveground, compared to previous studies. Moreover, bacterial community composition is not only influenced by the physical and chemical properties of the soil, but also by altitude [82] and temperature and rainfall [83], and the dominating factors may be different in different regions or at different scales.

4.3. Changes in the Structure of Soil Bacterial Networks

Analysis of ecological network structure can reveal the relationship between the complexity and stability of ecosystems [84]. In accordance with our first hypothesis, our study showed that nitrogen addition changed potential bacterial competition and network complexity (Figure 5, Table 2). Compared with CK, LN and HN increased the average clustering coefficient of the network but decreased the average geodesic distance of the network. This implied that the application of NH4+ might remit the competition within bacteria. The nitrogen addition enhanced the efficiency of inter-species transport, which rapidly distributed the effects of disturbance across the network and contributed to a reduction in network stability [85,86]. NH4+ input changed the total links and average degree of the bacterial network, indicating that nitrogen addition altered the network complexity. It may be that the nitrogen application treatment selectively impacted the reproduction and growth of soil microorganisms, causing changes in the interactions between microorganisms and thereby affecting the complexity of the network structure. Among these, positive links dominated in the three treatments, implying that mutual cooperation more than competitive exclusion played an important role in the microbial community. The LN reduced the interactions between species, either by mutual cooperation or competitive exclusion. The amount of NH4+ may be insufficient to significantly alleviate the limit of soil nitrogen use by the soil bacteria. Consistent with Sha et al. [49], who found that nitrogen addition reduced positive mutualism and depressed coprological interactions between soil microbial species, Yu et al. 2019 showed that continued urea addition reduces the redundancy of species interactions. Under conditions of limited N content, most of the nitrogen applied in soil may be absorbed by vegetation rather than by soil bacteria [49,87] because vegetation may be more competitive than soil bacteria. The HN increased mutual cooperation between species but decreased interspecies competition. The application of NH4+ provided an enriched nitrogen source for bacteria, which reduced competition for limited resources and promoted cooperation among species. Consistent with Li et al. [13], who investigated the effect of nitrogen fertilization on the network structure of soil microorganisms, this study found that excess N addition enriched a number of mutualistic microbes and increased the percentage of positive links. However, the input of NH4+ decreased the number of modules and negative links in the bacterial network, indicating that nitrogen addition reduced the stability of the soil bacterial community [88,89]. Exogenous nitrogen inputs may cluster species that are sensitive to nitrogen fertilization and thereby have an impact on the bacterial community [13]. Similar results were obtained by Yao et al. [58], who found that high nitrogen addition reduced the stability of the microbial systems. The decrease in stability may be attributed to a loss of diversity and changes in community composition caused by nitrogen enrichment effects [90].

5. Conclusions

Our study investigated the effects of more than five years ammonium addition on soil bacterial diversity and composition in the D. angutifolia Wetland of Sanjiang Plain. High nitrogen addition significantly reduced the heterogeneity of the soil bacterial community and had no significant effect on the richness, whereas low nitrogen addition did not alter the soil bacterial diversity. NH4+ addition changed the bacterial composition and the community network complexity and stability in the D. angutifolia Wetland of Sanjiang Plain. SMC, NO3, and pH were key soil parameters driving the composition of soil bacterial community structures. These results provide a theoretical basis for studying the impacts of future climate change on the diversity and functional stability of wetland ecosystems.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/microorganisms11082033/s1.

Author Contributions

X.W., M.W. and X.S. did the experiments, analyzed the data, and prepared this original manuscript; X.S., B.F. and M.L. revised this manuscript and edited the linguistics; Y.L. and R.Z. helped to analyze the data and do the field experiment; H.N. and M.L. designed this experiment. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Natural Sciences Foundation of Heilongjiang Province (LH2020C088); the Heilongjiang Province Postdoctoral Research Start-up Fund Project (LBH-Q21167); the Outstanding Youth Foundation of Heilongjiang University (JCL202006); Heilongjiang Provincial Ecological Environmental Protection Research Project (HST2022ST008); Basic Scientific Research project of Higher Education Institutions of Heilongjiang Province and the China Scholarship Council Visiting Scholar Program (201908230401).

Data Availability Statement

All raw sequences have been deposited into an NCBI Sequence Read Archive with the BioProject accession: PRJNA880534.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationship that could have appeared to influence the work reported in this paper.

References

  1. Galloway, J.N.; Dentener, F.J.; Capone, D.G.; Boyer, E.W.; Howarth, R.W.; Seitzinger, S.P.; Asner, G.P.; Cleveland, C.C.; Green, P.A.; Holland, E.A.; et al. Nitrogen cycles: Past, present, and future. Biogeochemistry 2004, 70, 153–226. [Google Scholar] [CrossRef]
  2. Hu, Y.; Chen, M.; Yang, Z.; Cong, M.; Zhu, X.; Jia, H. Soil microbial community response to nitrogen application on a swamp meadow in the arid region of central asia. Front. Microbiol. 2022, 12, 797306. [Google Scholar] [CrossRef] [PubMed]
  3. Zeng, J.; Liu, X.; Song, L.; Lin, X.; Zhang, H.; Shen, C.; Chu, H. Nitrogen fertilization directly affects soil bacterial diversity and indirectly affects bacterial community composition. Soil Biol. Biochem. 2016, 92, 41–49. [Google Scholar] [CrossRef]
  4. Knelman, J.E.; Nemergut, D.R. Changes in community assembly may shift the relationship between biodiversity and ecosystem function. Front. Microbiol. 2014, 5, 424. [Google Scholar] [CrossRef] [PubMed]
  5. Averill, C.; Waring, B. Nitrogen limitation of decomposition and decay: How can it occur? Glob. Chang. Biol. 2018, 24, 1417–1427. [Google Scholar] [CrossRef] [Green Version]
  6. Huang, L.B.; Bai, J.H.; Wen, X.J.; Zhang, G.L.; Zhang, C.D.; Cui, B.S.; Liu, X.H. Microbial resistance and resilience in response to environmental changes under the higher intensity of human activities than global average level. Glob. Chang. Biol. 2020, 26, 2377–2389. [Google Scholar] [CrossRef]
  7. Carreiro, M.M.; Sinsabaugh, R.L.; Repert, D.A.; Parkhurst, D.F. Microbial Enzyme Shifts Explain Litter Decay Responses to Simulated Nitrogen Deposition. Ecology 2000, 81, 2359. [Google Scholar] [CrossRef]
  8. Song, C.C.; Liu, D.Y.; Yang, G.S.; Song, Y.Y.; Mao, R. Effect of nitrogen addition on decomposition of Calamagrostis angustifolia litters from freshwater marshes of Northeast China. Ecol. Eng. 2011, 37, 1578–1582. [Google Scholar] [CrossRef]
  9. Song, C.C.; Wang, L.L.; Tian, H.Q.; Liu, D.Y.; Lu, C.Q.; Xu, X.F.; Zhang, L.H.; Yang, G.S.; Wan, Z.M. Effect of continued nitrogen enrichment on greenhouse gas emissions from a wetland ecosystem in the Sanjiang Plain, Northeast China: A 5-year nitrogen addition experiment. J. Geophys. Res. Biogeosci. 2013, 118, 741–751. [Google Scholar] [CrossRef]
  10. Bardgett, R.D.; Wardle, D.A. Aboveground-Belowground Linkages: Biotic interactions, Ecosystem Processes, and Global Change; Oxford University Press: New York, NY, USA, 2010. [Google Scholar]
  11. Peralta, R.M.; Ahn, C.; Gillevet, P.M. Characterization of soil bacterial community structure and physicochemical properties in created and natural wetlands. Sci. Total Environ. 2013, 443, 725–732. [Google Scholar] [CrossRef]
  12. Fierer, N.; Jackson, R.B. The diversity and biogeography of soil bacterial communities. Proc. Natl. Acad. Sci. USA 2006, 103, 626–631. [Google Scholar] [CrossRef] [PubMed]
  13. Li, B.B.; Roley, S.S.; Duncan, D.S.; Guo, J.; Quensen, J.F.; Yu, H.-Q.; Tiedje, J.M. Long-term excess nitrogen fertilizer increases sensitivity of soil microbial community to seasonal change revealed by ecological network and metagenome analyses. Soil Biol. Biochem. 2021, 160, 108349. [Google Scholar] [CrossRef]
  14. Xie, L.; Zhang, Q.; Cao, J.; Liu, X.; Xiong, D.; Kong, Q.; Yang, Y. Effects of Warming and Nitrogen Addition on the Soil Bacterial Community in a Subtropical Chinese Fir Plantation. Forests 2019, 10, 861. [Google Scholar] [CrossRef] [Green Version]
  15. Jia, M.; Gao, Z.; Gu, H.; Zhao, C.; Liu, M.; Liu, F.; Xie, L.; Wang, L.; Zhang, G.; Liu, Y.; et al. Effects of precipitation change and nitrogen addition on the composition, diversity, and molecular ecological network of soil bacterial communities in a desert steppe. PLoS ONE 2021, 16, e248194. [Google Scholar] [CrossRef] [PubMed]
  16. Ma, X.; Wang, T.; Shi, Z.; Chiariello, N.R.; Docherty, K.; Field, C.B.; Gutknecht, J.; Gao, Q.; Gu, Y.; Guo, X.; et al. Long-term nitrogen deposition enhances microbial capacities in soil carbon stabilization but reduces network complexity. Microbiome 2022, 10, 112. [Google Scholar] [CrossRef]
  17. Xu, L.; Zhang, B.; Peng, X.; Zhang, X.; Sun, B.; Sun, H.; Jiang, C.; Zhou, S.; Zeng, X.; Bai, Z.; et al. Dynamic variations of microbial community structure in Myriophyllum aquaticum constructed wetlands in response to different NH4+-N concentrations. Process Biochem. 2020, 93, 55–62. [Google Scholar] [CrossRef]
  18. Li, Y.; Tian, D.; Wang, J.; Niu, S.; Tian, J.; Ha, D.; Qu, Y.; Jing, G.; Kang, X.; Song, B. Differential mechanisms underlying responses of soil bacterial and fungal communities to nitrogen and phosphorus inputs in a subtropical forest. PeerJ 2019, 7, e7631. [Google Scholar] [CrossRef]
  19. Xu, R.H.; Tan, M.; Zhu, J.F.; Liu, Z.H. The response of soil microbial diversity to the difference of N deposition concentration in alpine wetland. J. Biol. 2021, 38, 75–81. [Google Scholar]
  20. Lu, G.; Xie, B.; Cagle, G.A.; Wang, X.; Han, G.; Wang, X.; Hou, A.; Guan, B. Effects of simulated nitrogen deposition on soil microbial community diversity in coastal wetland of the yellow river delta. Sci. Total Environ. 2021, 757, 143825. [Google Scholar] [CrossRef]
  21. Zhang, R.T.; Liu, Y.N.; Zhong, H.X.; Chen, X.W.; Sui, X. Effects of simulated nitrogen deposition on the soil microbial community diversity of a Deyeuxia angustifolia wetland in the Sanjiang plain, northeastern China. Ann. Microbiol. 2022, 72, 11. [Google Scholar] [CrossRef]
  22. Liu, X.J.; Zhang, Y.; Han, W.X.; Tang, A.H.; Shen, J.L.; Cui, Z.L.; Vitousek, P.; Erisman, J.W.; Goulding, K.; Christie, P.; et al. Enhanced nitrogen deposition over China. Nature 2013, 494, 459–462. [Google Scholar] [CrossRef] [PubMed]
  23. Fang, H.J.; Cheng, S.L.; Yu, G.R.; Zheng, J.J.; Zhang, P.L.; Xu, M.J.; Li, Y.N.; Yang, X.M. Responses of CO2 efflux from an alpine meadow soil on the Qinghai Tibetan Plateau to multi-form and low-level N addition. Plant Soil. 2012, 351, 177–190. [Google Scholar] [CrossRef]
  24. Fang, H.; Cheng, S.; Yu, G.; Cooch, J.; Wang, Y.; Xu, M.J.; Li, L.S.; Dang, X.S.; Li, Y.N. Low-level nitrogen deposition significantly inhibits methane uptake from an alpine meadow soil on the Qinghai-Tibetan Plateau. Geoderma 2014, 213, 444–452. [Google Scholar] [CrossRef]
  25. Kou, L.; Zhang, X.; Wang, H.; Yang, H.; Zhao, W.; Li, S. Nitrogen additions inhibit nitrification in acidic soils in a subtropical pine plantation: Effects of soil pH and compositional shifts in microbial groups. J. For. Res. 2018, 30, 669–678. [Google Scholar] [CrossRef]
  26. Matson, P.A.; McDowell, W.H.; Townsend, A.R.; Vitousek, P.M. The globalization of N deposition: Ecosystem consequences in tropical environments. Biogeochemistry 1999, 46, 67–83. [Google Scholar] [CrossRef]
  27. Recous, S.; Machet, J.M.; Mary, B. The partitioning of fertilizer-N between soil and crop-comparison of ammonium and nitrate applications. Plant Soil. 1992, 144, 101–111. [Google Scholar] [CrossRef]
  28. Zhang, Y.T.; Shen, H.; He, X.H.; Thomas, B.W.; Lupwayi, N.Z.; Hao, X.Y.; Thomas, M.C.; Shi, X.J. Fertilization shapes bacterial community structure by alteration of soil pH. Front. Microbiol. 2017, 8, 1325. [Google Scholar] [CrossRef] [Green Version]
  29. Liu, C.; Dong, Y.; Sun, Q.; Jiao, R. Soil bacterial community response to short-term manipulation of the nitrogen deposition form and dose in a Chinese fir plantation in southern China. Water Air Soil Poll. 2016, 227, 447. [Google Scholar] [CrossRef]
  30. Geng, J.; Fang, H.; Cheng, S.; Pei, J. Effects of N deposition on the quality and quantity of soil organic matter in a boreal forest: Contrasting roles of ammonium and nitrate. Catena 2021, 198, 104996. [Google Scholar] [CrossRef]
  31. Weng, X.H.; Sui, X.; Liu, Y.N.; Yang, L.B.; Zhang, R. Effect of nitrogen addition on the carbon metabolism of soil microorganisms in a Calamagrostis angustifolia wetland of the Sanjiang Plain, northeastern China. Ann. Microbiol. 2022, 72, 18. [Google Scholar] [CrossRef]
  32. Lei, K.; Zhang, M.X. The wetland resources in China and the conservation advices. Wetl. Sci. 2005, 3, 81–86. [Google Scholar]
  33. He, D.M.; Ruan, H.H. Long term effect of land reclamation from lake on chemical composition of soil organic matter and its mineralization. PLoS ONE 2014, 9, e99251. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Wang, Y.; Liu, J.; He, L. Responses of carbon dynamics to nitrogen deposition in typical freshwater wetland of Sanjiang Plain. J. Chem. 2014, 2014, 603948. [Google Scholar] [CrossRef] [Green Version]
  35. Wen, Z.; Xu, W.; Li, Q.; Han, M.J.; Tang, A.; Zhang, Y.; Luo, X.S.; Shen, J.L.; Wang, W.; Li, K.H.; et al. Changes of nitrogen deposition in China from 1980 to 2018. Environ. Int. 2020, 144, 106022. [Google Scholar] [CrossRef]
  36. Sui, X.; Zhang, R.T.; Liu, Y.N.; Xu, N.; Ni, H.W. Influence of simulation nitrogen deposition on soil microbial functional diversity of Calamagrostis angustifolia wetland in Sanjiang Plain. Acta. Agrestia Sin. 2016, 24, 1226–1233. [Google Scholar]
  37. Luan, Z.Q.; Deng, W.; Bai, J.H.; Zhu, B.G. Protection of Honghe national reserve wetland habitat. Res. Soil Water Conserv. 2003, 10, 154–157. [Google Scholar]
  38. Jones, D.L.; Willett, V.B. Experimental evaluation of methods to quantify dissolved organic nitrogen (DON) and dissolved organic carbon (DOC) in soil. Soil Biol. Biochem. 2005, 38, 991–999. [Google Scholar] [CrossRef]
  39. Lane, D. 16S/23S rRNA sequencing. In Nucleic Acid Techniques in Bacterial Systematics; Stackebrandt, E., Goodfellow, M., Eds.; John Wiley and Sons: New York, NY, USA, 1991. [Google Scholar]
  40. Parks, D.H.; Tyson, G.W.; Hugenholtz, P.; Beiko, R.G. STAMP: Statistical analysis of taxonomic and functional profiles. Bioinformatics 2014, 30, 3123–3124. [Google Scholar] [CrossRef] [Green Version]
  41. Harpole, W.S.; Tilman, D. Grassland species loss resulting from reduced niche dimension. Nature 2007, 446, 791–793. [Google Scholar] [CrossRef]
  42. Zhang, X.; Han, X. Nitrogen deposition alters soil chemical properties and bacterial communities in the inner mongolia grassland. J. Environ. Sci. 2012, 24, 1483–1491. [Google Scholar] [CrossRef]
  43. Weng, X.H.; Sui, X.; Li, M.S.; Zhang, R.; Yang, L.B. Effect of simulated nitrogen deposition on soil microbial carbon metabolism in Calamagrostis angustifolia wetland in Sanjiang Plain. Huan Jing Ke Xue 2022, 43, 4674–4683. [Google Scholar]
  44. Wang, C.; Liu, D.; Bai, E. Decreasing soil microbial diversity is associated with decreasing microbial biomass under nitrogen addition. Soil Biol. Biochem. 2018, 120, 126–133. [Google Scholar] [CrossRef]
  45. Li, Z.M.; Shen, J.P.; Zhang, L.M.; Liu, G.P.; Bai, P.M.; He, J.Z. Effect of stimulated nitrogen deposition on the bacterial community structure of semiarid temperate grassland. Environ. Sci. 2018, 39, 5665–5671. [Google Scholar]
  46. Frey, B.; Carnol, M.; Dharmarajah, A.; Brunner, I.; Schleppi, P. Only minor changes in the soil microbiome of a sub-alpine forest after 20 years of moderately increased nitrogen loads. Front. For. Glob. Chang. 2020, 3, 77. [Google Scholar] [CrossRef]
  47. Deng, J.J.; Zhou, Y.B.; Yin, Y.; Wei, Y.W.; Qin, S.J.; Zhu, W.X. Soil bacterial community structure characteristics in coniferous forests of Montane Regions of eastern Liaoning Province, China. Acta. Ecol. Sin. 2019, 39, 997–1008. [Google Scholar]
  48. Sui, X.; Zhang, R.; Frey, B.; Yang, L.; Liu, Y.; Ni, H.; Li, M.H. Soil physicochemical properties drive the variation in soil microbial communities along a forest successional series in a degraded wetland in northeastern China. Ecol. Evol. 2021, 11, 2194–2208. [Google Scholar] [CrossRef]
  49. Sha, M.; Xu, J.; Zheng, Z.; Fa, K. Enhanced atmospheric nitrogen deposition triggered little change in soil microbial diversity and structure in a desert ecosystem. Glob. Ecol. Conserv. 2021, 31, e01879. [Google Scholar] [CrossRef]
  50. Gavrichkova, O.; Kuzyakov, Y. Ammonium versus nitrate nutrition of Zea mays and Lupinus albus: Effect on root-derived CO2 efflux. Soil Biol. Biochem. 2008, 40, 2835–2842. [Google Scholar] [CrossRef]
  51. Thirukkumaran, C.M.; Parkinson, D. Microbial respiration, biomass, metabolic quotient and litter decomposition in a lodgepole pine forest floor amended with nitrogen and phosphorous fertilizers. Soil Biol. Biochem. 2000, 32, 59–66. [Google Scholar] [CrossRef]
  52. Tischner, R. Nitrate uptake and reduction in higher and lower plants. Plant Cell Environ. 2000, 23, 1005–1024. [Google Scholar] [CrossRef]
  53. Wang, C.Y.; Lv, Y.N.; Liu, X.Y.; Wang, L. Ecological effects of atmospheric nitrogen deposition on soil enzyme activity. J. For. Res. 2013, 24, 109–114. [Google Scholar] [CrossRef]
  54. Xie, F.; Ma, A.; Zhou, H.; Liang, Y.; Yin, J.; Ma, K.; Zhuang, X.L.; Zhuang, G.Q. Revealing fungal communities in alpine wetlands through species diversity, functional diversity and ecological network diversity. Microorganisms. 2020, 8, 632. [Google Scholar] [CrossRef] [PubMed]
  55. Fierer, N.; Bradford, M.A.; Jackson, R.B. Toward an ecological classification of soil bacteria. Ecology. 2007, 88, 1354–1364. [Google Scholar] [CrossRef]
  56. Fierer, N.; Lauber, C.L.; Ramirez, K.S.; Zaneveld, J.; Bradford, M.A.; Knight, R. Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients. ISME J. 2012, 6, 1007–1017. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Dai, Z.; Su, W.; Chen, H.; Barberan, A.; Zhao, H.; Yu, M.; Yu, L.; Brookes, P.C.; Schadt, C.W.; Chang, S.X.; et al. Long-term nitrogen fertilization decreases bacterial diversity and favors the growth of actinobacteria and proteobacteria in agro-ecosystems across the globe. Glob. Chang. Biol. 2018, 24, 3452–3461. [Google Scholar] [CrossRef] [PubMed]
  58. Yao, M.J.; Rui, J.P.; Li, J.B.; Dai, Y.M.; Bai, Y.F.; Heděnec, P.; Wang, J.M.; Zhang, S.H.; Pei, K.Q.; Liu, C.; et al. Rate-specific responses of prokaryotic diversity and structure to nitrogen deposition in the Leymus chinensis steppe. Soil Biol. Biochem. 2014, 79, 81–90. [Google Scholar] [CrossRef]
  59. Lauber, C.L.; Hamady, M.; Knight, R.; Fierer, N. Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. Appl. Environ. Microb. 2009, 75, 5111–5120. [Google Scholar] [CrossRef] [Green Version]
  60. Lauber, C.L.; Strickland, M.S.; Bradford, M.A.; Fierer, N. The influence of soil properties on the structure of bacterial and fungal communities across land-use types. Soil Biol. Biochem. 2008, 40, 2407–2415. [Google Scholar] [CrossRef]
  61. Wei, Z.W.; Li, Y.Y.; Jiang, W.; Liao, X.R. Diversity of Acidobacteria in rhizosphere soils of common trees in Wuxi. Chin. J. Ecol. 2018, 37, 2649–2656. [Google Scholar]
  62. Rousk, J.; Bååth, E.; Brookes, P.C.; Lauber, C.L.; Lozupone, C.; Caporaso, J.G.; Knight, R.; Fierer, N. Soil bacterial and fungal communities across a pH gradient in an arable soil. ISME J. 2010, 4, 1340–1350. [Google Scholar] [CrossRef]
  63. Ramirez, K.S.; Lauber, C.L.; Knight, R.; Bradford, M.A.; Fierer, N. Consistent effects of nitrogen fertilization on soil bacterial communities in contrasting systems. Ecology 2010, 91, 3463–3470. [Google Scholar] [CrossRef]
  64. Zhang, C.; Zhang, X.Y.; Zou, H.T.; Kou, L.; Yang, Y.; Wen, X.F.; Li, S.G.; Wang, H.M.; Sun, X.M. Contrasting effects of ammonium and nitrate additions on the biomass of soil microbial communities and enzyme activities in subtropical China. Biogeosciences 2017, 14, 4815–4827. [Google Scholar] [CrossRef] [Green Version]
  65. Breuker, A.; Stadler, S.; Schippers, A. Microbial community analysis of deeply buried marine sediments of the New Jersey shallow shelf (IODP Expedition 313). FEMS Microbiol. Ecol. 2013, 85, 578–592. [Google Scholar] [CrossRef] [Green Version]
  66. Lovisa, B.R.; Philip, H.; Tyson, G.W.; Blackall, L.L. Filamentous Chloroflexi (green non-sulfur bacteria) are abundant in wastewater treatment processes with biological nutrient removal. Microbiology 2002, 148, 2309. [Google Scholar]
  67. Degens, B.P.; Schipper, L.A.; Sparling, G.P.; Vojvodic-Vukovic, M. Decreases in organic C reserves in soils can reduce the catabolic diversity of soil microbial communities. Soil. Biol. Biochem. 2000, 32, 189–196. [Google Scholar] [CrossRef]
  68. Zhang, N.; Wan, S.; Guo, J.; Han, G.; Gutknecht, J.; Schmid, B.; Yu, L.; Liu, W.; Bi, J.; Wang, Z.; et al. Precipitation modifies the effects of warming and nitrogen addition on soil microbial communities in northern Chinese grasslands. Soil Biol. Biochem. 2015, 89, 12–23. [Google Scholar] [CrossRef] [Green Version]
  69. Keeler, B.L.; Hobbie, S.E.; Kellogg, L.E. Effects of long-term nitrogen addition on microbial enzyme activity in eight forested and grassland sites: Implications for litter and soil organic matter decomposition. Ecosystems 2009, 12, 1–15. [Google Scholar] [CrossRef]
  70. Liu, W.X.; Zhang, Z.; Wan, S.Q. Predominant role of water in regulating soil and microbial respiration and their responses to climate change in a semiarid grassland. Glob. Chang. Biol. 2009, 15, 184–195. [Google Scholar] [CrossRef]
  71. Zak, D.R.; Kling, G.W. Microbial community composition and function across an arctic tundra landscape. Ecology 2006, 87, 1659–1670. [Google Scholar] [CrossRef] [PubMed]
  72. Anderson, C.R.; Hamonts, K.; Clough, T.J.; Condron, L.M. Biochar does not affect soil N-transformations or microbial community structure under ruminant urine patches but does alter relative proportions of nitrogen cycling bacteria. Agric. Ecosyst. Environ. 2014, 191, 63–72. [Google Scholar] [CrossRef]
  73. Du, S.Y.; Yu, M.; Liu, F.H.; Xiao, L.L.; Zhang, H.X.; Tao, J.; Gu, W.; Gu, J.Y.; Chen, X. Effect of facility management regimes on soil bacterial diversity and community structure. Chin. J. Eco-Agric. 2017, 25, 1615–1625. [Google Scholar]
  74. Fei, Y.C.; Ye, Y.Q.; Zheng, H.; Lu, J.; You, Y.F.; Huang, Y.; Chen, A.L.; Cao, G.Q. Effect of C/N ratio under the control of exogenous nitrogen on the bacterial community structure in litters of Cunninghamia lanceolate plantation under different undergrowth vegetation management measures. Acta Ecol. Sin. 2021, 41, 2011–2023. [Google Scholar]
  75. Yao, F.; Yang, S.; Wang, Z.; Wang, X.; Ye, J.; Wang, X.; DeBruyn, J.M.; Feng, X.; Jiang, Y.; Li, H. Microbial taxa distribution is associated with ecological trophic cascades along an elevation gradient. Front. Microbiol. 2017, 8, 2071. [Google Scholar] [CrossRef] [Green Version]
  76. Xu, R.Y.; Zuo, M.X.; Yuan, Y.L.; Sun, J.; Lu, Y.S.; Jie, K.Z.; Xu, P.Z. Effects of nitrogen fertilizer dosage optimization on nitrogen uptake content and utilization efficiency and microbial function genes of nitrogen cycle in sweet corn. J. South. Agr. 2020, 51, 2919–2926. [Google Scholar]
  77. Polcyn, W.; Luciński, R. Aerobic and anaerobic nitrate and nitrite reduction in free-living cells of Bradyrhizobium sp. (Lupinus). FEMS Microbiol. Lett. 2003, 226, 331–337. [Google Scholar] [CrossRef] [Green Version]
  78. Rivas, R.; Martens, M.; Lajudie, P.; De Willems, A. Multilocus sequence analysis of the genus Bradyrhizobium. Syst. Appl. Microbiol. 2009, 32, 101–110. [Google Scholar] [CrossRef]
  79. Chen, C.L.; Wu, M.N.; Wei, W.X. Effect of long-term application of nitrogen fertilizer on the diversity of nitrifying genes (amoA and hao) in Paddy soil. Environ. Sci. 2011, 32, 1489–1496. [Google Scholar]
  80. Khan, A.; Jiang, H.; Bu, J.; Adnan, M.; Gilanni, S.W.; Hussain, M.A.; Zhang, M.Q. Untangling the Rhizosphere Bacterial Community Composition and Response of Soil Physiochemical Properties to Different Nitrogen Applications in Sugarcane Field. Fron. Microbiol. 2022, 13, 856078. [Google Scholar] [CrossRef]
  81. Wan, Y.; Li, W.; Wang, J.; Shi, X. Bacterial diversity and community in response to long-term nitrogen fertilization gradient in citrus orchard soils. Diversity 2021, 13, 282. [Google Scholar] [CrossRef]
  82. Shen, C.; Ni, Y.; Liang, W.; Wang, J.; Chu, H. Distinct soil bacterial communities along a small-scale elevational gradient in alpine tundra. Front. Microbiol. 2015, 6, 582. [Google Scholar] [CrossRef] [Green Version]
  83. Zhang, Y.; Dong, S.K.; Gao, Q.Z.; Liu, S.L.; Zhou, H.K.; Ganjurjav, H.; Wang, X.X. Climate change and human activities altered the diversity and composition of soil microbial community in alpine grasslands of the Qinghai-Tibetan Plateau. Sci. Total Environ. 2016, 562, 353–363. [Google Scholar] [CrossRef] [PubMed]
  84. Montoya, J.M.; Pimm, S.L.; Solé, R.V. Ecological networks and their fragility. Nature 2006, 442, 259–264. [Google Scholar] [CrossRef]
  85. Helbing, D. Globally networked risks and how to respond. Nature 2013, 497, 51–59. [Google Scholar] [CrossRef]
  86. Yuan, M.M.; Guo, X.; Wu, L.W.; Zhang, Y.; Xiao, N.J.; Ning, D.L.; Shi, Z.; Zhou, X.S.; Wu, L.Y.; Yang, Y.F. Climate warming enhances microbial network complexity and stability. Nat. Clim. Chang. 2021, 11, 343–348. [Google Scholar] [CrossRef]
  87. Han, X. Effects of Soil Nitrogen to Plant—Microbial on Nitrogen Competition in Temperate Forest; Beijing Forestry University: Beijing, China, 2015. [Google Scholar]
  88. Coyte, K.Z.; Schluter, J.; Foster, K.R. The ecology of the microbiome: Networks, competition, and stability. Science 2015, 350, 663–666. [Google Scholar] [CrossRef] [PubMed]
  89. Zhou, H.; Gao, Y.; Jia, X.; Wang, M.; Ding, J.; Cheng, L.; Bao, F.; Wu, B. Network analysis reveals the strengthening of microbial interaction in biological soil crust development in the Mu Us Sandy Land, northwestern China. Soil Biol. Biochem. 2020, 144, 107782. [Google Scholar] [CrossRef]
  90. Yu, Z.; Hu, X.; Wei, D.; Liu, J.; Zhou, B.; Jin, J.; Liu, X.; Wang, G. Long-term inorganic fertilizer use influences bacterial communities in mollisols of northeast China based on high-throughput sequencing and network analyses. Arch. Agron. Soil Sci. 2019, 65, 1331–1340. [Google Scholar] [CrossRef]
Figure 1. Effects of N application on soil microbial community α-diversity. (a) Shannon index analysis of bacteria. (b) Simpson index analysis of bacteria. (c) Chao index analysis of bacteria. (d) OTU numbers of bacteria. Lowercase letters indicate significant differences (p < 0.05). The asterisk indicates the average of each group. CK: control; LN: low ammonium; HN: high ammonium.
Figure 1. Effects of N application on soil microbial community α-diversity. (a) Shannon index analysis of bacteria. (b) Simpson index analysis of bacteria. (c) Chao index analysis of bacteria. (d) OTU numbers of bacteria. Lowercase letters indicate significant differences (p < 0.05). The asterisk indicates the average of each group. CK: control; LN: low ammonium; HN: high ammonium.
Microorganisms 11 02033 g001
Figure 2. Non-metric multidimensional scaling (NMDS) of bacterial communities under different N addition levels. CK: control; LN: low ammonium; HN: high ammonium.
Figure 2. Non-metric multidimensional scaling (NMDS) of bacterial communities under different N addition levels. CK: control; LN: low ammonium; HN: high ammonium.
Microorganisms 11 02033 g002
Figure 3. Relative abundance of bacterial phyla. (a) Top 10 phyla and genera, (b) top 20 genera under different N addition levels. “Other” of the legend in (a) includes 34 phyla except top 10 bacterial phyla. “Other” of the legend in (b) includes 53 genera except after top 20 bacterial genera. CK: control; LN: low ammonium; HN: high ammonium.
Figure 3. Relative abundance of bacterial phyla. (a) Top 10 phyla and genera, (b) top 20 genera under different N addition levels. “Other” of the legend in (a) includes 34 phyla except top 10 bacterial phyla. “Other” of the legend in (b) includes 53 genera except after top 20 bacterial genera. CK: control; LN: low ammonium; HN: high ammonium.
Microorganisms 11 02033 g003
Figure 4. The cladogram of the soil bacteria community under different N addition levels (LDA score = 3.5). CK: control; LN: low ammonium; HN: high ammonium.
Figure 4. The cladogram of the soil bacteria community under different N addition levels (LDA score = 3.5). CK: control; LN: low ammonium; HN: high ammonium.
Microorganisms 11 02033 g004
Figure 5. Network structures of bacteria under different N addition levels. Each node signifies an OTU. The size of the node represents the degree of node. Colors of the nodes indicate different phyla. The values in the circle represent the OTU numbers of the corresponding phylum. Colors of the lines represent the strength of the correlation between OTUs. A red line indicates a negative interaction between two individual nodes, while a blue line indicates a negative interaction. CK: control; LN: low ammonium; HN: high ammonium.
Figure 5. Network structures of bacteria under different N addition levels. Each node signifies an OTU. The size of the node represents the degree of node. Colors of the nodes indicate different phyla. The values in the circle represent the OTU numbers of the corresponding phylum. Colors of the lines represent the strength of the correlation between OTUs. A red line indicates a negative interaction between two individual nodes, while a blue line indicates a negative interaction. CK: control; LN: low ammonium; HN: high ammonium.
Microorganisms 11 02033 g005
Figure 6. Redundancy analysis (RDA) of the relationships between soil properties and bacterial community structure under N different addition. CK: control; LN: low ammonium; HN: high ammonium.
Figure 6. Redundancy analysis (RDA) of the relationships between soil properties and bacterial community structure under N different addition. CK: control; LN: low ammonium; HN: high ammonium.
Microorganisms 11 02033 g006
Figure 7. Pearson correlations of bacterial abundant phyla (a), and genera (b), with soil properties. * indicates 0.5 level; ** indicates 0.01 level; *** indicates 0.001 level. The color gradient (pink to yellow) represents the correlation degree from higher to lower. Pink represents positive correlation and yellow represents negative correlation. CK: control; LN: low ammonium; HN: high ammonium.
Figure 7. Pearson correlations of bacterial abundant phyla (a), and genera (b), with soil properties. * indicates 0.5 level; ** indicates 0.01 level; *** indicates 0.001 level. The color gradient (pink to yellow) represents the correlation degree from higher to lower. Pink represents positive correlation and yellow represents negative correlation. CK: control; LN: low ammonium; HN: high ammonium.
Microorganisms 11 02033 g007
Table 1. Physicochemical characteristics of the soil samples for the different N treatments.
Table 1. Physicochemical characteristics of the soil samples for the different N treatments.
TreatmentsSMCpHNO3
(mg/kg)
NH4+
(mg/kg)
DOC
(mg/kg)
DON
(mg/kg)
TN
(%)
SOC
(%)
CK0.43 ± 0.06 b4.41 ± 0.43 a23.83 ± 6.34 b21.15 ± 3.45 a3897.88 ± 625.37 b299.74 ± 130.71 b0.57 ± 0.23 b5.34 ± 2.23 a
LN0.48 ± 0.01 ab4.19 ± 0.06 ab23.27 ± 3.49 b5.62 ± 0.63 b2571.62 ± 294.54 c520.08 ± 166.01 ab0.68 ± 0.05 ab4.53 ± 0.42 a
HN0.58 ± 0.05 a4.02 ± 0.11 b40.60 ± 6.17 a4.72 ± 0.06 b5004.81 ± 331.11 a719.68 ± 156.87 a1.18 ± 0.29 a6.86 ± 2.00 a
Note: Each treatment was performed in triplicate. The data are expressed as the mean ± standard deviation; lowercase letters indicate significant differences (p < 0.05). CK: control; LN: low ammonium; HN: high ammonium. SMC: soil moisture content; DOC: dissolved organic carbon; DON: dissolved organic nitrogen; TN: total nitrogen; SOC: soil organic carbon.
Table 2. Parameters of bacterial network under different N addition levels.
Table 2. Parameters of bacterial network under different N addition levels.
TreatmentOTU NumbersTotal LinksPositive LinkNegative LinkAverage DegreeAverage Clustering CoefficientAverage Geodesic DistanceNetwork DensityModularity (No. of Modules)
CK46511,9026637526525.60.78411.7230.0550.812 (11)
LN46410,9225809511323.50.79111.6040.0510.824 (10)
HN45812,5427329521327.40.81211.4090.0600.779 (9)
CK: control; LN: low ammonium; HN: high ammonium.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Weng, X.; Wang, M.; Sui, X.; Frey, B.; Liu, Y.; Zhang, R.; Ni, H.; Li, M. High Ammonium Addition Changes the Diversity and Structure of Bacterial Communities in Temperate Wetland Soils of Northeastern China. Microorganisms 2023, 11, 2033. https://doi.org/10.3390/microorganisms11082033

AMA Style

Weng X, Wang M, Sui X, Frey B, Liu Y, Zhang R, Ni H, Li M. High Ammonium Addition Changes the Diversity and Structure of Bacterial Communities in Temperate Wetland Soils of Northeastern China. Microorganisms. 2023; 11(8):2033. https://doi.org/10.3390/microorganisms11082033

Chicago/Turabian Style

Weng, Xiaohong, Mingyu Wang, Xin Sui, Beat Frey, Yingnan Liu, Rongtao Zhang, Hongwei Ni, and Maihe Li. 2023. "High Ammonium Addition Changes the Diversity and Structure of Bacterial Communities in Temperate Wetland Soils of Northeastern China" Microorganisms 11, no. 8: 2033. https://doi.org/10.3390/microorganisms11082033

APA Style

Weng, X., Wang, M., Sui, X., Frey, B., Liu, Y., Zhang, R., Ni, H., & Li, M. (2023). High Ammonium Addition Changes the Diversity and Structure of Bacterial Communities in Temperate Wetland Soils of Northeastern China. Microorganisms, 11(8), 2033. https://doi.org/10.3390/microorganisms11082033

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop