Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Microbiol., 16 July 2025

Sec. Terrestrial Microbiology

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1602633

This article is part of the Research TopicMicrobial Solutions for Restoring Depleted SoilsView all 7 articles

Organic material mulching regulated core microbial groups to promote soil carbon and nitrogen cycling and improve faba bean productivity under a triple-cropping system in purple soil hilly region of southwest China

Ke RenKe Ren1Wenfeng SongWenfeng Song1Zehui WeiZehui Wei2Lixia SongLixia Song1Ming LiuMing Liu1Yuling ZhouYuling Zhou1Yuzhuo ZhenYuzhuo Zhen1Xinyao WuXinyao Wu1Kaiyuan GuKaiyuan Gu1Khanom Simarani
Khanom Simarani2*Longchang Wang
Longchang Wang1*
  • 1College of Agronomy and Biotechnology Southwest University Engineering Research Center of South Upland Agriculture, Ministry of Education, Chongqing, China
  • 2Division of Microbiology, Institute of Biological Sciences, Faculty of Science, Universiti Malaya, Kuala Lumpur, Malaysia

Background: Intensive agricultural production leads to the reduction of soil carbon (C) and nitrogen (N) reserves, and organic material mulching (OMM) can improve microbial community structure and promote C and N accumulation. The multi-cropping system based on legumes can provide abundant organic mulching material and improve soil quality, but the mechanism by which OMM provides ecological benefits via C and N cycling in this system is still unclear.

Methods: In this study, a field experiment of organic mulch under the triple-cropping system of faba bean (Vicia faba L.)-corn (Zea mays L.)-sweet potato (Ipomoea batatas Lam.) was performed. Four treatments were set up: without organic material mulching (CK), straw mulching (S), milk vetch (Astragalus sinicus L.) mulching (M), and straw and milk vetch mulching (SM). The dynamic changes in soil aggregates, C and N contents, enzyme activities, microbial communities and faba bean productivity under organic mulching were studied for 2 years (from October 2020 to May 2022).

Results: The results demonstrated that supplementation of OMM (S, M, and SM) significantly improved the stability of soil aggregates, the nutrient (C and N) contents, and the activities of acquiring enzymes compared to CK. OMM promoted the establishment of key microbial communities dominated by Actinobacteria, Bacteroidetes, Ascomycota, and Basidiomycota. Predicted functional profiles based on PICRUSt and FUNGuild analyses suggest possible upregulation of up-regulated genetic information processing, metabolism, and organismal systems functional pathways. Additionally, the enhancement of soil C and N cycling efficiency may be associated with an increase in the proportion of saprotrophs and symbiotrophs. Both the random forest model (RFM) and partial least squares path model (PLS-PM) demonstrated that OMM increased faba bean productivity by improving soil microbial diversity and the efficiency of C cycling.

Conclusion: This study highlighted that OMM could promote C and N cycling by regulating core microbial groups, thereby improving the faba bean productivity in dryland of purple soil hilly region of southwest China.

1 Introduction

Agroindustrial waste such as crop straw and legume green manure are high-quality organic mulch materials, which play an important role in improving soil ecological environment, soil quality, and crop yield (Chen S. H. et al., 2022). Crop straw to corn (Zea mays L.), faba bean (Vicia faba L.) and rape (Brassica campestris L.) is much abundant as a resource of renewable cellulose in farmland resources and yet it is underutilized (Yu S. et al., 2022). Straw is rich in cellulose, hemicellulose, lignin, and ash, as well as a rich source of carbon (C) (Wang et al., 2021). According to the agricultural production statistics of the Food and Agriculture Organization (FAO) (FAOSTAT, 2021) and the crop straw index (Zhang et al., 2021), the total global straw output in 2020 was about 5.849 billion tons, with China accounting for 17.85% which is the largest straw producer in the world. In 2021, the return rate of crop straw in China was about 61.82%, and organic material mulching (OMM) holds immense potential for scientific research (Yang et al., 2022). Since crop straws are characterized by a component of high carbon (C) to nitrogen (N) ratio (C:N), their return to the field leads to a competition between soil microorganisms and crops for mineral N, resulting in high greenhouse gas emissions (Hassan et al., 2024; Xin et al., 2024). Further, it can lead to a significant imbalance in the soil microbial community structure (Xie et al., 2022b). Therefore, determining the approach to optimize the interactions among soil microbes, soil quality, and crops after crop straw return is a crucial prerequisite to ensure the scientific utilization of straw resources. Milk vetch (Astragalus sinicus L.) is a high-quality legume green fertilizer, which can contribute to supplementing the nutrient requirements of cultivated land and replace part of fertilizer usage upon its return to the field (Ma et al., 2020). Moreover, it enriches the diversity of microbial community structure (Xia et al., 2023). The most important benefit of milk vetch is that it can reduce soil C:N (Fang et al., 2024), which helps in alleviating the increase in C:N after the return of straw to the field, and also promotes the release of mineral N by microbial decomposition of organic materials (Xie et al., 2022a). However, the mineralization of straw and milk vetch is a relatively slow process in the soil (Xu et al., 2021). Currently, there are limited long-term dynamic studies to understand the effects of mulching different crop straws in combination with milk vetch on soil microorganisms and soil quality.

Purple soil region of southwest China is a mountainous and hilly area, with high temperatures, high humidity, low levels of sunlight, and regular human activities (Tan et al., 2021), which has resulted in a fragile ecological environment and lack of cultivated land resources (Chu et al., 2020). At the same time, the extensive area under spring and summer grain crops in purple soil hilly region of southwestern China and the fallow time after harvest leads to a low multiple cropping index, resulting in a significant waste of natural resources (Chu et al., 2020). The multi-cropping system, commonly practiced in the purple soil hilly region of southwest China (Wu et al., 2023), effectively utilizes light, heat, water, and limited cultivated land resources to enhance farmland ecosystem diversity and ensure food production (He et al., 2023; Nasar et al., 2023). However, a multi-cropping system requires large amounts of nutrients for land cultivation. In addition, complex issues such as a wide range of crop options, greater niche competition, unreasonable fertilization, and extensive tillage practices are associated with this system (Chen et al., 2022). The selection of core crops based on the local conditions is vital for optimizing the soil microecological environment and improving soil quality. Faba bean is an important crop in the multi-cropping system in purple soil hilly region of southwest China, where the area under faba bean production is the largest in China (Zhou et al., 2022). The cultivation of faba bean provides multiple benefits such as improvement in the multiple cropping index, leguminous crop-based N fixation, enhancing the growth of numerous beneficial rhizosphere microorganisms, and stabilizing the soil ecosystem (Karkanis et al., 2018). More importantly, the faba bean as the major crop of the multi-cropping system can improve the soil quality of the fragile agricultural land and enhance the cropping index in purple soil hilly region of southwest China (Jithesh et al., 2024; Nasar et al., 2024). On the other hand, the multi-cropping system has significant advantages such as the production of straws, which are rich sources for OMM. Therefore, it is essential to investigate the multi-cropping system with the faba bean as the core crop for improving the soil ecosystem and stabilizing grain production in the purple soil hilly region of southwest China.

Soil microorganisms are key biological factors driving biotic and abiotic processes in farmland ecosystems (Barq et al., 2021). They play a vital role in soil N cycling (Zhang et al., 2022), contribute to collaborative crop resilience processes (Grobelak et al., 2018), and maintain ecosystem stability (Bhaduri et al., 2022). Straw mulching can regulate the changes in soil microbial community structure and function (Xu et al., 2023). The variations in organic materials composition and C:N ratio of different crop straws is the major reason for the changes in the soil microecological environment (Chen L.M. et al., 2022). Su et al. (2020) demonstrated the presence of copiotrophic bacteria mainly in the soil under corn straw return, while oligotrophic bacteria were mainly found in the soil under wheat straw return. Mulching and returning the crop straw to the field can enhance microorganisms that play key roles in N fixation, nitrification, and denitrification (Dou et al., 2023), stimulate the activity of sulfur-reducing bacteria, and reduce soil acidification (Zhong et al., 2021). In particular, milk vetch can effectively improve the N-use efficiency of soil microorganisms and maintain a suitable C:N ratio under crop straw mulching (Li et al., 2021). Gao et al. (2021) found that the return of milk vetch to rice soil had a significant impact on the community and function of soil bacteria, increasing the abundance of Actinobacteria and Firmicutes, and also enhancing nutrient cycling. All studies mentioned above indicated that the impacts of diverse organic materials inputs on soil microbial community structure and function are complex. Further research on the ecological benefits provided by the soil microorganisms under the multi-cropping system with different organic mulching is of great significance for the sustainable development of agricultural production in purple soil hilly region of southwest China.

According to previous studies, the ecological mechanism that explains the role of OMM in dynamically driving the soil nutrient cycling and changes in the microbial community structure that affect faba bean yield in the context of the triple cropping model in purple soil hilly region of southwest China is still unclear. Therefore, in this study, a field experiment was designed for two consecutive years in typical purple soil farmland in Chongqing, China. It was assumed that mulching straw combined with milk vetch onto the field has the potential to gradually improve the physicochemical properties and the diversity of microbial community in the soil, which leads to an increase in faba bean productivity. The objectives of this study were (1) to investigate the dynamic effects of straw and milk vetch mulching on the soil physicochemical properties and faba bean productivity; (2) to explore the dynamic changes in the stability of core microbial functions and molecular ecological networks (MENs) in faba bean soil following straw and milk vetch mulching; (3) to elucidate the ecological mechanism of increase in faba bean productivity associated with the changes in soil physicochemical properties and microbial communities driven by straw and milk vetch mulching. Briefly, this study under conservation tillage with OMM can provide a theoretical basis and practical guidance for the multi-cropping system in purple soil hilly region of southwest China to alleviate farmland soil erosion, improve ecological benefits, and stabilize sustainable grain production.

2 Materials and methods

2.1 Site description and experimental design

The experiment was conducted in the National Purple Soil Fertility and Fertilizer Benefit Monitoring Base (29°51′N, 106°27′E, and altitude 244 m) of Southwest University, Beibei District, Chongqing, China from October 2020 to May 2022. The experimental site has a subtropical humid climate with an average annual temperature of 18.3°C and an average annual rainfall of 1,105 mm. The soil type is purple soil, which is derived from the taupe purple sandy mudstone of the Jurassic Shaximiao Formation (J2s) (Zuo et al., 2020). The physicochemical properties of the soil in the top layer (0–15 cm) were pH—6.54, bulk density—1.1 g·cm−3, soil organic carbon (SOC)—13.50 g·kg−1, total nitrogen (TN)—1.18 g·kg−1, total phosphorus (TP)—0.9 g·kg−1, total potassium (TK)—20.38 g·kg−1, alkaline hydrolyzable nitrogen (AN)—101.68 g·kg−1, available phosphorus (AP)—74.98 g·kg−1, and available potassium (AK)—243.13 g·kg−1. The slope of the test field is shallow and the soil fertility is relatively uniform. This study was based on the faba bean-corn-sweet potato intercropping model in the multi-cropping system in dry land. In this intercropping model, the faba bean variety is the major local variety of small green faba bean, the corn variety is Zhengda 999, and the sweet potato (Ipomoea batatas Lam.) variety is Yuhongxin 98. The experiment was characterized by a single-factor randomized block design, and a total of four treatments was set up, namely, without organic material mulching (CK), straw mulching (S), milk vetch (Astragalus sinicus L.) mulching (M), and straw and milk vetch mulching (SM) (Supplementary Table S1). The treatments were replicated for three times, and each plot was 8 m long, 4 m wide, and 32 m2 in area. The crop establishment and management plans are presented in Figure 1A.

