Study area
We carried out this research work in Gansu Province (China) in Minqin County, (103°07′ 16″E, 38°37′10″N) situated on the Shiyang River Basin’s lower reaches in the Hexi Corridor, in an area bound by the Tenggeli and Badanjaran Deserts (Feng et al. 2011). In this area a typical arid continental climate prevails, whose mean annual air temperature and precipitation is respectively 7.8℃ and 113.2 mm (most of the latter falling between July to September), and an average annual evaporation totaling 2646 mm. The uneven rainfall distribution causes intense seasonal droughts in spring and/or winter. According to the, The research site has an irrigated desert soil (Gansu Provincial Soil Survey Office. 1992), this not unlike similar to Anthropic Camborthids according to Soil Taxonomy (Jiang et al. 2017). Feng et al. (2011) describe other relevant environmental conditions and agriculture production in this study area.
Design of the experimental and soil sampling
In 2015 the experiment was installed, by flat planting the field with spring maize (Zea mays L.) early on in April. The crop was allowed to grow until October’s end, when it was harvested. During each winter, all plots in the field were left fallow. When the experiment was begun, the soil pH was 8.63, and it 0.29 ms cm–1 specific conductance, 9.80 g kg–1 organic matter, 4.84 g kg–1 of total N, 0.40 mg kg–1 of available P, and 144.14 mg kg–1 of available K. An RBD (randomized block design) was used for this maize experiment, in which tillage was applied in four treatments, each having three replicates (i.e., plots). The tillage systems mainly differed in the cultivation method practiced each autumn after harvesting the main crop: 1) conventional sub-tillage (ST): moldboard plow down to 30 cm; 2) fold-tillage with sub-soiling (FT): chisel plowing to 20 cm; 3) conservation minimal tillage (MT): no autumn cultivation and rotary tillage to 10 cm; 4) conservation no-tillage (NT) (Fig. 1). Supplemental irrigation of all maize plots was applied using basin irrigation at the seedling, jointing, pre-heading, silking, flowering, and grain-filling stages, respectively. Plots were 30 m × 20 m, separated by concrete borders.
The plow layer (0–20 cm depth) was sampled for soil, by using10 cm diameter auger inserted at five position (these chosen randomly) within given crop row, about 20 cm from maize stems during their harvest in October 2021; these soil sampled were then mixed to into a composite sample at the plot level. These freshly collected soil samples were stored in containers made of hard plastic; to safeguard their soil structure from destruction while in transit to laboratory, each container was filled with buffer foam pieces.
Soil aggregate fractionation
To isolate the aggregates, a modified dry-sieving method was used (Dorodnikov et al. 2009) to minimize the disturbance of soil microbes. In a well-ventilated place, sample was dried to ca. 10% gravimetric soil water content (‘optimal’ moisture), to limit the mechanical stress “to induce maximum, brittle failure along natural planes of weakness” (Dorodnikov et al. 2009). That is, along their natural breakpoints, individual soil clods were gently separated, then sieved through 8-mm mesh to remove any visible stones (and other debris). Then, each soil sample was divided into two parts: the first was for bulk soil analysis and aggregate fractionation, by using the following sieving procedure. Around 1 kg soil per plot was positioned over two sieves (2.0 and 0.25 mm mesh) and gently shaken thrice, for 2 min each time. Soil material that was left on the sieves was collected and respectively designated as large macro-aggregates (i.e., > 2 mm) and small macro-aggregates (i.e., 0.25–2 mm). Any remaining soil which passed through the 0.25-mm sieve was then collected, this being the micro-aggregate particles (hereon ‘micro-aggregates’) (i.e., < 0.25 mm). Next, we weighed each soil aggregate fraction and divided each into two portions (1) for analysis of available nutrients, microbial community biomass, diversity and composition, and soil enzymatic activity assays; and (2) for air drying prior to its analysis of soil properties (chemical).
To obtain the soil aggregate distribution, the respective proportions corresponding to each soil aggregate fraction were calculated, with the mean weight diameter (MWD) used an index of aggregate stability: MWD = \(\sum _{i}^{3}xi\times wi,\) where i denotes the number count of aggregate size fractions; xi denotes the mean diameter (mm) of the aggregates in any size fraction; and wi denotes the weight of the aggregates in that particular size fraction, expressed as a proportion of the total soil dry weight (Somasundaram et al. 2017).
Chemical and biological properties of soil
Soil chemical analysis
To determined soil moisture, the difference in mass before and after drying at 105℃ to constant weight was used. To quantify soil organic C (SOC), the Walkley and Black dichromate oxidation method was used with 1.08 as the factor. For total N (TN), it was measured by the K2SO4–CuSO4–Se (100:10:1) distillation method; for total P (TP), the molybdate colorimetric method after perchloric acid digestion and ascorbic acid reduction was used. After first soaking them in 2 mol L–1 KCl, soil samples were shaken (200 rpm, 1 h) and analyzed for their NH4+-N and NO3–-N by an auto-flow injection system (Auto-Analyzer AA3, Germany). To quantify the available P (AP) in soil, we applied the Olsen-P method (Yuan et al. 2016). To determine the microbial biomass C (MBC) or N (MBN) in soil, the difference in values between the C or N extracted via chloroform fumigated vis-à-vis non-fumigated soil samples was obtained, to which the KEC factor (0.45) or KEN (0.57) was respectively applied. Briefly, 20 g fresh soil was fumigated for 24 h at 25°C with CHCl3 and then added to 80 mL of 0.5 mol L–1 K2SO4, shaken (200 rpm, 1 h) and filtered after removing the fumigant. A 20 g sample of non-fumigated soil was extracted in tandem. To determine the OC and TN contents of each extract a total C analyzer was used (TOC-L CPH Basic Analyzer System, Japan). In the non-fumigated samples, their organic C concentration corresponded to the dissolved organic C (DOC) (Wallenstein et al. 2006).
Soil enzyme assays
Four hydrolytic soil enzymes were quantified for their respectively activity. These consisted of two enzymes participating in C-acquisition (β-1,4-glucosidase and cellobiohyrolase), one enzyme participating in N-acquisition (β-1,4-N-acetylglucosaminidase), and one enzyme participating in P-acquisition (alkaline phosphatase) (Sinsabaugh et al. 2009). We also quantified the activities of two oxidases, phenoloxidase and peroxidase, which enable the degradation of recalcitrant organic C (Sinsabaugh 2010). To measure hydrolytic enzymes, a fluorescence microplate assay was employed, in which the substrates were labeled with 4-methylumbelliferone (MUB) (Bell et al. 2013). To do this, 1.25 g of fresh soil and 125 mL of 50 mmol L–1 Tris buffer solution (pH 8.30) were mixed at 25°C for 20 min by an oscillator, this then serving as a soil buffer suspension. Next, 400 µL of that suspension and 100 µL of substrate solution (200 µmol L–1) were pipetted into a 96-well microplate. Meanwhile, 400 µL of each sample mixture, a negative control, the reference standard, a blank control, and a quench control were pipetted into a black 96-well plate and incubated at 25°C for 3 h in the dark; to end the reaction, 20 µL of NaOH was added. Activity levels were measured fluorometrically, by a fluorescence plate reader (Varioskan LUX, ThermoFisher Scientific, USA) at excitation and emission of 365 and 450 nm, respectively, and expressed in units of substrate converted (µmol) per mL of sample (µmol g− 1 h− 1). Oxidases were measured spectrophotometrically in a clear 96-well microplate using the substrate L-3, 4-dihydroxyphenylalanine (DOPA). To represent the potential capacity of these C-, N-, and P-acquisition activities, their corresponding enzymes were grouped accordingly, and their data transformed (normalized) as follows: respectively, ln(β-glucosidase + cellobiohyrolase + phenoloxidase + peroxidase), ln(β-N-acetylglucosaminidase), and ln(phosphatase) (Sinsabaugh et al. 2009).
Soil microbial community composition
Soil DNA extraction and sequencing
From each soil sample (0.5 g) its microbial DNA was extracted by using the HiPure Soil DNA Kits (Magen, Guangzhou, China). The amplicons were assessed using 2% agarose gel electrophoresis and purified with an AxyPrep DNA Gel Extraction Kit (Axygen Bioscience, Union City, USA). From the bacterial 16S rDNA target region of the rRNA gene, the V3–V4 hypervariable region was PCR-amplified (conditions: 95°C for 2 min, then 35 cycles of 95°C for 30 s, 60°C for 45 s, and 72°C for 90 s; ending with a final extension 72°C for 10 min), using the primers 515F: GTGYCAGGMGCCGCGGTAA; 806 R: GGACTACNVGGGTWTCTAAT. To amplify the fungal ITS2 region, the primers ITS3-KYO2: GATGAAGAACGYAGYRAA and ITS4: TCCTCCGCTTATTGATATGC were used under the same condition. Each PCR reaction was performed in a 50-µL reaction volume that included TransGen High-Fidelity PCR SuperMix (TransGen Biotech, Beijing, China), 0.2 µM each F/R primer, and 5 ng of template DNA. To detect amplicons, 2% agarose gels were used, followed by purification with the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA). To generate the sequencing libraries, SMRTbell TM Template Prep Kit (PacBio, Menlo Park, CA, USA) was used; to assess library quality, we used the Qubit 3.0 Fluorometer (ThermoFisher Scientific, USA) and FEMTO Pulse system (Agilent Technologies, Santa Clara, CA, USA). All the obtained sequence data has been deposited into the CNGB Sequence Archive (CNSA) of China National GeneBank DataBase (CNGBdb), where it can be found under accession number CNP0003819.
Bioinformatics analysis
To merge the paired-end clean reads (raw) FLASH (v 1.2.11) tool (Magoč and Salzberg 2011) was used with 10-bp minimal overlap and 2% mismatch error rate at most. Noisy reads sequences were screened by prescribed filtering conditions to obtain a subset of clean reads high in quality, as described by Bokulich et al. (2013). The retained clean reads were clustered into operational taxonomic units (OTUs) have ≥ 97% similarity in identity, by implementing the UPARSE pipeline (v 9.2.64) (Edgar 2013). Any chimeric sequences present were found and excised by running the UCHIME algorithm (Edgar et al. 2011). Bacteria were identified taxonomically by the RDP classifier (Wang et al. 2007) according to the SILVA database (Prüesse et al. 2007), and the fungi were classified using the UNITE database (Nilsson et al. 2018) with the BLAST tool. A suite of alpha diversity metrics, namely the observed number of species (richness), Pielou’s evenness, in addition to the well-known Shannon–Wiener index, as well as the Chao1 estimator, were derived for the communities of bacteria or fungi, using QIIME. Bioinformatics analyses were carried out in R software (v 2.2.1).
Statistical analyses
To assess the tillage intensity (regime) effect on aggregate stability (MWD), properties of soil, biomass of microbial, enzymatic activity and stoichiometry, aggregate size fraction, and microbial diversity, we used ANOVA (one-way model) followed by pairwise comparisons using Tukey’s HSD test. In addition, the single effects of each size fraction of soil aggregate on enzymatic activity and stoichiometry as well as soil properties under each tillage treatment were assessed by univariate ANOVAs (one-way). Statistical significance was designated using an alpha = 0.05 (i.e., P ˂ 0.05). All above ANOVAs and multiple mean comparisons were implemented in Genstat 18 software (VSN International Ltd., Rothamsted, UK); to draw the figures, we used Origin 9.2 software (OriginLab 2016, USA). Principal component analysis (PCA) was applied to based on OUT levels data, to examine in what way the structures of bacterial and fungal communities were related or differed between aggregate size fractions and tillage management regimes.