- 1College of Land and Environment, Shenyang Agricultural University, Shenyang, China
- 2National Engineering Laboratory for Efficient Utilization of Soil and Fertilizer Resources, Tai’an, China
- 3Northeast Scientific Observation Station of Corn Nutrition and Fertilization of Ministry of Agriculture, Shenyang, China
- 4Department of Biotechnology, University of Sargodha, Sargodha, Pakistan
The objective of this study was to find out changes in ammonia oxidation microorganisms with respect to fertilizer as investigated by real-time polymerase chain reaction and high-throughput sequencing. The treatments included control (CK); chemical fertilizer nitrogen low (N) and high (N2); nitrogen and phosphorus (NP); nitrogen phosphorus and potassium (NPK) and organic manure fertilizer (M); MN; MN2; MNPK. The results showed that long-term fertilization influenced soil fertility and affected the abundance and community of ammonia-oxidizing microorganisms by changing the physical and chemical properties of the soil. The abundance and community structure of ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) was influenced by soil organic carbon, total nitrogen, total soil phosphorus, available phosphorus, available potassium, and soil nitrate. Soil environmental factors affected the nitrification potential by affecting the structure of ammonia-oxidizing microorganisms; specific and rare AOA and AOB rather than the whole AOA or AOB community played dominant role in nitrification.
Introduction
Agricultural nitrogen cycling is a crucial part in global nitrogen cycling, and ammoxidation is one of critical process of agricultural nitrogen cycling, which turns ammonia to nitrite and then to nitrate. Ammonia-oxidizing organisms, limiting the rate, play an important role in the first step of nitrification. As Könneke et al. (2005) found that archaea could oxidize ammonia to nitrite, researchers studied the comparison between ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) about which one was more important under different environmental conditions (Adair and Schwartz, 2008; Chen et al., 2008, 2015; Li et al., 2010; Tao et al., 2017). However, it is controversial whether AOA or AOB play a dominant role under different environmental conditions (Li et al., 2010; Liu et al., 2015; Gan et al., 2016), which may be due to that ammonia-oxidizing microorganisms community structure were affected by physiochemical properties of soil (Li et al., 2010; Huang et al., 2013). In a previous study, AOA might be dominant within nitrogen cycling under extreme environmental conditions (Erguder et al., 2009); AOA rather than AOB could adapt to oxygenated/hypoxic water conditions (Liu et al., 2015). However, AOB rather than AOA dominated in potential nitrification in other environment (Dang et al., 2018). Different NH4+-N concentration and different levels of total carbon and total nitrogen in soils influenced the composition and AOA and AOB abundance in various environmental conditions (Huang et al., 2013; Li and Gu, 2013; Wang et al., 2013).
Long-term fertilization, a good way to investigate agricultural soil, changed physicochemical properties and microbial community structure, such as soil pH, dissolved organic carbon (DOC), available phosphorus, etc. (Luo et al., 2015). Crop rotation could also change the physiochemical properties of soil and microbial community structure (Belay et al., 2002; Bhattacharyya et al., 2010). Whether fertilization changed the ammonia-oxidizing microorganisms’ community structure and quantities, which ammonia-oxidizing microorganisms play a dominant role in the nitrification process under the condition of long-term rotation and fertilization is still unclear. In this study, we evaluated the effects of fertilization on the changing of AOA and AOB because of a long-term (38 years) rotation and fertilization experiment. Our objectives were (i) to determine whether AOA and AOB had different responses to variations under the condition of long-term fertilization and rotation, (ii) to identify which taxa play a more crucial role during nitrification process, and (iii) to assess the relative contribution of those variables to changes in AOA and AOB activities and function.
Materials and Methods
Experimental Site and Soil Sampling
This 38 years fertilization experiment was established in 1979; the field soil properties and block design were described by Hua et al. (2020). The field experiment was conducted under a rotation of maize–maize–soybean. The crop-growing season started at the end of April and ended at the end of September. In this study, we have chosen nine treatments as follows: CK (no fertilizer), N (mineral nitrogen fertilizer), N2 (high mineral nitrogen fertilizer), NP (N and phosphate fertilizer), NPK (N, P, and potassium fertilizer), and M (manure), MN, MN2, MNPK. N was applied in the form of urea, whereas P was calcium superphosphate, and K was potassium sulfate, and the manure fertilizer used was pig manure; the specific application rates are detailed in Supplementary Table S1. Soil samples were collected from surface soil layer (0–20 cm) and plow pan soil layer (20–40 cm) on each plot on October 7, 2016. Three soil cores (5 cm diameter) were taken from each layer of soil in each plot, respectively. All samples were shifted through a 2.0 mm sieve. Some part of fresh soil samples was used for DNA extraction, whereas some part was used to determine the soil moisture content, and the rest was air-dried and ground for physiochemical properties determination.
Chemical Analysis
Soil properties determined was soil pH, DOC, soil organic carbon (SOC), total nitrogen (TN), total phosphorus (TP), total potassium (TK), and available phosphorus (AP) and potassium (AK), as described earlier (Luo et al., 2015). The fresh soil samples were used to determine soil ammonium (NH4+-N) and nitrate (NO3–-N), extracted with 1 M KCl and determined by a continuous flow analyzer (AutoAnalyzer 3; SEAL Analytical Ltd., Germany). Soil potential nitrification rate (PNR) was determined using the method of Kandeler (1995).
DNA Extraction and Quantitative Polymerase Chain Reaction Amplification
DNA was extracted from soil (0.250 g) using the PowerSoilTM DNA Isolation Kit (MoBio Laboratories Inc., United States). Purified DNA extracts were stored at −20°C until further determination. Amplification of archaeal and bacterial amoA genes was performed using the primer, which is listed in Supplementary Table S2. The reaction mixture (20 μL) contained 1 μL of DNA template (20–40 ng per reaction), 0.8 μL of forward and reverse primers (10 μM), respectively, and 10 μL of TB Green Premix Ex Taq (2X) 10 μL (Takara Biomedical Technology Co. Ltd., Beijing, China). The quantitative polymerase chain reaction (qPCR) amplification procedure of amoA gene was performed as prime denaturation at 95°C for 30 s, followed with 40 cycles of denaturation for 5 s at 95°C, annealing for 30 s at 55°C, and extension for 30 s at 72°C. The qPCR amplification procedure of arch-amoA gene was performed as prime denaturation at 94°C for 30 s, followed with 45 cycles of denaturation for 30 s at 94°C, annealing for 60 s at 55°C, and extension for 60 s at 72°C.
High-Throughput Sequencing
The pure PCR product equimolar and paired end was sequenced according to the standard protocols by Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China) (Zhang et al., 2020). The raw amoA gene sequencing reads were demultiplexed, quality-filtered by fast p version 0.20.0 (Chen et al., 2018), and merged by FLASH version 1.2.7 (Magoč and Salzberg, 2011) with the same following criteria as described by Zhang et al. (2020).
UPARSE version 7.1 clustering was used to identify and remove chimeric sequences of operation classification unit (OTUs) with 97% similarity truncation (Edgar, 2013). RDP version 2.2 was used for classification analysis of amoA database, and the confidence threshold was 0.7 (Wang et al., 2007).
Statistical Analysis
Redundancy analysis (RDA) (CANOCO 4.5 for Windows) was used to evaluate the relationships among high-throughput sequencing results, quantity PCR results, and physicochemical properties of fertilization treatments. The measured data were used to analyze the variation, using IBM SPSS Statistics 19.0 for Windows (IBM, Inc., Armonk, NY, United States). Analysis of variance and non-metric multidimensional scaling (NMDS) analysis (calculated by Qiime, and analysis and mapping by Vegan software package for R) were used to analyze the environmental factors and microbes index. A probability level of 5% was adopted for accepting or rejecting null hypotheses.
Results
Soil Physicochemical Properties
Long-term fertilization prominently changed the brown soil properties from both physical and chemical aspects (Table 1). Soil pH was varied from 5.03 to 6.95 and 5.78 to 6.86; SOC contents varied from 13.91 to 23.32 g ⋅ kg–1 and 10.42 to 20.51 g ⋅ kg–1; AP contents varied from 1.65 to 136.75 mg ⋅ kg–1 and 1.45 to 117.90 mg ⋅ kg–1; and NO3–-N contents varied from 2.65 to 27.32 mg ⋅ kg–1 and 1.50 to 19.80 mg ⋅ kg–1. PNR was varied from 0.86 to 5.06 mg ⋅ kg–1 ⋅ d–1 and 0.93 to 3.99 mg ⋅ kg–1 ⋅ d–1; DOC contents varied from 21.83 to 36.04 mg ⋅ kg–1 and 3.72 to 14.05 mg ⋅ kg–1 in surface arable soil layer (0–20 cm) and plow pan soil layer (20–40 cm), respectively. The application of manure significantly (P < 0.05) increased the soil pH, PNR, and the stocks of SOC, TN, TP, AP, AK, and NO3–-N of soil in both surface arable soil layer and plow pan soil layer. It had no remarkable difference (P > 0.05) in the stocks of TK and NH4+-N in all fertilization treatments.
Table 1. Physicochemical properties of soils in different fertilization treatments after 38 years of long-term fertilization.
In addition to organic fertilizer treatments, contents of SOC in chemical phosphorus fertilizer treatments (NP, NPK) were significantly higher (P < 0.05) than that in N and N2 treatments. The content of PNR in organic fertilizer treatments was significantly higher (P < 0.05) than that in chemical fertilizer treatments, and that in surface arable soil was significantly higher (P < 0.05) than that in plow pan soil. It was highest in M treatment in both soil layers. The stock of DOC in surface arable soil in the M treatment was the highest in all treatments, and contents of DOC in plow pan soil in N, MN, and MN2 treatments were the highest in all treatments.
Abundance and Community Structures of AOB and AOA in Soil
The aomA gene (AOB) quantified by qPCR ranging from 4.68 × 105 to 7.28 × 106 copy numbers per gram of dry soil in surface arable soil layer (Figure 1A) and 7.44 × 105 to 6.1 × 106 copy numbers per gram of dry soil in plow pan soil layer (Figure 1B) in this experiment. The abundance of amoA gene in both soil layers in organic fertilizer treatments was prominently higher (P < 0.05) than that in chemical fertilizer treatments. Except for N2 and MNPK treatments, the abundance of amoA gene in surface arable soil layer was higher than that in plow pan soil layer.
Figure 1. The abundance and community structure in species level of ammonia-oxidizing bacteria in different treatments in surface arable soil (A) and plow pan soil (B), respectively. The mean value ± SD (LSD method, and p < 0.05, n = 3).
Community structures of AOB were similar in both of two soil layers under long-term fertilization (Figures 1A,B). The abundance of uncultured bacterium MF323203.1 (MNPK 23.49 and 50.90%, MN 24.99 and 10.21%, MN2 30.86 and 8.91%, M 11.52 and 5.92% in surface arable soil layer and plow pan soil layer, respectively), uncultured bacterium MF323620.1 (MNPK 22.25 and 81.09%, MN 14.22 and 3.88%, MN2 28.83 and 9.14%, M 13.95 and 3.40% in surface arable soil layer and plow pan soil layer, respectively), and Nitrosospira PJA1 (MNPK 19.00 and 74.35%, MN 22.97 and 9.78%, MN2 37.87 and 10.68%, M 3.27 and 0.76% in surface arable soil layer and plow pan soil layer, respectively) in organic fertilizer treatments was prominently higher (P < 0.05) than that in other treatments. The abundance of unclassified Nitrosomonadaceae in NP (27.70 and 7.67% in surface arable soil layer and plow pan soil layer, respectively), NPK (46.70 and 10.24% in surface arable soil layer and plow pan soil layer, respectively), and Nitrosospira.L115 in NP (27.34% in surface arable soil layer) and NPK (25.24% in surface arable soil layer) treatments were significantly higher (P < 0.05) than that in other treatments. The abundance of Nitrosomonadaceae (40.78 and 40.02% in surface arable soil layer and plow pan soil layer, respectively) was the highest in N2 treatment. The abundance of uncultured Nitrosomonadales bacterium was the highest in M treatment (87.62 and 62.64% in surface arable soil layer and plow pan soil layer, respectively). The abundance of Nitrosospira.Ka3, uncultured AOB and Nitrosospira.Nsp12 in fertilizer treatments were prominently higher (P < 0.05) than that in CK treatment (Nitrosospira.Ka3 0.62 and 5.82%, uncultured AOB 0.73 and 0.07% in surface arable soil layer and plow pan soil layer, respectively). Unclassified Nitrosospira in surface arable soil layer was different from that in plow pan soil layer. The quantities of unclassified Nitrosospira were dominantly higher in NP (11.52%), NPK (20.01%) treatments, and organic fertilizer treatments (MNPK 23.11%, MN 11.86%, MN2 25.40%, M 3.83%) than that in other treatments in surface arable soil, whereas the quantities of unclassified Nitrosospira were significantly higher (P < 0.05) in organic fertilizer treatments (MNPK 69.69%, MN 5.09%, MN2 11.79%, M 5.68%) than that in chemical fertilizer treatments in plow pan soil.
The arch-amoA gene (AOA) quantified by qPCR ranging from 9.82 × 105 to 6.49 × 106 copy numbers per gram of dry soil in surface arable soil layer (Figure 2A) and 2.07 × 106–9.09 × 106 copy numbers per gram of dry soil in plow pan soil layer (Figure 2B). The abundance of arch-amoA gene of soil in MN2 treatment was significantly higher (P < 0.05) than that in other treatments in surface arable soil layer, and the abundance of arch-amoA gene of soil in M treatment was significantly higher (P < 0.05) than that in other treatments in plow pan soil layer. There was no prominent difference (P > 0.05) in arch-amoA gene abundance in surface arable soil layer of N, NPK, MN, and MNPK treatments. The abundance of arch-amoA gene in NP, N2, and CK treatments were significantly lower (P < 0.05) than that in other treatments in surface arable soil. The abundance of arch-amoA gene in M, MNPK, N2, and N treatments were significantly higher than CK treatment, and CK treatment had no significant difference with MN and MN2 treatments in plow pan soil. Except for MN2 and NPK treatments, the abundance of arch-amoA gene in surface arable soil layer was significantly higher (P < 0.05) than that in plow pan soil layer.
Figure 2. The abundance and community structure in species level of ammonia-oxidizing archaea in different treatments in surface arable soil (A) and plow pan soil (B), respectively. The mean value ± SD (LSD method, and p < 0.05, n = 3).
Community structures of AOA were similar in both of two soil layers under long-term fertilization (Figures 2A,B). The abundance of Crenarchaeota in M (0.37 and 0.56% in surface arable soil layer and plow pan soil layer, respectively) and CK (1.28 and 0.93% in surface arable soil layer and plow pan soil layer, respectively) treatments was dominantly lower (P < 0.05) than that in other treatments; MN2 treatment (27.69%) was the highest than other treatments in surface arable soil, and N2 treatment (31.58%) was the highest treatment than other treatments in plow pan soil. It indicated that excessive nitrogen was good to Crenarchaeota, and the increment of Crenarchaeota increased with the increase of chemical nitrogen fertilizer; in opposite, organic manure alone inhibited Crenarchaeota. The abundance of uncultured AOA in NPK (22.22 and 9.11% in surface arable soil layer and plow pan soil layer, respectively) and M (19.31 and 34.69% in surface arable soil layer and plow pan soil layer, respectively) treatments was dominantly higher (P < 0.05) than that in other treatments. The abundance of Nitrososphaera.JG1 (MNPK 13.75 and 27.33%, MN 18.69 and 27.70%, MN2 38.52 and 10.08%, M 22.21 and 25.37% in surface arable soil layer and plow pan soil layer, respectively) and Thaumarchaeota (MNPK 20.50 and 17.97%, MN 20.44 and 12.88%, MN2 12.68 and 8.18%, and M 44.44 and 56.92% in surface arable soil layer and plow pan soil layer, respectively) in organic fertilizer treatments was significantly higher (P < 0.05) than that in chemical fertilizer treatments. The abundance of uncultured ammonia-oxidizing Crenarchaeote and uncultured Crenarchaeote were different between two soil layers. The abundance of uncultured ammonia-oxidizing Crenarchaeote in organic manure treatments (MNPK 9.41%, MN 13.88%, MN2 22.87%, M 15.64%) was significantly higher (P < 0.05) than that in other treatments in surface arable soil, whereas that in M (15.21%) and CK (25.67%) treatments was significantly higher (P < 0.05) as compared to other treatments in plow pan soil. The abundance of uncultured Crenarchaeota in NPK (19.75%), N (15.89%), MN (12.37%), and MN2 (17.13%) was significantly higher (P < 0.05) than that in other treatments in surface arable soil, whereas that in CK treatment (25.95%) was significantly higher (P < 0.05) with respect to other treatments in plow pan soil.
Community Structures of AOB and AOA in Soil Under Relative Abundance
The relative abundance of AOB community structures was different between surface arable soil with plow pan soil under long-term fertilization (Supplementary Figures S1A,B). The relative quantity of Nitrosospira.Ka3 in CK treatment was prominently lower than that in other treatments in surface arable soil (Supplementary Figure S1A). The relative abundance of unclassified Nitrosomonadaceae, unclassified environmental samples Nitrosomonadaceae, and Nitrosospira.L115 in chemical treatments was prominently higher than that in other treatments, and the relative abundance of unclassified environmental samples Nitrosomonadaceae in N and N2 treatments was prominently higher as compared to other treatments. The relative abundance of uncultured bacterium MF323203.1 in MNPK, MN, and MN2 treatments was significantly higher than that in other treatments. It indicated that application of organic in combination with chemical fertilizer could increase the relative quantity of uncultured AOB. The relative abundance of unclassified Nitrosomonadaceae, unclassified environmental samples Nitrosomonadaceae in CK treatment, was significantly higher than that in other treatments in plow pan soil (Supplementary Figure S1B). The relative abundance of Nitrosospira.Ka3 in chemical fertilizer treatments was significantly higher than that in other treatments. The relative abundance of uncultured Nitrosomonadales bacterium and Nitrosospira.Nsp12 in M treatment was prominently higher than that in other treatments. The relative abundance of unclassified Nitrosospira, Nitrosospira.PJA1 in N, N2, M, and CK treatments was significantly lower among all other treatments. The relative abundance of AOA community structures was similar in both soil layers under long-term fertilization (Supplementary Figures S2A,B). The relative abundance of Crenarchaeota in M and CK treatments was dominantly lower than that in other treatments. Crenarchaeota in chemical fertilizer treatments was significantly higher than that in organic manure treatments; N2 treatment was higher than N treatment; MN2 treatment was higher than MN treatment. The relative abundance of uncultured ammonia-oxidizing Crenarchaeote and uncultured Crenarchaeote in CK treatment was significantly higher than that in other treatments. Uncultured Crenarchaeote of relative abundance in chemical fertilizer treatments was higher than that in organic fertilizer treatments. The relative quantity of Thaumarchaeota in M treatment was higher than that in other treatments, whereas the relative abundance of Nitrososphaera.JG1 in organic manure treatments was significantly higher than that in chemical fertilizer treatments. It indicated that application of organic manure fertilizer increased the quantity.
Shannon Index of AOA and AOB
Shannon index of AOB (Figure 3A) in the N2 treatment was significantly lower (P < 0.05) than other treatments in both of two soil layers. Shannon index of AOB was higher in surface arable soil than that in plow pan soil. There was no significant difference (P > 0.05) between organic fertilizer treatments and CK treatment in Shannon index of AOB in plow pan soil, whereas the Shannon index of AOB was higher in organic manure treatments than that in chemical fertilizer treatments.
Figure 3. Shannon index of ammonia-oxidizing bacteria (A) and archaea (B) in different treatments in surface arable soil and plow pan soil, respectively. The mean value ± SD (n = 3). Different letters represent significant differences among fertilization treatments (LSD method, and p < 0.05). UP represented surface arable soil (0–20 cm), and Down represented plow pan soil (20–40 cm).
Shannon index of AOA (Figure 3B) in N and N2 treatments was significantly lower (P < 0.05) than that in other treatments in surface arable soil layer, whereas Shannon index of AOA in N, N2, and MN2 treatments in plow pan soil layer was significantly lower (P < 0.05) than that in other treatments. Shannon index of AOA in plow pan soil layer was higher than that in surface arable soil layer.
NMDS Analysis of AOA and AOB
In order to determine whether long-term fertilization and soil layers were related to changes in soil ammonia-oxidizing microbial community structure, we used NMDS based on Bray–Curtis dissimilarities to analyze the overall archaeal and bacterial community structural changes in Figure 4. The NMDS ordinations showed that the AOA (Figure 4A) community structure was divided into two parts: the organic fertilizer treatment and chemical fertilizer treatments. The community structure of AOB (Figure 4B) had a significant correlation with the soil layer, and the community structure of AOB was divided into two parts: the surface arable soil and the plow pan soil.
Figure 4. Non-metric multidimensional scaling (NMDS) of AOA (A) and AOB (B) in two soil layers under different fertilization treatments. UP represented surface arable soil (0–20 cm) and Down represented plow pan soil (20–40 cm). Points of different colors or shapes represent samples of different groups. The closer the two sample points are, the more similar the species composition of the two samples. The abscissa and ordinate indicate the relative distance and have no practical meaning. Stress: Test the pros and cons of NMDS analysis results. It is generally believed that when stress is less than 0.2, it can be represented by a two-dimensional dot diagram of NMDS, which has a certain explanatory meaning; when stress is less than 0.1, it can be considered as a good ranking.
The Relationship Among Soil Physicochemical Properties, Abundance, and Community Structure of AOA and AOB
The quantities of AOB in both soil layers (Figure 5) and AOA in surface arable soil layer (Figure 5A) had significant positive correlation (P < 0.05) with stocks of SOC, TN, TP, AP, AK, and NO3–-N under long-term fertilization, and the quantities of AOB in surface soil (Figure 5A) had significant positive correlation (P < 0.05) with PNR, whereas the quantity of AOA had significant positive correlation (P < 0.05) with PNR and NO3–-N in plow pan soil layer and had significant negative correlation with NH4+-N.
Figure 5. Redundancy analysis (RDA) of AOB and AOA in diversity and quantity in surface arable soil layer (A) and plow pan soil layer (B). RDA correlation plot showed changes in ammonia oxidizers diversity and quantity explained by soil environmental factors. SOC, soil organic carbon; TN, total nitrogen; TP, total phosphorus; TK, total potassium; AP, available phosphorus; AK, available potassium; NO3–-N, nitrate nitrogen; PNR, potential nitrite rate; SMC, soil moisture carbon.
Shannon index of AOA had no relationship with environmental factors in both two soil layers (Figure 5), whereas the Shannon index of AOB had significant positive correlation (P < 0.05) with stocks of SOC, TN, and AK in both two soil layers. Unclassified Nitrosomonadaceae, Nitrosomonadaceae, and Nitosospira.L115 had no significant correlation (P > 0.05) with environmental factors, whereas uncultured bacterium MF323203.1, uncultured Nitrosomonadaceae bacterium, and Nitrosospira.Nsp12 had significant positive correlation (P < 0.05) with PNR and stocks of SOC, TN, TP, and NO3–-N in surface arable soil layer (Figure 6A). Uncultured bacterium MF323203.1 had prominently positive correlation (P < 0.05) with TK, AP, and AK, whereas Nitrosospira.PJA1 had prominent correlation (P < 0.05) with SOC, TK, AP, and AK, and unclassified Nitrosospira only had significant correlation (P < 0.05) with AK. Nitrosospira.Ka3, uncultured bacterium MF323203.1, and environmental samples of Nitrosomonadaceae had no correlation (P < 0.05) with environmental factors. Uncultured Nitrosomonadaceae bacterium and unclassified Nitrosospira had significant positive correlation (P < 0.05) with PNR and stocks of TP, AP, and NO3–-N in plow pan soil layer (Figure 6B). Uncultured Nitrosomonadaceae bacterium had prominent correlation with soil pH. Uncultured bacterium MF323203.1, uncultured bacterium MF323620.1, and unclassified Nitrosospira were significantly correlated with SOC, TP, TK, AP, and AK, whereas Nitrosospira.PJA1 had significant correlation with pH, SOC, and TP.
Figure 6. Redundancy analysis (RDA) of AOB taxa on surface arable soil (A) and plow pan soil (B), respectively, and of AOA taxa on surface arable soil (C) and plow pan soil (D). RDA correlation plot showed variance in taxa of ammonia-oxidation microorganisms explained by soil environmental factors. SOC, soil organic carbon; TN, total nitrogen; TP, total phosphorus; TK, total potassium; AP, available phosphorus; AK, available potassium; NO3–-N, nitrate nitrogen; PNR, potential nitrite rate; DOC, dissolved organic carbon; SMC, soil moisture carbon.
Crenarchaeota, unclassified AOA, uncultured Crenarchaeota, and uncultured archaea had no significant correlation (P < 0.05) with the environmental factors, whereas Nitrososphaera.JG1 and Thaumarchaeota had significant positive correlation (P < 0.05) with PNR, soil pH, and stocks of SOC, TN, TP, TK, and NO3–-N in surface arable soil layer (Figure 6C). Unclassified AOA, uncultured ammonia-oxidizing Crenarchaeota, Nitrososphaera.JG1, and Thaumarchaeota had significant positive correlation (P < 0.05) with PNR and soil pH in plow pan soil layer, whereas Nitrososphaera.JG1 and Thaumarchaeota had significant positive correlation (P < 0.05) with soil pH and stocks of SOC, TN, TP, AP, AK, and NO3–-N in plow pan soil layer (Figure 6D).
Discussion
Long-term application of manure with chemical fertilizer could effectively increase soil nutrient content, improve soil physicochemical properties, and increase crop yields, thereby influencing microorganisms (Oehl et al., 2004; Shi et al., 2008; Feitosa de Souza et al., 2016; Lin et al., 2018), whereas long-term application of chemical nitrogen fertilizer alone was likely to decrease soil pH and microbial biomass and damage soil microbial structure (Sun et al., 2017; Huang et al., 2018; Li et al., 2019). Our experiment showed the same results, and the phosphorus fertilizer application could increase organic matters of soil (Table 1). Long-term application of manure with chemical fertilizer effectively slowed down soil acidification and increased SOC, TN, TP, AP, and NO3–-N in both soil layers, which might be due to the application of manure that could effectively improve physicochemical property of soil.
Our experiment showed that application of organic fertilizer could significantly increase the abundance of amoA gene in both soil layers (Figure 1). This result was consistent with results reported as the abundance of AOB in organic fertilizer treatment was higher than that in chemical fertilizer (Li et al., 2017, 2018), which was due to the enrichment of soil nutrients caused the increase of bacterial abundance. Our result has supported the above conjecture, and RDA in abundance of amoA gene with environment showed that the abundance of amoA gene was significant positively correlated with the stocks of SOC, TN, TP, AP, AK, and NO3–-N in both layers of soil (Figure 5).
The results of absolute abundance (Figures 1, 2) and relative abundance (Supplementary Figures S1, S2) analysis were different from each other on community structure of ammonia-oxidizing microorganisms. We adopted the combined determination results of absolute abundance with the proportion of relative abundance of each taxon in different treatments for subsequent analysis in order to get the analysis results closer to the real situation.
The application of organic fertilizer, balanced application of chemical fertilizer, and rational application of nitrogen fertilizer could improve the abundance of arch-amoA gene in surface arable soil (Figure 2A), whereas the application of organic fertilizer and nitrogen fertilizer could keep the abundance of arch-amoA gene constant in plow pan soil layer (Figure 2B). The abundance of arch-amoA gene was higher in plow pan soil layer than that in surface arable soil layer. This might be due to the oxygen content in surface arable soil layer being higher than that in plow pan soil layer, and parts of AOA belong to anaerobic heterotrophic bacteria under our experimental condition. Previous studies have confirmed that AOA could adapt to lower oxygen conditions (Liu et al., 2015), and AOB quantities were more easily decreased along with the sediment depth (Wang et al., 2014).
The results of NMDS (Figure 4) showed that the community structure of AOB was divided into two parts in all treatments according to soil layers. This result indicated that the taxa of AOB in different soil layers were significantly different. On the one hand, soil nutrient of plow pan soil layer was significantly lower than that of surface arable soil layer; on the other hand, dissolved oxygen could greatly affect the growth of bacteria (Canziani et al., 2006), whereas AOB were significantly increased in less bioturbated zone (Chen and Gu, 2017). Therefore, community structure of AOB in plow pan soil layer was significantly different from that in surface arable soil layer due to less soil nutrient, oxygen, and disturbance in plow pan soil layer.
The results of NMDS (Figure 4A) showed that the community structure of AOA was divided into two parts in both soil layers due to application of organic fertilizer. This result indicated that AOA had different requirements for nutrient forms. This was consistent with the result reported for readily mineralized organic matter, which was important for AOA (Jiang et al., 2012), and addition of glucose could influence the community of AOA in Beijing Miyun Reservoir (Wang et al., 2015).
PNR was the ability of soil to oxidize ammonium nitrogen to nitrate nitrogen. In our experiment, long-term application of organic fertilizer was more conducive to improving soil PNR; this might be due to that organic fertilizer enhanced the PNR, whereas chemical fertilizer suppressed PNR (Fan et al., 2011). Han et al. (2016) found that PNR was significant positively correlated with abundance of AOB but not AOA. Some studies (Qin et al., 2012; Xu et al., 2012) found that AOA, but not AOB, was significantly correlated with PNR in acidic soils of slope land, and Nguyen et al. (2018) found that PNR was significantly correlated with both AOA and AOB. However, we found that PNR was significantly correlated with the abundance of amoA of AOB in surface arable soil layer (Figure 5A) and significantly correlated with the abundance of AOA in plow pan soil layer (Figure 5B). The abundance of uncultured bacterium MF323203.1, uncultured Nitrosomonadaceae bacterium, and Nitrosospira.Nsp12 had significant correlation with PNR and were mostly influenced by the soil phosphorus (Figures 6A,B), whereas the stock of P was significantly higher in surface arable soil than that in plow pan soil (Table 1). Thus, we hypothesized that this might be related to the beneficial effects of phosphorus on the growth and reproduction of these AOB. The abundance of unclassified AOA, uncultured ammonia-oxidizing Crenarchaeota, Nitrososphaera.JG1, and Thaumarchaeota had significant correlation with PNR and was prominently higher in manure fertilizer treatments than that in chemical fertilizer treatments (Figures 2A,B). NMDS ordinations (Figure 4A) showed that AOA was divided into two parts by organic fertilizer. Thus, we hypothesized that abundance of these AOA could be increased by complex changes in the soil physical and chemical properties, which were due to organic fertilizer. These results might be due to that uncultured bacterium MF323203.1, uncultured Nitrosomonadaceae bacterium, Nitrosospira.Nsp12, Nitrososphaera.JG1, and Thaumarchaeota played a dominant role in PNR (Figure 6). However, the changes of abundance of uncultured bacterium MF323203.1, uncultured Nitrosomonadaceae bacterium, and Nitrosospira.Nsp12 that were positively correlated with PNR were consistent with the total abundance of AOB. The changes of abundance of Nitrososphaera.JG1 and Thaumarchaeota that had positive correlation with PNR were different from the total abundance of AOA in surface arable soil layer. The changes of abundance of unclassified AOA, uncultured ammonia-oxidizing Crenarchaeota, Nitrososphaera.JG1, and Thaumarchaeota were positively correlated with PNR and consistent with the total abundance of AOA. The changes of abundance of uncultured Nitrosomonadaceae bacterium and Nitrosospira.NSP12 were positively correlated with PNR and different from the total abundance of AOB in plow pan soil layer. PNR was correlated with unclassified AOA, uncultured ammonia-oxidizing Crenarchaeota, Nitrososphaera.JG1, Thaumarchaeota (AOA) and uncultured bacterium MF323203.1, uncultured Nitrosomonadaceae bacterium, and Nitrosospira.Nsp12 (AOB). PNR was also positively related to the abundance of AOB in surface arable layer and AOA in plow pan layer macroscopically. The essence was that these taxa affected the ammonia oxidation under long-term fertilization in brown soil. It might be a good explanation for some studies suggesting that AOA or AOB played a dominant role in various environments.
Conclusion
In conclusion, long-term fertilization changed soil PNR by influencing specific ammonia-oxidizing microorganisms, which are influenced by soil physiochemical properties in the brown soil. Only specific and rare AOA and AOB instead of the whole AOA or AOB community played a dominant role in nitrification in brown soil under long-term fertilization.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
FC, XH, and PL conceived, designed the study, and wrote the manuscript. JY, MI, SZ, NA, and JD collected and analyzed the data. All authors read and approved the final manuscript.
Funding
This work was supported by the National Key Research and Development Program of China (Grant No. 2018YFD0300308) and the National Natural Science Foundation of China (Grant No. 31972511).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.622454/full#supplementary-material
References
Adair, K. L., and Schwartz, E. (2008). Evidence that ammonia-oxidizing archaea are more abundant than ammonia-oxidizing bacteria in semiarid soils of northern Arizona. Microb. Ecol. 56, 420–426. doi: 10.1007/s00248-007-9360-9
Belay, A., Claassens, A., and Wehner, F. (2002). Effect of direct nitrogen and potassium and residual phosphorus fertilizers on soil chemical properties, microbial components and maize yield under long-term crop rotation. Biol. Fertil. Soils 35:420427. doi: 10.1007/s00374-002-0489-x
Bhattacharyya, R., Prakash, V., and Kundu, S. (2010). Long term effects of fertilization on carbon and nitrogen sequestration and aggregate associated carbon and nitrogen in the Indian sub-Himalayas. Nutr. Cycl. Agroecosyst. 86:116. doi: 10.1007/s10705-009-9270-y
Canziani, R., Emondi, V., Garavaglia, M., Malpei, F., Pasinetti, E., and Buttiglieri, G. (2006). Effect of oxygen concentration on biological nitrification and microbial kinetics in a cross-flow membrane bioreactor (MBR) and moving-bed biofilm reactor (MBBR) treating old landfill leachate. J. Membr. Sci. 286:202212. doi: 10.1016/j.memsci.2006.09.044
Chen, J., and Gu, J. (2017). Faunal burrows alter the diversity, abundance, and structure of AOA, AOB, anammox and n-Damo communities in coastal mangrove sediments. Microb. Ecol. 74:140156. doi: 10.1007/s00248-017-0939-5
Chen, J., Zhang, H., Liu, W., Lian, J., Ye, W., and Shen, W. (2015). Spatial distribution patterns of ammonia-oxidizing archaea abundance in subtropical forests at early and late successional stages. Sci. Rep. 5:16587.
Chen, S., Zhou, Y., Chen, Y., and Gu, J. (2018). fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884–i890. doi: 10.1093/bioinformatics/bty560
Chen, X. P., Zhu, Y. G., Xia, Y., Shen, J. P., and He, J. Z. (2008). Ammonia-oxidizing archaea: important players in paddy rhizosphere soil? Environ. Microbiol. 10, 1978–1987. doi: 10.1111/j.1462-2920.2008.01613.x
Dang, C., Liu, W., Lin, Y., Zheng, M., Jiang, H., Chen, Q., et al. (2018). Dominant role of ammonia-oxidizing bacteria in nitrification due to ammonia accumulation in sediments of Danjiangkou reservoir, China. Appl. Microbiol. Biotechnol. 102:33993410. doi: 10.1007/s00253-018-8865-0
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604
Erguder, T. H., Boon, N., Wittebolle, L., Marzorati, M., and Verstraete, W. (2009). Environmental factors shaping the ecological niches of ammonia-oxidizing archaea. FEMS Microbiol. Rev. 33, 855–869. doi: 10.1111/j.1574-6976.2009.00179.x
Fan, F., Yang, Q., Li, Z., Wei, D., Cui, X. A., and Liang, Y. (2011). Impacts of organic and inorganic fertilizers on nitrification in a cold climate soil are linked to the bacterial ammonia oxidizer community. Microb. Ecol. 62, 982–990.
Feitosa, de Souza, T. A., Rodrígues, A. F., and Marques, L. F. (2016). The trend of soil chemical properties, and rapeseed productivity under different long-term fertilizations and stubble management in a Ferralsols of Northeastern Brazil. Org. Agric. 7, 353–363. doi: 10.1007/s13165-016-0164-4
Gan, X. H., Zhang, F. Q., Gu, J. D., Guo, Y. D., Li, Z. Q., Zhang, W. Q., et al. (2016). Differential distribution patterns of ammonia-oxidizing archaea and bacteria in acidic soils of nanling national nature reserve forests in subtropical China. Antonie Van Leeuwenhoek 109, 237–251. doi: 10.1007/s10482-015-0627-8
Han, J., Shi, J., Zeng, L., Xu, J., and Wu, L. (2016). Impacts of continuous excessive fertilization on soil potential nitrification activity and nitrifying microbial community dynamics in greenhouse system. J. Soils Sediments 17, 471–480. doi: 10.1007/s11368-016-1525-z
Hua, W., Luo, P., An, N., Cai, F., Zhang, S., Chen, K., et al. (2020). Manure application increased crop yields by promoting nitrogen use efficiency in the soils of 40-year soybean-maize rotation. Sci. Rep. 10:14882. doi: 10.1038/s41598-020-71932-9
Huang, L., Dong, H., Wang, S., Huang, Q., and Jiang, H. (2013). Diversity and abundance of ammonia-oxidizing archaea and bacteria in diverse chinese paddy soils. Geomicrobiol. J. 31, 12–22. doi: 10.1080/01490451.2013.797523
Huang, X., Liu, Y., Li, Y., Guo, P., Fang, X., and Yi, Z. (2018). Foliage application of nitrogen has less influence on soil microbial biomass and community composition than soil application of nitrogen. J. Soils Sediments 19, 221–231. doi: 10.1007/s11368-018-2027-y
Jiang, X., Liu, W., Liu, Q., Jia, Z., Wright, A. L., and Cao, Z. (2012). Soil N mineralization, nitrification and dynamic changes in abundance of ammonia-oxidizing bacteria and archaea along a 2000 year chronosequence of rice cultivation. Plant Soil 365, 59–68. doi: 10.1007/s11104-012-1377-2
Kandeler, E. (1995). “Potential nitrification,” in Methods in Soil Biology, eds F. Schinner, R. Öhlinger, E. Kandeler, and A. Margesin (Berlin; Heidelberg: Springer), 146149.
Könneke, M., Bernhard, A. E., de la Torre, J. R., Walker, C. B., Waterbury, J. B., and Stahl, D. A. (2005). Isolation of an autotrophic ammonia-oxidizing marine archaeon. Nature 437, 543–546. doi: 10.1038/nature03911
Li, D., Chen, L., Xu, J., Ma, L., Olk, D. C., Zhao, B., et al. (2018). Chemical nature of soil organic carbon under different long-term fertilization regimes is coupled with changes in the bacterial community composition in a Calcaric Fluvisol. Biol. Fertil. Soils 54, 999–1012. doi: 10.1007/s00374-018-1319-0
Li, F., Chen, L., Zhang, J., Yin, J., and Huang, S. (2017). Bacterial community structure after long-term organic and inorganic fertilization reveals important associations between soil nutrients and specific taxa involved in nutrient transformations. Front. Microbiol. 8:187.
Li, M., and Gu, J. D. (2013). Community structure and transcript responses of anammox bacteria, AOA, and AOB in mangrove sediment microcosms amended with ammonium and nitrite. Appl. Microbiol. Biotechnol. 97, 9859–9874. doi: 10.1007/s00253-012-4683-y
Li, M., Cao, H., Hong, Y., and Gu, J. D. (2010). Spatial distribution and abundances of ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) in mangrove sediments. Appl. Microbiol. Biotechnol. 89, 1243–1254. doi: 10.1007/s00253-010-2929-0
Li, N., Kumar, P., Lai, L., Abagandura, G. O., Kumar, S., Nleya, T., et al. (2019). Response of soil greenhouse gas fluxes and soil properties to nitrogen fertilizer rates under camelina and carinata nonfood oilseed crops. BioEnergy Res. 12, 524–535. doi: 10.1007/s12155-019-09987-4
Lin, Y., Slessarev, E. W., Yehl, S. T., D’Antonio, C. M., and King, J. Y. (2018). Long-term nutrient fertilization increased soil carbon storage in california grasslands. Ecosystems 22, 754–766. doi: 10.1007/s10021-018-0300-y
Liu, S., Hu, B., He, Z., Zhang, B., Tian, G., Zheng, P., et al. (2015). Ammonia-oxidizing archaea have better adaptability in oxygenated/hypoxic alternant conditions compared to ammonia-oxidizing bacteria. Appl. Microbiol. Biotechnol. 99, 8587–8596. doi: 10.1007/s00253-015-6750-7
Luo, P., Han, X., Wang, Y., Han, M., Shi, H., Liu, N., et al. (2015). Influence of long-term fertilization on soil microbial biomass, dehydrogenase activity, and bacterial and fungal community structure in a brown soil of northeast China. Ann. Microbiol. 65, 533–542. doi: 10.1007/s13213-014-0889-9
Magoč, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507
Nguyen, L. T. T., Osanai, Y., Anderson, I. C., Bange, M. P., Braunack, M., Tissue, D. T., et al. (2018). Impacts of waterlogging on soil nitrification and ammonia-oxidizing communities in farming system. Plant Soil 426, 299–311. doi: 10.1007/s11104-018-3584-y
Oehl, F., Sieverding, E., Mader, P., Dubois, D., Ineichen, K., Boller, T., et al. (2004). Impact of long-term conventional and organic farming on the diversity of arbuscular mycorrhizal fungi. Oecologia 138, 574–583. doi: 10.1007/s00442-003-1458-2
Qin, H., Yuan, H., Zhang, H., Zhu, Y., Yin, C., Tan, Z., et al. (2012). Ammonia-oxidizing archaea are more important than ammonia-oxidizing bacteria in nitrification and NO3–N loss in acidic soil of sloped land. Biol. Fert. Soils 49, 767–776. doi: 10.1007/s00374-012-0767-1
Rotthauwe, J.-H., and Witzel, K.-P. (1997). The ammonia monooxygenase structural gene amoA as a functional marker: molecular fine-scale analysis of natural ammonia- oxidizing populations. Appl. Environ. Microbiol. 63, 4704–4712. doi: 10.1128/aem.63.12.4704-4712.1997
Shi, Z. H., Chen, L. D., Cai, C. F., Li, Z. X., and Liu, G. H. (2008). Effects of long-term fertilization and mulch on soil fertility in contour hedgerow systems: a case study on steeplands from the three gorges area, China. Nutr. Cycling Agroecosyst. 84, 39–48. doi: 10.1007/s10705-008-9223-x
Sun, Y. F., Shen, J. P., Zhang, C. J., Zhang, L. M., Bai, W. M., Fang, Y., et al. (2017). Responses of soil microbial community to nitrogen fertilizer and precipitation regimes in a semi-arid steppe. J. Soils Sediments 18, 762–774. doi: 10.1007/s11368-017-1846-6
Tao, R., Wakelin, S. A., Liang, Y., and Chu, G. (2017). Response of ammonia-oxidizing archaea and bacteria in calcareous soil to mineral and organic fertilizer application and their relative contribution to nitrification. Soil Biol. Biochem. 114, 20–30. doi: 10.1016/j.soilbio.2017.06.027
Wang, C., Zhu, G., Wang, W., and Yin, C. (2013). Preliminary study on the distribution of ammonia oxidizers and their contribution to potential ammonia oxidation in the plant-bed/ditch system of a constructed wetland. J. Soils Sediments 13, 1626–1635. doi: 10.1007/s11368-013-0750-y
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/aem.00062-07
Wang, X., Wang, C., Bao, L., and Xie, S. (2015). Impact of carbon source amendment on ammonia-oxidizing microorganisms in reservoir riparian soil. Ann. Microbiol. 65, 1411–1418. doi: 10.1007/s13213-014-0979-8
Wang, Z., Wang, Z., Huang, C., and Pei, Y. (2014). Vertical distribution of ammonia-oxidizing archaea (AOA) in the hyporheic zone of a eutrophic river in North China. World J. Microbiol. Biotechnol. 30, 1335–1346. doi: 10.1007/s11274-013-1559-y
Xu, Y., Yu, W., Ma, Q., and Zhou, H. (2012). Responses of bacterial and archaeal ammonia oxidisers of an acidic luvisols soil to different nitrogen fertilization rates after 9 years. Biol. Fert. Soils 48, 827–837. doi: 10.1007/s00374-012-0677-2
Keywords: ammonia-oxidizing archaea, ammonia-oxidizing bacteria, potential nitrification rate, brown soil, long-term fertilization
Citation: Cai F, Luo P, Yang J, Irfan M, Zhang S, An N, Dai J and Han X (2021) Effect of Long-Term Fertilization on Ammonia-Oxidizing Microorganisms and Nitrification in Brown Soil of Northeast China. Front. Microbiol. 11:622454. doi: 10.3389/fmicb.2020.622454
Received: 28 October 2020; Accepted: 22 December 2020;
Published: 04 February 2021.
Edited by:
Anne E. Taylor, Oregon State University, United StatesReviewed by:
Christoph Mueller, University of Giessen, GermanyCuijing Zhang, Shenzhen University, China
Copyright © 2021 Cai, Luo, Yang, Irfan, Zhang, An, Dai and Han. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Peiyu Luo, ibtyoufe@syau.edu.cn; Xiaori Han, hanxiaori@163.com