Figure 1
Diagram depicting crop rotations and growth cycles over time. Section A shows different crop combinations (CK, S, M, SM) in late October and late May. Section B illustrates a timeline from 2020 to 2022 with symbols representing the sowing and growth of faba beans, maize, Chinese milk vetch, and sweet potatoes. A legend describes symbols for seeds, plants, sowing, and mulching.

Figure 1. Mulching with organic materials (A) and faba bean-corn-sweet potato three-ripe preparation succession scheme (B). CK, without organic material mulching; S, straw mulching; M, milk vetch mulching; SM, straw and milk vetch mulching.

Faba bean was sown in late October each year, with two rows planted in each strip (row spacing of 30 cm, hole spacing of 30 cm, and hole depth of 5 cm). Subsequently, the seedlings were thinned and patched after emergence and harvested in early May of the following year. Corn was sown in late March each year, with two rows planted in each strip (row spacing 65 cm, plant spacing 30 cm, and seeding depth 15 cm), and harvested in late July. Sweet potato seedlings were transplanted using the direct insertion method in late May each year, with three rows planted in each strip (row spacing of 65 cm, hole spacing of 25 cm, 16 holes in each row, and 2 plants in each hole depth of 10 cm), and harvested in late October. Milk vetch was sown in uniform lines with bare seeds at a rate of 75 kg·ha−1, and intercropped with faba beans. The seeds were planted in late October each year, and the waxy layer on the surface of the seed coat was removed before sowing. The crop planting sequence and returning operation are presented in Figure 1B.

During the growth phase of each crop, compound fertilizer (N:P2O5:K2O = 16:16:16) was applied once as the base fertilizer at a rate of 200 kg·ha−1 and also applied to the holes at the time of sowing or transplantation. In late March every year, urea at a rate of 115 kg·ha−1 was applied as the flower and pod fertilizer for faba beans, and 0.2% potassium dihydrogen phosphate solution was sprayed on the leaf surface 2–3 times. In early May, 270 kg·ha−1 of urea was applied to corn plantation during the jointing stage and 150 kg·ha−1 of urea was applied to sweet potato plantation during the growth phase. Other field management practices applicable for each crop are the same as local conventional cultivation methods. Fresh samples of straw were randomly collected after crop harvesting in the current season and also during the blooming period of milk vetch. The samples were initially oven-dried at 105°C for 30 min and then further dried to constant weight at 65°C. After grinding, the nutrient contents of each sample of crop straw and milk vetch were analyzed (Supplementary Table S2).

2.2 Plant and soil sampling and analysis

Field soil samples were collected from 0 to 15 cm top layer using the point-centered quarter sampling method with a soil drill from May 1 to 3 in 2021 and 2022. According to the point-centered quarter sampling method, 1.5 kg of samples were collected and stored in a sealed bag. All the samples were immediately transported to the laboratory and kept at 4°C. The collected soil samples were screened using a 2 mm sieve to remove plant residues, fine roots, stones, and other debris. Subsequently, each sample was divided into three parts. One part of the fresh soil sample was stored at −80°C for examination of soil microorganisms. The other part of the fresh soil sample was stored at 4°C for the determination of soil enzyme activity. The remaining part of the soil sample was air-dried for soil aggregate measurement and chemical characterization.

2.3 Soil aggregate fractionation

The primary methods for evaluating the stability of soil aggregates are the wet-dry sieve method. The dry sieving method assesses the mechanical stability of aggregates, while the wet sieving method evaluates their water stability.

Separation of dry aggregates: the aggregate samples underwent size classification through dry sieving using an aggregate sieving apparatus (TPF-100, Shanghai Tomos Scientific Instrument Co., Ltd.). Approximately 200 g of soil samples were processed through a trio of sieves featuring different mesh sizes to obtain four distinct aggregate size fractions. The particles experienced oscillation for 10 min at a frequency of 200 cycles per minute with an amplitude of 2 mm, followed by weighing each particle size. The soil collected on the 2,000 μm, 250 μm, and 106 μm sieves was carefully gathered and categorized as large macro-aggregates (>2,000 μm), small macro-aggregates (250–2,000 μm), micro-aggregates (106–250 μm), and silt and clay particles (<106 μm), respectively (Wang et al., 2016).

Separation of wet aggregates: A mixed soil sample weighing 50 g was prepared according to the proportions of soil particle sizes determined by the dry sieving technique and placed on top of the sieves (dimensions: 2,000 μm, 250 μm, 106 μm). Deionized water was used to reach the water level, ensuring that it sufficiently submerged the sieves. The setup was then allowed to stand in the bucket for 10 min. Following this, the samples underwent oscillation at a frequency of 30 times per minute with an amplitude of 2 mm for a duration of 10 min. The aggregates from each sieve were subsequently collected by rinsing with deionized water and were then moved into aluminum boxes for drying in a thermostatically controlled electric oven set at 105°C.

To determine the distribution of soil aggregates, the proportion of each fraction based on dry weight was calculated. Additionally, metrics such as macro-aggregate content (R0.25), percentage of aggregate destruction (PAD), mean weight diameter (MWD), geometric mean diameter (GMD), and fractal dimensions (Ds) were derived as indicators of aggregate stability, following the Equations 15 outlined by Tagar et al. (2020) and Zhu et al. (2020).

R 0.25 = D R 0.25 / DM T × 100 %     (1)
PAD = ( D R 0.25 W R 0.25 ) / D R 0.25 × 100 %     (2)
MWD = i = 1 n X i × M i     (3)
GMD = exp [ i = 1 n M i × Ln X i i = 1 n M i ]     (4)
log 10 { M ( δ < X i ¯ ) M 0 } = ( 3 Ds ) log 10 X i ¯ X max     (5)

Where, D R 0.25 is the greater than 0.25 mm (250 μm) mechanical-stable aggregates mass. DM T is the mechanical-stable aggregate mass. W R 0.25 is the greater than 0.25 mm (250 μm) water-stable aggregates mass. Xi is the mean diameter (mm) of size i aggregates in any aggregate size fraction, and Mi is the weight of size i aggregates in that aggregate size fraction as a proportion of the dry weight of soil. M ( δ < X i ¯ ) is the cumulative soil aggregates mass (g) less than the average of two adjacent aggregate sizes, M0 is the sum of the mass of soil aggregates of each aggregate size (g), X i ¯ is the mean particle diameter (mm) of size i aggregates and is considered as the arithmetic mean of the upper and smaller sieve sizes, Xmax is the mean diameter of the aggregate with the maximum size, and Ds is the fractal dimension of soil aggregates (dimensionless).

2.4 Chemical analyses

Soil water content (SWC) was determined from mass loss after drying at 105°C for 24 h following the method by Zhao et al. (2023). The chemical properties of soil were measured as described by Bao (1999). Soil pH was measured by potentiometric method (soil: water ratio of 1:5), soil organic carbon (SOC) was determined by potassium dichromate capacity method, soil total nitrogen (TN) was determined by micro-Kjeldahl method, soil total phosphorus (TP) was determined by H2SO4–HClO4 digestion and molybdenum-antimony anti-colorimetric method, and soil total potassium (TK) was NaOH fusion-flame photometry.

Soil hydrolytic nitrogen (AN) was determined by alkaline diffusion method, and soil nitrate nitrogen (NO3–N) and ammonium nitrogen (NH4+–N) were determined by continuous flow analyzer method (Zhang et al., 2015). Soil available phosphorus (AP) was determined by sodium bicarbonate extraction-molybdenum-antimony resistance colorimetric method (Ahmed et al., 2020). Soil available potassium (AK) was determined by ammonium acetate extraction-flame photometer (Lu et al., 2017).

2.5 Soil enzyme assays

Extracellular enzyme activities were determined using an enzyme-linked immunosorbent assay (ELISA) kit (Jiangsu Meimian Industrial Co., Ltd., Yancheng, China) following the manufacturer’s instructions. Activities of six hydrolytic soil enzymes were measured, including two enzymes involved in C-acquisition: β-1,4-glucosidase, β-1,4-xylosidase; two enzymes involved in N-acquisition: urease, L-leucine aminopeptide; and one enzyme involved in P-acquisition: acid phosphatase. In addition, we measured activities of one oxidase involved in the degradation of recalcitrant organic C: peroxidase (Tian et al., 2022). All the above activities were represented by μmol·g−1·h−1 dry soil. To represent the general potential C, N, and P-acquiring enzyme activity, the C, N, and P acquiring enzymes were grouped and normalized as ln (β-1,4-glucosidase + β-1,4-xylosidase), ln (Urease + L-leucine aminopeptide), and ln (acid phosphatase), respectively (Sinsabaugh et al., 2009).

2.6 Extracting DNA from soil, Illumina MiSeq sequencing, and microbial co-occurrence network construction

Sample preparation and DNA extraction consist of the following steps: Soil specimens were flash-frozen in liquid nitrogen and stored at −80°C to prevent nucleic acid degradation. Total genome DNA from soil samples was extracted by the CTAB method (Kamdem et al., 2023). DNA concentration and purity were monitored on 1% (w/v) agarose gels. According to the concentration, DNA was diluted to l μg·μL−1 using sterile water. 16S rRNA/18S rRNA/ITS genes of distinct regions (16S V4/16S V3/16S V3-V4/16S V4-V5, 18S V4/18S V9, ITS1/ITS2, Arc V4) were amplified by specific primer 16S rRNA (515F: GTGCCAGCMGCCGCGGTAA and 907R: CCGTCAATTCCTTTGAGTTT) and ITS (ITS1-1F-F: CTTGGTCATTTAGAGGAAGTAA and ITS1-1F-R: GCTGCGTTCTTCATCGATGC) with the barcode. All PCR reactions were carried out with 15 μL of Phusion® High-Fidelity PCR Master Mix (New England Biolabs), 0.2 μM of forward and reverse primers, and about 10 ng template DNA. Thermal cycling consisted of initial denaturation at 98°C for 1 min, followed by 30 cycles of denaturation at 98°C for 10 s, annealing at 50°C for 30 s, and elongation at 72°C for 30 s, finally extension at 72°C for 5 min. Reactions were performed in triplicate under optimized cycling conditions (30–35 cycles) to minimize amplification bias. Amplicons were purified using AMPure XP beads and quantified via Agilent Bioanalyzer. The library was constructed using the TruSeq® DNA PCR-Free Sample Preparation Kit. The constructed libraries were quantified using Qubit and quantitative PCR (Q-PCR). Once the libraries were validated, they were sequenced on the NovaSeq 6000 platform. Finally, the libraries underwent sequencing on the Illumina NovaSeq platform. The raw sequence data reported in this paper have been deposited in the Genome Sequence Archive in the National Genomics Data Center, Beijing Institute of Genomics (China National Center for Bioinformation), Chinese Academy of Sciences, under accession project numbers CRA011413 and CRA011419 for 16S rRNA and ITS rRNA genes, respectively, which are publicly accessible at https://bigd.big.ac.cn/gsa (Edgar, 2004; Edgar, 2013; Quast et al., 2013). More details can be found in Appendix A. Supplementary material (Guimerà and Amaral, 2005; Tolosana-Delgado, 2008; Zhou et al., 2011; Deng et al., 2012; Luo et al., 2014; Ju and Zhang, 2015; Wu et al., 2021; Mao et al., 2022).

