Next Article in Journal
Comparative Pathogenicity of Aeromonas spp. in Cultured Red Hybrid Tilapia (Oreochromis niloticus × O. mossambicus)
Previous Article in Journal
Biological Crusts to Increase Soil Carbon Sequestration: New Challenges in a New Environment
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biology
Volume 10
Issue 11
10.3390/biology10111191
biology-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nitrogen and Sulfur Additions Improved the Diversity of nirK- and nirS-Type Denitrifying Bacterial Communities of Farmland Soil

by
Xuan Chen
1,
Hui Wei
1,2,3,4 and
Jiaen Zhang
1,2,3,4,5,*
1
College of Natural Resources and Environment, South China Agricultural University, Guangzhou 510642, China
2
Guangdong Engineering Research Center for Modern Eco-Agriculture and Circular Agriculture, South China Agricultural University, Guangzhou 510642, China
3
Guangdong Provincial Key Laboratory of Eco-Circular Agriculture, South China Agricultural University, Guangzhou 510642, China
4
Key Laboratory of Agro-Environment in the Tropics, Ministry of Agriculture and Rural Affairs, South China Agricultural University, Guangzhou 510642, China
5
Guangdong Laboratory for Lingnan Modern Agriculture, South China Agricultural University, Guangzhou 510642, China
*
Author to whom correspondence should be addressed.
Biology 2021, 10(11), 1191; https://doi.org/10.3390/biology10111191
Submission received: 6 October 2021 / Revised: 6 November 2021 / Accepted: 8 November 2021 / Published: 16 November 2021
(This article belongs to the Section Microbiology)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

In this study, we investigated the effects of combined N and S additions on soil properties and N-cycling microorganisms. We found that soil treated with N and S additions had higher Shannon diversity and Chao 1 richness for the nirK- and nirS-related OTUs. We also found that, in the nirK-type denitrifying community, the addition of N and S promoted increases in the relative abundance of Bradyrhizobium, and in terms of the nirS-type denitrifying community, the N and S additions increased the relative abundance of Thiobacillus. Moreover, the N and S additions increased the activity of nitrate reductase and nitrite reductase. Most importantly, the nirK-type denitrifying community demonstrated a higher sensitivity to N and S additions. We expect that this research will be useful for the study of the nitrogen cycle in soil.

Abstract

Anthropogenic nitrogen (N) and sulfur (S) deposition can change above- and belowground biodiversity, including soil microbial diversity. The diversity of denitrifying microorganisms is of great significance to the calculation of the global nitrogen cycle and nitrogen flux. For a long time, nirK and nirS have been used as the functional genes to study denitrifying microorganisms, and have gradually become molecular markers for studying the composition and diversity of denitrifying bacteria. Here, three-time exposures to N and S applications (7, 30, and 60 days), were independently established. Additionally, the abundance, diversity, and structure of nirK- and nirS-type denitrifying communities were examined by sequencing analyses in response to three treatments, namely, N and S (TN/S), sodium chloride (TNaCl) and deionized water (pH = 7.0) (CK). Our results suggest that TN/S led to higher electrical conductivity (EC), total nitrogen (TN), total organic carbon (TOC), nitrate nitrogen (NO3−-N), ammonium nitrogen (NH4+-N), and lower pH compared with TNaCl and CK, which affected the diversity of nirK- and nirS-type denitrifying bacterial communities. We also observed that the nirK-type denitrifying community demonstrated a higher sensitivity to N and S additions. Overall, our results are important for the understanding of nitrogen in soil and N2O emissions.

1. Introduction

In recent years, the rapid development of industry and agriculture has led to an increase in global atmospheric N deposition [1,2,3,4,5,6,7]. Over the past two centuries, anthropogenic N input in ecosystems has increased by nearly 10 times [8], and it is expected to double by 2050 [2]. The imbalance of the natural N cycle has become increasingly serious. Overall, the impact of N deposition exists on a global scale and has a profound impact on ecosystems [9,10], especially on terrestrial ecosystems [11,12]. At the same time, N deposition also affects other ecological types, such as oceans [13], lakes [14], and seaport systems [15]. The dramatic increase in N deposition leads to changes in the N cycle [16], which has adverse effects on ecosystems, such as the eutrophication and acidification of ecosystems [17], and changes in microbial community structures and diversity in soils [18].
Soil microorganisms are important parts of the ecosystem [19], and they play an important role in material circulation, litter decomposition, and maintaining the stability of the ecosystem [20,21]. Therefore, it is valuable to explore the effect of N deposition on soil microorganisms. The current research on the effects of N deposition on soil microorganisms mainly focuses on microbial biomass, community structure and diversity, and microbial functional diversity. Studies have shown that the community structure and diversity of soil microorganisms respond differently under different N deposition conditions. Sui et al. [22] found that N deposition increased the diversity of bacteria in the Sanjiang Plain, and that the bacterial diversity was the highest when the level of N was low, indicating that the impact of N deposition on the increase in bacterial diversity had a certain concentration limit. The addition of N can also change the abundance of different phyla of soil microorganisms. Fierer et al. [23] and Ramirez et al. [24] found that copiotrophic taxa (including members of the Proteobacteria and Bacteroidetes phyla) typically increased in relative abundance in the high-N plots, with oligotrophic taxa (Acidobacteria and Verrucomicrobia) exhibiting the opposite pattern. Moreover, Cederlund et al. [25] have shown that N fertilization reduced the relative abundances of Acidobacteria, Gemmatimonadetes, and Verrucomicrobia. In addition, other studies have shown that N deposition affects the diversity of microbial functional genes by changing soil microbial activity. Yu et al. [26] reported that for the soil microbial gram aerobic bacteria, the metabolic function was enhanced significantly under different N additions, while the metabolic activity of soil microbial fungi and yeast was unaffected by different N additions. Wu et al. [27] also found that N deposition enhanced the metabolic function of microorganisms.
Sulfur, another nutrient element, also plays an important role in the nutrient cycle of the ecosystem. However, atmospheric sulfur deposition is the main factor in soil acidification in terrestrial ecosystems, which causes the release of base ions, nutrient imbalance, and has an impact on the structure and function of soil biological communities [28]. Studies have shown that soil acidification not only reduces the number of microorganisms but also inhibits their growth and activity [29,30]. Wang et al. [31] confirmed that Actinomycetes and Bacteria were greatly reduced and their activity was decreased due to soil acidification, while fungi were less affected by soil acidification. Acidification also affects the reproduction of microbial populations. Sierra et al. [32] indicated that it is not conducive to the survival of Pseudomonas and Corynebacterium. In addition, Wang et al. [33] found that the relative abundance of Lysobacter and Flavobacterium was significantly increased, while the relative abundance of Massilia was significantly reduced. Overall, high levels of N and S emissions have caused a large amount of acid deposition, but previous studies paid more attention to the negative effects of acid rain, while potential fertilization effects of acid rain-induced N and S inputs have rarely been studied.
Denitrification is an important link in the complete cycle of nitrogen [34]. Heterotrophic denitrification comprises sequential reduction steps from NO3− to NO2−, NO, N2O, and finally to elemental N2 [35]. Specific enzymes, namely NO3− reductase (Nar), NO2− reductase (Nir), NO reductase (Nor), and N2O reductase (Nos) regulate the rate and product stoichiometry of heterotrophic denitrification [36,37]. From an ecological point of view, the denitrification process is a self-balancing mechanism that can protect the system from instability and imbalance [38]. Moreover, denitrification is one of the main pathways of nitrogen loss [39]. Marcel et al. [40] showed that nitrogen which was recycled through denitrification accounted for 52% of the total global nitrogen input to the soil, while the denitrification loss of nitrogen fertilizer accounted for 12% of nitrogen fertilizer applications [41]. In addition, denitrification is the main source of greenhouse gases, such as N2O [42]. In 1990, the application of chemical nitrogen fertilizer resulted in about 1.5 × 106 t N2O-N, accounting for 44% of the N2O-N release into the atmosphere by human activities [3], but by 2020, the figure in China had reached 4.30 times that of 1990 [42]. Overall, the denitrification process has a significant impact on the nitrogen cycle, fertilizer loss, and the environment. In-depth research into this is helpful for us to further quantify nitrogen turnover, predict the fate of nitrogen fertilizer in farmland ecosystems, formulate sewage treatment measures, and deal with global warming caused by the greenhouse effect.
The reduction of nitrite to NO is the most important step in denitrification (Figure 1), which is catalyzed by Nir. There are two types of enzymes involved in this process. One is the soluble copper-containing enzyme (Cu-Nir), which is encoded by the nirK gene. The other is cytochrome reductase (cd1-Nir), which is encoded by the nirS gene. These two enzymes perform the same function, but feature different structures and catalytic sites and generally do not exist in the same denitrification species at the same time. They represent two denitrification groups with different ecological characteristics and occupy different niches [43]. Generally, nirS-type denitrifying bacteria feature higher gene abundance and dominate the environment, while nirK-type denitrifying bacteria exist in a wider microbial population and cover more diverse taxonomic units [39]. Studies have shown that nirK-type denitrifying bacteria exhibit more powerful habitat selection in response to different environmental gradients [43,44]. NirK and nirS were the first functional genes used to study denitrifying microorganisms [45], and have gradually become molecular markers for studying the composition and diversity of denitrifying bacteria. Studies have shown that nitrogen addition significantly increased the abundance of the NO2− reducing genes (nirK and nirS) in Chinese fir plantations [46]. There is also a study indicating that long-term fertilizer applications had a significant impact on the size of the nirK community, but had little impact on the nirS community [47]. Previous studies have shown that N deposition levels [48] and N forms [49] can affect the abundance and community compositions of N-cycle functional genes, whereas few studies have addressed the influence of S additions and, in particular, the interaction between S and N on the paddy soil microbial communities involved in N cycling. Therefore, it is necessary to better understand how mixed additions of N and S influence the N cycling process in soil ecosystems.
Additionally, in complex and variable soil environments, the coexistence of nirS- and nirK-type denitrifying communities leads to cooperation and competition, which are factors determining community variation [50,51]. Therefore, in this regard, the potential ecological roles of microorganisms may be revealed by the network analyses used to study the co-existence within functional groups [52,53]. However, for now, few such studies have focused on functionally equivalent genes [51]. Thus, the study of functional equivalent genes may help to improve our understanding of soil denitrifying microbial communities.
In this study, we investigated the effects of combined N and S additions on soil properties and N-cycling microorganisms. Based on the molar ratio of S to N in acid rain of 2.75:1 [54], we conducted this experiment to study how N and S additions affect denitrifying microbial communities and enzyme activities in the paddy soil. We hypothesized that: (1) the soil treated with N and S additions would have higher Shannon diversity and Chao 1 richness for nirK- and nirS-related operational taxonomic units (OTUs); (2) at the genus level, in the case of the nirK-type denitrifying community, the application of N and S could promote increases in the relative abundance of Bradyrhizobium, while in the nirS-type denitrifying community, the application of N and S might increase the relative abundance of Thiobacillus; (3) N and S additions would increase the activity of nitrate reductase and nitrite reductase. We expected that we could use this study to better evaluate the impact of N and S deposition on ecosystems, and that it might provide an effective reference for the study of the nitrogen cycle.

