Chunhua Jia1,2†
Guixiang Zhou2†
Ling Ma1,2
Xiuwen Qiu3
Jiabao Zhang1,2*
Jingkuan Wang1
Congzhi Zhang2
Lin Chen2Donghao Ma2Zhanhui Zhao4Zaiqi Xue2- 1Northeast Key Laboratory of Conservation and Improvement of Cultivated Land (Shenyang), Ministry of Agriculture and Rural Affairs, College of Land and Environment, Shenyang Agricultural University, Shenyang, China
- 2Fengqiu Experimental Station of National Ecosystem Research Network of China, State Key Laboratory of Soil and Sustainable Agriculture, Institute of Soil Science, Chinese Academy of Sciences, Nanjing, China
- 3College of Landscape Architecture, Jiangsu Vocational College of Agriculture and Forestry, Jurong, China
- 4School of Geomatics and Urban Spatial Informatics, Henan University of Urban Construction, Pingdingshan, Henan, China
Introduction: Ammonia oxidizing archaea (AOA) and ammonia oxidizing bacteria (AOB) have been proven to be key microorganisms driving the ammonia oxidation process. However, under different fertilization practices, there is a lack of research on the impact of interaction between predators and AOA or AOB on nitrogen cycling at the multi-trophic level.
Methods: In this study, a network-oriented microscopic culture experiment was established based on four different long-term fertilization practices soils. We used the nitrification inhibitors 2-phenyl-4,4,5,5-tetramethylimidazoline-1-oxide-3-oxyl (PTIO) and 3, 4-Dimethylpyrazole phosphate (DMPP) inhibited AOA and AOB, respectively, to explore the impact of interaction between protists and AOA or AOB on nitrogen transformation.
Results: The results showed that long-term nitrogen application promoted the potential nitrification rate (PNR) and nitrous oxide (N2O) emission, and significantly increased the gene abundance of AOB, but had no obvious effect on AOA gene abundance. DMPP significantly reduced N2O emission and PNR, while PTIO had no obvious effect on them. Accordingly, in the multi-trophic microbial network, Cercozoa and Proteobacteria were identified as keystone taxa of protists and AOB, respectively, and were significantly positively correlated with N2O, PNR and nitrate nitrogen. However, Nitrososphaerota archaeon as the keystone species of AOA, had an obvious negative linkage to these indicators. The structural equation model (SEM) showed that AOA and AOB may be competitors to each other. Protists may promote AOB diversity through direct trophic interaction with AOA.
Conclusion: The interaction pattern between protists and ammonia-oxidizing microorganisms significantly affects potential nitrification rate and N2O emission, which has important implications for soil nitrogen cycle.
1 Introduction
Soil nitrogen cycle is one of the cores of elemental cycling in soil ecosystems, which mainly includes four processes, biological nitrogen fixation, ammonification, nitrification and denitrification (Wang et al., 2021). Nitrification is an intermediate link between nitrogen fixation and denitrification (Lehnert et al., 2018). Nitrification determines the effective utilization of nitrogen by plants, but also is directly related to a series of ecological and environmental problems such as soil acidification, and greenhouse gas nitrous oxide (N2O) release (Aryal et al., 2022). The ammonia oxidation process is the first and rate-limiting step of nitrification, and is the central link of the global nitrogen cycle. It is mainly driven by ammonia oxidizing archaea (AOA) and ammonia oxidizing bacteria (AOB) (Zhao et al., 2020; Kaur-Bhambra et al., 2021). Previous studies of ammonia-oxidizing microorganisms affecting nitrogen transformation mainly focused on effects of abiotic factors at the same trophic level. For example, AOB was significantly correlated with potential nitrification rate (PNR) and dominated nitrification rather than AOA in calcareous soils (Zou et al., 2022). Urea nitrogen in paddy soil promoted N2O production by AOB more than ammonium nitrogen, and the N2O production by AOA was more sensitive to ammonium nitrogen (Fu et al., 2020). However, the effects of microbial interactions on nitrogen cycle under different fertilization practices and at the multi-trophic level of predation relationships (e.g., protists) have rarely been reported. Therefore, it is necessary to add discriminative inhibitors to silence the activity of AOA and AOB, in order to better explain the trophic interactions between predators and ammonia-oxidizing microorganisms.
2-phenyl-4,4,5,5-tetramethylimidazoline-1-oxide-3-oxyl (PTIO) is a stable and water-soluble organic free radical that can combine with NO to form NO2, and it is usually used as a nitric oxide scavenger (Lin et al., 2023). It has been reported that low concentrations of PTIO have a strong inhibitory effect on AOA rather than AOB (Sun et al., 2022). For example, under the condition of 100 μM PTIO, PTIO can inhibit the ammonia oxidation process by inhibiting AOA (Duan et al., 2018; Choi et al., 2023), and thus it was used to selectively inhibit the activity of AOA in the presence of AOB (Kits et al., 2019). 3,4-dimethylpyrazole phosphate (DMPP) is a highly soluble compound with high stability and persistent suppression effects (Li et al., 2023). Studies have shown that DMPP effectively inhibited the metabolic activity of AOB rather than AOA, which in turn led to a reduction of nitrification rate (Lei et al., 2022). Furthermore, in calcareous soils where the nitrification process was regulated by AOB, the addition of DMPP significantly suppressed N2O emissions (Bai et al., 2020). In general, these findings demonstrated that PTIO and DMPP can serve as effective inhibitors of AOA and AOB.
To investigate the effects of biological factors on nitrification process at the multi-trophic level, we focused on trophic interactions between protists and ammonia-oxidizing microorganisms. Protists, as predators of soil micro-food webs, can control the abundance and function of microbial species by preying on bacteria, fungi, or nematodes (Geisen et al., 2015; Geisen, 2016). Thus, the nutrient turnover rate of soil food web is increased, which ultimately regulates the decomposition of organic matter and affects nutrient cycling. Predatory behavior of protists leads to interactions between organisms within the soil microbiome to influence the response of the protist community to nitrogen application practices (Gao et al., 2019). Protists resorbed growth-limiting nutrients by feeding on bacterial cells in oligotrophic environments and ultimately stimulated the proliferation of primary producers and other bacteria (Chen et al., 2021). In addition, the interaction between protists and other microorganisms can promote the efficient utilization of organic nitrogen by AM fungal hyphae. Due to the fact that AM fungal hyphae acquire nitrogen more efficiently when both protists and prokaryotes are present in the organic nitrogen zone (Rozmoš et al., 2022). Thus, indirect (top-down control) and direct (heterotrophic) effects of protists on soil nutrient cycling perhaps increase the availability of nitrogen in the soil (Koller et al., 2013) and there is increasing attention to the role of protists in nutrient cycling. Although both protists and ammonia-oxidizing microorganisms have important implications for the soil nitrogen cycle, the effects of trophic interactions between protists and AOA or AOB on nitrification process remain unknown.
In general, understanding the effects of abiotic and biotic factors on nitrification is critically important for the regulation of soil nitrogen cycling processes. This study conducted microcosm experiments using four different long-term fertilization soils, using high-throughput sequencing technology and incorporating structural equation modeling (SEM). The aim was to (1) assess the efficacy of fertilization practices and ammonia-oxidizing microorganisms in influencing N2O emissions and PNR. (2) Reveal the trophic interactions between protists and AOA or AOB and their effects on nitrification process at the multi-trophic level.
2 Materials and methods
2.1 Field site and sampling
In June 2022, the soil was sampled in a long-term field experiment at Fengqiu Agricultural Ecological Experimental Station of Chinese Academy of Sciences (35°01′N, 114°32′E). The experiment has been conducted since 2010. The basic soil properties at the beginning of the experiment were shown in Table 1. This study collected soil samples from four fertilization treatments: NN (No straw and no nitrogen fertilization), CF (No straw and conventional fertilization), SM (Straw mulching and no nitrogen fertilization), SMCF (Straw mulching and conventional fertilization). NN and SM treatments were applied without nitrogen and only phosphorus and potassium fertilizers were applied. Conventional fertilization of CF and SMCF treatments indicated that both nitrogen, phosphorus and potassium fertilizers were applied. In all fertilizer treatments, the total amount of nitrogen (pure nitrogen) applied to wheat throughout the whole growth period was 210 kg ha–1, and the phosphorus and potassium fertilizers were 157 kg ha–1 for P2O5 and 105 kg ha–1 for K2O, respectively. Straw mulching was done by cutting corn straw into 20–50 mm strips after the previous season’s corn harvest, and then evenly covering the soil surface with straw after tilling (at a depth of about 15 cm). All the crop straw from each plot was returned to the same plot. The soil samples (0–20 cm depth) were collected using the five-point diagonal sampling method. Subsequently, the samples were put in sterile bags and temporarily stored in foam boxes with dry ice and taken back to the laboratory at the station within 2 h. Adequate dry ice was used to ensure that the soil samples were in an environment below 0°C. All the soil samples were passed through a 2 mm diameter mesh sieve to remove plant debris and stones, and the remaining fine root and straw residues were manually removed. The soil samples were stored at 4°C before constructing the micro world.
Table 1. Basic soil properties of the field.
2.2 Soil microcosm experiment
Soil samples were pre-incubated at 28°C for 7 days with moisture adjusted to 50% of water holding capacity. A 21-day microcosm experiment was established by placing 20 g (equivalent dry weight) of pre-incubated soil in 120 ml serum bottles. The following treatments were set up for each fertilized soil: (i) water only (CK), (ii) urea only (U), (iii) urea + PTIO (UP), and (iv) urea + DMPP (UD). Urea was added at 200 mg N kg–1 dry soil and DMPP was added at 1.5% of urea N (Fan et al., 2019; Yin et al., 2021) and the concentration of PTIO was set to 100 μM, (Jung et al., 2014). Three replications were set up for all treatments and a total of 144 microcosms were established. Gas samples (20 ml) were collected from the top with a syringe at 1, 2, 3, 4, 7, 10, 14, and 21 days after incubation. N2O concentrations were analyzed by gas chromatograph (Agilent 7890, Agilent Technologies, Santa Clara, CA, USA). To maintain aerobic conditions, all remaining bottles were vented for 30 min after each sampling and then resealed. Soil ammonium nitrogen (NH4+-N) and nitrate nitrogen (NO3–-N) contents were determined by destructive collection of three bottles from each treatment on days 0, 7, 14, and 21 of incubation. A portion of the last sampled soil was stored at −80°C for DNA extraction.
2.3 DNA extraction, real-time fluorescence quantitative PCR analysis and high-throughput sequencing
Total DNA was extracted from 0.5 g of fresh Soil using the Fast DNA® Spin Kit for Soil (MP Biomedicals, Inc.). The extracted DNA concentration and purity were examined by NanoDrop2000 (Thermo Fisher Scientific, Waltham, MA, USA). We used primers TAReuk454FWD1 (CCAGCASCYGCGGTA ATTCC)/TAReukREV3 (ACTTTCGTTCTTGATYRA) (Stoeck et al., 2010), Arch-amoA26F (GACTACATATTCTACACWGACT GGGC)/Arch-amoA417R (GGTGTCATATATTGGAGGCAACGT TGG) and amoA1F (GGGGTTTCTACTGGTGGT)/amoA2R(CCC CTCKGSAAAGCCTTCTTC) (Xia et al., 2011) amplified protists, AOA and AOB, respectively. AOA and AOB were analyzed on the qTOWER3/G amplification apparatus (Analytik Jena AG, Germany) with three replicates per sample. The standard curve is constructed from 8 series of diluents containing known copy number plasmids of the target gene. The reaction system for qPCR was 20 μl, including 1 μl of DNA template, 10 μl of TB Green® Premix Ex Taq TM, 0.25 μl forward and 0.25 μl reverse primers, and 8.5 μl of sterilized double-distilled water. Negative control was performed using sterilized distilled water instead of the template DNA sample as the reaction template.
Amplicon sequence variants (ASVs) were generated from raw sequence data using DADA2 (Callahan et al., 2016). Briefly, after removing the primer and adapt sequences using “cutadapt” (Martin, 2011), the sequences were filtered and trimmed using “filterAndTrim” with the following parameters: maxN = 0, maxEE = c(2,2), truncQ = 2. Chimeric sequences were detected and removed using “removeBimeraDenovo” based on the “consensus” method. The taxonomy assignment for protist was performed with the “assignTaxonomy” function (minimal bootstrap value 50) based on the PR2 database (V4.14.0) (Guillou et al., 2012). Taxonomy of the amoA ASVs was assigned using BLASTn based on the nucleotide database in NCBI with the default algorithm parameters (expect threshold was 0.05, word size was 28, match/mismatch scores were 1, -2, gap costs was linear, and e-value was 0.00001) (Kõljalg et al., 2005). The unassigned ASVs were filter for further study. The ASVs of 18S rRNA assigned as Rhodophyta, Streptophyta, Metazoa, Fungi, and unidentified Opisthokonta were removed for protist community analysis. All high-throughput sequencing steps were performed by Guangdong Magi Gene Technology Co., Ltd. (Guangzhou, China).
2.4 Statistical analysis
One-way analysis of variance (ANOVA) in SPSS Statistics 21 was used to evaluate the effects of different treatments on N2O emission and copy number of ammonia-oxidizing microorganisms. The “vegan” package in R v4.2.1 software was used to conduct principal coordinate analysis (PCoA) (Anderson and Willis, 2003). Pearson correlation analysis was conducted using R v4.2.1 software packages reshape2 and psych to examine the correlation between key taxa and soil physicochemical properties, as well as between key taxa. A heat map was generated using the “pheatmap” package.
The co-occurrence network analysis based on Spearman’s correlation matrix was calculated in the “Hmisc” package in R v4.2.1 (Wang et al., 2022). ASVs with a sum of relative abundance higher than 90% were selected, and then two correlations between all ASVs were calculated by combining protists (744 PASVs), AOA (105 ASVs) and AOB (67 BASVs) into one relative abundance table. To distinguish different microorganisms, we used PASV to represent the ASV of protists and BASV to represent the ASV of AOB. A correlation threshold higher than 0.6 and a P-value lower than 0.01 are considered to be stable co-occurrence relationships. Adjustments were made by using test corrections from the Benjamini-Hochberg procedure to reduce the chances of obtaining false positive results (Benjamini and Hochberg, 1995). The co-occurrence network was visualized and analyzed by “Gephi” software (Bastian et al., 2009). Subsequently, the connections within modules (Zi) and between modules (Pi) were calculated using the package “igraph” in R v4.2.1. According to the topological characteristics of nodes, node attributes were classified into four types, including: Module hubs, nodes with high connectivity within a module (Zi > 2.5 and Pi < 0.62). Connectors, nodes with high connectivity between two modules (Zi < 2.5 and Pi > 0.62). Network hubs, nodes with high connectivity throughout the network (Zi > 2.5 and Pi > 0.62), and Peripherals, nodes that do not have high connectivity both within and between modules (Zi < 2.5 and Pi < 0.62) (Deng et al., 2016). Moreover, according to Banerjee et al. (2018) and Shi et al. (2020), module hubs connectors, and network hubs were defined as keystone species in this study. Among them, the taxa of module hubs generally mediate interactions between species within the module and have significant implications for module function by regulating energy and material exchange within the module (Guimera and Nunes Amaral, 2005; Toju et al., 2018). Therefore, this study mainly focused on module hubs.
The SEM quantified the direct and indirect effects of soil physicochemical properties, protists, AOA, and AOB on PNR and N2O. Soil physicochemical properties (soil NO3–-N and NH4+-N) and microbial diversity (protists, AOA and AOB diversity) were represented by PCoA data of corresponding indicators. The maximum likelihood estimation method was applied to fit the model, and the parameters were consistent with relatively low chi-square and P > 0.05, goodness fit index (GFI) > 0.90 and root mean square error of approximation (RMSEA) < 0.05 (Grace and Keeley, 2006). All SEM analyses were performed using AMOS 21.
3 Results
3.1 Changes in N2O emission rates and PNR
Nitrous oxide emission rate was significantly affected by soil fertilization practices and nitrification inhibitors. For the urea (alone) treatments, the N2O emission rate of the long-term no-nitrogen treatments (NN, SM) peaked on day 4 (Figures 1A, C), and the long-term nitrogen treatments (CF, SMCF) reached the maximum rates on day 2 (Figures 1B, D). Meanwhile, long-term nitrogen and straw treatment (SMCF) had the highest N2O emission rate (Figure 1D). N2O emissions ranged from 18.84∼64.32 pmol g–1 h–1 under different fertilization practices, thus, long-term nitrogen and straw application promoted N2O emission by 3.4 times. The N2O emission rate was reduced at different degrees by the application of nitrification inhibitors. In all soils, N2O emission rates of U and UP treatments rapidly increased during the first week of incubation. During the whole incubation period, the UD treatment significantly decreased the N2O emission rate, and the emission trend was essentially consistent with CK. The range of N2O emissions under nitrification inhibitor treatments was 1.03∼64.32 pmol g–1h–1, which reduced N2O emissions by 62.4 times. In contrast, nitrification inhibitors were much more effective on N2O emissions than fertilizer treatments. Additionally, compared with U treatment, NO3–-N increased slowly in the UD treatment, and the UP treatment reduced NH4+-N content (Supplementary Figure 1).
Figure 1. Dynamic of N2O emission rates in the four fertilization practices from NN (A), CF (B), SM (C) and SMCF (D) under four treatments, including CK, U, UP and UD during the experimental incubation. Vertical bars indicate standard errors of the means (n = 3). NN, no straw and no nitrogen fertilization; CF, no straw and conventional fertilization; SM, straw mulching and no nitrogen fertilization; SMCF, straw mulching and conventional fertilization. CK, water only; U, urea only; UP, urea + PTIO; UD, urea + DMPP.
The effect of soil fertilization practices on PNR was significant only in the CK treatment, i.e., the long-term nitrogen application treatments (CF and SMCF) obviously increased PNR under the condition of no exogenous urea addition. Moreover, long-term nitrogen and straw application promoted PNR by 7.5 times. However, the addition of urea during incubation reduced the differences of soil fertilization practices. The U and UP treatments had the highest PNR in all soils. Compared with U treatment, the PNR was significantly increased in UP treatment on day 7 of NN and CF treatments, while there was a decreasing trend on the day 14 and 21. The PNR of UP treatment decreased with the extension of incubation time compared to the U treatment in SM treatment. The UD treatment significantly decreased PNR, with a trend consistent with CK or even obviously lower than CK (Table 2). PNR was reduced by 80 times after treatment with nitrification inhibitors. In general, the effect of fertilization practices on PNR was much lower than that of nitrification inhibitors.
Table 2. PNR during incubation of four nitrification inhibitors treatments (CK, U, UP, UD) within four fertilization practices (NN, CF, SM, SMCF).
3.2 Changes in gene abundance and analysis of microbial community structure
To explore the effects of fertilization practices and nitrification inhibitors on the activity of ammonia-oxidizing microorganisms, we analyzed the changes of AOA and AOB gene abundance. For AOA, the highest gene abundance of all soil treatments was found in the NN treatment on day 7. Moreover, the UP treatment reduced AOA gene abundance in all soils. Unexpectedly, obviously higher AOA gene abundance was observed in UD and CK treatments than in U treatment (Figures 2A–D). For AOB, CF treatment had the highest gene abundance on day 7. It is noteworthy that the gene abundance of SMCF treatment was lower than that of CF treatment. The effect of nitrification inhibitors on the AOB gene abundance under different fertilization practices showed the same trend. AOB gene abundance was significantly improved in U treatment and obviously decreased in UD treatment for all soils (Figures 2E–H). It appeared that AOA was better suited to oligotrophic environments, while AOB was more active in nitrogen-rich environments.
Figure 2. Dynamic of the AOA (A–D) and AOB (E–H) gene copies in the four fertilization practices from NN (A,E), CF (B,F), SM (C,G), and SMCF (D,H) under four treatments, including CK, U, UP, and UD during the experimental incubation. Vertical bars indicate standard errors of the means, different lowercase letters represent significant differences within groups (n = 3). NN, no straw and no nitrogen fertilization; CF, no straw and conventional fertilization; SM, straw mulching and no nitrogen fertilization; SMCF, straw mulching and conventional fertilization; CK, water only; U, urea only; UP, urea + PTIO; UD, urea + DMPP.
We further analyzed the response of microbial community structure to fertilization practices and nitrification inhibitors. Soil fertilization practices and nitrification inhibitor treatments had different effects on the community structure of protist, AOA, and AOB. Notably, all microbial communities in NN and SMCF treatments showed obvious separation, suggesting that long-term nitrogen and straw application significantly changed the microbial community structure (Figure 3). For protists, CK treatment was significantly detached from other treatments, while there was no obvious variation in community structure between U, UP, and UD treatments. For AOA, significant changes in community structure occurred in the UP treatment. Compared with nitrification inhibitor treatments, different fertilization soils had a greater impact on AOB community structure. For example, the AOB community structure in soils with long-term non-nitrogen application was clearly apart from soils with long-term nitrogen application.
Figure 3. The PCoA analysis of protists, AOA and AOB within four fertilization practices (NN, CF, SM, SMCF) under four treatments, including CK, U, UP, and UD. NN, no straw and no nitrogen fertilization; CF, no straw and conventional fertilization; SM, straw mulching and no nitrogen fertilization; SMCF, straw mulching and conventional fertilization; CK, water only; U, urea only; UP, urea + PTIO; UD, urea + DMPP.
3.3 Microbial interaction networks and key taxa analysis under nitrification inhibitor treatments
In order to determine the impact of nitrification inhibitors on the symbiotic relationship of soil microorganisms, the co-occurring networks of protists, AOA, and AOB were constructed for CK, U, UP, and UD treatments, respectively, and the topological parameters of each network were calculated (Figure 4A). Compared to CK, the network nodes, edges, average degree, clustering coefficient and network density were increased for all urea-added treatments. Notably, all topological properties of UD treatment were reduced, resulting in a lower complexity of the network (Supplementary Table 1). Meanwhile, the direct microbial interaction showed the lowest positive correlation ratio (79.92%) under UD treatment. Nodes with higher within module connectivity (Module hubs, Zi > 2.5 and Pi < 0.62) were identified as key taxa in the network. These taxa may play an important role in maintaining the stability of network structure. There were no module hubs in CK treatment (Figure 4B), while UP and UD treatments had more key taxa than U treatment. The key taxa identified in the U, UP, UD treatments were mainly Cercozoa, Lobosa, Sagenista, Crenarchaeota and Proteobacteria (Supplementary Table 2). Moreover, the relative abundances of key taxa were also highly correlated with physicochemical properties associated with nitrogen transformation (Figure 5).
Figure 4. The co-occurring network of protist, AOA and AOB under four treatments, including CK, U, UP and UD (A). Classification of nodes to identify keystone species within the networks (B). NN, no straw and no nitrogen fertilization; CF, no straw and conventional fertilization; SM, straw mulching and no nitrogen fertilization; SMCF, straw mulching and conventional fertilization; CK, water only; U, urea only; UP, urea + PTIO; UD, urea + DMPP.
Figure 5. Pearson correlations between key taxa and physicochemical properties (A) and within key taxa (B). *P < 0.05, **P < 0.01, ***P < 0.001. NN, no straw and no nitrogen fertilization; CF, no straw and conventional fertilization; SM, straw mulching and no nitrogen fertilization; SMCF, straw mulching and conventional fertilization; CK, water only; U, urea only; UP, urea + PTIO; UD, urea + DMPP.
N2O, PNR, NO3–-N, and NH4+-N are indicators closely related to the nitrification process. However, in contrast to the other indicators, NH4+-N was mainly negatively correlated with key species of microorganisms. Among the key taxa of the protists, the proportion of key taxa with significant positive correlations with N2O, PNR and NO3–-N was 50%, while the proportion of negative correlations was only 25%. The key taxa of protists with positive correlations were PASV1204, PASV634, PASV23, and PASV393, classified at the phylum level as Cercozoa, Lobosa and Sagenista. For AOA, only ASV56 was strongly linked to the four physicochemical properties, but the correlation was mainly negative. For AOB, although there were only two key taxa, BASV1 was positively associated with both N2O and PNR, and belonged to the Proteobacteria at the phylum level (Figure 5A). Additionally, the proportion of key species of AOA that were significantly negatively related to PASV1204, PASV634, PASV23, and PASV393 was 11%, and the proportion of positively correlated key species was 50% for AOB. For AOA, only ASV56 was significantly negatively linked to key species of protists, and BASV1 of AOB had the opposite pattern to ASV56 (Figure 5B).
3.4 Effects of microbial interactions on PNR and N2O under nitrification inhibitor treatments
The SEM was constructed to further clarify the relationship between soil properties, microorganisms, PNR and N2O under different nitrification inhibitor treatments (Figure 6). In all treatments, AOA had a direct negative effect on protists, conversely, AOB had a direct positive impact on protists. In CK and U treatment, the influences of microorganisms on PNR and N2O were mostly insignificant (Figures 6A, B). Microorganisms were more closely related to each other and to PNR and N2O with the addition of nitrification inhibitors (Figures 6C, D). In terms of standardized total effects, the addition of nitrification inhibitors resulted in a greater contribution of AOA and AOB on protist compared to the U treatment. To be specific, the interaction between AOB and protists increased the contribution to N2O after inhibiting AOA, the suppression of AOB enhanced the contribution of interaction between AOA and protists to PNR (Supplementary Figure 2). Overall, when AOA was suppressed, only N2O was significantly affected by interaction between AOB and protist. The interaction between AOA and protists obviously affected both N2O and PNR with high contributions when AOB was inhibited.
Figure 6. Structural equation model (SEM) analysis of the relationships between soil properties, protist, AOA, AOB and N2O and PNR under CK (A), U (B), UP (C), UD (D) treatments. The arrow width is proportional to the strength of the relationship. Green edge: positive correlation; Red edge: negative correlation. Adjacent numbers in the same direction as the arrow indicate the correlation coefficient. Paths with non-significant coefficients are shown as gray lines. The significance levels are indicated by *P < 0.05, **P < 0.01, ***P < 0.001. NN, no straw and no nitrogen fertilization; CF, no straw and conventional fertilization; SM, straw mulching and no nitrogen fertilization; SMCF, straw mulching and conventional fertilization; CK, water only; U, urea only; UP, urea + PTIO; UD, urea + DMPP.
4 Discussion
4.1 Effects of soil fertilization practices and nitrification inhibitors on N2O emission and PNR
The highest N2O emissions occurred in treatment with long-term nitrogen and straw application (Figure 1), which is consistent with recent studies (Xia et al., 2018). Long-term nitrogen application stimulates N2O emissions by affecting biological processes such as nitrification and denitrification in soil (Wu et al., 2018; Bai et al., 2020). Meanwhile, straw application can alleviate the limitation of soil N2O emission by providing carbon and nitrogen sources. In particular, straw addition enhances the availability of unstable carbon, which provided energy for the growth of N2O-producing microorganisms (e.g., AOA and AOB) in the soil (Wu et al., 2020).
As is well known, N2O is an important byproduct of nitrification or denitrification (Ma et al., 2012). The ammonia oxidation process driven by AOA and AOB is the initial and rate limiting step of nitrification (Hu and He, 2017; Lin et al., 2020). Our study confirmed that N2O emission was slightly reduced when AOA was inhibited and significantly reduced when AOB was suppressed (Figure 1). It indicated that AOB may be the main driver of N2O emission during the nitrification process in this study. Firstly, nitrification denitrification and incomplete NH2OH oxidation are the causes of N2O production by AOB in alkaline soils (Wu et al., 2018). Secondly, AOA seems to lack nitrification-denitrification ability due to the absence of NO reductase (Walker et al., 2010). Therefore, the effect of suppressing AOA on N2O emissions was not obvious. Except for N2O, PNR is also an important chemical index to characterize the nitrification of aerobic ammoxidation microorganisms, which can be used to quantify nitrification potential. Our study revealed that long-term nitrogen application significantly increased PNR, while PNR was reduced sharply after AOB was suppressed (Table 2). On the one hand, long-term nitrogen application provided sufficient substrate for nitrifying microorganisms and increased microbial activity (Wang et al., 2019), which in turn improved PNR. On the other hand, suppression of AOB reduced PNR probably owing to the fact that AOB is the major contributor to nitrification potential in neutral and alkaline soils (Ren et al., 2023). Thus, inhibition of AOB by DMPP indirectly suppressed PNR.
4.2 Effects of soil fertilization practices and nitrification inhibitors on AOA and AOB gene abundance
In this study, the long-term no nitrogen treatment had the highest AOA gene abundance (Figure 2A). This can be explained by the fact that AOA is dominant in nitrogen-limited environment (Wang et al., 2015), and long-term non-nitrogen soils provide favorable conditions for AOA to compete for ecological niches. In contrast to AOA, the peak of AOB gene abundance appeared in long-term nitrogen treatment (CF). It is due to that AOB is suited to survive in alkaline and nitrogen-enriched environments, and the application of nitrogen can provide high soil NH4+-N concentration for AOB growth (Verhamme et al., 2011; Sun et al., 2023). Additionally, we found that AOB gene abundance was lower in SMCF treatment than in CF treatment (Figure 2H). The incorporation of straw increased soil C/N, thereby delaying the nitrification process (Elrys et al., 2021; Wang et al., 2023), and thus the decrease of AOB gene abundance is understandable.
The addition of nitrification inhibitors had different effects on AOA and AOB. Our results showed that PTIO significantly inhibited the AOA gene abundance (Figures 2A–D). This may be explained by the high sensitivity of AOA to PTIO (Martens-Habbena et al., 2015). Interestingly, AOA gene abundance was obviously increased after suppressing AOB in long-term nitrogen fertilization soils (Figures 2A–D). This was inconsistent with the results of most studies that revealed AOA populations were unchanged or decreased (Di and Cameron, 2011). Firstly, there may be substrate competition between AOA and AOB, and the suppression effect of DMPP may free AOA from the competition of AOB (Fan et al., 2019). Secondly, DMPP, as an organic compound, may serve as a useful C-substrate to promote AOA growth (Florio et al., 2014). Meanwhile, the application of DMPP strongly inhibited the AOB gene abundance (Figures 2E–H) in all soils. It was worth noting that AOA activity was not affected by DMPP (Figures 2A–D). In fact, the different responses of AOA and AOB to DMPP inhibition may be attributed to different nitrifying enzyme systems. The AOA amoA sequence is more similar to the gene encoding bacterial micromethane monooxygenase (pMMO), resulting in functional differences between the AMO of AOA and AOB (Jung et al., 2014). The nitrification inhibitor DMPP decreases AMO activity, whereas MMO is unaffected (Weiske et al., 2001). Therefore, DMPP specifically inhibited AOB activity.
4.3 Interactions between protists and AOA or AOB
In contrast to AOA and AOB, protists are significantly affected by nitrogen fertilizer rather than inhibitors. Protists are the predators of soil micro-food web and the response to nitrogen are sensitive (Zhao et al., 2021). To explore the effect of interactions between protists and AOA or AOB on nitrogen transformation, a symbiotic network between them was constructed. We found that the addition of urea increased the complexity of the soil microbiome network and caused closer network connectivity (Supplementary Table 1). As nitrogen addition increased the availability of resources and food in the soil, it enhanced the interactions within the microbiome to improve resource utilization efficiency (Shi et al., 2016; Yu et al., 2018). Remarkably, the network complexity was lower when AOB was suppressed by DMPP. This may be due to the key taxa of AOB was inhibited, which can maintain network stability and increase network complexity. It was also confirmed by the fact that AOB was not identified in the keystone species treated with UD (Figure 4).
Nodes with high connectivity within the module (Zi > 2.5 and Pi < 0.62) were identified as key taxa. These taxa may function through energy and substance exchange within the module (Deng et al., 2016). When they are removed or disappeared from a particular ecosystem, it can lead to drastic changes in that ecosystem (Ikegami, 2005). In this study, the key species of protists after inhibitor addition mainly included PASV1204, PASV634, PASV23, and PASV393 (Supplementary Table 2), all of which were significantly positively correlated with N2O, PNR and NO3–-N. We further classified the key species and determined that they belong to Cercozoa, Lobosa and Sagenista at the phylum level. However, only Cercozoa was the protist species that appeared consistently since suppressing AOA or AOB, which was classified as generalist in microbial networks (Barberán et al., 2012). It emphasized the important role in microbial symbiosis. Cercozoa belongs to the predatory protist in the functional classification. It generally affects soil nutrient mineralization by top-down regulation of microbial biomass and activity and by affecting the acquisition of soil nitrogen and carbon enzyme activities (Wei et al., 2021; Wu et al., 2022). Additionally, BASV1, belonging to Proteobacteria at the phylum level, is the only AOB key species identified after inhibiting AOA. It is strongly positively associated with N2O, PNR, and NO3–-N, and key species of protist (Figure 5). It is well known that Proteobacteria is a key denitrification species with the potential to produce and reduce N2O (Xiang et al., 2023).
For AOA, on the one hand, N2O and PNR showed no significant changes or increased trends after AOA was inhibited. Meanwhile, we found that only ASV56 was significantly negatively correlated with both N2O and PNR among the key species of AOA. On the other hand, only ASV56 was negatively linked to keystone species of protists. Furthermore, the SEMs of four treatments showed a significant negative correlation between AOA and protists. The standardized total effects also indicated high negative contributions of AOA to protists. Therefore, we speculated that among the key species of AOA, ASV56 may play a major role in microbial interactions in this study. Since the key taxa of AOA was not identified, the relevant functions were analyzed by NCBI similar sequence alignment (Supplementary Table 3). By comparison, the key AOA species may be Nitrososphaerota archaeon, which can obtain energy during ammonia oxidation and utilize carbon dioxide as a chemoautotrophic carbon source (Wei et al., 2023). The interactions between protists and AOA or AOB may provide substrates and motivation for the soil nitrification process and stimulate nitrogen transformation. Specifically, Cercozoa had a higher C/N than sub-trophic level organisms, and they excreted excess nitrogen making it available to other microorganisms (Jousset, 2017). Meanwhile, Proteobacteria and Nitrososphaerota archaeon can also use the nitrogen secreted by protists to promote PNR and the N2O emission. Recent studies have found that the symbiotic relationship between Cercozoa and Proteobacteria was closely related to soil nitrogen transformation process (Xiao et al., 2022). The predator-prey interactions between Cercozoa and Thaumarchaeota increased nutrient turnover (Gu et al., 2023). These all support the above view.
The SEM analysis also provided a basis for the above argument. Protists were found to be the highly negatively correlated with AOA and significant positively associated with AOB. The addition of nitrification inhibitors resulted in a strong negative linkage between AOA and AOB. At the same time, standardized total effects indicated that AOA or AOB was closely related to protists. When one of AOA and AOB was inhibited, the other interacted with protists to increase contributions on N2O or PNR. Since the fact that both AOA and AOB use ammonia nitrogen as energy substrates, there is a direct nutritional competition between them (French et al., 2021; Sun et al., 2022). When one of the microorganisms was suppressed, the interaction between the other and protists may be more closely. Thus, the addition of nitrification inhibitors may promote the interactions between ammonia-oxidizing microorganisms and protists, and the trophic interaction may stimulate nitrogen turnover (Gao et al., 2019). Therefore, the addition of nitrification inhibitors is a noteworthy agricultural strategy to improve the soil nitrogen cycle.
5 Conclusion
Both abiotic and biotic factors are important factors affecting nitrogen transformation indicators N2O and PNR. The highest N2O emissions in different soil fertilization practices were found in long-term nitrogen and straw application treatment, and the PNR was also significantly increased in the treatment of long-term nitrogen application. In addition, our results demonstrated that the addition of PTIO and DMPP significantly inhibited the activity of AOA and AOB, respectively, and AOB was the main microorganism driving N2O emissions and stimulating PNR in alkaline fluvo-aquic soil. Finally, the relationship between keystone taxa and N2O and PNR in the network directly reflected the relevance of corresponding microorganisms with these indicators. Protists were significantly negatively associated with AOA and positively linked to AOB. There may be competition between AOA and AOB, and protists may promote the diversity of AOB through direct trophic interaction with AOA, thereby affecting the nitrogen transformation process. In the future, we can regulate the interaction between protists and ammonia oxidizing microorganisms by modulating key taxa of protists, AOA and AOB, and finally achieve the goal of managing soil nitrogen cycle in the field.
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 below: NCBI–PRJNA1035458, PRJNA1035811, and PRJNA1035828.
Author contributions
CJ: Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing—original draft, Writing—review and editing. GZ: Conceptualization, Validation, Supervision, Writing—review and editing. LM: Writing—review and editing. XQ: Writing—review and editing. JZ: Writing—review and editing, Conceptualization, Funding acquisition, Project administration, Resources, Validation. JW: Methodology, Writing—review and editing. CZ: Methodology, Writing—review and editing. LC: Resources, Writing—review and editing. DM: Resources, Writing—review and editing. ZZ: Data curation, Writing—review and editing. ZX: Data curation, Writing—review and editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Key Research Development Program Projects (2022YFD1500203 and 2022YFD1500401), the Strategic Research and Consulting Project of Chinese Academy of Engineering (2023-JB-06-05), the China Agriculture Research System (CARS-03-15 and CARS-52), the National Natural Science Foundation of China (41967002 and 42277336), the Natural Science Foundation of Jiangsu Province (BK20221561), and the Key Science and Technology Project of Jiangsu Province (BE2019387), the Strategic Priority Research Program of Chinese Academy of Sciences (XDA28020203).
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.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1337507/full#supplementary-material
References
Anderson, M. J., and Willis, T. J. (2003). Canonical analysis of principal coordinates: A useful method of constrained ordination for ecology. Ecology 84, 511–525.
Aryal, B., Gurung, R., Camargo, A. F., Fongaro, G., Treichel, H., Mainali, B., et al. (2022). Nitrous oxide emission in altered nitrogen cycle and implications for climate change. Environ. Pollut. 314:120272. doi: 10.1016/j.envpol.2022.120272
Bai, J., Qiu, S., Jin, L., Wei, D., Xu, X., Zhao, S., et al. (2020). Quantifying soil N pools and N2O emissions after application of chemical fertilizer and straw to a typical chernozem soil. Biol. Fertil. Soils 56, 319–329. doi: 10.1007/s00374-019-01422-2
Banerjee, S., Schlaeppi, K., and van der Heijden, M. G. A. (2018). Keystone taxa as drivers of microbiome structure and functioning. Nat. Rev. Microbiol. 16, 567–576. doi: 10.1038/s41579-018-0024-1
Barberán, A., Bates, S. T., Casamayor, E. O., and Fierer, N. (2012). Using network analysis to explore co-occurrence patterns in soil microbial communities. ISME J. 6, 343–351. doi: 10.1038/ismej.2011.119
Bastian, M., Heymann, S., and Jacomy, M. (2009). “Gephi: An open source software for exploring and manipulating networks,” in Proceedings of the International AAAI Conference on Web and Social Media, Vol. 3, (Palo Alto, CA: AAAI Press), 361–362. doi: 10.1609/icwsm.v3i1.13937
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B-Stat. Methodol. 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Chen, B., Xiong, W., Qi, J., Pan, H., Chen, S., Peng, Z., et al. (2021). Trophic interrelationships drive the biogeography of protistan community in agricultural ecosystems. Soil Biol. Biochem. 163:108445. doi: 10.1016/j.soilbio.2021.108445
Choi, E., Chaudhry, S. I., and Martens-Habbena, W. (2023). Role of nitric oxide in hydroxylamine oxidation by ammonia-oxidizing bacteria. Appl. Environ. Microbiol. 89:e0217322. doi: 10.1128/aem.02173-22
Deng, Y., Zhang, P., Qin, Y. J., Tu, Q. C., Yang, Y. F., He, Z. L., et al. (2016). Network succession reveals the importance of competition in response to emulsified vegetable oil amendment for uranium bioremediation. Environ. Microbiol. 18, 205–218. doi: 10.1111/1462-2920.12981
Di, H. J., and Cameron, K. C. (2011). Inhibition of ammonium oxidation by a liquid formulation of 3,4-Dimethylpyrazole phosphate (DMPP) compared with a dicyandiamide (DCD) solution in six New Zealand grazed grassland soils. J. Soils Sediments 11, 1032–1039. doi: 10.1007/s11368-011-0372-1
Duan, P., Wu, Z., Zhang, Q., Fan, C., and Xiong, Z. (2018). Thermodynamic responses of ammonia-oxidizing archaea and bacteria explain N2O production from greenhouse vegetable soils. Soil Biol. Biochem. 120, 37–47. doi: 10.1016/j.soilbio.2018.01.027
Elrys, A. S., Wang, J., Metwally, M. A., Cheng, Y., Zhang, J. B., Cai, Z. C., et al. (2021). Global gross nitrification rates are dominantly driven by soil carbon-to-nitrogen stoichiometry and total nitrogen. Glob. Change Biol. 27, 6512–6524. doi: 10.1111/gcb.15883
Fan, X., Yin, C., Chen, H., Ye, M., Zhao, Y., Li, T., et al. (2019). The efficacy of 3, 4-dimethylpyrazole phosphate on N2O emissions is linked to niche differentiation of ammonia oxidizing archaea and bacteria across four arable soils. Soil Biol. Biochem. 130, 82–93. doi: 10.1016/j.soilbio.2018.11.027
Florio, A., Clark, I. M., Hirsch, P. R., Jhurreea, D., and Benedetti, A. (2014). Effects of the nitrification inhibitor 3,4-dimethylpyrazole phosphate (DMPP) on abundance and activity of ammonia oxidizers in soil. Biol. Fertil. Soils 50, 795–807. doi: 10.1007/s00374-014-0897-8
French, E., Kozlowski, J. A., and Bollmann, A. (2021). Competition between ammonia-oxidizing archaea and bacteria from freshwater environments. Appl. Environ. Microbiol. 87:e01038-21. doi: 10.1128/AEM.01038-21
Fu, Q., Xi, R., Zhu, J., Hu, H., Xing, Z., and Zuo, J. (2020). The relative contribution of ammonia oxidizing bacteria and archaea to N2O emission from two paddy soils with different fertilizer N sources: A microcosm study. Geoderma 375:114486. doi: 10.1016/j.geoderma.2020.114486
Gao, Z., Karlsson, I., Geisen, S., Kowalchuk, G., and Jousset, A. (2019). Protists: Puppet masters of the rhizosphere microbiome. Trends Plant Sci. 24, 165–176. doi: 10.1016/j.tplants.2018.10.011
Geisen, S. (2016). The bacterial-fungal energy channel concept challenged by enormous functional versatility of soil protists. Soil Biol. Biochem. 102, 22–25. doi: 10.1016/j.soilbio.2016.06.013
Geisen, S., Rosengarten, J., Koller, R., Mulder, C., Urich, T., and Bonkowski, M. (2015). Pack hunting by a common soil amoeba on nematodes. Environ. Microbiol. 17, 4538–4546. doi: 10.1111/1462-2920.12949
Grace, J. B., and Keeley, J. E. (2006). A structural equation model analysis of postfire plant diversity in California shrublands. Ecol. Appl. 16, 503–514.
Gu, S., Zhou, X., Yu, H., Yan, H., Wang, Y., Liu, Y., et al. (2023). Microbial and chemical fertilizers for restoring degraded alpine grassland. Biol. Fertil. Soils 59, 911–926. doi: 10.1007/s00374-023-01759-9
Guillou, L., Bachar, D., Audic, S., Bass, D., Berney, C., Bittner, L., et al. (2012). The Protist Ribosomal Reference database (PR2): A catalog of unicellular eukaryote small sub-unit rRNA sequences with curated taxonomy. Nucleic Acids Res. 41, D597–D604. doi: 10.1093/nar/gks1160
Guimera, R., and Nunes Amaral, L. A. (2005). Functional cartography of complex metabolic networks. Nature 433, 895–900. doi: 10.1038/nature03288
Hu, H. W., and He, J. Z. (2017). Comammox-a newly discovered nitrification process in the terrestrial nitrogen cycle. J. Soils Sediments 17, 2709–2717. doi: 10.1007/s11368-017-1851-9
Ikegami, T. (2005). Neutral phenotypes as network keystone species. Popul. Ecol. 47, 21–29. doi: 10.1007/s10144-005-0208-1
Jousset, A. (2017). “Application of protists to improve plant growth in sustainable agriculture,” in Rhizotrophs: Plant Growth Promotion to bioremediation, ed. S. Mehnaz (Singapore: Springer), 263–273. doi: 10.1007/978-981-10-4862-3_13
Jung, M. Y., Well, R., Min, D., Giesemann, A., Park, S. J., Kim, J. G., et al. (2014). Isotopic signatures of N2O produced by ammonia-oxidizing archaea from soils. ISME J. 8, 1115–1125. doi: 10.1038/ismej.2013.205
Kaur-Bhambra, J., Wardak, D. L., Prosser, J. I., and Gubry-Rangin, C. (2021). Revisiting plant biological nitrification inhibition efficiency using multiple archaeal and bacterial ammonia-oxidising cultures. Biol. Fertil. Soils 58, 241–249. doi: 10.1007/s00374-020-01533-1
Kits, K. D., Jung, M. Y., Vierheilig, J., Pjevac, P., Sedlacek, C. J., Liu, S., et al. (2019). Low yield and abiotic origin of N2O formed by the complete nitrifier Nitrospira inopinata. Nat. Commun. 10:1836. doi: 10.1038/s41467-019-09790-x
Kõljalg, U., Larsson, K. H., Abarenkov, K., Nilsson, R. H., Alexander, I. J., Eberhardt, U., et al. (2005). UNITE: A database providing web-based methods for the molecular identification of ectomycorrhizal fungi. New Phytol. 166, 1063–1068. doi: 10.1111/j.1469-8137.2005.01376.x
Koller, R., Rodriguez, A., Robin, C., Scheu, S., and Bonkowski, M. (2013). Protozoa enhance foraging efficiency of arbuscular mycorrhizal fungi for mineral nitrogen from organic matter in soil to the benefit of host plants. New Phytol. 199, 203–211. doi: 10.1111/nph.12249
Lehnert, N., Dong, H. T., Harland, J. B., Hunt, A. P., and White, C. J. (2018). Reversing nitrogen fixation. Nat. Rev. Chem. 2, 278–289. doi: 10.1038/s41570-018-0041-7
Lei, J., Fan, Q., Yu, J., Ma, Y., Yin, J., and Liu, R. (2022). A meta-analysis to examine whether nitrification inhibitors work through selectively inhibiting ammonia-oxidizing bacteria. Front. Microbiol. 13:962146. doi: 10.3389/fmicb.2022.962146
Li, Z., Xu, P., Han, Z., Wu, J., Bo, X., Wang, J., et al. (2023). Effect of biochar and DMPP application alone or in combination on nitrous oxide emissions differed by soil types. Biol. Fertil. Soil. 59, 123–138. doi: 10.1007/s00374-022-01688-z
Lin, Y., Duan, C., Fan, J., Hu, H. W., He, Z. Y., Ye, G., et al. (2023). Nitrification inhibitor 1-octyne inhibits growth of comammox Nitrospira but does not alter their community structure in an acidic soil. J. Soils Sediments 23, 989–997. doi: 10.1007/s11368-022-03367-w
Lin, Y., Ye, G., Ding, W., Hu, H. W., Zheng, Y., Fan, J., et al. (2020). Niche differentiation of comammox Nitrospira and canonical ammonia oxidizers in soil aggregate fractions following 27-year fertilizations. Agric. Ecosyst. Environ. 304:107147. doi: 10.1016/j.agee.2020.107147
Ma, Y., Wang, J., Zhou, W., Yan, X., and Xiong, Z. (2012). Greenhouse gas emissions during the seedling stage of rice agriculture as affected by cultivar type and crop density. Biol. Fertil. Soils 48, 589–595. doi: 10.1007/s00374-011-0656-z
Martens-Habbena, W., Qin, W., Horak, R. E. A., Urakawa, H., Schauer, A. J., Moffett, J. W., et al. (2015). The production of nitric oxide by marine ammonia-oxidizing archaea and inhibition of archaeal ammonia oxidation by a nitric oxide scavenger. Environ. Microbiol. 17, 2261–2274. doi: 10.1111/1462-2920.12677
Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 17, 10–12. doi: 10.14806/ej.17.1.200
Ren, B., Wang, W., Shen, L., Yang, W., Yang, Y., Jin, J., et al. (2023). Nitrogen fertilization rate affects communities of ammonia-oxidizing archaea and bacteria in paddy soils across different climatic zones of China. Sci. Total Environ. 902:166089. doi: 10.1016/j.scitotenv.2023.166089
Rozmoš, M., Bukovská, P., Hršelová, H., Kotianová, M., Dudáš, M., Gančarčíková, K., et al. (2022). Organic nitrogen utilisation by an arbuscular mycorrhizal fungus is mediated by specific soil bacteria and a protist. ISME J. 16, 676–685. doi: 10.1038/s41396-021-01112-8
Shi, S., Nuccio, E. E., Shi, Z. J., He, Z., Zhou, J., and Firestone, M. K. (2016). The interconnected rhizosphere: High network complexity dominates rhizosphere assemblages. Ecol. Lett. 19, 926–936. doi: 10.1111/ele.12630
Shi, Y., Delgado-Baquerizo, M., Li, Y., Yang, Y., Zhu, Y.-G., Peñuelas, J., et al. (2020). Abundance of kinless hubs within soil microbial networks are associated with high functional potential in agricultural ecosystems. Environ. Int. 142:105869. doi: 10.1016/j.envint.2020.105869
Stoeck, T., Bass, D., Nebel, M., Christen, R., Jones, M. D., Breiner, H. W., et al. (2010). Multiple marker parallel tag environmental DNA sequencing reveals a highly complex eukaryotic community in marine anoxic water. Mol. Ecol. 19, 21–31. doi: 10.1111/j.1365-294X.2009.04480.x
Sun, D., Tang, X., Li, J., Liu, M., Hou, L., Yin, G., et al. (2022). Chlorate as a comammox Nitrospira specific inhibitor reveals nitrification and N2O production activity in coastal wetland. Soil Biol. Biochem. 173:108782. doi: 10.1016/j.soilbio.2022.108782
Sun, R., Wang, D., Guo, Z., Hua, K., Guo, X., Chen, Y., et al. (2023). Combined application of organic manure and chemical fertilizers stabilizes soil N-cycling microflora. Soil Ecol. Lett. 5:220165. doi: 10.1007/s42832-022-0165-z
Toju, H., Tanabe, A. S., and Sato, H. (2018). Network hubs in root-associated fungal metacommunities. Microbiome 6, 1–16. doi: 10.1186/s40168-018-0497-1
Verhamme, D. T., Prosser, J. I., and Nicol, G. W. (2011). Ammonia concentration determines differential growth of ammonia-oxidising archaea and bacteria in soil microcosms. ISME J. 5, 1067–1071. doi: 10.1038/ismej.2010.191
Walker, C. B., de la Torre, J. R., Klotz, M. G., Urakawa, H., Pinel, N., Arp, D. J., et al. (2010). Nitrosopumilus maritimus genome reveals unique mechanisms for nitrification and autotrophy in globally distributed marine crenarchaea. Proc. Natl. Acad. Sci. U. S. A. 107, 8818–8823. doi: 10.1073/pnas.0913533107
Wang, B., Zhao, J., Guo, Z., Ma, J., Xu, H., and Jia, Z. (2015). Differential contributions of ammonia oxidizers and nitrite oxidizers to nitrification in four paddy soils. ISME J. 9, 1062–1075. doi: 10.1038/ismej.2014.194
Wang, G., George, T. S., Pan, Q., Feng, G., and Zhang, L. (2022). Two isolates of Rhizophagus irregularis select different strategies for improving plants phosphorus uptake at moderate soil P availability. Geoderma 421:115910. doi: 10.1016/j.geoderma.2022.115910
Wang, J., Wang, J., Rhodes, G., He, J. Z., and Ge, Y. (2019). Adaptive responses of comammox Nitrospira and canonical ammonia oxidizers to long-term fertilizations: Implications for the relative contributions of different ammonia oxidizers to soil nitrogen cycling. Sci. Total Environ. 668, 224–233. doi: 10.1016/j.scitotenv.2019.02.427
Wang, Q., Han, Y., Lan, S., and Hu, C. (2021). Metagenomic insight into patterns and mechanism of nitrogen cycle during biocrust succession. Front. Microbiol. 12:633428. doi: 10.3389/fmicb.2021.633428
Wang, S., Zhai, L., Guo, S., Zhang, F., Hua, L., and Liu, H. (2023). Returned straw reduces nitrogen runoff loss by influencing nitrification process through modulating soil C: N of different paddy systems. Agric. Ecosyst. Environ. 354:108438. doi: 10.1016/j.agee.2023.108438
Wei, T. S., Gao, Z. M., Gong, L., Li, Q. M., Zhou, Y. L., Chen, H. G., et al. (2023). Genome-centric view of the microbiome in a new deep-sea glass sponge species Bathydorus sp. Front. Microbiol. 2023:1078171. doi: 10.3389/fmicb.2023.1078171
Wei, X., Ge, T., Wu, C., Wang, S., Mason-Jones, K., Li, Y., et al. (2021). T4-like phages reveal the potential role of viruses in soil organic matter mineralization. Environ. Sci. Technol. 55, 6440–6448. doi: 10.1021/acs.est.0c06014
Weiske, A., Benckiser, G., Herbert, T., and Ottow, J. (2001). Influence of the nitrification inhibitor 3,4-dimethylpyrazole phosphate (DMPP) in comparison to dicyandiamide (DCD) on nitrous oxide emissions, carbon dioxide fluxes and methane oxidation during 3 years of repeated application in field experiments. Biol. Fertil. Soils 34, 109–117. doi: 10.1007/s003740100386
Wu, B., Zhou, L., Liu, S., Liu, F., Saleem, M., Han, X., et al. (2022). Biogeography of soil protistan consumer and parasite is contrasting and linked to microbial nutrient mineralization in forest soils at a wide-scale. Soil Biol. Biochem. 165:108513. doi: 10.1016/j.soilbio.2021.108513
Wu, D., Zhao, Z. C., Han, X., Meng, F. Q., Wu, W. L., Zhou, M. H., et al. (2018). Potential dual effect of nitrification inhibitor 3, 4-dimethylpyrazole phosphate on nitrifier denitrification in the mitigation of peak N2O emission events in North China Plain cropping system. Soil Biol. Biochem. 121, 147–153. doi: 10.1016/j.soilbio.2018.03.010
Wu, L., Hu, R., Tang, S., Shaaban, M., Zhang, W., Shen, H., et al. (2020). Nitrous oxide emissions in response to straw incorporation is regulated by historical fertilization. Environ. Pollut. 266:115292. doi: 10.1016/j.envpol.2020.115292
Xia, L., Lam, S. K., Wolf, B., Kiese, R., Chen, D., and Butterbach-Bahl, K. (2018). Trade-offs between soil carbon sequestration and reactive nitrogen losses under straw return in global agroecosystems. Glob. Change Biol. 24, 5919–5932. doi: 10.1111/gcb.14466
Xia, W., Zhang, C., Zeng, X., Feng, Y., Weng, J., Lin, X., et al. (2011). Autotrophic growth of nitrifying community in an agricultural soil. ISME J. 5, 1226–1236. doi: 10.1038/ismej.2011.5
Xiang, H., Hong, Y., Wu, J., Wang, Y., Ye, F., Ye, J., et al. (2023). Denitrification contributes to N2O emission in paddy soils. Front. Microbiol. 14:1218207. doi: 10.3389/fmicb.2023.1218207
Xiao, D., He, X., Zhang, W., Hu, P., Sun, M., and Wang, K. (2022). Microbiological mechanism underlying vegetation restoration across climatic gradients in a karst ecosystem. Land Degrad. Dev. 33, 3245–3259. doi: 10.1002/ldr.4385
Yin, C., Fan, X., Chen, H., Jiang, Y., Ye, M., Yan, G., et al. (2021). 3, 4-Dimethylpyrazole phosphate is an effective and specific inhibitor of soil ammonia-oxidizing bacteria. Biol. Fertil. Soil 57, 753–766. doi: 10.1007/s00374-021-01565-1
Yu, J., Deem, L. M., Crow, S. E., Deenik, J. L., and Penton, C. R. (2018). Biochar application influences microbial assemblage complexity and composition due to soil and bioenergy crop type interactions. Soil Biol. Biochem. 117, 97–107. doi: 10.1016/j.soilbio.2017.11.017
Zhao, J., Bello, M. O., Meng, Y., Prosser, J. I., and Gubry-Rangin, C. (2020). Selective inhibition of ammonia oxidising archaea by simvastatin stimulates growth of ammonia oxidising bacteria. Soil Biol. Biochem. 141:107673. doi: 10.1016/j.soilbio.2019.107673
Zhao, Z., He, J., Geisen, S., Han, L., Wang, J., Shen, J., et al. (2021). Protist communities are more sensitive to nitrogen fertilization than other microorganisms in diverse agricultural soils. Microbiome 7:33. doi: 10.1186/s40168-019-0647-0
Keywords: protists, ammonia-oxidizing bacteria (archaea), keystone taxa, predatory relationship, N2O emission, potential nitrification rate
Citation: Jia C, Zhou G, Ma L, Qiu X, Zhang J, Wang J, Zhang C, Chen L, Ma D, Zhao Z and Xue Z (2024) The addition of discrimination inhibitors stimulations discrimination potential and N2O emissions were linked to predation among microorganisms in long term nitrogen application and straw returning systems. Front. Microbiol. 14:1337507. doi: 10.3389/fmicb.2023.1337507
Received: 13 November 2023; Accepted: 19 December 2023;
Published: 09 January 2024.
Edited by:
Xingjia Xiang, Anhui University, ChinaReviewed by:
Chenhao Zhao, Yangzhou University, ChinaWu Xiong, Nanjing Agricultural University, China
Wenchen Song, Minzu University of China, China
Copyright © 2024 Jia, Zhou, Ma, Qiu, Zhang, Wang, Zhang, Chen, Ma, Zhao and Xue. 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: Jiabao Zhang, jbzhang@issas.ac.cn
†These authors have contributed equally to this work