- 1Liaoning Key Laboratory of Development and Utilization for Natural Products Active Molecules, School of Chemistry and Life Science, Anshan Normal University, Anshan, Liaoning, China
- 2Key Laboratory of Vegetation Ecology, Jilin Songnen Grassland Ecosystem National Observation and Research Station, Institute of Grassland Science, Northeast Normal University, Ministry of Education, Changchun, China
- 3Remote Sensing Laboratories, Department of Geography, University of Zürich, Winterthurerstrasse, Zürich, Switzerland
- 4Fredericton Research and Development Centre, Agriculture and Agri-Food Canada, Fredericton, NB, Canada
Introduction: Changes in plant diversity and increased atmospheric nitrogen deposition independently influence soil nitrogen cycling in terrestrial ecosystems. However, the interactive effects of plant diversity and nitrogen deposition on soil nitrogen cycling multifunctionality (NCMF) in grassland ecosystems remain poorly understood.
Methods: We conducted a fully factorial microcosm experiment to quantify the responses and underlying mechanism of soil NCMF to nitrogen addition (0, 5, and 10 g N m−2 yr.−1) and plant diversity gradients (1, 3, and 6 species).
Results: Our results revealed a significant interactive effect between plant diversity and nitrogen addition on soil NCMF. Specifically, high plant diversity alleviated the negative effects of nitrogen addition on soil NCMF. The addition of nitrogen reduced the soil pH, which imposed microbial stress by limiting carbon availability. In contrast, higher plant diversity increased soil organic matter via below-ground carbon inputs, thereby reducing the soil carbon limitation of microorganims and enhancing the soil NCMF.
Discussion: Overall, our findings suggest that maintaining or enhancing plant diversity in grasslands could be a key strategy to mitigate the adverse effects of atmospheric nitrogen deposition on soil nitrogen cycling, highlighting the crucial role of plant diversity in regulating ecosystem nutrient cycling under global change.
1 Introduction
Nitrogen is a critical nutrient for both plants and soil microorganisms (Wang et al., 2023). Nitrogen cycling, as an essential component of biogeochemical cycles, plays a pivotal role in supporting organic growth and maintaining ecosystem stability (Chen and Chen, 2021), and has been demonstrated to be strongly linked to changes in plant diversity and nitrogen deposition (Hallin et al., 2009; Maestre et al., 2011). However, previous studies have primarily examined the individual effects of plant diversity and nitrogen deposition on nitrogen cycling (Lu et al., 2010; Byrnes et al., 2014), and their interactive effects remain unclear, hampering our ability to understand the responses and adaptations of ecosystem nutrient cycling to plant diversity and nitrogen deposition.
Many studies have shown that plant diversity and nitrogen deposition significantly affect soil nitrogen cycling functions, but the results are inconsistent. For example, plant diversity is positively correlated with soil ammonium nitrogen availability and nitrogen mineralization rates, whereas soil nitrate nitrogen availability tends to decrease with increasing plant diversity (Oelmann et al., 2011). Moreover, plant species mixtures can increase the abundance of soil nitrifiers and denitrifiers (Shah et al., 2013). Additionally, previous studies have shown that nitrogen addition enhances nitrification by increasing the abundance of ammonia-oxidizing bacteria (AOB) and the Nitrospira genus (Ma et al., 2016; Yang et al., 2020). However, other studies have reported the neutral impacts of nitrogen fertilization on nitrogen mineralization rates and the inhibitory effects on denitrification (e.g., narG gene abundance) (Tang et al., 2016). Furthermore, leucine aminopeptidase (LAP) and β-1,4-N-acetylglucosaminidase (NAG), the hydrolytic enzymes related to nitrogen cycling, have been found to decrease under high nitrogen addition (Cui et al., 2020). These results indicate that the response of individual nitrogen cycling functions to plant diversity and nitrogen deposition significantly differs and that there is a certain trade-off, because soil nitrogen cycling functionality is inherently linked to multiple ecosystem processes rather than isolated functions (Cheng et al., 2022; Wang et al., 2023; Yang et al., 2023). Therefore, assessing single functions cannot fully capture how soil nitrogen cycling responds to environmental changes in real-world scenarios (Cui et al., 2020). This highlights the importance of exploring the relationship between soil nitrogen cycling multifunctionality (NCMF) and environmental conditions in grassland ecosystems.
Increasing attention has been given to the response of soil NCMF to changes in plant diversity and nitrogen deposition. Both plant diversity and nitrogen addition have been shown to independently alter soil NCMF (Cui et al., 2020; Zhang et al., 2021; Freitag et al., 2023). For example, greater plant diversity has been associated with greater soil NCMF, likely due to increased resource heterogeneity, niche differentiation, and complementary resource use among plant species, increasing ecosystem stability and productivity (Byrnes et al., 2014). Similarly, diverse plant communities can also support a wider range of microbial taxa involved in nitrogen cycling, such as nitrifiers, denitrifiers, and nitrogen-fixing bacteria (Chen et al., 2020), thereby enhancing soil NCMF. Nitrogen addition has been shown to increase soil NCMF (Cui et al., 2020) because it provides a readily available nutrient source for soil microbes, stimulating their activity and growth (Li et al., 2020b). However, excessive nitrogen addition can lead to soil acidification, changes in microbial community composition, and potential declines in microbial diversity (Ma et al., 2016), which may negatively affect soil NCMF over time. Nevertheless, the interactive effects of plant diversity and nitrogen addition on soil NCMF remain unclear.
The interactive effects of plant diversity and nitrogen deposition on soil NCMF have not been fully explored. Although separate studies on plant diversity and nitrogen deposition suggest that both factors may influence NCMF through distinct mechanisms (Fornara and Tilman, 2008; Cui et al., 2020; Delgado-Baquerizo et al., 2020; Li et al., 2020b; Freitag et al., 2023). The interactive effects of plant diversity and nitrogen addition on soil NCMF are complex and may depend on several factors, including the magnitude of nitrogen input, the composition of the plant community, and the initial soil nutrient status. The combination of high plant diversity and nitrogen addition may lead to synergistic effects on soil NCMF. For example, diverse plant communities may be better able to utilize added nitrogen, leading to increased soil microbial activity and increased NCMF (Shah et al., 2013; Meyer et al., 2016). Increased nitrogen deposition alters soil environmental conditions, such as the available nitrogen content and soil pH, which directly affect microbial diversity and activity (Wang et al., 2011; Wu et al., 2024), thereby influencing nitrogen cycling processes. However, the presence of diverse plant species may also buffer against the negative effects of nitrogen addition, such as soil acidification, by maintaining a more stable soil environment (Turnbull et al., 2012; van der Plas, 2019; Xu et al., 2022). Higher plant diversity has been associated with increased soil microbial biomass and activity due to greater plant-derived carbon inputs and expanded microbial niches (García-Palacios et al., 2010; Chen et al., 2020), which can enhance both mineralization rates and immobilization rates of soil nitrogen (Fornara and Tilman, 2008; Meyer et al., 2016). Additionally, increased soil carbon resulting from increased plant diversity can meet the energy demands of carbon-intensive microbial processes and nitrogen cycling functional gene expression, thereby accelerating soil nitrogen turnover (Shu et al., 2011; Cui et al., 2021). Nitrogen addition may reduce the positive effects of plant diversity on soil NCMF. For instance, nitrogen addition leads to the dominance of a few plant species that outcompete others, which could reduce plant diversity and, in turn, negatively impact soil microbial diversity and function (Chen et al., 2020). Furthermore, complementarity in resource use and or facilitative interactions were the main drivers of increased productivity at higher levels of species richness (van Ruijven et al., 2005). Higher plant diversity communities have varied root architectures and growth forms, which could lead to more efficient exploration of soil volume (Xu et al., 2022). In contrast, nitrogen addition would homogenize soil conditions and reduce the spatial heterogeneity of soil resources (He et al., 2021), potentially reducing the positive effects of plant diversity on soil NCMF. Nevertheless, it remains unclear whether and to what extent plant diversity modulates the effects of nitrogen addition on soil NCMF.
In this study, we conducted a two-year microcosm experiment to investigate the effects and potentially microbial mechanisms of plant diversity and nitrogen addition on soil NCMF. We propose three hypotheses: (1) Plant diversity and nitrogen addition interactively affect soil NCMF (Cui et al., 2020; Delgado-Baquerizo et al., 2020; Li et al., 2020a,b; Freitag et al., 2023); (2) higher nitrogen addition increases soil microbial carbon limitation by reducing soil pH and thus affects soil NCMF (Ma et al., 2016; Yang et al., 2020); and (3) higher plant diversity positively influences soil organic matter and reduces microbial carbon limitation, thereby significantly increasing soil NCMF (Fornara and Tilman, 2008; Meyer et al., 2016; Chen et al., 2020).
2 Materials and methods
2.1 Study region
We conducted a microcosm experiment at the Jilin Songnen Grassland Ecosystem National Observation and Research Station, which is located in the Songnen grassland, China. The region is characterized by a meadow steppe dominated by Leymus chinensis, with widespread salinized soils. The climate is a continental monsoon climate, with a mean annual temperature ranging from 4.6 to 6.4°C and an average annual precipitation of 470.9 mm, 70% of which occurs between June and August (Cui et al., 2020).
2.2 Preparation of the plant species pool
We surveyed local plant species in the Songnen grassland and selected six native species that frequently occur in the area: Leymus chinensis, Hierochloe glabra, Lespedeza daurica, Vicia amoena, Carex duriuscula, and Kalimeris integrifolia. These species were used to prepare aseptic seedlings. The sterilized seeds were sown in plastic boxes (20 cm internal diameter, 30 cm depth) and placed in a climate-controlled greenhouse at Northeast Normal University (25°C, 70% relative humidity, 16-h daylight at 5000 lux illumination). Seedlings approximately 5 cm in height were selected and transplanted into microcosms (30 cm internal diameter, 40 cm depth). The homogenized meadow soil used in this study was collected from the Songnen grassland and sieved to remove roots and stones. To minimize interactions with the soil seed bank, the top 5 cm of soil was removed (García-Palacios et al., 2010).
2.3 Experimental design
The experiment comprised 120 pots arranged in 24 combinations, representing three nitrogen addition levels (N0: 0, N1: 5, and N2: 10 g NH4NO3–N m−2 yr.−1) and three levels of plant diversity (S1: 1 species, S3: 3 species, and S6: 6 species), with five replicates per combination. As soil NCMF tended to saturate at N rates ≥ 10 g m−2 yr.−1 in this grassland, nitrogen addition was applied at above three levels (Cui et al., 2020). We cultivated monocultures of all six plant species, as well as three-species combinations and mixtures of all six species, resulting in eight distinct plant compositions. Theoretically, there are 20 potential combination models for selecting three species from the six available; however, the experimental workload would be prohibitively large. Therefore, we randomly selected one representative combination model—L. chinensis + L. daurica + V. amoena—to exemplify the three-species combinations. While this selection does not encompass all possible tri-species combinations, it effectively represents the ecological responses associated with other plant assemblages across different families, genera, or functional groups.
Each pot was filled with 21.5 kg of homogenized meadow soil. The pH of the homogenized soil was 8.14, and the NH4+-N, NO3−-N, and total carbon (TC) concentrations were 2.53 mg kg−1, 0.21 mg kg−1, and 4.07 g kg−1, respectively. The soil was irrigated to adjust the gravimetric water content to match field conditions. Twelve seedlings selected were planted in each pot, ensuring an equal likelihood of neighboring species across microcosms (Fukami et al., 2001). Pots were randomly arranged in an open-ended greenhouse at the field station. The microcosms were buried with their top edges 3 cm above ground level to prevent surface runoff and left for 30 days to allow plant community establishment. Deionized water was applied weekly to ensure plant survival. In late September, the microcosms were mown to simulate local farming practices, and weeds were removed regularly throughout the experiment.
Nitrogen treatments were applied via the addition of chemically pure NH4NO3 (35% nitrogen content). The annual fertilization amount was divided into monthly applications during the growing season (mid-June to August) over 2 years. Each NH4NO3 portion was dissolved in deionized water and evenly sprayed into the pots. The control pots received equivalent amounts of deionized water to account for the water addition effect. In the first year, no plants or soils were harvested to allow for transient responses of the soil and plant properties. Sampling occurred at peak plant biomass in the second year.
2.4 Plant sample and soil collection and analyses
The plant and soil samples were collected in August of the second year (peak growing season). The above-ground biomass (AGB) and below-ground biomass (BGB) were harvested from each pot. The plant tissues were heated at 105°C for 30 min to halt metabolic activity, then dried at 65°C and weighed (Rillig et al., 2002).
Soil samples were collected from the upper 10 cm at five points in each pot and mixed into a composite sample. The soil was sieved through a 2 mm mesh to remove roots and divided into three parts for analysis. One part was stored at 4°C for measuring the soil water content (SWC), available nitrogen (NO₃−–N and NH₄+–N), soil net nitrogen mineralization rate (Rm), and net nitrogen nitrification rate (Rn). The second part was air-dried to analyze the soil pH, total nitrogen (TN), soil organic matter (SOM), and TC. The third part was stored at −80°C for qPCR and enzyme activity assays (LAP, NAG, alkaline phosphatase: ALP, and β-1,4-glucosidase: βG).
The soil pH was measured via a PHS-3E glass electrode (Leichi, Shanghai, China) in a 1:5 soil-to-water suspension. The SWC was determined by drying the samples at 105°C to a constant weight. Available nitrogen was quantified via a continuous flow analyzer (Futura, Alliance-AMS, France). Rm and Rn were measured during aerobic incubation (Hart et al., 1994). TN and TC were analyzed via an elemental analyzer (Isoprime 100, Isoprime Ltd., Manchester, United Kingdom). Soil organic carbon was analyzed using a total organic carbon analyzer (vario TOC cube, Elementar, Hanau, Germany) through dry combustion. According to Van Bemmelen factor, soil organic matter content (SOM) = soil organic carbon content × 1.724. Microbial biomass carbon was determined with the chloroform fumigation–extraction method (Cui et al., 2020). Enzyme activities were measured via a microplate fluorometric assay (TECAN Infinite F200, Tecan Group, Switzerland) (Cui et al., 2020).
2.5 DNA extraction and qPCR analysis
Soil DNA was extracted using the Power Soil DNA Isolation Kit (MoBio Laboratories, San Diego, CA, United States) according to the manufacturer’s instructions. The purified DNA was stored at −80°C after purification with a DNA purification kit and verification via 2% agarose gel electrophoresis. The abundance of 16S rRNA and 18S rRNA genes was determined using quantitative PCR (qPCR). PCR amplifications were conducted in triplicate 20 μL reaction mixtures, each containing 4 μL of 5 × FastPfu Buffer, 2 μL of dNTPs, 0.8 μL of each primer, 0.4 μL of FastPfu Polymerase, 10 ng of template DNA, and ultra-pure water to adjust the volume to 20 μL. The PCR protocol was as follows: initial denaturation at 95°C for 2 min, followed by 27 cycles of 95°C for 30 s, 55°C for 30 s, and 72°C for 40 s, with a final extension at 72°C for 10 min. The PCR products were then excised from 2% agarose gels and purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, United States) according to the manufacturer’s guidelines. The quantification of the abundances of nitrogen cycling microbial genes was performed by a StepOne™ Real-Time PCR System (Applied Biosystems, CA, United States) with a 20 μL reaction volume, including 2 μL DNA templates, 0.4 μL forward and reverse primers, and 10 μL Fast qPCR Master Mix (BBI, Canada). Tenfold serial dilutions of the linearized plasmid DNA were used to establish a standard curve for each gene. Copy numbers per ng DNA were transformed into gene copy numbers per gram of soil using the DNA amount retrieved per gram of soil (Ma et al., 2016). Standard curves were generated via tenfold serial dilutions of linearized plasmid DNA, and gene copy numbers were normalized per gram of soil. Meanwhile, the microbial community structure was assessed by the ratio of fungal gene copies and bacterial gene copies.
2.6 Soil microbial carbon limitation
Soil microbial carbon limitation was assessed using vector analysis of enzymatic stoichiometry (Chen et al., 2021). The vector length, calculated as , represents the degree of carbon limitation, with a greater vector length indicating a stronger limitation.
2.7 Assessing soil NCMF
Six soil variables related to nitrogen storage and cycling, including soil nifH, AOB amoA, nirK, nirS, Rm, and Rn, were measured in this study (Maestre et al., 2012; Delgado-Baquerizo et al., 2020; Osburn et al., 2021; Xu et al., 2021). These variables represent the processes of nitrogen cycling and available nutrient supply (Wang et al., 2023). The soil NCMF was calculated using an averaging approach, which has been widely used in multifunctionality analyses (Manning et al., 2018; Cui et al., 2020; Xu et al., 2021; Wang et al., 2023). To obtain the average soil NCMF index, we tested all variables for normal distribution using the Shapiro–Wilk test prior to analysis, and logarithm or square root transformation was performed when necessary. Nine soil variables were standardized and normalized using z-score transformations, and the average of these transformed values was calculated as the final result (Maestre et al., 2012; Valencia et al., 2018).
2.8 Statistical analyses
Two-way ANOVA was used to assess the individual and combined effects of plant diversity and nitrogen addition on ecosystem attributes (BGB, AGB, soil pH, SWC, TC, SOM, NO3−-N, NH4+-N, TN, the gene abundances of nifH, AOB amoA, nirK, nirS, NAG, LAP, ALP, βG, Rm, Rn, bacterial abundance, fungal abundance, fungi to bacteria ratio, microbial biomass carbon, microbial carbon limitation, and soil NCMF). Microbial gene abundance was log10 transformed before two-way ANOVA. p < 0.05 was considered to identify statistically significant differences. Tukey’s HSD post-hoc tests were used to evaluate the effects of nitrogen addition and plant diversity on soil NCMF. A piecewise structural equation model was used to test the direct and indirect effects of plant diversity and nitrogen addition on soil NCMF. The model assumes that plant diversity, nitrogen addition, and their interaction alter soil microbial carbon limitation levels via changes in soil pH, BGB, bacterial abundance, and SOM content and that they ultimately influence the soil NCMF. We constructed a prior model based on the known effects and potential relationships (Supplementary Figure S1). Fisher’s C statistic and the Akaike information criteria (AIC) were used to assess the goodness of-fit of the model (Shipley, 2013). All statistical analyses were carried out using R software (4.1.2, R Development Core Team, 2015).
3 Results
3.1 Effects of plant diversity and nitrogen addition on plants and soil properties
Two-way ANOVA revealed that both plant diversity and nitrogen addition significantly influenced SWC, TC, NH4+–N, and TN concentrations, AGB and BGB (Supplementary Table S1). Additionally, the interaction between plant diversity and nitrogen addition significantly affected on the soil pH, SWC, and TC, soil NO3−–N, NH4+–N, and TN concentrations (Supplementary Table S1). Compared with the N0 treatment, the N2 treatment significantly reduced the soil pH under the S3 and S6 treatments but significantly increased the SWC (Table 1). Furthermore, compared with the N0 treatment, the N1 and N2 treatments significantly elevated the soil TC and BGB under the S3 and S6 treatments (Table 1). Under N1 condition, S6 treatment significantly reduced soil pH than the S1 treatment (Table 1). Under N0 and N2 condition, S6 treatment significantly increased SWC than the S1 treatment (Table 1). Under N0, N1 and N2 condition, S6 treatment significantly increased SOM, TC, and BGB than the S1 treatment (Table 1). Under the S1, S3, and S6 conditions, the NO3−–N concentrations in the N2 treatment were significantly greater than those in the N0 treatment (Table 1). Under N0, N1 and N2 condition, S6 treatment significantly increased the soil NO3−–N concentrations than the S1 treatment (Table 1). Under the N2 condition, S6 treatment significantly increased the soil NH4+–N concentrations than the S1 treatment (Table 1).