2. Materials and Methods

2.1. Site Description and Experimental Design

The experimental soils were collected from an ecological farm (23°06′ N, 113°20′ E) of South China Agricultural University in Guangzhou City, Guangdong Province, China. The climate in this area is classified as a subtropical monsoon climate, according to the Köppen–Geiger climate classification. In addition, the average annual air temperature and precipitation are 22 °C and 1800 mm, respectively. The soil is classified as Rhodudult [55]. The soil pH was 6.2, with an organic matter content of 19.64 g·kg−1, a TN of 0.98 g·kg−1, and an EC of 73.33 μS·cm−1. Furthermore, the other chemical properties of the soil, including contents of the available phosphorous (AP), the available potassium (AK), NH4+-N, and NO3−-N, were 40.88, 144.76, 10.70, and 28.82 mg·kg−1, respectively.
Three-time exposures to N and S applications (7, 30, and 60 days) were independently established, and three treatments with five replications of each were set up, including TN/S, TNaCl, and CK. Among them, CK was set as the control. N was added as sodium nitrate (NaNO3) and S was added as sodium sulfate (Na2SO4). To observe the effects of Na input, we set up an extra TNaCl treatment with an input of the same amount of Na as in the TN/S treatment. NaNO3, Na2SO4, and NaCl were purchased from Guangzhou Congyuan Instrument Ltd., Guangzhou, China. According to the characteristics of the N and S input associated with acid rain in Guangzhou, a mixed solution of NaNO3 and Na2SO4 was prepared with a molar ratio of S to N of 2.75:1 [54], that is, the concentrations of NaNO3 and Na2SO4 were 0.60 mol·L−1 and 1.06 mol·L−1, respectively, and they were sealed and stored for reserve. Before the experiment began, each conical flask (200 mL) was filled with 150 g of soil and 27 mL of treatment solution, respectively. All the soil samples were placed in a temperature-controlled (28 °C) incubator (RXZ-500A) in the dark and the soil moisture was controlled at 60% of the water holding capacity throughout the incubation period. At the end of each incubation time, 3 g of soil from each sample were stored at −20 °C for DNA extraction and the remaining sample was kept at 4 °C for further analysis of the soil’s physicochemical properties and enzyme activities.

2.2. Soil Sampling and Physicochemical Analysis

The NO3−-N and NH4+-N were measured using the methods described by Bao [56]. The EC and pH were measured in a 1:5 (w:v) aqueous extract, with a conductivity meter and pH meter (Crison mod. 2001, Barcelona, Spain), respectively. The TN and TOC concentrations were determined using an elemental analyzer (Vario TOC cube). The nitrate reductase assays were carried out according to the method described by Li et al. [57]. One unit of enzyme activity (U) was defined as the number of milligrams of nitrite produced per g of soil sample in 24 h at 30 °C. The nitrite reductase assays were determined by the Griess reagent colorimetric method [58]. One unit of enzyme activity (U) was defined as the number of milligrams of the difference between the content of nitrite nitrogen before and after the enzymatic reaction in per g of soil sample in 24 h at 30 °C.

2.3. Microbial DNA Extraction, nirK, and nirS Genes Amplification and Sequencing

The batch tests samples obtained on days 7, 30, and 60 were used for DNA extraction in triplicate with the OMEGA soil DNA extraction kit (USA), following the manufacturer’s instructions. This study used 1% agarose gel electrophoresis to determine the quality and concentration of the extracted DNA, which was finally checked with a NanoDrop ND-2000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA) and stored at −20 °C until use. The primers nirK1aCuF (ATCATGGTSCTGCCGCG) and nirKR3CuR (GCCTCGATCAGRTTGTGGTT) were used to amplify the hypervariable region fragments of nirK genes [59], and the primers cd3aF (GTSAACGTSAAGGARACSGG) and R3cdR (GASTTCGGRTGSGTCTTGA) were used to amplify the hypervariable region fragments of nirS genes [57]. The 20-μL PCR reaction consisted of a 4 μL of 5 × FastPfu Buffer, 2 μL of 2.5 mmol·L dNTPs, 0.8 μL of each primer (5 μM, Forward and Reverse), 0.4 μL of FastPfu Polymerase, 0.2 μL of BSA and 10 ng of template DNA. The thermal cycling programs were as follows: initial denaturation at 95 °C for 3 min; 40 cycles of 95 °C for 30 s, 60 °C for 30 s, 72 °C for 45 s; and a final extension at 72 °C for 10 min. The PCR products were purified using an AxyPrep DNA Gel Extraction Kit (Axygen Bio-Sciences, Union City, CA, USA), and then the purified products of nirS and nirK were pooled in equimolar quantities. Finally, according to the standard protocols issued by Majorbio Bio-pharm Technology Co., Ltd. (Shanghai, China), paired-end sequencing (2 × 250) was performed using the Illumina MiSeq platform (Illumina, San Diego, CA, USA). QIIME was used to filter the raw data and screen out sequences shorter than 150 bp, low-quality sequences, and mismatched sequences, and the filtered data were normalized at the same depth by Mothur software. In addition, through the UPARSE pipeline 7.1 [60], all the optimized sequences were clustered into OTUs according to the principle of 97% similarity clustering. An average of 194,400 nirK and 214,806 nirS effective sequence reads were obtained from each sample, and the average lengths of nirK and nirS were 450 bp and 400 bp, separately.

2.4. Statistical Analysis

According to the identification results of the OTU division and classification level, the specific composition of each sample at each classification level (kingdom, phylum, class, order, family, genus, species) was obtained by QIIME version 1.8.0. Moreover, two α diversity indices (Chao 1 index and Shannon index) for each sample were calculated by QIIME version 1.8.0. One-way ANOVA was performed using SPSS16.0. The similarity analysis (ANOSIM) and non-metric multidimensional scaling (NMDS: based on Bray–Curtis distance) of the nirK- and nirS-type denitrifying communities between different treatments were carried out through the “Vegan” package, and the data were visualized by the “ggplot2” package in R software (Version 4.0.5). In addition, CANNOCO 5.0 software was used for redundancy analysis (RDA) of dominant denitrifying microorganisms and soil physicochemical properties. Meanwhile, the Spearman test was used to analyze the correlation between them, and the correlation heat map was drawn by using R software.