2.7 Grain yield and the economic value of faba bean

The yield of faba bean were measured at the maturity stage. The faba bean seeds were collected separately from each plot, and the weight of the seeds after the dry pod threshing was measured and recorded as the grain yield of the corresponding plot. The economic value was calculated based on the current year’s price of faba beans.

2.8 Data analysis and statistics

The data obtained in this study were maintained in Microsoft Excel 2016. Statistical analyses were performed using R software (Version 4.2.1). Data were subjected to Levene’s test for homogeneity of variance and the Shapiro–Wilk test for normality. One-way ANOVA analysis was performed to compare the effects of soil physicochemical properties, soil biological characteristics, and economic properties among different treatments (p = 0.05). Qiime software (Version 1.9.1) was employed to calculate Chao1, Shannon, Simpson, and Ace indices. The differences in the alpha diversity index between groups were further analyzed using R software. Two-way ANOVA and Bonferroni test were performed to examine the significant differences in soil aggregate stability, chemical composition, enzyme activity, and microbial diversity across different years and among diverse treatments (p = 0.05). Further details of data analyses can be found in Appendix A of Supplementary material.

3 Results

3.1 Impacts of organic mulch materials on the soil physicochemical properties and enzymatic activities

OMM treatment (S, M, and SM), as well as the year of plantation and their interactions, had significant effects on soil aggregate stability, nutrients content, enzymatic activity, and enzyme stoichiometric ratios (p < 0.05) (Supplementary Tables S3, S4; Figure 2 and Supplementary Figure S1). In 2022, compared to CK, soil aggregate stability of OMM treatment significantly improved (p < 0.05). In 2021 and 2022, compared to CK, the soil organic carbon (SOC), total nitrogen (TN) content, and the activities of soil C and N acquiring enzymes of OMM treatment increased. Soil macro-aggregate content (R0.25) and SOC increased with the increase in mulching duration. All the results mentioned above indicate that the parameters increased with SM treatment compared to other OMM treatments.

Figure 2
Nine-panel composite bar chart with data on enzyme and nutrient ratios for years 2021 and 2022. Panels A, B, and C show levels of C, N, and P acquiring enzymes, respectively. Panels D, E, and F display the C:N, C:P, and N:P enzyme ratios. Panels G, H, and I present the C:N, C:P, and N:P nutrient ratios. Each panel presents statistical results for treatment, year, and their interaction with annotations like significant (**) and non-significant (NS). The legend distinguishes data by year, with light brown for 2021 and green for 2022.

Figure 2. Enzyme activities, enzymatic stoichiometry, and soil nutrient stoichiometry under treatments using different organic mulching materials. (A) Activities of normalized C-acquiring enzymes. (B) Activities of normalized N-acquiring enzymes. (C) Activities of normalized P-acquiring enzymes. (D) C:N acquiring enzyme activity ratio. (E) C:P acquiring enzyme activity ratio. (F) N:P acquiring enzyme activity ratio. (G) C:N ratio. (H) C:P ratio. (I) N:P ratio. Different lowercase letters indicate significant differences between different treatments of the same year. Different capital letters indicate significant differences between different years of the same treatment. CK, without organic material mulching; S, straw mulching; M, milk vetch mulching; SM, straw and milk vetch mulching.

3.2 Effects of organic mulching materials on soil microbial community structure

In 2021 and 2022, the relative abundance of the dominant bacterial communities at the phylum level, i.e., Proteobacteria and Acidobacteria, among treatments demonstrated no significant difference (Figure 3A). However, the relative abundance of Actinobacteria and Acidobacteria in OMM treatment soils was higher than that in CK, i.e., an increase of 66.67 and 150%, respectively in 2022 (p < 0.05). Further, in 2021 and 2022, the relative abundance of dominant fungal communities at the phylum level, i.e., Ascomycota, Basidiomycota, Glomeromycota, and Mortierellomycota among treatments demonstrated no significant difference (Figure 3B).

Figure 3
(A) Stacked bar chart showing relative abundance of bacterial phyla over different samples. (B) Stacked bar chart showing relative abundance of fungal phyla over different samples. (C) CAP analysis of bacterial communities with variance explained and statistical significance. (D) CAP analysis of fungal communities with variance explained and statistical significance. (E) Circular phylogenetic tree of bacterial taxa for 2021. (F) Circular phylogenetic tree of bacterial taxa for 2022. (G) Circular phylogenetic tree of fungal taxa for 2021. (H) Circular phylogenetic tree of fungal taxa for 2022. Treatments include CK, S, M, and SM.

Figure 3. Differences in the microbial communities between CK and organic mulch treatments. Relative abundance of soil bacteria (A) and fungi (B) at the phylum level, CPCoA of soil bacteria (C) and fungi (D). Biomarkers of cultivated ginseng soil microbiota at different organic mulch treatments using the LEfSe method. Only results with |LDA| > 4 and p < 0.05 according to Tukey’s honest significant difference test are shown. (E,F) Bacterial biomarkers in the soil samples obtained from fields in 2021 and 2022, respectively. (G,H) Fungal biomarkers in the soil samples obtained from 2021 and 2022, respectively. Different lowercase letters indicate significant differences between different treatments of the same year. Different capital letters indicate significant differences between different years of the same treatment. CK, without organic material mulching; S, straw mulching; M, milk vetch mulching; SM, straw and milk vetch mulching.

The results of bacterial α-diversity indicated that in 2021, the observed species, Shannon index, and Simpson index in the SM treatment were significantly higher than those in the other treatments (Supplementary Table S5). In 2022, the observed species and Shannon indices in OMM treatment soils significantly increased by 7.7–18.85% and 6.05–8.48%, respectively (p < 0.05), when compared with CK. Moreover, the bacterial α diversity indices increased with the increase in mulching time. The fungal α diversity results demonstrated no significant difference among treatments in 2021 (Supplementary Table S6). However, in 2022, except for the Shannon index, other fungal α diversity indices in OMM treatment soils were higher than those in CK. The above results indicated greater improvement in the SM treatment when compared with other OMM treatments. The constrained principal coordinate analysis (CPCoA) of β diversity results for 2021 and 2022 (Figures 3C,D) demonstrated significant differences in the distribution of bacterial and fungal communities between CK and OMM treatments (p < 0.05), which corresponded to 57.4 and 44.1% of the total changes in the community, respectively.

The results of this study indicate that the predominant bacterial groups at the phylum level in CK were Gemmatimonadetes, Acidobacteria, and Planctomycete; and at the class level were Gemmatimonadetes, Subgroup 6, and Planctomycetacia (Figures 3E,F; Supplementary Figure S2A; Supplementary Table S9). The predominant bacterial groups at the phylum level in OMM treatment soil were Proteobacteria, Actinobacteria, and Bacteroidetes; and at the class level were Gammaproteobacteria, Actinobacteria, and Bacteroidia. In 2021 and 2022, the predominant fungal groups at the phylum level in CK were Ascomycota and Glomeromycota, and at the class level were Paraglomeromycetes and Saccharomycetes. Similarly, the predominant fungal groups in OMM treatment soil at the phylum level were Ascomycota and Basidiomycota, and at the class level were Sordariomycetes, Dothideomycetes, and Agaricomycetes. In addition, the major bacteria in the SM treatment soil at the class level were Anaerolineae and Verrucomicrobiae, and the major fungus in SM treatment soil at the class level was Mortierellomycetes (Figures 3G,H; Supplementary Figure S2B; Supplementary Table S9).

3.3 Effects of organic mulching materials on molecular ecological network characteristics of soil microbial communities

A comparative analysis of the molecular ecological network between CK and OMM treatments revealed improvements in the size and complexity of the microbial network under OMM treatments, particularly in terms of total nodes, total links, average degree, and total modularity number. In addition, the microbial molecular ecological network characteristics treated by OMM in 2022 were superior to those in 2021 (Figure 4; Supplementary Tables S7, S8).

Figure 4
Network graphs showing modules of bacterial and fungal communities in 2021 and 2022. Each graph represents modules with various colors, indicating different groupings. The top half displays bacteria modules for 2021 and 2022 with varying numbers, while the bottom half shows fungi modules for the same years. Labels such as CK, S, M, and SM indicate different conditions. A color legend identifies module numbers from zero to ten.

Figure 4. Co-occurrence networks of soil bacteria and fungi. (A–D) Bacteria in CK and S, M, and SM treatment fields in 2021. (E–H) Bacteria in CK and S, M, and SM treatment fields in 2022. (I–L) Fungi in CK and S, M, and SM treatment fields in 2021. (M–P) Fungi in CK and S, M, and SM treatment fields in 2022. CK, without organic material mulching; S, straw mulching; M, milk vetch mulching; SM, straw and milk vetch mulching.

The key soil microbial species in CK and OMM treatments were found in the connectors and module hubs of the MENs (Supplementary Figure S3; Supplementary Table S12) of soil microbial communities. In 2021 and 2022, the key bacterial species in CK at the phylum level were Gemmatimonadetes, Proteobacteria, and Acidobacteria, and at the class level were Gemmatimonadetes, Alphaproteobacteria, and Acidobacteriia. The key bacterial species in OMM treatment soils at the phylum level were Proteobacteria and Acidobacteria and at the class level were Alphaproteobacteria, Gammaproteobacteria, and Acidobacteriia. Further, in 2021 and 2022, the key fungus in CK at the phylum level was Ascomycota, and at the class level was Sordariomycetes. The key fungi in the OMM treatment soils at the phylum level were Ascomycota and Basidiomycota, and at the class level were Sordariomycetes, Eurotiomycetes, and Agaricomycetes. Specifically, the fungus in SM treatment at the class level was Glomeromycetes.

3.4 Relationship among soil environmental factors, microbial community structure, and soil microbial predicted functions

A significant linear relationship was observed between environmental factors and bacterial community structure in 2021 (p = 0.002; Supplementary Table S15). The redundancy analysis (RDA) of the bacterial community structure revealed that environmental factors accounted for 63.72% of the total variance. The C-acquiring enzyme activity (variation accounting for 15.48%, p = 0.002) and N-acquiring enzyme activity (variation accounting for 13.62%, p = 0.001) were the major environmental factors influencing the bacterial communities. Both these factors had a positive correlation with Actinobacteria. Additionally, the SM had a strong positive correlation with C-acquiring enzyme activity, N-acquiring enzyme activity, and Actinobacteria (Figure 5A). In 2022, the RDA of the bacterial community structure revealed that the environmental factors accounted for 72.94% of the total variance. C-acquiring enzyme activity (variation accounting for 32.9%, p = 0.121) and SOC content (variation accounting for 14.9%, p = 0.311) were the major environmental factors influencing the bacterial communities. These factors were positively correlated with Actinobacteria, Acidobacteria, and Bacteroidetes. OMM had a strong positive correlation with C acquiring enzyme activity, SOC content, Actinobacteria, Acidobacteria, and Bacteroidetes (Figure 5B).

Figure 5
Four redundancy analysis plots (A-D) displaying changes in bacterial and fungal communities in 2021 and 2022. Panels (A) and (B) show bacteria, while (C) and (D) illustrate fungi. Symbols represent different treatments: CK (red triangles), S (blue triangles), M (green triangles), and SM (purple triangles). Arrows indicate factors like total nitrogen (TN), total potassium (TK), soil organic carbon (SOC), and various nutrient-acquiring enzymes interacting with microbial groups such as Proteobacteria, Acidobacteria, Ascomycota, and Glomeromycota. Each plot details changes in community composition based on RDA1 and RDA2 axes percentages.