Table 1. Effects of plant diversity (S) and nitrogen addition (N) on plant and soil properties in semi-arid grassland mesocosms.
3.2 Effects of plant diversity and nitrogen addition on soil extracellular and microbial carbon limitation
Two-way ANOVA showed that plant diversity and nitrogen addition significantly influenced the activities of LAP and ALP (Supplementary Table S1). Additionally, the interaction between plant diversity and nitrogen addition significantly affected on the activities of NAG, LAP, βG, ALP, and microbial carbon limitation (Supplementary Tables S1, S2). Under the S1 condition, the N2 treatment significantly reduced the activities of NAG and ALP. However, under the S3 and S6 conditions, N2 significantly increased their activity (Figures 1A,B). Under the S1, S3, and S6 conditions, the activities of βG and ALP in the N2 treatment were significantly greater than those in the N0 treatment (Figures 1C,D). Additionally, under the S1 treatment, the N2 treatment significantly increased soil microbial carbon limitation compared with the N0 treatment, but this effect was not observed under the S6 treatment (Figure 1E). Under N1 and N2 condition, S3 and S6 treatment significantly increased the activities of NAG than the S1 treatment (Figure 1A). Under N0, N1 and N2 condition, S3 and S6 treatment significantly increased the activities of LAP than the S1 treatment (Figure 1B). Under N0 and N2 condition, S6 treatment significantly increased the activities of βG than the S1 treatment (Figure 1C). Under N0, N1 and N2 condition, S6 treatment significantly increased the activities of ALP than the S1 treatment (Figure 1D). Under N1 and N2 condition, S3 and S6 treatment significantly reduced the soil microbial carbon limitation levels than the S1 treatment (Figure 1E).

