- 1Heilongjiang Black Soil Conservation and Utilization Research Institute, Harbin, China
- 2The Centre for Ion Beam Bioengineering Green Agriculture, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei, China
- 3Science Island Branch, Graduate School of USTC, Hefei, China
- 4Keshan Branch of Heilongjiang Academy of Agricultural Sciences, Qiqihar, China
Conservation tillage (CT) is an important agronomic measure that facilitates soil organic carbon (SOC) accumulation by reducing soil disturbance and plant residue mulching, thus increasing crop yields, improving soil fertility and achieving C neutrality. However, our understanding of the microbial mechanism underlying SOC fraction accumulation under different tillage practices is still lacking. Here, a 6-year in situ field experiment was carried out to explore the effects of CT and traditional tillage (CK) practices on SOC fractions in an eolian sandy soil. Compared with CK, CT increased the particulate OC (POC) content in the 0–30 cm soil layer and the mineral-associated OC (MAOC) content in the 0–20 cm soil layer. Moreover, tillage type and soil depth had significant influences on the bacterial, fungal and protistan community compositions and structures. The co-occurrence network was divided into 4 ecological modules, and module 1 exhibited significant correlations with the POC and MOC contents. After determining their topological roles, we identified the keystone taxa in the network. The results indicated that the most common bacterial taxa may result in SOC loss due to low C use efficiency, while specific fungal (Cephalotrichum) and protistan (Cercozoa) species could facilitate SOC fraction accumulation by promoting macroaggregate formation and predation. Therefore, the increase in keystone fungi and protists, as well as the reduction in bacteria, drove module 1 community function, which in turn promoted SOC sequestration under CT. These results strengthen our understanding of microbial functions in the accrual of SOC fractions, which contributes to the development of conservation agriculture on the Northeast China Plain.
1 Introduction
Soils are the largest carbon (C) pool in the global terrestrial system and contain more than 2,500 Gt of C (Banerjee et al., 2016). Global soil organic carbon (SOC) dynamics immensely influence soil productivity, greenhouse gas emissions and C neutrality (Tang et al., 2019; Duan et al., 2023). In agroecosystems, tillage practices are regarded as crucial agronomic regimes that mediate SOC sequestration and depletion processes (Topa et al., 2021). Traditional agricultural practices involve frequent tillage to accomplish achieve high crop productivity (Zhang et al., 2019). However, high-intensity tillage not only decreases crop productivity but also leads to reduced soil fertility and sustainability and results in the loss of other agroecosystem services, causing soil erosion, water shortages or biodiversity decline (Raus et al., 2016). Hence, abundant attempts have been made to transition from traditional tillage to conservation tillage (CT) to increase the soil C stock of agroecosystems in recent years. In particular, CT practices, such as no-tillage or reduced-tillage, may minimize the degree and frequency of tillage passes and maintain an adequate soil surface covered with residues to reduce soil physical disturbance and increase the soil C sink capacity (Pearsons et al., 2023). Nevertheless, contrasting tillage practices cause changes in resource availability in the topsoil and subsoil, which leads to differences in SOC formation (Angers and Eriksen-Hamel, 2008). Therefore, elucidating the SOC sequestration mechanisms that occur under different tillage practices and at different soil depths is crucial for maintaining soil health and facilitating agroecosystem services.
Overall, the input of exogenous organic materials (i.e., manure and crop residues) is an essential prerequisite for SOC accumulation (Lehmann and Kleber, 2015). Many expert researchers have confirmed that the mechanisms underlying SOC accumulation are commonly attributed to physical protection by aggregates and chemical stabilization by soil minerals (Six et al., 2004). Correspondingly, semidecomposed exogenous large organic biopolymers are readily encapsulated by aggregates and form particulate organic carbon (POC); as biopolymers further decompose, the C monomers tend to be adsorbed by mineral surfaces and become mineral-associated organic carbon (MAOC) (Bastian et al., 2009; Herath et al., 2014). There is nearly a consensus regarding the disruption of topsoil aggregates due to frequent tillage, which leads to the loss of POC (Hewins et al., 2017). These conclusions also reflect the potential of using CT to increase SOC sequestration in agroecosystems. Several previous studies also suggested that traditional tillage practices involving straw return can transport crop residues to the subsoil, thus contributing to the accumulation of SOC in the subsoil to some extent (Zhang et al., 2013). This contradiction also illustrates the complexity of enhancing soil fertility through tillage practices. Achankeng and Cornelis (2023) found that climate, soil texture, rotation pattern, and crop type are all crucial factors to be considered under different tillage treatments in a meta-analysis, which increases the challenge for researchers in optimizing tillage practices. Therefore, to date, we still lack a comprehensive understanding of the direct associations between tillage practices and SOC fractions.
Soil microorganisms play an important role in SOC formation (Cotrufo et al., 2013; Sarker et al., 2018). Generally, when plant residues are applied, soil animal- and meso-fauna-driven fragmentation constitute the first stage of straw degradation (Gessner et al., 2010). Subsequently, bacteria and fungi successively regulate further degradation processes due to changes in C and N availability in the substrate (Wang et al., 2021). In the early stage of decomposition, adequate amounts of labile straw C and nitrogen can sustain bacterial proliferation (Huang et al., 2017). With continuous plant decomposition, the microbial community synchronously undergoes succession. Fungi may be the dominant decomposers due to their potent ability to utilize recalcitrant straw components (Clemmensen et al., 2015). Therefore, many scholars consider the ratio of bacteria to fungi to be an important indicator of the straw decomposition process (Li J. W. et al., 2020; Zhao et al., 2021). Moreover, microbial diversity and module community was the key drivers of SOC turnover. Previous study found that bacterial, fungal protistan richness was significantly correlated with carbon use efficiency, microbial biomass carbon, microbial respiration and growth rate, which changed the SOC turnover process (Ma et al., 2024). Soil microbial module community also played an important role in influencing SOC. Numerous studies have demonstrated that the core microbial module community was involved in maintaining the stability of soil microbial function, promoting soil nutrient cycling and SOC accumulation; while the keystone species were the core to achieve them (Shi et al., 2020; Jiao et al., 2022). Therefore, keystone taxa-driven the changes of microbial module communities and diversity are pivotal factors leading to SOC turnover. In particular, soil bacterial and fungal communities are extremely sensitive to tillage practices. Li Y. et al. (2020) reported that, compared with traditional tillage, minimum tillage increases fungal biomass and bacterial diversity, which may further influence residue decomposition. However, as major members of the soil microbiome, protists drive plant residue decomposition, and microbial community regulation has rarely been included in microbiome analyses associated with SOC fraction turnover (Geisen and Bonkowski, 2018). Specific protozoan taxa participate in aggregate formation and SOC turnover. According to the report of Pellegrino et al. (2021), CT increased the abundance of Alveolata and Cercozoa, which contributed to SOC accumulation by reshaping soil aggregates. Additionally, the top-down control of protists in the soil microfood web demonstrated great potential for influencing SOC turnover (Gao et al., 2019). Therefore, a thorough empirical understanding of the microbial roles (including bacterial, fungal and protist roles) in SOC fraction sequestration under different tillage practices has not been achieved.
To bridge these gaps, eolian sandy soil located in the Northeast Plain, the largest grain-producing area in China, was selected as the research object. In recent years, the Northeast Plain has been facing continuous depletion of SOC stocks since the 1980s, when straw return was widely implemented (Wang et al., 2018). Therefore, we conducted a 5-year in situ field experiment to reveal the effect of tillage practices on the microbial community, SOC fraction and maize yield at different soil depths. The soil samples were collected from the 0–50 cm soil profile under CT and traditional tillage practices. In this study, we determined the SOC fraction content and microbial community and attempted to explain the potential relationships between them. We hypothesized that (1) the SOC fraction content and microbial traits exhibit distinct responses to traditional and CT practices at different depths and that (2) specific microbial taxa may be involved in the turnover of SOC fractions.
2 Materials and methods
2.1 Experimental site
The experimental field was located in Dulbert Mongolian Autonomous County (46°54′N, 124°26′E), Daqing city, Heilongjiang Province, which has a semiarid monsoon continental climate. The mean annual precipitation and temperature are 400 mm and 5.6°C, respectively. According to the USDA soil taxonomy, the soils in the area are carbonate meadow soils. The basic nutrient contents of the soil before the experiment were as follows: 0.53 g kg−1 total nitrogen; 60.95 mg kg−1 alkali-hydrolyzable nitrogen; 60.22 mg kg−1 available phosphorus; 44.51 mg kg−1 available potassium; and 9.52 g kg−1 SOC with a pH of 5.56. The cropping system used was single spring maize (Zea mays L.).
2.2 Field trial design and soil sampling
The experimental trial was set up in 2017 in accordance with a randomized complete block design with three replicates. Each field plot was covering an area of 64 m2 (4 m × 16 m). Before the experiment, all the plots were treated with N-P-K fertilizers, and straw was removed via shallow tillage to 25 cm. Traditional tillage (CK): the plots were plowed with large machinery, and a five-share turning plow tilled the soil to 25 cm after harvest at October (amount of maize stover return about 7,200 kg ha−1); and CT: no-tillage with 100% straw mulch after harvest at October (amount of maize stover mulch about 7,600 kg ha−1). A straw crusher was used to crush the straw into fragments with lengths less than 10 cm before mulching. Except for the field surface drilling of maize in October, the no-tillage plots remained undisturbed, and maize straw was evenly distributed over the field surface after harvest every year. Chemical fertilizers were applied in May, and the N-P2O5-K2O application rates ranged from 180–115–75 kg hm−2. All the other normal management practices were consistent between the treatments during the experiment.
Soil profiles were excavated at a depth of 50 cm in October 2022 in each replicate plot. Soil samples were collected at depths of 0–10, 10–20, 20–30, 30–40, and 40–50 cm. Three soil samples were collected from each plot. The three soil samples were placed into the same sterile plastic bag and mixed to create one composite sample. All the samples were immediately transported to the laboratory in an incubator with ice packs. Each soil sample was divided into two parts: one part was stored at −80°C for DNA extraction, and the remaining part was air-dried for use in the additional chemical analyses. The basic chemical properties of the soils under different tillage practices and at different depths in 2022 are shown in Supplementary Figure S1.
2.3 Basic soil properties, SOC fractionation, and soil aggregate isolation
The basic chemical properties of the soil were measured using the method described by Lu (2000). The soil pH was measured at a 1:2.5 soil: water ratio for 30 min. The SOC concentration was determined using K2Cr2O7 digestion, and total nitrogen was determined by the Kjeldahl method. Available phosphorus and potassium were determined using molybdenum blue colorimetric and flame photometry methods, respectively.
SOC fractionation was determined using a method described by Yu et al. (2017). First, the SOC was further fractionated into POC and MAOC. Generally, 5.0 g (dry weight) of soil was dispersed by adding 30 mL of 0.5% sodium hexametaphosphate solution and centrifuging at 200 r min−1 for 18 h. Thereafter, the POC and MAOC fractions were obtained by passing the samples through 53-μm filters. All the fractions were dried (50°C), and POC and MAOC were measured using K2Cr2O7 digestion.
The wet sieving method was used to separate the water-stable aggregates (Six et al., 2002). First, the soil samples were gently broken apart into small pieces along natural break points, and the fragmented samples were subsequently placed on top of a 0.25-mm sieve and soaked in deionized water for 5 min. The samples were subsequently separated at an amplitude of 3 cm and a frequency of 30 cycles per min for a duration of 2 min by a wet-sieving apparatus. Afterward, the aggregate subsamples above each sieve were obtained as follows: macroaggregates (>0.25 mm), microaggregates (0.053–0.25 mm), and silt and clay fractions (<0.053 mm). After wet sieving, all the aggregates were immediately oven-dried at 60°C and weighed.
The mean weight diameter (MWD) was used to describe the aggregate stability and was calculated by the following formula:
where Xi represents the average diameter of each aggregate size and Wi represents the proportion of each aggregate weight relative to the total sample weight after wet sieving. The upper limit of the macroaggregate diameter was 2 mm.
2.4 DNA extraction and 16S, ITS and 18S amplification and sequencing
Total DNA was extracted from 0.5 g of soil using a Fast DNA Spin Kit for Soil (MP Biomedicals, CA, United States) in accordance with the manufacturer’s instructions. Each treatment contained three replicates. The extracted DNA samples were stored at −80°C for molecular analysis.
High-throughput sequencing was performed using the Illumina MiSeq sequencing platform (Illumina, Inc.). Both the forward and reverse primers were tagged with adapter and linker sequences, and 8-bp barcode oligonucleotides were added to distinguish the amplicons derived from different soil samples.
The primers 515F (5′-GTGCCAGCMGCCGCGGTAA-3′) and 907R (5′-CCGTCAATTCMTTTRAGTTT-3′) were chosen to amplify the 16S rRNA genes in the V4–V5 hypervariable region. PCR was conducted in a 50-μL reaction mixture containing 27 μL of ddH2O, 2 μL (5 μM) of each forward/reverse primer, 2.5 μL (10 ng) of template DNA, 5 μL (2.5 mM) of deoxynucleoside triphosphates, 10 μL of 5× Fastpfu buffer, 0.5 μL of bovine serum albumin, and 1 μL of TransStart Fastpfu polymerase (TransGen, Beijing, China). The PCR procedure was 94°C for 5 min; 30 cycles of 94°C for 30 s, 52°C for 30 s and 72°C for 30 s, followed by 72°C for 10 min (Biddle et al., 2008).
The fungal ITS1 region was amplified using the primer pair ITS1F (CTTGGTCATTTAGAGGAAGTAA)/ITS2 (GCTGCGTTCTTCATCGATGC). The 50-μL reaction mixture contained 1 μL (30 ng) of DNA, 4 μL (1 μM) each of forward/reverse primer, 25 μL of PCR Master Mix, and 16 μL of ddH2O. PCR amplification was conducted at 98°C for 3 min, followed by 30 cycles of 98°C for 45 s, 55°C for 45 s, and 72°C for 45 s, with a final extension at 72°C for 7 min (Ghannoum et al., 2010).
The eukaryotic V4 region was amplified using the primer pair V4_1f (CCAGCASCYGCGGTAATWCC)/TAReukREV3 (ACTTTCGTTCTTGATYRA). PCR was performed in a 20 μL volume consisting of 4 μL of 5× reaction buffer, 2 μL of dNTPs (2.5 mM), 0.8 μL of each primer (10 μM), 0.4 μL of FastPfu Polymerase, 10 ng of DNA template, and ddH2O to reach the final volume. PCR amplification was conducted at 95°C for 5 min, followed by 30 cycles of 95°C for 30 s, 55°C for 30 s, and 72°C for 45 s, with a final extension at 72°C for 10 min. To construct the protistan amplicon sequence variant (ASV) table, we removed sequences belonging to Rhodophyta, Streptophyta, Metazoa, and Fungi (Stoeck et al., 2010).
Raw Illumina amplicon reads were processed using the QIIME2 Core 2019.7 distribution. The Divisive Amplicon Denoising Algorithm (DADA2) pipeline implemented in the QIIME 2 platform was used to conduct sequence quality control, which included quality filtering reads, denoising reads, merging forward and reverse reads, removing chimeric reads, and assigning reads to ASVs. The Silva 138 Bacterial 16S rRNA gene database, the UNITE Fungal ITS database, and the Protist Ribosomal Reference (PR2) database v4.14.0 were used to classify the representative sequences of ASVs (Ghannoum et al., 2010). All singletons and nonfungal ASVs were removed, and each sample was rarefied to 30,000, 36,000, and 44,000 sequences for the bacterial, fungal and eukaryotic diversity analysis. The alpha diversity and Bray–Curtis distances for principal coordinate analysis of the soil microbial community were calculated after all the samples were rarefied to the same sequencing depth.
2.5 Statistical analysis
Crop yield and soil biochemical and other relevant properties under different tillage practices were subjected to the chi-square test for independence of variance. Significant differences were determined by one-way analysis of variance (ANOVA) based on the post hoc Tukey test at the 5% level. Prior to ANOVA, the normality and homogeneity of variance were tested by the Kolmogorov–Smirnov test and Levene’s test, respectively. If normality was not met, log or square-root transformation was carried out. One-way ANOVA was performed using SPSS 21.0 (SPSS, Inc., Chicago, IL, United States).
Principal component analysis (PCA) was used to determine and evaluate the changes in the community structure of the soil microbiome via the R (ver 4.2.3) package “vegan.” To characterize the patterns of soil microbial interactions, we constructed a co-occurrence network with the “igraph” and “WGCNA” R packages. We constructed microbial networks using bacteria, fungi and protists with relative abundances greater than 0.01%; screened nodes with Pearson’s correlations greater than 0.6 and p < 0.05; performed modular analysis based on the connectivity between nodes; visualized the network using Gephi (ver. 0.9.2); and calculated information on network topological features. The within-cluster connectivity (Zi) and among-cluster connectivity (Pi) of different clusters were calculated using the R packages “reshape2,” “igraph,” “ggrepel,” “dplyr,” and “Rcpp” and filtered for peripherals (Zi ≤ 2.5, Pi ≤ 0.62), connectors (Zi ≤ 2.5, Pi > 0.62), cluster hubs (Zi > 2.5, Pi ≤ 0.62), and network hubs (Zi > 2.5, Pi > 0.62) (Deng et al., 2012). ASVs in module hubs, connectors and network hubs may be regarded as the microbial keystone taxa of network systems (Deng et al., 2012). An interactive platform “Gephi” (default parameters set) was used to identify the modules (ecological clusters) of soil taxa strongly interacting with each other.
Linear regressions between SOC fractions and the microbial community modules were conducted to determine the relationships between microbial communities and SOC fraction contents using Origin 2018. The microbial module community variation data were quantified by the PCA 1 axis. Heatmaps were constructed to reveal the potential associations between keystone taxa richness and SOC fraction content via the “heatmap.2” function in the R package “ggplots.”
3 Results
3.1 Crop yields, aggregate size distributions, and SOC fractions
Overall, CT had a positive effect on crop yield and SOC fraction content. After 5 consecutive years of different tillage practices, the maize yield, aggregate stability and SOC fraction content changed significantly (Supplementary Figures S2, S3); Figure 1). However, no significant changes were found before 2021. Compared with those in CK, the maize yields in CT significantly increased in 2021 and 2022 (Supplementary Figure S2, p < 0.05).
Figure 1. Soil organic carbon (A), particulate organic carbon (B), and mineral-associated organic carbon (C) contents at various depths under different soil tillage treatments. *p < 0.05; CT, conservation tillage; CK, traditional tillage.
Different tillage practices changed the aggregate size distribution and stability (Supplementary Figure S3). CT significantly increased the proportion of macroaggregates at the 0–30 cm soil depth (p < 0.05) and decreased the proportion of microaggregates at the 0–20 cm soil depth (p < 0.05). However, no changes were observed in the silt and clay fractions under the different treatments. Additionally, the MWD was greater under CT than CK at 0–30 cm depth, while there were no significant changes at 40–50 cm depth.
SOC, POC, and MAOC contents were also affected by tillage practice (Figure 1). Generally, CT significantly increased the SOC, POC, and MAOC contents in the topsoil (0–20 cm layer, p < 0.05, except for MAOC in the 10–20 cm layer). Moreover, the POC content was significantly greater under CT than CK. No significant changes were observed in the 30–50 cm layer under the different tillage practices.
3.2 Microbial community, co-occurrence network, and keystone taxa
PCA was used to evaluate the changes in the soil microbial community under the different tillage practices (Supplementary Figures S4, S5). The results indicated that the soil communities changed significantly under the different tillage practices and at different soil depths (except for tillage practices on the protistan community). Generally, the effect of soil depth on microbial communities was observed mainly along the PCA 1 axis, while the effect of tillage practices on microbial communities was observed mainly along the PCA 2 axis.
Although soil microbial community compositions were changed after different tillage practices and depths, Proteobacteria, Acidobacteriota, Actinobacteriota, and Gemmatimonadota were the main phyla of bacteria, contributing more than 60% of the total bacterial abundance (Supplementary Figure S6A). Generally, with the increase of depth, the relative abundance of Proteobacteria decreased gradually. Ascmycota, Basidiomycota, and Mortierellomycota were the main phyla of fungi, contributing almost 80% of the total fungal abundance (Supplementary Figure S6B); while protist were composed mainly of Intramacronucleata, Cercozoa, Chlorophyta, and Apicomplexa (Supplementary Figure S6C).
A co-occurrence network was constructed to reveal the connections between specific microbial species (Figure 2). We found four dominant ecological modules (Figure 2A). Modules 1, 2, 3, and 4 contained 225, 182, 160, and 124 nodes, respectively. Among the four modules, protists and fungi accounted for the greatest proportion of microbial species in module 1, while bacteria accounted for the highest proportion of microbial species in module 3 (Figure 2B). Additionally, the percentage of intraspecies edges was greater in module 3 than in the other modules, while module 1 contained more interspecies edges than did the other modules (Figure 2C).
Figure 2. Co-occurrence network analysis of bacterial, fungal and protistan ASVs under different tillage practices and at different soil depths. (A) Multitrophic network including multiple ecological modules. The colors of the nodes represent different ecological modules; the percentages of bacterial, fungal and protistan ASVs (B); and the intraspecies and interspecies relationships (C) in each module.
ZP plots were constructed to identify the topological roles of each node in the co-occurrence network. A total of 30 microbial taxa (including 13 bacteria, 11 fungi and 6 protists) were detected as keystone species (Figure 3). Nineteen keystone taxa belonged to module 1, and modules 2 and 3 each contained 2 keystone taxa. The information for the selected keystone taxa is displayed in Supplementary Table S1. The bacterial keystone species mainly belonged to Proteobacteria (7 taxa), the fungal keystone species mainly belonged to Ascomycota (6 taxa), and the protist keystone species mainly belonged to Cercozoa (2 taxa).
Figure 3. ZP plot showing the distribution of ASVs based on their module-based topological roles. The topological role of each ASV was determined according to the scatter plot of within-module connectivity (Z) and among-module connectivity (P).
3.3 Relationships between microbial traits and SOC fractions
To determine the potential relationships between the microbial communities of specific modules and SOC fractions, we constructed correlations between the microbial community module connections and the SOC fractions. Figure 4 shows that there were significant correlations between microbial module communities and POC content (R2 = 0.74 for module 1 and R2 = 0.68 for module 2) and MAOC content (R2 = 0.51 for module 1 and R2 = 0.44 for module 2). However, no significant associations were detected between the SOC fractions and the other microbial module communities.
Figure 4. Links between the soil community of each module within the co-occurrence network and the SOC fraction content. The links between soil microbial community of module 1-4 with particulate organic carbon were sown in A-D; the links between soil microbial community of module 1-4 with mineral-associated organic carbon were sown in E-H.
Heatmaps revealed close associations between the richness of keystone taxa and SOC fractions (Figure 5). Overall, the bacterial richness demonstrated significant negative correlations with the POC and MAOC contents (except for BASV 147) (Figure 5A). There were significant positive correlations between FASV945 richness and SOC fraction contents (POC and MAOC), as well as between FASV945 richness and MAOC content. FASV95 richness was negatively correlated with MAOC content (Figure 5B). Moreover, the richness of PASV45 and PASV17 was positively correlated with the POC content, and the richness of PASV45 was also positively correlated with the MAOC content (Figure 5C).
Figure 5. Relationships of SOC fraction contents with the richness of bacterial (A), fungal (B), and protistan (C) keystone species.
In the present study, we selected keystone species that were significantly associated with the SOC fraction content for further analysis. The richness of keystone taxa was sensitive to the different tillage practices (Supplementary Figure S7). Compared with CT, CK increased the richness of bacterial keystone taxa by 21.51–520.75%. The richness of BASV8256 decreased by 32.15% under CK compared with that under CT (Supplementary Figure S7A). Compared with CK, CT increased the richness of FASV945 and FASV95 by 58.81 and 42.81%, respectively, and decreased the richness of FASV946 by 82.61%. In addition, compared to that under CK, PASV17 and PASV45 richness increased by 97.79 and 12.97%, respectively, under CT (Supplementary Figure S7B).
4 Discussion
CT has been considered a sustainable technique for properly managing soil and hence maintaining agroecosystem services by minimizing tillage operations to effectively avoid water infiltration and erosion (Müller et al., 2009). In this study, we compared crop yields, SOC fraction contents and microbial traits under different tillage practices and investigated the relationships between SOC fractions and microbial functions. The results of this research strengthen our understanding of SOC accumulation under different tillage practices on the Northeast China Plain.
4.1 Response of SOC fractions and maize yields to different tillage practices
The effect of tillage practices on crop yields has been studied frequently; however, no consistent conclusions have been drawn. It was reported that no- or minimum-tillage practices led to a 0–30% reduction in yields in Europe, which was affected by crop type, tillage technique, soil texture and crop rotation (Alaoui et al., 2020). Another study revealed that shallow tillage (8 cm strip depth) achieved the greatest yields (Licht and Al-Kaisi, 2005). This is mainly because shallow tillage can be a neutral solution to the problem of late seed emergence due to no-tillage by reducing soil disturbance (Araya et al., 2021). Notably, based on a 17-year experiment, CT practices increased maize yields by 12.2 to 20.1% (Ren et al., 2024), which was consistent with our results (Supplementary Figure S2). This can be explained partly by the fact that the CT method is generally implemented with straw residue left on topsoil while minimizing soil disturbance, which favors soil nutrient accumulation, moisture retention and temperature increase, resulting in faster seed emergence (Licht and Al-Kaisi, 2005). In summary, although optimizing tillage practices requires consideration of factors such as crop rotation and soil texture, CT increases crop yield in maize monoculture systems in the eolian sandy soil of the Northeast China Plain.
SOC is the key to soil fertility and is sensitive to changes in tillage practices (Figure 1). Lehmann and Kleber (2015) proposed a soil continuum model indicating that the input of exogenous organic materials (such as plant residues) is a prerequisite for SOC accumulation. Thus, the finding that CT can improve topsoil OC fractions is no surprise. However, the POC and MOC contents were not consistent. No-tillage practices improved the POC content in the 0–30 cm soil layer. Undecomposed and semidecomposed plant residues are the “core” of POC, which is encapsulated by macroaggregates (Samson et al., 2020). Therefore, the formation of macroaggregates and the accumulation of POC are generally complementary. The data on the distribution of aggregate sizes in the present study also validate this view (Supplementary Figure S3). An increased proportion of macroaggregates provides physical protection for POC and thus favors POC accumulation. CK practices disturb the physical structure at 0–25 cm soil depths through frequent plowing, thus destroying the formation of macroaggregate structures (Jat et al., 2019).
When the large biopolymers in the residues further decomposed, the small biopolymers and C monomers (such as root exudates and microbial necromass) can be adsorbed to the soil mineral surface and become MAOC (Lehmann and Kleber, 2015). In addition to straw return, microbial and maize biomass are important factors that cannot be ignored. As a supplementary exogenous C source in the soil, straw inevitably increases MAOC content after further degradation of residues (Vogel et al., 2014). Posteriorly, CT decreases the soil structure distribution and increases the soil density, which promotes the growth of crop roots to a certain extent (Shi et al., 2012). In addition, the crop yield data imply a well-developed root system that has secreted more organic matter (Van den Putte et al., 2010). With respect to microbial biomass, minimum tillage and residue retention increase the soil microbial population size due to the adequate energy supply and appropriate stoichiometry and result in MAOC sequestration in topsoil under CT practices (Ren et al., 2024). In summary, soil microbiomes play an irreplaceable role in straw degradation and SOC turnover. Therefore, revealing the microbial mechanisms responsible for SOC fraction accumulation under different tillage practices is crucial for improving soil fertility.
4.2 Microbial keystone species-driven SOC fraction sequestration by regulating specific module community functions
CT changed the soil microbial diversity and community, which subsequently affected SOC formation and accumulation. We found that soil depth affected the soil bacterial, fungal and protistan communities much more than did tillage practice (Supplementary Figure S4). This is partly because the soil layer has the greatest influence on the changes in nutrient accessibility and availability rather than tillage practices (Kong et al., 2011). As soil depth increases, nutrient pools decrease, and mineral protection increases, leading to increased difficulty in nutrient acquisition by soil microorganisms (Modak et al., 2019). As a result, oligotrophs may become the dominant species in the community.
Microbial communities were classified into different functional modules by identifying soil taxa strongly interacting with each other (Figure 2), which can indicate important ecological processes, different niches, and habitat preferences. Each module in a network is considered a functional unit that conducts an identifiable task (Chen et al., 2019). In the present study, strong relationships were observed between the SOC fractions and the microbial community in modules 1 and 2, which indicated the potential function of SOC turnover (Figure 4). A previous study demonstrated that microbial keystone species have great explanatory power in terms of network (module) structure and function (Delgado-Baquerizo et al., 2018). Thus, it is necessary to explore the function of keystone species within modules 1 and 2. However, more than two-thirds of the identified keystone species belong to module 1 communities (19/30), while only two keystone species (2/30) belong to module 2 communities (Supplementary Table S1). Therefore, it is more meaningful to study the function of keystone species-driven module communities on SOC turnover in module 1 community.
Microbial keystone taxa also exhibited significant correlations with SOC fractions (Figure 5). Among the bacteria, Rhizobiales_Incertae_Sedis (BASV7693) and Reyranellaceae (BASV11161) within module 1 were identified as Proteobacteria, which encompasses typical copiotrophs with low C use efficiency that leads to straw C loss (Dove et al., 2021). A previous study demonstrated that Gaiellaceae (BASV7998) richness is recognized as an indicator of the carbon-to-nitrogen ratio due to its ability to utilize labile C (Duan et al., 2021). Therefore, the increase in the above bacterial taxa indicated negative effects on straw-derived OC accumulation. Moreover, only the Sphingomonadaceae richness exhibited a positive association with the POC and MAOC contents. This was likely because Sphingomonadaceae can consume various C sources and become major exopolysaccharide contributors, which provide a source for MAOC formation (Lan et al., 2022). Compared with bacteria, fungi generally exhibit greater C use efficiency and straw decomposition ability (Clemmensen et al., 2015). The richness of Cephalotrichum (FASV945) and Herpotrichiellaceae (FASV946) in module 1 was conducive to POC and MAOC sequestration. As typical saprophytic fungi, Cephalotrichum and Herpotrichiellaceae are often considered to play roles in straw degradation, pathogen control, and crop growth promotion and are considered important indicators of soil health (Zhang et al., 2022). Zhang et al. (2023) indicated that Cephalotrichum was enriched after straw addition and promoted straw decomposition and SOC accumulation in saline-alkaline soils. Most Cephalotrichum species are known for their saprotrophic function in decomposing plant materials (Woudenberg et al., 2017), corresponding with our results that the Cephalotrichum played a dominant role in POC and MAOC sequestration. In addition, higher abundance of Cephalotrichum led to the higher fungal diversity (Wang et al., 2023). Jin et al. (2022) found that Cephalotrichum exerted significant inhibitory effects on several pathogenic bacteria. Due to these strong abilities, Cephalotrichum abundance was also considered as the biomarker of soil health. Therefore, we speculated that the increase in the abundances of Cephalotrichum can accelerate straw degradation and that the early and late products favor the formation of POC and MAOC, respectively.
In addition, protists can influence SOC accumulation through direct or indirect effects (Pellegrino et al., 2021). Our results revealed that Cercozoa (PASV45 in module 1) may be a possible participant that was significantly positively correlated with POC (p < 0.05) and MAOC (p < 0.001). Pellegrino et al. (2021) reported that CT practices increase the abundance of Cercozoa partly via several abiotic factors, such as soil moisture, clay content and N availability. As expected, Cercozoa are important consumers of straw residues, which promotes the fragmentation of straw to facilitate further decomposition and is a prerequisite for organic carbon accumulation (Gessner et al., 2010). These authors also indicated that Cercozoa was the keystone taxon in macroaggregates and was positively correlated with SOC by promoting residue decomposition (Delgado-Baquerizo et al., 2020). Moreover, Cercozoa was also thought to have been an important driving force in the formation of macroaggregates by reshaping the pore sizes in the soil (Berisso et al., 2012). Therefore, according to the theory of interaction between SOC and aggregate structure, Cercozoa-driven straw fragments can be encapsulated by aggregates and become the core of aggregate formation. The formation of aggregates provides physical protection for POC, which is conducive to POC accumulation. This conclusion is also consistent with previous reports that Cercozoa are crucial microorganisms in macroaggregate taking part to long-term sequestration and storage of SOC (Pellegrino et al., 2021). Another study indicated that Cercozoa species exhibited the highest numbers of links with bacteria and fungi through the construction of co-occurrence networks, which implied that they were potentially vital to soil food webs (microbiome predation) (Kou et al., 2020). Cercozoa are phagotrophs that may consume Acidobacteria, Proteobacteria, and Ascomycota and consequently increase microbe-derived C. Furthermore, Ma et al. (2024) confirmed that protozoa can regulate microbial carbon use efficiency and SOC formation by regulating fungal, bacterial and keystone module communities through structural equation model analysis. Among these factors, the Cercozoa-driven protozoan community was the most influential factor. Accordingly, Cercozoa mediated POC and MAOC accumulation, mainly through macroaggregate formation and microbial necromass supply.
After comparing the richness of the selected keystone species, our results showed that C-accumulating microbes were enriched under CT. Specifically, compared with CK, CT decreased the abundances of most keystone bacterial taxa and increased the abundances of specific fungal and protistan species in module 1 of the network, which promoted the sequestration of SOC fractions by straw degradation, aggregate formation and predation effects. As a consequence, soil bacteria, fungi and protistan taxa all participate in SOC turnover under different tillage practices, while the appointed keystone species-driven community of module 1 facilitated POC and MOC accumulation under CT.
5 Conclusion
In this study, we demonstrated the associations between microbial keystone taxa and SOC fractions under different tillage practices. Compared with CK, continuous 6-year CT significantly increased maize yields, aggregate stability, and POC (0–30 cm) and MAOC (0–20 cm) contents. Tillage practice and soil depth both influence bacterial, fungal and protistan communities, which might change the turnover of SOC fractions. The co-occurrence network indicated that the connectivity of module 1 was significantly related to POC and MAOC contents CT increased the richness of specific fungal (Cephalotrichum) and protistan (Cercozoa) species and thus promoted SOC fraction accumulation through straw degradation, macroaggregate formation and predation effects. The selected bacterial taxa were enriched in the CK treatment and resulted in SOC loss due to low C use efficiency. Taken together, our results revealed that stimulating the function of keystone taxa can drive the function of the module 1 community in SOC accumulation under CT practices, which is beneficial for maintaining soil fertility and productivity in eolian sandy soils on the Northeast China Plain.
Data availability statement
Raw sequencing data were deposited in the NCBI Sequence Read Archive under accession number PRJNA1111948. The data of SOC, soil nutrients content and maize yields can be found in the article/supplementary material.
Author contributions
Y-mL: Data curation, Funding acquisition, Investigation, Methodology, Writing – original draft. Y-mW: Writing – review & editing, Data curation, Investigation. G-wQ: Writing – review & editing. H-jY: Writing – review & editing. F-mL: Writing – review & editing. G-lW: Funding acquisition, Supervision, Writing – review & editing. YD: Funding acquisition, Investigation, Methodology, Supervision, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was funded by the National Key Research and Development Program of China (2022YFD1500704, 2022YFD1500305, 2016YFD0300806, 2023YFD1901005), the Grant of the President Foundation of Hefei Institutes of Physical Science of Chinese Academy of Sciences (YZJJ2023QN37).
Conflict of interest
The 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/fmicb.2024.1394179/full#supplementary-material
References
Achankeng, E., and Cornelis, W. (2023). Conservation tillage effects on European crop yields: a meta-analysis. Field Crop Res. 298:108967. doi: 10.1016/j.fcr.2023.108967
Alaoui, A., Barao, L., Ferreira, C. S. S., Schwilch, G., Basch, G., Garcia-Orenes, F., et al. (2020). Visual assessment of the impact of agricultural management practices on soil quality. Agron. J. 112, 2608–2623. doi: 10.1002/agj2.20216
Angers, D. A., and Eriksen-Hamel, N. S. (2008). Full-inversion tillage and organic carbon distribution in soil profiles: a meta-analysis. Soil Sci. Soc. Am. J. 72, 1370–1374. doi: 10.2136/sssaj2007.0342
Araya, T., Gebremedhin, A., Baudron, F., Hailemariam, M., Birhane, E., Nyssen, J., et al. (2021). Influence of 9 years of permanent raised beds and contour furrowing on soil health in conservation agriculture based systems in Tigray region, Ethiopia. Land Degrad. Dev. 32, 1525–1539. doi: 10.1002/ldr.3816
Banerjee, S., Kirkby, C. A., Schmutter, D., Bissett, A., Kirkegaard, J. A., and Richardson, A. E. (2016). Network analysis reveals functional redundancy and keystone taxa amongst bacterial and fungal communities during organic matter decomposition in an arable soil. Soil Biol. Biochem. 97, 188–198. doi: 10.1016/j.soilbio.2016.03.017
Bastian, F., Bouziri, L., Nicolardot, B., and Ranjard, L. (2009). Impact of wheat straw decomposition on successional patterns of soil microbial community structure. Soil Biol. Biochem. 41, 262–275. doi: 10.1016/j.soilbio.2008.10.024
Berisso, F. E., Schjonning, P., Keller, T., Lamandé, M., Etana, A., de Jonge, L. W., et al. (2012). Persistent effects of subsoil compaction on pore size distribution and gas transport in a loamy soil. Soil Till. Res. 122, 42–51. doi: 10.1016/j.still.2012.02.005
Biddle, J. F., Fitz-Gibbon, S., Schuster, S. C., Brenchley, J. E., and House, C. H. (2008). Metagenomic signatures of the Peru margin subseafloor biosphere show a genetically distinct environment. Proc. Natl. Acad. Sci. U.S.A. 105, 10583–10588. doi: 10.1073/pnas.0709942105
Chen, F., Yan, G. Y., Xing, Y. J., Zhang, J. H., Wang, Q. G., Wang, H. L., et al. (2019). Effects of N addition and precipitation reduction on soil respiration and its components in a temperate forest. Agric. For. Meteorol. 271, 336–345. doi: 10.1016/j.agrformet.2019.03.021
Clemmensen, K. E., Finlay, R. D., Dahlberg, A., Stenlid, J., Wardle, D. A., and Lindahl, B. D. (2015). Carbon sequestration is related to mycorrhizal fungal community shifts during long-term succession in boreal forests. New Phytol. 205, 1525–1536. doi: 10.1111/nph.13208
Cotrufo, M. F., Wallenstein, M. D., Boot, C. M., Denef, K., and Paul, E. (2013). The microbial efficiency-matrix stabilization (MEMS) framework integrates plant litter decomposition with soil organic matter stabilization: do labile plant inputs form stable soil organic matter? Glob. Chang. Biol. 19, 988–995. doi: 10.1111/gcb.12113
Delgado-Baquerizo, M., Oliverio, A. M., Brewer, T. E., Benavent-González, A., Eldridge, D. J., Bardgett, R. D., et al. (2018). A global atlas of the dominant bacteria found in soil. Science 359:320-+. doi: 10.1126/science.aap9516
Delgado-Baquerizo, M., Reich, P. B., Trivedi, C., Eldridge, D. J., Abades, S., Alfaro, F. D., et al. (2020). Multiple elements of soil biodiversity drive ecosystem functions across biomes. Nat. Ecol. Evol. 4, 210–220. doi: 10.1038/s41559-019-1084-y
Deng, Y., Jiang, Y. H., Yang, Y. F., He, Z. L., Luo, F., and Zhou, J. Z. (2012). Molecular ecological network analyses. BMC Bioinformatics 13:113. doi: 10.1186/1471-2105-13-113
Dove, N. C., Torn, M. S., Hart, S. C., and Tas, N. (2021). Metabolic capabilities mute positive response to direct and indirect impacts of warming throughout the soil profile. Nat. Commun. 12:2089. doi: 10.1038/s41467-021-22408-5
Duan, Y., Chen, L., Li, Y. M., Li, J. Y., Zhang, C. Z., Ma, D. H., et al. (2023). Nitrogen input level modulates straw-derived organic carbon physical fractions accumulation by stimulating specific fungal groups during decomposition. Soil Tillage Res. 225:105560. doi: 10.1016/j.still.2022.105560
Duan, Y., Chen, L., Li, Y. M., Wang, Q. Y., Zhang, C. Z., Ma, D. H., et al. (2021). N, P and straw return influence the accrual of organic carbon fractions and microbial traits in a Mollisol. Geoderma 403:115373. doi: 10.1016/j.geoderma.2021.115373
Gao, Z. L., Karlsson, I., Geisen, S., Kowalchuk, G., and Jousset, A. (2019). Protists: puppet masters of the rhizosphere microbiome. Trends Plant Sci. 24, 165–176. doi: 10.1016/j.tplants.2018.10.011
Geisen, S., and Bonkowski, M. (2018). Methodological advances to study the diversity of soil protists and their functioning in soil food webs. Appl. Soil Ecol. 123, 328–333. doi: 10.1016/j.apsoil.2017.05.021
Gessner, M. O., Swan, C. M., Dang, C. K., Mckie, B. G., Bardgett, R. D., Wall, D. H., et al. (2010). Diversity meets decomposition. Trends Ecol. Evol. 25, 372–380. doi: 10.1016/j.tree.2010.01.010
Ghannoum, M. A., Jurevic, R. J., Mukherjee, P. K., Cui, F., Sikaroodi, M., Naqvi, A., et al. (2010). Characterization of the Oral fungal microbiome (Mycobiome) in healthy individuals. PLoS Pathog. 6:e1000713. doi: 10.1371/journal.ppat.1000713
Herath, H. M. S. K., Camps-Arbestain, M., Hedley, M., Van Hale, R., and Kaal, J. (2014). Fate of biochar in chemically- and physically-defined soil organic carbon pools. Org. Geochem. 73, 35–46. doi: 10.1016/j.orggeochem.2014.05.001
Hewins, D. B., Sinsabaugh, R. L., Archer, S. R., and Throop, H. L. (2017). Soil-litter mixing and microbial activity mediate decomposition and soil aggregate formation in a sandy shrub-invaded Chihuahuan Desert grassland. Plant Ecol. 218, 459–474. doi: 10.1007/s11258-017-0703-4
Huang, R., Lan, M. L., Liu, J., and Gao, M. (2017). Soil aggregate and organic carbon distribution at dry land soil and paddy soil: the role of different straws returning. Environ. Sci. Pollut. R. 24, 27942–27952. doi: 10.1007/s11356-017-0372-9
Jat, S. L., Parihar, C. M., Singh, A. K., Nayak, H. S., Meena, B. R., Kumar, B., et al. (2019). Differential response from nitrogen sources with and without residue management under conservation agriculture on crop yields, water-use and economics in maize-based rotations. Field Crop Res. 236, 96–110. doi: 10.1016/j.fcr.2019.03.017
Jin, Q., Zhang, Y., Ma, Y., Sun, H., Guan, Y., Liu, Z., et al. (2022). The composition and function of the soil microbial community and its driving factors before and after cultivation of Panax ginseng in farmland of different ages. Ecol. Indic. 145:109748. doi: 10.1016/j.ecolind.2022.109748
Kong, A. Y. Y., Scow, K. M., Córdova-Kreylos, A. L., Holmes, W. E., and Six, J. (2011). Microbial community composition and carbon cycling within soil microenvironments of conventional, low-input, and organic cropping systems. Soil Biol. Biochem. 43, 20–30. doi: 10.1016/j.soilbio.2010.09.005
Lan, J. C., Wang, S. S., Wang, J. X., Qi, X., Long, Q. X., and Huang, M. Z. (2022). The shift of soil bacterial community after afforestation influence soil organic carbon and aggregate stability in karst region. Front. Microbiol. 13:901126. doi: 10.3389/fmicb.2022.901126
Lehmann, J., and Kleber, M. (2015). The contentious nature of soil organic matter. Nature 528, 60–68. doi: 10.1038/nature16069
Li, J. W., Li, M. Y., Dong, L. B., Wang, K. B., Liu, Y. L., Hai, X. Y., et al. (2020). Plant productivity and microbial composition drive soil carbon and nitrogen sequestrations following cropland abandonment. Sci. Total Environ. 744:140802. doi: 10.1016/j.scitotenv.2020.140802
Li, Y., Zhang, Q. P., Cai, Y. J., Yang, Q., and Chang, S. X. (2020). Minimum tillage and residue retention increase soil microbial population size and diversity: implications for conservation tillage. Sci. Total Environ. 716:137164. doi: 10.1016/j.scitotenv.2020.137164
Licht, M. A., and Al-Kaisi, M. (2005). Strip-tillage effect on seedbed soil temperature and other soil physical properties. Soil Tillage Res. 80, 233–249. doi: 10.1016/j.still.2004.03.017
Lu, R. K. (2000). The analysis method of soil agricultural chemistry. Beijing: Chinese Agricultural Sciences and Technology Press (in Chinese).
Ma, L., Zhou, G., Zhang, J., Jia, Z., Zou, H., Chen, L., et al. (2024). Long-term conservation tillage enhances microbial carbon use efficiency by altering multitrophic interactions in soil. Sci. Total Environ. 915:170018. doi: 10.1016/j.scitotenv.2024.170018
Modak, K., Ghosh, A., Bhattacharyya, R., Biswas, D. R., Das, T. K., Das, S., et al. (2019). Response of oxidative stability of aggregate-associated soil organic carbon and deep soil carbon sequestration to zero-tillage in subtropical India. Soil Till. Res. 195:104370. doi: 10.1016/j.still.2019.104370
Müller, E., Wildhagen, H., Quintern, M., Hess, J., Wichern, F., and Joergensen, R. G. (2009). Spatial patterns of soil biological and physical properties in a ridge tilled and a ploughed Luvisol. Soil Till. Res. 105, 88–95. doi: 10.1016/j.still.2009.05.011
Pearsons, K. A., Omondi, E. C., Zinati, G., Smith, A., and Rui, Y. C. (2023). A tale of two systems: does reducing tillage affect soil health differently in long-term, side-by-side conventional and organic agricultural systems? Soil Till. Res. 226:105562. doi: 10.1016/j.still.2022.105562
Pellegrino, E., Piazza, G., Helgason, T., and Ercoli, L. (2021). Eukaryotes in soil aggregates across conservation managements: major roles of protists, fungi and taxa linkages in soil structuring and C stock. Soil Biol. Biochem. 163:108463. doi: 10.1016/j.soilbio.2021.108463
Raus, L., Jitareanu, G., Ailincai, C., Pârvan, L., and Topa, D. (2016). Impact of different soil tillage systems and Organo-mineral fertilization on physical properties of the soil and on crops yield in Pedoclimatical conditions of Moldavian plateau. Rom. Agric. Res. 33, 111–123,
Ren, Z. J., Han, X. J., Feng, H. X., Wang, L. F., Ma, G., Li, J. H., et al. (2024). Long-term conservation tillage improves soil stoichiometry balance and crop productivity based on a 17-year experiment in a semi-arid area of northern China. Sci. Total Environ. 908:168283. doi: 10.1016/j.scitotenv.2023.168283
Samson, M. E., Chantigny, M. H., Vanasse, A., Menasseri-Aubry, S., Royer, I., and Angers, D. A. (2020). Management practices differently affect particulate and mineral-associated organic matter and their precursors in arable soils. Soil Biol. Biochem. 148:107867. doi: 10.1016/j.soilbio.2020.107867
Sarker, J. R., Singh, B. P., Cowie, A. L., Fang, Y. Y., Collins, D., Dougherty, W. J., et al. (2018). Carbon and nutrient mineralisation dynamics in aggregate-size classes from different tillage systems after input of canola and wheat residues. Soil Biol. Biochem. 116, 22–38. doi: 10.1016/j.soilbio.2017.09.030
Shi, X. H., Yang, X. M., Drury, C. F., Reynolds, W. D., McLaughlin, N. B., and Zhang, X. P. (2012). Impact of ridge tillage on soil organic carbon and selected physical properties of a clay loam in southwestern Ontario. Soil Till. Res. 120, 1–7. doi: 10.1016/j.still.2012.01.003
Six, J., Bossuyt, H., Degryze, S., and Denef, K. (2004). A history of research on the link between (micro)aggregates, soil biota, and soil organic matter dynamics. Soil Till. Res. 79, 7–31. doi: 10.1016/j.still.2004.03.008
Six, J., Conant, R. T., Paul, E. A., and Paustian, K. (2002). Stabilization mechanisms of soil organic matter: implications for C-saturation of soils. Plant Soil 241, 155–176. doi: 10.1023/A:1016125726789
Stoeck, T., Bass, D., Nebel, M., Christen, R., Jones, M. D. M., Breiner, H. W., et al. (2010). Multiple marker parallel tag environmental DNA sequencing reveals a highly complex eukaryotic community in marine anoxic water. Mol. Ecol. 19, 21–31. doi: 10.1111/j.1365-294X.2009.04480.x
Tang, Q., Ti, C. P., Xia, L. L., Xia, Y. Q., Wei, Z. J., and Yan, X. Y. (2019). Ecosystem services of partial organic substitution for chemical fertilizer in a peri-urban zone in China. J. Clean. Prod. 224, 779–788. doi: 10.1016/j.jclepro.2019.03.201
Topa, D., Cara, I. G., and Jitareanu, G. (2021). Long term impact of different tillage systems on carbon pools and stocks, soil bulk density, aggregation and nutrients: a field meta-analysis. Catena 199:105102. doi: 10.1016/j.catena.2020.105102
Van den Putte, A., Govers, G., Diels, J., Gillijns, K., and Demuzere, M. (2010). Assessing the effect of soil tillage on crop growth: a meta-regression analysis on European crop yields under conservation agriculture. Eur. J. Agron. 33, 231–241. doi: 10.1016/j.eja.2010.05.008
Vogel, C., Mueller, C. W., Höschen, C., Buegger, F., Heister, K., Schulz, S., et al. (2014). Submicron structures provide preferential spots for carbon and nitrogen sequestration in soils. Nat. Commun. 5:2947. doi: 10.1038/ncomms3947
Wang, J., Li, H., Cheng, Z., Yin, F., Yang, L., and Wang, Z. (2023). Changes in soil bacterial and fungal community characteristics in response to long-term mulched drip irrigation in oasis agroecosystems. Agric. Water Manage. 279:108178. doi: 10.1016/j.agwat.2023.108178
Wang, C., Qu, L. R., Yang, L. M., Liu, D. W., Morrissey, E., Miao, R. H., et al. (2021). Large-scale importance of microbial carbon use efficiency and necromass to soil organic carbon. Glob. Chang. Biol. 27, 2039–2048. doi: 10.1111/gcb.15550
Wang, S. C., Zhao, Y. W., Wang, J. Z., Zhu, P., Cui, X., Han, X. Z., et al. (2018). The efficiency of long-term straw return to sequester organic carbon in Northeast China's cropland. J. Integr. Agric. 17, 436–448. doi: 10.1016/S2095-3119(17)61739-8
Woudenberg, J., Sandoval-Denis, M., Houbraken, J., Seifert, K. A., and Samson, R. A. (2017). Cephalotrichum and related synnematous fungi with notes on species from the built environment. Stud. Mycol. 88, 137–159. doi: 10.1016/j.simyco.2017.09.001
Yu, Z. H., Zhang, J. B., Zhang, C. Z., Xin, X. L., and Li, H. (2017). The coupling effects of soil organic matter and particle interaction forces on soil aggregate stability. Soil Till. Res. 174, 251–260. doi: 10.1016/j.still.2017.08.004
Zhang, Y., Li, X. J., Gregorich, E. G., McLaughlin, N. B., Zhang, X. P., Guo, Y. F., et al. (2019). Evaluating storage and pool size of soil organic carbon in degraded soils: tillage effects when crop residue is returned. Soil Till. Res. 192, 215–221. doi: 10.1016/j.still.2019.05.013
Zhang, S. X., Li, Q., Lü, Y., Zhang, X. P., and Liang, W. J. (2013). Contributions of soil biota to C sequestration varied with aggregate fractions under different tillage systems. Soil Biol. Biochem. 62, 147–156. doi: 10.1016/j.soilbio.2013.03.023
Zhang, X. P., Li, Q. L., Zhong, Z. K., Huang, Z. Y., Bian, F. Y., Yang, C. B., et al. (2022). Changes in soil organic carbon fractions and fungal communities, subsequent to different management practices in Moso bamboo plantations. J. Fungi (Basel) 8:8. doi: 10.3390/jof8060640
Zhang, L., Tang, C., Yang, J., Yao, R., Wang, X., Xie, W., et al. (2023). Salinity-dependent potential soil fungal decomposers under straw amendment. Sci. Total Environ. 891:164569. doi: 10.1016/j.scitotenv.2023.164569
Zhao, R. D., He, M., Jiang, C. L., Li, C. L., and Liu, F. (2021). Microbial community structure in rhizosphere soil rather than that in bulk soil characterizes aggregate-associated organic carbon under long-term forest conversion in subtropical region. Rhizosphere 20:100438. doi: 10.1016/j.rhisph.2021.100438
Keywords: conservation tillage, particulate organic carbon, mineral-associated organic carbon, microbial community, soil depth
Citation: Li Y-m, Wang Y-m, Qiu G-w, Yu H-j, Liu F-m, Wang G-l and Duan Y (2024) Conservation tillage facilitates the accumulation of soil organic carbon fractions by affecting the microbial community in an eolian sandy soil. Front. Microbiol. 15:1394179. doi: 10.3389/fmicb.2024.1394179
Edited by:
Peng Shi, Xi'an University of Technology, ChinaReviewed by:
Shankar Ganapathi Shanmugam, Mississippi State University, United StatesBaoku Shi, Northeast Normal University, China
Copyright © 2024 Li, Wang, Qiu, Yu, Liu, Wang and Duan. 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: Gen-lin Wang, wanggenlin2005@163.com; Yan Duan, duanyan@iim.ac.cn