- 1Ministry of Agriculture Key Laboratory of Plant Nutrition and Fertilizer, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing, China
- 2College of Forestry, Hebei Agricultural University, Baoding, China
- 3Institute of Plant Protection and Soil Fertility, Hubei Academy of Agricultural Sciences, Wuhan, China
Although soil microorganism is an active area of research, we are still in the early stages of understanding how living microorganisms influence the accumulations of soil microbial residues under different agricultural practices. Based on a 39-year fertilization experiment, we characterized the soil microbiota and correlated their compositions to soil microbial residues, which are indicated by amino sugars under a rice-wheat rotation. In the present study, fertilization regimes and crop season all exerted significant impacts on the compositions of soil microbial communities and their residues, although no significant difference in the microbial residues was found between soil depth (0–10 cm vs. 10–20 cm). Compared within fertilization regimes, the long-term fertilization, especially the application of organic manure, stimulated the accumulations of carbon (C) and nitrogen in soils and microbial residues. Upland soils in wheat season accumulated more microbial residues, particularly in fungal residues, than paddy soils in rice season. Our results suggested that the long-term application of organic manure favored the growth of soil microbial communities, and then increased the contents of microbial residues, particularly in fungal residues, leading to an enlargement of soil C pools. The keystone taxa Pseudaleuria identified by network analysis showed a significantly positive potential in soil C sequestration by increasing the accumulation of fungal residues. Thus, this study revealed the strong and close connections between microbial communities and their residues, and provided evidence about the critical role of keystone taxa in regulating C sequestration.
Introduction
Soil organic carbon (SOC), retained in the largest pool of terrestrial carbon, is critical in improving soil fertility, maintaining ecosystems sustainability, and mitigating climate change (Schimel and Schaeffer, 2012; Keenan and Williams, 2018). Recent developments have raised the attention in exploring the processes of SOC transformation and sequestration (Miltner et al., 2012; Liang et al., 2019), which can be influenced by various factors, e.g., environmental conditions, edaphic factors, microorganism, and human interference (Schimel and Schaeffer, 2012; Wiesmeier et al., 2019). Empirical studies have suggested the important role of microorganisms in SOC sequestration through microbial metabolism and their residues (Naylor et al., 2020). For instance, living soil microbes usually participate in SOC sequestration by converting plant detritus, root exudates, and other nutrients into microbial cellular components and byproducts, while dead microbes (microbial residues) generally have organic-mineral bonds that stabilize their C content (Liang et al., 2017). Current estimates suggested that soil microbial residues contribute approximately 30–60% to SOC, fungal residues C (FRC) contribute approximately 20–50%, and bacterial residues C (BRC) contribute approximately 10–30% (Liang et al., 2019). Hence, it is necessary to clarify the changes in microbial residues to further evaluate the soil C stocks.
A number of studies have focused on the mechanism of microbial residues affected by abiotic factors, e.g., temperature (Chen G. P. et al., 2020), precipitation (Fierer, 2017), soil pH (Hu et al., 2018), nutrient availability (He et al., 2011), and redox status (Xia et al., 2019), but the biotic factors in regulating the accumulations of microbial residues are still uncertain. Previous studies have reported that the microbial cell walls with a greater proportion of recalcitrant components, such as peptidoglycan in gram-positive bacteria and chitin in fungal cell walls, were considered to be more resistant to decomposition (Kögel-Knabner, 2002). Bacteria and fungi, as the widely spread microbiota, shared a microhabitat and shaped the biological communities above and below the Earth’s surface (Nazir et al., 2009; Frey-Klett et al., 2011), leading to a further change in organic C accumulation and stability. Although there is no strong difference between the turnover of bacterial and fungal residues (Joergensen, 2018), the interactions between living microorganism and their residues remain unclear and elusive, which needs further evidence.
The agroecosystem, which contained 16.32 ± 0.41 Pg C in croplands and accounted for 20.6% of the total terrestrial C pool in China (Tang et al., 2018), could regulate the global C balance through different agricultural practices, e.g., crop rotation (Schimel and Schaeffer, 2012), fertilization regimes (Ye et al., 2019), as well as soil depth (Roth et al., 2011). Soils in crop rotation, especially under the rice-wheat rotation, are typified by high heterogeneity in the distribution of soil nutrients and ecological niches, which can promote the changes in the microbial community biomass and compositions, as well as microbial functions in C cycling (Pan et al., 2003; Schimel and Schaeffer, 2012). A large-scale research on the distribution of microbial residues in paddy and upland soils throughout eastern China suggested that paddy soils were enriched with a greater proportion of bacterial residues, whereas upland soils were replenished by fungal residues (Chen et al., 2021). The long-term inorganic fertilization usually decreased the relative abundances of functional genes that participated in soil C and N cycling, while organic manure increased their abundances (Jiao et al., 2019) and the accumulations of microbial residues (Ye et al., 2019). Current studies also suggested a decrease in soil C accumulations in microbial residues with the increase in soil depth (Peixoto et al., 2020; Gies et al., 2021). Thus, based on the critical role of living and dead microorganisms in soil C cycling, understanding how agricultural practices impact the interactions between living microbiota and microbial residues is an important subject toward a more sustainable agroecosystem.
Amino sugars that can be hydrolyzed from microbial cell walls have been extensively used as proxies for soil microbial residues (Liang et al., 2019). Muramic acid (MurN) occurs exclusively in peptidoglycans in bacterial cell walls (Vollmer et al., 2008), glucosamine (GlcN) primarily originates from fungal chitin (Parsons, 1981), while galactosamine (GluN) originates from both bacterial and fungal cell wall (Engelking et al., 2007). Hence, we quantified the bacterial and fungal residues in soils based on the contents of MurN and GlcN, respectively (Appuhn and Joergensen, 2006; Joergensen, 2018). Moreover, the Python language package Networkx provided a useful tool for exploring and analyzing the complicated connections between soil variables and microbial communities.1 Therefore, we carried out this study on the compositions of bacterial and fungal communities and the content of C and N in the soil and microbial residues based on a 39-year fertilization experiment with a rice-wheat crop rotation. This study aimed to (1) explore how the long-term fertilization regimes, crop season, and soil depth affected the compositions of living and dead microorganisms and (2) evaluate the influence of fertilization regimes, crop season, and soil depth on soil C pools by mediating soil living and dead microorganisms. Considering the different fertilization regimes and heterogeneous microenvironment in the rice-wheat rotation, we hypothesized that (i) fertilization regimes, crop season, and soil depth exerted impacts on the soil microbial community compositions, and then altered the contents and components of microbial residues and (ii) keystone taxa influencing soil microbial residues also affected the size of soil C pools.
Materials and Methods
Study Site and Soil Sampling
A long-term field experiment, i.e., a rice-wheat cropping system, was established at Ezhou (114.7°E, 30.4°N), Hubei Province, China, in 1981. Wheat is sown in late November and harvested in late May, and rice is sown in early June and harvested in early October. According to the National Data Center for Meteorological Sciences,2 from 2019 to 2020, the average precipitation was 3.2 and 2.7 mm, and the average temperature was 27.3 and 12.5°C in rice and wheat season, respectively, per day. The soil is derived from yellow-brown paddy soil, classified as Cambisols (Nachtergaele, 2017). The soil before cultivation in 1981 had the following properties: pH 6.3, SOC 15.9 g kg–1, TN 1.8 g kg–1, available phosphate (AP) 5.0 mg kg–1, available potassium (AK), 98.5 mg kg–1. The experiment had a complete randomized design with different fertilization regimes and three replicates. Each plot measured 5 m × 8 m. The fertilization regimes in this study contained: (1) no fertilizer, NoF, (2) inorganic nitrogen, phosphate and potassium fertilizers, NPK, (3) manure, M, and (4) the application of NPK combined with manure, NPKM. We applied urea, superphosphate, and potassium chloride as N, P, and K fertilizers, respectively. The detailed application rates of fertilizers are displayed in Table 1. The organic manure was mainly composed of swine manure (69% moisture) and had the following properties: 262.2 g kg–1 SOC, 15.1 g kg–1 TN, 20.8 g kg–1 phosphorous pentoxide, and 13.6 g kg–1 potassium dioxide.
Soil Sampling and Analysis
Five intact soil cores were collected in October 2019 and May 2020, using stainless steel cores (5 cm inner diameter × 10 cm height) from a depth of 0–10 and 10–20 cm in each plot, which finally contained 48 samples (4 fertilization regimes × 2 crop season × 2 soil depth × 3 replicates). The fresh soil was transported to the laboratory immediately, sieved through a 2 mm mesh, and divided into three subsamples. One subsample was prepared for microbial DNA extraction, one was stored at 4°C for the analysis of extracellular enzyme activity, and the third was air-dried to analyze the soil physicochemical properties. SOC and TN were determined by an Elemental Analyzer (VARIO MACRO cube, Germany).
Soil Amino Sugar Analysis
The extraction of soil amino sugars, including GlcN, GalN, and MurN, was determined following Zhang and Amelung (1996). Briefly, a weighted air-dried soil sample (contained about 0.3 mg N) was hydrolyzed with 6 mol l–1 HCl at 105°C for 8 h, and then cooled to room temperature. The solution was mixed with 100 μl inositol and dried by a rotatory evaporator. The precipitate was dissolved in distilled water, adjusted to pH 6–6.8, centrifuged at 3,000 r min–1, and freeze-dried. Then, methanol was added to freeze-dried supernatants, centrifuged at 3,000 r min–1, transferred to vials, and dried by N2 at 45°C. Next, the residue was dissolved in 1 ml distilled water and 100 μl N-methylglucamine, mixed, and freeze-dried. The purified amino sugars were reacted with ethyl acetate and n-hexane. The amino sugar derivatives were separated on a gas chromatograph equipped with an HP-5 column (30 m × 0.25 mm × 0.25 mm) and quantified using a flame ionization detector (GC-7890B, Agilent, United States). The temperature program was set as follows: 120°C for 1 min, increased by 10°C min–1 to 250°C and held for 2.5 min, increased by 20°C min–1 to 270°C and held for 2 min. The injector temperature was 250°C and the split ratio was 10: 1. The peaks were identified by comparing sample retention times with those of pure amino sugar standards using the Chemstation software. Myo-inositol was added as the first standard after hydrolysis and N-methylglucamine was added as the second standard before the derivatization to assess the recovery of amino sugars.
The contents of amino sugars, microbial residue-C and N were calculated using the analytical equations developed by Appuhn and Joergensen (2006) and Joergensen (2018).
where mi represents the quality of inositol, and Ai and Ax present the peak area of inositol and amino sugars, respectively. In addition, 45 is the conversion coefficient of MurN to BRC, 6.67 is the conversion coefficient of MurN to bacterial residues-N (BRN), 179.17 is the molecular weight of GlcN, 9 is the conversion coefficient of fungal GlcN to FRC, and 1.4 is the conversion coefficient of GlcN to fungal residues-N (FRN).
Amplicon Sequencing and Bioinformatics
DNA was extracted from 500 mg soil using the FASTDNA Spin Kit (MP Biomedical, Santa Ana, CA, United States) according to the manufacturer’s protocols. The DNA concentration and quality were checked by NanoDrop (Thermo NanoDrop, 2000, United States). For bacteria, the V3-V4 hypervariable regions of the 16S rRNA gene were targeted using primers 515F (GTGCCAGCMGCCGCGGTAA) and 806R (GGACTACVSGGGTATCTAAT); for fungi, the ITS1 region was amplified using ITS1F (CTTGGTCATTTAGAGGAAGTAA) and ITS2R (GCTGCGTTCTTCATCGATGC). Sequencing was performed on the Illumina MiSeq platform. Sequenced paired-end reads were joined using FLASH3 and were filtered according to quality using Fastp tools.4 Briefly, sequences were joined (overlapping pairs) and grouped by samples following the barcodes and then the barcodes were removed. Then, sequences <50 bp or with ambiguous base calls were removed. Reads were denoised by unoise3, and any chimeric sequences were removed using the USEARCH tool based on silva_16s_v123. Sequences were then split into operational taxonomic units (OTUs) at 97% similarity using the UPARSE pipeline in USEARCH. Plastid and non-bacteria were removed by USEARCH based on silva_16s_v123. This involved quality filtering and 97% clustering of the ITS1 region as indicated above for the 16S processing, using the UNITE database (utax_reference_dataset_all_02.02.2019), for chimera removal, taxonomic identification, and non-fungal removal of representative OTUs. The bacterial and fungal OTU abundance tables were rarefied to an even depth of 30,000 reads for 16S and 25,000 for ITS sequences. The raw sequencing data were deposited in the NCBI Sequence Read Archive (SRA)5 with accession No. PRJNA796782.
Statistical Analysis
The data were processed using Microsoft Excel 2019 and the analysis of variance (ANOVA) was conducted by using SPSS 18.0 (IBM, United States), followed by a Fisher’s least significant difference (LSD) test at p < 0.05 to determine the effects of fertilization regimes, crop season, and soil depth on the contents of C and N in soil and microbial residues. PERMANOVA was conducted using the adonis () function with 9,999 permutations, variance portioning analysis (VPA) was determined using the varpart () function, and differences within-group centroids were calculated using the betadisper () function. All these functions were from the ‘‘vegan’’ package6 in R (version 3.6.1). Principal co-ordinates analysis (PCoA) was conducted to calculate the differences within microbial community compositions using the BetaDiv () function in “phyloseq” package, which was based on the Bray-Curits distance (Adonis test) with 999 permutations.
Partial least squares path model (PLS-PM) was performed using the “plspm” package (Sanchez et al., 2015) to further investigate the potential direct and indirect effects of inorganic and organic fertilization, crop season, and soil depth on soil C and N contents by influencing the compositions of microbial communities and their residues, which finally influence soil C pools. Bacterial community was represented by the sample values of PC1 and PC2 derived from PCoA analysis. Bacterial residues contained BRC and BRN, while fungal residues contained FRC and FRN. The quality of the PLS-PM was evaluated by examining the goodness of fit (GoF) index, in which >0.7 means an acceptable overall prediction performance of the model, and by calculating the coefficients of determination of the variables (R2). Then, we estimated the standardized total effects of all variables upon soil SOC and TN by using the “pls” package.
To explore the complex relationships between microbial communities and their residues, we further constructed the networks using the Python language package Networkx (see text footnote 1). Only OTUs with top 500 abundances were considered. Edges were represented as links between microbial OTUs and the soil variables, which were tested by the Pearson correlation analysis. Only edges with p < 0.05 and |ρ| > 0.5 were included in the final networks. The network properties, including node degree, degree centrality, closeness centrality, and betweenness centrality, were also calculated. A node with the highest score of closeness centrality, betweenness centrality, and degree centrality are indicated as keystone taxa (Chen et al., 2020). And the networks were visualized using the Matplotlib Python package.7
Some figures were plotted using OriginPro 9 (OriginLab, United States) and modified using Adobe Illustrator CS6 (Adobe, United States).
Results
Changes in the Contents of Soil C and N
Fertilization regimes and crop season exerted significant impacts on the contents of SOC and TN (p < 0.05), while soil depth did not exert any impact (Table 2, p > 0.05). For fertilization regimes, compared with the original soil, the long-term fertilization enlarged SOC and TN pools, with an increase of 3.8–31.4% and 1.5–19.2%, respectively, while the long-term non-fertilization significantly decreased, with an average decrease of 15.8 and 6.6%, respectively. Compared with soils in NPK, the long-term application of organic manure, especially NPKM, significantly increased SOC and TN contents by 42.2 and 24.6%, respectively. Compared with soils in rice season, a relatively higher TN content was observed in wheat fields, which increased by 12.7%.
For the contents of soil microbial residues, fertilization regimes and crop season provided significant contributions to their variations, especially crop rotation. Meanwhile, soil depth only exerted significant impacts on the contents of FRC and FRN. In detail, when compared with original soil, the long-term fertilization significantly increased the contents of BRC and BRN by 29.0–220.5% and 28.9–221.0%, respectively, while NPKM significantly increased the FRC and FRN by an average of 13.4 and 13.5%, respectively. Compared with NoF, NPKM significantly increased the accumulations of BRC, BRN, FRC, and FRN by 12.7–52.6%, 12.8–53.3%, 13.7–59.2%, and 13.6–59.3%, respectively. For crop season, the accumulations of microbial residues were relatively higher in the wheat fields, and the contents of BRC, BRN, FRC, and FRN were on average 24.2, 24.4, 33.7, and 33.8% higher than those in the rice fields. Compared with the upper soil layer, the contents of FRC and FRN were significantly higher in the lower layer with an increase of 33.7 and 33.8%, respectively.
Moreover, the ratios of microbial residues to SOC and TN can reflect the contributions of the bacterial or fungal residues to soil C and N pools. The ratios of FRC/SOC and FRN/TN accounted for 23.7–51.1% and 30.0–55.0%, respectively, which were significantly higher than the ratios of BRC/SOC (9.1–15.7%) and BRN/TN (10.9–19.9%).
Changes in the Living Soil Microbial Communities
Differences in the overall microbial community compositions based on the PCoA revealed that both bacterial and fungal community compositions demonstrated a significant structuring in response to fertilization regimes, crop season, and soil depth (Figure 1A and Table 3). In detail, soil depth made the highest contribution to the variation of bacterial community compositions, and organic manure contributed more to the variation of fungi (Figure 1B). Major bacterial and fungal phyla were observed in all treatments (Figure 1C). Within fertilization regimes, compared with NoF, NPKM significantly increased the relative abundances of Planctomycetota, Verrucomicrobia, and Ascomycota, with an increase of 22.6, 15.1, and 13.3%, respectively. On the contrary, NPKM significantly increased the relative abundances of Chloroflexi, Gemmatimonadetes, Basidiomycota, Glomeromycota, and Chytridiomycota, with a decrease of 24.2, 13.1, 108.9, and 10.7%, respectively. For crop season, only the relative abundance of Bacteroidetes was 20.7% higher in the wheat fields than that in the rice fields. Compared with the lower soil layer, the relative abundances of Actinobacteria, Firmicutes, Alphaproteobacteria, Gammaproteobacteria, Rozellomycota, Mortierellomycota, and Chytridiomycota were significantly higher in the upper layer, while the relative abundances of Chloroflexi, Bacteroidetes, Nitrospirae, and Glomeromycota were significantly lower in the upper layer.
Figure 1. Soil bacterial and fungal community compositions. (A) PCoA analysis of soil bacterial (upper) and fungal (lower) community compositions. (B) Variance partitioning analysis (VPA) of the bacterial (upper) and fungal (lower) communities by inorganic fertilization, organic manure, crop rotation, and soil depth. Values < 0 were not shown. *** indicated significance at p < 0.001. (C) Taxonomic profiles of bacterial and fungal communities at phylum level. Bacteria and fungi phyla with relative abundances lower than 1% were summarized with “others” (Bray-Curtis).
Table 3. Results of PERMANOVA and BETADISP testing for differences in multivariate dispersion between fertilization regimes, crop rotation, and soil depth in bacterial and fungal communities.
Regression of Bray-Curtis dissimilarity with soil variables revealed the influence of soil C and N on microbial community compositions (Table 4). For soil microbial communities in rice fields, the contents of soil TN, BRC, BRN, FRC, and FRN significantly correlated with bacterial and fungal community compositions, and SOC only affected soil fungal community compositions (R2 = 0.66, p < 0.001). In the wheat fields, FRC and FRN (bacteria R2 = 0.23, p < 0.05; fungi R2 = 0.26, p < 0.01) were significantly correlated with microbial community compositions, and fungi was also affected by SOC (R2 = 0.53, p < 0.001) and TN (R2 = 0.56, p < 0.001).
Correlations Between Living and Dead Microorganism
Considering the significant correlations between microbial community compositions and microbial residues, we constructed a PLS-PM to access the direct and indirect effects of fertilization regimes, crop season, and soil depth on the contents of SOC and TN by affecting living and dead microorganisms in the soil (Figure 2). We found that inorganic fertilization increased soil SOC contents directly and indirectly by positively correlating with soil bacterial community compositions and increasing the contents of microbial residues C and N. Both organic fertilization and microbial communities exerted indirect positive impacts on SOC and TN contents. Relative to bacterial residues, fungal residues made more contributions to the contents of SOC and TN.
Figure 2. The partial least squares path models illustrating the effects of fertilization regimes, crop rotation, and soil depth on soil C and N by influencing the distribution of living and dead microorganism. (A) Bacterial and fungal communities were represented by the PC1 and PC2 of bacterial and fungal PCoA, respectively. Bacterial residues contained bacterial residues C and N. The blue and red lines indicated significant positive and negative correlations, respectively, while gray dash lines indicated insignificant correlation. The data on lines indicated the coefficient, and R2 indicated the coefficients of determination of the variables. GoF, goodness of fit. The cumulative bar plot displayed the standardized total effects of variables on SOC (B) and TN (C). The asterisk indicated significant differences, *p < 0.05, **p < 0.01, ***p < 0.001. BPC1, sample values in PC1 axis of bacterial PCoA; BPC2, sample values in PC2 of bacterial PCoA; FPC1, values in PC1 of fungal PCoA; FPC2, values in PC2 of fungal PCoA; BRC, bacterial residues C; BRN, bacterial residues N; FRC, fungal residues C; FRN, fungal residues N.
To further investigate the correlations between soil-specific microorganisms and microbial residues, we conducted a network analysis (Figure 3). Results revealed that both C and N in soil and microbial residues showed significant links to specific microbial OTUs. For the bacterial community in the rice fields, most OTUs affiliated to the phyla Actinobacteria, Acidobacteria, Chloroflexi, and Gammaproteobacteria showed significantly positive correlations with microbial residues, while the OTUs affiliated to the phyla Nitrospirae and Myxococcota showed significantly negative correlations with the microbial residues. On the contrary, in the wheat fields, there were fewer correlations between bacterial OTUs and bacterial residues. OTUs affiliated to Acidobacteria and Methylomirabilota were positively correlated with the fungal residues, while OTUs affiliated with Actinobacteria, Gammaproteobacteria, Alphaproteobacteria, and Gemmatimonadetes were negatively correlated. Considering their relative abundances, only Actinobacteria, which correlated positively with fungal residues, showed relatively higher abundance (1.85–6.42%, Supplementary Figure 1). Among all the identified bacterial genera, Gaiella (0.28–0.93%) and Marmoricola (0.17–0.86%) showed relatively higher abundances. Specifically, fertilization regimes and soil depth affected the genus Gaiella, which showed relatively higher abundances in the upper soil layer with M and NPKM treatments (Supplementary Figure 2). For the fungal community, most OTUs affiliated with the class Eurotiomycetes presented negative correlations with bacterial residues, and OTUs from the class Sordariomycetes showed negative links to the fungal residues, including the genera Cosmospora, Nigrospora, and Zopfiella. All these genera showed the lowest relative abundances in NPKM treatment than other fertilization treatments (Supplementary Figure 3). However, based on their relative abundances, only the class Pezizomycetes, which correlated positively with the fungal residues, presented a relatively higher abundance (0.22–61.09%), which was primarily attributed to the genus Pseudaleuria.
Figure 3. Association networks of microbial OTUs and soil variables. Bacterial networks in rice (A) and wheat (C) season, and fungal networks in rice (B) and wheat (D) season. OTUs are represented by circles, and soil variables by rectangle. Nodes’ sizes reflect their degrees, and are colored based on their taxonomic information. Green and red edges indicate positive and negative correlations, respectively. Wider edges reflect stronger correlations. BPC1, sample values in PC1 axis of bacterial PCoA; BPC2, sample values in PC2 of bacterial PCoA; FPC1, values in PC1 of fungal PCoA; FPC2, values in PC2 of fungal PCoA; BRC, bacterial residues C; BRN, bacterial residues N; FRC, fungal residues C; FRN, fungal residues N. Abbreviations: see Figure 2.
Compared with fertilization regimes, crop season, and soil depth, only fertilization regimes significantly affected the relative abundance of Psedaleuria (Figures 4B, P < 0.001), which showed higher abundances in M and NPKM with an average abundance of 34.1%. Considering its high relative abundance and significant correlations with fungal residues, we constructed the PLS-PM (Figure 4A) to further investigate the possible effects on microbial residues caused by Psedaleuria. Results revealed that the relative abundance of Psedaleuria was significantly positively correlated with fungal residues, which then indirectly increased the stocks of SOC and TN (Figures 4C,D, p < 0.01).
Figure 4. The partial least squares path models illustrating the effects of fertilization regimes, crop rotation, and soil depth on soil C and N by influencing the distribution of Pseudaleuria and microbial residues (for more details, see Figure 2). (B) The relative abundance of Pseudaleuria in different treatments and its difference between fertilization regimes, crop season and soil depth were test by LSD t-test. The different letters indicate significant difference. The standardized total effects of parameters on SOC (C) and TN (D). Asterisks indicate significant differences. *p < 0.05, **p < 0.01, ***p < 0.001.
Discussion
Both empirical studies and conceptual models have revealed that the incorporation of microbial biomass components into soils via microbial residues is large, which suggests that microbial inputs play a greater role in soil C sequestration (Lehmann and Kleber, 2015; Liang et al., 2019). Soil microorganisms convert plant C into microbial cell components and produce by-products in the process of decomposing plant residues (Liang et al., 2017). And compared with other cellular components, i.e., metabolites, microbial cell wall components are more stable due to their high polymerization (Gunina et al., 2017). After the death of microorganisms, the remnants of dead microbial cells serve as power for these biogeochemical processes because their chemical constituents persist as one of the largest C pools (Liang et al., 2017; Sokol et al., 2022). Soil bacteria and fungi are the dominant decomposers in soil, and their distinct physiologies are likely to differentially influence the rates of soil C cycling (Waring et al., 2013). In the present study, the long-term fertilization significantly stimulated the accumulations of soil bacterial residues, especially the application of organic manure, which was consistent with previous findings (Ding and Han, 2013; Faust et al., 2017). This was primarily because the addition of organic manure increased the nutrient resources and promoted the growth of soil bacteria, which preferred labile substrate (Schneider et al., 2010; Zhang et al., 2015). In addition, the microorganism carried in organic manure could also promote the accumulations of microbial residues (Ye et al., 2019; Chen et al., 2020).
The distinct microhabitats within the rice-wheat rotation also affected the compositions of microbial communities and their residues. In the present study, the contents of microbial residues in the wheat fields were much higher than the contents in the rice fields, as well as the contributions of microbial residues to soil C and N pools, especially the fungal residues. This is because the anaerobic environment in paddy soils inhibits the growth and metabolism of the fungal community (Moritz et al., 2009). Some studies have suggested that the high oxygen availability in the wheat fields enhances the turnover of microbial biomass, and increases the recolonization and decomposition of soil organic matter (Roth et al., 2011), which is inconsistent with our result. However, some studies demonstrated that high turnover rate of microbial biomass increased the incorporation of soil C and N into microbial residues, and higher temperatures in the rice season also enhanced the decomposition of SOC (Ma et al., 2018; Ding et al., 2019).
Microbial residues made a great contribution to soil C pools, and fungal residues contributed more than bacterial residues, which account for about 40% in most soils. The possible reason for that might be as follows. (1) Compared with bacteria, fungi have higher rates of growth and C assimilation, leading to a relatively high efficiency in utilizing C-containing substrates (Schneider et al., 2012). (2) Fungi are the main decomposer of cellulose, hemicellulose, and lignin, and can poly-condense the products into more stable C substrates (Roth et al., 2011; Schneider et al., 2012). (3) The amounts of the fungi hyphae are relatively huge, and their cell wall components are resistant to degradation. (4) Bacteria have relatively faster metabolic rates, and their residues are mostly readily degradable and easily utilized by microorganisms (Joergensen, 2018). Moreover, we found that crop season strongly affected the contribution of fungal resides to soil C and N pools, possibly due to the aerobic microhabitat in the wheat fields, which favored the fungal growth. Although no significant difference in the amounts of microbial residues was found between soil depths, higher ratios of microbial residues to SOC were observed in the upper soil layer, which indicated that soil microorganisms contributed more to soil C pools in surface soil layer.
The contents and components of microbial residues are also influenced by the soil microbial communities. A previous study has demonstrated that compared with other components of soil organic matter, soil microorganisms preferentially decompose and metabolize microbial residues, resulting in a rapid turnover of microbial residues (Zeglin and Myrold, 2013). As the soil-dominant microorganisms, bacteria and fungi contribute to the soil microbial residues that largely depend on their divergent responses to soil environment (Xia et al., 2019). Therefore, differences in the soil microbial communities directly affect the contents and compositions of soil microbial residues. In the present study, our results revealed that the long-term fertilization, especially the application of organic manure, favored the growth of soil microbial communities, and then increased the contents of microbial residues, leading to an enlarged soil C and N storage. And compared with bacterial residues, fungal residues contributed more to soil C and N pools. This is because the components in the bacterial cell walls, such as peptidoglycan, are easily degraded into labile C substrates, which are unfavorable to soil C sequestration (Joergensen, 2018). The SOC and TN were intrinsically coupled with the accumulations of microbial residues (Roth et al., 2011), which also supported our result.
Considering the strong correlations between microbial community compositions and their residues, we identified that some OTUs with relatively high abundances respond strongly to the accumulations of microbial residues. The much more complicated and significant correlations between FRC, FRN, and microbial OTUs indicated the stronger relationships between soil microbial communities and fungal residues. Among all microorganisms, the OTUs belonging to Actinobacteria showed the highest relative abundances, possibly due to the relative higher abundances of Gaiella and Marmoricola. However, their functions in mediating microbial residues were significantly weaker than fungal communities, especially the phyla Ascomycota, which contained a large proportion of saprotrophic fungi (Gourmelon et al., 2016). Previous studies have demonstrated that saprotrophic fungi are specialized in decomposing recalcitrant, ligno-cellulose containing materials, which involved in soil C cycling (De Boer et al., 2005; Lehmann and Rillig, 2015). Our study revealed that the fungal OTUs correlated with fungal residues were mostly classified into saprotrophic fungi, including Zopfiella, Cosmospora, Nigrospora, and Pseudaleuria (Tedersoo et al., 2014). Furthermore, the highest relative abundance of Pseudaleuria in all fungal OTUs was relatively higher in soils treated with organic manure, which is consistent with a previous study (Xiang et al., 2020). This was possibly because of the critical role of saprophytic fungi Pseudaleuria functioning in SOC cycling by secreting extracellular enzymes (Xu et al., 2012). Considering the strong and positive correlations with fungal residues and its relatively higher abundance, we suggested that the application of organic manure favored the growth of Pseudaleuria, then stimulated the accumulations of fungal residues, and finally benefited soil C sequestration.
Conclusion
Microbial community compositions influence soil C sequestration by shaping the microbial set and by triggering changes in the concentrations and compositions of microbial residues. Our study found that the long-term fertilization, especially the application of organic manure combined with inorganic fertilizers, significantly increased the contents of C and N in soil and microbial residues. Compared with paddy soils in the rice season, upland soils in the wheat season accumulated more fungal residues by favoring the growth of fungi. Soil bacterial and fungal communities showed strong and complicated correlations with the contents of microbial residues, particularly the positive correlations between certain fungal taxa Pseudaleuria and fungal residues. These results suggested that the long-term application of organic manure significantly favored the growth of soil bacterial and fungal communities, and then stimulated the accumulations of microbial residues, particularly in fungal residues, leading to the enlargement of soil C and N pools. The quantitative impact of Pseudaleuria exerted on the contents and compositions of microbial residues needs further research.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
JC contributed to the sampling work and laboratory experiments, analyzed all the results of the experiments, and wrote the first draft of this manuscript. DS and HL contributed to the statistical analysis. DL helped to maintain the long-term experiment and support the sampling work. XW and JS revised this manuscript. WZ and GL guided and proofread the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the National Key Research and Development Program of China (2016YFD0200109) and the National Natural Science Foundation of China (31471943).
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.
Acknowledgments
We thank all colleagues who contributed their efforts to maintain this long-term experiment and support our sampling work.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.854216/full#supplementary-material
Supplementary Figure 1 | The relative abundance of microbial OTUs linked to microbial residues in the phylum/class level.
Supplementary Figure 2 | The relative abundance of microbial OTUs linked to microbial residues in the genus level.
Supplementary Figure 3 | The relative abundances of some genus under different fertilization regimes. Significance was test by ANOVA with LSD t-test.
Footnotes
- ^ https://github.com/networkx/networkx
- ^ http://data.cma.cn/
- ^ https://ccb.jhu.edu/software/FLASH/index.shtml
- ^ https://github.com/OpenGene/fastp
- ^ https://submit.ncbi.nlm.nih.gov/subs/sra
- ^ https://github.com/vegandevs/vegan
- ^ https://matplotlib.org/
References
Appuhn, A., and Joergensen, R. (2006). Microbial colonization of roots as a function of plant species. Soil Biol. Biochem. 38, 1040–1051. doi: 10.1016/J.SOILBIO.2005.09.002
Chen, G. P., Ma, S. H., Tian, D., Xiao, W., Jiang, L., Xing, A. J., et al. (2020). Patterns and determinants of soil microbial necromass from tropical to boreal forest. Soil Biol. Biochem. 151:108059. doi: 10.1016/j.soilbio.2020.108059
Chen, J., Guo, Q., Liu, D., Hu, C., Sun, J., Wang, X., et al. (2020). Composition, predicted functions, and co-occurrence networks of fungal and bacterial communities_Links to soil organic carbon under long-term fertilization in a rice-wheat cropping system. Eur. J. Soil Biol. 100:103226. doi: 10.1016/j.ejsobi.2020.103226
Chen, X., Hu, Y., Xia, Y., Zheng, S., Ma, C., Rui, Y., et al. (2021). Contrasting pathways of carbon sequestration in paddy and upland soils. Glob. Change Biol. 27, 2478–2490. doi: 10.1111/gcb.15595
De Boer, W., Folman, L. B., Summerbell, R. C., and Boddy, L. (2005). Living in a fungal world: impact of fungi on soil bacterial niche development. FEMS Microbiol. Rev. 29, 795–811. doi: 10.1016/J.FEMSRE.2004.11.005
Ding, X., Chen, S., Zhang, B., Liang, C., He, H., and Horwath, W. R. (2019). Warming increases microbial residue contribution to soil organic carbon in an alpine meadow. Soil Biol. Biochem. 135, 13–19. doi: 10.1016/J.SOILBIO.2019.04.004
Ding, X., and Han, X. (2013). Effects of long-term fertilization on contents and distribution of microbial residues within aggregate structures of a clay soil. Biol. Fertil. Soil. 50, 549–554. doi: 10.1007/s00374-013-0867-6
Engelking, B., Flessa, H., and Joergensen, R. G. (2007). Shifts in amino sugar and ergosterol contents after addition of sucrose and cellulose to soil. Soil Biol. Biochem. 39, 2111–2118. doi: 10.1007/s00374-016-1165-x
Faust, S., Heinze, S., Ngosong, C., Sradick, A., Oltmanns, M., Raupp, J., et al. (2017). Effect of biodynamic soil amendments on microbial communities in comparison with inorganic fertilization. Appl. Soil Ecol. 114, 82–89. doi: 10.1016/J.APSOIL.2017.03.006
Fierer, N. (2017). Embracing the unknown: disentangling the complexities of the soil microbiome. Nat. Rev. Microbiol. 15, 579–590. doi: 10.1038/nrmicro.2017.87
Frey-Klett, P., Burlinson, P., Deveau, A., Barret, M., Tarkka, M., and Sarniguet, A. (2011). Bacterial-fungal interactions: Hyphens between agricultural, clinical, environmental, and food microbiologists. Microbiol. Mol. Biol. Rev. 75, 583–609. doi: 10.1128/MMBR.00020-11
Gies, H., Hagedorn, F., Lupker, M., Montluçon, D., Haghipour, N., and Sophia Van, et al. (2021). Millennial-age glycerol dialkyl glycerol tetraethers (GDGTs) in forested mineral soils: 14C-based evidence for stabilization of microbial necromass. Biogeosciences 18, 189–205. doi: 10.5194/BG-18-189-2021
Gourmelon, V., Maggia, L., Powell, J. R., Gigante, S., Hortal, S., Gueunier, C., et al. (2016). Environmental and geographical factors structure soil microbial diversity in New Caledonian ultramafic substrates: a metagenomic approach. PLoS One 11:e0167405. doi: 10.1371/journal.pone.0167405
Gunina, A., Dippold, M., Glaser, B., and Kuzyakov, Y. (2017). Turnover of microbial groups and cell components in soil: C-13 analysis of cellular biomarkers. Biogeosci. 14, 271–283. doi: 10.5194/BG-14-271-2017
He, H., Zhang, W., Zhang, X., Xie, H., and Zhuang, J. (2011). Temporal responses of soil microorganisms to substrate addition as indicated by amino sugar differentiation. Soil Biol. Biochem. 43, 1155–1161. doi: 10.1016/J.SOILBIO.2011.02.002
Hu, Y. T., Zheng, Q., Zhang, S. S., Noll, L., and Wanek, W. (2018). Significant release and microbial utilization of amino sugars and D-amino acid enantiomers from microbial cell wall decomposition in soils. Soil Biol. Biochem. 123, 115–125. doi: 10.1016/j.soilbio.2018.04.024
Jiao, S., Yang, Y. F., Xu, Y. Q., Zhang, J., and Lu, Y. H. (2019). Balance between community assembly process mediates species coexistence in agricultural soil microbiomes across eastern China. ISME J. 14, 202–216. doi: 10.1038/s41396-019-0522-9
Joergensen, R. G. (2018). Amino sugars as specific indices for fungal and bacterial residues in soil. Biol. Fert. Soil. 54, 5959–5568. doi: 10.1007/s00374-018-1288-3
Keenan, T. K., and Williams, C. A. (2018). The terrestrial carbon sink. Nat. Annu. Rev. Environ. Resour. 43, 219–243. doi: 10.1146/annurev-environ-102017-030204
Kögel-Knabner, I. (2002). The macromolecular organic composition of plant and microbial residues as inputs to soil organic matter. Soil Biol. Biochem. 34, 139–162. doi: 10.1016/S0038-0717(01)00158-4
Lehmann, A., and Rillig, M. (2015). Understanding mechanisms of soil biota involvement in soil aggregation: a way forward with saprobic fungi? Soil Biol. Biochem. 88, 298–302.
Lehmann, J., and Kleber, M. (2015). The contentious nature of soil organic matter. Nature 528, 60–68. doi: 10.1038/nature16069
Liang, C., Amelung, W., Lehmann, J., and Kästner, M. (2019). Quantitative assessment of microbial necromass contribution to soil organic matter. Glob. Change Biol. 25, 3578–3590. doi: 10.1111/gcb.14781
Liang, C., Schimel, J., and Jastro, J. D. (2017). The importance of anabolism in microbial control over soil carbon storage. Nat. Microbiol. 2:17105. doi: 10.1038/nmicrobiol.2017.105
Ma, T., Zhu, S., Wang, Z., Chen, D., Dai, G., Feng, B., et al. (2018). Accumulation of nitrogen and microbial residues during 2000 years of rice paddy and non-paddy soil development in the Yangtze River Delta China. Glob. Change Biol. 17, 3405–3417. doi: 10.1111/J.1365-2486.2011.02500.X
Miltner, A., Bombach, P., Schmid-Brücken, B., and Kästner, M. (2012). SOM genesis: microbial biomass as a significant source. Biogeochem. 111, 41–55. doi: 10.1007/s10533-011-9658-z
Moritz, L., Liang, C., Waggi, R., Kitayama, K., and Balser, T. C. (2009). Vertical distribution and pools of microbial residues in tropical forest soils formed from distinct parent materials. Biogeochem. 92, 83–94. doi: 10.1007/S10533-008-9264-X
Nachtergaele, F. O. (2017). Classification Systems: FAO. Ref. Mod. Earth Syst. Environ. Sci. 2017:2. doi: 10.1016/B978-0-12-409548-9.10520-2
Naylor, D., Sadler, N., Bhattacharjee, A., Graham, E. B., Anderton, C. R., McClure, R., et al. (2020). Soil microbiomes under climate change and implications for carbon cycling. Nat. Annu. Rev. Environ. Resour. 45, 29–59. doi: 10.1146/annurev-environ-012320-082720
Nazir, R., Warmink, J. A., Boersma, H., and van Elsas, J. D. (2009). Mechanisms that promote bacterial fitness in fungal-affected soil microhabitats. FEMS Microbiol. Ecol. 71, 169–185. doi: 10.1111/j.1574-6941.2009.00807.x
Pan, G. X., Li, L. Q., Wu, L. S., and Zhang, X. H. (2003). Storage and sequestration potential of topsoil organic carbon in China’s paddy soils. Glob. Chang. Biol. 10, 79–92. doi: 10.1111/j.1365-2486.2003.00717.x
Parsons, J. W. (1981). Chemistry and distribution of amino sugars in soils and soil organisms. Soil Biochem. 5, 197–227. doi: 10.1128/AEM.00956-20
Peixoto, L., Elsgaard, L., Rasmussen, J., Kuzyakov, Y., Banfield, C. C., Dippold, M. A., et al. (2020). Decreased rhizodeposition, but increased microbial carbon stabilization with soil depth down to 3.6 m. Soil Biol. Biochem. 150:108008. doi: 10.1016/J.SOILBIO.2020.108008
Roth, P. J., Lehndorff, E., Cao, Z., Zhuang, S. Y., Bannert, A., Wissing, L., et al. (2011). Accumulation of nitrogen and microbial necromass during 2000 years of rice paddy and non-paddy soil development in the Yangtze River Delta, China. Glob. Change Biol. 17, 3405–3417. doi: 10.1111/j.1365-2486.2011.02500.x
Sanchez, G., Trinchera, L., and Russolollo, G. (2015). plspm: Tools for partial least squares path modeling (PLS-PM). R package version 0.4.9.
Schimel, J. P., and Schaeffer, S. M. (2012). Microbial control over carbon cycling in soil. Front. Microbiol. 3:348. doi: 10.3389/fmicb.2012.00348
Schneider, M., Scheel, T., Mikutta, R., van Hees, P., Kaiser, K., and Kalbitz, K. (2010). Sorptive stabilization of organic matter by amorphous Al hydroxide. Geochimica et Cosmochimica Acta 74, 1606–1619. doi: 10.1016/J.GCA.2009.12.017
Schneider, T., Keiblinger, K. M., Schmid, E., Sterflinger-Gleixner, K., Ellersdorfer, G., Roschitzki, B., et al. (2012). Who is who in litter decomposition? Metaproteomics reveals major microbial players and their biogeochemical functions. ISME J. 6, 1749–1762. doi: 10.1038/ismej.2012.11
Sokol, N., Slessarev, E., Marschmann, G., Nicolas, A., Blazewicz, S., Brodie, E., et al. (2022). Life and death in the soil microbiome: how ecological processes influence biogeochemistry. Nat. Rev. Micobiol. 2022:695. doi: 10.1038/s41579-022-00695-z
Tang, X. L., Zhao, X., Bai, Y. F., Tang, Z. Y., Wang, W. T., Zhao, Y. C., et al. (2018). Carbon pools in China’s terrestrial ecosystems: new estimates based on an intensive field survey. PNAS 115, 4021–4026. doi: 10.1073/pnas.1700291115
Tedersoo, L., Bahram, M., and Põlme, S. (2014). Global diversity and geography of soil fungi. Sci. 346:1256688. doi: 10.1126/science.1256688
Vollmer, W., Joris, B., Charlier, P., and Foster, S. (2008). Bacterial peptidoglycan (murein) hydrolases. FEMS Microbiol. Rev. 32, 259–286. doi: 10.1111/j.1574-6976.2007.00099.x
Waring, B., Averill, C., and Hawkes, C. (2013). Differences in fungal and bacterial physiology alter soil carbon and nitrogen cycling: insights from meta-analysis and theoretical models. Ecol. Lett. 16, 887–894. doi: 10.1111/ele.12125
Wiesmeier, M., Urbanski, L., Hobley, E., Lang, B., von Lützow, M., Martin-Spiotta, E., et al. (2019). Soil orgnic carbon storage as a key function of soils – A review of drivers and indicators at various scales. Geoderma 333, 149–162. doi: 10.1016/J.GEODERMA.2018.07.026
Xia, Y. H., Chen, X. B., Hu, Y. J., Zheng, S. M., Ning, Z., Guggenberger, G., et al. (2019). Contrasting contribution of fungal and bacterial necromass to organic carbon accumulation in paddy soils across eastern China. Biol. Fertil. Soils 55, 767–776. doi: 10.1007/s00374-019-01390-7
Xiang, X., Liu, J., Zhang, J., Li, D., Xu, C., and Kuzyakov, Y. (2020). Divergence in fungal abundance and community structure between soils under long-term mineral and organic fertilization. Soil Till. Res. 196:104491. doi: 10.1016/j.still.2019.104491
Xu, L., Ravnskov, S., Larsen, J., Nilsson, R., and Nicolaisen, M. (2012). Soil fungal community structure along a soil health gradient in pea fields examined using deep amplicon sequencing. Soil Biol. Biochem. 46, 26–32. doi: 10.1016/J.SOILBIO.2011.11.010
Ye, G. P., Lin, Y. X., Kuzyakov, Y., Liu, D. Y., Luo, D. Y., Lindsey, S., et al. (2019). Manure over crop necromass increases soil organic matter but decreases microbial necromass relative contribution in upland Ultisols: Results of a 27-year field experiment. Soil Biol. Biochem. 134, 15–24. doi: 10.1016/J.SOILBIO.2019.03.018
Zeglin, L., and Myrold, D. (2013). Fate of decomposed fungal cell wall materials in organic horizons of old-growth douglas-fir forest soils. Soil Sci. Soci. Am. J. 77, 489–500. doi: 10.2136/SSSAJ2012.0204
Zhang, Q., Zhou, W., Liang, G. Q., Sun, J. W., Wang, X. B., and He, P. (2015). Distribution of soil nutrients, extracellular enzyme activities and microbial communities across particle-size fractions in a long-term fertilizer experiment. Appl. Soil Ecol. 94, 59–71. doi: 10.1016/J.APSOIL.2015.05.005
Keywords: microbial community, long-term fertilization, soil amino sugars, network analysis, partial least squares path model
Citation: Chen J, Song D, Luan H, Liu D, Wang X, Sun J, Zhou W and Liang G (2022) Living and Dead Microorganisms in Mediating Soil Carbon Stocks Under Long-Term Fertilization in a Rice-Wheat Rotation. Front. Microbiol. 13:854216. doi: 10.3389/fmicb.2022.854216
Received: 13 January 2022; Accepted: 09 May 2022;
Published: 10 June 2022.
Edited by:
Leonardo Erijman, Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), ArgentinaReviewed by:
Eva L. M. Figuerola, Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), ArgentinaMilena Gonzalo, Université Bourgogne Franche-Comté, France
Copyright © 2022 Chen, Song, Luan, Liu, Wang, Sun, Zhou and Liang. 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: Guoqing Liang, liangguoqing@caas.cn; Wei Zhou, zhouwei02@caas.cn