Figure 1. Responses of NAG (A), LAP (B), βG (C), ALP activities (D) and soil microbial carbon limitation levels (E) to plant diversity and nitrogen addition in semi-arid grassland microcosm. Different lower- and upper-case letters indicate significant differences (p < 0.05) among different nitrogen addition and plant diversity treatment, respectively (Tukey’s HSD) in semi-arid grassland mesocosms. ns indicates that the difference is not significant. S1: 1 species, S3: 3 species, S6: 6 species. N0: without nitrogen addition, N1: 5 g N m−2 yr.−1, N2: 10 g N m−2 yr.−1.
3.3 Effects of plant diversity and nitrogen addition on soil microbial community properties
Two-way ANOVA revealed that both plant diversity and nitrogen addition significantly influenced soil bacterial abundance (Supplementary Table S1). Moreover, the interaction between plant diversity and nitrogen addition significantly affected on the soil bacterial abundance and microbial biomass carbon (Supplementary Table S1). There was no significant interaction between plant diversity and nitrogen addition on fungal abundance and the fungi to bacteria ratio (Supplementary Table S1). Compared with the N0 treatment, the N2 treatments significantly increased the soil bacterial abundance under the S1, S3 and S6 treatments (Figure 2A). Under N0 and N2 condition, S6 treatment significantly increased the soil bacterial abundance than the S1 treatment (Figure 2A). Compared with the N0 treatment, the N2 treatments have no significantly increased the fungal abundance, fungi to bacteria ratio, and microbial biomass carbon under the S1and S3 treatments (Figures 2B–D).

