2.1 Study site description
The present study was conducted in the core areas of Hui-river National Nature Reserve (latitude: 48° 10′ to 48° 57′ N, longitude: 118° 40′ to 119° 45′ E) (Fig. 1), Hulunbuir steppe, China. This site is influenced by the temperate continental monsoon climate and has the annual average precipitation and temperature of 375 mm and –2.4 ℃ to 2.2 ℃, respectively. The soils of Hui-river riparian zones are Fluvisols and Kastanozems (WRB 2015). In 2007, this river was included in the China Biosphere Reserve Network (CBRN), and become one of the 151 critical areas for biodiversity protection. Hui-river also provides the important breeding sites and habitats for migratory birds from East Asia‒Australia (https://www.iucnredlist.org). The riparian zone of Hui-river is less disturbed bceuase its location in the core areas of the Hui-river National Nature Reserve, and has high species diversity. The dominant plant species belong to Cyperaceae, Gramineae, and Asteraceae family (Li et al. 2019).
2.2 Experimental design and sampling process
Due to the obvious variations in submergence frequency in the transversal direction (i.e., perpendicular to the river channel), the soil properties and vegetation types sharply change in this direction. For example, the areas of riparian near the river are frequently submerged, and conducive to hygrophyte plant growth (Fig. 1). However, the areas of riparian far away the river are rarely submerged or unsubmerged, and dominanted by xerophyte. According to this natural vegetation type differences in the transverse direction of the riparian zones, three plant functional groups (PFGs) were selected within a 50 m × 200 m sample area in July 2021. They were hygrophyte, mesophyte, and xerophyte. In each PFG, three plant species were selected based on the method that dominance is greater than 50% (e.g., in the 2 m × 2 m sample quadrat, if the dominance of A species is greater than 50%, we will define this observed quadrat as the A plant patch, and select plant A as an investigate object) (Sawchik et al. 2003). In total, nine plant species were selected (Table 1). For each species, three plant patches were investigated as the repetition. And, in each observed plant patch, three plants were selected, totally 81 plants were investigated. Three soil samples were collected at the depth of 0-20 cm around each plant (the black points in Fig. 1) after harvesting the plant leaf samples in July 2021, totally 243 soil samples were harvested. These soils and plant leaf samples were immediately stored at 4 ℃ in cooler boxes and shipped to the laboratory for examining the soil C, N, P concentrations, soil pH, plant leaf C, N, P concentrations, soil microbial biomass C, N and soil enzymatic activity.
2.3 Sample analysis
The fresh leaves of each plant species were oven-dried at 105 ℃ for 20 minutes, and then at 60 ℃ for 48 hours. These dried leaf samples were ground and passed through 0.15 mm sieve for determining the nutrients of plant leaf–TOC, leaf–TN, and leaf–TP. Plant leaf–TOC was measured using the vario TOC (vario TOC select, Elementar Analysensysteme GmbH, Elementar-Straße 1, 63505 Langenselbold, Germany). Leaf–TN was measured by the Kjeldahl method, and the digestion solution was examined using the FOSS Kjeltec TM 8420 (Kjeltec TM 8420, FOSS. Foss Allé 1, DK-3400 Hillerød, Denmark) to examine the N concentrations (Bao 2000). Leaf–TP was determined by the H2SO4-HClO4 ammonium molybdate method (Bao 2000).
In the laboratory, each soil sample was divided into two parts. One was stored at 4 ℃ to determine soil MBC, MBN, and seven soil enzyme activities after passing through a 2 mm sieve. Another soil sample was air-dried to measure soil physicochemical properties. SOC was determined by the vario TOC. Soil TN and TP were measured using the Kjeldahl method and H2SO4-HClO4 method, respectively. Their digestion solutions were examined using the SEAL Auto Analyzer 3 (SEAL Analytical, Inc., 6501 W Donges Bay Rd, Mequon, WI 53092, United States) to determine the concentrations of N and P. Soil pH was measured in a 1: 2.5 soil-DI H2O suspension using a Sartorius pH meter (PB-10, Sartorius Corporate Administration GmbH, Goettingen, Germany).
Soil MBC and MBN were determined by the chloroform (CHCl3) fumigation-extraction method (Brookes et al. 1985). Briefly, two fresh soil samples were weighed (each of 20 g), one was fumigated with CHCl3 for 12 h at 25 ℃. Then, the fumigated soil was extracted with 80 ml 0.5 mol·L-1 K2SO4 and shaken for 30 min at 200 rpm on an oscillator. Another soil sample (non-fumigated) was directly extracted as discussed above. The supernatants of fumigated and non-fumigated soil samples were filtered to determine fumigated and non-fumigated soil MBC and MBN. The MBC and MBN in the filtrate were measured using the elementar vario TOC.
The potential activities of seven extracellular hydrolytic enzymes were determined by the standard fluorimetric techniques using a TECAN Spark® multimode microplate reader (TECAN Spark®, Tecan Group Ltd., Seestrasse 103, 8708 Männedorf, Switzerland ) setting excitation-wavelength to 365 nm and emission-wavelength to 450 nm (Bell et al. 2013, Li et al. 2020). These enzymes degrade the common constituent of organic matter and help microorganisms and plants to acquire C, N, and P nutrients (Li et al. 2020). They are C–acquisition enzymes:αG, βG, βDC, and βX, N–acquisition enzymes: NAG and LAP, and P–acquisition enzyme: aP.
2.4 Statistical analysis
We performed Shapiro-Wilk and Levene’s test to analyze normality and equal variances of our data by R ‘stats’ and ‘car’ packages, respectively, using R: The R Project for Statistical Computing (Bell et al. 2014, R-Core-Team 2013). To improve the homoscedasticity and normality for the subsequent analysis, some data were natural-log transformed such as enzyme C: N = ln (αG + βG + βDC + βX): ln (NAG + LAP), enzyme N: P = ln (NAG + LAP): ln (aP), and enzyme C: P = ln (αG + βG + βDC + βX) : ln (aP). An IBM–SPSS Statistics for Windows version 25.0 (IBM Corp., Armonk, NY, USA) was used to calculate the means and standard errors for all parameters in this study.
We used the distance-based redundancy analysis (dbRDA) as a statistical-criterion to evaluate the overall difference in the plant-soil-microbial C: N: P stoichiometric ratios among plant species and PFGs by the ‘vegan’ package in the R (R-Core-Team 2013). The dbRDA is an excellent nonlinear distance-metric model with the multidimensional resolution, which has been generally accepted by ecologists (Soong et al. 2016). This model is an ordination approach with 3 steps to assess the effect of the environmental variables (i.e. stoichiometric ratio) on the defined group (i.e. plant species, PFGs). Briefly, a dissimilarity or distance matrix was calculated for different plan species and PFGs using the Bray–Curtis dissimilarity (nonlinear) to obtain the dissimilarity or distance matrices for plant species and PFGs (Soong et al. 2016). In steps 2 and 3 of the dbRDA model, a Principal Component Analysis (PCA) was calculated based on the distance matrices, from which the eigenvalues were applied to a Redundancy Analysis (RDA) (Soong et al. 2016). The results of dbRDA were visualized using R: ‘ggplot2’ package.
One-way Analysis of Variance (ANOVA) was used to check the differences in plant leaf stoichiometry (leaf–TOC, leaf–TN, leaf–TP, leaf C: N, leaf N: P, and leaf C: P), soil stoichiometry (SOC, TN, TP, soil C: N, soil N: P, and soil C: P), soil microbial biomass stoichiometry (MBC, MBN, and MBC: MBN), and enzyme acquisition activities stoichiometry (activities of αG, βG, βDC, βX, NAG, LAP, and P, enzyme C: N, enzyme N: P, and enzyme C: P) at different plant species levels and PFGs levels, respectively, using the IBM–SPSS v. 25.0. Fisher’s least significant difference (LSD) test to compare the mean values for all univariate analysis. The results of one-way ANOVA were visualized by Origin Pro v. 2018 (Origin Pro for Windows, v. 2018, Origin Lab Corporation, Northampton, MA, USA).
To assess the relationships between consumer stoichiometry (plant leaf C: N: P) and their resource stoichiometry (i.e. soil C: N: P, MBC: MBN, and enzyme C: N: P) of the riparian zones in Hulunbuir steppe, a regression analysis was calculated in R. To adhere to the conventions of stoichiometric analysis and to normalize variance because of the data of enzymes stoichiometry has been natural-log transformed, other stoichiometric data were also natural-log transformed before performing the regression analysis (Sterner &Elser 2002). The results of regression analysis were visualized using R: ‘ggplot2’ package. Furthermore, the homeostatic index (H) with the corresponding r2 between consumer stoichiometry and resource stoichiometry was calculated by the nonlinear equation y = cx1/H as Bell et al. (2014) described using Origin Pro v. 2018. In this equation, y represents consumer, and x represents resources (Sterner &Elser 2002). H > 1 indicates a lower change in consumer stoichiometry than in resource stoichiometry, which implies there is homeostatic. However, H ≤ 1 indicates non-homeostatic relationships among consumer and resource stoichiometry because there are consistent changes or higher changes in consumer stoichiometry comparing with resource stoichiometry (Bell et al. 2014, Cleveland &Liptzin 2007, Sterner &Elser 2002). Through assessing the homeostasis relationship, we can indirectly evaluate the nutrients limitation among plant, soil, and microorganisms (Sperfeld et al. 2017).