Figure 5. Redundancy analysis (RDA) of the soil microbial community structure and physicochemical factors in the fields. (A) Bacteria in 2021, (B) Bacteria in 2022, (C) Fungi in 2021, (D) Fungi in 2022. CK, without organic material mulching; S, straw mulching; M, milk vetch mulching; SM, straw and milk vetch mulching.

In 2021 and 2022, a significant linear relationship was observed between the environmental factors and fungal community structure. The RDA of the fungal community structure demonstrated that the environmental factors accounted for 37.2% of the total variance in 2021. The C-acquiring enzyme activity (variation accounting for 10%, p = 0.059) and TN content (variation accounting for 10.7%, p = 0.085) were the major environmental factors influencing the fungal communities. Both these factors were positively correlated with Ascomycota (Figure 5C). In 2022, the RDA of the fungal community structure demonstrated that the environmental factors accounted for 47.82% of the total variance. Among them, TK content and C-acquiring enzyme activity were the major environmental factors influencing the fungal communities and were positively correlated with Basidiomycota. OMM had a strong positive correlation with TK content, C acquiring enzyme activity, SOC content, TN content, and Basidiomycota (Figure 5D).

In 2021 and 2022, when compared with CK, PICRUSt predicted that OMM significantly up-regulated 5 first-order biological metabolic pathways (Level 1) and 17 s-order biological metabolic pathways (Level 2) (p < 0.05). The 5 Level 1 pathways included genetic information processing (Level 2: genetic information processing unclassified, folding, sorting, degradation, translation and replication, and repair), metabolism (Level 2: glycan biosynthesis and metabolism, enzyme families, and energy metabolism), cellular processes (Level 2: cell motility and cell communication), environmental information processing (Level 2: signal transduction, metabolism of cofactors and vitamins, and nucleotide metabolism), and organismal systems (Level 2: digestive system, immune system, environmental adaptation, and sensory system) (p < 0.05) (Supplementary Table S13).

In 2021 and 2022, when compared with CK, FUNGuild predicted that OMM significantly increased the two trophic modes and 19 functional groups, and both saprotrophs and symbiotrophs significantly increased in the trophic modes (p < 0.05). Fungal parasite, wood saprotroph, undefined saprotroph, nematophagous, plant saprotroph, endophyte, litter saprotroph, soil saprotroph, and epiphyte were significantly up-regulated in the functional groups (p < 0.05) (Supplementary Table S14).

3.5 Mechanism underlying organic mulching materials for enhanced faba bean productivity

In 2021, the grain yield and economic value of the faba bean in the OMM treatment field increased by 40.94–73.87% when compared with CK field and the increments observed in the case of S and SM treatment fields were highly significant (p < 0.05). In 2022, except for the S treatment field, the grain yield and economic value of the faba bean in the OMM treatment field increased by 9.95–12.99% when compared with CK field. The grain yield and economic value in M treatment fields increased with the increase in mulching time (Figures 6A,B).

Figure 6
Bar graphs and path models show the effects of different treatments on faba bean productivity. (A) and (B) depict grain yield and economic value across treatments CK, S, M, and SM for 2021 and 2022. (C) and (D) illustrate the relationships between soil physical structure, nutrient levels, enzyme activity, and microbial diversity, highlighting the impact of organic material mulching versus no mulching on faba bean productivity, marked by positive and negative effects. Statistical measures are included with significance levels.

Figure 6. Faba bean productivity under CK and organic materials mulching treatments (A,B) and partial least squares path model (PLS-PM) illustrating the direct effects of (C) CK on the soil quality and grain yield of faba bean and (D) OMM (S, M, and SM) on the soil quality and grain yield of faba bean. Numbers on the arrowed lines and the thickness of the arrows indicate normalized path coefficients. Different lowercase letters indicate significant differences between different treatments of the same year. Different capital letters indicate significant differences between different years of the same treatment. CK, without organic material mulching; S, straw mulching; M, milk vetch mulching; SM, straw and milk vetch mulching.

The RFM analysis was performed using 17 abiotic and microbial variables (Supplementary Table S16). The results showed that no index had a significant effect on the grain yield of faba bean under the CK system. However, the changes in GMD, SOC, and Shannon index induced by OMM treatments had a significant positive effect on the grain yield of the faba bean (p < 0.05).

The goodness-of-fit (GOF) of the partial least squares path model (PLS-PM) for soil variables (abiotic and microbial) and faba bean grain yield was higher than 0.8 in CK as well as OMM treatments (Figures 6C,D), which confirms the reliability of the model. The PLS-PM indicated that predictor variables accounted for 68–89%, 67–69%, 57–82%, 82–90%, and 71–89% of the variations in the soil physical structure, soil nutrients, soil enzyme activity, soil microbial diversity, and faba bean productivity, respectively.

CK treatment significantly reduced the stability of soil physical structure (path coefficient = 0.89) and soil enzyme activity (path coefficient = 0.82), significantly increasing the grain yield of faba bean (path coefficient = 0.71) (p < 0.05). Further, soil enzyme activity significantly enhanced soil microbial diversity (path coefficient = 0.85), but significantly reduced the grain yield of faba bean (path coefficient = 0.75) (p < 0.05). The soil nutrients (path coefficient = 0.65) and soil microbial diversity (path coefficient = 0.63) had significant impacts on the productivity of the faba bean (p < 0.05) (Figure 6C).

OMM treatments significantly increased soil nutrients (path coefficient = 0.79), soil enzyme activity (path coefficient = 0.85), soil microbial diversity (path coefficient = 0.91), and faba bean grain yield (path coefficient = 1.42) (p < 0.05). Further, the increase in soil microbial diversity was associated with a significant increase in soil nutrients (path coefficient = 0.73) and faba bean productivity (path coefficient = 1.06) (p < 0.05).

4 Discussion

4.1 Mechanism underlying organic mulching materials for enhanced soil aggregates stability and carbon cycling efficiency of faba bean

The results of this study demonstrated that in 2021 and 2022, OMM treatments enhanced the SOC content and activities of the C acquisition enzyme when compared with CK. In 2022, OMM treatments significantly improved soil aggregate stability when compared with CK (p < 0.05). The stability of soil aggregates and SOC content increase with the increase in mulching time. These results indicated that the improvement in soil aggregates and C cycle by OMM treatments was a result of long-term accumulation. Returning organic materials to the field enhances soil aggregate stability and C cycling efficiency through interactions among aggregates, SOC content, C-acquisition enzymes, and microbial communities. Current research demonstrates that enhanced aggregate stability arises not only from organic-mineral coordination forming stable complexes, but also involves organic-induced adsorption of microbially-derived extracellular polymers (e.g., polysaccharides, proteins, humic substances), fungal hyphae, and colloidal materials, ultimately fostering microbial residue formation (Kan et al., 2022; Ai S. et al., 2023). Wang Y. et al. (2023) conducted a study for five consecutive years and reported that straw mulching can significantly improve the distribution and stability of soil aggregates. Kamran et al. (2021) performed long-term field experiments for 10 consecutive years and demonstrated that fresh milk vetch could significantly increase the distribution of large soil aggregates >0.25 mm.

The major reason for the improvement in soil aggregate stability through the return of organic materials to the field is the synergistic effect of abiotic and biotic activities in the soil, which is closely linked to the C cycle (Huang et al., 2017). In terms of abiotic components, SOC is a major contributor to several decomposition products that promote soil microbial activity, including carboxylic acids, alcohols, phenols, polysaccharides, and amines. These organic compounds bind with soil minerals and agglomerate to form microaggregates (Liu Y. et al., 2023). In this study, when compared with CK, OMM treatments effectively increased the SOC content (Supplementary Table S4). The RDA analysis demonstrated a positive correlation of OMM with SOC (Figures 5B,D). Numerous research studies have demonstrated that the return of straw and legume green fertilizer to the field has a significant influence on SOC and C mineralization (Hobbie and Hobbie, 2013; Dhaliwal et al., 2020). Soil biological activities manifest as variations in enzyme activity and microbial community structure, primarily derived from microorganisms, plants, animals, and organic residues (Sun et al., 2023). This study confirmed that the OMM treatments can promote the activity of key enzymes involved in the soil C cycle (Figure 5). Previous studies have demonstrated that the return of straw and legume green fertilizer to the field had a positive effect on enzyme activity associated with soil C cycle (Khan et al., 2020; Liu et al., 2022). Bowles et al. (2014) demonstrated that the application of legume green fertilizer can increase the activity of β-glucosidase and α-galactosidase, which are considered as key enzymes of the soil C cycle. Further, this study demonstrated the positive effect of OMM treatments on soil microbial diversity. Increased microbial diversity drives soil aggregation through mycelia production and polysaccharide secretion, enhancing aggregate formation and stabilization (Ren et al., 2018; Tsutsui et al., 2022). The microbial residues contain a large number of amino sugars, which are also important sources of SOC (Wang et al., 2021). In addition, the results of this study demonstrated a significant improvement in the relative abundance of Bacteroidetes and Actinobacteria by OMM treatments (Figure 3A). The metabolic energy of Bacteroidetes is derived from starch, cellulose, and other polysaccharide substances (Zheng et al., 2021). Bacteroidetes, as key decomposers mediating carbohydrate fermentation and nitrogenous organic compound utilization, drive straw decomposition processes (Tlaskal and Baldrian, 2021). Actinobacteria are capable of degrading cellulose and chitin, which are the major sources of soil nutrients (Ren et al., 2022). Further, this study revealed that the dominant bacterial groups in OMM treatment soil were Alphaproteobacteria and Gammaproteobacteria (Supplementary Figure S3; Supplementary Table S12). These two groups of bacteria can have a significant role in straw degradation and promote the accumulation of organic carbon (Ma et al., 2024). The RDA analysis confirmed a positive correlation of Bacteroidetes, Acidobacteria, Ascomycota, and Basidiomycota with C-acquiring enzyme activity and SOC (Figure 5).

4.2 Mechanism underlying organic mulching materials for enhanced soil nitrogen cycling efficiency of faba bean

OMM treatments elevated soil nitrogen content and N-acquisition enzyme activity relative to CK during 2021–2022, demonstrating enhanced coupling of soil carbon and nitrogen cycling. Straw serves as a key organic substrate providing C and energy for soil microorganisms, while enhancing the mineralizable organic N: available C ratio to stimulate soil N mineralization and nitrification (Chen et al., 2017). Several research studies have demonstrated that soil N mineralization and nitrification can increase with the increase in SOC content (Yu H. et al., 2022). This study revealed that the SOC, TN, and AN content in the S treatment field were higher than those in CK throughout the experiment in 2021 and 2022. The RDA confirmed a positive correlation between C and N contents and their enzyme activities, which revealed that SOC input could promote N cycling in dry farmlands to a certain extent (Supplementary Table S4; Figure 5). M treatment enhanced soil nitrogen cycling via milk vetch incorporation, a leguminous source supplying N–P–K-micronutrients and mediating symbiotic N fixation for nutrient enrichment (Ma et al., 2020). The results of this study demonstrated that S treatment could effectively improve soil N nutrient and N acquiring enzyme activity when compared with CK. In addition, RDA also confirmed a strong positive correlation of S with TN content and N-acquiring enzyme activity (Supplementary Table S4; Figure 5).