3. Results

3.1. Physico-Chemical Analysis of Soil Samples

Relative to the TNaCl and CK, we found that the soil physicochemical properties changed significantly at 7-day incubation and 60-day incubation in the TN/S (Table 1), including significant increases in the EC, TN, TOC, NO3−-N, and NH4+-N, and a significant reduction in the pH (p < 0.05). The 30-day incubation did not significantly change the investigated soil physicochemical properties (p > 0.05, Table 1).

3.2. Alpha-Diversity Indices of nirS- and nirK-Type Denitrifying Communities

In terms of the nirK gene, a total of 312 OTUs for species classification were obtained by clustering. On the 7th day, the total OTUs of TN/S, TNaCl, and CK were 191, 203, and 213, respectively. Similarly, on the 30th day and the 60th day, they were 221, 210, and 220 and 197, 219, and 228, respectively. Among them, there were 110 OTUs shared by all groups (Figure 2a). As for the nirS gene, a total of 751 OTUs for species classification were obtained by clustering. With seven-day incubation, the total OTUs of TN/S, TNaCl, and CK were 529, 490, and 561, respectively. Similarly, they were 533, 513, and 493 and 508, 508, and 513 with 30-day incubation and 60-day incubation, respectively. Of these, there were 259 OTUs common to all groups (Figure 2b). After OTU identification, the OTUs of each sample were divided into six classification levels, including phylum, class, order, family, genus, and species. The results showed that the nirK-type denitrifying microbial communities were more abundant. In terms of the nirK gene, on the 7th day, the numbers of OTUs at the phylum, class, and order levels under each treatment were the same, while the numbers of OTUs at the family, genus, and species levels were lower than those of the control, while after 30, 60 days of incubation, the numbers of OTUs in the soil treated with TN/S at the family, genus, and species levels were higher than those of the control (Figure 3a). In terms of the nirS gene, we found that there was a non-significant difference in the number of OTUs between the treatments at the phylum and class levels. However, a significant difference was observed, that is, compared with CK, TN/S increased the numbers of OTUs at the order, family, genus, and species levels in all three time periods (Figure 3b).
Compared with CK, TN/S resulted in a significantly higher Shannon diversity of nirS-related OTUs at 30-day incubation and 60-day incubation, and a significantly lower Chao 1 richness at 60-day incubation. Conversely, the TN/S-treated soils demonstrated higher Shannon diversity and Chao 1 richness for nirK-related OTUs relative to that of TNaCl and CK throughout the incubation period (Table 2).

3.3. Beta Diversity of nirS- and nirK-Type Denitrifying Microorganisms

NMDS analysis showed that the nirK-type community at 7-day incubation was separated from those in the other two periods (30-day incubation and 60-day incubation) (ANOSIM R = 0.540, p = 0.032) (Figure 4a). Moreover, in the nirS-type denitrifying community, we found that the microbial community changed significantly from the 7th day to the 60th day, and the similarity of the microbial community changed from small to large and then to small in different periods. However, such shifts were not statistically significant (ANOSIM R = 0.161, p = 0.201) (Figure 4b). Interestingly, we found a significant similarity in the community composition of soil nirS denitrifying microorganisms between TN/S and TNaCl at 7-day incubation. We also found a similarity between TN/S and CK at 60-day incubation, while there was a significant difference in the community composition of soil nirK denitrifying microorganisms between each treatment at each stage. The results indicated that treatment duration could affect the community composition of denitrifying microorganisms, and the nirK-type denitrifying community was more responsive to treatment duration.
Under different treatment conditions, the application of N and S changed the compositions of the nirK- and nirS-type denitrifying communities. At the genus level, Bradyrhizobium was the dominant nirK genus, accounting for 17.2–32.9% of the total nirK-type denitrifying community sequences in all soil samples (Figure 5a). Compared with TNaCl and CK, we found that in the TN/S treatment groups, the relative abundance of Bradyrhizobium significantly increased at 30-day incubation, but reduced at 60-day incubation. In the nirS-type denitrifying community, we observed that the predominant genus was Thiobacillus, accounting for 3.6–9.8% of the total nirS sequences (Figure 5b). Compared with the other two treatments, it is evident that the composition of the nirS-type denitrifying community was not sensitive to TN/S.

3.4. Enzyme Activities of Soil Samples

As expected, the nitrate reductase and the nitrite reductase were affected significantly by N and S additions (p < 0.05, Table 3), especially at 7-day incubation. We found that in the case of the nitrate reductase, it increased by 6.63% and 6.20%, respectively (p < 0.05), during that period when compared to the TNaCl and CK. Similarly, it increased by 34.71% and 37.98%, respectively (p < 0.05), over the same time for the nitrite reductase.

3.5. Correlations between Functional Genes and Soil Properties

A redundancy analysis (RDA) was used to explore the environmental factors affecting the composition of nirK- and nirS-type denitrifying microbial communities (Figure 6), and eight soil physical and chemical factors (pH, EC, TOC, NO3−N, and NH4+-N, nitrate reductase, nitrite reductase) affecting the microbial community were selected and analyzed with an OTUs table. In the case of the nirK-type denitrifying communities, pH and NO3−-N were significantly positively correlated with axis 1, while nitrite reductase was significantly negatively correlated (Figure 6a). In addition, according to the projection of physicochemical factors on the first axis, pH may have been the first influencing factor. Similarly, in the nirS-type denitrifying communities, TN and NH4+-N were significantly positively correlated with axis 1, while nitrite reductase was significantly negatively correlated (Figure 6b). Moreover, TN may have been the first influencing factor through the projection of physicochemical factors on the first axis.
In order to select the 50 most abundant bacterial genera, the Spearman correlation coefficient was used to analyze the correlation between soil physicochemical factors and the composition of soil nirK- and nirS-type denitrifying microbial communities in the three treatments (Figure 7). The results showed that there were differences in the correlation between the two genotypes of denitrifying microbial communities and soil physicochemical factors, but the difference was more obvious in the nirK genes. In the case of the nirK-type denitrifying communities, we found that only two genera were significantly associated with pH at 7-day incubation, namely Mesorhizobium and Afipia (Figure 7a) (p ≤ 0.001). At 30-day incubation, three genera were significantly related to pH, TN, NH4+-N, nitrate reductase, and nitrite reductase, respectively, namely Mesorhizobium, Bradyrhizobium, and Bosea (Figure 7a) (p ≤ 0.001); while at the 60-day incubation, four genera were significantly related to pH and nitrate reductase, respectively, namely Norank_p_environmental_samples_k_norank, Nitrosopira, Ensifer, and Bradyrhizobium (Figure 7a) (p ≤ 0.001). In terms of the nirS-type denitrifying communities, Thiobacillus was only significantly associated with pH at 7-day incubation. However, at 30-day incubation and 60-day incubation, Thiobacillus was also related to nitrate reductase, nitrite reductase, and EC, TN, NH4+-N, respectively (Figure 7b) (p ≤ 0.001).

4. Discussion

4.1. Impact of N and S Additions on the Community Structure and Diversity of Soil nirK- and nirS-Type Denitrifying Microorganisms

