Front. Plant Sci.Frontiers in Plant ScienceFront. Plant Sci.1664-462XFrontiers Media S.A.10.3389/fpls.2022.860590Plant ScienceOriginal ResearchN Addition Overwhelmed the Effects of P Addition on the Soil C, N, and P Cycling Genes in Alpine Meadow of the Qinghai-Tibetan PlateauXiaoJiannan1DongShikui2*ShenHao2LiShuai3WessellKelly4LiuShiliang1LiWei1ZhiYangliu1MuZhiyuan1LiHongbo51School of Environment, Beijing Normal University, Beijing, China2School of Grassland Sciences, Beijing Forestry University, Beijing, China3College of Resource and Environment, Shanxi Agricultural University, Taigu, China4Tompkin Cortland Community College, Ithaca, NY, United States5Institute of Environment and Sustainable Development in Agriculture, Chinese Academy of Agricultural Sciences, Beijing, China
Edited by: Huakun Zhou, Key Laboratory of Restoration Ecology in Cold Regions, Northwest Institute of Plateau Biology (CAS), China
Reviewed by: Jin Hua Li, Lanzhou University, China; Lu Wen, Inner Mongolia University, China; Xukun Su, Research Center for Eco-environmental Sciences (CAS), China
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.
Although human activities have greatly increased nitrogen (N) and phosphorus (P) inputs to the alpine grassland ecosystems, how soil microbial functional genes involved in nutrient cycling respond to N and P input remains unknown. Based on a fertilization experiment established in an alpine meadow of the Qinghai-Tibetan Plateau, we investigated the response of the abundance of soil carbon (C), N, and P cycling genes to N and P addition and evaluated soil and plant factors related to the observed effects. Our results indicated that the abundance of C, N, and P cycling genes were hardly affected by N addition, while P addition significantly increased most of them, suggesting that the availability of P plays a more important role for soil microorganisms than N in this alpine meadow ecosystem. Meanwhile, when N and P were added together, the abundance of C, N, and P cycling genes did not change significantly, indicating that the promoting effects of P addition on microbial functional genes abundances were overwhelmed by N addition. The Mantel analysis and the variation partitioning analysis revealed the major role of shoot P concentration in regulating the abundance of C, N, and P cycling genes. These results suggest that soil P availability and plant traits are key in governing C, N, and P cycling genes at the functional gene level in the alpine grassland ecosystem.
Nitrogen (N) is a common nutrient element limiting plant productivity in most terrestrial ecosystems (Vitousek and Howarth, 1991). Atmospheric N deposition has increased substantially over the last decades due to anthropogenic activities (Liu et al., 2015) and has dramatically altered the characteristics of soils and plants (Pan et al., 2014). For example, studies have shown soil acidification and eutrophication, reduced species diversity, and increased primary productivity under N deposition (Jing et al., 2016; She et al., 2018; Li et al., 2019). Phosphorus (P) is another key nutrient for plant growth and ecosystem health, but P has not increased by anthropogenic activity to the same degree as N (Penuelas et al., 2013). N deposition may result in the limitation of P by influencing the nutrient balance in the soil (Cui et al., 2020).
Soil microorganisms play a major role in driving carbon (C), N, and P biogeochemical cycles and in responses to environmental changes (Lynn et al., 2016; Smercina et al., 2021). Functional genes, encoding various enzymes, not only provide important information about the genetic potentials for specific metabolic processes in soil but also reflect the diversity and activity of microbial communities (He et al., 2012; Turner et al., 2019). Nutrient availability is a key player in shaping soil microbial community (Lauber et al., 2008). N and P addition that increase N and P availability are responsible for changing the abundance and the community structure of soil microorganisms (Wang et al., 2017). The response of soil microbes to N addition has been widely reported across ecosystems. Generally, N addition tends to decrease the abundance of microbial community due to reduced belowground C allocation under N addition (Frey et al., 2004; Bae et al., 2015). N addition decreased the gene abundance of ammoxidation (AOA) and N fixation (nifH) (He et al., 2010) and increased the abundance of narG, nirK, nirS, and nosZ genes (denitrification) (Kastl et al., 2015). In a temperate steppe ecosystem, 5-year of N addition enhanced the abundance of the AOB gene by 11-fold and decreased the abundance of nifH, chiA, AOA, nirS, nirK, and nosZ genes by 30–50% (Zhang et al., 2019). N addition (exceeding 40 kg N ha–1 year–1) notably enhanced the abundance of C degradation genes in an alpine steppe (Chen et al., 2020). Similar to N addition, P addition can alter the structure of microbial functional genes (Tan et al., 2013). Both increased and decreased soil microbial diversity under P fertilization have been reported in various ecosystems (Chhabra et al., 2013; Tan et al., 2013). These inconsistent results may be due to the fact that P input may not be the only factor responsible for changes in the microbial community and that other environmental factors such as soil characteristics and climate may also be involved in the process (Tan et al., 2013).
It is noteworthy that N and P additions not only put a direct influence on microbial community but also indirectly affected microbial abundances through vegetation (Li et al., 2017; She et al., 2018). Plant traits can alter soil abiotic factors that regulate the abundance, composition, and activity of soil microbial communities (Abalos et al., 2019). N and P additions can influence plant growth, leading to the alteration of plant residue, such as litter and root exudate, which directly affect the structure and the composition of the microbial community (She et al., 2018). For example, with the increase of belowground plant biomass after P addition, the plant C distribution to soil was enhanced, hence improving microbial biomass and microbe phenotype (Dong J. et al., 2020). Aboveground plant biomass indirectly stimulated the abundance of soil microbial functional genes by increasing soil dissolved organic carbon (Li et al., 2017). To date, the impacts of N and P additions, as well as their complex interactive influences on the microbial mediation of C, N, and P biogeochemical cycling in alpine grassland ecosystems, remain elusive, despite the fact that alpine grassland ecosystems are globally vital in providing ecosystems services and functions such as C/N storage, biodiversity conservation, water conservation, and nutrient cycling (Dong et al., 2007).
The Qinghai-Tibetan Plateau (QTP) is the largest plateau in China and the highest in the world, and it is described as a “Global Change Warning Area” and the “Global Climate Stove” due to its sensitivity to global change (Dong et al., 2007). This region has been experiencing a significant increase in N deposition, approximately 8 kg N ha–1 year–1 (Lü and Tian, 2007). Due to its low contents and decomposition rates, N is always considered to be the critical limiting factor of alpine ecosystems on the QTP (Li et al., 2019). However, it has also been reported that P, not N, is the limiting factor in the QTP’s alpine grassland ecosystem (Jing et al., 2016). Therefore, it has not been determined whether the alpine ecosystems on the QTP are limited by either N, P, or N and P. N and P inputs could affect ecosystem function significantly (Jing et al., 2016). Elevated N deposition has significantly influenced soil properties and vegetation characteristics of QTP’s alpine grasslands (Jiang et al., 2013; Ma et al., 2016; Li et al., 2019; Chen et al., 2020). P addition suppressed topsoil microbial activity and increased plant primary productivity in an alpine grassland (Jing et al., 2016). Alpine meadow, which is widely distributed on the QTP, accounts for over 40% of the QTP’s territory (Kang et al., 2007) and provides essential ecosystem functional services for millions of people living there and downstream (Dong et al., 2007). However, to our knowledge, few researchers have documented the response of soil C, N, and P cycling genes to N and P addition in QTP’s alpine meadow, and the environmental factors that shape these genes in soil have not yet been established. These uncertainties will restrict our ability to predict the response of the biogeochemistry cycle to future environmental changes in the QTP’s alpine ecosystems.
Herein, we measured the abundances of soil microbial functional genes for C fixation (smtA, rbcL, mct, pccA, acsE, acsA, aclB, frdA, korA, and accA), C degradation (amyA, apu, manA, abfA, CDH, chiA, exg, gam, glx, IsoP, lig, mnp, exoPG, and xylA), N cycling (amoA, amoB, hzsB, nirK, nirS, nosZ, nxrA, nifH, napA, gdh, and UreC), and P cycling (phnK, pqqC, bpp, phoD, phoX, ppx, gmGDH, and emGDH). We also measured plant biomass and nutrients and soil nutrients to analyze their relationships with C, N, and P cycling genes under N and P additions. We hypothesized that (1) combined N and P additions would promote the functional capabilities of C, N, and P cycling genes in this alpine meadow ecosystem by elevating nutrient availabilities; (2) its promotion effect would be greater than that of N alone addition due to the alleviation of P limitation; and (3) N and P additions not only put direct influence on microbial functional genes but also indirectly affect them through vegetation.
Materials and MethodsSite Description and Experimental Design
This study was carried out in Xihai Town of Haiyan County (36°56′N, 100°57′E, 3100 m ASL) of the Qinghai Province, China (Figure 1). The major vegetation in Haiyan County was classified as alpine meadow, and the dominant plant species include Agropyron cristatum, Poa pratensis, Elymus dahuricus, and Stipa capillata (Li et al., 2019). The mean annual temperature and mean annual precipitation at Haiyan County were 1.4°C and 330–370 mm, respectively. The mean daily temperature and precipitation during the growing season are shown in Supplementary Figure 1. The soil was classified as Mat-Gryic Cambisols (Chinese National Soil Survey), with sand, silt, and clay contents of 16.75, 58.7.3, and 24.55%, respectively.
The location of study site.
In May of 2019, the field experiment was established using a randomized block design (2 m × 5 m plots separated by a buffer of 1 m). In brief, 12 plots were assigned to four treatments with 3 replications. The four treatments were CK (no fertilizer), N (72 kg N hm–2), P (36 kg P hm–2), and NP (72 kg N hm–2 + 36 kg P hm–2). The N addition level (72 kg N ha–1 year–1) sufficiently simulated N critical loads (50 kg N ha–1 year–1 reported by Zong et al., 2016). N added as ammonium nitrate (NH4NO3) and P added as calcium superphosphate [Ca(H2PO4)2] were applied in May, July, and September each year since 2019. The same amount of N and P fertilizers were applied each time. NH4NO3 and Ca(H2PO4)2 were weighed, mixed with 0.5 kg sand, and evenly distributed to each plot by hand after sunset (for higher moisture).
Plant and Soil Sampling and Analyses
In mid-August 2020, aboveground biomass (AGB) and belowground biomass (BGB) were measured at the peak of plant growth. We clipped all aboveground plants at the ground level from a 0.25 m × 0.25 m quadrat in each plot and took 3 soil cores (3.5 cm diameter × 20 cm depth) in the clipped area. The soil cores were soaked in water to separate the roots using a 0.5 mm mesh. The collected plant samples were oven-dried at 65°C to constant weight and weighed as AGB and BGB.
Soil samples were collected in the clipped area using a 3.5 cm diameter corer at depths of 0–20 cm in mid-August 2020. Fresh soil was immediately processed to determine soil moisture, and the remaining soils of each sample were divided into three parts. One part was air-dried for the determination of total nutrients, another part was stored at 4°C for the determination of available nutrients, and the remaining part was stored at −80°C for DNA extraction.
An element analyzer (EA 3000, Italy) was used to measure total C and N concentration in plant shoot (Cshoot and Nshoot), root (Croot and Nroot), and soil (TC and TN), and an ICP-AES analyzer (Thermo-Jarrell Ash Corp., MA, United States) was used to measure total P concentration in plant shoot (Pshoot), root (Proot), and soil (TP). A flow injection auto-analyzer (AACE, Germany) was used to measure soil NH4+ and NO3–, and inductively coupled plasma spectrometers (ICP; SPECTRO ARCOS EOP, Germany) were used to measure the soil’s available phosphorus (AP).
Microbial DNA Extraction and Quantitative PCR Analysis
Soil DNA was extracted from 0.5 g soil using a FastDNA Spin kit for soil (MP Biomedicals, United States) by following the manufacturer’s instructions. DNA quality was assessed using a spectrophotometer (NanoDrop ND-1000, Thermo Scientific, Waltham, MA, United States). DNA concentrations were quantified with a QuantiFluor dsDNA system (Promega, Madison, WI, United States) using a microplate reader (Spectramax M5, United States). DNA was stored at −20°C for later use.
Furthermore, 48 primer sets targeting genes were used to characterize the soil DNA, including 24 genes involved in C cycling, 16 genes in N cycling, 8 genes in P cycling, and 1 reference gene (16S rRNA) (Supplementary Table 1). High-throughput qPCR was used in this study to detect the copy numbers of these genes. Triplicate amplicons were analyzed using the Wafergen SmartChip Real-time PCR system (Wafergen, Fremont, CA, United States) as follows: 95°C for 10 min, followed by 40 cycles of 95°C for 30 s, 60°C for 30 s, and 72°C for 30 s. The melting curve was automatically generated by WaferGen software. Results with multiple melting peaks or amplification efficiencies < 90% and > 110% were discarded by the SmartChip qPCR software. A threshold cycle (CT) < 31 was used for subsequent analysis. The gene absolute copy number was calculated according to Eq. 1, and the gene relative copy number was calculated according to Eq. 2, as described by Dong S. et al. (2020).
where CT is the threshold cycle, CNi is the copy number of one of the 48 genes, and CN16S is the copy number of the 16S rRNA gene.
Statistical Analysis
A one-way ANOVA was performed to identify significant differences in microbial functional genes, soil nutrients, and plant biomass and nutrients in response to N and P addition. A two-way ANOVA was performed to test the effects of N addition, P addition, and their interactions on soil C, N, and P cycling genes. Principal coordinates analysis (PCoA) based on the Bray–Curtis distances was applied to evaluate the differences of functional genes among experimental treatments. Non-parametric multivariate statistical analyses, including multi-response permutation procedures (MRPP), analysis of similarity (ANOSIM), and permutational multivariate analysis of variance (ADONIS) were used to determine the dissimilarity of microbial functional gene abundance among different treatments. Pearson correlations were performed to analyze the relationships between microbial functional genes and plant and soil nutrients. The Mantel test was conducted to examine the correlations between microbial functional genes and soil and plant variables. The variation partitioning analysis (VPA) was performed to determine the relationships between soil and plant variables and microbial functional genes. All statistical analyses were performed using the R Development Core Team. “vegan” and “corrplot” package (version 4.0.4).
ResultsImpacts of N and P Addition on Microbial Functional Genes
Phosphorus addition had significant effects on 39 of the 48 C, N, and P cycling genes (Supplementary Table 2). In contrast to CK, P alone addition significantly (p < 0.05) increased the abundance of 24 functional genes (korA, aclB, acsA, mct, smtA, gam, glx, IsoP, mnp, xylA, nifH, amoB, hzsB, niK3, nirS1, nirS2, nosZ1, nosZ2, napA, phnK, phoD, phoX, ppx, and emGDH) but had no significant effects (p > 0.05) on the other 24 functional genes. N and NP additions had insignificant (p > 0.05) effects on all functional genes (Figures 2–4). However, there were interactive effects between N and P on 35 of the 48 C, N, and P cycling genes (Supplementary Table 2).
Abundance of C fixation (A) and C degradation (B) genes among experimental treatments. Different letters indicate significant differences at a p < 0.05 level across experimental treatments. No letter indicates a non-significant difference (CK = control, N = addition of 72 kg N hm–2 year–1, P = addition of 36 kg P hm–2 year–1, NP = interaction of N and P).
Abundance of N cycling functional genes among experimental treatments. Different letters indicate significant differences at a p < 0.05 level across experimental treatments. No letter indicates a non-significant difference (CK = control, N = addition of 72 kg N hm–2 year–1, P = addition of 36 kg P hm–2 year–1, NP = interaction of N and P).
Abundance of P cycling functional genes among experimental treatments. Different letters indicate significant differences at a p < 0.05 level across experimental treatments. No letter indicates a non-significant difference (CK = control, N = addition of 72 kg N hm–2 year–1, P = addition of 36 kg P hm–2 year–1, NP = interaction of N and P).
The PCoA plot revealed a distinct separation in β diversity by N and P addition, with the first two axes explaining 66.76% (axis 1) and 20.95% (axis 2) (Figure 5A). The cluster analysis showed that C, N, and P cycling genes of P addition were clearly separated from those under other treatments (CK, N, and NP) (Figure 5B). Non-parametric multivariate dissimilarity tests (MRPP, Adonis, and ANOSIM) further indicated that the differences in the pattern of C, N, and P cycling genes were significant (Supplementary Table 3).
(A) Principal coordinates analysis (PCoA) and (B) cluster analysis of microbial functional gene composition (CK = control, N = addition of 72 kg N hm–2 year–1, P = addition of 36 kg P hm–2 year–1, NP = interaction of N and P).
Impacts of N and P Additions on Soil and Plant Variables
Compared with CK, N, and P additions had no significant effects on soil TC, and TN, NH4+, and AP. Soil total P (TP) and NO3– increased significantly under P addition (p < 0.05) (Figure 6).
Soil nutrients under different experimental treatments. Different letters indicate significant differences at a p < 0.05 level across experimental treatments. No letter indicates a non-significant difference (CK = control, N = addition of 72 kg N hm–2 year–1, P = addition of 36 kg P hm–2 year–1, NP = interaction of N and P).
N alone addition showed no significant effects on plant biomass and nutrients (p > 0.05, Figure 7). P alone addition significantly (p < 0.05) increased Pshoot. Combined N and P addition significantly (p < 0.05) increased Nshoot.
Plant biomass and nutrients under different experimental treatments. Different letters indicate significant differences at a p < 0.05 level across experimental treatments. No letter indicates a non-significant difference (CK = control, N = addition of 72 kg N hm–2 year–1, P = addition of 36 kg P hm–2 year–1, NP = interaction of N and P).
Relationship Between Microbial Functional Genes and Environmental Variables
Most C, N, and P cycling genes were positively correlated with Pshoot and negatively correlated with Croot (p < 0.05, Supplementary Figures 2–4). A few functional genes were positively correlated with NH4+, AP, TP, and Proot and were negatively correlated with NO3– and Nshoot (p < 0.05). Nroot, Cshoot, TN, and TC were not significantly correlated with all microbial functional genes (p > 0.05). It showed that P, whether in plants or in soil, was positively correlated with microbial functional genes (Supplementary Figures 2–4).
The Mantel test showed that the functional genes were significantly correlated with NO3–, Pshoot, and Nshoot (p < 0.05). Additionally, plant and soil variables were independently tested with the C, N, and P cycling genes. The C, N, and P cycling genes were all significantly (p < 0.05) correlated with NO3– and Pshoot, and the C cycling genes were also significantly correlated with Nshoot and Croot.
The results from VPA showed that, for all functional genes, plant variables (Pshoot and Nshoot) alone explained 43% of the variance, and the interaction between plant variables, and NO3– explained 27% of the variance (Figure 8A). For C cycling genes, plant variables (Pshoot, Nshoot, and Croot) alone accounted for 37% of the variance, while the interaction between plant variables and NO3– accounted for 30% of the variance (Figure 8B). For N cycling genes, plant variables (Pshoot) and NO3– explained 37 and 2% of the variance, respectively, and the interaction between plant variables and NO3– accounted for 19% of the variance (Figure 8C). For P cycling genes, plant variables (Pshoot) and NO3– contributed to 44 and 4% of the variance, respectively, and the interaction between plant variables and NO3– contributed to 21% of the variance (Figure 8D). The Mantel test and the VPA results showed that Pshoot and NO3– were the main factors influencing the C, N, and P cycling genes.
Variation partitioning analysis (VPA) of all functional genes (A), C cycling genes (B), N cycling genes (C), and P cycling genes (D) explained by soil and plant variables. Nshoot, N concentration of plant shoot; Pshoot, P concentration of plant shoot; Croot, C concentration of plant root.
DiscussionNeutral Effects of N Addition on Microbial Functional Genes
All C, N, and P cycling genes showed a consistent neutral response to N addition, which did not support our first hypothesis and was inconsistent with previous scholars’ results (Zhang et al., 2019; Chen et al., 2020; Xiao et al., 2020). For example, N addition significantly increased the abundance of nitrifiers and denitrifiers in a meadow steppe (Xiao et al., 2020) of C degradation genes in an alpine steppe (Chen et al., 2020), whereas the neutral effects of N addition on the abundance of nifH and nirS have also been reported (Sun et al., 2015). Zhang et al. (2019) reported that N addition decreased the abundance of nifH, chiA, AOA, nirS, nirK, and nosZ. These discrepancies may be attributed to the differences in climate, soil properties, and vegetation characteristics in specific sites. In this study, we proposed two possible explanations for the results of this study.
One non-exclusive explanation refers to the N condition of this alpine meadow. The neutral effects of N addition on microbial functional genes in this study reiterate that N may be not limiting factors in this alpine ecosystem, corroborating previous findings reported from another alpine ecosystem (Jing et al., 2016). The background of N deposition (8 kg N ha–1 year–1) at this region over the past decades may have provided adequate N for the alpine plants and soil microbes (Lü and Tian, 2007). Once N saturation is reached, that is, exceeds the threshold of demand by plants and soil microbes, the additional N input no longer increases soil microbial abundance and diversity (Zhang et al., 2017). The insignificant responses of plant biomass (AGB and BGB) (Figure 7) and soil microbial biomass (Xiao et al., 2022) to N addition in the present study further proved that this alpine meadow ecosystem may be not limited by N.
Another plausible explanation may be related to the short-term study durations. The short timescale (2 years) of this study may hinder our examination of significant changes in microbial functional genes. The timescale over which to detect the response of soil microbial community to nutrient input is important (Jing et al., 2016). Allison and Martiny (2008) have reported that an average of 8 years was required before the microbial community significantly responded to fertilization. Rinnan et al. (2007) also demonstrated that, in contrast to significant responses after 15 years of fertilization to a tundra heath, 5, 6, and 10 treatment years had slight effects on microbial community composition, indicating the need for at least 10 years for the appearance of significant responses. The present study began in May 2019 and therefore the results represent the first 2 years after N addition. Several more years might be needed before a notable shift in the microbial community after the given N addition, highlighting the necessity for a long-term experiment for the reliable prediction of N deposition effects.
Positive Effect of P Addition on Soil Microbial Functional Genes
Consistent with our first hypothesis, P alone addition strikingly increased the abundance of most microbial functional genes, indicating the particular sensitivity of the soil microbial community to P input. However, inconsistent with our second hypothesis, we observed non-significant effects of the combined addition of N and P on the abundance of microbial functional genes. The PCoA and the cluster analysis also showed that P treatment clustered together and distinguished from CK, N, and NP treatment, further confirming the limitation of P to soil microbial community in this alpine ecosystem. Similar influence and limitation of P on soil microbial communities in other ecosystems have also been observed (Jing et al., 2016; Wang et al., 2018; Li et al., 2020; Zou et al., 2020). Wang et al. (2018) reported that P addition increased the abundance of soil nifH gene in a fir plantation. Li et al. (2020) observed significant correlations between the abundance of C, N, and P cycling genes and the different levels of P amendments.
In alpine ecosystems, soil P is considered a limiting nutrient for plant growth due to the low decomposition rate (Jiang et al., 2013). P can be critical to soil microorganisms when it is limited (Zou et al., 2020). Appropriate P input can alleviate plants’ P limitation, promote plant growth, and then enhance rhizosphere secretions and C input to the soil, and the increased soil C from plant subsequently provide more substrates for soil microbes (Li et al., 2017; Jing et al., 2019). In this study, we found that P addition enhanced most of the C, N, and P cycling genes, indicating the important influence of P on the biogeochemical cycle in the alpine ecosystem. Tang et al. (2016) also reported that the soil’s available P was strongly correlated with all N cycling genes except amoA, AOB, and nirS, suggesting that P plays a major role in regulating soil N cycling microbe. P limitation may occur in the growth of microorganisms because P-rich RNA needs P to synthesize proteins (He and Dijkstra, 2015). We speculate that P additions alleviate the P limitations of soil microorganisms and provide a more suitable environment for them. However, P addition showed a non-significant influence on microbial functional genes abundance when N was added, implying that N addition overwhelmed the positive effects of P on the abundance of microbial functional genes.
Relationship Between Microbial Functional Genes and Environmental Variables
The structure and the function of the soil microbial community can be shaped by soil properties and plant traits (Li et al., 2017; Wang et al., 2017). In this study, the Mantel test showed that NO3– and Pshoot significantly correlated with the whole community structure of C, N, and P cycling genes. In general, N availability can regulate the abundance, composition, and functional diversity of soil microbial microorganisms (Frey et al., 2004), because it directly restricts microbial growth (Treseder, 2010). Cong et al. (2015) demonstrated the key role of soil available N in controlling soil microbial functional gene structure and metabolic potential in a tropical rainforest. Combined with the results of VPA, Pshoot appeared to be the main factor altering the microbial functional communities in this study, because the effects of NO3– alone were weak (Figure 8). In this study, Pshoot positively correlated with many C, N, and P cycling genes (Supplementary Figures 2–4). Based on this, we speculate that plant traits have indirect effects on soil microbial function genes. Cong et al. (2015) also reported that plant communities are major players in shaping microbial community structure. Moreover, a previous study in a semi-arid grassland ecosystem demonstrated the indirect effects of plant biomass on soil microbial functional genes (Li et al., 2017), supporting the strong links between plant traits and soil microbial community. We proposed that the increase in Pshoot induced by P addition would indirectly promote microbial functional genes through increasing rhizosphere secretion and C inputs to the soil, because plant aboveground traits, as well as symbiosis and exudation, can alter soil abiotic factors that regulate the abundance, composition, and activity of soil microbial communities (Abalos et al., 2019). Not only plant traits are important for the structure and function of soil microbial populations but the plant-microorganism interaction also can influence the composition of the microbial community.
Conclusion
This study examined the responses of the C, N, and P biogeochemical cycle to N and P additions at the functional gene level in the alpine meadow. The results showed that N and P additions had inconsistent effects on the abundances of microbial functional genes. N addition had no significant effects on the abundances of C, N, and P cycling genes, while P addition significantly increased most of them, suggesting that soil microbial communities were more sensitive to P availability than to N availability in this alpine ecosystem. Combined N and P additions obviously decreased the abundances of microbial functional genes compared to P alone addition, which may be due to the overwhelming effects of N on P. The results confirmed the key role of P in C, N, and P cycling processes. P addition may increase microbial functional genes by enhancing Pshoot because Pshoot explained a large part of the variance of microbial functional genes. In future studies, the important role of plant traits in regulating the biogeochemical cycle should be addressed under N and P additions in the alpine grassland ecosystem.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
JX: conceptualization, investigation, writing-original draft, and visualization. SD: supervision, resources, and funding acquisition. HS and SL: project administration and investigation. SL and WL: project administration, and writing-review. YZ and ZM: resources and methodology. KW and HL: writing-review and editing. All authors contributed to the article and approved the submitted version.
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.
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.
Funding
This study was funded by the Second Tibetan Plateau Scientific Expedition and Research Program (2019QZKK0307), the National Science Foundation of China (U20A2007-01), the Key R&D Program of Qinghai Province (2019-SF-145 and 2018-NK-A2), and the Qinghai innovation platform construction project (2017-ZJ-Y20).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.860590/full#supplementary-material
ReferencesAbalosD.GroenigenJ.PhilippotL.LubbersI. M.DeynG. (2019). Plant trait-based approaches to improve nitrogen cycling in agroecosystems.J. Appl. Ecol.562454–2466. 10.1111/1365-2664.13489AllisonS. D.MartinyJ. (2008). Resistance, resilience, and redundancy in microbial communities.Proc. Natl. Acad. Sci. U.S.A.10511512–11519. 10.1073/pnas.080192510518695234BaeK.FaheyT. J.YanaiR. D.FiskM. (2015). Soil nitrogen availability affects belowground carbon allocation and soil respiration in northern hardwood forests of new hampshire.Ecosystems181179–1191. 10.1007/s10021-015-9892-7ChenQ. Y.YuanY. L.HuY. L.JianW.GsD.RxA. (2020). Excessive nitrogen addition accelerates N assimilation and P utilization by enhancing organic carbon decomposition in a Tibetan alpine steppe.Sci. Total Environ.764:142848. 10.1016/j.scitotenv.2020.14284833268263ChhabraS.BrazilD.MorrisseyJ.BurkeJ.O’GaraF.DowlingD. N. (2013). Fertilization management affects the alkaline phosphatase bacterial community in barley rhizosphere soil.Biol. Fertility Soils4931–39. 10.1007/s00374-012-0693-2CongJ.LiuX.LuH.XuH.LiY.DengY. (2015). Available nitrogen is the key factor influencing soil microbial functional gene diversity in tropical rainforest.BMC Microbiol.15:167. 10.1186/s12866-015-0491-8CuiY. X.ZhangY. L.DuanC. J.WangX.ZhangX. C.JuW. L. (2020). Ecoenzymatic stoichiometry reveals microbial phosphorus limitation decreases the nitrogen cycling potential of soils in semi-arid agricultural ecosystems.Soil Tillage Res.197:104463. 10.1016/j.still.2019.104463DongJ.WangS.NiuH.CuiX.WangK. (2020). Responses of soil microbes and their interactions with plant community after nitrogen and phosphorus addition in a Tibetan alpine steppe.J. Soils Sediments202236–2247. 10.1007/s11368-020-02586-3DongS. K.GaoW. H.XuC. G.HouX. Y.LongR. J.KangM. Y. (2007). Farmer and professional attitudes to the large-scale ban on livestock grazing of grasslands in China.Environ. Conserv.34246–254.DongS.LiY.GanjurjavH.GaoQ.GaoX.ZhangJ. (2020). Grazing promoted soil microbial functional genes for regulating C and N cycling in alpine meadow of the Qinghai-Tibetan Plateau.Agric. Ecosyst. Environ.303:107111. 10.1016/j.agee.2020.107111FreyS. D.KnorrM.ParrentJ. L.SimpsonR. T. (2004). Chronic nitrogen enrichment affects the structure and function of the soil microbial community in temperate hardwood and pine forests.For. Ecol. Manag.196159–171. 10.1016/j.foreco.2004.03.018HeJ. Z.ShenJ. P.ZhangL. M.ZhuY. G.ZhengY. M.XuM. G. (2010). Quantitative analyses of the abundance and composition of ammonia-oxidizing bacteria and ammonia-oxidizing archaea of a Chinese upland red soil under long-term fertilization practices.Environ. Microbiol.92364–2374. 10.1111/j.1462-2920.2007.01358.x17686032HeM.DijkstraF. A. (2015). Phosphorus addition enhances loss of nitrogen in a phosphorus-poor soil.Soil Biol. Biochem.8299–106. 10.1016/j.soilbio.2014.12.015HeZ.DengY.ZhouJ. (2012). Development of functional gene microarrays for microbial community analysis.Curr. Opin. Biotechnol.2349–55. 10.1016/j.copbio.2011.11.00122100036JiangJ.ZongN.SongM.ShiP.MaW.FuG. (2013). Responses of ecosystem respiration and its components to fertilization in an alpine meadow on the Tibetan Plateau.Eur. J. Soil Biol.56101–106. 10.1016/j.ejsobi.2013.03.001JingH.ZhangP.LiJ.YaoX.LiuG.WangG. (2019). Effect of nitrogen addition on the decomposition and release of compounds from fine roots with different diameters: the importance of initial substrate chemistry.Plant Soil151:108041.JingX.YangX.RenF.ZhouH.ZhuB.HeJ. S. (2016). Neutral effect of nitrogen addition and negative effect of phosphorus addition on topsoil extracellular enzymatic activities in an alpine grassland ecosystem.Appl. Soil Ecol.107205–213. 10.1016/j.apsoil.2016.06.004KangL.HanX.ZhangZ.SunO. J. (2007). Grassland ecosystems in China: review of current knowledge and research advancement.Philos. Trans. R. Soc. Lond.362997–1008. 10.1098/rstb.2007.202917317645KastlE. M.Schloter-HaiB.BueggerF.SchloterM. (2015). Impact of fertilization on the abundance of nitrifiers and denitrifiers at the root–soil interface of plants with different uptake strategies for nitrogen.Biol. Fertil. Soils5157–64. 10.1007/s00374-014-0948-1LauberC. L.StricklandM. S.BradfordM. A.FiererN. (2008). The influence of soil properties on the structure of bacterial and fungal communities across land-use types.Soil Biol. Biochem.402407–2415. 10.1016/j.soilbio.2008.05.021LiH.BiQ.YangK.LassonS. B.ZhengB.CuiL. (2020). High starter phosphorus fertilization facilitates soil phosphorus turnover by promoting microbial functional interaction in an arable soil.J. Environ. Sci.94179–185. 10.1016/j.jes.2020.03.04032563482LiH.YangS.XuZ.YanQ.LiX.NostrandJ. V. (2017). Responses of soil microbial functional genes to global changes are indirectly influenced by aboveground plant biomass variation.Soil Biol. Biochem.10418–29. 10.1016/j.soilbio.2016.10.009LiJ.YangC.LiuX.ShaoX. (2019). Inconsistent stoichiometry response of grasses and forbs to nitrogen and water additions in an alpine meadow of the Qinghai-Tibet plateau.Agric. Ecosyst. Environ.279178–186. 10.1016/j.agee.2018.12.016LiuL.ZhangX.LuX. (2015). The composition, seasonal variation, and potential sources of the atmospheric wet sulfur (S) and nitrogen (N) deposition in the southwest of china.Environ. Sci. Pollut. Res. Int.236363–6375. 10.1007/s11356-015-5844-126620861LüC.TianH. (2007). Spatial and temporal patterns of nitrogen deposition in China: synthesis of observational data.J. Geophys. Res.112:D22S05.LynnT. M.GeT.YuanH.WeiX.WhiteleyA. S. (2016). Soil carbon-fixation rates and associated bacterial diversity and abundance in three natural ecosystems.Microb. Ecol.731–13. 10.1007/s00248-016-0890-x27838764MaW.JiangS.AssemienF.QinM.MaB.XieZ. (2016). Response of microbial functional groups involved in soil N cycle to N, P and NP fertilization in Tibetan alpine meadows.Soil Biol. Biochem.101195–206. 10.1016/j.soilbio.2016.07.023PanY.NorikoC.MattiasD. H.MendesL. W.HeinK.GeertsR. H. E. M. (2014). Impact of long-term N, P, K, and NPK fertilization on the composition and potential functions of the bacterial community in grassland soil.Fems Microbiol. Ecol.90195–205. 10.1111/1574-6941.1238425046442PenuelasJ.PoulterB.SardansJ.CiaisP.Van Der VeldeM.BoppL. (2013). Human-induced nitrogen–phosphorus imbalances alter natural and managed ecosystems across the globe.Nat. Commun.41–10. 10.1038/ncomms393424343268RinnanR.MichelsenA.BthE.JonassonS. (2007). Fifteen years of climate change manipulations alter soil microbial communities in a subarctic heath ecosystem.Glob. Change Biol.1328–39. 10.1111/j.1365-2486.2006.01263.xSheW.BaiY.ZhangY.QinS.WeiF.SunY. (2018). Resource availability drives responses of soil microbial communities to short-term precipitation and nitrogen addition in a desert shrubland.Front. Microbiol.9:186. 10.3389/fmicb.2018.00186SmercinaD. N.BaileyV. L.HofmockelK. S. (2021). Micro on a macroscale: relating microbial-scale soil processes to global ecosystem function.FEMS Microbiol. Ecol.97:fiab091. 10.1093/femsec/fiab09134223869SunR. B.GuoX. S.WangD. Z.ChuH. Y. (2015). Effects of long-term application of chemical and organic fertilizers on the abundance of microbial communities involved in the nitrogen cycle.Appl. Soil Ecol.95171–178. 10.1139/cjm-2018-068330901528TanH.BarretM.MooijM. J.RiceO.MorrisseyJ. P.DobsonA. (2013). Long-term phosphorus fertilisation increased the diversity of the total bacterial community and the phoD phosphorus mineraliser group in pasture soils.Biol. Fertil. Soils49661–672. 10.1007/s00374-012-0755-5TangY.ZhangX.LiD.WangH.ChenF.FuX. (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.103284–293. 10.1016/j.soilbio.2016.09.001TresederK. K. (2010). Nitrogen additions and microbial biomass: a meta-analysis of ecosystem studies.Ecol. Lett.111111–1120. 10.1111/j.1461-0248.2008.01230.x18673384TurnerS.MikuttaR.GuggenbergerG.SchaarschmidtF.SchippersA. (2019). Distinct pattern of nitrogen functional gene abundances in top- and subsoils along a 120,000-year ecosystem development gradient.Soil Biol. Biochem.132111–119. 10.1016/j.soilbio.2019.02.006VitousekP. M.HowarthR. W. (1991). Nitrogen limitation on land and in the sea: how can it occur?Biogeochemistry1387–115.WangQ.LiuY. R.ZhangC. J.ZhangL. M.HanL. L.ShenJ. P. (2017). Responses of soil nitrous oxide production and abundances and composition of associated microbial communities to nitrogen and water amendment.Biol. Fertil. Soils53601–611. 10.1007/s00374-017-1203-3WangQ.WangJ.LiY.ChenD.AoJ.ZhouW. (2018). Influence of nitrogen and phosphorus additions on N2-fixation activity, abundance, and composition of diazotrophic communities in a Chinese fir plantation.Sci. Total Environ.619-6201530–1537. 10.1016/j.scitotenv.2017.10.06429129329XiaoH.YangH.ZhaoM.MonacoT. A.WangD. (2020). Soil extracellular enzyme activities and the abundance of nitrogen-cycling functional genes responded more to N addition than P addition in an Inner Mongolian meadow steppe.Sci. Total Environ.759:143541. 10.1016/j.scitotenv.2020.14354133198996XiaoJ.DongS.ShenH.LiS.DingC. (2022). Phosphorus addition promotes nitrogen retention in alpine grassland plants while increasing N deposition.Catena210:105887. 10.1016/j.catena.2021.105887ZhangJ.AiZ.LiangC.WangG.XueS. (2017). Response of soil microbial communities and nitrogen thresholds of bothriochloa ischaemum to short-term nitrogen addition on the loess plateau.Geoderma308112–119. 10.1016/j.geoderma.2017.08.034ZhangX.JohnstonE. R.WangY.YuQ.HanX. (2019). Distinct drivers of core and accessory components of soil microbial community functional diversity under environmental changes.mSystems4:e00374–19. 10.1128/mSystems.00374-1931575666ZongN.ShiP.SongM.ZhangX.JingJ.XiC. (2016). Nitrogen critical loads for an alpine meadow ecosystem on the Tibetan plateau.Environ. Manag.57531–542. 10.1007/s00267-015-0626-626475686ZouY.NingD.HuangY.LiangY.ZhouJ. (2020). Functional structures of soil microbial community relate to contrasting N2O emission patterns from a highly acidified forest.Sci. Total Environ.725:138504. 10.1016/j.scitotenv.2020.13850432302850