However, in this study, SM treatment outperformed other treatments in improving the efficiency of C and N cycling. Fresh straw decomposition depletes soil available nitrogen via microbial mineralization, impeding overall organic materials decomposition through induced N limitation (Hicks et al., 2021). Although the simple application of milk vetch is beneficial for N storage, it is not beneficial for the long-term enhancement of SOC (Zhong et al., 2021). Research studies have demonstrated that organic materials with a low C:N ratio are conducive to the decomposition of straw C by soil microorganisms. The decomposition products bind with soil minerals and contribute to the formation of aggregates (Zhang et al., 2017). SM treatment stimulated SOC accumulation through enhanced fresh straw decomposition and alleviated microbial mineral N competition via legume application (Xie et al., 2022a). These results provided evidence that SM treatment enhances SOC, TN, AN, and the activities of C and N acquiring enzymes in the soil (Supplementary Table S4; Figure 2). Straw-milk vetch co-incorporation stimulated microbial community diversification, yielding abundant polysaccharides that drive aggregate cementation through metabolic processes and enhance macroaggregate formation for C–N adsorption sites (Kamran et al., 2021). The findings of this study demonstrated that the major microbial class in the SM-enriched soil was Anaerolineae (Supplementary Table S6), a group of Gram-negative bacteria that can effectively degrade humus in the soil (Zhang et al., 2025). Research studies have demonstrated that Anaerolineae can reduce N fertilizer application under the influence of straw returning to the field, thus enriching the soil with this bacterial class and increasing the accumulation of nutrients (Zhang et al., 2019; Podosokorskaya et al., 2013).

4.3 Regulation of core microbial groups by organic mulching materials

The soil microbial community is the major component of the soil that mediates the cycling of important nutrients such as C, N, and P through anabolic and catabolic processes (Liu S. et al., 2023). The results of this study indicated that the OMM treatments can increase the structural and functional diversity of the soil microbial communities. The addition of organic materials can provide necessary nutrients and energy for microbial activities, and also increase the diversity of soil bacteria and fungi (Adetunji et al., 2020; Ren et al., 2021). Long-term straw return stimulates soil bacterial/fungal proliferation and microbial community richness, with microbial-derived amino sugars enhancing SOC accumulation and aggregate formation (Li, 2022), which is consistent with the results of this study. Straw with low density and light texture can increase the accumulated temperature and water retention capacity of the soil (Feng et al., 2021). Straw-milk vetch co-incorporation enhances macroaggregate formation and stability by supplying microbial membrane C reservoirs and metabolic energy, while fostering microbial colonization through habitat optimization and species richness elevation (Asghar and Kataoka, 2022; Yuan et al., 2023). These studies explain the enrichment of soil microorganisms through OMM treatments in this research. Kim et al. (2020) reported that OMM is a sustainable method to enhance agricultural production and has a positive impact on soil microbial metabolism. PICRUSt and FUNGuild analyses revealed OMM treatment possible upregulation of functional pathways (genetic information processing, metabolism, organismal systems) enhancing soil metabolic potential, while elevated proportion of saprotrophs and symbiotrophs may be drove soil C–N cycling efficiency through synergistic functional contributions (Supplementary Tables S13, S14). These findings are consistent with the previous studies demonstrating that organic materials return to the field can enhance the soil microbial community structure, promote co-trophic bacteria, and contribute to soil sustainability (Zhao et al., 2017; Price et al., 2021). In addition, this study demonstrated that the positive impact of SM treatment on soil microbial diversity was significantly higher than that of S and M treatments. This may be due to the higher P acquiring enzyme activity and soil P content in SM treatment, and the lower C:P acquiring enzyme activity ratio and N:P acquiring enzyme activity ratio. The greater soil P metabolic efficiency resulting from SM treatment may increase P availability for microbial energy production and enhance the microbial community structure (Lv et al., 2022).

Microbial MENs construction elucidates soil microbial community architecture and interaction dynamics through topological feature analysis (e.g., node/link influence), identifying keystone taxa driving community structure–function relationships (Ai Z. M. et al., 2023). This study demonstrated that OMM treatments enhanced the topological metrics of microbial MENs, including total nodes, connections, average degree, and modularity (Supplementary Tables S7, S8). These findings indicate an increase in the complexity and stability of microbial system architecture. Straw return to field enhanced soil co-occurrence network complexity in Xu et al. (2023), manifested through elevated connectivity, keystone taxa abundance, and topological metrics (average degree, clustering coefficient) reinforcing microbial network stability. In addition, the dominant species and microbial groups often play a vital role in maintaining the stability of microbial community structure and function (Marshall and Lynch, 2020). In this study, OMM treatments led to the establishment of a microbial community with Actinobacteria, Bacteroidetes, Ascomycota, and Basidiomycota as the core microbial phyla. These microorganisms play an active role in soil nutrient cycling and maintain the large energy consumption of the microbial community. It is worth mentioning that the dominant microbes of CK and OMM treatment soils belong to the phylum Ascomycota, which plays a vital role in the decomposition of soil organic materials. However, it requires a large amount of N as an energy source (Lopez-Mondejar et al., 2018). OMM can provide more N to soil compared with CK, which is one of the main reasons for the complexity of soil microbial MENs following the application of OMM.

4.4 Enhanced faba bean productivity by organic mulching materials through the promotion of soil microbial diversity and carbon cycling

The results of this study demonstrated that the OMM treatments could improve the yield and economic value of faba bean, which is consistent with the conclusion that straw mulching and purple milk vetch return to the field could improve crop yield (Wang M. et al., 2023). Since the variations of yield and economic value are the same, this study analyzed the effects of OMM on the faba bean yield alone when compared to CK using RFM analysis. The analysis results indicated that GMD, SOC, and Shannon index had significant effects on the faba bean grain yield under OMM treatments. There was a gradual decrease in the influence of SOC, C acquiring enzyme activity, N acquiring enzyme activity, and Shannon and Simpson indices on the grain yield of faba bean. These results indicated that SOC and microbial diversity were the major driving factors for the enhanced faba bean grain yield under the OMM treatments. Several studies have demonstrated that straw mulching enhanced soil active C pool and SOC content, thereby providing favorable conditions for growth and development of crops (Liu D. et al., 2023; Wang M. et al., 2023). Gai et al. (2019) reported that straw return to the field could improve crop yield by augmenting the direct nutrient supply for increasing soil fertility. Soil biodiversity elevation stimulates beneficial microorganism proliferation and interpopulation coordination, thereby improving soil quality and environmental stress resilience as demonstrated by Stefan et al. (2021). Soil microbial diversity stimulates crop yield through functional enhancement, with meta-analysis by Shu et al. (2022) demonstrating significant positive correlations among diversity, community structure, functionality, and productivity. In this study, PLS-PM demonstrated that OMM treatments significantly improved the soil nutrient content, soil microbial diversity, and faba bean productivity. OMM treatments enhanced soil nutrient availability and faba bean yield via soil microbial diversity-driven mechanisms, establishing microbial diversity as the primary driver under such amendments. The yield of S and SM treatments decreased during 2021–2022, with S exhibiting significant interannual variation contrasting SM stability. This divergence likely originated from S-mediated straw return without N supplementation, elevating soil organic carbon but inducing long-term microbial-crop N competition that reduced yields. Conversely, SM with N supplementation mitigated such competition. Strategic integration of optimized N inputs with organic-mineral fertilizers is recommended to maximize ecological benefits (Li et al., 2024; Yu et al., 2024).

Corn straw mulching revealed superior efficacy in enhancing soil quality and crop yield versus granulator-produced pellet incorporation, offering novel strategies for optimizing straw incorporation practices as demonstrated by Fan and Wu (2020). No-tillage combined with straw returning during conservation tillage plays an important role in promoting soil aggregate stability, soil C and N cycling, and microbial diversity (Guo et al., 2016). Further studies can be conducted to explore the ecological benefits of OMM treatments combined with no-tillage mode on farmland soils and crops. Future OMM treatment research should prioritize multi-omics integration for key C/N metabolic gene mapping and soil metabolome profiling, addressing current methodological gaps. These studies will be essential for a systematic understanding of the potential role of OMM in enhancing the ecological benefits of soil and crops in farmlands.

5 Conclusion

In this study, organic mulching materials (OMM) significantly enhanced soil macroaggregate stability and stimulated carbon (C) and nitrogen (N) accrual through elevated enzymatic activities over time. OMM treatment induced microbial community diversification with complex co-occurrence networks. The dominant phyla (Actinobacteria, Bacteroidetes, Ascomycota, Basidiomycota) exhibited strong positive correlations with C and N contents, as well as with nutrient-acquiring enzyme activities. The functional prediction analysis conducted using PICRUSt and FUNGuild suggests that the enhancement of soil metabolic capacity may be related to the processing, metabolism, and improvement of biological system pathways associated with microbial functional gene information. Furthermore, the increase in the proportion of saprotrophs and symbiotrophs may have improved the efficiency of soil C and N cycling. These microbial-driven mechanisms, particularly diversity indices followed by soil organic content (SOC) and C-acquiring enzymes activities, directly enhanced faba bean productivity and economic value. Our findings establish straw-milk vetch co-mulching as a sustainable strategy for augmenting agroecosystem services through simultaneous soil health improvement and yield optimization.

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

KR: Formal analysis, Investigation, Software, Writing – original draft. WS: Formal analysis, Investigation, Methodology, Writing – review & editing. ZW: Data curation, Writing – review & editing. LS: Formal analysis, Visualization, Writing – review & editing. ML: Methodology, Writing – review & editing. YulZ: Methodology, Writing – review & editing. YuzZ: Project administration, Resources, Writing – review & editing. XW: Formal analysis, Supervision, Validation, Writing – review & editing. KG: Conceptualization, Data curation, Writing – review & editing. KS: Investigation, Writing – review & editing. LW: Conceptualization, Funding acquisition, Supervision, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was supported by the Natural Science Foundation Project of China (No. 31871583).

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.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

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.2025.1602633/full#supplementary-material

References

Adetunji, A. T., Ncube, B., Mulidzi, R., and Lewu, F. B. (2020). Management impact and benefit of cover crops on soil quality: a review. Soil Tillage Res. 204:104717. doi: 10.1016/j.still.2020.104717

Crossref Full Text | Google Scholar

Ahmed, K., Mohamed, B., and Cherki, G. (2020). Effect of soil available phosphorus levels on chickpea (Cicer arietinum L.) - rhizobia symbiotic association. Legum. Res. 43, 878–883. doi: 10.18805/lr-524

Crossref Full Text | Google Scholar

Ai, Z. M., Deng, Y., Li, X. H., Zhang, J. Y., Liu, H. F., Xu, H. W., et al. (2023). Effect of plant-soil feedback on soil microbial co-occurrence network depends on the stage of secondary succession. Rhizosphere 27:11. doi: 10.1016/j.rhisph.2023.100733

Crossref Full Text | Google Scholar

Ai, S., Meng, X., Zhang, Z., Li, R., Teng, W., Cheng, K., et al. (2023). Artificial humic acid regulates the impact of fungal community on soil macroaggregates formation. Chemosphere 332:138822. doi: 10.1016/j.chemosphere.2023.138822

PubMed Abstract | Crossref Full Text | Google Scholar

