- 1Key Laboratory of Farmland Eco-Environment of Hebei/ College of Resources and Environmental Sciences, Hebei Agricultural University, Baoding, China
- 2State Key Laboratory of North China for Crop Improvement and Regulation/North China Key Laboratory for Crop Germplasm Resources of Education Ministry, Hebei Agricultural University, Baoding, China
- 3Biology Institute, Hebei Academy of Sciences, Shijiazhuang, China
- 4Laboratory of Crop Germplasm Resources of Hebei Province/ College of Agronomy, Hebei Agricultural University, Baoding, China
- 5School of Environmental Science and Engineering, Hebei University of Science and Technology, Shijiazhuang, China
The application of straw-degrading microbes (SDMs) with straw returned to the field is an effective measure to improve soil quality, increase yield, and maintain soil microorganisms. However, the utilization of SDMs in winter in north China is limited by the poor effects at low temperatures. This study investigated the effects of a new compound SDM, including a novel low-temperature fungus Pseudogymnoascus sp. SDF-LT, on winter wheat yield, soil improvement, and soil microbial diversity. A 2-year field experiment was conducted in two different soil textures of wheat–maize rotation fields with full corn straw return and application of SDMs at an amount of 67.5 kg hm−2. After 2 years of continuous application of SDMs, the winter wheat yield increased significantly, reaching 9419.40 kg hm−2 in Ningjin (NJSDM) and 9107.25 kg hm−2 in Mancheng (MCSDM). The soil properties have been significantly improved compared with the single straw return group, especially the sandy loam soil, whose quality is relatively low. The analysis of soil microbial diversity showed that SDMs significantly reduced the Chao1, Shannon, Simpson, and observed species of the sandy loam soil in the MCSDM group. The Simpson and Shannon indexes of fungi diversity in the two experimental sites were significantly increased by SDMs. The negative correlation of fungi increased from 47.1 to 48.85% in the SDM groups. The soil-dominant microbes changed in the SDM groups, in which the interactions between microbes were enhanced. These results suggested that the SDMs changed the the soil microbial community structure and its diversity and complexity, which is beneficial for crop growth. Our study provided sufficient evidence for the utilization of low-temperature SDMs with straw return in cold winter, which plays a role in soil improvement, especially for low-quality soils, to increase crop yield.
Introduction
China is a large agricultural country that has produced more than 819 million tons of crop residues in 2021 (Zhao et al., 2019; Bai et al., 2023). The winter wheat–summer maize rotation system is the largest planting mode in north China that produces more than 500 million tons of straw annually (Islam et al., 2023). The straw resources are excessive with the increase in crop yield resulting in potential environmental problems unless the straws are recycled appropriately (Wang et al., 2022). Straw burning and disposal are forbidden in cases of adverse effects on agriculture and environment, which is in line with the concept of sustainable agricultural development (Hong et al., 2016). Therefore, the outstanding contribution of straw return to the agricultural ecological environment has gradually emerged.
Straw returning technology is currently the most effective and greenish way to utilize straw resources, which have been improved to increase soil organic carbon, maintain soil quality and health, promote crop yield, and improve plant immunity (Tang et al., 2019; Hao et al., 2022; Yue et al., 2023). Straw is rich in organic matter, nitrogen, phosphorus, potassium, and other trace elements, which makes it an ideal fertilizer for crop growth. However, the main components of straw, including lignin, cellulose, and hemicellulose, are all aromatic polymer compounds with complex three-dimensional structures that are difficult to degrade in the natural environment, resulting in limited effects of the application of in situ straw return (Bilal et al., 2017; Gong et al., 2020). Furthermore, the high C/N ratio of the returned straws could result in N immobilization, leading to a lack of available nitrogen (AN) in the soil, which needs to be compensated by additional N application, affecting crop growth in turn. Additionally, slowly degrading straws provide a more habitable place for pathogens and eggs, thus causing aggravation of crop pests and diseases (Li et al., 2016; Liu et al., 2016).
There are abundant microorganisms that can produce cellulase and degrade straws. Microorganisms play a crucial role in natural straw degradation. However, only 200 microorganisms have been found and utilized (Taha et al., 2015). Fungi are most widely applied among bacteria, fungi, and actinomycetes that could degrade cellulose. This is because of the ability of fungi to produce strong enzymes and deep penetration of the hyphae into the plant stratum corneum (Gao et al., 2020). The compound microbes commonly function better to degrade cellulose and lignin resulting in an increased straw degradation rate (Chen et al., 2019; Chu et al., 2021). Exogenous application of compound bacterial systems plays a positive role in the abundance and diversity of the soil-dominant bacteria, thereby changing the indigenous microbial community structure of bacteria and fungi (Chen et al., 2021; Nigussie et al., 2021). Previous research found that adding exogenous lignocellulose decomposition microbes changed the indigenous soil microbial community structure and increased crop yield, and it also increased the functional abundance of soil microorganisms (Phukongchai et al., 2022; Wang et al., 2023).
Most of the commercial straw-degrading microbes (SDMs) are functional at warm temperature. However, agricultural land is located broadly in the temperate zone that harbors a cold winter, during which the long-term low temperature of the earth is fatal for SDMs. Many SDMs cease growth and metabolism in winter, resulting in the loss of their straw degradation function. In contrast, the low-temperature straw decomposition bacteria are unconstrained by cold temperatures. Thereby, low-temperature SDMs fully utilize their superior cellulose digestion abilities under natural conditions, regardless of the effects of low and warm temperatures, and save energy and production costs during cultivation (Huang et al., 2020). Thus, we have previously isolated a fungus from the wheat–maize rotation areas that degraded corn straw effectively under 10°C. In this study, we aimed to investigate the effect of compound microbial systems on soil quality, wheat yield, and soil microbial community structure in different soil textures to support the further application of the low-temperature SDMs.
Materials and methods
Experimental site description
The experimental sites were located in Ningjin County (NJ) of Xintai City and Mancheng County (MC) of Baoding City in Hebei Province of China (37°36′–38°43′N, 114°43′–115°6′E). These regions had typical temperate continental monsoon climates. The average annual temperature is 12.82–12.90°C. The mean annual precipitation is 499.0–546.5 mm, mainly occurring in June, July, and August. The tested soil was classified into two types: medium loam soil (MLS) in NJ and sandy loam soil (SLS) in MC. The general chemical properties of soil are shown in Table 1.
Experimental design
This study was initiated in 2021 and conducted in a rotation system of winter wheat and summer maize, in which winter wheat was planted in October and harvested in June of the next year, and summer maize was directly seeded and grown from June to October. The basal compound fertilizer was applied before wheat sowing, which contained 280.5 kg hm−2 N, 132 kg hm−2 P2O5, and 60 kg hm−2 K2O. Topdressing was applied at the regreening stage of wheat with 375 kg hm−2 urea. The corn straws were chopped into pieces of 3–5 cm (length) before mulching on the surface of the soil, followed by fully retained by rotary tilling (depth of 0–10 cm). The planting density of wheat (variety No. 1 of Malan) was 262.50 kg hm−2 with a row spacing of 20 cm. The crop cultivation and farming management were in accordance with the local conventional cultivation practices. The straw return treatment without a straw-degrading agent served as a control check (CK), and the straw return with the application of the SDMs with an amount of 67.5 kg hm−2 was set as the experimental groups. The experiments were performed in NJ and MC in the same year, with a randomized block design with three replicates for three plots, and each plot area was 100 m2 (10.0 m × 10.0 m).
The SDMs were composed of Pseudogymnoascus sp. SDF-LT and Bacillus licheniformis. The effective number of viable microbes was 100 million g−1. The author previously isolated the fungus Pseudogymnoascus sp. SDF-LT with a strong cellulose degradation ability at 10°C, whose registration number is OM 304630 on NCBI.
Soil sampling and analysis
A previous study on the rate of straw degradation by the SDMs was performed in NJ with the string-bag method, which showed a higher straw degradation rate until the wheat filling stage (Supplementary Figure S1). Since the filling stage is important for wheat yield, the soil samples were collected immediately after the wheat filling stage. For each plot, five cores (0–20 cm in depth and 5 cm in diameter) were randomly selected to collect the soil samples. All soil samples were collected into separate sterile plastic bags and stored temporarily in an icebox before being examined in the lab. After removing leaf litter, rocks, and roots, the soil samples were divided into two subsamples: one was stored at −80°C for the determination of microbial diversity, and the other was air-dried for the determination of chemical properties.
Soil pH was determined with a water: soil ratio of 2.5:1 (V/W) using a pH meter (Luo et al., 2017). Soil organic matter (SOM) was determined using the volumetric method. Soil total nitrogen (TN) was measured by the micro-Kjeldahl method. Soil AN was determined by the alkali-hydrolysis diffusion method. Soil available phosphorous (AP) was extracted with 0.5 M NaHCO3 solution and analyzed by the Mo-Sb colorimetric method. Soil available potassium (AK) was determined by flame photometry (FP640, Shanghai, China).
During the winter wheat harvest, spike number, grain number, 1000-grain weight, and yield were measured.
High-throughput sequencing
The microbial DNA was extracted from the 12 soil samples using the OMEGA Soil DNA Kit (M5635-02, Omega Bio-Tek, Norcross, GA, United States), followed by the quantity and quality examination using a NanoDrop NC2000 spectrophotometer (ThermoFisher Scientific, Waltham, MA, United States) and agarose gel electrophoresis, respectively. The V3–V4 regions of the 16S rRNA genes of the soil bacteria were amplified using primers 338F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′-GGAC TACHVGGGTWTCTAAT-3′). Fungal ITS rRNA genes were amplified with primers ITS1F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and ITS2R (5′-GCTGCGTTCTTCATCGATGC-3′). Sample-specific 7-bp barcodes were incorporated into the primers for multiplex sequencing. The PCR amplicons were purified with Vazyme VAHTS DNA Clean Beads (Vazyme, Nanjing, China) and quantified using the Quant-iT PicoGreen dsDNA Assay Kit (Invitrogen, Carlsbad, CA, United States). The amplicons were then pooled in equal amounts and pair-end (2,250 bp) sequenced with NovaSeq 6,000 SP Reagent Kit (500 cycles) on the Illumina NovaSeq platform at Shanghai Personal Biotechnology Co., Ltd. (Shanghai, China).
Bioinformatic analysis
The raw sequence data were demultiplexed using the demux plugin, followed by cutting primers with the cutadapt plugin (Martin, 2011). Then the DADA2 plugin was applied for filtering, dereplication, chimera identification, and merging paired-end reads (Callahan et al., 2016) for further analysis. The alpha diversity indexes, including the Chao1 richness estimator, Shannon diversity, Simpson, and observed species, were calculated using the Amplicon Sequence Variant (ASV) table in QIIME2 and visualized as box plots. Beta diversity analysis was performed to investigate the structural variation of microbial communities across samples using Jaccard metrics (Jaccard, 1908) and visualized using principal coordinate analysis (PCoA). A Venn diagram was generated to visualize the shared and unique ASV among groups based on their occurrence, regardless of their relative abundance (Zaura et al., 2009). Taxa abundances at the ASV levels were statistically compared among samples or groups by MetagenomeSeq and visualized as Manhattan plots (Zgadzaj et al., 2016). Linear discriminant analysis effect size (LEfSe) was performed to detect differentially abundant taxa across groups using the default parameters (Segata et al., 2011). Co-occurrence network analysis was performed by SparCC analysis. The co-occurrence network was analyzed with the igraph package and visualized by Gephi 0.9.2 software.
Statistical analysis
All data were displayed as mean ± standard deviation (SD). Statistical analysis was conducted with IBM SPSS Statistics 22 software (SPSS IBM Corp). Data variance was evaluated by one-way analysis. Means were compared by Duncan’s new multiple range test at the 5% level. Relationships between soil microbial ecological diversity and soil physical and chemical properties were analyzed by Pearson analysis. A p-value of <0.05 was considered statistically significant.
Results
Effects of straw-degrading microbes on soil’s physical and chemical properties
The application of SDMs in different soil textures significantly promoted soil chemical properties (Figures 1A–E). Two years of continuous application of the SDMs promoted the content of organic matter, total N, alkaline hydrolysis N, available P, and available K by 26.39, 15.32, 9.66, 3.94, and 15.19%, respectively, in SLS. In MLS, the content of organic matter, total N, alkaline hydrolysis N, available P, and available K was increased by 6.82, 3.62, 7.88, 28.39, and 32.26%, respectively. Soil C and N are the crop nutrition sources that were increased much higher in SLS than in MLS, and the P and K that promote crop development increased much higher in MLS compared with SLS. These results revealed that the SDMs significantly improved the soil nutrients, especially the SLS. The soil bulk density of MLS and SLS with the application of SDMs was higher than that of CK in the second year, but the difference was not significant (Figure 1F). These results suggested that the continuous application of SDMs with straw return would improve soil quality, especially for barren soil.
Figure 1. Effects of SDMs on soil physical and chemical properties, including organic matter (A), total N (B), alkaline hydrolysis N (C), available P (D), available K (E), and soil bulk density (F). The different lowercase letters stand for significant differences at a p-value of <0.05. CK is the control that the straw returned without SDMs in Ningjin (NJCK) and Mancheng (MCCK). The SDM is the treatment that the straw returned with the application of SDM in Ningjin (NJSDM) and Mancheng (MCSDM).
Effects of straw-degrading microbes on wheat growth
Changes in soil chemical and biological properties can significantly affect crop yield. Thus, the winter wheat yield was evaluated in 2023. The yield traits of wheat were higher in the SDM groups in both of the regions, including the numbers of spikes, grains per spike, and weight of 1,000 grains (Supplementary Figure S2). Although the yield traits that measured were not significantly varied, the total yield of winter wheat was significantly higher with the application of SDMs in both farmlands (Figure 2). The yield of winter wheat was increased by 17.4 and 9.5%, respectively, in MC and NJ after the application of SDMs. Thus, the compound SDM was effective in promoting the yield of winter wheat which were grown on soils with different textures.
Figure 2. Effects of SDMs on the yield of winter wheat in the second year. Asterisks refer to significant differences at a p-value of < 0.05 level.
Effects of straw-degrading microbes on the α-diversity of soil microbial community
Plants, animals, and microbes are the main components of the agricultural ecosystem. In order to investigate the effects of microbes on wheat after the application of SDMs, we analyzed the soil microbial diversity through sequencing. We assessed a total of 461,206 and 374,949 sequences (average 115,302 and 93,737 per sample) for bacteria and fungi, respectively. The results of a-diversity analysis demonstrated that the bacterial diversity and richness of the soils in different places were significantly different from each other (Figure 3A). The application of SDMs in the SLSs significantly decreased the bacterial Chao1, Shannon, and Simpson indexes and observed species. In contrast, the opposite results were obtained from the soil fungi (Figure 3B). The Chao1, Shannon, and Simpson indexes and observed species of fungi were only a little higher in the SLSs applied with SDMs. The Shannon and Simpson indexes of fungi were significantly more abundant with the application of SDMs in NJ, but the Chao1 and observed species indexes were lower. These results suggested that soil texture could affect soil microbial diversity, and the effects of the application of SDMs on soil microbial diversity might be determined by soil texture.
Figure 3. α-diversity of soil bacteria and fungi among the treatments in two experimental sites in the second year. (A) α-diversity of soil bacteria. (B) α-diversity of soil fungi. The different lowercase letters stand for significant differences at a p-value of <0.05 level.
Effect of straw-degrading microbes on the β-diversity of soil microbial community
The PCoA analysis showed that the three repeated microbial communities in different treatments all converged in one area and separated from the other treatments, suggesting the consistency of the replications. The samples from NJ were distributed far from the MC samples, suggesting great differences between the two soils (Supplementary Figure S3). The total variance explained by bacterial (PCo1) and fungal (PCo2) communities was 42.4 and 65.4%, respectively, between the four experimental groups, which suggested a greater effect of the SDMs on fungal diversity. The bacterial and fungal communities showed greater differences in MC, which was consistent with the α-diversity analysis.
A Venn diagram analysis of OTUs by applying straw-degrading microbes
A Venn diagram was used to compare the OTUs in the four groups. A total of 2.72% (1348) of the bacterial OTUs and 5.25% (94) of the fungal OTUs were shared in the four experimental groups (Figure 4). The specific bacteria in each group took 77.38% (11611), 76.38% (11252), 79.42% (12645), and 76.66% (10766) of the NJCK, NJSDM, MCCK, and MCSDM OTUs, respectively (Supplementary Figures S3A,B). While the specific fungi in each group took 61.55% (325), 78.15% (273), 64.16% (471), and 64.51% (478) of the NJCK, NJSDM, MCCK, and MCSDM OTUs, respectively (Supplementary Figures S4C,D). These results showed that the unique bacterial OTU composition was reduced, while the unique fungal OTUs were increased by applying SDMs.
Microbial community structure and composition
Furthermore, the microbial community structure was investigated at the phylum level. The predominant bacterial phyla were Proteobacteria, Actinobacteriota, Acidobacteriota, Chloroflexi, and Gemmatimonadota, with a relative abundance above 5% (Figure 5A). Proteobacteria was the most dominant bacterial phylum in NJ. However, the most dominant bacterial phylum in MC was Actinobacteriota. The application of the SDMs increased the relative abundance of Actinobacteriota in SLS. The dominant fungal phylum was Ascomycota in all the experimental groups (Figure 5B). The SDMs increased the relative abundance of Mortierellomycota in NJ, while increased the relative abundance of Basidiomycota in MC.
Figure 5. Microbial composition at the phylum level. (A) The 16S rRNA gene-based bacterial community composition. (B) The ITS gene-based fungal community composition.
Subsequently, a genus-level comparison was made between the microbial communities of the experimental groups. The genera type with relative abundance above 1% was the same, while the ratio of each genus varied between the two experimental places. The bacterial genus Vicinamibacteraceae was commonly abundant in all the groups with the highest relative abundance above 5% (Figure 6A). The application of SDMs could increase the relative abundance of Vicinamibacteraceae in MLS, but decrease in SLS. Furthermore, the application of SDMs could also decrease the relative abundance of Sphingomonas. At the genus level of fungus, the preponderant fungal genus was quite different in the two experimental places (Figure 6B). In NJ, Botryotrichum was the most dominant and SDMs increased the relative abundance of Schizothecium to lower the ratio of the dominant fungal genus. In MC, the SDMs greatly increased the abundance of the genus Subulicystidium (from 6.55 to 12.38%), Cladorrhinum (from 3.20 to 6.77%), and Apodus (from 1.23 to 7.72%). The above results suggested different variations in the bacterial and fungal community structures in different soil textures after applying SDMs.
Figure 6. Microbial composition at the genus level. (A) The 16S rRNA gene-based bacterial community composition. (B) The ITS gene-based fungal community composition. The top 20 genera of bacteria and fungi are shown in different colors.
Effect of straw-degrading microbes on the compositions of bacterial and fungal biomarkers
LEfSe was applied to identify the bacterial and fungal community biomarkers (p < 0.05). Apart from the differences at phylum and genus levels, there were also other groups of microorganisms at the class, order, family, and genus levels (bacteria and fungi) that were enriched in the four experimental groups (Figure 7). Gemmatimonadetes (from class to order) were significantly enriched in NJCK. Two groups of Alphaproteobacteria and Gammaproteobacteria were significantly enriched in NJSDM. Three groups of bacteria, namely, Myxococcota, Vicinamibacteria (from class to order), and Polyangia, were significantly enriched in MCCK. Actinobacteriota (from phylum to class) were significantly enriched in MCSDM (Figure 7A). Four groups of fungi were significantly enriched in NJCK, namely, Chaetomiaceae, Botryotrichum, Sordariales, and Ascomycota. Sordariomycetes, Hypocreales, Nectriaceae, and Schizothecium were significantly enriched in NJSDM. Four fungal groups were enriched in MCCK, namely, Cucurbitariaceae, Pyrenochaetopsis, Penicillium, and Agaricales. Fungal groups, such as Agaricomycetes, Basidiomycota, and Lasiosphaeriaceae, were significantly enriched in MCSDM (Figure 7B).
Figure 7. Cladograms plotted from the LEfSe analysis show significant differences (p < 0.05) in the relative abundance of 16S rRNA gene-based bacterial taxa (A) and ITS gene-based fungal taxa (B).
Effect of straw-degrading microbes on the co-occurrence network
The symbiotic network relationship between the CK and SDM groups was obtained based on Spearman’s correlation coefficient. We found that bacteria dominated the symbiotic relationship among soil microbes and were the main components of the symbiotic systems (Figure 8). The dominant bacteria and fungi of CK groups were Vicinamibacteraceae and Botryotrichum, while the dominant bacteria of SDM groups was Sphingomonas. The correlation of bacteria decreased from 52.9 to 51.15%, while the correlation of fungi increased from 47.1 to 48.85%. The number of edges that indicate the complexity of microbial interaction increased by 40.65% in the SDMs compared with CK. The disease bacterium, Fusarium, displayed a more active correlation with the other microbes, indicating a limitation of its pathogenicity. These results suggested that the application of SDMs with straw changed the soil-dominant microbes and rebalanced their interactions, thus increasing their complexity, which is beneficial for crop growth.
Figure 8. Co-occurrence network analysis of soil bacterial and fungal communities in the CK groups (A) and the SDM groups (B) of the two experimental sites. The size of each node is proportional to the relative abundance of the genus. Each connection shown has a Spearman’s coefficient of >0.6 and a p-value of <0.05. The red color represents bacteria, and the red line shows a positive correlation. The green color represents fungi, and the green line shows a negative correlation.
Correlation analysis between soil microbial diversity and soil physicochemical properties
Correlation analysis was performed between soil microbial diversity and environmental factors. From the soil bacterial correlation analysis (Table 2), we found that OM, TN, AN, and AK were negatively correlated with the Chao1, Shannon, and Simpon indexes and observed species. Significantly negative correlations were found between pH and AK with Shannon and Simpon indexes (p < 0.05). Furthermore, TN and AN were extremely significantly negatively correlated with Shannon and Simpon indexes (p < 0.01). In contrast, AP was positively correlated with Shannon and Simpon indexes. From the soil fungal correlation analysis (Table 3), OM was positively correlated with Chao 1, Shannon, and Simpon indexes and observed species. AP was extremely significantly positively correlated with Chao 1, Shannon, and Simpon indexes and observed species (p < 0.01). These results revealed that the changes in soil nutrients after the application of SDMs significantly affected the soil microbial richness and diversity.
Table 2. Correlation analysis between soil bacterial ecological diversity and soil chemical properties.
Table 3. Correlation analysis between soil fungal ecological diversity and soil chemical properties.
Discussion
Soil is a very diverse system consisting of various organic and inorganic components. One of the most essential constituents in the soil is microorganisms, which can degrade several compounds, resulting in a highly dynamic soil environment. Straw decomposition is a process of nutrient release, organic C mineralization, and soil organic C balance mediated by SDMs (Zhao et al., 2019; Mei et al., 2020), which promotes soil fertility, changes the number of soil microorganisms, and affects plant growth (Han and He, 2010; Said-Pullicino et al., 2014). Our research confirmed that applying SDMs significantly improved wheat yield, especially the wheat yield increased by 17.4% in the SLS of MC (Figure 2).
The SDMs promoted straw decomposition (Supplementary Figure S1), which improved soil chemical properties (Figure 1). The straw’s nutrients are immobilized, but the decomposition of the straw releases C and N into the soil. Organic matter and N determine soil fertility. The significant increase in the amount of SOM, total N, and available N revealed that the application of SDMs increased soil nutrient condition, which could significantly affect both crop yield and soil respiration (Chen et al., 2017; Liang et al., 2019, 2021). Phosphorus (P) is vital for plant growth and production, as it is required as an important part of the compounds of ATP (the energy carrier). P is necessary for starch and fat metabolism, cellular respiration, N uptake, and carbohydrate transportation (Adnan et al., 2020). Thus, the optimum level of P is important to obtain maximum yield (Cui et al., 2022). Studies proved that straw returning can increase the soil phosphorus supply capacity by reducing soil phosphorus adsorption and increasing soil phosphorus desorption capacity (Ding et al., 2024). Potassium (K) is another mineral element that is essential for plant growth and development and increases wheat yield (Zhan et al., 2016). Long-term straw returning is beneficial for soil K balance and K supply ability (Zhang et al., 2021). The straws are good sources of organic matter, plants N, P, and K. The straw degradation rate increased fast with the SDMs during the wheat jointing stage and the filling stage (Supplementary Figure S1), which not only provided nutrients for wheat but also improved yield traits at the main developmental stages, resulting in a significant increase in wheat yield.
Straw decomposition affects soil microorganisms whose diversity and community structure influence key soil chemical characteristics (Zhang et al., 2023). Different levels of straw degradation with different nutrient release capacities and metabolites create different soil biological environments (Kandasamy et al., 2019). Our results were consistent with this that SDM groups increased the fungi composition and diversity (Figures 3, 4). The coexistence of bacteria and fungi in soil forms a complex system of species interactions. During the decomposition process of straw, the dominant microbial communities in the early and late stages of decomposition are different (Bastian et al., 2009). Previous research showed that bacteria are the dominant microbes in the early stage of straw decomposition, while fungi will become dominant in the soil after active straw degradation (Bastian et al., 2009). In this study, the appropriate amount of SDMs was applied considering environmental safety so that the exogenous microbes were unable to ‘overturn’ the microbial community (Ho and Ko, 1985), which is consistent with the fact that Pseudogymnoascus was not dominant in the soil (Figure 6). The degradation of the straw by Pseudogymnoascus gradually changed the soil microbial community, in which fungal richness increased (Figure 3; Supplementary Figure S4). The bacterial diversity in SLS in MC was significantly affected by SDMs, but there was little difference in MLS in NJ. The fungal diversity only significantly varied in MLS in NJ. The same effects of different composition agents were found on SLS and MLS under straw return (Sa et al., 2020). Thus, it was speculated that the effects of SDMs on soil microbial diversity were related to the composition of straw decomposing agents and soil type.
The dominant bacterial phyla in the soil were Proteobacteria, Actinobacteriota, and Acidobacteria (Figure 5A), which are commonly found in soils worldwide (Maestre et al., 2015). The dominant fungal groups were Ascomycota and Basidiomycota (Figure 5B), which prevail in terrestrial ecosystems (Tedersoo et al., 2014). The relative abundance of the bacterial genus Sphingomonas in the Proteobacteria phylum increased with the application of SDMs (Figure 6). Sphingomonas can metabolize a wide variety of carbon sources and clean up toxic substances in soil. Some possess the characteristics of nitrogen fixation and denitrification (Yang et al., 2016). The members of Actinobacteriota produce extracellular hydrolytic enzymes to decompose organic matter and enhance soil carbon cycling. The relative abundance of dominant fungal phyla in MCSDM increased significantly. The ability of cellulose and lignin degradation is particularly common in Ascomycota and Basidiomycota (Hannula et al., 2010; Zhang et al., 2019; Wu et al., 2021; Zhao et al., 2021). Basidiomycota is positively correlated with soil C and N (Zhu et al., 2022). The changes in bacterial and fungal community structure supported the higher increase of soil C and N resulting in a more effective improvement of wheat yield in MC.
Combining the LEfSe analysis with microbial community composition, we found that Alphaproteobacteria and Gammaproteobacteria were significantly enriched in NJSDM, and Actinobacteriota was significantly enriched in MCSDM (Figure 7A), which may be caused by the soil particle size and nutrient contents. The growth of microorganisms is regulated by soil fraction size and nutrients (Sessitsch et al., 2001; Lu et al., 2020; Xia et al., 2020). The effects of soil type on soil microbial community structure are determined by different soil environmental factors. Our results displayed different variations of microbial community structure related to soil nutrients (Tables 2, 3). The soil pH is considered the main factor driving the changes in soil microbial diversity (Pennanen, 2001; Oehl et al., 2010). The vital activities of the microbial community are supported by soil nutrients. Total N and available N are negatively correlated with microbial diversity, while available P is positively correlated with soil microbial diversity (Tables 2, 3), which was supported by previous research (Yue et al., 2023). Additionally, the plant root exudates are decomposed and absorbed by soil microorganisms, which accelerate soil nutrient accumulation. The microbial community is essential for soil quality.
The microorganism is a part of the soil ecosystem, whose interaction is crucial for determining the function and stability of the entire ecological network. The key species with specific functions in the soil were changed by SDM (Figure 8B). The disease bacterium, Fusarium, displayed a more active correlation with the other microorganisms, indicating a limitation of its pathogenicity. The interactions with Sphingomonas were more complex. The clustering coefficient, graph density, and average degree of SDMs were all higher than those of CK, suggesting a more complex microbial network. Furthermore, the application of SDMs with straw return increased specific fungal OTUs and decreased specific bacterial OTUs (Figure 4; Supplementary Figure S4). This shift revealed a higher ratio of fungi/bacteria, suggested as an indicator for a more sustainable agroecosystem where organic matter decomposition and N mineralization dominate the nutrient supply for crop growth (De Vries et al., 2006), which is beneficial for crop reproduction. Additionally, a large proportion of Basidiomycota can form mycorrhizas with plant roots in the soil (Kumar and Atri, 2018). The increase in this phylum indicated interactions with wheat roots to promote wheat growth and development.
Conclusion
The application of SDMs could significantly promote the soil’s chemical properties, increase wheat yield, and promote the microorganism environment in the soil. The density and complexity of the symbiotic network of SDMs were increased, and the community structure was more stable. Furthermore, the SDMs were found to be more effective in SLS, revealing a potential utilization for straw return in winter.
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 at: https://www.ncbi.nlm.nih.gov/, OM 304630.
Author contributions
YuH: Conceptualization, Funding acquisition, Investigation, Methodology, Writing – original draft, Writing – review & editing. YY: Data curation, Formal analysis, Investigation, Writing – original draft, Writing – review & editing. YM: Formal analysis, Writing – review & editing. XZ: Formal analysis, Writing – review & editing. QZ: Formal analysis, Writing – review & editing. MM: Conceptualization, Supervision, Writing – review & editing. YaH: Conceptualization, Supervision, Writing – review & editing. ZP: Conceptualization, Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was supported by the Key R&D Projects of Hebei Sciences (24305), National Key R&D Program Projects (2023YFD2301500), and the Key R&D Projects in Hebei Province (21327310D).
Acknowledgments
We are grateful to Ru Zhao for their technical help and data analysis.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1391632/full#supplementary-material
SUPPLEMENTARY FIGURE S1 | The degradation rate of corn straw at different developmental stages of wheat. CK is the control that the straw returned without straw-degrading microbes. SDM stands for experimental group that the straw returned with application of the straw-degrading microbes.
SUPPLEMENTARY FIGURE S2 | Effects of SDM on traits of wheat yield including spike number (A), number of grains per spike (B), and 1000-grain weight (C).
SUPPLEMENTARY FIGURE S3 | Effect of SDM on β-diversity of soil microbial community. Principal coordinate analysis (PCoA) of the bacterial 16S rRNA gene (A) and fungal ITS gene (B). Each colored dot represents the sample.
SUPPLEMENTARY FIGURE S4 | Venn diagrams of OTU distribution of the 16S rRNA gene (A,B) and the ITS gene (C,D) among CK and SDM treatments in two experimental sites.
References
Adnan, M., Hussain, M., Anjum, M. Z., Rehman, F. U., and Ahmad, R. (2020). Role of phosphorous in wheat production: a review. Int. J. Agric. Sci. Res. 27, 29–43. doi: 10.1051/agro:2006011
Bai, W. L., Zhang, L., Yan, L., Wang, X., and Zhou, Z. (2023). Crop straw resource utilization as pilot policy in China: An event history analysis. Int. J. Environ. Res. Public Health 20:3939. doi: 10.3390/ijerph20053939
Bastian, F., Bouziri, L., Nicolardot, B., and Ranjard, L. (2009). Impact of wheat straw decomposition on successional patterns of soil microbial community structure. Soil Biol. Biochem. 41, 262–275. doi: 10.1016/j.soilbio.2008.10.024
Bilal, M., Asgher, M., Iqbal, H. N., Hu, H. B., and Zhang, X. H. (2017). Bio-based degradation of emerging endocrine-disrupting and dye-based pollutants using cross-linked enzyme aggregates. Environ. Sci. Pollut. Res. Int. 24, 7035–7041. doi: 10.1007/s11356-017-8369-y
Callahan, B. J., Mcmurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J., and Holmes, S. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Chen, Y. H., Li, S., Liu, N., He, H., Cao, X. Y., Lv, C., et al. (2021). Effects of different types of microbial inoculants on available nitrogen and phosphorus, soil microbial community, and wheat growth in high-P soil. Environ. Sci. Pollut. Res. Int. 28, 23036–23047. doi: 10.1007/s11356-020-12203-y
Chen, K. J., Tang, J. C., Xu, B. H., Lan, S. L., and Cao, Y. (2019). Degradation enhancement of rice straw by co-culture of Phanerochaete chrysosporium and Trichoderma viride. Sci. Rep. 9:19708. doi: 10.1038/s41598-019-56123-5
Chen, Z. M., Xu, Y. H., Fan, J. L., Yu, H. Y., and Ding, W. X. (2017). Soil autotrophic and heterotrophic respiration in response to different N fertilization and environmental conditions from a cropland in Northeast China. Soil Biol. Biochem. 110, 103–115. doi: 10.1016/j.soilbio.2017.03.011
Chu, X. D., Awasthi, M., Liu, Y. Y. N., Cheng, Q., Qu, J. B., and Sun, Y. (2021). Studies on the degradation of corn straw by combined bacterial cultures. Bioresour. Technol. 320:124174. doi: 10.1016/j.biortech.2020.124174
Cui, H. X., Luo, Y. L., Chen, J., Jin, M., Li, Y., and Wang, Z. L. (2022). Straw return strategies to improve soil properties and crop productivity in a winter wheat-summer maize cropping system. Eur. J. Agron. 133:126436. doi: 10.1016/j.eja.2021.126436
De Vries, F. T., Hoffland, E., Van Eekeren, N., Brussaard, L., and Bloem, J. (2006). Fungal/bacterial ratio in grasslands with contrasting nitrogen management. Soil Biol. Biochem. 38, 2092–2103. doi: 10.1016/j.soilbio.2006.01.008
Ding, J. K., Xu, B., Ye, X. X., Sun, R. B., Xiong, Q. Z., Wen, Y. C., et al. (2024). Effects of long-term different straw returning measures on soil phosphorus supply capacity. J. Anhui Agric. Univ. 51, 110–116. doi: 10.13610/j.cnki.1672-352x.20240313.010
Gao, X. A., Wang, X., Xie, J., Wang, F. H., Gong, Y. X., Guan, F. C., et al. (2020). Research Progress on low temperature straw-degrading compound microbial agent. Biol. Bull. 36, 144–150. doi: 10.13560/j.cnki.biotech.bull.1985.2019-0888
Gong, X. J., Zou, H. T., Qian, C. R., Yu, Y., Hao, Y. B., Li, L., et al. (2020). Construction of in situ degradation bacteria of corn straw and analysis of its degradation efficiency. Ann. Microbiol. 70:62. doi: 10.1186/s13213-020-01601-9
Han, W., and He, M. (2010). The application of exogenous cellulase to improve soil fertility and plant growth due to acceleration of straw decomposition. Bioresour. Technol. 101, 3724–3731. doi: 10.1016/j.biortech.2009.12.104
Hannula, S. E., De Boer, W., and Van Veen, J. A. (2010). In situ dynamics of soil fungal communities under different genotypes of potato, including a genetically modified cultivar. Soil Biol. Biochem. 42, 2211–2223. doi: 10.1016/j.soilbio.2010.08.020
Hao, X. X., Han, X. Z., Wang, S. Y., and Li, L. J. (2022). Dynamics and composition of soil organic carbon in response to 15 years of straw return in a Mollisol. Soil Tillage Res. 215:105221. doi: 10.1016/j.still.2021.105221
Ho, W. C., and Ko, W. H. (1985). Soil microbiostasis: effects of environmental and edaphic factors. Soil Biol. Biochem. 17, 167–170. doi: 10.1016/0038-0717(85)90110-5
Hong, J. L., Ren, L. J., Hong, J. M., and Xu, C. Q. (2016). Environmental impact assessment of corn straw utilization in China. J. Clean. Prod. 112, 1700–1708. doi: 10.1016/j.jclepro.2015.02.081
Huang, Y. L., Huang, Y. Y., Ma, H. Y., Xu, B. X., Jia, Z. H., Song, S. S., et al. (2020). Fungi with corn straw decomposing characterization at low-temperature: screening and application in straw returning field. Chin Agric Sci Bull. 36, 53–60. doi: 10.11924/j.issn.1000-6850.casb20200200102
Islam, M. U. I., Jiang, F. H., Halder, M., Liu, S., and Peng, X. H. (2023). Impact of straw return combined with different fertilizations on soil organic carbon stock in upland wheat and maize croplands in China: a meta-analysis. Crop Environ. 2, 233–241. doi: 10.1016/j.crope.2023.10.003
Jaccard, P. (1908). Nouvelles recherches sur la distribution florale. Bull. Soc Vaud. Sci. Nat. 44, 223–270. doi: 10.5169/seals-268384
Kandasamy, S., Liu, E., Patterson, G., Ali, S., and Lazarovits, G. (2019). Introducing key microbes from high productive soil transforms native soil microbial community of low productive soil. Microbiol. Open. 8:e895. doi: 10.1002/mbo3.895
Kumar, J., and Atri, N. S. (2018). Studies on Ectomycorrhiza: An appraisal. Bot. Rev. 84, 108–155. doi: 10.1007/s12229-017-9196-z
Li, T., He, C. E., Ge, X. Y., and Ouyang, Z. (2016). Responses of soil mineral N contents, enzyme activities and crop yield to different C/N ratio mediated by straw retention and N fertilization. Chin. J. Eco Agric. 24, 1633–1642. doi: 10.13930/j.cnki.cjea.160357
Liang, G. P., Cai, A. D., Wu, H. J., Wu, X. P., Houssou, A. A., Ren, C. J., et al. (2019). Soil biochemical parameters in the rhizosphere contribute more to changes in soil respiration and its components than those in the bulk soil under nitrogen application in croplands. Plant Soil 435, 111–125. doi: 10.1007/s11104-018-3886-0
Liang, G. P., Wu, X. P., Cai, A. D., Dai, H. C., Zhou, L. P., Cai, D. X., et al. (2021). Correlations among soil biochemical parameters, crop yield, and soil respiration vary with growth stage and soil depth under fertilization. Agron. J. 113, 2450–2462. doi: 10.1002/agj2.20699
Liu, T., Chen, X. Y., Hu, F., Ran, W., Shen, Q. R., Li, H. X., et al. (2016). Carbon-rich organic fertilizers to increase soil biodiversity: evidence from a meta-analysis of nematode communities. Agric. Ecosyst. Environ. 232, 199–207. doi: 10.1016/j.agee.2016.07.015
Lu, P. N., Bainard, L. D., Ma, B., and Liu, J. H. (2020). Bio-fertilizer and rotten straw amendments alter the rhizosphere bacterial community and increase oat productivity in a saline-alkaline environment. Sci. Rep. 10:19896. doi: 10.1038/s41598-020-76978-3
Luo, X. S., Han, S., Lai, S. S., Huang, Q. Y., and Chen, W. L. (2017). Long-term straw returning affects Nitrospira-like nitrite oxidizing bacterial community in a rapeseed-rice rotation soil. J. Basic Microbiol. 57, 309–315. doi: 10.1002/jobm.201600400
Maestre, F. T., Delgado-Baquerizo, M., Jeffries, T. C., Eldridge, D. J., Ochoa, V., Gozalo, B., et al. (2015). Increasing aridity reduces soil microbial diversity and abundance in global drylands. Proc. Natl. Acad. Sci. U S A. 112, 15684–15689. doi: 10.1073/pnas.1516684112
Martin, M. (2011). CUTADAPT removes adapter sequences from high-throughput sequencing reads. EMBnet J. 17, 10–12. doi: 10.14806/ej.17.1.200
Mei, J. F., Shen, X. B., Gang, L., Xu, H. J., Wu, F. F., and Sheng, L. Q. (2020). A novel lignin degradation bacteria-Bacillus amyloliquefaciens SL-7 used to degrade straw lignin efficiently. Bioresour. Technol. 310:123445. doi: 10.1016/j.biortech.2020.123445
Nigussie, A., Dume, B., Ahmed, M., Mamuye, M., Ambaw, G., Berhiun, G., et al. (2021). Effect of microbial inoculation on nutrient turnover and lignocellulose degradation during composting: a meta-analysis. Waste Manag. 125, 220–234. doi: 10.1016/j.wasman.2021.02.043
Oehl, F., Laczko, E., Bogenrieder, A., Stahr, K., Bösch, R., Van Der Heijden, M., et al. (2010). Soil type and land use intensity determine the composition of arbuscular mycorrhizal fungal communities. Soil Biol. Biochem. 42, 724–738. doi: 10.1016/j.soilbio.2010.01.006
Pennanen, T. (2001). Microbial communities in boreal coniferous forest humus exposed to heavy metals and changes in soil pH—a summary of the use of phospholipid fatty acids, biolog® and 3H- thymidine incorporation methods in field studies. Geoderma 100, 91–126. doi: 10.1016/s0016-7061(00)00082-3
Phukongchai, W., Kaewpradit, W., and Rasche, F. (2022). Inoculation of cellulolytic and ligninolytic microorganisms accelerates decomposition of high C/N and cellulose rich sugarcane straw in 580 tropical sandy soils. Appl. Soil Ecol. 172:104355. doi: 10.1016/j.apsoil.2021.104355
Sa, R. L., Yang, H. S., Tai, J. C., Gao, J. L., and Li, Y. (2020). Effect of straw maturing agents on fungal diversity in soil with different textures under returned straw conditions. Chin. J. Eco Agric. 28, 1061–1071. doi: 10.13930/j.cnki.cjea.190900
Said-Pullicino, D., Cucu, M. A., Sodano, M., Birk, J. J., Glaser, B., and Celi, L. (2014). Nitrogen immobilization in paddy soils as affected by redox conditions and rice straw incorporation. Geoderma 228-229, 44–53. doi: 10.1016/j.geoderma.2013.06.020
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60
Sessitsch, A., Weilharter, A., Gerzabek, M. H., Kirchmann, H., and Kandeler, E. (2001). Microbial population structures in soil particle size fractions of a long-term fertilizer field experiment. Appl. Environ. Microbiol. 67, 4215–4224. doi: 10.1128/AEM.67.9.4215-4224.2001
Taha, M., Shahsavari, E., Al-Hothaly, K., Mouradov, A., Smith, A. T., Ball, A. S., et al. (2015). Enhanced biological straw saccharification through coculturing of lignocellulose-degrading microorganisms. Appl. Biochem. Biotechnol. 175, 3709–3728. doi: 10.1007/s12010-015-1539-9
Tang, L. L., Nie, S. R., Li, W. H., Fan, C., Wang, S. Q., Wu, F. Z., et al. (2019). Wheat straw increases the defense response and resistance of watermelon monoculture to Fusarium wilt. BMC Plant Biol. 19:551. doi: 10.1186/s12870-019-2134-y
Tedersoo, L., Bahram, M., Polme, S., Koljalg, U., Yorou, N. S., Wijesundera, R., et al. (2014). Fungal biogeography. Global diversity and geography of soil fungi. Science 346:1256688. doi: 10.1126/science.1256688
Wang, Y. R., Kang, Z. C., Zhu, G. P., Wang, Y., Qi, G. Q., and Yu, H. W. (2022). Optimization of low temperature resistant corn stalk degrading bacterial community and its effect. Acta Agric. Zhejianggensis 34, 2720–2727. doi: 10.3969/j.issn.1004-1524.2022.12.15
Wang, Y. L., Zhang, X. L., Lou, Z. X., An, X. Y., Xue, L., Jiang, X. B., et al. (2023). The effects of adding exogenous lignocellulose degrading bacteria during straw incorporation in cold regions on degradation characteristics and soil indigenous bacteria communities. Front. Microbiol. 14:1141545. doi: 10.3389/fmicb.2023.1141545
Wu, X. F., Yang, J. J., Ruan, H., Wang, S. N., Yang, Y. R., Naeem, I., et al. (2021). The diversity and cooccurrence network of soil bacterial and fungal communities and their implications for a new indicator of grassland degradation. Ecol. Indic. 129:107989. doi: 10.1016/j.ecolind.2021.107989
Xia, Q., Rufty, T., and Shi, W. (2020). Soil microbial diversity and composition: links to soil texture and associated properties. Soil Biol. Biochem. 149:107953. doi: 10.1016/j.soilbio.2020.107953
Yang, H., Hu, J., Long, X., Liu, Z., and Rengel, Z. (2016). Salinity altered root distribution and increased diversity of bacterial communities in the rhizosphere soil of Jerusalem artichoke. Sci. Rep. 6:20687. doi: 10.1038/srep20687
Yue, X. S., Liu, X., Wang, F., Shen, C. W., and Zhang, Y. (2023). Contrasting effects of organic materials versus their derived biochars on maize growth, soil properties and bacterial community in two type soils. Front. Microbiol. 14:1174921. doi: 10.3389/fmicb.2023.1174921
Zaura, E., Keijser, B. J., Huse, S. M., and Crielaard, W. (2009). Defining the healthy "core microbiome" of oral microbial communities. BMC Microbiol. 9:259. doi: 10.1186/1471-2180-9-259
Zgadzaj, R., Garrido-Oter, R., Jensen, D. B., Koprivova, A., Schulze-Lefert, P., and Radutoiu, S. (2016). Root nodule symbiosis in Lotus japonicus drives the establishment of distinctive rhizosphere, root, and nodule bacterial communities. Proc. Natl. Acad. Sci. U. S. A. 113, E7996–E8005. doi: 10.1073/pnas.1616564113
Zhan, A., Zou, C. Q., Ye, Y. L., Liu, Z. H., Cui, Z. L., and Chen, X. P. (2016). Estimating on-farm wheat yield response to potassium and potassium uptake requirement in China. Field Crop Res. 191, 13–19. doi: 10.1016/j.fcr.2016.04.001
Zhang, S., Han, S., Gao, J., Yu, X., and Hu, S. (2023). Low-temperature corn straw-degrading bacterial agent and moisture effects on indigenous microbes. Appl. Microbiol. Biotechnol. 107, 5241–5255. doi: 10.1007/s00253-023-12644-8
Zhang, Z. Y., Liu, D. B., Wu, M. Q., Xia, Y., Zhang, F. L., and Fan, X. P. (2021). Long-term straw returning improve soil K balance and potassium supplying ability under rice and wheat cultivation. Sci. Rep. 11:22260. doi: 10.1038/s41598-021-01594-8
Zhang, B. L., Wu, X. K., Tai, X. S., Sun, L. L., Wu, M. H., Zhang, W., et al. (2019). Variation in Actinobacterial community composition and potential function in different soil ecosystems belonging to the arid Heihe River basin of Northwest China. Front. Microbiol. 10:2209. doi: 10.3389/fmicb.2019.02209
Zhao, Y. C., Fu, W. J., Hu, C. W., Chen, G. Q., Xiao, Z. W., Chen, Y. R., et al. (2021). Variation of rhizosphere microbial community in continuous mono-maize seed production. Sci. Rep. 11:1544. doi: 10.1038/s41598-021-81228-1
Zhao, S. C., Qiu, S. J., Xu, X. P., Ciampitti, I. A., Zhang, S. Q., and He, P. (2019). Change in straw decomposition rate and soil microbial community composition after straw addition in different long-term fertilization soils. Appl. Soil Ecol. 138, 123–133. doi: 10.1016/j.apsoil.2019.02.018
Keywords: straw return, straw-degrading microbes, soil textures, yield, soil microbial community
Citation: Huang Y, Yan Y, Ma Y, Zhang X, Zhao Q, Men M, Huang Y and Peng Z (2024) The effect of low-temperature straw-degrading microbes on winter wheat growth and soil improvement under straw return. Front. Microbiol. 15:1391632. doi: 10.3389/fmicb.2024.1391632
Edited by:
Febri Doni, Padjadjaran University, IndonesiaReviewed by:
Ahmad Ali Pourbabaee, University of Tehran, IranMia Miranti, Padjadjaran University, Indonesia
Copyright © 2024 Huang, Yan, Ma, Zhang, Zhao, Men, Huang and Peng. 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: Mingxin Men, bWVubWluZ3hpbkBzaW5hLmNvbQ==; Yali Huang, aHVhbmd5YWxpMjI5MUAxNjMuY29t; Zhengping Peng, cGVuZ3poZW5ncGluZ0BzaW5hLmNvbQ==
†These authors share first authorship