Figure 2. Responses of soil bacterial abundance (A), fungal abundance (B), fungi to bacteria ratio (C), and microbial biomass carbon (D) to plant diversity and nitrogen addition in semi-arid grassland microcosm. Different lower- and upper-case letters indicate significant differences (p < 0.05) among different nitrogen addition and plant diversity treatment, respectively (Tukey’s HSD) in semi-arid grassland mesocosms. ns indicates that the difference is not significant. S1: 1 species, S3: 3 species, S6: 6 species. N0: without nitrogen addition, N1: 5 g N m−2 yr.−1, N2: 10 g N m−2 yr.−1.
3.4 Effects of plant diversity and nitrogen addition on soil nitrogen cycling functions and NCMF
Two-way ANOVA showed that both plant diversity and nitrogen addition significantly influenced the abundances of nifH, nirK, nirS (Supplementary Table S1). Additionally, the interaction between plant diversity and nitrogen addition significantly affected on the abundances of nifH, AOB amoA, nirK, nirS, Rm, Rn, and soil NCMF (Supplementary Tables S1, S2). The abundance of nifH was significantly lower in the N2 treatment than in the N0 and N1 treatments (Figure 3A). Under the S1 condition, the N2 treatment significantly reduced the abundance of AOB amoA, nirK, nirS, and soil NCMF (Figures 3B–D, 4). However, under the S3 and S6 conditions, N2 significantly increased their abundance and soil NCMF (Figure 4). Under the N2 condition, S3 and S6 treatment significantly increased the abundances of nifH, AOB amoA, nirK, nirS, Rm, and Rn than the S1 treatment (Figures 3A–F). Under N1 and N2 condition, S3 and S6 treatment significantly increased the soil NCMF than the S1 treatment (Figure 4).

Figure 3. Responses of nifH (A), AOB amoA (B), nirK (C), nirS (D), Rm (E), and Rn (F) to plant diversity and nitrogen addition in semi-arid grassland microcosm. Different lower- and upper-case letters indicate significant differences (p < 0.05) among different nitrogen addition and plant diversity treatment, respectively (Tukey’s HSD) in semi-arid grassland mesocosms. ns indicates that the difference is not significant. S1: 1 species, S3: 3 species, S6: 6 species. N0: without nitrogen addition, N1: 5 g N m−2 yr.−1, N2: 10 g N m−2 yr.−1.

Figure 4. Responses of soil NCMF index to nitrogen addition and plant diversity. Different lower- and upper-case letters indicate significant differences (p < 0.05) among different nitrogen addition and plant diversity treatment, respectively (Tukey’s HSD) in semi-arid grassland mesocosms. S1: 1 species, S3: 3 species, S6: 6 species. N0: without nitrogen addition, N1: 5 g N m−2 yr.−1, N2: 10 g N m−2 yr.−1.
3.5 Indirect effects of plant diversity and nitrogen addition on soil NCMF
Piecewise structural equation modeling (SEM) demonstrated that soil NCMF was indirectly influenced by plant diversity and nitrogen addition through pathways involving BGB, soil pH, SOM, bacterial abundance, and microbial carbon limitation (Figure 5A). Plant diversity alleviated the negative effects of nitrogen addition on soil NCMF by increasing BGB and SOM, which reduced microbial carbon limitation and mitigated the pH reduction caused by nitrogen addition (Figure 5A). SOM displayed largest positive effects on soil NCMF (Figure 5B). However, soil microbial carbon limitation exhibited a highest negative effect on soil NCMF (Figure 5B).

Figure 5. The piecewise structural equation model was used to test the direct and indirect causality between plant diversity, nitrogen addition, and soil NCMF (A). The corresponding values of the solid line arrows are the width of normalized path system arrows reflect the size of normalized path coefficient. The blue and red arrows show significant positive and negative correlations, respectively (*p < 0.05, **p < 0.01 and ***p < 0.001). Dashed lines show non-significant relationships. The values above each variable represents the explanatory degree (R2) of each variable in the model. Standardized total effects (direct plus indirect effects) derived from the piecewise structural equation model (B).
4 Discussion
Our findings indicate that plant diversity and nitrogen addition interactively influence soil NCMF (Table 2), suggesting that the impact of nitrogen addition depends on the diversity of the plant community (Figure 4). Specifically, nitrogen addition significantly increased the soil NCMF under the S3 and S6 conditions but not under the S1 condition (Figure 4). This implies that greater plant diversity can mitigate the adverse effects of nitrogen addition on soil NCMF. In more diverse plant communities, nitrogen addition disrupts soil nitrogen cycling less severely than it does in less diverse communities.