Asghar, W., and Kataoka, R. (2022). Green manure incorporation accelerates enzyme activity, plant growth, and changes in the fungal community of soil. Arch. Microbiol. 204:7. doi: 10.1007/s00203-021-02614-x

PubMed Abstract | Crossref Full Text | Google Scholar

Bao, S. D. (1999). Soil agricultural chemistry analysis. Beijing, China: Agriculture Press.

Google Scholar

Barq, M. G., Hassan, M. M., Yasmin, H., Shahzad, A., Malik, N. H., Lorenz, N., et al. (2021). Variation in archaeal and bacterial community profiles and their functional metabolic predictions under the influence of pure and mixed fertilizers in paddy soil. Saudi J. Biol. Sci. 28, 6077–6085. doi: 10.1016/j.sjbs.2021.08.091

PubMed Abstract | Crossref Full Text | Google Scholar

Bhaduri, D., Sihi, D., Bhowmik, A., Verma, B. C., Munda, S., and Dari, B. (2022). A review on effective soil health bio-indicators for ecosystem restoration and sustainability. Front. Microbiol. 13:25. doi: 10.3389/fmicb.2022.938481

Crossref Full Text | Google Scholar

Bowles, T. M., Acosta-Martínez, V., Calderón, F., and Jackson, L. E. (2014). Soil enzyme activities, microbial communities, and carbon and nitrogen availability in organic agroecosystems across an intensively-managed agricultural landscape. Soil Biol. Biochem. 68, 252–262. doi: 10.1016/j.soilbio.2013.10.004

Crossref Full Text | Google Scholar

Chen, S., Ding, X. Q., Zhu, Z. K., Wang, J., Peng, P. Q., Ge, T. d., et al. (2017). Effect of straw application on the dynamics of exogenous nitrogen and microbial activity in paddy soil. Environ. Sci. 38, 1613–1621. doi: 10.13227/j.hjkx.201609219

Crossref Full Text | Google Scholar

Chen, W., Li, P., Li, F., Xi, J., and Han, Y. (2022). Effects of tillage and biochar on soil physiochemical and microbial properties and its linkage with crop yield. Front. Microbiol. 13:929725. doi: 10.3389/fmicb.2022.929725

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, L. M., Sun, S. L., Yao, B., Peng, Y. T., Gao, C. F., Qin, T., et al. (2022). Effects of straw return and straw biochar on soil properties and crop growth: a review. Front. Plant Sci. 13:12. doi: 10.3389/fpls.2022.986763

Crossref Full Text | Google Scholar

Chen, S. H., Xiang, X. L., Ma, H. L., Penttinen, P., Zheng, T., Huang, X. L., et al. (2022). Response of soil bacterial communities in wheat rhizosphere to straw mulching and N fertilization. Front. Microbiol. 13:11. doi: 10.3389/fmicb.2022.982109

Crossref Full Text | Google Scholar

Chu, L., Sun, T., Wang, T., Li, Z., and Cai, C. (2020). Temporal and spatial heterogeneity of soil erosion and a quantitative analysis of its determinants in the three gorges reservoir area, China. Int. J. Environ. Res. Public Health 17:8486. doi: 10.3390/ijerph17228486

PubMed Abstract | Crossref Full Text | Google Scholar

Deng, Y., Jiang, Y.-H., Yang, Y., He, Z., Luo, F., and Zhou, J. (2012). Molecular ecological network analyses. BMC Bioinf. 13:113. doi: 10.1186/1471-2105-13-113

PubMed Abstract | Crossref Full Text | Google Scholar

Dhaliwal, S. S., Naresh, R. K., Gupta, R. K., Panwar, A. S., Mahajan, N. C., Singh, R., et al. (2020). Effect of tillage and straw return on carbon footprints, soil organic carbon fractions and soil microbial community in different textured soils under rice-wheat rotation: a review. Rev. Environ. Sci. Biotechnol. 19, 103–115. doi: 10.1007/s11157-019-09520-1

Crossref Full Text | Google Scholar

Dou, Y., Wen, M., Yang, C., Zhao, F., Ren, C., Zhang, N., et al. (2023). Effects of straw and plastic film mulching on microbial functional genes involved in soil nitrogen cycling. Front. Microbiol. 14:1205088. doi: 10.3389/fmicb.2023.1205088

PubMed Abstract | Crossref Full Text | Google Scholar