Denitrifying microorganisms widely exist among bacteria and archaea, and denitrification is also found in some fungal mitochondria [37]. In this study, N and S applications were found to have a significant effect on soil properties, which in turn contributed to promoting changes in the abundance and community compositions of nirK- and nirS-type denitrifying bacteria. Through the present research, we found that the application of N and S increased the content of TOC in the soil, which provided essential nutrients for the increase of bacterial biomass [62]. Therefore, we found that bacteria were mainly involved in denitrification in the three treatments, and many unidentified gene sequences were also found. As for whether there are denitrifying fungi in the soil, further analysis and the identification of the nirK and nirS gene sequences of these unknown species are also needed. Among the identified nirK and nirS denitrifying microorganisms, Proteobacteria was the most prominent among the three treatments, and fertilization significantly increased its relative abundance. It was confirmed that Proteobacteria can use difficult biodegradable carbon sources in acidic environments and decompose them into small molecular substances to provide nutrients for other microorganisms [63]. This is consistent with the results of this study. N and S applications decreased the pH of the soil, which created a favorable acidic environment for Proteobacteria. Thus, the abundance of Proteobacteria was higher under fertilization. It is worth noting that the abundance of the nirK gene was significantly correlated with soil pH. As mentioned above, soil pH is a key factor affecting the responses of denitrifier genes to N fertilization, which may be because pH is closely related to the metabolic substrate content of dissolved organic carbon and nitrate [49,64,65,66] and directly or indirectly affects denitrifying communities. Other associated factors include NO3−-N. As one of the initial substrates of denitrification, NO3−-N is an important factor in controlling the denitrification of nirK-type denitrifying bacteria [67]. At the genus level, in terms of the nirK gene, our results demonstrated that Bradyrhizobium was the core genus, which was consistent with the results of Fan et al. [68]. Previous studies reported that the dominant genera, including Bradyrhizobium and Burkholderia were facultative nitrogen-fixing bacteria, which were found in Alpha-Proteobacteria and Beta-Proteobacteria, respectively, and they were usually classified as eutrophic organisms. These advantageous eutrophic organisms can use large amounts of inorganic substances, such as NO3−-N, as energy sources to reproduce rapidly in a nutrient-rich environment [69,70,71]. The results of this study also supported the hypothesis that the application of N increased the content of NO3−-N in the soil compared with the other two treatments, which were more conducive to the growth of Bradyrhizobium. As for the nirS gene, our results indicated that Thiobacillus, a typical sulfur-oxidizing bacteria, was the dominant genus. In our study, the addition of S led to an increase in the SO42− content in the soil, with nitrite gradually replaced by sulfate, resulting in a significant increase in the relative abundance of Thiobacillus [72]. Additionally, in our study, TN and NH4+-N were significantly related to nirS, but not correlated with nirK denitrifying bacteria. Previous studies demonstrated that the composition and abundance of denitrifying bacteria could be affected by indirectly creating an environment favorable to denitrifying bacteria [49,73,74,75]. Other studies also reported that the abundance and community structure of soil bacteria in farmland were affected by many factors. Different water management systems caused differences in pH, inorganic nitrogen content, microbial biomass carbon content, microbial biomass nitrogen content, soluble organic nitrogen content, and soluble organic carbon content. All of these factors may lead to differences in soil bacterial abundance and community structure [76,77,78,79].
The results of this study indicated that the addition of N and S significantly increased the alpha diversity index of nirK-type denitrifying communities, while the nirS-OTU-related Shannon and Chao 1 indices were less affected, suggesting that alpha diversity was not necessarily affected by the change in population sizes and compositions of denitrifying communities [80,81]. In conclusion, compared with the change in nirS-type denitrifying communities at 60-day incubation, the nirK-type denitrifying communities responded more strongly to the application of N and S. These observations confirmed that nirK- and nirS-type denitrifying communities respond differently to N and S applications, and it is clear that the most sensitive is the nirK-type denitrifying community, as was revealed in previous studies [82,83].

4.2. Impact of Treatment Duration on the Diversity of Soil nirK- and nirS-Type Denitrifying Microorganisms

The NMDS analysis showed that treatment duration could affect the community composition of denitrifying microorganisms. This may be because with the increase in time, the contents of carbon and nitrogen change, which can affect the diversity of soil-denitrifying microorganisms [84]. On the other hand, pH is also one of the factors influencing microbial diversity [35,85]. When the soil pH value exceeds a certain range (niche), the net growth of a single taxon that cannot survive decreases, which may alter the competitive outcome [45,86]. However, in this study, the pH is within the optimal range of all microorganisms (5.9–6.8), so there are no extreme conditions. At present, there are few studies on the effect of treatment duration on the diversity of soil nirK- and nirS-type denitrifying substances, so further study is needed in the future.

5. Conclusions

To summarize, genes for denitrification are ubiquitous in soil microorganisms. The abundance of nirK and nirS genes was linked to soil pH, NO3−-N concentration, TN concentration, and NH4+-N concentration. Through the quantitation of nirS and nirK, the addition of N and S had an effect on soil denitrification. Therefore, further study of the application of N and S will contribute to our knowledge of how the soil community composition affects N2O flux. The difference between nirK and nirS abundance may provide a more accurate result than using single genes.

Author Contributions

J.Z. contributed to the conception of the study and supervised this study. X.C. collected and carried out the data analysis, made charts, and wrote the manuscript. H.W. provided constructive help and comments for this study. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number U1701236, the Key Joint Team Project of Guangdong Laboratory for Lingnan Modern Agriculture, grant number NT2021010, and the Science and Technology Planning Project of Guangdong Province, grant number 2019B030301007.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data availability statement

All the data generated or analyzed during this study are included in this published article, and the data are available from the corresponding author on reasonable request.

Acknowledgments

