- 1Institute of Microbiology, Xinjiang Academy of Agricultural Sciences, Urumqi, Xinjiang, China
- 2Key Laboratory of Agricultural Environment in Northwest Oasis of Ministry of Agriculture and Countryside, Urumqi, Xinjiang, China
- 3Xinjiang Laboratory of Special Environmental Microbiology, Urumqi, Xinjiang, China
- 4Institute of Soil, Fertilizer and Agricultural Water Conservation, Xinjiang Academy of Agricultural Sciences, Urumqi, Xinjiang, China
Introduction: The soil microbial community plays an important role in modulating cotton soil fertility. However, the effects of chemical fertilizer combined with organic fertilizer on soil chemical properties, microbial community structure, and crop yield and quality in arid areas are still unclear. This study aimed to explore the effects of different organic fertilizers on soil microbial community structure and diversity and cotton growth and yield.
Methods: High-throughput sequencing was used to study the soil bacteria and fungi in different growth stages of cotton. The field fertilization experiment had five treatments.
Results: The results indicated that the treatments of chemical fertilizer reduction combined with organic fertilizer significantly increased soil available nitrogen and phosphorus in cotton field. There were significant differences in the abundance of the bacterial and fungal communities in the dominant phyla among the treatments. At the phyla level, there were not significantly different in the diversity of bacteria and fungi among treatments. There were significant differences in the composition and diversity of bacterial and fungal communities during the entire cotton growth period (p = 0.001). The rhizosphere bacterial and fungal community structure was significantly affected by soil TK, NH4+, AK, TP, AN, and NO3−. The different fertilization treatments strongly influenced the modular structure of the soil bacterial and fungal community co-occurrence network. A reduction in chemical fertilizer combined with organic fertilizer significantly improved cotton stem diameter and seed yield, and the effect of the biological organic fertilizer on plant growth and yield formation was greater than that of ordinary organic fertilizer.
Discussion: This study provide a scientific and technical basis for the establishment of environmentally friendly green fertilization technology for cotton in arid areas and the promotion of sustainable development of cotton industry.
1 Introduction
Xinjiang is the largest commodity cotton base in China (Niu et al., 2021). Applying chemical fertilizer increases the number of bolls per plant, reduces the bud and boll abscission rate, and increases cotton yields (Zheng et al., 2015). However, the long-term use of chemical fertilizers can lead to a deterioration in soil physicochemical properties, resulting in agricultural non-point source pollution, which can lead to sustained increases in fertilizer use and water, air, and soil pollution (Zhou et al., 2011). Excessive nitrogen fertilizer inputs not only waste resources, but also lead to environmental pollution and soil degradation (Chen J. S. et al., 2021).
It has become an urgent challenge in agricultural production to identify reasonable fertilization measures that reduce the harm caused by excessive soil fertilizer applications. Chemical fertilizer reduction combined with organic fertilizer is an environmental protection fertilization technology (Ye et al., 2020). Appropriate organic fertilizer substitution promotes crop yield, improves soil nutrients and protects soil ecology (Ding et al., 2021; Li et al., 2021a,b; Xiao et al., 2022).
Organic fertilizer combined with chemical fertilizer can reduce nutrient losses in cotton fields (Adeli et al., 2008), improve the fertilizer utilization rate, and increase soil organic matter content (Mensik et al., 2018; Wang et al., 2022). This fertilization method activate soil nutrients, improve soil microbial community structure and diversity (Guo et al., 2019; Zhang et al., 2019; Dong et al., 2020), and change rhizosphere soil enzyme activity. It enhance soil microbial activity, improve soil fertility, and reduce soil salinization and soil-borne diseases (Halihashi et al., 2020; Zhang et al., 2021). Then it improve crop root activity (Wang et al., 2022), antioxidant enzyme activity (Tang et al., 2018), and leaf pigment content (Wang N. et al., 2020; Wang Y. et al., 2020; Wang Q. Q. et al., 2020). It also increases water and fertilizer conservation and ultimately improves crop yield and stress resistance.
One study reported that 80% conventional fertilization combined with moderate amino acid, fulvic acid, and biogas slurry inputs produced the highest cotton yield and benefit (Li et al., 2019), and was significantly higher than that produced by conventional fertilization (Li et al., 2019). Previous studies have shown that when the amount of organic fertilizer was the same, and there was no significant difference in cotton yield between different types of organic fertilizer treatments (Tao et al., 2017; Li et al., 2019; Sun et al., 2020). The application of 70% chemical fertilizer combined with green manure did not reduce the nutrient accumulation, seed cotton yield and fertilizer utilization rate of cotton plant (Li et al., 2019). Furthermore, the agronomic efficiency of chemical fertilizer was significantly higher than that for the 100 and 85% chemical fertilizer treatments, and the seed cotton yield increased by nearly 30% (Sun et al., 2015; Li et al., 2021a,b). The improvement in phosphorus use efficiency after applying biological organic fertilizer is better than that achieved by applying ordinary organic fertilizer. The replacement of 40% chemical fertilizer with bio-organic fertilizer significantly increased cashmere length and yield (Sun et al., 2020). Chemical fertilizer combined with organic fertilizer also significantly increased the number of bacteria and actinomycetes, decreased the number of fungi, and increased urease, catalase, sucrose, and alkaline phosphatase activities (Wang N. et al., 2020; Wang Y. et al., 2020; Wang Q. Q. et al., 2020). Reducing nitrogen fertilization by 30% can reduce soil electrical conductivity and improve soil nutrient content, which will ultimately increase cotton yield and quality (Zhu et al., 2020). It can be seen that chemical fertilizer reduction combined with organic fertilizer is economically feasible. Few studies have investigated the effect of replacing chemical fertilizer with organic fertilizer on microorganisms and cotton growth traits and yield.
In this study, ordinary organic fertilizer, humic acid urea, biological bacterial fertilizer substitution, and conventional fertilization were taken as the research objects. Field experiments were used to investigate the effects of different fertilization regimes on cotton yield, soil nutrients, and the microbial community. In addition, the scientific and economic efficiency of combining chemical fertilizer reduction and organic fertilizer applications was also evaluated. The underlying hypotheses are that (I) the reduction of chemical fertilizer combined with organic fertilizer increases the content of soil nutrient in cotton fields; and (II) affects the community structures of soil microorganism over the whole cotton growth period.
2 Materials and methods
2.1 Field site and experimental design
The experiment was conducted at Xinjiang Academy of Agricultural Sciences in Korla City, Xinjiang, in 2020. The experimental area was located on Baotouhu Farm, Xinjiang Academy of Agricultural Sciences, Korla City, southern Xinjiang (E 85°52′, N 41°41′). The experimental area has a typical arid climate, the average annual rainfall is 56.20 mm. The tested soil was a medium fertility sandy loam soil, the soil organic matter in the 0–30 cm of soil above the plow layer was 10.24 g/kg, available nitrogen was 48.78 mg/kg, available phosphorus was 20.36 mg/kg, available potassium was 139.00 mg/kg, and the pH was 8.20.
The experiment consisted of five treatments (Supplementary Table S1): T1: no fertilizer application; T2: conventional fertilization (CF; NPK dosage: 714 kg/ha: 357 kg N/ha, 207 kg P/ha, and 150 kg K/ha); T3: 60% CF + 12,000 kg/ha organic fertilizer; T4: 46% CF + 428.4 kg/ha humic acid urea; and T5: 73% CF + 225 kg/ha bio-organic fertilizer. Each treatment was replicated three times and the area of each experimental plot was 40.5 m2. Organic fertilizer refers to sheep manure. The nutrient contents of sheep manure were N 0.192%, P2O5 1.16%, K2O 0.82%, humic acid 19.04%, organic matter 56.05%, water content 51.50%. Humic acid urea is a urea containing 0.12% humic acid. The nutrient contents of Bacillus velezensis were N 0.357%, P2O5 0.09%, K2O 0.04%, humic acid 2.00%, organic matter 1.94%. The bio-organic fertilizer is a bacterial fertilizer produced by fermentation of Bacillus velezensis BHZ-29, which we isolated from cotton fields (Zhang et al., 2018).
The cotton variety was Xinluzhong 66. The common organic fertilizer was used as a base fertilizer and applied to the 0–20 cm soil layer before sowing. The chemical fertilizers, humic acid urea, and bio-organic fertilizer were drip applied with water six times throughout the cotton growing season according to their fertilizer requirements (Niu et al., 2021).
2.2 Soil sampling and analysis of physicochemical properties
Soil samples from the 0–20 cm soil layer were collected by soil drilling at the cotton seedling stage (B), bud stage (M), flowering stage (H), boll opening stage (T) in 2020. There were three replicates per treatment, a five-spot sampling method for each treatment. Each replicate collected 500 g soil samples. Plant residues and stones were removed from the soil samples using a 2 mm mesh sieve. Some of the samples were placed in a 4°C refrigerator for the determination of soil microbial diversity, whereas the others were stored under natural air drying to determine the soil physicochemical properties. The SPAD (Soil and Plant Analyzer Development) value of cotton leaves was measured at bud stage, flowering stage and boll stage (Chen et al., 2011). The agronomic traits, yield, and yield related factors were measured at the cotton boll opening stage and the harvest stage. Ten cotton plants were randomly selected from each plot before boll opening. The number of fruit branches per plant, plant height, and boll number per plant were measured. The yield was measured after drying and the seed cotton yield for the plot was calculated. The lint yield and lint percentage were calculated after rolling and the cotton quality index was determined by the Cotton Quality Supervision and Testing Center of the Ministry of Agriculture and Rural Areas (Wan et al., 2018; Chen et al., 2020).
Soil moisture content (WCR) was determined by the drying method. Soil organic matter (OM) was determined using potassium dichromate dilution heat method, total nitrogen (TN) was quantified by a Kjeldahl nitrogen analyzer (multi N/C 2100 S, Analytik Jena, Jena, Germany), alkaline nitrogen (AN) was measured using the diffusion method, soil ammonium nitrogen (NH4+) and nitrate (NO3-) was determined by flow analyzer, total phosphorus (TP) was determined by the acid-soluble-molybdenum-antimony colorimetric method, effective phosphorus (AP) was determined using sodium bicarbonate extraction-molybdenum antimony antispectrophotometry, total potassium (TK) was measured using an NaOH alkali fusion-flame photometer, and available potassium (AK) was determined using the ammonium acetate extraction-flame photometric method (Ma et al., 2022). Each analysis was repeated three times.
2.3 Extraction and sequencing of soil microbial DNA
The soil microbial DNA was extracted according to the Power Soil DNA Isolation Kit instructions (MoBio Laboratories, Carlsbad, CA, United States). The extracted DNA was detected by 1% agarose gel electrophoresis and spectrophotometry, and the qualified samples were stored at −20°C until needed. Primers 338F (5’-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5’-GGACTACHVGGGTWTCTAAT-3′) were used to amplify the V3–V4 region of bacterial 16S rRNA gene (Jones et al., 2021) and the primers for ITS gene sequencing of the fungi were ITS1F (5’-CTTGGTCATTTAGAGAAGTAA-3′) and ITS2 (5’-TGCGTTCTTCATCGATGC-3′) (Li S. et al., 2020; Li Y. et al., 2020). The above primers with barcode sequences were synthesized for the PCR amplification procedure. The PCR products were detected by 1.5% agarose gel electrophoresis, purified using a TIANgel nucleic acid purification kit (Tiangen, Beijing, China) and then used to construct a microbial diversity sequencing library. Paired-end sequencing was performed using the Illumina MiSeq high-throughput sequencing platform by Shanghai Majorbio Bio-pharm Technology Co., Ltd. (Shanghai, China).
2.4 Analysis method for the sequencing data
To improve the accuracy and reliability of the information analysis, the offline data were split into samples according to their Barcode sequence by QIIME1 (v 1.8.0) software. The data were filtered and spliced by Pear (v 0.9.6) software and the chimera sequences are removed using the Uchime method according to the Gold Database. Finally, the Vsearch (v 2.7.1) software UPARSE algorithm was used to cluster the high-quality sequences based on 97% consistency and the Silva128 database was used to annotate the OTUs (Operational Taxonomic Units) using the RDP Classifier algorithm. Then, the community composition of each sample was analyzed at the phylum, class, order, family, genus, and species levels to obtain the species composition and relative abundance at each taxonomic level (Shi et al., 2022). The α-diversity of each sample was calculated using Mothur 1.45.3 (Barouillet et al., 2022) and Excel 2007 (Microsoft, Redmond, WA, United States) software and included the Chao1 value, the Shannon index, and the Simpson index. Rarefaction curves were obtained from the ratio of the number of OTUs in the sample to the effective reads.
2.5 Data processing
Microsoft Excel 2007 was used for data pre-processing and presentation, and SPSS 20.0 (IBM Corp, Armonk, NY, United States) software was used for the one-way analysis of variance (ANOVA) and the Duncan method (α = 0.05) for multiple comparisons and difference testing. The data are shown as the mean ± standard error. A principal coordinate analysis (PCoA) and a permutation multi-factor analysis of variance (Adonis) based on the sample Bray-Curtis distance (OTU level) were performed by R Studio (version 4.0.3), and a dominant OTU heat map and a correlation heat map were constructed. Canoco 5.1 was used to draw a redundancy analysis diagram (redundancy analysis, RDA) (Dalu et al., 2021). Finally, the Spearman’s correlation coefficient was used to calculate the correlations among OTUs (relative abundance ≥0.1%) and OTUs with correlation coefficients of r ≥ 0.6 and p < 0.01 were selected as data sources to construct the co-occurrence network for soil bacteria using Gephi (version 9.2) (Tang et al., 2018).
3 Results
3.1 Soil properties in response to chemical fertilizer reduction combined with organic fertilizer
It can be seen from Supplementary Table S2 that the treatments of chemical fertilizer reduction combined with organic fertilizer increased the contents of soil organic matter (OM), total nitrogen (TN), total phosphorus (TP), total potassium (TK), available nitrogen (AN), available phosphorus (AP) and available potassium (AK) in cotton field. Compared with T1, T3, T4 and T5 treatments significantly increased the content of soil available nitrogen and available phosphorus in cotton fields. The content of soil available nitrogen and available phosphorus in T3 treatment was the highest, reaching 32.33 mg / kg and 32.35 mg/kg. The trend of soil TN, TP, AN and AK in different fertilization treatments was T3 > T5 > T2 > T4 > T1. The content of SOM in each growth period of T3, T4 and T5 treatments showed an increasing trend. Except for T5, the total nitrogen of each treatment increased first and then decreased, and the total nitrogen of each treatment at flowering stage was the highest. Total phosphorus changed little in each period, and total potassium showed an increasing trend. Available nitrogen, phosphorus and potassium showed a downward trend.
3.2 Soil microbial alpha diversity
The 16S rRNA gene and ITS sequencing results showed that 2,616,391 and 3,812,536 original sequences were obtained from 60 samples. The samples contained 34,107–73,318 bacterial and 38,032–74,915 fungal sequences. The minimum numbers of sequences (34,107 and 38,032, respectively) were used as the sampling depth for the bacteria and fungi. The number of OTUs after flattening was 2,685–3,576 and 183–372 for bacteria and fungi, respectively (Table 1). There was no significant difference in the Chao1 indexes, Shannon indexes, and Simpson indexes for the rhizosphere bacterial and fungal communities among treatments (Table 2), indicating that chemical fertilizer reduction combined with organic fertilizer did not significantly change soil microbial α-diversity. The ACE, Chao and Shannon indexes of soil bacteria in the boll opening stage decreased significantly, and the ACE, Chao and Simpson indexes of soil fungi in the flowering and boll opening stages increased significantly.
3.3 Soil microbial community composition and structure
Among the five fertilization treatments, organic fertilizer combined with chemical fertilizer treatments (T3, T5) did not change the community structure of soil bacteria and fungi, but affected the bacterial abundance at different species classification levels (Figure 1). Among the five fertilization treatments, Proteobacteria, Actinobacteriota, Chloroflexi, Acidobacteriota, Gemmatimonadota, Firmicutes, Bacteroidota, and Ascomycota, Mortierellomycota, Basidiomycota, and Chytridiomycota were the dominant phyla in the bacterial and fungal communities, respectively (Figures 1A,D). Alphaproteobacteria, Actinobacteria, Gammaproteobacteria, Acidimicrobiia, Anaerolineae, Vicinamibacteria, Chloroflexia, Gemmatimonadetes, Bacilli, Dehalococcoidia, Thermoleophilia, Bacteroidia. KD4-96 was the dominant class of bacterial community (Figures 1B,E). The class with relative abundances greater than 2% in the T1, T2, T3, T4, and T5 treatments accounted for 73.11, 74.62, 73.11, 74.48, 73.76, 96.19, 95.91, 97.02, 96.69, and 92.94% of the total bacterial and fungal class, respectively (Figures 1B,E). With regards to the dominant genera of the bacterial and fungal community (relative abundance >1%), the T3 and T5 treatments increased the relative abundance of norank_f__norank_o__Vicinamibacterales, norank_f__Vicinamibacteraceae, norank_f__norank_o__norank_c__MB-A2-108 and Preussia compared to T1 and T2, but the norank_f__Gemmatimonadaceae and Mortierella relative abundances decreased (Figures 1C,F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1295722/fmicb-14-1295722-HTML/image_m/fmicb-14-1295722-g001.jpg)
Figure 1. Relative abundances of bacterial (A–C) and fungal (D–F) taxa at the phylum, class, genus level, respectively. T1: no fertilizer application (CK); T2: Conventional fertilization (CF; NPK dosage: 714 kg ha-1: 357 kg N ha-1, 207 kg P ha-1, and 150 kg K ha-1); T3: 60% CF + 12,000 kg ha organic fertilizer (CFO; NPK dosage: 428.4 kg ha-1: 214.2 kg N ha-1, 124.2 kg P ha-1, and 90 kg K ha-1); T4: 46% CF + 428.4 kg ha humic acid urea (CFH; NPK dosage: 714 kg ha-1: 357 kg N ha-1, 207 kg P ha-1, and 150 kg K ha-1); T5: 73% CF + 225 kg ha bio-organic fertilizer (CFB; NPK dosage: 617.55 kg ha-1: 260.55 kg N ha-1, 207 kg P ha-1, and 150 kg K ha-1).
3.4 Dynamics of the soil microbial communities
The PCoA analysis showed that there are seasonal differences in the bacterial and fungal community structure among the samples (Figures 2A,B). The PC1 axis for bacteria and fungi explained 23.97 and 28.26% of the difference, respectively, and the PC2 axis explained 19.34 and 17.95%, respectively. Together, they explained 43.31 and 46.21% of the variability, respectively (Figures 2A,B). Figure 2 also showed that the dispersion degree of the samples at the bud, flowering, and boll stages was low. However, the dispersion degree was high for the samples at the seedling stage. There were significant differences in the community structures of the bacteria and fungi over the whole cotton growth period. The bacteria samples were highly clustered on the PC1 axis, which showed that the rhizosphere bacteria were highly similar in the samples, but the fungi were very different at the flowering stage.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1295722/fmicb-14-1295722-HTML/image_m/fmicb-14-1295722-g002.jpg)
Figure 2. Principal coordinate analysis (PCoA) plots of bacterial (A) and fungal (B) community composition at the OTU level. T1, T2, T3, T4, and T5 are as defined in the Figure 1 legend. B: seedling stage: M: bud stage: H: flowering stage: T: boll opening stage.
The cluster analysis (Supplementary Figure S1) showed that there were significant differences in bacterial and fungal community composition among the samples at the bud, flowering, and boll stages, while the samples at the seedling stage were clustered on the same branch. The growth period had significant effects on the bacterial floras Actinobacteria, Gammaproteobacteria, Acidimicrobiia, Chloroflexia, Gemmatimonadetes, Dehalococcoidia, and Bacteroidia, and the fungal floras Leotiomycetes, Sordariomycetes, Mortierellomycetes, Eurotiomycetes, and Pezizomycetes were dominant in the five treatments (Supplementary Figure S2).
3.5 Correlation between bacterial and fungal abundance and soil properties
The Spearman’s correlation heatmap analysis indicated that the soil physicochemical factors affecting the soil bacterial and fungal communities were TK, NH4, AK, TP, and AN and TP, AN, NO3, and TK, respectively (Supplementary Figure S3). The effect of soil physicochemical factors on Bacillus was weak and only TK, NH4, TP, and AN were significant factors. Pseudounas was sensitive to most soil physicochemical factors, except for NH4+ and NO3-. The physical and chemical factors OM and NH4+ had significant effects on Microascus, and WCR and NO3 had significant effects on Fusarium, and NO3-had significant effects on Tricharina (Supplementary Figures S3A,B).
The RDA analysis showed that soil physicochemical factors explained 35.50% of the differences in the community structure of the dominant bacteria and fungi, and that RDA1 and RDA2 explained 29.09, 21.97 and 6.41%, 9.77% of the bacterial and fungal differences, respectively (Figure 3). Total phosphate (R2 = 0.75, p = 0.0000), TK (R2 = 0.47, p = 0.0000), and AN (R2 = 0.59, p = 0.0000) for the bacterial community and TP (R2 = 0.17, p = 0.0012), TK (R2 = 0.10, p = 0.0139), and AN (R2 = 0.15, p = 0.0022) for the fungal community reached significant levels (Figures 3A,B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1295722/fmicb-14-1295722-HTML/image_m/fmicb-14-1295722-g003.jpg)
Figure 3. Redundancy analysis (RDA) of bacterial (A) and fungal (B) communities with soil chemical properties. WCR, electrical conductivity; OM, soil organic matter, NO3, electrical conductivity, NH4, soil organic matter.TN. total nitrogen: TP. total phosphorus: TK, total potassium; AN, total nitrogen, AP. total phosphorus: AK, total potassium. T1, T2, T3, T4, and T5 are as defined in the Figure 1 legend. B: seedling stage; M: bud stage, H: flowering stage. T. boll opening stage.
3.6 Soil microbial network structure and composition analysis
The soil bacterial and fungal co-occurrence networks under the T1, T2, T3, T4, and T5 treatments were composed of 251, 253, 247, 254, 361, and 84, 85, 87, 79, 84 nodes, respectively, and 3,698, 3,403, 3,071, 3,776, and 4,291, 308, 337, 332, 305, 297 had highly significantly positively correlated edges, with average degrees of 28.337, 26.901, 24.866, 29.732, 32.881, and 7.333, 7.929, 7.632, 7.722, 7.071, respectively (Figures 4A,B). The modular indexes were 3.249, 5.458, 2.793, 3.498, 2.955, and 1.798, 0.638, 0.757, 0.316, 1.965, respectively. A modularity index >0.4 indicates that the co-occurrence network has a modular structure.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1295722/fmicb-14-1295722-HTML/image_m/fmicb-14-1295722-g004.jpg)
Figure 4. The co-occurrence network of soil bacterial and fungal community composition in a cotton field with different fertilization treatments on the phylum level. (A) bacterial community, (B) fungal community. Each node denotes a bacterial or a fungal OTU (defined at a 97% similarity level); each edge linking two nodes represents a positive (pink line) or negative (black line) relationship. OTUS are colored by different phylum. The size of each node is proportional to the number of connections. A connection between two nodes is a statistically significant (p < 0.01) and strong (r > 0.60) correlation. The percentage of positive links in every network: A 90.08%, B 73.87%.
Figure 4A also indicated that the bacterial co-occurrence networks for T1, T2, and T4 were similar, whereas the bacterial co-occurrence networks for T5 were relatively complex. The fungal co-occurrence networks for T1, T4, and T5 were similar, as were the fungal co-occurrence networks for T2 and T3 (Figure 4B).
Supplementary Table S3 showed that the proportions of the bacterial and fungal groups in the co-occurrence network were different in each treatment. The Proteobacteria and Actinobacteriota proportions were the largest in the bacterial co-occurrence network and the Ascomycota and Basidiomycota proportions were largest in the fungal co-occurrence network. The bacterial groups in each treatment were basically the same, but the number of fungal groups was lower in T3. In addition. Proteobacteria, Actinobacteriota, Chloroflexi, Bacteroidota, Acidobacteriota, Firmicutes, and Myxococcota were the dominant genera in the soil bacterial co-occurrence network and Pseudeurotium, Metarhizium, Talaromyces, Trichoderma, Preussia, Chaetomium, and Mortierella were the dominant genera in the soil fungal co-occurrence network.
3.7 Cotton agronomic traits, quality and yield
The T5 treatment produced the highest yield (Table 3). The four fertilization treatments significantly increased the boll number per plant and the seed cotton yield compared to T1 (Table 3). Compared to T2, the stem diameter and seed cotton yield significantly increased under the chemical fertilizer reduction combined with organic fertilizer treatments, but there was no significant difference in plant height, leaf number, SPAD, single boll weight, and lint percentage. This showed that organic fertilizer promoted cotton plant growth and boll formation, but had no significant effect on cotton yield per plant. The raw cotton yields of the chemical fertilizer reduction combined with organic fertilizer treatments were in the order T5 > T4 > T3 and the differences among the three treatments were significant. The raw cotton yields of T4 and T5 were 5.27 and 12.10% higher than those of T2, respectively, when T2 was reduced by 50 and 54%, combined with 428.4 kg hm−2 humic acid urea (T4) and 225 L hm-2 Bacillus velezensis fertilizer (T5), respectively, and there was no significant effect on cotton quality.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1295722/fmicb-14-1295722-HTML/image_m/fmicb-14-1295722-t003.jpg)
Table 3. Effects of different fertilizer treatments on agronomic characters, quality and yield of cotton.
4 Discussion
4.1 Effects of fertilization on soil physicochemical properties
The proportion of soil available nutrients compared to the total soil nutrients is one of the important indicators used to evaluate the effectiveness of soil nutrients (Zhao and Huang, 2022). Reasonable combinations of organic and inorganic fertilizers can improve the soil physicochemical properties. In this study, we have demonstrated that reductions in chemical fertilizer combined with organic fertilizer could increase soil nutrients, such as AP, SOM and TN, which was consistent with previous studies on agricultural soils (Lu et al., 2009; Hossain et al., 2021; Yang W. N. et al., 2022). The increase in AP content by the fertilizer combined with organic fertilizer suggests that the organic fertilizer provided a large amount of soil organic matter (SOM), and then, organic matter decomposes under the action of microorganisms to produce a large number of organic acids, which can activate soil phosphorus, thereby increasing the soil AP content (Alori et al., 2017).
4.2 Effects of fertilization on the diversity and composition of soil microbial community
In microbial diversity studies, the greater the Shannon index, the higher the microbial community evenness; the higher the Chao1 and Ace indices, the higher the richness of the microbial communities (Hartmann et al., 2006). We also found that chemical fertilizer reduction combined with organic fertilizer increased soil bacterial richness and diversity, and reduced fungal richness and diversity (Table 2), which was consistent with previous studies that showed increased richness and diversity of the bacterial community, and decreased richness and diversity of the fungal community (Liu et al., 2021; Jin et al., 2022).
Microorganisms are an important part of the soil ecosystem. Their community composition and diversity show certain dynamic changes at different growth stages and under different nutritional statuses, fertilization measures, and tillage modes (Yang and Crowley, 2000; Wang et al., 2016; Liu et al., 2020; Iqbal et al., 2022). They are important indicators for evaluating soil fertility and are related to the occurrence of plant soil-borne diseases, continuous cropping obstacles, and soil acidification and salinization (Yao and Wu, 2010; Shen et al., 2016; Gao et al., 2021; Song et al., 2022). We found that there was little differentiation in the composition of bacterial or fungal communities among different fertilization treatments (Figures 2A,B). The reason might be that the soil organic matter content and microbial community diversity of long-term continuous cropping cotton field were low, and the organic matter content in fertilization treatment was not high.
4.3 Correlation between microorganism and soil physicochemical properties
Soil can directly provide the necessary environmental conditions for the survival of microorganisms (Santoyo et al., 2017). Organic fertilizer application can not only increase the total amount of soil nutrients such as soil total organic carbon, total nitrogen, and total phosphorus but can also increase the content of available soil nutrients (alkali nitrogen, available phosphorus, available potassium) (Liu et al., 2014; Qaswar et al., 2020). TP, TK, AN is significantly associated with the bacterial and fungal communities and is one of the main factors promoting changes in the soil microbial community (Dong et al., 2014; Hu et al., 2019; Chen Y. et al., 2021; Li et al., 2021a,b; Zhao et al., 2021). Proteobacteria and Ascomycota is a dominant taxon in the bacterial and fungal community, separately and is positively correlated with the soil nitrogen, phosphorus and potassium pool; therefore, TP, TK, AN has been proven to drive changes in the soil bacterial and fungal community (Landesman and Dighton, 2010; Zhong et al., 2010; Ducousso-Détrez et al., 2022). In the present study, TP, TK, and AN were found to be the key determinants driving changes in the bacterial and fungal community in this study (Figures 3A,B; Supplementary Figures S3A,B), which was consistent with previous studies on soils (Landesman and Dighton, 2010; Zhong et al., 2010; Ducousso-Détrez et al., 2022).
Another study reported that SOM were the major driver on structuring soil microbial communities across land uses, while soil bacterial communities were more sensitive to variations in SOM and geochemical characteristics compared with fungi (Bahadori et al., 2022). This is inconsistent with our results and suggests that the effect of chemical fertilizer reduction combined with organic fertilizer on microbial community structure is not the only factor. This makes it difficult to assess the full effect of their nutritional potential on the soil microbial community in a short-term study period (Zhang et al., 2021). Therefore, the short-term effect of chemical fertilizer reduction combined with organic fertilizer on soil microbial community may be smaller than chemical fertilizer.
4.4 Co-occurrence networks for soil microorganisms in cotton fields
The network complexity became gradually more complicated as soil nutrient levels increased (Liao et al., 2020; Wang et al., 2021; Ye et al., 2021), indicating that the increase in the co-occurrence network complexity in the T5 treatment was closely related to the increases in SOM, TN, TP, and AP contents (Supplementary Table S2). The proportion of Proteobacteria in the T5 co-occurrence network was higher than that for the other treatments, but the Patescibacteria, Chloroflexi, and Planctomycetota proportions were lower than in the other treatments in our study (Supplementary Table S3). Proteobacteria are mainly distributed in the upper humus layer of farmland soil or the rhizosphere and are enriched in environments with high levels of active organic carbon (Pinhassi and Berman, 2003; Sagova-Mareckova et al., 2016). In contrast, Chloroflexi are oligotrophic bacteria with slow growth characteristics (Feng et al., 2022) and are ubiquitous in nutrient-poor soils (Wang N. et al., 2020; Wang Y. et al., 2020; Wang Q. Q. et al., 2020). When the substrate concentration of the microbial environment increased, the nutrient-rich bacteria replaced the oligotrophic bacteria and colonize the nutrient-rich environment, which strongly suggested that the T5 treatment improved the growth of nutrient-rich bacteria and restricted the growth of oligotrophic bacteria (Wang N. et al., 2020; Wang Y. et al., 2020; Wang Q. Q. et al., 2020).
The dominant groups are considered to be important drivers of microbial community structure and function (Banerjee et al., 2018). Zhang et al. (2017) found that Acinetobacter and Pseudomonas were the core functional bacteria in acid red soil paddy fields and that their relative abundances increased by 5.7–10 times over 31 years, indicating that bacterial communities experienced ecological succession, and that dominant bacteria occupied specific niches and had specific functions. In this study, Proteobacteria, Actinobacteriota, Chloroflexi, and Pseudeurotium, and Metarhizium and Talaromyces were the core genera of soil bacteria and fungi in the co-occurrence networks.
4.5 Effects of fertilization on agronomic traits and cotton yield
Combined applications of chemical fertilizer and organic fertilizer can regulate the release and intensity of soil and fertilizer nutrients, which means that crops can obtain stable and balanced amounts of nutrients at all growth stages (Zhao et al., 2016; Muhammad et al., 2020). This study showed that the chemical fertilizer reduction combined with organic fertilizer treatments increased cotton yields by different degrees, which is consistent with previous studies on cotton field soils (Bai et al., 2014). Previous studies have shown that the ratio of chemical fertilizer to organic fertilizer is related to soil fertility and climatic conditions (Ye et al., 2020). When soil fertility is high, increasing the proportion of organic fertilizer can promote reproduction by microorganisms, improve soil structure, and increase crop yields (Ye et al., 2020). Bioorganic fertilizers can replace 23–52% of nitrogen fertilizers without causing a loss of yield (Rose et al., 2014).
In this study, the effect of fertilizer reduction combined with bacterial fertilizer on yield increase was clear, indicating that the nutrient release law for organic fertilizer was consistent with the nutrient demand law for cotton growth. At the same time, crop yield was significantly correlated with soil nutrients, microorganisms, and related enzyme activities (Jiang et al., 2017). Therefore, reductions in chemical fertilizer and reasonable applications of organic fertilizer could effectively control the number of soil microorganisms, improved soil enzyme activity and soil fertility. At present, soil nutrient imbalances in cotton fields caused by unreasonable fertilization is common in Xinjiang (Gong et al., 2012; Yang J. Y. et al., 2022). The fertilization structure can be adjusted and optimized to reduce chemical fertilizer application rates through the combined use of chemical fertilizer and organic fertilizer. In this study, reducing chemical fertilizer by 30% combined with 12,000 kg·hm−2 common organic fertilizer or 225 kg·hm−2 bio-organic fertilizer produced the highest yields, but the continuous application of organic fertilizer changed soil fertility. Therefore, future research should investigate the most appropriate proportion of organic fertilizer to apply to a field when combined with a chemical fertilizer application.
5 Conclusion
This study indicated that the reduction of chemical fertilizer combined with organic fertilizer significantly increased the content of soil available nitrogen and phosphorus in cotton fields, and that total and available nitrogen, phosphorus and potassium contents in the chemical fertilizer reduction combined with organic fertilizer treatments were basically stable in the cotton fields. The application of chemical fertilizer reduction combined with organic fertilizer significantly affected the community structures of the bacteria and fungi over the whole cotton growth period, without significantly changing the soil microbial alpha diversity. The different fertilization treatments strongly influenced the modular structure of the soil bacterial and fungal community co-occurrence network. A reduction in chemical fertilizer combined with organic fertilizer significantly improved cotton stem diameter and seed yield, and the effect of the biological organic fertilizer on plant growth and yield formation was greater than that of ordinary organic fertilizer.
Data availability statement
The datasets presented in this study can be found in online repositories. The sequences from rhizosphere bacterial samples are available at NCBI SRA: BioProject PRJNA876230.
Author contributions
YS: Formal analysis, Investigation, Writing – original draft. XN: Funding acquisition, Investigation, Resources, Writing – original draft, Writing – review & editing. BC: Investigation, Resources, Writing – original draft. SP: Investigation, Resources, Writing – review & editing. HM: Data curation, Formal analysis, Investigation, Writing – original draft. PL: Data curation, Resources, Writing – review & editing. GF: Formal analysis, Investigation, Writing – review & editing. XM: Formal analysis, Software, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Open Fund Project of Key Laboratory of Agricultural Environment of Northwest Oasis, Ministry of Agriculture and rural Affairs (XBLZ-20205), the National Nature Science Foundation of China (31860024), the Project of Renovation Capacity Building for the Young Sci-Tech Talents Sponsored by Xinjiang Academy of Agricultural Sciences (xjnkq-2019014), the Major Scientific and Technological Project of Xinjiang Uygur Autonomous Region of China (2020A01002), the National Key Research and Development Program of China (2017YFD0201903, 2018YFD0800804).
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.2023.1295722/full#supplementary-material
References
Adeli, A., Shankle, M. W., Tewolde, H., Sistani, K. R., and Rowe, D. E. (2008). Nutrient dynamics from broiler litter applied to no-till cotton in an upland soil. Agron. J. 100, 564–570. doi: 10.2134/agronj2007.0224
Alori, E. T., Glick, B. R., and Babalola, O. O. (2017). Microbial phosphorus solubilization and its potential for use in sustainable agriculture. Front. Microbiol. 8:971. doi: 10.3389/fmicb.2017.00971
Bahadori, M., Wang, J. T., Shen, J. P., Lewis, S., Rezaei Rashti, M., and Chen, C. (2022). Soil organic matter and geochemical characteristics shape microbial community composition and structure across different land uses in an Australian wet tropical catchment. Land Degrad. Dev. 33, 817–831. doi: 10.1002/ldr.4174
Bai, L., Li, J. H., Chu, G. X., and Hua, R. (2014). Effects of organic manure application combined with chemical fertilizers on nutrients absorption and nitrogen efficiency of cotton in Xinjiang China. Agric. Res. Arid Areas 32, 143–148.
Banerjee, S., Schlaeppi, K., and Van Der Heijden, M. G. A. (2018). Keystone taxa as drivers of microbiome structure and functioning. Nat. Rev. Microbiol. 16, 567–576. doi: 10.1038/s41579-018-0024-1
Barouillet, C., Vasselon, V., Keck, F., Millet, L., Etienne, D., Galop, D., et al. (2022). Paleoreconstructions of ciliate communities reveal long-term ecological changes in temperate lakes. Sci. Rep. 12, 1–12. doi: 10.1038/s41598-022-12041-7
Chen, Y., Li, S., Liu, N., He, H., Cao, X., 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. 28, 23036–23047. doi: 10.1007/s11356-020-12203-y
Chen, B. Y., Ma, X. W., Yang, T., Yao, Q. Q., Xue, X. R., Wang, B., et al. (2011). Research on the relationship between the accuracy of SPAD value and sample number in growing period of cotton. Sci. Agric. Sin. 44, 4748–4755. doi: 10.3864/j.issn.0578-1752.2011.22.024
Chen, J. S., Sun, W. H., Wang, G. S., Abubakar, S. A., and Gao, Y. (2021). Effects of different nitrogen application rates on soil water stable aggregates and N(2)O emission in winter wheat field. Chin. J. Appl. Ecol. 32, 3961–3968. doi: 10.13287/j.1001-9332.202111.024
Chen, L. Y., Zhang, L. F., Lu, Z. Y., Xian, F., Zhang, J. Z., and Zhang, X. Q. (2020). Effect of pipeline buried depth and water amount of trace irrigation on cotton quality and yield. Water Saving Irrigation 3, 7–11. doi: 10.3969/j.issn.1007-4929.2020.03.002
Dalu, T., Dlamini, P., Wasserman, R. J., Mokgoebo, M. J., Mutshekwa, T., Dondofema, F., et al. (2021). Effects of environmental variables on littoral macroinvertebrate community assemblages in subtropical reservoirs. Chem. Ecol. 37, 419–436.
Ding, W. T., Fang, J. J., Wu, X. P., Zhang, J. Z., Zhang, J. M., Zhang, J. Z., et al. (2021). Effects of different ratios of organic fertilizers instead of chemical fertilizers on black soil microbiological properties and spring wheat yield and quality. Soil Fert. Sci. China 2, 44–52. doi: 10.11838/sfsc.1673-6257.20452
Dong, W. Y., Zhang, X. Y., Dai, X. Q., Fu, X. L., Yang, F. T., Liu, X. Y., et al. (2014). Changes in soil microbial community composition in response to fertilization of paddy soils in subtropical China. Appl. Soil Ecol. 84, 140–147. doi: 10.1016/j.apsoil.2014.06.007
Dong, W., Zhang, Q., Luo, T., and Wang, H. P. (2020). Effects of continuous application of different organic fertilizers on soil quality. Chin. Agric. Sci. Bull. 36, 106–110.
Ducousso-Détrez, A., Fontaine, J., Lounes-Hadj Sahraoui, A., and Hijri, M. (2022). Diversity of phosphate chemical forms in soils and their contributions on soil microbial community structure changes. Microorganisms 10:609. doi: 10.3390/microorganisms10030609
Feng, H., Pan, H., Li, C., and Zhuge, Y. (2022). Microscale heterogeneity of soil bacterial communities under long-term fertilizations in fluvo-aquic soils. Soil Ecol. Lett. 4, 337–347. doi: 10.1007/s42832-021-0121-3
Gao, Z., Hu, Y., Han, M., Xu, J., Wang, X., Liu, L., et al. (2021). Effects of continuous cropping of sweet potatoes on the bacterial community structure in rhizospheric soil. BMC microbiol. 21, 1–13.
Gong, L., Ran, Q. Y., and Han, L. (2012). Enzymatic activities and related physicochemical properties in soils of continuous cropping cotton fields within a typical oasis in upper reaches of Tarim River. Bull. Soil Water Conserv. 32, 36–42.
Guo, Z., Han, J., Li, J., Xu, Y., and Wang, X. (2019). Effects of long-term fertilization on soil organic carbon mineralization and microbial community structure. PLoS One 14:e0211163. doi: 10.1371/journal.pone.0211163
Halihashi, Y. B. T., Li, Q. J., and Zhang, Y. (2020). Effects of cheimcal feritilzers appilcaiton combined with manure on yield and nutrient absorption in cotton field. Xinjiang Agric. Sci. 57, 1049–1056. doi: 10.6048/j.issn.1001-4330.2020.06.008
Hartmann, M., Hartmann, M., and Widmer, F. (2006). Community structure analyses are more sensitive to differences in soil bacterial communities than anonymous diversity indices. Appl. Environ. Microbiol. 72, 7804–7812. doi: 10.1128/AEM.01464-06
Hossain, M. E., Mei, X., Zhang, W., Dong, W., Yan, Z., Liu, X., et al. (2021). Substitution of chemical fertilizer with organic fertilizer affects soil Total nitrogen and its fractions in northern China. Int. J. Environ. Res. Public Health 18:12848. doi: 10.3390/ijerph182312848
Hu, L., Zi, H., Wu, P., Wang, Y., Lerdau, M., Wu, X., et al. (2019). Soil bacterial communities in grasslands revegetated using Elymus nutans are largely influenced by soil pH and total phosphorus across restoration time. Land Degrad. Dev. 30, 2243–2256. doi: 10.1002/ldr.3414
Iqbal, A., Ali, I., Yuan, P., Khan, R., Liang, H., Wei, S., et al. (2022). Combined application of manure and chemical fertilizers alters soil environmental variables and improves soil fungal community composition and Rice grain yield. Front. Microbiol. 13:856355. doi: 10.3389/fmicb.2022.856355
Jiang, R., Tang, L., Li, M., and Xu, Z. (2017). Response of greenhouse soil microbial community and enzyme activities to combined application of chemical fertilizer reduction with organic fertilizer. Chin. J. Soil Sci. 48, 639–646. doi: 10.19336/j.cnki.trtb.2017.03.19
Jin, N., Jin, L., Wang, S., Li, J., Liu, F., Liu, Z., et al. (2022). Reduced chemical fertilizer combined with bioorganic fertilizer affects the soil microbial community and yield and quality of lettuce. Front. Microbiol. 13:863325. doi: 10.3389/fmicb.2022.863325
Jones, B., Goodall, T., George, P. B., Gweon, H. S., Puissant, J., Read, D. S., et al. (2021). Beyond taxonomic identification: integration of ecological responses to a soil bacterial 16S rRNA gene database. Front. Microbiol. 12:682886. doi: 10.3389/fmicb.2021.682886
Landesman, W. J., and Dighton, J. (2010). Response of soil microbial communities and the production of plant-available nitrogen to a two-year rainfall manipulation in the New Jersey pinelands. Soil Biol. Biochem. 42, 1751–1758. doi: 10.1016/j.soilbio.2010.06.012
Li, S., Deng, Y., Wang, Z., Zhang, Z., Kong, X., Zhou, W., et al. (2020). Exploring the accuracy of amplicon-based internal transcribed spacer markers for a fungal community. Mol. Ecol. Resour. 20, 170–184. doi: 10.1111/1755-0998.13097
Li, Y., Li, Q. J., Zhang, Y., Halihashi, Y. B. T., and Wang, A. L. (2020). Effects of liquid organic fertilizers with reduced N and P fertilizer application on soil nutrient of drip irrigated cotton field under plastic film mulching. Agric. Res. Arid Areas. 38, 99–104. doi: 10.7606/j.issn.1000-7601.2020.02.14
Li, H., Qiu, Y., Yao, T., Han, D., Gao, Y., Zhang, J., et al. (2021a). Nutrients available in the soil regulate the changes of soil microbial community alongside degradation of alpine meadows in the northeast of the Qinghai-Tibet plateau. Sci. Total Environ. 792:148363. doi: 10.1016/j.scitotenv.2021.148363
Li, H., Yang, S., Semenov, M. V., Yao, F., Ye, J., Bu, R., et al. (2021b). Temperature sensitivity of SOM decomposition is linked with a K-selected microbial community. Glob. Chang. Biol. 27, 2763–2779. doi: 10.1111/gcb.15593
Li, Y., Zhang, Y., Halihashi, Y. B. T., and Li, Q. J. (2019). Effects of chemical fertilizer reduction application combined with liquid organic fertilizer on nutrient absorption and yield of cotton (Gossypium hirsutum L.). Xinjiang Agric. Sci. 56, 137–145. doi: 10.6048/j.issn.1001-4330.2019.01.017
Liao, H., Zhang, Y., Wang, K., Hao, X., Chen, W., and Huang, Q. (2020). Complexity of bacterial and fungal network increases with soil aggregate size in an agricultural Inceptisol. Appl. Soil Ecol. 154:103640. doi: 10.1016/j.apsoil.2020.103640
Liu, Z., Liu, J., Yu, Z., Yao, Q., Li, Y., Liang, A., et al. (2020). Long-term continuous cropping of soybean is comparable to crop rotation in mediating microbial abundance, diversity and community composition. Soil Tillage Res. 197:104503. doi: 10.1016/j.still.2019.104503
Liu, K., Yang, Y. L., Yang, J. Z., Wang, Z. X., Tian, Y. Y., Yang, B. K., et al. (2021). Effects of reduced chemical fertilizer combined with organic fertilizer on bacterial community diversity of Tobacco rhizosphere soil in field. Southwest China J. Agri. Sci. 34, 2191–2196. doi: 10.16213/j.cnki.scjas.2021.10.018
Liu, P. P., Zhou, Y., Fu, G. X., Gao, X., Zhang, P., Zhang, L., et al. (2014). Effects of fertilization on crop yield and soil phosphorus availablity based on the returning straw. Nanjing Nongye Daxue Xuebao 37, 27–33. doi: 10.7685/j.issn.1000-2030.2014.05.005
Lu, S. J., Han, X. Z., Zhag, D., Wang, F. J., and Zhu, W. W. (2009). Effect of long-term fertilization on microbial biomass C and N in soybean rhizosphere. Soybean Sci. 28, 495–498.
Ma, L., Song, D., Liu, M., Li, Y., and Li, Y. (2022). Effects of earthworm activities on soil nutrients and microbial diversity under different tillage measures. Soil Tillage Res. 222:105441. doi: 10.1016/j.still.2022.105441
Mensik, L., Hlisnikovsky, L., Pospisilova, L., and Kunzova, E. (2018). The effect of application of organic manures and mineral fertilizers on the state of soil organic matter and nutrients in the long-term field experiment. J. Soils Sediments 18, 2813–2822. doi: 10.1007/s11368-018-1933-3
Muhammad, Q., Huang, J., Waqas, A., Li, D., Liu, S., Zhang, L., et al. (2020). Yield sustainability, soil organic carbon sequestration and nutrients balance under long-term combined application of manure and inorganic fertilizers in acidic paddy soil. Soil Tillage Res. 198:104569. doi: 10.1016/j.still.2019.104569
Niu, X. X., Pu, S., Wu, X. L., Li, P., Yang, T., Ma, H. H., et al. (2021). Effects of the substitution for chemical fertilizer with organic manure on the cotton growth character and yield. Xinjiang Agric. Sci. 58, 2043–2048. doi: 10.6048/j.issn.1001-4330.2021.11.009
Pinhassi, J., and Berman, T. (2003). Differential growth response of colony-forming α-and γ-proteobacteria in dilution culture and nutrient addition experiments from Lake Kinneret (Israel), the eastern Mediterranean Sea, and the Gulf of Eilat. Appl. Environ. Microbiol. 69, 199–211. doi: 10.1128/AEM.69.1.199-211.2003
Qaswar, M., Chai, R., Ahmed, W., Jing, H., Han, T., Liu, K., et al. (2020). Partial substitution of chemical fertilizers with organic amendments increased rice yield by changing phosphorus fractions and improving phosphatase activities in fluvo-aquic soil. J. Soils Sediments 20, 1285–1296. doi: 10.1007/s11368-019-02476-3
Rose, M. T., Phuong, T. L., Nhan, D. K., Cong, P. T., Hien, N. T., and Kennedy, I. R. (2014). Up to 52% N fertilizer replaced by biofertilizer in lowland rice via farmer participatory research. Agron. Sustain. Dev. 34, 857–868. doi: 10.1007/s13593-014-0210-0
Sagova-Mareckova, M., Zadorova, T., Penizek, V., Omelka, M., Tejnecky, V., Pruchova, P., et al. (2016). The structure of bacterial communities along two vertical profiles of a deep colluvial soil. Soil Biol. Biochem. 101, 65–73. doi: 10.1016/j.soilbio.2016.06.026
Santoyo, G., Pacheco, C. H., Salmeron, J. H., and Leon, R. H. (2017). The role of abiotic factors modulating the plant-microbe-soil interactions: toward sustainable agriculture. A review. Span. J. Agric. Res. 15:e03R01. doi: 10.5424/SJAR/2017151-9990
Shen, W., Ni, Y., Gao, N., Bian, B., Zheng, S., Lin, X., et al. (2016). Bacterial community composition is shaped by soil secondary salinization and acidification brought on by high nitrogen fertilization rates. Appl. Soil Ecol. 108, 76–83.
Shi, K., Shao, H., Han, C. X., and Toshmatov, Z. (2022). Diversity of the rhizosphere soil Fungi of the invasive plant (Solanum rostratum Dunal) and the Allelopathic potential of their secondary metabolites. Chin. J. Soil Sci. 53, 548–557. doi: 10.19336/j.cnki.trtb.2021090202
Song, L., Wang, J., Pan, J., Yan, Y., and Niu, S. (2022). Chronic nitrogen enrichment decreases soil gross nitrogen mineralization by acidification in topsoil but by carbon limitation in subsoil. Geoderma 428:116159.
Sun, G. L., Feng, K. Y., Zhao, X. X., Wang, Z. B., and Li, Y. B. (2020). Effects of organic manure replacing chemical fertilizer on cotton growth, development and yield. Xinjiang Agric. Sci. 57, 762–769. doi: 10.6048/j.issn.1001-4330.2020.04.023
Sun, W. Y., Yin, H. J., Xu, J. X., Li, Z. J., Tian, C. Y., and Tang, J. W. (2015). The study of winter green manure on cotton production and soil fertility. Acta Agric. Boreali-Sin. 30, 129–132. doi: 10.7668/hbnxb.2015.S1.023
Tang, H. M., Xiao, X. P., Li, C., Tang, W. G., Guo, L. J., Wang, K., et al. (2018). Effects of different long-term fertilization managements on the physiological characteristics of leaves and yield of rice in double cropping paddy field. J. China Agric. Univ. 23, 60–71. doi: 10.11841/j.issn.1007-4333.2018.11.06
Tao, R., Wakelin, S. A., Liang, Y., and Chu, G. (2017). Organic fertilization enhances cotton productivity, nitrogen use efficiency, and soil nitrogen fertility under drip irrigated field. Agron. J. 109, 2889–2897. doi: 10.2134/agronj2017.01.0054
Wan, H. L., Liu, P. C., Liu, L. T., Zhang, Y. J., Liu, Y. C., Bai, Z. Y., et al. (2018). Effect of moderate drought in the early stage on cotton yield, Fiber quality and water use efficiency. Cotton Sci. 30, 464–472. doi: 10.11963/1002-7807.whlshc.20181113
Wang, N., Feng, K. Y., Nan, H. Y., and Zhang, T. H. (2022). Effects of combined application of organic fertilizer and chemical fertilizer on root characteristics and yield of cotton under different water conditions. Chin. Agric. Sci. 55, 2187–2201. doi: 10.3864/j.issn.0578-1752.2022.11.009
Wang, E., Lin, X., Tian, L., Wang, X., Ji, L., Jin, F., et al. (2021). Effects of short-term rice straw return on the soil microbial community. Agriculture 11:561. doi: 10.3390/agriculture11060561
Wang, N., Nan, H. Y., and Feng, K. Y. (2020). Effects of reduced chemical fertilizer with organic fertilizer application on soil microbial biomass, enzyme activity and cotton yield. Chin. J. Appl. Ecol. 31, 173–181. doi: 10.13287/j.1001-9332.202001.022
Wang, Y., Osman, J. R., and DuBow, M. S. (2020). Bacterial communities on the surface of the mineral sandy soil from the desert of Maine (USA). Curr. Microbiol. 77, 1429–1437. doi: 10.1007/s00284-020-01946-z
Wang, J., Xue, C., Song, Y., Wang, L., Huang, Q., and Shen, Q. (2016). Wheat and rice growth stages and fertilization regimes alter soil bacterial community structure, but not diversity. Front. Microbiol. 7:1207. doi: 10.3389/fmicb.2016.01207
Wang, Q. Q., Zhu, D. Q., Wang, C., Niu, T. H., Mao, Z. Y., and Xie, J. M. (2020). Effects of chemical fertilizer reduction and bio-organic fertilizer on photosynthesis, yield and quality of garlic seedlings. Gansu Nongye Daxue Xuebao 55, 69–75. doi: 10.13432/j.cnki.jgsau.2020.04.010
Xiao, X., Li, J., Lyu, J., Feng, Z., Zhang, G., Yang, H., et al. (2022). Chemical fertilizer reduction combined with bio-organic fertilizers increases cauliflower yield via regulation of soil biochemical properties and bacterial communities in Northwest China. Front. Microbiol. 13:922149. doi: 10.3389/fmicb.2022.922149
Yang, C. H., and Crowley, D. E. (2000). Rhizosphere microbial community structure in relation to root location and plant iron nutritional status. Appl. Environ. Microb. 66, 345–351. doi: 10.1128/AEM.66.1.345-351.2000
Yang, J. Y., Xu, Y. M., and Sun, J. S. (2022). Effects of long-term fertilization on nutrients in Irrigated Desert soil of eastern Xinjiang. Xinjiang Agric. Sci. 59, 433–440. doi: 10.6048/j.issn.1001-4330.2022.02.021
Yang, W. N., Yu, L., Luo, D. H., Xiong, Z. Y., Wang, Y. Y., Xu, M., et al. (2022). Effect of combined application of biochar with chemical fertilizer and organic fertilizer on soil phosphatase activity and microbial community. Environ. Sci. 43, 540–549. doi: 10.13227/j.hjkx.202105279
Yao, H., and Wu, F. (2010). Soil microbial community structure in cucumber rhizosphere of different resistance cultivars to fusarium wilt. FEMS Microbiol. Ecol. 72, 456–463. doi: 10.1111/j.1574-6941.2010.00859.x
Ye, Z., Li, J., Wang, J., Zhang, C., Liu, G., and Dong, Q. G. (2021). Diversity and co-occurrence network modularization of bacterial communities determine soil fertility and crop yields in arid fertigation agroecosystems. Biol. Fert. Soils 57, 809–824. doi: 10.1007/s00374-021-01571-3
Ye, L., Zhao, X., Bao, E., Li, J., Zou, Z., and Cao, K. (2020). Bio-organic fertilizer with reduced rates of chemical fertilization improves soil fertility and enhances tomato yield and quality. Sci. Rep. 10:177. doi: 10.1038/s41598-019-56954-2
Zhang, T., Kong, Y., Xiu, W. M., Li, G., Zhao, J. N., Yang, D. L., et al. (2019). Effects of fertilization treatments on soil microbial community characteristics under the wheat-maize rotation system in fluvo-aquic soil region in North China. Ecol. Environ. Sci. 28, 1159–1167. doi: 10.16258/j.cnki.1674-5906.2019.06.011
Zhang, T., Li, X. Y., Yang, H. M., Chu, M., and Shi, Y. W. (2018). Isolation, screening and identification of antagonistic bacteria against Verticillium dahliae Kleb. Xinjiang 45, 2418–2428. doi: 10.13344/j.microbiol.china.180003
Zhang, X., Zhang, R., Gao, J., Wang, X., Fan, F., Ma, X., et al. (2017). Thirty-one years of rice-rice-green manure rotations shape the rhizosphere microbial community and enrich beneficial bacteria. Soil Biol. Biochem. 104, 208–217. doi: 10.1016/j.soilbio.2016.10.023
Zhang, M., Zhang, X., Zhang, L., Zeng, L., Liu, Y., Wang, X., et al. (2021). The stronger impact of inorganic nitrogen fertilization on soil bacterial community than organic fertilization in short-term condition. Geoderma 382:114752. doi: 10.1016/j.geoderma.2020.114752
Zhao, W., and Huang, L. M. (2022). Stoichiometric characteristics and influencing factors of soil nutrients under different land use types in an alpine mountain region. Acta Ecol. Sin. 42, 4415–4427. doi: 10.5846/stxb202105061179
Zhao, J., Li, Y., Ran, W., Zhang, R. F., Shen, B., and Shen, Q. R. (2016). Effects of organic manure partial substitution for chemical fertilizer on crop yield and soil microbiome in a rice-wheat cropping system. Nanjing Nongye Daxue Xuebao 39, 594–602. doi: 10.7685/jnau.201603016
Zhao, Y., Mao, X., Zhang, M., Yang, W., Di, H. J., Ma, L., et al. (2021). The application of Bacillus Megaterium alters soil microbial community composition, bioavailability of soil phosphorus and potassium, and cucumber growth in the plastic shed system of North China. Agric. Ecosyst. Environ. 307:107236. doi: 10.1016/j.agee.2020.107236
Zheng, J. C., Yan, M. M., Zhang, J. S., Shi, H. L., and Wang, J. (2015). Effects of nitrogen fertilizer top dressed earlier under the shade on the yield and the fiber quality of cotton. Xinjiang Agric. Sci. 52, 1982–1988. doi: 10.6048/j.issn.1001-4330.2015.11.004
Zhong, W., Gu, T., Wang, W., Zhang, B., Lin, X., Huang, Q., et al. (2010). The effects of mineral fertilizer and organic manure on soil microbial community and diversity. Plant Soil 326, 511–522. doi: 10.1007/s11104-009-9988-y
Zhou, J., Gao, Z., Ma, H., and Chen, W. (2011). “Prevention and study on agricultural non-point source pollution in China” in Advances in electrical engineering and electrical machines (Berlin, Heidelberg: Springer), 565–569.
Keywords: cotton, Illumina MiSeq, organic fertilizer, reduced chemical fertilizer, soil microbial community, yield
Citation: Shi Y, Niu X, Chen B, Pu S, Ma H, Li P, Feng G and Ma X (2023) Chemical fertilizer reduction combined with organic fertilizer affects the soil microbial community and diversity and yield of cotton. Front. Microbiol. 14:1295722. doi: 10.3389/fmicb.2023.1295722
Edited by:
Jian-Wei Guo, Kunming University, ChinaReviewed by:
Nagaraju Yalavarthi, National Bureau of Agriculturally Important Microorganisms (ICAR), IndiaWei Jiang, Western Sydney University, Australia
Copyright © 2023 Shi, Niu, Chen, Pu, Ma, Li, Feng and Ma. 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: XingWang Ma, maxw2001@163.com
†These authors have contributed equally to this work