Edgar, R. C. (2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797. doi: 10.1093/nar/gkh340

PubMed Abstract | Crossref Full Text | Google Scholar

Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604

PubMed Abstract | Crossref Full Text | Google Scholar

Fan, W., and Wu, J. (2020). Short-term effects of returning granulated straw on soil microbial community and organic carbon fractions in dryland farming. J. Microbiol. 58, 657–667. doi: 10.1007/s12275-020-9266-5

PubMed Abstract | Crossref Full Text | Google Scholar

Fang, C., Li, P., Zhang, J., Lu, Y., Tang, Y., Tu, N., et al. (2024). Soil cd bioavailability response characteristics to microbes in paddy fields with co-incorporation of milk vetch, rice straw and amendments. Sci. Total Environ. 935:173306. doi: 10.1016/j.scitotenv.2024.173306

PubMed Abstract | Crossref Full Text | Google Scholar

FAOSTAT (2021). Crops and livestock products. Rome, Italy: Food and Agriculture Organization of the United Nations (FAO).

Google Scholar

Feng, W., Yang, F., Cen, R., Liu, J., Qu, Z., Miao, Q., et al. (2021). Effects of straw biochar application on soil temperature, available nitrogen and growth of corn. J. Environ. Manag. 277:111331. doi: 10.1016/j.jenvman.2020.111331

PubMed Abstract | Crossref Full Text | Google Scholar

Gai, X., Liu, H., Liu, J., Zhai, L., Wang, H., Yang, B., et al. (2019). Contrasting impacts of long-term application of manure and crop straw on residual nitrate-N along the soil profile in the North China plain. Sci. Total Environ. 650, 2251–2259. doi: 10.1016/j.scitotenv.2018.09.275

PubMed Abstract | Crossref Full Text | Google Scholar

Gao, S., Cao, W., Zhou, G., and Rees, R. M. (2021). Bacterial communities in paddy soils changed by milk vetch as green manure: a study conducted across six provinces in South China. Pedosphere 31, 521–530. doi: 10.1016/s1002-0160(21)60002-4

Crossref Full Text | Google Scholar

Grobelak, A., Kokot, P., Hutchison, D., Grosser, A., and Kacprzak, M. (2018). Plant growth-promoting rhizobacteria as an alternative to mineral fertilizers in assisted bioremediation - sustainable land and waste management. J. Environ. Manag. 227, 1–9. doi: 10.1016/j.jenvman.2018.08.075

PubMed Abstract | Crossref Full Text | Google Scholar

Guimerà, R., and Amaral, L. A. N. (2005). Cartography of complex networks: modules and universal roles. J. Stat. Mech. P02001, 1–13. doi: 10.1088/1742-5468/2005/02/p02001

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, L. J., Lin, S., Liu, T. Q., Cao, C. G., and Li, C. F. (2016). Effects of conservation tillage on topsoil microbial metabolic characteristics and organic carbon within aggregates under a rice (Oryza sativa L.) –wheat (Triticum aestivum L.) cropping system in Central China. PLoS One 11:e0146145. doi: 10.1371/journal.pone.0146145

PubMed Abstract | Crossref Full Text | Google Scholar

Hassan, M. U., Guoqin, H., Arif, M. S., Mubarik, M. S., Tang, H., Xu, H., et al. (2024). Can urea-coated fertilizers be an effective means of reducing greenhouse gas emissions and improving crop productivity? J. Environ. Manag. 367:121927. doi: 10.1016/j.jenvman.2024.121927

PubMed Abstract | Crossref Full Text | Google Scholar

He, C., Zhou, B., Wang, H., Wei, Y., and Huang, J. (2023). A first-year maize/cassava relay intercropping system improves soil nutrients and changes the soil microbial community in the symbiotic period. Front. Microbiol. 14:1087202. doi: 10.3389/fmicb.2023.1087202

PubMed Abstract | Crossref Full Text | Google Scholar

Hicks, L. C., Lajtha, K., and Rousk, J. (2021). Nutrient limitation may induce microbial mining for resources from persistent soil organic materials. Ecology 102:e03328. doi: 10.1002/ecy.3328

PubMed Abstract | Crossref Full Text | Google Scholar

Hobbie, J. E., and Hobbie, E. A. (2013). Microbes in nature are limited by carbon and energy: the starving-survival lifestyle in soil and consequences for estimating microbial rates. Front. Microbiol. 4:324. doi: 10.3389/fmicb.2013.00324

PubMed Abstract | Crossref Full Text | Google Scholar

Huang, R., Lan, M., 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. Res. 24, 27942–27952. doi: 10.1007/s11356-017-0372-9

PubMed Abstract | Crossref Full Text | Google Scholar

Jithesh, T., James, E. K., Iannetta, P. P. M., Howard, B., Dickin, E., and Monaghan, J. M. (2024). Recent progress and potential future directions to enhance biological nitrogen fixation in faba bean (Vicia faba L.). Plant-Environ. Interact. 5:e10145. doi: 10.1002/pei3.10145

PubMed Abstract | Crossref Full Text | Google Scholar

Ju, F., and Zhang, T. (2015). Bacterial assembly and temporal dynamics in activated sludge of a full-scale municipal wastewater treatment plant. ISME J. 9, 683–695. doi: 10.1038/ismej.2014.162

PubMed Abstract | Crossref Full Text | Google Scholar

Kamdem, C. N., Fogue, P. S., Tiofack, A. A. Z., Mewamba, E. M., Womeni, H. M., Koffi, M., et al. (2023). Assessment of cetyl-trimethyl-ammonium bromide (CTAB) based method for the extraction of soil-transmitted helminth dnas from stools for molecular dagnostic of soil-transmitted helminthiasis. J. Microbiol. Methods 204:106661. doi: 10.1016/j.mimet.2022.106661

PubMed Abstract | Crossref Full Text | Google Scholar

Kamran, M., Huang, L., Nie, J., Geng, M., Lu, Y., Liao, Y., et al. (2021). Effect of reduced mineral fertilization (NPK) combined with green manure on aggregate stability and soil organic carbon fractions in a fluvo-aquic paddy soil. Soil Tillage Res. 211:105005. doi: 10.1016/j.still.2021.105005

Crossref Full Text | Google Scholar

Kan, Z. R., Liu, W. X., Liu, W. S., Lal, R., Dang, Y. P., Zhao, X., et al. (2022). Mechanisms of soil organic carbon stability and its response to no-till: a global synthesis and perspective. Glob. Change Biol. 28, 693–710. doi: 10.1111/gcb.15968

PubMed Abstract | Crossref Full Text | Google Scholar

Karkanis, A., Ntatsi, G., Lepse, L., Fernandez, J. A., Vagen, I. M., Rewald, B., et al. (2018). Faba bean cultivation - revealing novel managing practices for more sustainable and competitive European cropping systems. Front. Plant Sci. 9:1115. doi: 10.3389/fpls.2018.01115

PubMed Abstract | Crossref Full Text | Google Scholar

Khan, M. I., Gwon, H. S., Alam, M. A., Song, H. J., Das, S., and Kim, P. J. (2020). Short term effects of different green manure amendments on the composition of main microbial groups and microbial activity of a submerged rice cropping system. Appl. Soil Ecol. 147:103400. doi: 10.1016/j.apsoil.2019.103400

Crossref Full Text | Google Scholar

Kim, N., Zabaloy, M. C., Guan, K., and Villamil, M. B. (2020). Do cover crops benefit soil microbiome? A meta-analysis of current research. Soil Biol. Biochem. 142:107701. doi: 10.1016/j.soilbio.2019.107701

Crossref Full Text | Google Scholar

Li, J. (2022). Effects and the microbial mechanism of straw deep-burying to soil nitrogen use. [Dissertation]. Shenyang: Shenyang Agricultural University.

Google Scholar

Li, P., Jia, L., Chen, Q., Zhang, H., Deng, J., Lu, J., et al. (2024). Adaptive evaluation for agricultural sustainability of different fertilizer management options for a green manure-maize rotation system: impacts on crop yield, soil biochemical properties and organic carbon fractions. Sci. Total Environ. 908:168170. doi: 10.1016/j.scitotenv.2023.168170

PubMed Abstract | Crossref Full Text | Google Scholar

Li, L., Konkel, J., Jin, V. L., and Schaeffer, S. M. (2021). Conservation management improves agroecosystem function and resilience of soil nitrogen cycling in response to seasonal changes in climate. Sci. Total Environ. 779:146457. doi: 10.1016/j.scitotenv.2021.146457

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, D., Song, C., Xin, Z., Fang, C., Liu, Z., and Xu, Y. (2023). Agricultural management strategies for balancing yield increase, carbon sequestration, and emission reduction after straw return for three major grain crops in China: a meta-analysis. J. Environ. Manag. 340:117965. doi: 10.1016/j.jenvman.2023.117965

Crossref Full Text | Google Scholar

Liu, Y., Wang, P., and Wang, J. (2023). Formation and stability mechanism of soil aggregates: progress and prospect. Acta Pedol. Sin. 60, 627–643. doi: 10.11766/trxb202112180686

Crossref Full Text | Google Scholar

Liu, B., Xia, H., Jiang, C., Riaz, M., Yang, L., Chen, Y., et al. (2022). 14 year applications of chemical fertilizers and crop straw effects on soil labile organic carbon fractions, enzyme activities and microbial community in rice-wheat rotation of middle China. Sci. Total Environ. 841:156608. doi: 10.1016/j.scitotenv.2022.156608

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, S., Xu, G., Chen, H., Zhang, M., Cao, X., Chen, M., et al. (2023). Contrasting responses of soil microbial biomass and extracellular enzyme activity along an elevation gradient on the eastern Qinghai-Tibetan plateau. Front. Microbiol. 14:974316. doi: 10.3389/fmicb.2023.974316

PubMed Abstract | Crossref Full Text | Google Scholar

Lopez-Mondejar, R., Brabcova, V., Stursova, M., Davidova, A., Jansa, J., Cajthaml, T., et al. (2018). Decomposer food web in a deciduous forest shows high share of generalist microorganisms and importance of microbial biomass recycling. ISME J. 12, 1768–1778. doi: 10.1038/s41396-018-0084-2

PubMed Abstract | Crossref Full Text | Google Scholar

Lu, D., Li, C., Sokolwski, E., Magen, H., Chen, X., Wang, H., et al. (2017). Crop yield and soil available potassium changes as affected by potassium rate in rice-wheat systems. Field Crop Res. 214, 38–44. doi: 10.1016/j.fcr.2017.08.025

Crossref Full Text | Google Scholar

Luo, C., Rodriguez-R, L. M., Johnston, E. R., Wu, L., Cheng, L., Xue, K., et al. (2014). Soil microbial community responses to a decade of warming as revealed by comparative metagenomics. Appl. Environ. Microbiol. 80, 1777–1786. doi: 10.1128/aem.03712-13

PubMed Abstract | Crossref Full Text | Google Scholar

Lv, F. L., Zheng, W., Zhai, B. N., and Li, Z. Y. (2022). Cover cropping and chemical fertilizer seasonally mediate microbial carbon and phosphorus metabolisms in an apple orchard: evidence from the enzymatic stoichiometry method. Appl. Soil Ecol. 178:104579. doi: 10.1016/j.apsoil.2022.104579

Crossref Full Text | Google Scholar

Ma, L., Li, R., Luan, H., Tang, J., Wang, L., Guo, T., et al. (2024). Impacts of long-term different fertilization regimes on microbial utilization of straw-derived carbon in greenhouse vegetable soils: Insights from its ecophysiological roles and temperature responses. Front. Plant Sci. 15:1486817. doi: 10.3389/fpls.2024.1486817

Crossref Full Text | Google Scholar

Ma, Q., Li, J., Aamer, M., and Huang, G. (2020). Increasing methane (CH4) emissions and altering rhizosphere microbial diversity in paddy soil by combining Chinese milk vetch and rice straw. PeerJ 8:e9653. doi: 10.7717/peerj.9653

PubMed Abstract | Crossref Full Text | Google Scholar

Mao, X., Yang, Y., Guan, P., Geng, L., Ma, L., Di, H., et al. (2022). Remediation of organic amendments on soil salinization: focusing on the relationship between soil salts and microbial communities. Ecotoxicol. Environ. Saf. 239:113616. doi: 10.1016/j.ecoenv.2022.113616

PubMed Abstract | Crossref Full Text | Google Scholar

Marshall, C. B., and Lynch, D. H. (2020). Soil microbial and macrofauna dynamics under different green manure termination methods. Appl. Soil Ecol. 148:103505. doi: 10.1016/j.apsoil.2020.103505

Crossref Full Text | Google Scholar

Nasar, J., Ahmad, M., Gitari, H., Tang, L., Chen, Y., and Zhou, X. B. (2024). Maize/soybean intercropping increases nutrient uptake, crop yield and modifies soil physio-chemical characteristics and enzymatic activities in the subtropical humid region based in Southwest China. BMC Plant Biol. 24:434. doi: 10.1186/s12870-024-05061-0

PubMed Abstract | Crossref Full Text | Google Scholar

Nasar, J., Zhao, C. J., Khan, R., Gul, H., Gitari, H., Shao, Z., et al. (2023). Maize-soybean intercropping at optimal N fertilization increases the N uptake, N yield and N use efficiency of maize crop by regulating the N assimilatory enzymes. Front. Plant Sci. 13:1077948. doi: 10.3389/fpls.2022.1077948

PubMed Abstract | Crossref Full Text | Google Scholar

Podosokorskaya, O. A., Bonch-Osmolovskaya, E. A., Novikov, A. A., Kolganova, T. V., and Kublanov, I. V. (2013). Ornatilinea apprima gen. nov., sp nov., a cellulolytic representative of the class Anaerolineae. Int. J. Syst. Evol. Microbiol. 63, 86–92. doi: 10.1099/ijs.0.041012-0

PubMed Abstract | Crossref Full Text | Google Scholar

Price, G. W., Langille, M. G. I., and Yurgel, S. N. (2021). Microbial co-occurrence network analysis of soils receiving short- and long-term applications of alkaline treated biosolids. Sci. Total Environ. 751:141687. doi: 10.1016/j.scitotenv.2020.141687

PubMed Abstract | Crossref Full Text | Google Scholar

Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219

PubMed Abstract | Crossref Full Text | Google Scholar

Ren, T., Feng, H., Xu, C., Xu, Q., Fu, B., Azwar, E., et al. (2022). Exogenous application and interaction of biochar with environmental factors for improving functional diversity of rhizosphere's microbial community and health. Chemosphere 294:133710. doi: 10.1016/j.chemosphere.2022.133710

PubMed Abstract | Crossref Full Text | Google Scholar

Ren, S., Guo, X., Lu, A., Guo, X., Wang, Y., Sun, G., et al. (2018). Effects of co-composting of lincomycin mycelia dregs with furfural slag on lincomycin degradation, maturity and microbial communities. Bioresour. Technol. 265, 155–162. doi: 10.1016/j.biortech.2018.05.087

PubMed Abstract | Crossref Full Text | Google Scholar

Ren, J., Liu, X., Yang, W., Yang, X., Li, W., Xia, Q., et al. (2021). Rhizosphere soil properties, microbial community, and enzyme activities: short-term responses to partial substitution of chemical fertilizer with organic manure. J. Environ. Manag. 299:113650. doi: 10.1016/j.jenvman.2021.113650

PubMed Abstract | Crossref Full Text | Google Scholar

Shu, X., He, J., Zhou, Z., Xia, L., Hu, Y., Zhang, Y., et al. (2022). Organic amendments enhance soil microbial diversity, microbial functionality and crop yields: a meta-analysis. Sci. Total Environ. 829:154627. doi: 10.1016/j.scitotenv.2022.154627

PubMed Abstract | Crossref Full Text | Google Scholar

Sinsabaugh, R. L., Hill, B. H., and Shah, J. J. F. (2009). Ecoenzymatic stoichiometry of microbial organic nutrient acquisition in soil and sediment. Nature 462, 795–U117. doi: 10.1038/nature08632

Crossref Full Text | Google Scholar

Stefan, L., Hartmann, M., Engbersen, N., Six, J., and Schoeb, C. (2021). Positive effects of crop diversity on productivity driven by changes in soil microbial composition. Front. Microbiol. 12:660749. doi: 10.3389/fmicb.2021.660749

Crossref Full Text | Google Scholar

Su, Y., Yu, M., Xi, H., Lv, J., Ma, Z., Kou, C., et al. (2020). Soil microbial community shifts with long-term of different straw return in wheat-corn rotation system. Sci. Rep. 10:6360. doi: 10.1038/s41598-020-63409-6

PubMed Abstract | Crossref Full Text | Google Scholar

Sun, L., Lin, C., Zhou, Z., Zhang, F., Xu, M., Jiao, P., et al. (2023). Characteristics of organic pollutants and their effects on the microbial composition and activity in the industrial soils of Pearl River Delta, China. Ecotoxicol. Environ. Saf. 256:114844. doi: 10.1016/j.ecoenv.2023.114844

PubMed Abstract | Crossref Full Text | Google Scholar

Tagar, A. A., Adamowski, J., Memon, M. S., Minh Cuong, D., Mashori, A. S., Soomro, A. S., et al. (2020). Soil fragmentation and aggregate stability as affected by conventional tillage implements and relations with fractal dimensions. Soil Tillage Res. 197:104494. doi: 10.1016/j.still.2019.104494

Crossref Full Text | Google Scholar

Tan, Y., Chen, H., Xiao, W., Meng, F., and He, T. (2021). Influence of farmland marginalization in mountainous and hilly areas on land use changes at the county level. Sci. Total Environ. 794:149576. doi: 10.1016/j.scitotenv.2021.149576

PubMed Abstract | Crossref Full Text | Google Scholar

Tian, S., Zhu, B., Yin, R., Wang, M., Jiang, Y., Zhang, C., et al. (2022). Organic fertilization promotes crop productivity through changes in soil aggregation. Soil Biol. Biochem. 165:108533. doi: 10.1016/j.soilbio.2021.108533

Crossref Full Text | Google Scholar

Tlaskal, V., and Baldrian, P. (2021). Deadwood-inhabiting bacteria show adaptations to changing carbon and nitrogen availability during decomposition. Front. Microbiol. 12:685303. doi: 10.3389/fmicb.2021.685303

PubMed Abstract | Crossref Full Text | Google Scholar

Tolosana-Delgado, R. (2008). Compositional data analysis in a nutshell. Göttingen, Germany: University of Gottingen on-line reference.

Google Scholar

Tsutsui, S., Hatano, T., Funada, R., and Kaneko, S. (2022). Microorganisms capable of producing polysaccharides from d-xylose. J. Appl. Glycosci. 69, 83–89. doi: 10.5458/jag.jag.JAG-2022_0008

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Q. Y., Hu, B., and Yu, H. W. (2016). Adsorption behaviors of fungicide-derived copper onto various size fractions of aggregates from orchard soil. Environ. Sci. Pollut. Res. 23, 24983–24990. doi: 10.1007/s11356-016-7743-5

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, M., Liu, Z., Zhai, B., Zhu, Y., and Xu, X. (2023). Long-term straw mulch underpins crop yield and improves soil quality more efficiently than plastic mulch in different maize and wheat systems. Field Crop Res. 300:109003. doi: 10.1016/j.fcr.2023.109003

Crossref Full Text | Google Scholar

Wang, Y., Pang, J., Zhang, M., Tian, Z., Wei, T., Jia, Z., et al. (2023). Is adding biochar be better than crop straw for improving soil aggregates stability and organic carbon contents in film mulched fields in semiarid regions? –evidence of 5-year field experiment. J. Environ. Manag. 338:117711. doi: 10.1016/j.jenvman.2023.117711

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, C., Qu, L., Yang, L., Liu, D., Morrissey, E., Miao, R., et al. (2021). Large-scale importance of microbial carbon use efficiency and necromass to soil organic carbon. Glob. Change Biol. 27, 2039–2048. doi: 10.1111/gcb.15550

PubMed Abstract | Crossref Full Text | Google Scholar

Wu, B., Wang, P., Devlin, A. T., Chen, L., Xia, Y., Zhang, H., et al. (2021). Spatial and temporal distribution of bacterioplankton molecular ecological networks in the Yuan river under different human activity intensity. Microorganisms 9:1532. doi: 10.3390/microorganisms9071532

PubMed Abstract | Crossref Full Text | Google Scholar

Wu, Y., Wang, E., Gong, W., Xu, L., Zhao, Z., He, D., et al. (2023). Soybean yield variations and the potential of intercropping to increase production in China. Field Crop Res. 291:108771. doi: 10.1016/j.fcr.2022.108771

Crossref Full Text | Google Scholar

Xia, M., Ma, X., Liu, J., Wu, M., Li, Z., and Liu, M. (2023). Potential effect of key soil bacterial taxa on the increase of rice yield under milk vetch rotation. Front. Microbiol. 14:1150505. doi: 10.3389/fmicb.2023.1150505

PubMed Abstract | Crossref Full Text | Google Scholar

Xie, Z., Shah, F., and Zhou, C. (2022a). Combining rice straw biochar with leguminous cover crop as green manure and mineral fertilizer enhances soil microbial biomass and rice yield in South China. Front. Plant Sci. 13:778738. doi: 10.3389/fpls.2022.778738

PubMed Abstract | Crossref Full Text | Google Scholar

Xie, Z., Zou, H., Zheng, Y., and Fu, S. F. (2022b). Improving anaerobic digestion of corn straw by using solid-state urea pretreatment. Chemosphere 293:133559. doi: 10.1016/j.chemosphere.2022.133559

Crossref Full Text | Google Scholar

Xin, J., Yan, L., and Cai, H. (2024). Response of soil organic carbon to straw return in farmland soil in China: a meta-analysis. J. Environ. Manag. 359:121051. doi: 10.1016/j.jenvman.2024.121051

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, Z., Sun, R., He, T., Sun, Y., Wu, M., Xue, Y., et al. (2023). Disentangling the impact of straw incorporation on soil microbial communities: enhanced network complexity and ecological stochasticity. Sci. Total Environ. 863:160918. doi: 10.1016/j.scitotenv.2022.160918

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, P., Zhu, J., Wang, H., Shi, L., Zhuang, Y., Fu, Q., et al. (2021). Regulation of soil aggregate size under different fertilizations on dissolved organic materials, cellobiose hydrolyzing microbial community and their roles in organic materials mineralization. Sci. Total Environ. 755:142595. doi: 10.1016/j.scitotenv.2020.142595

Crossref Full Text | Google Scholar

Yang, W., Li, X., and Zhang, Y. (2022). Research progress and the development trend of the utilization of crop straw biomass resources in China. Front. Chem. 10:904660. doi: 10.3389/fchem.2022.904660

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, S., Gan, M., Chen, Y., Hu, Z., Xie, Y., and Feng, Q. (2022). Fabrication of lignin-containing cellulose bio-composite based on unbleached corncob and wheat straw pulp. Int. J. Biol. Macromol. 208, 741–747. doi: 10.1016/j.ijbiomac.2022.03.192

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, H., Sui, Y., Chen, Y., Bao, T., and Jiao, X. (2022). Soil organic carbon mineralization and its temperature sensitivity under different substrate levels in the mollisols of northeast China. Life 12:712. doi: 10.3390/life12050712

Crossref Full Text | Google Scholar

Yu, K., Xiao, S., Shen, Y., Liu, S., and Zou, J. (2024). Enhanced carbon sinks following double-rice conversion to green manure-rice cropping rotation systems under optimized nitrogen fertilization in Southeast China. Agric. Ecosyst. Environ. 362:108845. doi: 10.1016/j.agee.2023.108845

Crossref Full Text | Google Scholar

Yuan, L., Gao, Y., Mei, Y., Liu, J., Kalkhajeh, Y. K., Hu, H., et al. (2023). Effects of continuous straw returning on bacterial community structure and enzyme activities in rape-rice soil aggregates. Sci. Rep. 13:2357. doi: 10.1038/s41598-023-28747-1

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, N., Bai, L., Wei, X., Li, T., Tang, Y., Wen, J., et al. (2025). Effects of organic material addition on carbon cycling and soil fertility in paddy soil. J. Environ. Manag. 379:124898. doi: 10.1016/j.jenvman.2025.124898

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, M., Cheng, G., Feng, H., Sun, B., Zhao, Y., Chen, H., et al. (2017). Effects of straw and biochar amendments on aggregate stability, soil organic carbon, and enzyme activities in the loess plateau, China. Environ. Sci. Pollut. Res. 24, 10108–10120. doi: 10.1007/s11356-017-8505-8

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, S., Fang, Y., and Xi, D. (2015). Adaptation of micro-diffusion method for the analysis of 15N natural abundance of ammonium in samples with small volume. Rapid Commun. Mass Spectrom. 29, 1297–1306. doi: 10.1002/rcm.7224

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, J., Ke, W., Liu, J., Wang, L., Chen, H., Peng, T., et al. (2019). Influence of water controlling depth on soil microflora and bacterial community diversity in paddy soil. Chin. J. Eco-Agric. 27, 277–285. doi: 10.13930/j.cnki.cjea.180373

Crossref Full Text | Google Scholar

Zhang, R., Li, Y., Zhao, X., Degen, A. A., Lian, J., Liu, X., et al. (2022). Fertilizers have a greater impact on the soil bacterial community than on the fungal community in a sandy farmland ecosystem, Inner mongolia. Ecol. Indic. 140:108972. doi: 10.1016/j.ecolind.2022.108972

Crossref Full Text | Google Scholar

Zhang, X. Q., Wang, Z. F., Can, M. Y., Bai, H., and Ta, N. (2021). Analysis of yield and current comprehensive utilization of crop straws in China. J. China Agric. Univ. 26, 30–41. doi: 10.11841/j.issn.1007-4333.2021.09.04

Crossref Full Text | Google Scholar

Zhao, J., Ni, T., Xun, W., Huang, X., Huang, Q., Ran, W., et al. (2017). Influence of straw incorporation with and without straw decomposer on soil bacterial community structure and function in a rice-wheat cropping system. Appl. Microbiol. Biotechnol. 101, 4761–4773. doi: 10.1007/s00253-017-8170-3

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, Y., Wang, H., Guo, T., Li, Z., Mi, W., and Cao, Z. (2023). Response of soil C-, N-, and P- acquisition enzymes to moisture pulses in desert grassland to shrubland state transition. Sci. Total Environ. 861:160569. doi: 10.1016/j.scitotenv.2022.160569

PubMed Abstract | Crossref Full Text | Google Scholar

Zheng, R., Cai, R., Liu, R., Liu, G., and Sun, C. (2021). Maribellus comscasis sp. Nov., a novel deep-sea Bacteroidetes bacterium, possessing a prominent capability of degrading cellulose. Environ. Microbiol. 23, 4561–4575. doi: 10.1111/1462-2920.15650

PubMed Abstract | Crossref Full Text | Google Scholar

Zhong, C., Liu, Y., Xu, X., Yang, B., Aamer, M., Zhang, P., et al. (2021). Paddy-upland rotation with Chinese milk vetch incorporation reduced the global warming potential and greenhouse gas emissions intensity of double rice cropping system. Environ. Pollut. 276:116696. doi: 10.1016/j.envpol.2021.116696

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, J., Deng, Y., Luo, F., He, Z., and Yang, Y. (2011). Phylogenetic molecular ecological network of soil microbial communities in response to elevated CO2. MBio 2, e00122–00111. doi: 10.1128/mBio.00122-11

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, Y., Yao, M., Miao, Y., Jin, J., Gu, C., Zhao, N., et al. (2022). The developmental of fresh broad bean industry in China. J. Agric. For. 12, 80–84. doi: 10.11923/j.issn.2095-4050.cjas2021-0137

Crossref Full Text | Google Scholar

Zhu, M., Yang, S., Ai, S., Ai, X., Jiang, X., Chen, J., et al. (2020). Artificial soil nutrient, aggregate stability and soil quality index of restored cut slopes along altitude gradient in Southwest China. Chemosphere 246:125687. doi: 10.1016/j.chemosphere.2019.125687

PubMed Abstract | Crossref Full Text | Google Scholar

Zuo, F. L., Li, X. Y., Yang, X. F., Wang, Y., Ma, Y. J., Huang, Y. H., et al. (2020). Soil particle-size distribution and aggregate stability of new reconstructed purple soil affected by soil erosion in overland flow. J. Soils Sediments 20, 272–283. doi: 10.1007/s11368-019-02408-1

Crossref Full Text | Google Scholar

Keywords: straw mulch, milk vetch, soil carbon cycle, molecular ecological networks, faba bean productivity

Citation: Ren K, Song W, Wei Z, Song L, Liu M, Zhou Y, Zhen Y, Wu X, Gu K, Simarani K and Wang L (2025) Organic material mulching regulated core microbial groups to promote soil carbon and nitrogen cycling and improve faba bean productivity under a triple-cropping system in purple soil hilly region of southwest China. Front. Microbiol. 16:1602633. doi: 10.3389/fmicb.2025.1602633

Received: 31 March 2025; Accepted: 07 July 2025;
Published: 16 July 2025.

Edited by:

Samina Mehnaz, Forman Christian College, Pakistan

Reviewed by:

Asma Imran, National Institute for Biotechnology and Genetic Engineering, Pakistan
Muhammad Farhan Ul Haque, University of the Punjab, Pakistan

Copyright © 2025 Ren, Song, Wei, Song, Liu, Zhou, Zhen, Wu, Gu, Simarani and Wang. 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: Longchang Wang, d2FuZ2xjMjAwM0AxNjMuY29t; Khanom Simarani, aGFub21fc3NAdW0uZWR1Lm15

Disclaimer: 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.