Authors are deeply indebted to all the colleagues who helped us during the sampling process, and we thank the reviewers for their valuable comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kaiser, J. Environmental Policy: The Other Global Pollutant: Nitrogen Proves Tough to Curb. Science 2001, 294, 1268–1269. [Google Scholar] [CrossRef] [PubMed]
  2. Galloway, J.N.; Dentener, F.J.; Capone, D.G.; Boyer, E.; Howarth, R.; Seitzinger, S.P.; Asner, G.; 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]
  3. Galloway, J.N.; Cowling, E.B. Reactive Nitrogen and the World: 200 Years of Change. Ambio 2002, 31, 64–71. [Google Scholar] [CrossRef]
  4. Liu, X.; Duan, L.; Mo, J.; Du, E.; Shen, J.; Lu, X.; Zhang, Y.; Zhou, X.; He, C.; Zhang, F. Nitrogen Deposition and its Ecological Impact in China: An overview. Environ. Pollut. 2011, 159, 2251–2264. [Google Scholar] [CrossRef] [PubMed]
  5. Liu, X.; Zhang, Y.; Han, W.; Tang, A.; Shen, J.; Cui, Z.; 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]
  6. Gu, B.; Ju, X.; Chang, J.; Ge, Y.; Vitousek, P.M. Integrated Reactive Nitrogen Budgets and Future Trends in China. Proc. Natl. Acad. Sci. USA 2015, 112, 8792–8797. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Song, L.; Kuang, F.; Skiba, U.; Zhu, B.; Liu, X.; Levy, P.; Dore, A.; Fowler, D. Bulk Deposition of Organic and Inorganic Nitrogen in Southwest China from 2008 to 2013. Environ. Pollut. 2017, 227, 157–166. [Google Scholar] [CrossRef] [PubMed]
  8. Sutton, M.A.; Bleeker, A. Environmental Science: The Shape of Nitrogen to come. Nature 2013, 494, 435–437. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Stevens, C.J.; Dise, N.B.; Mountford, J.O.; Gowing, D.J. Impact of Nitrogen Deposition on the Species Richness of Grasslands. Science 2004, 303, 1876–1879. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Clark, C.M.; Tilman, D. Loss of Plant Species after Chronic Low-Level Nitrogen Deposition to Prairie Grasslands. Nature 2008, 451, 712–715. [Google Scholar] [CrossRef]
  11. Asner, G.P.; Townsend, A.R.; Riley, W.J.; Matson, P.A.; Neff, J.C.; Cleveland, C.C. Physical and Biogeochemical Controls over Terrestrial Ecosystem Responses to Nitrogen Deposition. Biogeochemistry 2001, 54, 1–39. [Google Scholar] [CrossRef]
  12. Matson, P.; Lohse, K.A.; Hall, S.J. The Globalization of Nitrogen Deposition: Consequences for Terrestrial Ecosystems. Ambio 2002, 31, 113–119. [Google Scholar] [CrossRef] [PubMed]
  13. Zhang, Y.; Yu, Q.; Ma, W.; Chen, L. Atmospheric Deposition of Inorganic Nitrogen to the Eastern China Seas and its Implications to Marine Biogeochemistry. J. Geophys. Res. Space Phys. 2010, 115. [Google Scholar] [CrossRef]
  14. Hobbs, W.O.; Lafrancois, B.M.; Stottlemyer, R.; Toczydlowski, D.; Engstrom, D.R.; Edlund, M.B.; Almendinger, J.E.; Strock, K.E.; VanderMeulen, D.; Elias, J.E.; et al. Nitrogen Deposition to Lakes in National Parks of the Western Great Lakes Region: Isotopic Signatures, Watershed Retention, and Algal Shifts. Glob. Biogeochem. Cycles 2016, 30, 514–533. [Google Scholar] [CrossRef] [Green Version]
  15. Pakeman, R.J.; Alexander, J.; Brooker, R.; Cummins, R.; Fielding, D.; Gore, S.; Hewison, R.; Mitchell, R.; Moore, E.; Orford, K.; et al. Long-term Impacts of Nitrogen Deposition on Coastal Plant Communities. Environ. Pollut. 2016, 212, 337–347. [Google Scholar] [CrossRef] [PubMed]
  16. Fowler, D.; Coyle, M.; Skiba, U.; Sutton, M.A.; Cape, J.; Reis, S.; Sheppard, L.J.; Jenkins, A.; Grizzetti, B.; Galloway, J.N.; et al. The Global Nitrogen Cycle in the Twenty-First Century. Philos. Trans. R. Soc. B Biol. Sci. 2013, 368, 20130164. [Google Scholar] [CrossRef] [PubMed]
  17. Bouwman, A.F.; Van Vuuren, D.P.; Derwent, R.G.; Posch, M. A Global Analysis of Acidification and Eutrophication of Terrestrial Ecosystems. Water Air Soil Pollut. 2002, 141, 349–382. [Google Scholar] [CrossRef]
  18. Bowman, W.D.; Cleveland, C.C.; Halada, L.; Hreško, J.; Baron, J.S. Negative Impact of Nitrogen Deposition on Soil Buffering Capacity. Nat. Geosci. 2008, 1, 767–770. [Google Scholar] [CrossRef]
  19. Ding, S.Y.; She, J.Y.; Yang, Q.P.; Wang, S.L. Effects of Thinning and Pruning on Soil Microbial Biomass Carbon and Soil Enzyme Activities in Chinese Fir Plantation. J. Cent. South Univ. For. Technol. 2015, 35, 75–79. [Google Scholar]
  20. Kazda, M. Indications of Unbalanced Nitrogen Nutrition of Norway Spruce Stands. Plant Soil 1990, 128, 97–101. [Google Scholar] [CrossRef]
  21. Lovett, G.M.; Reiners, W.A.; Olson, R.K. Cloud Droplet Deposition in Subalpine Balsam Fir Forests: Hydrological and Chemical Inputs. Science 1982, 218, 1303–1304. [Google Scholar] [CrossRef]
  22. 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]
  23. Fierer, N.; Lauber, C.L.; Ramirez, K.; Zaneveld, J.; Bradford, M.; Knight, R. Comparative Metagenomic, Phylogenetic and Physiological Analyses of Soil Microbial Communities across Nitrogen Gradients. ISME J. 2012, 6, 1007–1017. [Google Scholar] [CrossRef] [Green Version]
  24. Ramirez, K.S.; Lauber, C.L.; Knight, R.; Bradford, M.; Fierer, N. Consistent Effects of Nitrogen Fertilization on Soil Bacterial Communities in Contrasting Systems. Ecology 2010, 91, 3463–3470. [Google Scholar] [CrossRef] [PubMed]
  25. Cederlund, H.; Wessén, E.; Enwall, K.; Jones, C.; Juhanson, J.; Pell, M.; Philippot, L.; Hallin, S. Soil Carbon Quality and Nitrogen Fertilization Structure Bacterial Communities with Predictable Responses of Major Bacterial Phyla. Appl. Soil Ecol. 2014, 84, 62–68. [Google Scholar] [CrossRef]
  26. Yu, P.Y.; Zhu, F.; Wang, Z.Y.; Yan, W.D.; Su, S.F.; Li, T.P. Effects of Nitrogen Addition on Metabolic Function of Microbial Community in Red Soil of Cinnamomum camphora Forest. J. Cent. South Univ. For. Technol. 2013, 33, 70–74. [Google Scholar]
  27. Wu, S.Q.; Wang, C.Z.; Li, M.S. On Soil Functional Diversity of Native Coastal Wetland under Simulated Nitrogen Deposition. Soils 2017, 49, 1153–1158. [Google Scholar]
  28. Zhang, A.; Olatunji, O.A.; Tariq, A.; Li, T.; Wang, R.; Jiang, Y. Sulfur Deposition Changed the Community Structure of Soil Nematodes by Affecting Omnivores-Predators. Sci. Total Environ. 2021, 771, 144912. [Google Scholar] [CrossRef]
  29. Piertri, J.C.A.; Brookes, P.C. Relationships between Soil pH and Microbial Properties in a UK Arable Soil. Soil Biol. Biochem. 2008, 40, 1856–1861. [Google Scholar]
  30. Guo, L.L.; Yuan, Z.G.; Zhu, W.W.; Yi, Z.X.; Tu, N.M. Research Progress in Soil Acidification Effect on Soil Biological Characteristics. Hunan Agric. Sci. 2014, 24, 30–32. [Google Scholar]
  31. Wang, F.G.; Song, L.; Feng, Y.; Hong, Y.C.; Cui, D.J.; Yuan, Y.B. Characteristics of Soil Microbiology in Different Planting-Life Orchard Acid Soils. Chin. J. Soil Sci. 2011, 42, 46–50. [Google Scholar]
  32. Sierra, J.; Noel, C.; Dufour, L.; Ozier-Lafontaine, H.; Welcker, C.; Desfontaines, L. Mineral Nutrition and Growth of Tropical Maize as Affected by Soil Acidity. Plant Soil 2003, 252, 215–226. [Google Scholar] [CrossRef]
  33. Wang, X.T.; Lan, X.F.; An, W.L.; Xu, X.P.; Wang, W.Q. Effect of Simulated Acid Rain on Paddy Soil Bacterial Abundance and Diversity in Fuzhou Plain. China Environ. Sci. 2019, 39, 1237–1244. [Google Scholar]
  34. Wang, H.T.; Zheng, T.L.; Yang, X.R. Molecular Ecology Research Progress for Soil Denitrification and Research Status for its Incluencing Factors. J. Agro-Environ. Sci. 2013, 32, 1915–1924. [Google Scholar]
  35. Senbayram, M.; Budai, A.; Bol, R.; Chadwick, D.; Marton, L.; Gündogan, R.; Wu, D. Soil NO3− Level and O2 Availability are Key Factors in Controlling N2O Reduction to N2 following Long-Term Liming of an Acidic Sandy Soil. Soil Biol. Biochem. 2019, 132, 165–173. [Google Scholar] [CrossRef]
  36. Zhang, Z.; Wang, S.J.; Li, J.H.; Cao, R.; Chen, M.K.; Li, S.H. Response of Soil Readily Oxidizable Carbon to Community Succession of Xishuangbanna Tropical Forests. Acta Ecol. Sin. 2019, 39, 6257–6263. [Google Scholar]
  37. Zumft, W.G. Cell biology and Molecular Basis of Denitrification. Microbiol. Mol. Biol. Rev. 1997, 61, 533–616. [Google Scholar] [PubMed]
  38. Ding, B.J.; Zhang, H.; Luo, W.Q.; Sun, S.Y.; Cheng, F.; Li, Z.K. Nitrogen Loss Through Denitrification, Anammox and Feammox in a Paddy Soil. Sci. Total Environ. 2021, 773, 145601. [Google Scholar] [CrossRef] [PubMed]
  39. Zhang, J.; Lin, X.-G.; Yin, R. Advances in Functional Gene Diversity of Microorganism in Relation to Soil Nitrogen Cycling. Chin. J. Eco-Agric. 2009, 17, 1029–1034. [Google Scholar] [CrossRef]
  40. Marcel, M.M.K.; Hannah, K.M.; Boran, K. The Microbial Nitrogen-Cycling Network. Nat. Rev. Microbiol. 2018, 16, 263–276. [Google Scholar]
  41. Ding, H.; Wang, Y.S. Influence of Herbicides on Denitrification Loss of Nitrogen Fertilizer and Emission of N2O. Chin. Environ. Sci. 2004, 24, 596–599. [Google Scholar]
  42. Zhang, Y.S.; Ding, H.; Qin, S.J. Progress in the Studies of Nitrogen Denitrification and N2O Emission in Agroecosystem. Chin. Agric. Sci. Bull. 2010, 26, 253–259. [Google Scholar]
  43. Jones, C.M.; Hallin, S. Ecological and Evolutionary Factors Underlying Global and Local Assembly of Denitrifier Communities. ISME J. 2010, 4, 633–641. [Google Scholar] [CrossRef] [PubMed]
  44. Penton, C.R.; Johnson, T.A.; Quensen, J.F.; Iwai, S.; Cole, J.R.; Tiedje, J.M. Functional Fenes to Assess Nitrogen Cycling and Aromatic Hydrocarbon Degradation: Primers and Processing Matter. Front. Microbiol. 2013, 4, 279. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Braker, G.; Zhou, J.; Wu, L.; Devol, A.H.; Tiedje, J.M. Nitrite Reductase Genes (nirK and nirS) as Functional Markers to Investigate Diversity of Denitrifying Bacteria in Pacific Northwest Marine Sediment Communities. Appl. Environ. Microbiol. 2000, 66, 2096–2104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Tang, Y.; Zhang, X.; Li, D.; Wang, H.; Chen, F.; Fu, X.; Fang, X.; Sun, X.; Yu, G. Impacts of Nitrogen and Phosphorus Additions on the Abundance and Community Structure of Ammonia Oxidizers and Denitrifying Bacteria in Chinese Fir Plantations. Soil Biol. Biochem. 2016, 103, 284–293. [Google Scholar] [CrossRef]
  47. Chen, Z.; Luo, X.; Hu, R.; Wu, M.; Wu, J.; Wei, W. Impact of Long-Term Fertilization on the Composition of Denitrifier Communities Based on Nitrite Reductase Analyses in a Paddy Soil. Microb. Ecol. 2010, 60, 850–861. [Google Scholar] [CrossRef] [PubMed]
  48. Bárta, J.; Melichová, T.; Vaněk, D.; Picek, T.; Šantrůčková, H. Effect of pH and Dissolved Organic Matter on the Abundance of nirK and nirS Denitrifiers in Spruce Forest Soil. Biogeochemistry 2010, 101, 123–132. [Google Scholar] [CrossRef]
  49. Szukics, U.; Hackl, E.; Zechmeister-Boltenstern, S.; Sessitsch, A. Contrasting Response of Two Forest Soils to Nitrogen Input: Rapidly Altered NO and N2O Emissions and nirK Abundance. Biol. Fertil. Soils 2009, 45, 855–863. [Google Scholar] [CrossRef]
  50. Zhongjun, J.; Angel, R.; Veraart, A.; Daebeler, A.; Jia, Z.; Kim, S.Y.; Kerckhof, F.-M.; Boon, N.; Bodelier, P. Biotic Interactions in Microbial Communities as Modulators of Biogeochemical Processes: Methanotrophy as a Model System. Front. Microbiol. 2016, 7, 1285. [Google Scholar] [CrossRef] [Green Version]
  51. Huang, R.; Zeng, J.; Zhao, D.; Yong, B.; Yu, Z. Co-Association of Two nir Denitrifiers under the Influence of Emergent Macrophytes. Microb. Ecol. 2020, 80, 809–821. [Google Scholar] [CrossRef]
  52. Deng, Y.; Jiang, Y.-H.; Yang, Y.; He, Z.; Luo, F.; Zhou, J. Molecular Ecological Network Analyses. BMC Bioinform. 2012, 13, 113. [Google Scholar] [CrossRef] [Green Version]
  53. Barberán, A.; Bates, S.T.; Casamayor, E.O.; Fierer, N. Using Network Analysis to Explore Co-Occurrence Patterns in Soil Microbial Communities. ISME J. 2012, 6, 343–351. [Google Scholar] [CrossRef] [Green Version]
  54. Liu, Z.F.; Liu, L.L. Study on Increasing Resistance of Rice with Salicylic Acid. Guangdong Sci. Technol. News 2016, 2, 4–15. [Google Scholar]
  55. Luo, G.B. FAO World Soil Legend Classification System Revision. Soil Prog. 1988, 6, 22–27. [Google Scholar]
  56. Bao, S. Soil Agro-Chemistrical Analysis, 1st ed.; China Agriculture Press: Beijing, China, 2000; pp. 1–495. [Google Scholar]
  57. Di, H.J.; Cameron, K.; Podolyan, A.; Robinson, A. Effect of Soil Moisture Status and a Nitrification Inhibitor, Dicyandiamide, on Ammonia Oxidizer and Denitrifier Growth and Nitrous Oxide Emissions in a Grassland Soil. Soil Biol. Biochem. 2014, 73, 59–68. [Google Scholar] [CrossRef]
  58. Li, Z.G.; Luo, Y.M.; Teng, E. Soil and Environmental Microbiological Research, 1st ed.; Science Press: Beijing, China, 2008; pp. 50–100. [Google Scholar]
  59. Caporaso, J.C.; Christian, L.L.; Walters, W.A.; Berg-Lyons, D.; Huntley, J.; Fierer, N.; Owens, S.M.; Betley, J.; Fraser, L.; Bauer, M.; et al. Ultra-High-Throughput Microbial Community Analysis on the lllumina HiSeq and MiSeq Platforms. ISME J. 2012, 6, 1621–1624. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Edgar, R.C. UPARSE: Highly Accurate OTU Sequences from Microbial Amplicon Reads. Nat. Methods 2013, 10, 996–998. [Google Scholar] [CrossRef] [PubMed]
  61. Hou, S.; Ai, G.; Zhou, W.; Liang, G.; He, P. Structure and assembly cues for rhizospheric nirK- and nirS-type denitrifier communities in long-term fertilized soils. Soil Biol. Biochem. 2018, 119, 32–40. [Google Scholar] [CrossRef]
  62. Liu, C.H.; Wu, D.M.; Liu, Y.H.; Chen, H.; Shen, B.G.; Jang, Z.K.; Liu, X.F. Effects of Nitrogen Deposition on Soil Organic Carbon and Soil Microbial Communities in a Natural Castanopsis Carlesii Forest. For. Res. 2021, 34, 42–49. [Google Scholar]
  63. Campbell, B.J.; Polson, S.; Hanson, T.; Mack, M.C.; Schuur, E.A.G. The Effect of Nutrient Deposition on Bacterial Communities in Arctic Tundra Soil. Environ. Microbiol. 2010, 12, 1842–1854. [Google Scholar] [CrossRef]
  64. Herold, M.B.; Giles, M.E.; Alexander, C.J.; Baggs, E.M.; Daniell, T.J. Variable Response of nirK and nirS Containing Denitrifier Communities to Long-Term pH Manipulation and Cultivation. FEMS Microbiol. Lett. 2018, 365, fny035. [Google Scholar] [CrossRef]
  65. Ouyang, Y.; Evans, S.E.; Friesen, M.; Tiemann, L.K. Effect of Nitrogen Fertilization on the Abundance of Nitrogen Cycling Genes in Agricultural Soils: A Meta-Analysis of Field Studies. Soil Biol. Biochem. 2018, 127, 71–78. [Google Scholar] [CrossRef]
  66. Yang, Y.; Zhao, J.; Jiang, Y.; Hu, Y.; Zhang, M.; Zeng, Z. Response of Bacteria Harboring nirS and nirK Genes to Different N Fertilization Rates in an Alkaline Northern Chinese Soil. Eur. J. Soil Biol. 2017, 82, 1–9. [Google Scholar] [CrossRef]
  67. Xie, Z.; LE Roux, X.; Wang, C.; Gu, Z.; An, M.; Nan, H.; Chen, B.; Li, F.; Liu, Y.; Du, G.; et al. Identifying Response Groups of Soil Nitrifiers and Denitrifiers to Grazing and Associated Soil Environmental Drivers in Tibetan Alpine Meadows. Soil Biol. Biochem. 2014, 77, 89–99. [Google Scholar] [CrossRef]
  68. Fan, K.; Delgado-Baquerizo, M.; Guo, X.; Wang, D.; Wu, Y.; Zhu, M.; Yu, W.; Yao, H.; Zhu, Y.-G.; Chu, H. Suppressed N Fixation and Diazotrophs after Four Decades of Fertilization. Microbiome 2019, 7, 1–10. [Google Scholar] [CrossRef] [PubMed]
  69. Zhu, J.; Liu, H.; Wu, B.K.; Yu, F.; Liu, Z.Y.; Jin, T. Effects of Integrated Rice-Crayfish Farming System on Community Structure and Diversity of nirK Denitrification Microbe in Paddy Soils. Chin. J. Eco-Agric. 2018, 26, 1324–1332. [Google Scholar]
  70. Bremer, C.; Braker, G.; Matthies, D.; Reuter, A.; Engels, C.; Conrad, R. Impact of Plant Functional Group, Plant Species, and Sampling Time on the Composition of nirK-type Denitrifier Communities in Soil. Appl. Environ. Microbiol. 2007, 73, 6876–6884. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Wang, M.; Wu, J.; Zhou, T.; Liang, Y.; Zheng, L.X.; Sun, Y.X. Effects of Copper and Florfenicol on nirS- and nirK-type Denitrifier Communities and Related Antibiotic Resistance in Vegetable Soils. Ecotoxicol. Environ. Saf. 2021, 213, 112011. [Google Scholar] [CrossRef] [PubMed]
  72. Qin, Y.; Wei, Q.; Zhang, Y.; Li, H.; Jiang, Y.; Zheng, J. Nitrogen Removal from Ammonium- and Sulfate-Rich Wastewater in an Upflow Anaerobic Sludge Bed Reactor: Performance and Microbial Community Structure. Ecotoxicology 2021, 30, 1719–1730. [Google Scholar] [CrossRef] [PubMed]
  73. Enwall, K.; Philippot, L.; Hallin, S. Activity and Composition of the Denitrifying Bacterial Community Respond Differently to Long-Term Fertilization. Appl. Environ. Microbiol. 2005, 71, 8335–8343. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Henderson, S.L.; Dandie, C.E.; Patten, C.L.; Zebarth, B.J.; Burton, D.L.; Trevors, J.T.; Goyer, C. Changes in Denitrifier Abundance, Denitrification Gene mRNA Levels, Nitrous Oxide Emissions, and Denitrification in Anoxic Soil Microcosms Amended with Glucose and Plant Residues. Appl. Environ. Microbiol. 2010, 76, 2155–2164. [Google Scholar] [CrossRef] [Green Version]
  75. Yin, C.; Fan, F.; Song, A.; Cui, P.; Li, T.; Liang, Y. Denitrification Potential under Different Fertilization Regimes is Closely Coupled with Changes in the Denitrifying Community in a Black Soil. Appl. Microbiol. Biotechnol. 2015, 99, 5719–5729. [Google Scholar] [CrossRef] [PubMed]
  76. Liu, Y.L.; Zhang, B.; Hu, F.; Qiao, J.; Zhang, W.J. Carbon and Nitrogen Mineralization of Paddy Soils as Affected by Wetting and Drying. Soils 2008, 40, 554–560. [Google Scholar]
  77. Tang, H.; Shen, J.L.; Zhang, Y.Z.; Liu, J.Y.; Wang, C. Effect of Rice Straw Incorporation and Water Management on Soil Microbial Biomass Carbon, Nitrogen and Dissolved Organic Carbon, Nitrogen in a Rice Paddy Field. J. Soil Water Conserv. 2013, 27, 240–246. [Google Scholar]
  78. Karasu, A.; Dogan, M.Z. The Effect of Bacterial Inoculation and Different Nitrogen Dose on Yield and Yield Components of Some Chickpea Genotypes (Cicer arietinum L.). Afr. J. Biotechnol. 2009, 8, 59–64. [Google Scholar]
  79. Rinklebe, J.; Langer, U. Microbial Diversity in Three Flood Plain Soils at the Elbe River (Germany). Soil Biol. Biochem. 2006, 38, 2144–2151. [Google Scholar] [CrossRef]
  80. Chen, Z.; Hou, H.; Zheng, Y.; Qin, H.; Zhu, Y.; Wu, J.; Wei, W. Influence of Fertilization Regimes on a nosZ-Containing Denitrifying Community in a Rice Paddy Soil. J. Sci. Food Agr. 2012, 92, 1064–1072. [Google Scholar] [CrossRef] [PubMed]
  81. Tao, R.; Wakelin, S.; Liang, Y.; Hu, B.; Chu, G. Nitrous Oxide Emission and Denitrifier Communities in Drip-Irrigated Calcareous Soil as Affected by Chemical and Organic Fertilizers. Sci. Total Environ. 2018, 612, 739–749. [Google Scholar] [CrossRef]
  82. Delmont, T.O.; Francioli, D.; Jacquesson, S.; Laoudi, S.; Mathieu, A.; Nesme, J.; Ceccherini, M.T.; Nannipieri, P.; Simonet, P.; Vogel, T.M. Microbial Community Development and Unseen Diversity Recovery in Inoculated Sterile Soil. Biol. Fertil. Soils 2014, 50, 1069–1076. [Google Scholar] [CrossRef]
  83. Chen, C.; Zhang, J.; Lu, M.; Qin, C.; Chen, Y.; Yang, L.; Huang, Q.; Wang, J.; Shen, Z.; Shen, Q. Microbial Communities of an Arable Soil Treated For 8 Years with Organic and Inorganic Fertilizers. Biol. Fertil. Soils 2016, 52, 455–467. [Google Scholar] [CrossRef]
  84. Fang, L.-N.; Yang, X.-D.; Du, J. Effects of Land Use Pattern on Soil Microbial Biomass Carbon in Xishuangbanna. J. Appl. Ecol. 2011, 22, 837–844. [Google Scholar]
  85. Zhou, Z.; Wang, C.; Luo, Y. Meta-Analysis of the Impacts of Global Change Factors on Soil Microbial Diversity and Functionality. Nat. Commun. 2020, 11, 3072. [Google Scholar] [CrossRef] [PubMed]
  86. Enwall, K.; Throbäck, I.N.; Stenberg, M.; Söderström, M.; Hallin, S. Soil Resources Influence Spatial Patterns of Denitrifying Communities at Scales Compatible with Land Management. Appl. Environ. Microb. 2010, 76, 2243–2250. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Biology 10 01191 g001 550
