- 1State Key Laboratory of Urban and Regional Ecology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences (CAS), Beijing, China
- 2College of Resources and Environment, University of Chinese Academy of Sciences (CAS), Beijing, China
- 3School of Geography, Nanjing Normal University, Nanjing, China
- 4State Key Laboratory of Vegetation and Environmental Change, Institute of Botany, Chinese Academy of Sciences (CAS), Beijing, China
- 5Jiangsu Center for Collaborative Innovation in Geographical Information Resource Development and Application, Nanjing, China
- 6Zhongke Clean Soil (Guangzhou) Technology Service Co., Ltd., Guangzhou, China
- 7Suzhou Station of Farmland Quality Protection, Suzhou, China
- 8Key Laboratory of Virtual Geographical Environment, Ministry of Education, Nanjing Normal University, Nanjing, China
The microbial groups of nitrogen fixers, ammonia oxidizers, and denitrifiers play vital roles in driving the nitrogen cycle in grassland ecosystems. However, the understanding of the abundance and distribution of these functional microorganisms as well as their driving factors were limited mainly to topsoil. In this study, the abundances of nitrogen functional genes (NFGs) involved in nitrogen fixation (nifH), ammonia oxidation (amoA), and denitrification (nirK, nirS, and nosZ) were investigated in both topsoil (0–10 cm, soil layer with concentrated root) and subsoil (30–40 cm, soil layer with spare root) of three grassland habitats in northern China. The abundance of NFGs decreased with soil depth except for the archaeal amoA gene and the distribution of nifH, archaeal amoA, nirK, and nirS gene was significantly impacted by grassland habitats. Moreover, the distribution of NFGs was more responsive to the vertical difference than horizontal spatial heterogeneity. Redundancy analysis revealed that the distribution pattern of overall NFGs was regulated by grassland habitats, and these regulations were more obvious in the subsoil than in the topsoil. Variance partitioning analysis further indicated that soil resource supply (e.g., organic matter) may control the vertical distribution of NFGs. Taken together, the findings in this study could fundamentally improve our understanding of the distribution of N cycling-associated microorganisms across a vertical scale, which would be useful for predicting the soil N availability and guiding the soil N management in grassland ecosystems.
Introduction
Nitrogen (N) cycling is considered a vital biogeochemical cycle on Earth as it controls the availability of nitrogenous nutrients and biological productivity in natural ecosystems (Vitousek and Howarth, 1991; Nelson et al., 2016). The underlying microbe-mediated functions are associated with the availability of nitrogenous nutrient supply and, thus, influence plant-based primary production ecosystems (Wang H. et al., 2017; Song et al., 2020; Wang et al., 2021; Zhou et al., 2021). A better understanding of the distribution of N functional microbes remains constrained.
Soil N cycle serves as an important component of the global N cycle, and the transformation among various N forms is largely accomplished by the enzymes encoded by several N functional genes (NFGs) harbored in diverse soil microbes (Robertson and Groffman, 2007; Hu et al., 2015). For instance, the nitrogen fixation (nifH) gene encodes the nitrogenase and catalyzes the conversion of N2 to NH4+, which, thereby, provides a primary source of N to natural ecosystems (Levy-Booth et al., 2014). Nitrification is a two-step process linking the reduced and oxidized forms of N, in which ammonia oxidation (amoA) is thought to be a rate-limiting step and mediated by the ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) with the monooxygenase encoded by the gene amoA (Prosser and Nicol, 2008; Levy-Booth et al., 2014). Increasing evidence suggested that the NFGs abundances were correlated with the rate of the N cycling processes to some extent, and some of their ratios were considered as a useful indicator for the potential of N turnover capacity (Prosser and Nicol, 2012; Carey et al., 2016). Based on path analysis, AOA/AOB and nirK/S and nosZ gene abundances were the most important variables to predict nitrification and denitrification potential rates across terrestrial ecosystems (Xiong et al., 2012). Regional-scale experiments offered direct evidence of the relationship between N cycling process rate and N functional gene abundance (Chen et al., 2018). Moreover, the relative abundance of NFGs, such as AOA/AOB and (AOA + AOB)/(nirK + nirS) quantified as N turnover potential capacity, represented dominance in ammonia oxidization and NO3-N leaching, respectively (Tang et al., 2016, 2018). Therefore, better understanding the abundance and distribution of NFGs will offer an insight to know the soil N availability for a given natural ecosystem and thus to guide for predicting plant primary productivity.
Over the past few decades, spatial patterns and drivers of soil microbial communities have attracted extensive interest at various spatial scales (Fierer et al., 2010; Serna-Chavez et al., 2013; Fierer, 2017). These studies demonstrated that microbial taxa inhabiting varied surface soil habitats exhibit biogeographic patterns across broad geographical space (Chu et al., 2010; Brockett et al., 2012). It is believed that environmental factors, such as altitude, longitude, and latitude, usually covary with the changes in the geographical distance (Fierer and Jackson, 2006; Griffiths et al., 2011; Jia et al., 2017). Such horizontal differences are often characterized by heterogeneity of soil type, climate, vegetation, etc., and, of course, pronouncedly shape the distribution of microorganisms (Zhang et al., 2006; Hu et al., 2014; Chen et al., 2018). Moreover, the geographical distribution of microbial populations was also driven by geographical isolation even with similar characteristics and common ancestry (O’Malley, 2008). A previous study found that there is habitat specificity for microbial diversity and is no considerable dispersal limitation in bacteria distribution at a large scale in China (Wang X. B. et al., 2017). Although these findings have drastically improved our understanding of the pattern of microbes, we are still not sure if there is a horizontal difference in N-cycling microbes over hundreds of kilometers. Such knowledge is needed to improve our understanding of the N transformation processes in the grassland ecosystem especially under the ground of global climate change.
It is well recognized that both the biotic and abiotic factors could impact the abundance and distribution of N-cycling microbial communities (Hayden et al., 2010; Stone et al., 2015; Chen et al., 2018). For instance, Wei et al. (2017) and Cheng et al. (2019) exhibited that plant traits such as plant type, species richness, and biomass were important predictors of abundance and activity of ammonia-oxidizing and denitrifying communities in the grassland ecosystem. In contrast, numerous studies indicated that abiotic factors such as precipitation, soil pH, and organic carbon were pivotal to the proliferation of N-cycling microorganisms (Chen et al., 2015; Barrett et al., 2016; Zhou et al., 2021). Moreover, abiotic factors, especially precipitation, could indirectly affect the soil N-cycling microbes by impacting the accumulation of grassland biomass (Prober et al., 2015; Wang X. B. et al., 2017). However, these reports are restricted to the N-cycling microbes residing in near-surface soil while there was still less knowledge about their counterpart microbes inhabiting in deeper soil depths. Especially, Wakelin et al. (2008) provided evidence where soil C is strongly linked to diazotrophic ecology, pointing out that the quality and quantity of soil C influence the N-fixing bacteria. In addition, Chen et al. (2018) exhibited that abiotic and plant factors greatly predict the abundance and activity of ammonia-oxidizing and denitrifying communities across the Qinghai-Tibetan plateau.
Soil depth is regarded as an important factor in determining the distribution of microbial communities because of the environmental gradients along with the soil vertical profiles (Bai et al., 2017). Particularly, the oxygen level, moisture, nutrient availability, etc., sharply decrease with soil depth, and the soil layer of >30 cm is commonly characterized as the oligotrophic environment (Whitman et al., 1998; Eilers et al., 2012). Especially, more than 80% of roots distribute in the top 30 cm of soil in the grassland ecosystem (Yang et al., 2010). In fact, root vertical distribution could impact the microbial populations between soil depths mainly by providing different quantity and type of organic carbons and nutrients within a soil profile via root exudates and/or root litter decomposition (Peng et al., 2017; Shi et al., 2018). Therefore, the microbial biomass and activity considerably declined with the increasing depths as a consequence, which might lead to a reduction in population size as well as changes in the composition and diversity of microbial communities (Babujia et al., 2010; Bai et al., 2015). For example, researchers observed that both AOA and AOB rapidly decreased between 0 and 30 cm depth and then tended to mildly reduce below 30 cm soil depth (Banning et al., 2015; Tang et al., 2018). It has been widely accepted that AOA is preferred to adaptation in poor nutrient conditions and low soil pH environments (Prosser and Nicol, 2012), making the AOA less responsive to soil depth (Wang et al., 2014). In addition, different taxa might respond dissimilarly to the changes in abiotic factors along with soil profile (Tang et al., 2018). Until present, we still lack a comprehensive understanding of the N-turnover functional gene in the grassland ecosystem, which impedes our ability to predict the key services provided by the grassland ecosystem.
Grasslands are the dominant landscapes in China, which cover 41.7% of the national territory and are mainly distributed in the northern temperate zone, ranging from temperate grasslands to alpine grasslands (Kang et al., 2007). To control the precipitation variable, sampling sites were designed along the 400 mm isohyet in the grasslands of northern China to evaluate the distribution difference and N-transformation potential in the grasslands of Northern China (Wang et al., 2018). In this study, we use quantitative polymerase chain reaction of samples taken at two soil depths, namely, topsoil at 0–10 cm depth and subsoil at 30–40 cm depth from the representative steppe, and compare the difference between soil depths over hundreds of kilometers. In this study, we use the ratio of nifH/nosZ and (AOA + AOB)/(nirK + nirS) as the proxy of potential microbial N storage and NO3–-N leaching. We hypothesized that the abundance of NFGs varied with geographical distance and soil depth. It is expected that vertical difference and horizontal geographical differences in environmental factors would shape the distinct N-cycling microbial community. We also hypothesized that the major driving factors shifted from topsoil to subsoil, and we expected that soil resources prevailed over soil edaphic parameters in driving the variation of N-cycling microbes residing in the subsoil.
Materials and Methods
Site Description and Soil Sample Collection
Four study sites with three representative grassland habitat types were located in the transition zone along 400 mm isohyet from the meadow steppe, typical steppe to the alpine meadow (Figure 1). Two sites were located on the Inner Mongolian Plateau, with meadow steppe dominated by Leymus chinensis in Eerguna [called EEGN at 520 m average sea level (a.s.l.)] and the other was typical steppe dominated by Stipa grandis in Xilinhot (called XLHT at 1250 m a.s.l.). The other two sites, namely Menyuan (MY, 3200 m a.s.l.) and Naqu (NQ, 4611 m a.s.l.), were located on the Qinghai-Tibet Plateau with the alpine meadow dominated by Kobresia humilis and Potentilla nivea. The aboveground biomass in the EEGN, XLHT, MY, and NQ sites was 112, 103, 350, and 116 g m–2, respectively; belowground biomass was 711, 2512, 3000, and 1509 g m–2, respectively (Table 1). According to the root distribution characteristics of plants, we collected soil samples from a depth of 0–10 cm (topsoil) and 30–40 cm (subsoil) at each study site in the plant growing season in September 2016. At each study site, five 1 m × 1 m subplots were established, one at each corner and one in the center of a 10 m × 10 m plot. The soil samples collected from the same layer at each plot were mixed into single composite samples, which were passed through 2 mm sieves and removed root materials. In total, 24 sieved soils (4 study sites × 3 replicates × 2 soil depths) were processed and stored as appropriate for physicochemical, biological, and molecular analyses.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792002/fpls-12-792002-HTML-r1/image_m/fpls-12-792002-g001.jpg)
Figure 1. Locations of four selected grassland sites along 400 mm isohyet from the inner Mongolian Plateau and Qinghai-Tibet Plateau in China. The site abbreviations are defined in Table 1.
Soil Physicochemical Analyses
Soil pH was measured using a glass electrode meter (S200 K, Mettler-Toledo International Inc., Shanghai, China) with a soil/water ratio of 1:2.5. Soil organic matter (OM) content was determined by H2SO4-K2Cr2O7 oxidation, while the total N (TN) content was analyzed by semimicro Kjeldahl digestion. Dissolved organic carbon (DOC) was measured in a 4:1 water/soil extract using an Analyzer Multi N/C (Analytic Jena, Jena, Germany). Soil NH4+-N and NO3–-N were extracted using 2 M KCl with a soil/solution ratio of 1:5 and measured using a continuous-flow analyzer (Skalar, Breda, Netherlands). The moisture content was determined gravimetrically by oven drying at 105°C to constant mass.
DNA Extraction and Real-Time PCR Assay
Soil genomic DNA was extracted, quality-controlled, and quantified as described previously by Zhao et al. (2018). After extraction, the DNA samples were stored at −20°C for subsequent molecular analysis. Real-time PCR was performed on a CFX-96 thermocycler (Bio-Rad Laboratories Inc., CA, United States) to enumerate the abundances of microbial NFGs that involved in nitrogen fixation (nifH), ammonia oxidation (archaeal amoA and bacterial amoA), and denitrification (nirK, nirS, and nosZ) using the primer pairs listed in Table 2. The reaction mixture contained 10 μl of Premix Ex Taq™ (Takara Bio Inc., Kyoto, Japan), 2 μl of DNA extract (20 ng μl–1), 1 μl of each primer (10 μM), and 6 μl of sterile distilled water. The PCR procedure consisted of 2 min at 95°C, followed by 39 cycles of 10 s at 95°C, 20 s at 58°C, and 30 s at 72°C. Standard curves were plotted using 10-fold serial dilutions of the plasmid DNA inserted with the fragment of interest, following a previously established procedure (Zhao et al., 2016), and the amplification efficiencies ranged from 88.3 to 97.0% for different genes. The melting curve analysis was performed to ensure the amplification specificity, and the copy numbers were log10-transformed to normalize the values prior to statistical analysis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792002/fpls-12-792002-HTML-r1/image_m/fpls-12-792002-t002.jpg)
Table 2. Summary of the genes investigated, the enzymes they encoded, and their functions in the nitrogen cycle.
Statistical Analyses
The PASW Statistics 18 (SPSS Inc., Chicago, United States) and R (version 3.3.1) (R Core Team, 2012) were used to perform statistical analysis, and P-value of <0.05 was considered statistically significant. Data were tested for normality (Shapiro–Wilk normality test) and homogeneity of variance and were logarithm transformed when necessary to meet the criteria for a normal distribution. Multiple analysis of variance (MANOVA) was used to determine the effects of geographical location and soil depth on soil properties and NFGs abundance. For a given soil depth, the significant difference among sampling sites was tested using one-way ANOVA followed by LSD post hoc test. For a given site, two-tailed unpaired t-tests were used to compare each measurement between soil depths. The abundance of overall NFGs and soil properties was visualized with redundancy analysis (RDA). The Monte Carlo test (999 permutations) was further used to assess the effects of edaphic factors on the distribution of NFGs. Variance partitioning analysis was also performed to determine the relative contributions of soil edaphic parameters (soil pH and moisture) and soil sources (OM, DOC, TN, and NH4+-N content) on the variation in NFGs distribution in both soil depths.
Results
Soil Physicochemical Properties
Overall, physicochemical properties of the topsoil or subsoil differed significantly (P < 0.05) among the selected grassland sites except for the soil NO3–-N content in the subsoil. However, the different physicochemical properties varied among these four sampling sites at both topsoil and subsoil (Table 3). For instance, soil OM, TN, and NH4+-N contents were highest in the topsoil of the MY site, followed by the NQ, EEGN, and XLHT sites, while the highest NO3–-N content was observed in the EEGN site, followed by the NQ, XLHT, and MY sites. In the subsoil, soil OM and TN content were highest in the MY site, followed by the EEGN, XLHT, and NQ sites, whereas the highest NH4+-N and NO3–-N contents were observed in the EEGN site, with the lowest value being in the XLHT and MY sites, respectively (Table 3). Moreover, most tested soil physicochemical properties showed a decreasing trend from topsoil to the subsoil, whereas soil pH increased (P < 0.05) with the soil depth, except for the EEGN site (Table 3). The results of MANOVA showed that all soil edaphic properties were affected (P < 0.05) by geographical location, soil depth, and the geographical location × soil depth interaction term except for the soil DOC content (Table 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792002/fpls-12-792002-HTML-r1/image_m/fpls-12-792002-t003.jpg)
Table 3. Physicochemical characteristics of soils from different natural grassland sites and soil depths.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792002/fpls-12-792002-HTML-r1/image_m/fpls-12-792002-t004.jpg)
Table 4. Overall effects of the sampling site and soil depth on the soil physicochemical characteristics analyzed by two-way ANOVA.
Abundance of Nitrogen Fixers and Ammonia Oxidizers
The abundance of nitrogen fixers decreased with the increase in soil depth at all the grassland sites except for the MY site, ranging from 1.29 × 107 to 3.01 × 108 copies⋅g–1 dry soil, and significant differences between topsoil and subsoil were observed in the EEGN and NQ sites (Figure 2A). Moreover, the population of nitrogen fixers differed significantly among the sampling sites in the topsoil, being higher at the EEGN and XLHT sites as compared with the MY site, while no significant difference was observed in the subsoil. The distribution of nitrogen fixers was affected (P < 0.05) by the sampling site and soil depth (Figure 2A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792002/fpls-12-792002-HTML-r1/image_m/fpls-12-792002-g002.jpg)
Figure 2. Abundances of nifH (A), archaeal amoA (B), bacterial amoA (C), nirK (D), nirS (E), and nosZ (F) genes at soil depths of 0–10 and 30–40 cm at four grassland sites. Different uppercase and lowercase letters above the boxes indicate significant differences at P < 0.05 between sites according to LSD post hoc test in the topsoil and subsoil, respectively. An asterisk above the boxes indicates a significantly different level at P < 0.05 according to the student t-test for a given site. An asterisk next to the number indicates a significantly different level at P < 0.05 according to multiple analyses of variance. SS, sampling site; SD, soil depth; SS × SD, sampling site × soil depth. *P < 0.05; **P < 0.01; ***P < 0.001. The site abbreviations are defined in Table 1.
Unlike the results of nitrogen fixers, the population size of AOA increased with the increasing soil depth in all sampling sites except for the EEGN site, ranging from 6.16 × 108 to 6.27 × 109 copies⋅g–1 dry soil, although there was no significant difference between soil depth for a given site (Figure 2B). Furthermore, significant differences were observed among sampling sites in both topsoil and subsoil. For instance, the abundance of AOA was higher in the EEGN and XLHT sites than that in the MY and NQ sites in the topsoil. In contrast, the abundance of AOB showed a declining trend from topsoil to subsoil at all the grassland sites, ranging from 2.52 × 105 to 6.29 × 106 copies⋅g–1 dry soil (Figure 2C). A significant difference between soil depth was observed only in the XLHT site. The abundance of AOB changed considerably (P < 0.05) among the sampling sites in the topsoil, whereas there was no significant difference in the subsoil. In addition, MANOVA results revealed that the distribution of AOA was affected by the sampling site while the distribution of AOB was impacted by soil depth (Figures 2B,C).
Abundance of Nitrite Reductase (nirS and nirK) and Nitrous Oxide Reductase (nosZ) Genes
The abundance of nirK- and nirS-type nitrite reductase gene shifted from 3.49 × 106 to 3.85 × 107 and from 1.01 × 106 to 2.45 × 107 copies⋅g–1 dry soil, respectively, in all soil samples. They consistently decreased from topsoil to subsoil in all the sampling sites except for nirS-type nitrite reductase gene at the MY site. These declines were significant (P < 0.05) in the EEGN and MY sites for nirK-type nitrite reductase gene and the EEGN and NQ sites for nirS-type nitrite reductase gene, respectively (Figures 2D,E). Significant differences in abundance were observed in the subsoil of nirK-type nitrite reductase gene and the topsoil of nirS-type nitrite reductase gene among the sampling sites. Both sampling sites and soil depth significantly affected the distributions of nirK- and nirS-type nitrite reductase gene (Figures 2D,E). The abundance of nitrous oxide reductase (nosZ) gene ranged from 1.81 × 107 to 2.03 × 108 copies g–1 dry soil in all soil samples and accounted for the largest proportion (43.9–88.4%) of denitrifies. It declined considerably (P < 0.05) with the increasing soil depth in all the grassland sites, whereas there were no significant differences among the sites for both soil layers (Figure 2F). Thus, only soil depth was the significant (P < 0.001) factor influencing the distribution of nitrous oxide reductase (nosZ) gene.
Potential Nitrogen-Cycling Capacity
The topsoil N-storage potential at the EEGN and XLHT sites was significantly (P < 0.05) higher than that at the MY and NQ sites. The N-storage potential increased with the increase in soil depth at all the grassland sites, and the significant differences between topsoil and subsoil were observed at the MY and NQ sites (Figure 3A). The NO3–-N leaching potential in the EEGN and XLHT topsoil was almost two times higher than that at the MY and NQ sites, and the NO3–-N leaching potential significantly (P < 0.05) decreased with soil depth at all the sites except for MY (Figure 3B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792002/fpls-12-792002-HTML-r1/image_m/fpls-12-792002-g003.jpg)
Figure 3. Microbial N storage potential (A) and nitrate leaching potential (B) at soil depths of 0–10 and 30–40 cm at four grassland sites. Values (means ± SD, n = 3) within the same soil depth followed by different letters are significantly different at P < 0.05 according to LSD post hoc test. An asterisk above the bars indicates a significantly different level at P < 0.05 according to the student t-test for a given sampling site. *P < 0.05; **P < 0.01. The site abbreviations are defined in Table 1.
Relationships Between Overall Nitrogen Functional Genes and Soil Edaphic Factors
The Monte Carlo permutation test showed a significant (F = 2.71, P = 0.004) and positive correlation between the ordination of overall NFGs and soil edaphic properties. The RDA biplot revealed that the ordination of the overall NFGs was separated from the soil layers in the first axis, and the pattern of ordination in the topsoil was considerably distinct from that in the subsoil (Figure 4 and Supplementary Figure 1). In topsoil, the ordinations of overall NFGs at the MY and NQ sites were similar, with clear separation from that at the EEGN sites and a small overlap in the XLHT site (Figure 4A). The differences in the ordination of overall NFGs were enlarged in the subsoil. Especially, the ordinations of NFGs at the MY and NQ sites were grouped together, and they were distinctly separated from that at the EEGN and XLHT sites, with the EEGN and XLHT sites being separated from each other (Figure 4B). The effect of individual soil factors on the ordination of NFGs is shown by the direction and the length of the vectors. Soil moisture, OM, DOC, TN, and NH4+-N contents were closely (P < 0.05) correlated with the ordination of the overall NFGs in topsoil, whereas no significant correlations were observed between the NFGs ordination and individual soil property in the subsoil (Figure 4 and Supplementary Figure 1). Variance partitioning analysis showed the relative contributions of soil physical parameters and soil resource supply on NFGs distribution (Figure 4). The soil edaphic parameters and soil resource variables explained 23.9 and 34.1% of the observed variation in NFGs distribution in the topsoil, while these variables explained 9.7 and 26.2% of the observed variation, respectively (Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792002/fpls-12-792002-HTML-r1/image_m/fpls-12-792002-g004.jpg)
Figure 4. Redundancy analysis of overall nitrogen functional genes and soil characteristics for individual samples at a soil depth of 0–10 cm (A) and 30–40 cm (B), respectively. The site abbreviations are defined in Table 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792002/fpls-12-792002-HTML-r1/image_m/fpls-12-792002-g005.jpg)
Figure 5. Variance partitioning analysis (VPA) of the effects of soil physical parameter and soil resource supply on the variation in the NFGs distribution at a soil depth of 0–10 cm (A) and 30–40 cm (B). Circles on the edges of the triangle show the percentage of variation explained by each variable.
Discussion
Changes in the Distribution of Nitrogen-Cycling Microbes Across Three Grassland Habitats
We assumed that the effect of the geographical distance on N-cycling microorganisms was embodied in the effects of temperature and/or precipitation (Griffiths et al., 2011; Brockett et al., 2012; Yang et al., 2014). As a result, our experiment was designed along with 400 mm isohyet and thus eliminated differences in the environmental factors relating to geographical location (Jia et al., 2017). In this study, geographical location considerably impacted the distribution of N-fixer, AOA, nirK-type, and nirS-type denitrifiers, suggesting that a geographical distance over hundreds of kilometers might influence the distribution pattern of N-cycling microorganisms (De Vries et al., 2012; Chu et al., 2016). This result might be attributed to heterogeneous climate, vegetation, and soil type caused by geographical distances between every pair of sites ranging from 762 to 3069 km (Supplementary Table 1). As per the principle of microbial biogeography, everything is everywhere, but the environment selects (de Wit and Bouvier, 2006), distribution of free-living microbial taxa exhibited dissimilarity with increasing geographical distance (Chu et al., 2016; Azziz et al., 2017; Escalas et al., 2021). The AOA gene abundance was significantly influenced by the geomorphic zones at the landscape scale (Hayden et al., 2010; Song et al., 2020; Zhou et al., 2021), which was consistent with this study. Previous studies also found that soil type exerted an effect on the abundance of the nirK-type and nirS-type denitrifies from two geographical zones (Azziz et al., 2017). Notably, the abundance of AOB and nosZ genes was less responsive to geographical location (Figure 2), suggesting that dispersal limitation might slightly influence the distribution of AOB and nosZ-type denitrifier in the study area (Chu et al., 2016; Wang X. B. et al., 2017). The ratio of nifH/nosZ and (AOA + AOB)/(nirK + nirS) in the EEGN and XLHT topsoil is significantly higher than that in the MY and NQ topsoil, suggesting that the temperate meadow might have a higher potential N fixation and NO3–-N leaching than alpine meadow (Tang et al., 2018). This result might be explained by the fact that microbes related energy-consume N process habitat preference and, thus, microbes in an alpine meadow with lower temperature may limit N fixation and NO3–-N leaching capacity than that in temperature meadow (Chen et al., 2018; Frindte et al., 2019).
Shifts in the Abundance of Nitrogen Functional Genes Between Soil Depths
In general, subsoil differed in a range of edaphic factors from topsoils (Fierer et al., 2003; Rumpel and Kögel-Knabner, 2011; Schnecker et al., 2015), such as soil pH, moisture, oxygen level, and nutrient availability and consequently led to a shift in the abundance of soil microorganisms (Eilers et al., 2012; Wang et al., 2014; Ying et al., 2017). Alternatively, the distribution of AOA was not affected by the soil depth, thus indicating that AOA preferred to adaptation in lower soil depth with poor nutrient condition (Prosser and Nicol, 2012; Sun et al., 2019; Wang et al., 2019). It is also likely that the root system of vegetation cannot reach a certain depth of soil, supporting a depth-related root disturbance on the soil microbial community (Haichar et al., 2008; Shi et al., 2018). The relatively low N-cycling gene in the alpine meadow might be attributed to the fact that alpine grasslands have a relatively much shallower root distribution than temperate grassland (Yang et al., 2010). These observations together indicate that the N processes and soil functionality decline along with soil profiles, which is consistent with the results of previous studies (Bai et al., 2017; Jones et al., 2018; Tang et al., 2018), demonstrating that the soil functionality was greatly affected by soil depths (Yao et al., 2018; Xu et al., 2020). Researchers have found that nirK and nirS-type denitrifiers preferred to reside in surface soil with high organic C (Bru et al., 2011; Petersen et al., 2012; Prosser and Nicol, 2012). Soil C source and C/N ratio decreased from topsoil with dense root to subsoil with spare root, resulting in the reduction of nirK and nirS-type denitrifiers, and thus decrease in NO3–-N leaching potential between soil depths (Fierer et al., 2003; Tang et al., 2018). However, an increase in N-storage potential with increasing soil depths was observed in the MY and NQ sites, which might be attributed to permafrost-affected deeper alpine meadow soil (Chen et al., 2018).
Soil Resource Supply Drives the Distribution of Nitrogen Functional Genes
As shown by RDA, the abundances in NFGs at the NQ and MY sites were grouped together, and they were distinctly separated from the EEGN and MY sites, respectively, in the subsoil, suggesting their distinct N-cycling microbial pattern in the subsoil from kilometers apart. These results are supported by previous findings that the abundance and distribution of dominant taxonomic groups varied across habitats in the Chinese northern dryland (Wang X. B. et al., 2017). These findings together indicate that the distribution of N-cycling microbes is structured by habitat-associated vegetation type in the subsoil. In this study, the overall distribution of NFGs depicted a habitat-specific pattern, reflecting their ecological characteristics such as physiological capabilities or habitat preferences (Fierer et al., 2007; Pointing et al., 2009). We found that variation in the distribution of NFGs in the subsoil was primarily related to soil resource supply, but not to the soil edaphic parameters. The clear and broad consensus is that soil microbes are limited by soil resources including C and N availability (Fierer et al., 2003; Hu et al., 2014; Chen et al., 2016). Recently, Chen et al. (2018) examined the geographical pattern of the functional N gene using intensive soil sampling across the Tibetan alpine grassland and exhibited that soil resource was the main predictor of N-cycling microbial community, while plant traits affected them via altering the SOC availability (Li et al., 2017). Likewise, we observed that soil resource was an important predictor in the topsoil. Moreover, the soil resource explained more variation in N-cycling microbial community in deeper soil depth, indicating that the distribution of N-cycling microbes residing in the deeper soil depths was associated with soil resource supply (Yang et al., 2013; Liu et al., 2020). This result was supported by the finding that the vertical distribution of N-cycling microorganisms was largely attributed to the decline in soil C availability. Barrett et al. (2016) reported that the C amendment resulted in the increase of denitrifies abundance and this impact changed with soil depth due to the decline in C availability. In addition, the interaction of plant and edaphic factors in temperate grasslands affects the microbial community distribution and the effect depends on the soil depth (Yao et al., 2018). Hou et al. (2018) also presented direct evidence that root exudates can act as inducible C sources for heterotrophic denitrifying bacteria. Thus, the plant–soil interaction (i.e., rhizosphere effect and litter decomposition) involved in the distribution of NFGs should be further studied to reveal how the plant traits regulate the abundance of NFGs.
Conclusion
In general, the abundance of NFGs gradually decreased with soil depth in the three grassland habitats. We also found that the distribution of free N-cycling microbial taxa exhibited dissimilarity at the landscape scale. In addition, NFGs distribution was more responsive to the vertical difference than horizontal spatial heterogeneity except for the archaeal amoA gene. The variation in N-transformation potential among the three grassland habitats may be connected with the energy-consuming N-cycling microorganism habitat preference. The distribution of overall NFGs exerted a habitat-specific pattern and the driving factors of NFGs shifted from the topsoil to the subsoil. Moreover, soil resource supply controlled the vertical distribution of NFGs rather than soil edaphic parameters.
Uncertainties and Perspectives
The first uncertainty was the sampling rule along 400 mm isohyet, which neglects various environmental factors across the large geographical distance. It is also difficult to completely disentangle the vertical N-cycling ecological processes by only determining the NFGs abundance in limited sampling sites. Intensive sampling would be helpful to improve the understanding of active functional groups within the surface and subsoils and determine how these functional groups change with the soil profile. Second, we assume the copy number of a gene provides information about the potential N store and N release capacity of soil microorganisms in grassland, which may overestimate the nitrogen turnover potential. Besides, our research takes no account of other N-cycling processes, such as the mineralization of fungi. We do not provide a complete insight into measuring the capacity of microbial N turnover. Molecular analyses can be integrated with gene expression and rhizosphere isotopic tracing techniques to measure actual N storage and turnover potentials in future research.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author Contributions
JBZ, XF, JZ, and TW devised the study. YL carried out the experiments, analyzed the data, and wrote the first draft of the manuscript. QC, HC, and HD assisted with the experiments. All authors contributed to the preparation of the manuscript.
Funding
This work was financially supported by the National Natural Science Foundation of China (41977084 and 41977032), the Key-Area Research and Development Program of Guangdong Province (2020B0202010006), and the Opening Project of Zhejiang Key Laboratory of Urban Environmental Processes and Pollution Control (NUEORS202002).
Conflict of Interest
JZ, JBZ, and ZC were employed by the company Zhongke Clean Soil (Guangzhou) Technology Service Co., Ltd.
The remaining 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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.792002/full#supplementary-material
References
Azziz, G., Monza, J., Etchebehere, C., and Irisarri, P. (2017). nirS- and nirK-type denitrifier communities are differentially affected by soil type, rice cultivar and water management. Eur. J. Soil Biol. 78, 20–28. doi: 10.1016/j.ejsobi.2016.11.003
Babujia, L. C., Hungria, M., Franchini, J. C., and Brookes, P. C. (2010). Microbial biomass and activity at various soil depths in a Brazilian oxisol after two decades of no-tillage and conventional tillage. Soil Biol. Biochem. 42, 2174–2181.
Bai, R., Wang, J.-T., Deng, Y., He, J.-Z., Feng, K., and Zhang, L.-M. (2017). Microbial community and functional structure significantly varied among distinct types of paddy soils but responded differently along gradients of soil depth layers. Front. Microbiol. 8:945. doi: 10.3389/fmicb.2017.00945
Bai, R., Xi, D., He, J.-Z., Hu, H.-W., Fang, Y.-T., and Zhang, L.-M. (2015). Activity, abundance and community structure of anammox bacteria along depth profiles in three different paddy soils. Soil Biol. Biochem. 91, 212–221. doi: 10.1016/j.soilbio.2015.08.040
Banning, N. C., Maccarone, L. D., Fisk, L. M., and Murphy, D. V. (2015). Ammonia-oxidising bacteria not archaea dominate nitrification activity in semi-arid agricultural soil. Sci. Rep. 5:11146. doi: 10.1038/srep11146
Barrett, M., Khalil, M. I., Jahangir, M. M. R., Lee, C., Cardenas, L. M., Collins, G., et al. (2016). Carbon amendment and soil depth affect the distribution and abundance of denitrifiers in agricultural soils. Environ. Sci. Pollut. Res. 23, 7899–7910. doi: 10.1007/s11356-015-6030-1
Braker, G., Fesefeldt, A., and Witzel, K. P. (1998). Development of PCR primer systems for amplification of nitrite reductase genes (nirK and nirS) to detect denitrifying bacteria in environmental samples. Appl. Environ. Microbiol. 64, 3769–3775.
Brockett, B. F. T., Prescott, C. E., and Grayston, S. J. (2012). Soil moisture is the major factor influencing microbial community structure and enzyme activities across seven biogeoclimatic zones in western Canada. Soil Biol. Biochem. 44, 9–20. doi: 10.1016/j.soilbio.2011.09.003
Bru, D., Ramette, A., Saby, N. P., Dequiedt, S., Ranjard, L., Jolivet, C., et al. (2011). Determinants of the distribution of nitrogen-cycling microbial communities at the landscape scale. ISME J. 5, 532–542. doi: 10.1038/ismej.2010.130
Carey, C. J., Dove, N. C., Beman, J. M., Hart, S. C., and Aronson, E. L. (2016). Meta-analysis reveals ammonia-oxidizing bacteria respond more strongly to nitrogen addition than ammonia-oxidizing archaea. Soil Biol. Biochem. 99, 158–166. doi: 10.1016/j.soilbio.2016.05.014
Chen, D., Cheng, J., Chu, P., Hu, S., Xie, Y., Tuvshintogtokh, I., et al. (2015). Regional-scale patterns of soil microbes and nematodes across grasslands on the Mongolian plateau: relationships with climate, soil, and plants. Ecography 38, 622–631. doi: 10.1111/ecog.01226
Chen, Y., Kou, D., Li, F., Ding, J., Yang, G., Fang, K., et al. (2018). Linkage of plant and abiotic properties to the abundance and activity of N-cycling microbial communities in Tibetan permafrost-affected regions. Plant Soil 434, 453–466. doi: 10.1007/s11104-018-3845-9
Chen, Y. L., Ding, J. Z., Peng, Y. F., Li, F., Yang, G. B., Liu, L., et al. (2016). Patterns and drivers of soil microbial communities in Tibetan alpine and global terrestrial ecosystems. J. Biogeogr. 43, 2027–2039.
Cheng, Y., Wang, J., Chang, S. X., Cai, Z., Müller, C., and Zhang, J. (2019). Nitrogen deposition affects both net and gross soil nitrogen transformations in forest ecosystems: a review. Environ. Pollut. 244, 608–616. doi: 10.1016/j.envpol.2018.10.054
Chu, H., Fierer, N., Lauber, C. L., Caporaso, J. G., Knight, R., and Grogan, P. (2010). Soil bacterial diversity in the Arctic is not fundamentally different from that found in other biomes. Environ. Microbiol. 12, 2998–3006. doi: 10.1111/j.1462-2920.2010.02277.x
Chu, H., Sun, H., Tripathi, B. M., Adams, J. M., Huang, R., Zhang, Y., et al. (2016). Bacterial community dissimilarity between the surface and subsurface soils equals horizontal differences over several kilometers in the western Tibetan Plateau. Environ. Microbiol. 18, 1523–1533. doi: 10.1111/1462-2920.13236
De Vries, F. T., Manning, P., Tallowin, J. R. B., Mortimer, S. R., Pilgrim, E. S., Harrison, K. A., et al. (2012). Abiotic drivers and plant traits explain landscape-scale patterns in soil microbial communities. Ecol. Lett. 15, 1230–1239.
de Wit, R., and Bouvier, T. (2006). ‘Everything is everywhere, but, the environment selects’; what did Baas Becking and Beijerinck really say? Environ. Microbiol. 8, 755–758. doi: 10.1111/j.1462-2920.2006.01017.x
Eilers, K. G., Debenport, S., Anderson, S., and Fierer, N. (2012). Digging deeper to find unique microbial communities: the strong effect of depth on the structure of bacterial and archaeal communities in soil. Soil Biol. Biochem. 50, 58–65. doi: 10.1016/j.soilbio.2012.03.011
Escalas, A., Paula, F. S., Guilhaumon, F., Yuan, M., Yang, Y., Wu, L., et al. (2021). Macroecological distributions of gene variants highlight the functional organization of soil microbial systems. ISME J. doi: 10.1038/s41396-021-01120-8
Fierer, N. (2017). Embracing the unknown: disentangling the complexities of the soil microbiome. Nat. Rev. Microbiol. 15, 579–590. doi: 10.1038/nrmicro.2017.87
Fierer, N., Bradford, M. A., and Jackson, R. B. (2007). Toward an ecological classification of soil bacteria. Ecology 88, 1354–1364. doi: 10.1890/05-1839
Fierer, N., and Jackson, R. B. (2006). The diversity and biogeography of soil bacterial communities. Proc. Natl. Acad. Sci. U.S.A. 103:626. doi: 10.1073/pnas.0507535103
Fierer, N., Schimel, J. P., and Holden, P. A. (2003). Variations in microbial community composition through two soil depth profiles. Soil Biol. Biochem. 35, 167–176.
Fierer, N., Strickland, M. S., Liptzin, D., Bradford, M. A., and Cleveland, C. C. (2010). Global patterns in belowground communities. Ecol. Lett. 12, 1238–1249. doi: 10.1111/j.1461-0248.2009.01360.x
Francis, C. A., Roberts, K. J., Beman, J. M., Santoro, A. E., and Oakley, B. B. (2005). Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean. Proc. Natl. Acad. Sci. U.S.A. 102:14683. doi: 10.1073/pnas.0506625102
Frindte, K., Pape, R., Werner, K., Loffler, J., and Knief, C. (2019). Temperature and soil moisture control microbial community composition in an arctic-alpine ecosystem along elevational and micro-topographic gradients. ISME J. 13, 2031–2043. doi: 10.1038/s41396-019-0409-9
Griffiths, R. I., Thomson, B. C., James, P., Bell, T., Bailey, M. J., and Whiteley, A. S. (2011). The bacterial biogeography of British soils. Environ. Microbiol. 13, 1642–1654.
Haichar, F. Z., Marol, C., Berge, O., Rangel-Castro, J. I., Prosser, J. I., Balesdent, J., et al. (2008). Plant host habitat and root exudates shape soil bacterial community structure. ISME J. 2, 1221–1230. doi: 10.1038/ismej.2008.80
Hayden, H. L., Drake, J., Imhof, M., Oxley, A. P. A., Norng, S., and Mele, P. M. (2010). The abundance of nitrogen cycle genes amoA and nifH depends on land-uses and soil types in South-Eastern Australia. Soil Biol. Biochem. 42, 1774–1783. doi: 10.1016/j.soilbio.2010.06.015
Hou, S., Ai, C., Zhou, W., Liang, G., and He, P. (2018). Structure and assembly cues for rhizospheric nirK- and nirS-type denitrifier communities in long-term fertilized soils. Soil Biol. Biochem. 119, 32–40. doi: 10.1016/j.soilbio.2018.01.007
Hu, H. W., Chen, D., and He, J. Z. (2015). Microbial regulation of terrestrial nitrous oxide formation: understanding the biological pathways for prediction of emission rates. FEMS Microbiol. Rev. 39, 729–749. doi: 10.1093/femsre/fuv021
Hu, Y., Xiang, D., Veresoglou, S. D., Chen, F., Chen, Y., Hao, Z., et al. (2014). Soil organic carbon and soil structure are driving microbial abundance and community composition across the arid and semi-arid grasslands in northern China. Soil Biol. Biochem. 77, 51–57. doi: 10.1016/j.soilbio.2014.06.014
Jia, J., Feng, X., He, J.-S., He, H., Lin, L., and Liu, Z. (2017). Comparing microbial carbon sequestration and priming in the subsoil versus topsoil of a Qinghai-Tibetan alpine grassland. Soil Biol. Biochem. 104, 141–151. doi: 10.1016/j.soilbio.2016.10.018
Jones, D. L., Magthab, E. A., Gleeson, D. B., Hill, P. W., Sánchez-Rodríguez, A. R., Roberts, P., et al. (2018). Microbial competition for nitrogen and carbon is as intense in the subsoil as in the topsoil. Soil Biol. Biochem. 117, 72–82. doi: 10.1016/j.soilbio.2017.10.024
Kang, L., Han, X., Zhang, Z., and Sun, O. J. (2007). Grassland ecosystems in China: review of current knowledge and research advancement. Philos. Trans. R. Soc. B Biol. Sci. 362, 997–1008. doi: 10.1098/rstb.2007.2029
Levy-Booth, D. J., Prescott, C. E., and Grayston, S. J. (2014). Microbial functional genes involved in nitrogen fixation, nitrification and denitrification in forest ecosystems. Soil Biol. Biochem. 75, 11–25. doi: 10.1016/j.soilbio.2014.03.021
Li, H., Yang, S., Xu, Z., Yan, Q., Li, X. 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
Liu, Y., Ma, W., Kou, D., Niu, X., Wang, T., Chen, Y., et al. (2020). A comparison of patterns of microbial C:N:P stoichiometry between topsoil and subsoil along an aridity gradient. Biogeosciences 17, 2009–2019. doi: 10.5194/bg-17-2009-2020
Nelson, M. B., Martiny, A. C., and Martiny, J. B. H. (2016). Global biogeography of microbial nitrogen-cycling traits in soil. Proc. Natl. Acad. Sci. U.S.A. 113, 8033–8040. doi: 10.1073/pnas.1601070113
O’Malley, M. A. (2008). ‘Everything is everywhere: but the environment selects’: ubiquitous distribution and ecological determinism in microbial biogeography. Stud. Hist. Philos. Biol. Biomed. Sci. 39, 314–325. doi: 10.1016/j.shpsc.2008.06.005
Peng, Y., Guo, D., and Yang, Y. (2017). Global patterns of root dynamics under nitrogen enrichment. Glob. Ecol. Biogeogr. 26, 102–114. doi: 10.1111/geb.12508
Petersen, D. G., Blazewicz, S. J., Firestone, M., Herman, D. J., Turetsky, M., and Waldrop, M. (2012). Abundance of microbial genes associated with nitrogen cycling as indices of biogeochemical process rates across a vegetation gradient in Alaska. Environ. Microbiol. 14, 993–1008. doi: 10.1111/j.1462-2920.2011.02679.x
Pointing, S. B., Chan, Y., Lacap, D. C., Lau, M. C., Jurgens, J. A., and Farrell, R. L. (2009). Highly specialized microbial diversity in hyper-arid polar desert. Proc. Natl. Acad. Sci. U.S.A. 106, 19964–19969. doi: 10.1073/pnas.0908274106
Poly, F., Ranjard, L., Nazaret, S., Gourbiere, F., and Monrozier, L. J. (2001). Comparison of nifH gene pools in soils and soil microenvironments with contrasting properties. Appl. Environ. Microbiol. 67, 2255–2262. doi: 10.1128/AEM.67.5.2255-2262.2001
Prober, S. M., Leff, J. W., Bates, S. T., Borer, E. T., Firn, J., Harpole, W. S., et al. (2015). Plant diversity predicts beta but not alpha diversity of soil microbes across grasslands worldwide. Ecol. Lett. 18, 85–95. doi: 10.1111/ele.12381
Prosser, J. I., and Nicol, G. W. (2008). Relative contributions of archaea and bacteria to aerobic ammonia oxidation in the environment. Environ. Microbiol. 10, 2931–2941. doi: 10.1111/j.1462-2920.2008.01775.x
Prosser, J. I., and Nicol, G. W. (2012). Archaeal and bacterial ammonia-oxidisers in soil: the quest for niche specialisation and differentiation. Trends Microbiol. 20, 523–531. doi: 10.1016/j.tim.2012.08.001
R Core Team (2012). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Robertson, G. P., and Groffman, P. M. (2007). “13 - Nitrogen transformations,” in Soil Microbiology, Ecology and Biochemistry, Third Edn, ed. E. A. Paul (San Diego: Academic Press), 341–364.
Rosch, C., and Bothe, H. (2005). Improved assessment of denitrifying, N2-fixing, and total-community bacteria by terminal restriction fragment length polymorphism analysis using multiple restriction enzymes. Appl. Environ. Microbiol. 71, 2026–2035. doi: 10.1128/AEM.71.4.2026-2035.2005
Rumpel, C., and Kögel-Knabner, I. (2011). Deep soil organic matter—a key but poorly understood component of terrestrial C cycle. Plant Soil 338, 143–158. doi: 10.1007/s11104-010-0391-5
Schnecker, J., Wild, B., Takriti, M., Eloy Alves, R. J., Gentsch, N., Gittel, A., et al. (2015). Microbial community composition shapes enzyme patterns in topsoil and subsoil horizons along a latitudinal transect in Western Siberia. Soil Biol. Biochem. 83, 106–115. doi: 10.1016/j.soilbio.2015.01.016
Serna-Chavez, H. M., Fierer, N., and van Bodegom, P. M. (2013). Global drivers and patterns of microbial abundance in soil. Global Ecol. Biogeogr. 22, 1162–1172. doi: 10.1111/geb.12070
Shi, S., Herman, D. J., He, Z. Pett-Ridge, J., Wu, L., Zhou, J., et al. (2018). Plant roots alter microbial functional genes supporting root litter decomposition. Soil Biol. Biochem. 127, 90–99. doi: 10.1016/j.soilbio.2018.09.013
Song, Y., Liu, C., Wang, X., Ma, X., Jiang, L., Zhu, J., et al. (2020). Microbial abundance as an indicator of soil carbon and nitrogen nutrient in permafrost peatlands. Ecol. Indic. 115:106362. doi: 10.1016/j.ecolind.2020.106362
Stephen, J. R., Kowalchuk, G. A., Bruns, M.-A. V., McCaig, A. E., Phillips, C. J., Embley, T. M., et al. (1998). Analysis of β-subgroup proteobacterial ammonia oxidizer populations in soil by denaturing gradient gel electrophoresis analysis and hierarchical phylogenetic probing. Appl. Environ. Microbiol. 64, 2958–2965. doi: 10.1128/AEM.64.8.2958-2965.1998
Stone, M. M., Kan, J., and Plante, A. F. (2015). Parent material and vegetation influence bacterial community structure and nitrogen functional genes along deep tropical soil profiles at the Luquillo Critical Zone Observatory. Soil Biol. Biochem. 80, 273–282. doi: 10.1016/j.soilbio.2014.10.019
Sun, R., Myrold, D. D., Wang, D., Guo, X., and Chu, H. (2019). AOA and AOB communities respond differently to changes of soil pH under long-term fertilization. Soil Ecol. Lett. 1, 126–135. doi: 10.1007/s42832-019-0016-8
Tang, Y., Yu, G., Zhang, X., Wang, Q., Ge, J., 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
Throbäck, I. N., and Enwall, K. Jarvis, Å, and Hallin, S. (2004). Reassessing PCR primers targeting nirS, nirK and nosZ genes for community surveys of denitrifying bacteria with DGGE. FEMS Microbiol. Ecol. 49, 401–417. doi: 10.1016/j.femsec.2004.04.011
Vitousek, P. M., and Howarth, R. W. (1991). Nitrogen limitation on land and in the sea: how can it occur? Biogeochemistry 13, 87–115. doi: 10.1007/BF00002772
Wakelin, S. A., Macdonald, L. M., Rogers, S. L., Gregg, A. L., Bolger, T. P., and Baldock, J. A. (2008). Habitat selective factors influencing the structural composition and functional capacity of microbial communities in agricultural soils. Soil Biol. Biochem. 40, 803–813. doi: 10.1016/j.soilbio.2007.10.015
Wang, D., Rui, Y., Ding, K., Cui, X., Hao, Y., Tang, L., et al. (2018). Precipitation drives the biogeographic distribution of soil fungal community in Inner Mongolian temperate grasslands. J. Soils Sediments 18, 222–228. doi: 10.1007/s11368-017-1727-z
Wang, H., Bu, L., Song, F., Tian, J., and Wei, G. (2021). Soil available nitrogen and phosphorus affected by functional bacterial community composition and diversity as ecological restoration progressed. Land Degrad. Dev. 32, 183–198. doi: 10.1002/ldr.3707
Wang, H., Li, X., Li, X., Li, X., Wang, J., and Zhang, H. (2017). Changes of microbial population and N-cycling function genes with depth in three Chinese paddy soils. PLoS One 12:e0189506. doi: 10.1371/journal.pone.0189506
Wang, J., Zhang, L., Lu, Q., Raza, W., Huang, Q., and Shen, Q. (2014). Ammonia oxidizer abundance in paddy soil profile with different fertilizer regimes. Appl. Soil Ecol. 84, 38–44. doi: 10.1016/j.apsoil.2014.06.009
Wang, X. B., Lu, X. T., Yao, J., Wang, Z. W., Deng, Y., Cheng, W. X., et al. (2017). Habitat-specific patterns and drivers of bacterial beta-diversity in China’s drylands. ISME J. 11, 1345–1358. doi: 10.1038/ismej.2017.11
Wang, Z., Cao, Y., Wright, A. L., Shi, X., and Jiang, X. (2019). Different ammonia oxidizers are responsible for nitrification in two neutral paddy soils. Soil Tillage Res. 195:104433. doi: 10.1016/j.still.2019.104433
Wei, X., Reich, P. B., Hobbie, S. E., and Kazanski, C. E. (2017). Disentangling species and functional group richness effects on soil N cycling in a grassland ecosystem. Glob. Change Biol. 23, 4717–4727. doi: 10.1111/gcb.13757
Whitman, W. B., Coleman, D. C., and Wiebe, W. J. (1998). Prokaryotes: the unseen majority. Proc. Natl. Acad. Sci. U.S.A. 95, 6578–6583. doi: 10.1073/pnas.95.12.6578
Xiong, J., Liu, Y., Lin, X., Zhang, H., Zeng, J., Hou, J., et al. (2012). Geographic distance and pH drive bacterial distribution in alkaline lake sediments across Tibetan plateau. Environ. Microbiol. 14, 2457–2466. doi: 10.1111/j.1462-2920.2012.02799.x
Xu, T., Chen, X., Hou, Y., and Zhu, B. (2020). Changes in microbial biomass, community composition and diversity, and functioning with soil depth in two alpine ecosystems on the Tibetan plateau. Plant Soil 459, 137–153. doi: 10.1007/s11104-020-04712-z
Yang, Y., Fang, J., Ma, W., Guo, D., and Mohammat, A. (2010). Large-scale pattern of biomass partitioning across China’s grasslands. Glob. Ecol. Biogeogr. 19, 268–277. doi: 10.1111/j.1466-8238.2009.00502.x
Yang, Y., Gao, Y., Wang, S., Xu, D., Yu, H., Wu, L., et al. (2014). The microbial gene diversity along an elevation gradient of the Tibetan grassland. ISME J. 8, 430–440. doi: 10.1038/ismej.2013.146
Yang, Y., Wu, L., Lin, Q., Yuan, M., Xu, D., Yu, H., et al. (2013). Responses of the functional structure of soil microbial community to livestock grazing in the Tibetan alpine grassland. Glob. Change Biol. 19, 637–648. doi: 10.1111/gcb.12065
Yao, X., Zhang, N., Zeng, H., and Wang, W. (2018). Effects of soil depth and plant–soil interaction on microbial community in temperate grasslands of northern China. Sci. Total Environ. 630, 96–102. doi: 10.1016/j.scitotenv.2018.02.155
Ying, J., Li, X., Wang, N., Lan, Z., He, J., and Bai, Y. (2017). Contrasting effects of nitrogen forms and soil pH on ammonia oxidizing microorganisms and their responses to long-term nitrogen fertilization in a typical steppe ecosystem. Soil Biol. Biochem. 107, 10–18. doi: 10.1016/j.soilbio.2016.12.023
Zhang, Y., Li, D., Wang, H., Xiao, Q., and Liu, X. (2006). Molecular diversity of nitrogen-fixing bacteria from the Tibetan Plateau, China. FEMS Microbiol. Lett. 260, 134–142. doi: 10.1111/j.1574-6968.2006.00317.x
Zhao, J., Ni, T., Li, J., Lu, Q., Fang, Z., Huang, Q., et al. (2016). Effects of organic–inorganic compound fertilizer with reduced chemical fertilizer application on crop yields, soil biological activity and bacterial community structure in a rice–wheat cropping system. Appl. Soil Ecol. 99, 1–12. doi: 10.1016/j.apsoil.2015.11.006
Zhao, J., Zhou, X., Jiang, A., Fan, J., Lan, T., Zhang, J., et al. (2018). Distinct impacts of reductive soil disinfestation and chemical soil disinfestation on soil fungal communities and memberships. Appl. Microbiol. Biotechnol. 102, 7623–7634. doi: 10.1007/s00253-018-9107-1
Keywords: grassland ecosystem, nitrogen functional gene, soil depth, real-time PCR, abundance
Citation: Liu Y, Chi Q, Cheng H, Ding H, Wen T, Zhao J, Feng X, Zhang J, Cai Z and Liu G (2022) Comparative Microbial Nitrogen Functional Gene Abundances in the Topsoil vs. Subsoil of Three Grassland Habitats in Northern China. Front. Plant Sci. 12:792002. doi: 10.3389/fpls.2021.792002
Received: 09 October 2021; Accepted: 09 December 2021;
Published: 14 January 2022.
Edited by:
Shiliang Liu, Beijing Normal University, ChinaReviewed by:
Roberta Pastorelli, Council for Agricultural and Economics Research (CREA), ItalyDima Chen, China Three Gorges University, China
Jiapeng Wu, Guangzhou University, China
Copyright © 2022 Liu, Chi, Cheng, Ding, Wen, Zhao, Feng, Zhang, Cai and Liu. 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: Teng Wen, wenteng@njnu.edu.cn; Jun Zhao, junzhao37@njnu.edu.cn