Site description
The study site is a typical semiarid grassland representative of the Eurasian steppe17, located in the Xilin River Basin, Inner Mongolia Autonomous Region of China, close to the Inner Mongolia Grassland Ecosystem Research Station (IMGERS, 43°38′ N, 116°42′ E). Mean annual precipitation is 346 mm, with 60–80% of precipitation falling in the growing season (May to September). Mean annual temperature is 0.3 °C, with mean monthly temperatures ranging from −21.6 °C in January to 19.0 °C in July4. The topography at our experimental site consists of two landscape units (i.e. flat block and sloped block), with elevation ranging from 1200 to 1280 m above sea level, and slopes less than 5°18,19. The soil is classified as dark chestnut (Calcic Chernozem, ISSS Working Group RB, 1998) derived from aeolian sediments18,20. The soil substrate is dominated by sandy loam and loamy sand with more than 50% being fine sand and silt21. At the beginning of the experiment, soil organic carbon and total nitrogen contents were higher in the flat block than in the sloped block (Table 1). Plant species richness and above-ground biomass were also greater in the flat block than in the sloped block, although species composition in terms of relative biomass of common species did not differ between the two systems (Table 1). Leymus chinensis (perennial rhizomatous grass) and Stipa grandis (perennial bunchgrass) are the dominant species in the study area, together accounting for more than 70% of community aboveground biomass. Other dominant species include Cleistogenes squarrosa, Agropyron cristatum, Achnatherum sibiricum, and Carex korshinskyi.
Study design
The experimental area was used for moderate sheep grazing (1.5–3 ewes ha−1 year−1) by local herdsmen until 2003. Afterwards, grass swards recovered for two years before the experiment started20,22. At the end of the growing season in 2004, prior to beginning the experiment, swards in the entire area were cut to 3–5 cm in stubble height23. The experiment was established in June 2005 with split plots in a randomized complete block design (Fig. 1). The study area included two blocks (i.e., flat and sloped blocks), with each block further divided into seven plots. We included flat and sloped blocks because our project was designed to assess the impacts of grazing at spatial scales that are both relevant to land management and that can capture ecosystem and landscape-scale effects of grazing24. It is unrealistic to conduct such a study in an area with no variation in topography. Grazing intensity was randomly assigned to the plots, and each plot was divided into two subplots. The grazing or mowing management regime was randomly assigned into each subplot23. In the grazing regime, there were seven levels of grazing intensity (GI: 0, 1.5, 3.0, 4.5, 6.0, 7.5 and 9.0 sheep ha−1), and sheep grazed in the subplots continuously from June to September each year25. The ungrazed plots (0 sheep ha−1) had no sheep grazing for 12 years. Each subplot was 2 ha, except the subplot with 1.5 sheep ha−1, which was enlarged to 4 ha to ensure a minimum herd of six sheep per subplot. In the mowing regime, mowing was done once a year in the middle of August. Plant and soil microbial community data was collected in late July and early August 2017, after 12 years of grazing and mowing treatments.
Illustration of the grazing experiment design. G: grazing regime, M: mowing regime.
Plant community surveys
For each subplot, we randomly laid out ten 1 m × 1 m quadrats at least five meters from the edge of each plot to avoid edge effects. In each quadrat, plant species were identified, and the abundance of each species was counted by bunches (bunchgrasses) or stems (rhizomatous grasses). For each species, five individuals were randomly chosen to measure plant height and the average height of all species was used as plant canopy height. Plant canopy coverage was measured visually.
Soil sampling
For each quadrat, three soil cores (3 cm diameter, 10 cm depth) were collected, and soil was passed through a 2 mm sieve to form one composite soil sample per quadrat. Sieved soil was then divided into three subsamples. One subsample was air-dried for the analysis of soil pH, soil organic carbon (SOC), total nitrogen (TN), and total phosphorus (TP). The second fresh subsample was used for the analysis of microbial community structure and microbial biomass. The third subsample was stored at -20 °C prior to being used for microbial sequencing analysis.
Soil physical and chemical properties
To evaluate soil compaction, we measured soil hardness by using a Yamanaka-style soil hardness tester (Fujiwara Scientific Co., Japan). Soil moisture content was measured by using 10 g of moist soil that was oven-dried at 105 °C for 24 h. Soil pH was measured in a 1:2.5 soil:water suspension using a pH meter (FE20-FiveEasy, Mettler-Toledo, Switzerland).
We measured SOC content with the Walkley-black method, soil TN content by the micro-Kjeldahl digestion, followed by colorimetric determination with a 2300 Kjeltec Analyzer Unit, and soil TP content was by the H2SO4-HClO4 fusion method using a 6505 UV spectrophotometer26.
Soil microbial community structure
Microbial community structure was assessed using phospholipid fatty acids (PLFAs), as described by Bossio and Scow27. First, lipids were extracted from 10 g of fresh soil using a buffer (CHCL3:CH3OH:K2HPO4 = 1:2:0.8, v:v:v). Second, the fatty acid methyl esters (FAMEs) were separated, quantified and identified using a gas chromatograph system (Agilent 7890, Santa Clara, USA) and a MIDI Sherlock Microbial Identification System (MIDI Inc., Newark, USA). Peak areas were converted to nmol g−1 dry soil using the internal standard, methylnon-adecanoate (C19:0). Third, the specific microbial groups were identified according to their representative markers. Specifically, G+ bacteria correspond to iso-, anteiso- and 10Me-branched PLFAs; G- bacteria correspond to monounsaturated and cyclopropyl PLFAs; arbuscular mycorrhizal fungi (AMF) use 16:1ω5c as representative marker; saprotrophic fungi (SF) use 18:1ω9c, 18:2ω6c and 18:3ω6c as representative markers28,29,30. The 12:0, 14:0, 15:0, 16:0, 17:0, 18:0 PLFAs were general markers present in all microorganisms30,31. Bacterial PLFAs included G+ and G− bacteria PLFAs. Fungal PLFAs included arbuscular mycorrhizal and saprotrophic fungi PLFAs. Total microbial PLFAs were the sum of bacterial, fungal, and general PLFAs.
Soil microbial biomass carbon (MBC), nitrogen (MBN), and phosphorus (MBP) were measured using the chloroform-extraction method32,33. For MBC and MBN, two fresh soil samples were used for the analysis. One sample was placed in a chloroform steam bath for 24 h and another sample was kept non-fumigated. Then, organic C and total N were extracted by shaking two soil samples in 0.5 M K2SO4 for 1 h and filtering through a Whatman No. 1 filter paper (9 cm in diameter). The filtered extracts were measured with a total organic carbon (TOC) analyzer (Elementar vario TOC, Hanau, Germany). Microbial biomass P was measured using a similar method as for MBC and MBN except that P was extracted by 0.5 M NaHCO3 and then measured with a UV Spectrometer (6505 spectrometer, Jenway, Stone, UK).
DNA extraction and sequencing
We mixed ten soil samples of each plot to form one composite sample for DNA extraction and sequencing. Total genomic DNA was extracted from 0.5 g soil using a FastDNA Spin kit (MP Biomedical, Santa Ana, California, USA). The DNA quality was checked by 1% agarose gel electrophoresis and quantity was determined with a NanoDrop 2000 UV-vis spectrophotometer (Thermo Scientific, Wilmington). Bacterial 16 S rRNA genes were amplified with PCR primers 338 F (5′- ACTCCTACGGGAGGCAGCAG-3′) and 806 R (5′-GGACTACHVGGGTWTCTAAT-3′). Fungal internal transcribed spacer (ITS) rRNA genes were amplified with PCR primers ITS1F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and ITS2 (5′-GCTGCGTTCTTCATCGATGC-3′)34,35. The resulting PCR products were extracted from a 2% agarose gel and further purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA) and quantified using QuantiFluor™-ST (Promega, USA). Purified amplicons were pooled in equimolar concentrations and paired-end sequenced for high-throughput 16 S rRNA or ITS rRNA gene sequencing on an Illumina Hiseq. 2500 platform (Illumina Inc., USA) according to the standard protocols by Novogene Technology Co., Ltd. Operational taxonomic units (OTUs) were clustered with 97% similarity cut-off using UPARSE (version 7.1 https://drive5.com/uparse/), and chimeric sequences were identified and removed using UCHIME. Silva and Unite databases were used as references for bacteria and fungi, respectively34,35.
Source: Ecology - nature.com