Figure 1. Functional genes for nitrification and denitrification. Abbreviations: NO: nitric oxide; NO2−: nitrite ion; NO3−: nitrate ion; NH4+: ammonium ion; NH2OH: hydroxylamine; N2: nitrogen; DNRA: dissimilatory nitrate reduction to ammonium. Dotted lines indicate possible reaction process.
Figure 1. Functional genes for nitrification and denitrification. Abbreviations: NO: nitric oxide; NO2−: nitrite ion; NO3−: nitrate ion; NH4+: ammonium ion; NH2OH: hydroxylamine; N2: nitrogen; DNRA: dissimilatory nitrate reduction to ammonium. Dotted lines indicate possible reaction process.
Biology 10 01191 g001
Biology 10 01191 g002 550
Figure 2. Venn diagram of sample OTUs quantity. Different numbers after the treatments indicate the treatment duration. Lower-case letters (a and b) represent different genes. Treatments: TN/S: NaNO3 + Na2SO4; TNaCl: NaCl; CK: control. There were 110 OTUs shared by all groups (a). There were 259 OTUs common to all groups (b).
Figure 2. Venn diagram of sample OTUs quantity. Different numbers after the treatments indicate the treatment duration. Lower-case letters (a and b) represent different genes. Treatments: TN/S: NaNO3 + Na2SO4; TNaCl: NaCl; CK: control. There were 110 OTUs shared by all groups (a). There were 259 OTUs common to all groups (b).
Biology 10 01191 g002
Biology 10 01191 g003 550
Figure 3. The number of OTUs of the sample at different classification levels. Lower-case letters (a and b) represent different genes. Treatments: TN/S: NaNO3 + Na2SO4; TNaCl: NaCl; CK: control. In terms of the nirK gene, on the 7th day, the numbers of OTUs at the phylum, class, and order levels under each treatment were the same, while the numbers of OTUs at the family, genus, and species levels were lower than those of the control, while after 30, 60 days of incubation, the numbers of OTUs in the soil treated with TN/S at the fam-ily, genus, and species levels were higher than those of the control (a). In terms of the nirS gene, we found that there was a non-significant difference in the number of OTUs between the treatments at the phylum and class levels. However, a significant difference was observed, that is, compared with CK, TN/S increased the numbers of OTUs at the order, family, genus, and species levels in all three time periods (b).
Figure 3. The number of OTUs of the sample at different classification levels. Lower-case letters (a and b) represent different genes. Treatments: TN/S: NaNO3 + Na2SO4; TNaCl: NaCl; CK: control. In terms of the nirK gene, on the 7th day, the numbers of OTUs at the phylum, class, and order levels under each treatment were the same, while the numbers of OTUs at the family, genus, and species levels were lower than those of the control, while after 30, 60 days of incubation, the numbers of OTUs in the soil treated with TN/S at the fam-ily, genus, and species levels were higher than those of the control (a). In terms of the nirS gene, we found that there was a non-significant difference in the number of OTUs between the treatments at the phylum and class levels. However, a significant difference was observed, that is, compared with CK, TN/S increased the numbers of OTUs at the order, family, genus, and species levels in all three time periods (b).
Biology 10 01191 g003
Biology 10 01191 g004 550
Figure 4. Non-metric multidimensional scaling analysis of nirK- (a) and nirS-type (b) denitrifying communities under different treatments. Treatments: TN/S: NaNO3 + Na2SO4; TNaCl: NaCl; CK: control.
Figure 4. Non-metric multidimensional scaling analysis of nirK- (a) and nirS-type (b) denitrifying communities under different treatments. Treatments: TN/S: NaNO3 + Na2SO4; TNaCl: NaCl; CK: control.
Biology 10 01191 g004
Biology 10 01191 g005 550
Figure 5. The community compositions of nirK- (a) and nirS-type (b) denitrifying communities under different treatments. Treatments: TN/S: NaNO3 + Na2SO4; TNaCl: NaCl; CK: control. At the genus level, a large number of nirK reads (0.3–20.8%) and nirS reads (1.4–48.8%) are defined as “unidentified” OTUs, suggesting that there remains a large degree of nirK and nirS diversity, which is not reflected in the current genome database [61].
Figure 5. The community compositions of nirK- (a) and nirS-type (b) denitrifying communities under different treatments. Treatments: TN/S: NaNO3 + Na2SO4; TNaCl: NaCl; CK: control. At the genus level, a large number of nirK reads (0.3–20.8%) and nirS reads (1.4–48.8%) are defined as “unidentified” OTUs, suggesting that there remains a large degree of nirK and nirS diversity, which is not reflected in the current genome database [61].
Biology 10 01191 g005
Biology 10 01191 g006 550
Figure 6. Redundancy analysis (RDA) of soil properties and functional gene abundances from different treatments. Lower-case letters (a and b) represent different genes. Treatments: TN/S: NaNO3 + Na2SO4; TNaCl: NaCl; CK: control. In the case of the nirK-type denitrifying communities, pH and NO3−-N were signifi-cantly positively correlated with axis 1, while nitrite reductase was significantly nega-tively correlated (a). Similarly, in the nirS-type denitrifying communities, TN and NH4+-N were significantly positively correlated with axis 1, while nitrite reductase was significantly negatively correlated (b).
Figure 6. Redundancy analysis (RDA) of soil properties and functional gene abundances from different treatments. Lower-case letters (a and b) represent different genes. Treatments: TN/S: NaNO3 + Na2SO4; TNaCl: NaCl; CK: control. In the case of the nirK-type denitrifying communities, pH and NO3−-N were signifi-cantly positively correlated with axis 1, while nitrite reductase was significantly nega-tively correlated (a). Similarly, in the nirS-type denitrifying communities, TN and NH4+-N were significantly positively correlated with axis 1, while nitrite reductase was significantly negatively correlated (b).
Biology 10 01191 g006
Biology 10 01191 g007 550
Figure 7. Correlation heatmap of soil physicochemical factors and nirK- (a) and nirS-type (b) denitrifying communities in different periods. The R values are shown in different colors, and *** indicates p ≤ 0.001.
Figure 7. Correlation heatmap of soil physicochemical factors and nirK- (a) and nirS-type (b) denitrifying communities in different periods. The R values are shown in different colors, and *** indicates p ≤ 0.001.
Biology 10 01191 g007
Table
Table 1. Physicochemical properties of soil samples.
Table 1. Physicochemical properties of soil samples.
DaysTreatmentspHEC (μS·cm−1)TN (g·kg−1)TOC (g·kg−1)NO3−-N (mg·kg−1)NH4+-N (mg·kg−1)
7TN/S5.9 + 0.0 b115 + 8.9 a1.14 + 0.28 a28 + 0.6 a28 + 0.6 a11 + 0.2 a
TNaCl6.2 + 0.0 a82 + 4.6 b1.03 + 0.02 b20 + 0.5 b17 + 0.1 b7 + 0.3 b
CK6.2 + 0.1 a70 + 3.5 b1.05 + 0.04 b20 + 0.8 b16 + 0.4 c7 + 0.2 b
30TN/S6.8 + 0.1 a128 + 7.5 a1.09 + 0.07 a22 + 0.3 a44 + 0.7 a14 + 0.6 a
TNaCl6.7 + 0.1 a137 + 3.6 a1.11 + 0.08 a22 + 0.7 a45 + 0.9 a14 + 0.5 a
CK6.8 + 0.1 a130 + 7.5 a1.15 + 0.09 a22 + 0.4 a45 + 0.6 a15 + 0.7 a
60TN/S6.3 + 0.0 b266 + 2.1 a1.14 + 0.02 a23 + 0.6 a51 + 0.8 a6 + 0.1 a
TNaCl6.8 + 0.0 a220 + 8.5 b0.90 + 0.02 b17 + 0.4 b48 + 0.6 b5 + 0.1 b
CK6.9 + 0.1 a224 + 8.1 b0.96 + 0.03 c17 + 0.1 b48 + 0.5 b5 + 0.0 b
All values are presented as mean ± standard error (n = 5); different letters in the same column under the same period indicate significant differences between treatments (p < 0.05). Treatments: TN/S: NaNO3 + Na2SO4; TNaCl: NaCl; CK: control.
Table
Table 2. Alpha diversity index of soil nirK- and nirS-type denitrifying microorganisms.
Table 2. Alpha diversity index of soil nirK- and nirS-type denitrifying microorganisms.
Gene TypesDaysSampleShannonChao 1
nirK7TN/S2.77 ± 0.01 a246 ± 1.8 a
TNaCl2.56 ± 0.00 c219 ± 4.4 c
CK2.73 ± 0.00 b237 ± 5.4 b
30TN/S3.38 ± 0.01 a274 ± 8.5 a
TNaCl3.07 ± 0.01 b242 ± 10.8 b
CK3.03 ± 0.01 c238 ± 8.4 b
60TN/S3.07 ± 0.01 a262 ± 3.0 a
TNaCl2.79 ± 0.01 c254 ± 5.3 b
CK3.03 ± 0.01 b228 ± 2.7 c
nirS7TN/S4.49 ± 0.01 b619 ± 4.8 a
TNaCl4.34 ± 0.00 c620 ± 9.6 a
CK4.51 ± 0.01 a620 ± 12.4 a
30TN/S4.60 ± 0.01 a619 ± 9.2 a
TNaCl4.61 ± 0.01 a623 ± 13.1 a
CK4.55 ± 0.00 b619 ± 7.2 a
60TN/S4.53 ± 0.01 a602 ± 11.2 b
TNaCl4.54 ± 0.01 a596 ± 7.3 b
CK4.51 ± 0.00 b615 ± 7.7 a
All values are presented as mean ± standard error (n = 5); different letters in the same column under the same gene and the same period indicate significant differences between treatments (p < 0.05). Treatments: TN/S: NaNO3 + Na2SO4; TNaCl: NaCl; CK: control.
Table
Table 3. Enzyme activities of soil samples.
Table 3. Enzyme activities of soil samples.
EnzymesTreatments7 days30 days60 days
Nitrate reductaseTN/S0.54 + 0.02 a0.54 + 0.00 a0.54 + 0.01 a
TNaCl0.49 + 0.01 b0.51 + 0.00 b0.51 + 0.01 b
CK0.50 + 0.02 b0.51 + 0.01 b0.51 + 0.01 b
Nitrite reductaseTN/S1.53 + 0.00 a1.53 + 0.01 a1.39 + 0.01 a
TNaCl1.14 + 0.01 b1.12 + 0.00 b1.03 + 0.00 b
CK1.11 + 0.01 c1.12 + 0.01 b1.01 + 0.01 c
All values are presented as mean ± standard error (n = 5), enzyme activity is expressed in U·g−1, different letters in the same column under the same enzyme indicate significant differences between treatments (p < 0.05). Treatments: TN/S: NaNO3 + Na2SO4; TNaCl: NaCl; CK: control.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, X.; Wei, H.; Zhang, J. Nitrogen and Sulfur Additions Improved the Diversity of nirK- and nirS-Type Denitrifying Bacterial Communities of Farmland Soil. Biology 2021, 10, 1191. https://doi.org/10.3390/biology10111191

AMA Style

Chen X, Wei H, Zhang J. Nitrogen and Sulfur Additions Improved the Diversity of nirK- and nirS-Type Denitrifying Bacterial Communities of Farmland Soil. Biology. 2021; 10(11):1191. https://doi.org/10.3390/biology10111191

Chicago/Turabian Style

Chen, Xuan, Hui Wei, and Jiaen Zhang. 2021. "Nitrogen and Sulfur Additions Improved the Diversity of nirK- and nirS-Type Denitrifying Bacterial Communities of Farmland Soil" Biology 10, no. 11: 1191. https://doi.org/10.3390/biology10111191

APA Style

Chen, X., Wei, H., & Zhang, J. (2021). Nitrogen and Sulfur Additions Improved the Diversity of nirK- and nirS-Type Denitrifying Bacterial Communities of Farmland Soil. Biology, 10(11), 1191. https://doi.org/10.3390/biology10111191

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