Table 2. Two-way ANOVA testing the effects of plant diversity (S), nitrogen addition (N), and their interaction (S × N) on soil NCMF and soil microbial carbon limitation in semi-arid grassland microcosm.
4.1 Effects of plant diversity and nitrogen addition on soil NCMF
This study showed that nitrogen addition decreased the soil NCMF under low-diversity conditions (Figure 4); however, under high-diversity conditions, nitrogen addition significantly increased the soil NCMF, and the negative effects resulting from nitrogen addition in the low-diversity treatments were counteracted (Figure 4). The SEM clearly explained these effects (Figure 5A). Nitrogen enrichment can substantially alter soil properties, such as soil acidification (Table 1) and losses of base cations (Matson et al., 1999; Peñuelas et al., 2013). Furthermore, soil acidification-induced reductions in soil nitrogen cycling microbial activity and organic matter decomposition rates might overwhelm the positive effects of increased plant productivity under nitrogen addition (Song et al., 2021). Furthermore, the Songnen grassland in the semi-arid area was greatly limited by nitrogen (Shi et al., 2021); thus, improved nitrogen availability and supplies might have shifted soil nitrogen cycling microbes from nitrogen limitation to carbon limitation (Zhang et al., 2021). Finally, these decreases in nitrogen cycling function gene abundances have been known to be drivers of the decline in the soil nitrogen transformation rate (Ning et al., 2015; Tang et al., 2018), which likely caused the observed decrease in soil NCMF. This result is consistent with a previous study reporting that soil acidification may impede the N cycling process with nitrogen deposition (Tian and Niu, 2015).
As we hypothesized, the responses of the soil NCMF to nitrogen deposition were modified by plant diversity (Figures 4, 5). The SEM results showed that high plant diversity could offset the high nitrogen-induced decrease in soil NCMF through direct and indirect effects compared with low plant diversity (Figure 5), suggesting that plant diversity could significantly maintain and regulate nitrogen cycling in a grassland ecosystem. The direct pathway might be regulated by the changes in nitrogen concentration in roots, thus directly affecting soil nitrogen transformations because of the positive correlation between the nitrogen concentrations in roots and the rate of soil nitrogen turnover (Mueller et al., 2013).
The increase in plant below-ground biomass may be another possible mechanism by which high plant diversity offsets the negative effect of nitrogen addition on soil NCMF. It is well known that the input of plant below-ground biomass is an essential source of the grassland soil carbon pool (Reich et al., 2001; Li et al., 2017; Li et al., 2021). The present study found that high plant species diversity increased below-ground biomass, which is consistent with the previous results that there was a positive correlation between plant diversity and productivity (Reich et al., 2001). The increased soil carbon resources provided rich substrates and energy for microorganism growth and basal metabolism so that the abundance of nitrogen cycling genes and the extracellular enzyme activities associated with nitrogen cycling significantly increased, thus offsetting the negative effect of high nitrogen addition on soil NCMF under low diversity. The enhancement of soil enzyme activity in this study could be partly attributed to the increase in plant below-ground biomass (Table 1). Plant roots are not only the primary organs for nutrient uptake but also significant sources of soil extracellular enzymes. Under higher plant diversity communities, the root systems become more complex and extensive, leading to a higher input of root exudates into the soil (Reich et al., 2001; Li et al., 2017). These root exudates, which include organic acids, sugars, and amino acids, can directly promote the activity of soil extracellular enzymes (Yuan et al., 2020; Liu et al., 2021). Therefore, the impact of plant diversity on soil NCMF is not solely by the secretion of extracellular enzymes by microorganisms but also by the mediation of plant roots and their exudates. The relative abundance of corresponding nitrogen cycling functional genes has been demonstrated to regulate nitrogen transformation processes (Li et al., 2020a). For example, we found that nitrogen addition strongly increased the abundances of the AOB amoA, nirK, and nirS genes under high plant diversity conditions (Figure 3). Ammonia oxidation is the first rate-limiting step of nitrification that involves converting from ammonia to nitrite (Tang et al., 2016). In this study, the increased AOB amoA genes might accelerate the net nitrogen nitrification rate (Kowalchuk and Stephen, 2001; Osburn et al., 2021) and then enhance the soil NCMF (Figure 4C). These results supported our hypothesis, suggesting that higher plant diversity might benefit grassland ecosystem multifunctionality (Wang et al., 2019; Moi et al., 2021).
4.2 Interactive effects of plant diversity and nitrogen addition on soil NCMF
The observed interaction can be attributed to several underlying mechanisms. First, nitrogen enrichment can significantly alter soil properties (Peñuelas et al., 2013; Elrys et al., 2021), such as by inducing soil acidification (Table 1). A lower soil pH can directly inhibit the growth and activity of soil bacteria, which are often less tolerant of acidic conditions than fungi are (Zhang et al., 2008). Acidic conditions may favor fungal dominance over bacterial dominance (Hu et al., 2024). While fungi are generally more acid tolerant, they may decompose organic matter more slowly than bacteria do, potentially reducing overall nutrient cycling rates (Wang et al., 2011; Wang et al., 2023). Key nitrogen-cycling microbes, such as nitrifiers, are particularly sensitive to soil pH changes, and their reduced activity under acidic conditions can slow nitrogen transformation processes, such as nitrification and denitrification (Wu et al., 2024). Soil acidification can decrease the solubility of soil organic carbon, decreasing its availability to microbes, which can exacerbate microbial carbon limitations (Figure 1E), particularly in ecosystems where carbon availability is already low (Hallin et al., 2009; Cui et al., 2020; Cui et al., 2021; Song et al., 2021). Carbon-limited microbes may exhibit reduced enzymatic activity, slowing the decomposition of soil organic matter and the release of nutrients (Shu et al., 2011; Kivlin and Treseder, 2013). In some cases, high nitrogen availability can stimulate microbial activity, leading to faster decomposition of organic matter (Fornara and Tilman, 2008). However, carbon inputs (e.g., BGB) do not keep pace, which can deplete soil carbon stocks over time. In the nitrogen-limited Songnen grasslands of semi-arid regions (Shi et al., 2021), improved nitrogen availability may shift soil microbial communities from nitrogen limitation to carbon limitation (Zhang et al., 2008). Under soil carbon limitation, microbes may allocate more resources to acquire carbon (e.g., producing extracellular enzymes such as βG) (Figure 1C), potentially altering nutrient cycling dynamics (Hu et al., 2024).
Second, plant diversity may counteract microbial carbon limitation by increasing BGB (Figure 5). Diverse root structures and carbon-rich root exudates, along with increased organic matter from root turnover, provide additional carbon sources for soil microbes, alleviating carbon limitations (Yuan et al., 2020; Liu et al., 2021). Higher plant diversity leads to increased root turnover, contributing to SOM accumulation and providing a continuous supply of carbon for microbial communities (Figure 5) (Freitag et al., 2023). Furthermore, diverse plant communities may utilize nitrogen more efficiently, reducing soil acidification. Increased SOM not only supplies carbon but also enhances soil structure, water retention (Table 1), and microbial habitats (Chen et al., 2020; Cheng et al., 2022). This enhanced microbial activity can improve organic matter decomposition, increasing the amount of carbon released into the soil (Chen and Chen, 2021). With sufficient carbon, microbes can process nitrogen more efficiently (Tang et al., 2016), potentially increasing soil net nitrogen mineralization and nitrification rates. In diverse plant communities, species may utilize nitrogen in different ways or at different times, leading to more efficient nitrogen use overall (Maestre et al., 2011; Zhang et al., 2021; Xu et al., 2022). This reduces excess soil nitrogen, mitigating the negative effects of nitrogen addition. For example, deep-rooted and shallow-rooted plants may access nitrogen from different soil layers, preventing nitrogen leaching (Oelmann et al., 2011). The insurance hypothesis further suggests that diverse communities are more resilient to stressors such as excess nitrogen, as compensatory mechanisms among species stabilize ecosystem functions, including nitrogen cycling (Xu et al., 2022). Diverse plant communities may also increase soil microbial abundance (Figures 2A,B), fostering robust nitrogen cycling processes that can accommodate added nitrogen without becoming unbalanced (Lucas-Borja and Delgado-Baquerizo, 2019; Wang et al., 2023). This creates a positive feedback loop, as healthy nitrogen cycling supports plant growth. The buffering effect of plant diversity may depend on functional diversity rather than species richness alone (Liu et al., 2021). For example, legumes, which fix atmospheric nitrogen, may interact differently with added nitrogen than nonlegumes do. Future research should explore these interactions further (Suter et al., 2015; Singh et al., 2018).
4.3 Implications for ecosystem management and climate change mitigation
These findings have significant implications for ecosystem management and climate change mitigation. In the context of global nitrogen deposition and the increasing use of nitrogen fertilizers, understanding the role of plant diversity in regulating soil NCMF is essential for predicting soil nitrogen dynamics and formulating sustainable land management strategies. Our results demonstrate that maintaining or enhancing plant diversity in terrestrial ecosystems can serve as a buffer against the adverse effects of nitrogen addition on soil NCMF. This study not only underscores the interactive effects of plant diversity and nitrogen addition on soil NCMF but also elucidates the underlying mechanisms driving these interactions. These insights highlight the critical need to incorporate plant diversity into strategies aimed at managing soil nitrogen dynamics amidst changing nitrogen regimes. Future research should focus on investigating the long-term impacts of plant diversity on soil NCMF, as well as the potential feedback mechanisms between soil nitrogen cycling and ecosystem productivity across diverse environmental scenarios. Such efforts will further refine our understanding of ecosystem resilience and inform more effective management practices in the face of global environmental change.
5 Conclusion
Our study provides robust evidence that plant diversity plays a pivotal role in mitigating the adverse effects of nitrogen addition to soil NCMF. Specifically, we demonstrate that the interactive effects of plant diversity and nitrogen addition on soil NCMF are mediated through a combination of indirect pathways, including changes in BGB, SOM, soil pH, bacterial abundance, and microbial carbon limitation. These mechanisms collectively highlight the multifaceted ways in which plant diversity buffers the negative impacts of nitrogen enrichment on soil nitrogen cycling processes. These findings provide novel insights into the mechanisms by which plant diversity and nitrogen addition influence soil NCMF in semi-arid grasslands, offering valuable implications for the understanding and sustainable management of grassland ecosystems under global climate change.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
ZL: Resources, Visualization, Data curation, Project administration, Validation, Conceptualization, Investigation, Supervision, Writing – original draft, Methodology, Software, Formal analysis, Writing – review & editing. XL: Supervision, Writing – review & editing. MZ: Writing – review & editing. BZ: Writing – review & editing. TZ: Writing – review & editing. CW: Writing – review & editing. FX: Resources, Visualization, Validation, Funding acquisition, Writing – review & editing, Supervision.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This research was financially supported by the National Natural Science Foundation of China (31670524).
Acknowledgments
We are deeply grateful to Xiaote Li and Yu Bi for their invaluable assistance in collecting soil and plant samples during field investigations.
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.
Generative AI statement
The author(s) declare that no Gen AI was used in the creation of this manuscript.
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.2025.1596515/full#supplementary-material
References
Byrnes, J. E. K., Gamfeldt, L., Isbell, F., Lefcheck, J. S., Griffin, J. N., Hector, A., et al. (2014). Investigating the relationship between biodiversity and ecosystem multifunctionality: challenges and solutions. Methods Ecol. Evol. 5, 111–124. doi: 10.1111/2041-210x.12143
Chen, X., and Chen, H. Y. H. (2021). Plant mixture balances terrestrial ecosystem C:N:P stoichiometry. Nat. Commun. 12:4562. doi: 10.1038/s41467-021-24889-w
Chen, X., Chen, H. Y. H., Searle, E. B., Chen, C., and Reich, P. B. (2020). Negative to positive shifts in diversity effects on soil nitrogen over time. Nature Sustain. 4, 225–232. doi: 10.1038/s41893-020-00641-y
Chen, H., Tang, J., Sun, X., Ma, K., Chen, H., and Li, D. (2021). Topography modulates effects of nitrogen deposition on microbial resource limitation in a nitrogen-saturated subtropical forest. Forest Ecosystems 8, 914–922. doi: 10.1186/s40663-021-00341-9
Cheng, B., Liu, H., Bai, J., and Li, J. (2022). Soil fungal composition drives ecosystem multifunctionality after long-term field nitrogen and phosphorus addition in alpine meadows on the Tibetan plateau. Plan. Theory 11:2893. doi: 10.3390/plants11212893
Cui, H., Sun, W., Delgado-Baquerizo, M., Song, W., Ma, J.-Y., Wang, K., et al. (2020). Phosphorus addition regulates the responses of soil multifunctionality to nitrogen over-fertilization in a temperate grassland. Plant Soil 473, 73–87. doi: 10.1007/s11104-020-04620-2
Cui, H., Sun, W., Delgado-Baquerizo, M., Song, W., Ma, J. Y., Wang, K., et al. (2021). Cascading effects of N fertilization activate biologically driven mechanisms promoting P availability in a semi-arid grassland ecosystem. Funct. Ecol. 35, 1001–1011. doi: 10.1111/1365-2435.13773
Delgado-Baquerizo, M., Reich, P. B., Trivedi, C., Eldridge, D. J., Abades, S., Alfaro, F. D., et al. (2020). Multiple elements of soil biodiversity drive ecosystem functions across biomes. Nature Ecol. Evolution 4, 210–220. doi: 10.1038/s41559-019-1084-y
Elrys, A. S., Ali, A., Zhang, H., Cheng, Y., Zhang, J., Cai, Z. C., et al. (2021). Patterns and drivers of global gross nitrogen mineralization in soils. Glob. Chang. Biol. 27, 5950–5962. doi: 10.1111/gcb.15851
Fornara, D. A., and Tilman, D. (2008). Plant functional composition influences rates of soil carbon and nitrogen accumulation. J. Ecol. 96, 314–322. doi: 10.1111/j.1365-2745.2007.01345.x
Freitag, M., Hölzel, N., Neuenkamp, L., van der Plas, F., Manning, P., Abrahão, A., et al. (2023). Increasing plant species richness by seeding has marginal effects on ecosystem functioning in agricultural grasslands. J. Ecol. 111, 1968–1984. doi: 10.1111/1365-2745.14154
Fukami, T., Naeem, S., and Wardle, D. A. (2001). On similarity among local communities in biodiversity experiments. Oikos 95, 340–348. doi: 10.1034/j.1600-0706.2001.950216.x
García-Palacios, P., Maestre, F. T., and Gallardo, A. (2010). Soil nutrient heterogeneity modulates ecosystem responses to changes in the identity and richness of plant functional groups. J. Ecol. 99, 551–562. doi: 10.1111/j.1365-2745.2010.01765.x
Hallin, S., Jones, C. M., Schloter, M., and Philippot, L. (2009). Relationship between N-cycling communities and ecosystem functioning in a 50-year-old fertilization experiment. ISME J. 3, 597–605. doi: 10.1038/ismej.2008.128
Hart, S. C., Stark, J. M., Davidson, E. A., and Firestone, M. K. (1994). Nitrogen mineralization, immobilization, and nitrification.Methods of Soil Analysis: Part 2 Microbiological and Biochemical Properties, 5, 985–1018.
He, K., Huang, Y., Qi, Y., Sheng, Z., and Chen, H. (2021). Effects of nitrogen addition on vegetation and soil and its linkages to plant diversity and productivity in a semi-arid steppe. Sci. Total Environ. 778:146299. doi: 10.1016/j.scitotenv.2021.146299
Hu, Z., Delgado-Baquerizo, M., Fanin, N., Chen, X., Zhou, Y., Du, G., et al. (2024). Nutrient-induced acidification modulates soil biodiversity-function relationships. Nat. Commun. 15:2858. doi: 10.1038/s41467-024-47323-3
Kivlin, S. N., and Treseder, K. K. (2013). Soil extracellular enzyme activities correspond with abiotic factors more than fungal community composition. Biogeochemistry 117, 23–37. doi: 10.1007/s10533-013-9852-2
Kowalchuk, G. A., and Stephen, J. R. (2001). Ammonia-oxidizing bacteria: A model for molecular microbial ecology. Ann. Rev. Microbiol. 55, 485–529. doi: 10.1146/annurev.micro.55.1.485
Li, G. Y., Cai, J. T., Song, X. X., Pan, X. B., Pan, D. F., Jiang, S. C., et al. (2021). Herbivore grazing mitigates the negative effects of nitrogen deposition on soil organic carbon in low-diversity grassland. J. Appl. Ecol. 59, 483–491. doi: 10.1111/1365-2664.14066
Li, J., Wang, G., Yan, B., and Liu, G. (2020a). The responses of soil nitrogen transformation to nitrogen addition are mainly related to the changes in functional gene relative abundance in artificial Pinus tabulaeformis forests. Sci. Total Environ. 723:137679. doi: 10.1016/j.scitotenv.2020.137679
Li, H., Yang, S., Xu, Z. W., Yan, Q. Y., Li, X. B., van Nostrand, J. D., et al. (2017). Responses of soil microbial functional genes to global changes are indirectly influenced by aboveground plant biomass variation. Soil Biol. Biochem. 104, 18–29. doi: 10.1016/j.soilbio.2016.10.009
Li, J., Yang, C., Zhou, H., and Shao, X. (2020b). Responses of plant diversity and soil microorganism diversity to water and nitrogen additions in the Qinghai-Tibetan plateau. Global Ecol. Conservation 22:e01003. doi: 10.1016/j.gecco.2020.e01003
Liu, X., Shi, X., and Zhang, S. (2021). Soil abiotic properties and plant functional diversity co-regulate the impacts of nitrogen addition on ecosystem multifunctionality in an alpine meadow. Sci. Total Environ. 780:146476. doi: 10.1016/j.scitotenv.2021.146476
Lu, M., Yang, Y., Luo, Y., Fang, C., Zhou, X., Chen, J., et al. (2010). Responses of ecosystem nitrogen cycle to nitrogen addition: a meta-analysis. New Phytol. 189, 1040–1050. doi: 10.1111/j.1469-8137.2010.03563.x
Lucas-Borja, M. E., and Delgado-Baquerizo, M. (2019). Plant diversity and soil stoichiometry regulates the changes in multifunctionality during pine temperate forest secondary succession. Sci. Total Environ. 697:134204. doi: 10.1016/j.scitotenv.2019.134204
Ma, W., Jiang, S., Assemien, F., Qin, M., Ma, B., Xie, Z., et al. (2016). Response of microbial functional groups involved in soil N cycle to N, P and NP fertilization in Tibetan alpine meadows. Soil Biol. Biochem. 101, 195–206. doi: 10.1016/j.soilbio.2016.07.023
Maestre, F. T., Castillo-Monroy, A. P., Bowker, M. A., and Ochoa-Hueso, R. (2011). Species richness effects on ecosystem multifunctionality depend on evenness, composition and spatial pattern. J. Ecol. 100, 317–330. doi: 10.1111/j.1365-2745.2011.01918.x
Maestre, F. T., Quero, J. L., Gotelli, N. J., Escudero, A., Ochoa, V., Delgado-Baquerizo, M., et al. (2012). Plant species richness and ecosystem multifunctionality in global drylands. Science 335, 214–218. doi: 10.1126/science.1215442
Manning, P., van der Plas, F., Soliveres, S., Allan, E., Maestre, F. T., Mace, G., et al. (2018). Redefining ecosystem multifunctionality. Nature Ecol. Evolution 2, 427–436. doi: 10.1038/s41559-017-0461-7
Meyer, S. T., Ebeling, A., Eisenhauer, N., Hertzog, L., Hillebrand, H., Milcu, A., et al. (2016). Effects of biodiversity strengthen over time as ecosystem functioning declines at low and increases at high biodiversity. Ecosphere 7:e01619. doi: 10.1002/ecs2.1619
Moi, D. A., Romero, G. Q., Antiqueira, P. A. P., Mormul, R. P., de Mello, F. T., and Bonecker, C. C. (2021). Multitrophic richness enhances ecosystem multifunctionality of tropical shallow lakes. Funct. Ecol. 35, 942–954. doi: 10.1111/1365-2435.13758
Matson, P. A., McDowell, W. H., Townsend, A. R., and Vitousek, P. M. (1999). The globalization of N deposition: ecosystem consequences in tropical environments. Biogeochemistry, 46, 67–83. doi: 10.1007/Bf01007574
Mueller, K. E., Hobbie, S. E., Tilman, D., and Reich, P. B. (2013). Effects of plant diversity, N fertilization, and elevated carbon dioxide on grassland soil N cycling in a long-term experiment. Glob. Chang. Biol. 19, 1249–1261. doi: 10.1111/gcb.12096
Ning, Q. S., Gu, Q., Shen, J. P., Lv, X. T., Yang, J. J., Zhang, X. M., et al. (2015). Effects of nitrogen deposition rates and frequencies on the abundance of soil nitrogen-related functional genes in temperate grassland of northern China. J. Soils Sediments 15, 694–704. doi: 10.1007/s11368-015-1061-2
Oelmann, Y., Buchmann, N., Gleixner, G., Habekost, M., Roscher, C., Rosenkranz, S., et al. (2011). Plant diversity effects on aboveground and belowground N pools in temperate grassland ecosystems: development in the first 5 years after establishment. Glob. Biogeochem. Cycles 25:869. doi: 10.1029/2010gb003869
Osburn, E. D., Badgley, B. D., Strahm, B. D., Aylward, F. O., and Barrett, J. E. (2021). Emergent properties of microbial communities drive accelerated biogeochemical cycling in disturbed temperate forests. Ecology 102:e03553. doi: 10.1002/ecy.3553
Peñuelas, J., Poulter, B., Sardans, J., Ciais, P., van der Velde, M., Bopp, L., et al. (2013). Human-induced nitrogen–phosphorus imbalances alter natural and managed ecosystems across the globe. Nat. Commun. 4:2934. doi: 10.1038/ncomms3934
R Development Core Team. (2015). R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna.
Reich, P. B., Knops, J., Tilman, D., Craine, J., Ellsworth, D., Tjoelker, M., et al. (2001). Plant diversity enhances ecosystem responses to elevated CO2 and nitrogen deposition. Nature 411:824. doi: 10.1038/35081122
Rillig, M., Wright, S., and Eviner, V. (2002). The role of arbuscular mycorrhizal fungi and glomalin in soil aggregation: comparing effects of five plant species. Plant Soil 238, 325–333. doi: 10.1023/A:1014483303813
Shah, V., Le Roux, X., Schmid, B., Poly, F., Barnard, R. L., Niklaus, P. A., et al. (2013). Soil environmental conditions and microbial build-up mediate the effect of plant diversity on soil nitrifying and denitrifying enzyme activities in temperate grasslands. PLoS One 8:e61069. doi: 10.1371/journal.pone.0061069
Shi, Y., Wang, J., Ao, Y., Han, J., Guo, Z., Liu, X., et al. (2021). Responses of soil N2O emissions and their abiotic and biotic drivers to altered rainfall regimes and co-occurring wet N deposition in a semi-arid grassland. Glob. Chang. Biol. 27, 4894–4908. doi: 10.1111/gcb.15792
Shipley, B. (2013). The AIC model selection method applied to path analytic models compared using a d‐separation test. Ecology, 94, 560–564. doi: 10.1890/12-0976.1
Shu, W., Pablo, G. P., Jun, Y., and Danfeng, H. (2011). Abundance and diversity of nitrogen-fixing bacteria in rhizosphere and bulk paddy soil under different duration of organic management. World J. Microbiol. Biotechnol. 28, 493–503. doi: 10.1007/s11274-011-0840-1
Singh, A. K., Rai, A., Banyal, R., Chauhan, P. S., and Singh, N. (2018). Plant community regulates soil multifunctionality in a tropical dry forest. Ecol. Indic. 95, 953–963. doi: 10.1016/j.ecolind.2018.08.030
Song, J., Xia, J., Hui, D., Zheng, M., Wang, J., Ru, J., et al. (2021). Plant functional types regulate non-additive responses of soil respiration to 5-year warming and nitrogen addition in a semi-arid grassland. Funct. Ecol. 35, 2593–2603. doi: 10.1111/1365-2435.13902
Suter, M., Connolly, J., Finn, J. A., Loges, R., Kirwan, L., Sebastià, M. T., et al. (2015). Nitrogen yield advantage from grass–legume mixtures is robust over a wide range of legume proportions and environmental conditions. Glob. Chang. Biol. 21, 2424–2438. doi: 10.1111/gcb.12880
Tang, Y. Q., Yu, G. R., Zhang, X. Y., Wang, Q. F., Ge, J. P., and Liu, S. (2018). Changes in nitrogen-cycling microbial communities with depth in temperate and subtropical forest soils. Appl. Soil Ecol. 124, 218–228. doi: 10.1016/j.apsoil.2017.10.029
Tang, Y., Zhang, X., Li, D., Wang, H., Chen, F., Fu, X., et al. (2016). Impacts of nitrogen and phosphorus additions on the abundance and community structure of ammonia oxidizers and denitrifying bacteria in Chinese fir plantations. Soil Biol. Biochem. 103, 284–293. doi: 10.1016/j.soilbio.2016.09.001
Tian, D. S., and Niu, S. L. (2015). A global analysis of soil acidification caused by nitrogen addition. Environ. Res. Lett. 10, 1714–1721. doi: 10.1088/1748-9326/10/2/024019
Turnbull, L. A., Levine, J. M., Loreau, M., Hector, A., and Lambers, J. H. R. (2012). Coexistence, niches and biodiversity effects on ecosystem functioning. Ecol. Lett. 16, 116–127. doi: 10.1111/ele.12056
Valencia, E., Gross, N., Quero, J. L., Carmona, C. P., Ochoa, V., Gozalo, B., et al. (2018). Cascading effects from plants to soil microorganisms explain how plant species richness and simulated climate change affect soil multifunctionality. Glob. Chang. Biol. 24, 5642–5654. doi: 10.1111/gcb.14440
van Ruijven, J., and Berendse, F. (2005). Diversity-productivity relationships: initial effects, long-term patterns, and underlying mechanisms. Proc Natl Acad Sci U S A. 102, 695–700. doi: 10.1073/pnas.0407524102
van der Plas, F. (2019). Biodiversity and ecosystem functioning in naturally assembled communities. Biol. Rev. 94, 1220–1245. doi: 10.1111/brv.12499
Wang, L., Delgado-Baquerizo, M., Wang, D. L., Isbell, F., Liu, J., Feng, C., et al. (2019). Diversifying livestock promotes multidiversity and multifunctionality in managed grasslands. Proc. Natl. Acad. Sci. USA 116, 6187–6192. doi: 10.1073/pnas.1807354116
Wang, P., Edwards, I. P., Zak, D. R., Kellner, H., Eisenlord, S. D., and Pregitzer, K. S. (2011). Simulated atmospheric N deposition alters fungal community composition and suppresses Ligninolytic gene expression in a northern hardwood Forest. PLoS One 6:e20421. doi: 10.1371/journal.pone.0020421
Wang, C., Li, X. A., Hu, Y. X., Zheng, R. L., and Hou, Y. H. (2023). Nitrogen addition weakens the biodiversity-multifunctionality relationships across soil profiles in a grassland assemblage. Agric. Ecosyst. Environ. 342:108241. doi: 10.1016/j.agee.2022.108241
Wu, J., Wang, H., Li, G., and Chen, N. (2024). Effects of nitrogen deposition on soil nitrogen fractions and enzyme activities in wet meadow of the Qinghai-Tibet plateau. Sci. Rep. 14:31848. doi: 10.1038/s41598-024-83285-8
Xu, Y., Dong, K., Jiang, M., Liu, Y., He, L., Wang, J., et al. (2022). Soil moisture and species richness interactively affect multiple ecosystem functions in a microcosm experiment of simulated shrub encroached grasslands. Sci. Total Environ. 803:149950. doi: 10.1016/j.scitotenv.2021.149950
Xu, H., Yu, M., and Cheng, X. (2021). Abundant fungal and rare bacterial taxa jointly reveal soil nutrient cycling and multifunctionality in uneven-aged mixed plantations. Ecol. Indic. 129:107932. doi: 10.1016/j.ecolind.2021.107932
Yang, Y., Chai, Y., Xie, H., Zhang, L., Zhang, Z., Yang, X., et al. (2023). Responses of soil microbial diversity, network complexity and multifunctionality to three land-use changes. Sci. Total Environ. 859:160255. doi: 10.1016/j.scitotenv.2022.160255
Yang, K., Luo, S., Hu, L., Chen, B., Xie, Z., Ma, B., et al. (2020). Responses of soil ammonia-oxidizing bacteria and archaea diversity to N, P and NP fertilization: relationships with soil environmental variables and plant community diversity. Soil Biol. Biochem. 145:107795. doi: 10.1016/j.soilbio.2020.107795
Yuan, Z., Ali, A., Ruiz-Benito, P., Jucker, T., Mori, A. S., Wang, S., et al. (2020). Above- and below-ground biodiversity jointly regulate temperate forest multifunctionality along a local-scale environmental gradient. J. Ecol. 108, 2012–2024. doi: 10.1111/1365-2745.13378
Zhang, N., Wan, S., Li, L., Bi, J., Zhao, M., and Ma, K. (2008). Impacts of urea N addition on soil microbial community in a semi-arid temperate steppe in northern China. Plant Soil 311, 19–28. doi: 10.1007/s11104-008-9650-0
Keywords: nitrogen deposition, plant diversity, soil microbial carbon limitation, soil nitrogen cycling multifunctionality, global change
Citation: Li Z, Liu X, Zhang M, Zhang B, Zhang T, Wagg C and Xing F (2025) High plant diversity alleviates the negative effects of nitrogen deposition on soil nitrogen cycling multifunctionality. Front. Microbiol. 16:1596515. doi: 10.3389/fmicb.2025.1596515
Edited by:
Wenguang Sun, University of Nebraska-Lincoln, United StatesReviewed by:
Leilei Ding, Guizhou Academy of Agricultural Sciences, ChinaJunxiao Pan, Chinese Academy of Sciences (CAS), China
Tingting Ren, Nanjing Forestry University, China
Copyright © 2025 Li, Liu, Zhang, Zhang, Zhang, Wagg and Xing. 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: Fu Xing, eGluZ2Y1MjJAMTI2LmNvbQ==