These authors contributed equally to this work.
Microbial stoichiometry and its potential driving factors play crucial roles in understanding the balance of chemical elements in ecological interactions and nutrient limitations along the aridity gradient. However, little is known about the variation in these features along the aridity gradient due to the lack of comprehensive field investigations. Additionally, previous studies focused on the topsoil (0–10 or 0–20 cm); however, the minimum sampling depth for topsoil could impact the results of the vertical distribution of microbial stoichiometry. In the present study, we measured the variation in microbial stoichiometry, examined the major influential factors (climatic, edaphic, and biotic factors) along an aridity gradient, and determined whether the sampling depth affected microbial C : N : P stoichiometry. From the topsoil to the subsoil, the microbial C : N, C : P, and N : P ratios varied from 6.59 to 6.83, from 60.2 to 60.5, and from 9.29 to 8.91, respectively. Only the microbial C : N ratio significantly increased with soil depth. The microbial C : N ratio significantly increased with increasing aridity in both topsoil and subsoil, while the microbial N : P ratio decreased along the aridity gradient only for the topsoil. This result implied that drought-stimulated microbes tend to be more N conservative, especially those in topsoil. Among all the factors, the soil organic carbon (SOC) content and the fungi-to-bacteria ratio exerted the largest influence on the microbial C : N, C : P, and N : P ratios at both soil depths, implying that the substrate supply and microbial structure together controlled the microbial stoichiometry. The results also revealed that the aridity index (AI) and plant aboveground biomass (AGB) exerted negative impacts on the microbial C : N ratio at both soil depths, and the effects of AI decreased in the subsoil. The results of this study suggested that the flexibility of the microbial N : P ratio should be considered when establishing the sampling depth for microbial stoichiometry study.
Ecological stoichiometry is a powerful tool for understanding the balance of chemical elements required by organisms and the functions of ecosystems (Elser et al., 2000a; Sterner and Elser, 2002). C, N, and P are regarded as critical elements in global biogeochemical cycling, and C : N : P stoichiometry in soil microorganisms offers essential insight into the nutrient limitations of the organisms and communities in an ecosystem (Manzoni et al., 2010; McGroddy et al., 2004).
A few studies have addressed the pattern of microbial stoichiometry along latitudinal (Cleveland and Liptzin, 2007; Li et al., 2015; Xu et al., 2013) or environmental gradients (Li and Chen, 2004; Li et al., 2012, 2015). For example, Cleveland and Liptzin (2007) analyzed microbial stoichiometry at the global scale and showed an increasing trend in the microbial N : P ratio with higher latitude. However, Li et al. (2015) summarized the data and found that the microbial N : P ratio decreased with latitude. Undoubtedly, there is uncertainty in the values of microbial C : N, C : P, and N : P ratios in global studies, and the variations in these patterns might have been partially caused by the different methods that were used in the various studies (Xu et al., 2013; Chen et al., 2016). Furthermore, less exploration of soil microbial stoichiometry along an aridity gradient at the regional scale impedes our ability to disentangle the trend of the changes in microbial stoichiometry amid climate changes. Climate change, such as global warming, is increasing the degree of aridity in drylands owing to the decreased precipitation and increased evaporation (Wang et al., 2014; Li et al., 2017). Given this background, key ecosystem processes that are regulated by soil microbes, such as soil respiration and nutrient mineralization, may be dramatically impacted by the increased degree of aridity, especially in fragile areas of arid and semiarid ecosystems (Delgado-Baquerizo et al., 2013; Chen et al., 2014). Therefore, we conducted a field investigation across a 2100 km climatic transect in the Inner Mongolian grasslands to determine how the microbial C : N, C : P, and N : P ratios were affected by changing environmental conditions.
Previous studies have also shown that a variety of abiotic factors impact microbial C : N : P stoichiometry (Cleveland and Liptzin, 2007; Manzoni et al., 2010; Hartman and Richardson, 2013). For instance, control experiments have found that warming can indirectly affect the turnover of microbial biomass N by stimulating soil respiration (Veraart et al., 2011; Butterbach-Bahl et al., 2013). Climate could exert an influence on microbial stoichiometry through changes to the microenvironment, such as soil moisture and temperature, and it could also impact the availability of substrates in the soil (Nielsen et al., 2009). Moreover, edaphic variables, such as soil organic carbon (SOC; Maria et al., 2014; Chen et al., 2016) and soil texture (Li et al., 2015), could be associated with nutrient mineralization and availability, thus influencing the C : N : P stoichiometry in microbial biomass (Griffiths et al., 2012). A labeled incubation experiment showed that the mineralization of organic P was mainly driven by the microbial C demands in P-poor soils (Aponte et al., 2010; Heuck et al., 2015). In addition, microbial C : N, C : P, and N : P ratios were also affected by biotic factors such as plant productivity and the composition of the microbial community (Fanin et al., 2013; Chen et al., 2016). Generally, fungi exhibit a higher C : N ratio than bacteria (Strickland and Rousk, 2010); thus, a shift in the fungi-to-bacteria ratio is expected to result in microbial stoichiometry changes (Li et al., 2012; Heuck et al., 2015). However, those findings were based on literature analyses or small-scale experiments, and the variations in microbial C : N, C : P, and N : P ratios at the regional scale have rarely been assessed systematically, and the drivers of these variations need to be addressed more specifically with appropriate experimental designs. Moreover, most research has focused on the top 10 cm of soil, which often has high C availability and nutrient contents. It can be assumed that the effects of potential driving factors exhibit minimal differentiation at deeper soil depths. However, soil at a deeper depth might contain microbial communities that are specialized for their environment, and their functions may differ from the functions of the communities in the topsoil (Fritze et al., 2000; Blume et al., 2002). Certainly, the drivers that are responsible for the variations in microbial C : N, C : P, and N : P ratios in deeper soil remain poorly understood. Such knowledge of the nature of soil microbial stoichiometry is fundamental for understanding ecosystem function, especially at the 10–20 cm soil depth, which remains highly uncertain in published studies.
Substrates for microorganisms, such as available nutrients and water, decline exponentially with depth, and the top 20 cm of soil accumulates the greatest amount of microbial biomass, thereby attracting the attention of most researchers (Fierer et al., 2003; Xu et al., 2013). Soil at a 0–20 cm depth was regarded as the topsoil in some studies, while other researchers divided the soil from 0 to 20 cm into different soil depths to explore the vertical differences between these depths (Aponte et al., 2010; Peng and Wang, 2016). However, most studies used 0–10 cm as the topsoil to facilitate sampling and comparative research (Li and Chen, 2004; Cleveland and Liptzin, 2007; Chen et al., 2016). The depth of topsoil varies among studies, and sampling depth can therefore have impacts on the study of the vertical patterns in soil microbial stoichiometry (Tischer et al., 2014). Given that soil represents a highly heterogeneous environment, especially in terms of site-specific soil development history, it is difficult to draw general conclusions (Xu et al., 2013; Camenzind et al., 2018). In addition, if a large difference existed between the soil at 0–10 cm and that at 10–20 cm, microbial stoichiometry would be underestimated due to the ambiguous limitation of topsoil (Tischer et al., 2014). To identify the soil depth that is appropriate for sampling and to improve the understanding of topsoil research at a global scale, we designed a study that divided the topsoil into 0–10 and 10–20 cm depths to compare the differences in microbial stoichiometry at the regional scale.
In Inner Mongolia grasslands, the aridity exhibits a gradient that increases from northeast to southwest (aridity index, calculated as precipitation divided by potential evapotranspiration, ranges from 0.16 to 0.54), thus providing an ideal platform to better estimate the patterns and drivers of microbial C : N : P stoichiometry along an aridity gradient (Chen et al., 2014; Li et al., 2017). In this study, we aim to access the effect of soil depth on soil microbial C : N : P stoichiometry along the aridity gradient. We hypothesized that the microbial C : N and C : P ratios decrease and the microbial N : P ratio increases with temperature (Cleveland and Liptzin, 2007; Li et al., 2015; Xu et al., 2013), and the microbial C : N and C : P ratios decrease and the microbial N : P ratio increases with decreasing aridity index (Wang et al., 2014; Li et al., 2017). In addition, the identification of soil depth for vertical study is different in some published literature (Li and Chen, 2004; Aponte et al., 2010; Tischer et al., 2014; Peng and Wang, 2016). We predicted that variation in bacterial and fungal taxa between soil depths might contribute to the shifts in C : N : P stoichiometry, especially in the N : P ratio (Mouginot et al., 2014; Camenzind et al., 2018). Therefore, we focus on (i) the effects of potential driving factors on microbial C : N, C : P, and N : P ratios in topsoil and subsoil and (ii) the response of the microbial C : N, C : P, and N : P ratios to soil depth.
This study was performed across the Inner Mongolian temperate grassland,
which is a central part of the Eurasian steppe. The study area is located at
39.2–49.6
Along this transect, a total of 58 sites that were slightly disturbed by
humans and domestic animals were sampled, including 10 in the meadow steppe,
28 in the typical steppe, and 20 in the desert steppe (Fig. S1 in the Supplement). Five
1 m
The aridity index was extracted from the Global Aridity Index
(Global-Aridity) dataset, which provides high-resolution (30 arcsec or
Microbial biomass carbon (MBC) and microbial biomass nitrogen (MBN) were
determined following the chloroform fumigation–
The C : N, C : P, and N : P ratios in the soil microbial biomass were
The results indicate well-constrained relationships among C, N, and P in soil
microbial biomass (Fig. S2). The soil microbial biomass C : N, C : P, and N : P
ratios varied by an order of magnitude. Significantly different water
content, soil bulk density, sand percentages, and SOC content were found
between soil depths (
Basic characteristics of study sites across the Inner Mongolian
grassland at different soil depths. Different letters indicate significant
differences between soil depths on
Basic information of study sites.
Note that data represent the means, with minimum and maximum values in parentheses. MAT, mean annual temperature; MAP, mean annual precipitation; AGB, aboveground biomass.
The microbial biomass C, N, and P concentrations and microbial C : N : P stoichiometric ratios across the Inner Mongolian grassland at different soil depths.
Note that different letters indicate significant differences between soil depths
based on log10-transformed data (paired
The results revealed a significant positive relationship between the aridity index (AI) and
the microbial C : N ratio (topsoil,
Relationships between the C : N, C : P, and N : P ratios in soil
microbial biomass and aridity index
Relationships between the C : N, C : P, and N : P ratios in the soil
microbial biomass and AGB
The final SEM adequately fit the data, as shown by several robust
goodness-of-fit measures (
The structural equation model (SEM) shows the direct and indirect
influences of various ecological factors on the microbial C : N
As stated in our hypothesis, the increase in the microbial C : N ratio and the decrease in the microbial N : P ratio that were found along a temperature gradient in this study are in accordance with the findings of Li et al. (2015) and Chen et al. (2016), who reported similar variations in microbial stoichiometry along latitudinal gradient. Temperature drives the variation in the growth of the microbial community, as high growth rates at low latitudes require high RNA contents, causing the N : P ratio to decline (Chadwick et al., 1999; Kooijman et al., 2009; Xu et al., 2013). In addition, we observed that the microbial C : N ratio significantly increased with increasing aridity index, while the microbial N : P ratio decreased with increasing aridity index, indicating that drought (decreasing aridity index) affects ecological stoichiometry by mediating the growth rate of microorganisms in semiarid regions (Elser et al., 2000b; Peng and Wang, 2016). Wang et al. (2014) and Kou et al. (2018) both noted that increasing aridity reduced the soil microbial abundance in drylands, and a decreased growth rate in dry areas might result in decreased allocation to P-rich ribosomal RNA (and thus higher C : P and N : P ratios) (Maestre et al., 2015). Additionally, microbial C : N ratio decreased with the decreasing aridity index, which serves as a protective mechanism as microbes decrease their nitrogen use efficiency (NUE, the ratio of N invested in growth over total N uptake) and tend to be more N conservative under drier climatic conditions (Mooshammer et al., 2014; Delgado-Baquerizo et al., 2017). Moreover, under drier climate conditions, the soil microbial communities shift from acting as r-strategists (fast-growing copiotrophs) to acting as K-strategists (slow-growing oligotrophs), as microorganisms with K-strategies have lower nutrient demands (N and P) and growth rates, invest more nutrients into extracellular enzymes to gain limited nutrients and thus have higher cellular C : N : P ratios than r-strategists (Fierer et al., 2007, 2010).
The microbial C : N ratio demonstrated an increasing trend with increasing latitude, in contrast to the decreasing trend that was demonstrated for the microbial N : P ratio. Such results paralleled the results of studies on ecological stoichiometry, which revealed that the C : N ratio of microorganisms increased with latitude, while the N : P ratio decreased with latitude, suggesting increasing N limitations in microorganism ecosystems in high-latitude areas (Li et al., 2015; Chen et al., 2016). The regional-scale microbial stoichiometry followed the global-scale stoichiometry patterns that were observed for plant leaves (Reich and Oleksyn, 2004; Yuan et al., 2011), litter (McGroddy et al., 2004), and soil (Sardans et al., 2012), conforming to the substrate age hypothesis, which predicts young soils to be N-limited, whereas old soils tend to be P-limited (Walker and Syers, 1976; Vitousek et al., 2010). Our study further illustrated the latitudinal pattern of microbial stoichiometry and first attempted to examine the variation in microbial stoichiometry along an aridity gradient at the regional scale.
Among the ecological factors examined, our study found that the patterns of microbial C : N and C : P ratios were associated with SOC content and the F : B ratio, suggesting that the available C and microbial community structure together regulated the shift in microbial stoichiometry. If the environmental parameters were considered individually, SOC content was found to be significantly positively related to the microbial C : N and C : P ratios, which is consistent with the results observed from global data analysis, suggesting that SOC content may control microbial stoichiometry by mediating the substrate stoichiometry, e.g., the soil C : N and C : P ratios (Hartman and Richardson, 2013; Maria et al., 2014; Mooshammer et al., 2014). In deeper soil, microbial metabolic processes are limited by C availability and energy (C), such as denitrification and P mineralization, thus resulting in the effect of SOC on microbial C : N : P stoichiometry in subsoil (Fierer et al., 2003; Peng and Wang, 2016; Camenzind et al., 2018).
SEM also illustrated that the microbial community structure is an important feature in determining microbial stoichiometry. The F : B ratio has recently been found to have a vital influence on the patterns of microbial C : N and N : P ratios in soil at a large scale (Chen et al., 2016). An experiment indicated that fungi have lower resource requirements and higher C : N and C : P ratios than bacteria; thus, the shift in the F : B ratio impacted microbial C : N : P stoichiometry (Mouginot et al., 2014). In our study, the lower F : B ratio might have led to a shift in the microbial nutrient stoichiometry at deeper soil depths (Tischer et al., 2014). Overall, the SEM highlighted the important role of the C availability and microbial community structure in driving the variations in microbial C : N, C : P, and N : P ratios at both soil depths.
Moreover, AGB and AI also exerted direct influences on microbial C : N or C : P ratios, and those impacts mainly acted in the topsoil but were weaker in the subsoil. The climate imposes important controls on both the plant community and the microbial taxa along with their interactions with soil nutrients (Chadwick et al., 1999; Vitousek et al., 2010; Reich and Oleksyn, 2004). In particular, drier weather conditions, and the decreasing aridity index, could affect the growth and productivity of plants and then shape a shift in vegetation types along this grassland transect (Jaleel et al., 2009; Cherwin and Knapp, 2012). At the same time, the meadow steppe ecosystem with high productivity maintained relatively high soil C and N contents, which resulted in high C : P and N : P ratios in these regions; thus, plant productivity exerted a positive influence on microbial C : N (Aponte et al., 2010; Manzoni et al., 2010). Because of the vertical distribution of the influence of the AI, the effect decreased with soil depth.
Interestingly, our results revealed that the microbial N : P ratio was mainly impacted by the F : B ratio and SOC content, while the sand percentage and SOC content had direct negative effects on the ratio in the subsoil, suggesting the flexibility of microbial stoichiometry in response to distinct resource supplies between topsoil and subsoil (Peng and Wang, 2016). The soil depth affected the microbial biomass N and P, which decreased nearly 2-fold from the topsoil to the subsoil (Table 2). However, the results showed that the N and P cycles responded asymmetrically to soil depth, which might be attributed to the high variability in P availability (Li et al., 2015; Zechmeister-Boltenstern et al., 2016). Generally, P is mostly derived from parent material, while N is mainly a biological element (Vitousek and Farrington, 1997; Vitousek et al., 2010). Therefore, it is believed that P variations regulate large-scale patterns in microbial N : P stoichiometry and nutrient use strategies (Heuck et al., 2015; Camenzind et al., 2018). With a high proportion of sand, the soil becomes porous, which may lead to increased leaching of available P to deeper soil depths (Otten et al., 1999; Achbergerová and Nahálka, 2011). Similarly, P leaching caused by weathering led to a shift in the N : P ratio in the soil, and a vertical study found a high variation in the N : P ratio between soil depths across a large scale (Tian et al., 2010). The high variability of the N : P ratio in soil and soil microbial biomass therefore indicates that the N : P ratio could be an indicator of the ecosystem nutrient status at deeper soil depths (Cleveland and Liptzin, 2007; Li et al., 2015).
The results showed significant differences in the water content and SOC content between topsoil and subsoil, suggesting that the resource supplies between topsoil and subsoil were significantly different. We also observed that the microbial C : N, C : P, and N : P ratios varied between topsoil and subsoil, and a significant difference was found in the microbial C : N ratio. Those results indicated that the flexibility of the microbial stoichiometry responds to different resource supplies between soil depths (Tian et al., 2010; Peng and Wang, 2016). Similar findings were found in the top 16 cm of soil in a Mediterranean oak forest (0–8 and 8–16 cm), where the microbial nutrient ratios (C : N, C : P, and N : P) varied with soil depths (Aponte et al., 2010). Tischer et al. (2014) sampled the top 20 cm of soil (0–5, 5–10, 10–20 cm) and observed that the microbial C : N ratio changed with soil depth. Moreover, sampling to a depth of 10 cm showed a significant difference in the microbial N : P ratio (Tischer et al., 2014). The detection of the differences in the microbial N : P ratio in our study depended strongly on the sampling depth, suggesting that the microbial N : P ratio might provide insight into the nature of ecosystem nutrient limitations in a vertical study (Cleveland and Liptzin, 2007; Fierer et al., 2010). In addition, SEM also showed that the microbial N : P ratio was controlled by multiple potential driving factors at different soil depths, indicating that a 0–10 cm or shallower sampling interval should be used when studying the vertical patterns in the microbial N : P ratio.
The first uncertainty was related to the determination of fungal and bacterial biomasses by PLFA markers, which have limited targets for fungi and bacteria. This uncertainty should be noted when interpreting the results in the present study. Methodological advances in sequencing approaches might be used to more accurately index the microbial community and reveal insights into the regulation of microbial C : N : P stoichiometry in distinct soil microbial taxa or functional groups. Second, the theory used to construct the model is another source of uncertainty; the theory related to the drivers of microbial stoichiometry used in this study was mostly derived from a literature review and summarized data. In future research, more control experiments with the manipulation of C availability, especially at deeper soil depths, would further improve our understanding of the changes in microbial stoichiometry and nutrient limitations under the impacts of global change.
The ratios of C, N, and P in the microbial biomass were 6.59 : 60.2 : 9.29 in the topsoil, which deviated from the 6.83 : 60.5 : 8.91 ratio in the subsoil. Moreover, significant differences were found in the microbial C : N and N : P ratios between topsoil and subsoil, indicating that the flexibility of microbial stoichiometry should be considered for vertical study. In addition, the trend of microbial C : N ratio increasing with aridity index is consistent with the perspective that microbes mediate their nitrogen use efficiency and tend to be more N conservative under drier climatic conditions. The microbial N : P ratio trend along the aridity gradient was consistent with the growth rate hypothesis that a decreased growth rate in dry areas results in decreased allocation to P-rich ribosomal RNA and thus a higher N : P ratio. These findings confirmed the importance of SOC content, the microbial structure, and soil texture in shaping the pattern of microbial stoichiometry in semiarid grassland systems. The influence of ecological factors decreased from topsoil to subsoil, and the climatic, edaphic, and biotic factors declined as well. Overall, these results illustrated N and P limitation in microbial biomass at deeper soil depths along the aridity gradient and limited responses to ecological factors in the subsoil.
Data for aridity index are available on the global aridity and PET database at
The supplement related to this article is available online at:
HH, WM, and YY devised the study. YL carried out the experiment and data analyses. DK, YC, and DC assisted with the data analyses and interpretation. XN, TW, XZ, MZ, and HB assisted with the experiment. All authors contributed to the preparation of the paper.
The authors declare that they have no conflict of interest.
We thank the sampling team members of Institute of Botany, Chinese Academy of Science for assistance with the field data collection.
This research has been supported by the National Basic Research Program of China (grant nos. 2015CB954201 and 2014CB954303) and the National Natural Science Foundation of China (grant no. 31470498).
This paper was edited by Yakov Kuzyakov and reviewed by three anonymous referees.