- 1College of Horticulture, Gansu Agricultural University, Lanzhou, China
- 2Gansu Provincial Key Laboratory of Arid Land Crop Science, Gansu Agricultural University, Lanzhou, China
Reducing chemical fertilizers in combination with bio-organic fertilizers can limit the use of chemical fertilizers while maintaining soil fertility. However, the effects of combined fertilization on soil chemical properties, microbial community structure, and crop yield and quality are unknown. Using high-throughput sequencing, we conducted field experiments using lettuce plants subjected to five fertilization treatments: chemical fertilizer with conventional fertilization rate (CK), chemical fertilizer reduction by 30% + 6,000 kg ha–1 bio-organic fertilizer (T1), chemical fertilizer reduction by 30% + 9,000 kg ha–1 bio-organic fertilizer (T2), chemical fertilizer reduction by 40% + 6,000 kg ha–1 bio-organic fertilizer (T3), and chemical fertilizer reduction by 40% + 9,000 kg ha–1 bio-organic fertilizer (T4). Compared with CK, the T1–T4 had significantly higher soil pH and soil organic matter (SOM) and showed increased richness and diversity of the bacterial community, and decreased richness and diversity of the fungal community. Principal coordinate analysis evidenced that the bacterial and fungal communities of CK and T1–T4 were distinctly separated. The Kruskal-Wallis H-test demonstrated that the fungal community was more sensitive than the bacterial community to chemical fertilizer reduction combined with bio-organic fertilizer. Among the soil chemical parameters measured, only TN (total nitrogen) was significantly correlated with bacterial and fungal community composition. The T1 and T2 increased lettuce yield. Moreover, T1–T4 characterized reduced nitrate content and increased levels of soluble sugars and vitamin C in lettuce. Overall, the combined application of reduced chemical fertilizer and bio-organic fertilizer effectively improved soil fertility, microbial community structure, and lettuce yield and quality. These findings have valuable implications for vegetable safety and long-term environmental sustainability.
Introduction
Fertilization is a globally used management strategy for improving soil quality, and the application of chemical fertilizers is an important determinant of high agricultural production (Yang et al., 2019a). In recent decades, the growth of global crop yields has largely depended on heavy investments in chemical fertilizers (Geng et al., 2019), and farmers manage their farmlands with high fertilization rates to maintain soil productivity. To date, large amounts of chemical fertilizers have been applied to farmlands to obtain higher crop yields (Gruzdeva et al., 2007). However, the amount of chemical fertilizers used is often far greater than that required by crops (Zhu et al., 2019). Excess nutrients in the soil may cause potential damage to the soil, including soil acidification, a decline in soil organic matter (SOM), and a sharp decrease in soil biodiversity (Lv et al., 2020; Wu et al., 2020). Furthermore, it can cause a wide range of environmental problems in marine, freshwater, and terrestrial ecosystems (Tilman et al., 2011). The SOM plays a major role in regulating atmospheric CO2 concentrations and maintaining soil fertility and productivity (Lin et al., 2019). Moreover, organic fertilization is the most effective technique for increasing SOM content and improving soil fertility (Liang et al., 2014; Wang et al., 2017b). Numerous studies have shown that organic fertilization reduces soil’s bulk density and increases hydraulic conductivity and soil organic carbon (SOC) content, which forms the basis of soil fertility (Miller et al., 2008; Bandyopadhyay et al., 2010; Victoria et al., 2012). The application of organic fertilizers also increases soil pH and helps alleviate soil acidification (Afreh et al., 2018).
The development of techniques to improve crop yields, the quality of agricultural products, and the sustainability of agricultural production have attracted increasing attention. Belay et al. (2002) reported that long-term inorganic fertilization increased the levels of total organic carbon in soil and significantly increased grain yield in maize. However, the excessive use of chemical fertilizers causes the accumulation of nitrate in vegetable products, leading to a decline in food safety and quality. The negative effects of excessive use of chemical fertilizers can be minimized by reducing the use of chemical fertilizers and combining them with organic fertilizers, as products from organic agricultural systems have generally better nutritional properties (Luthria et al., 2010; Vallverdú-Queralt et al., 2012; Oliveira et al., 2013). Lettuce (Lactuca sativa var. angustana) is a popular edible vegetable with stems and leaves. China is the largest producer of lettuce, accounting for approximately half of the world’s total lettuce production and cultivation, mainly for stem lettuce products (Nie et al., 2017). The growth cycle of lettuce is relatively short (approximately 90 days). Therefore, although organic fertilizers are more environmentally friendly, nutrient release from organic fertilizers is too slow to support the production of lettuce in a short time (Xiao et al., 2017). The nutrient content of chemical fertilizers is higher than that of organic fertilizers, but it can easily cause soil degradation and environmental pollution. Several studies have found substantive evidence to show that combining chemical and organic fertilizers can improve the quality of a variety of fruits (Wu et al., 2021a). Caris-Veyrat et al. (2004) reported that the levels of carotenoids, polyphenols, and vitamin C were higher in organically cultivated tomatoes than in traditional agricultural tomatoes. In addition, the use of inorganic-organic compound fertilizers not only reduces the use of chemical fertilizers, but also improves the long-term efficiency and sustainability of agroecosystems (Geng et al., 2019).
Soil microorganisms play a crucial role in energy flow and nutrient cycling in soil, and participate in the decomposition of organic matter, degradation of xenobiotics, soil carbon sequestration, and prevention of crop diseases (Morris and Blackwood, 2015; Tao et al., 2015). Bacteria and fungi are indicator species of the soil environment (Hermans et al., 2017; Wolińska et al., 2021), and are conducive to the improvement of the soil structure and plant growth. Fungi can grow symbiotically with crops to form mycorrhizae, which are of major significance for promoting plant growth and maintaining agroecosystem stability (Brito et al., 2012). Bacteria and fungi also have a cooperative relationship in soil environments that lack nutrients; that is, fungal hyphae are beneficial for improving and balancing nutrients in the soil, thereby promoting the growth and metabolism of bacteria (Strickland and Rousk, 2010). Soil microbial biomass and diversity are potentially valuable indicators of soil quality and are sensitive to changes in soil nutrients, pH, and organic matter content (Zhong et al., 2010). Microorganisms can quickly respond to changes in environmental factors by modifying microbial diversity and community composition (Krashevska et al., 2015; Wolińska et al., 2017). Soil microbial communities are greatly affected by anthropogenic activities such as agricultural intensification and fertilization. Moreover, nutrient sources and application rates determine nutrient availability in soil; therefore, the application of different nutrient sources may cause significant shifts in the predominant microbial taxa (Hartmann et al., 2015; Sun et al., 2015). Previous studies have reported changes in soil microbial communities after fertilization. Based on 27 years of experiments, Lin et al. (2019) reported that organic fertilizers had a strong impact on soil microbial community structure and led to an increase in the abundance of Bacillales, Gaiellales, and Pezizales. Regardless of whether chemical fertilizers are used, organic fertilization has a positive effect on bacterial and fungal diversity (Kamaa et al., 2011). The combined application of organic and inorganic fertilizers can accelerate microbial growth, change the soil microbial community structure, and improve enzyme activity (Lazcano et al., 2013). Therefore, the use of organic fertilizers to maintain the stability of underground ecosystem function has attracted increased attention.
Previous studies have primarily focused on long-term chemical or organic fertilization alone, and few studies have investigated the use of bio-organic fertilizers, combined with reduced chemical fertilizers, especially in lettuce. As a result, the variation in yield, quality, and soil bacterial and fungal community diversity in organic-inorganic compound fertilization systems remains unclear. As discussed above, excessive chemical fertilization has negative effects on the environment, whereas organic fertilization has beneficial effects on crop’s growth and quality. Exploring how the soil’s bacterial and fungal communities respond to the agricultural fertilization system would be conducive to the development of sustainable agriculture. Our objectives were as follows: (i) to evaluate how different degrees of reduction of chemical fertilizers combined with bio-organic fertilizers affected crop yield and quality, soil nutrients, soil microbial (bacterial and fungal) diversity and composition, and soil microbial community structure; (ii) to assess the relationship between soil chemical properties and microbial communities. The results of this study provide scientific support for the establishment of an environmentally friendly fertilization technology and effectively promote the sustainable development of agriculture.
Materials and Methods
Field Site and Experiment Design
The study was conducted in Tianzhu, Gansu, China (37° 09′ N, 102° 99′ E), in May 2019 and 2020. The altitude of this region is approximately 2,630 m, the annual average temperature is −2°C, and the average annual precipitation is 400–450 mm. At the study location, the soil type was loam, the terrain was flat, and the soil fertility was uniformly moderate. Lettuce seeds (cv. “Yong’an lettuce—-Zilong lettuce”) were sown at Fuzhou Kexiang Seed Industry Co., Ltd., and thinned thereafter. The conventional fertilizers used in this study were urea (N ≥ 46%), calcium superphosphate (P2O5 ≥ 16%), and potassium sulfate (K2O ≥ 52%). The bio-organic fertilizer used in the experiment (organic matter content ≥ 40%, Bacillus subtilis, and Pseudomonas stutzeri with effective viable count ≥ 20 million CFU/g) was purchased from Gansu Lvneng Agricultural Technology Co., Ltd. The calcium superphosphate and bio-organic fertilizers are both applied as basal fertilizers. In this study, 30% of the urea and potassium sulfate were used as basal fertilizers. Moreover, 70% of the urea and potassium sulfate were top-dressed twice; of this, 30% was applied in the seedling stage and 40% in the rhizome expansion stage.
Five treatments were used in this experiment: (i) chemical fertilizer with conventional fertilization, as used by local farmers (CK; NPK dosage: 965.2 kg ha–1: 314.3 kg N ha–1, 478.5 kg P ha–1, and 172.4 kg K ha–1); (ii) chemical fertilizer reduction by 30% + 6,000 kg ha–1 bio-organic fertilizer (T1; NPK dosage: 675.6 kg ha–1: 293.9 kg N ha–1, 117.2 kg P ha–1, and 264.5 kg K ha–1); (iii) chemical fertilizer reduction by 30% + 9,000 kg ha–1 bio-organic fertilizer (T2; NPK dosage: 675.6 kg ha–1: 293.9 kg N ha–1, 117.2 kg P ha–1, and 264.5 kg K ha–1); (iv) chemical fertilizer reduction by 40% + 6,000 kg ha–1 bio-organic fertilizer (T3; NPK dosage: 579.1 kg ha–1: 251.9 kg N ha–1, 100.5 kg P ha–1, and 226.7 kg K ha–1); and (v) chemical fertilizer reduction by 40% + 9,000 kg ha–1 bio-organic fertilizer (T4; 579.1 kg ha–1: 251.9 kg N ha–1, 100.5 kg P ha–1, and 226.7 kg K ha–1). Each treatment included three replicates, and the area of each experimental plot was 80 m2. We used the ridge, double row, and semi-film cultivation modes for the experiments. The cultivation density was 57,750 plants ha–1, the ridge width was 50 cm, the furrow width was 40 cm, and the plant spacing was 30 cm. We ensured that all field management practices (other than fertilization) were consistent across treatments.
Soil Sampling and Analysis of Soil Chemical Properties
After lettuce harvesting in August 2020, soil samples were collected from the plow layer (0–20 cm) of each replicate plot for each treatment using a stainless-steel auger. Each soil sample was a composite of five randomly collected soil cores that were pooled in a sterile plastic bag and transported to the laboratory on ice. The soil samples were sieved using a 2-mm mesh sieve and divided into two subsamples before thorough homogenization. One part was air-dried to determine the physical and chemical properties of the soil, and the other part was stored at −80°C for DNA extraction.
The soil pH and electrical conductivity (EC) were measured after shaking the soil-water suspension (1:5 wt/vol) for 30 min. The pH of the filtrate was determined using a glass electrode (PHS-3E, Shanghai Jingke, China), and the EC was measured by inserting a conductivity meter (DSJ-308A, Shanghai Jingke) into the filtrate. The SOM content was measured using a titration method after oxidation with K2Cr2O7. The soil’s total nitrogen (TN), total phosphorus (TP), and total potassium (TK) were determined using the H2SO4-H2O2 wet digestion method. The TN was measured using the Kjeldahl method using a fully automatic Kjeldahl K1100F apparatus (Jinan Hanon Instruments Company, Jinan, China). The TP was measured using the Mo-Sb colorimetric method and analyzed using a UV-1780 spectrophotometer [Shimadzu Instruments (Suzhou) Co., Ltd., Suzhou, China]. The TK was analyzed using a ZEEnit 700P atomic absorption spectrometer (Analytik Jena, Germany). All samples were tested in triplicates.
Soil DNA Extraction and PCR Amplification
The genomic DNA of the soil microbial community was extracted from 0.5 g soil samples using the E.Z.N.A.® soil DNA Kit (Omega Bio-tek, GA, United States) according to the manufacturer’s instructions. The DNA extract was analyzed on a 1% agarose gel, and DNA concentration and purity were determined using a NanoDrop ND-2000 UV-VIS spectrophotometer (Thermo Scientific, Wilmington, DE, United States). The hypervariable V3–V4 region of the bacterial 16S rRNA gene was amplified with the primer pairs 338F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTA AT-3′) (Chen et al., 2016) using an ABI GeneAmp 9700 PCR thermocycler (Applied Biosystems, CA, United States). The internal transcribed spacer (ITS) region of fungi was amplified using the primers ITS1F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and ITS2R (5′-GCTGCGTTCTTCATCGATGC-3′) (Blaalid et al., 2013) on the same thermal cycler. The PCR was conducted using 20 μl reaction mixtures containing the template DNA (10 ng), 4 μL of 5 × TransStart FastPfu buffer (TransGen Biotech, Beijing, China), 2 μL of 2.5 mM dNTPs, 0.8 μL each of the forward and reverse primers (5 μm), 0.4 μl of TransStart FastPfu DNA Polymerase, and ddH2O. The PCR conditions were as follows: initial denaturation at 95°C for 3 min, 27 cycles of denaturation at 95°C for 30 s, annealing at 55°C for 30 s, and extension at 72°C for 45 s, a single extension at 72°C for 10 min, and storage at 4°C. The PCRs were performed in triplicate. The PCR product was extracted from a 2% agarose gel, purified using the AxyPrep DNA gel extraction kit (Axygen Biosciences, CA, United States) according to the manufacturer’s instructions, and quantified using a Quantus Fluorometer (Promega, WA, United States).
Illumina Miseq Sequencing and Sequence Processing
Purified amplicons were pooled in equimolar amounts and paired-end sequenced (2 × 300 bp) on an Illumina MiSeq PE300 platform (Illumina, CA, United States) by Shanghai Majorbio Bio-pharm Technology Co., Ltd. (Shanghai, China). The raw sequences were analyzed using the Fastp software (version 0.19.6)1 (Chen et al., 2018) for quality control. The paired-end reads were merged using the FLASH software (version 1.2.11)2 (Magoč and Salzberg, 2011). The reads were dereplicated, sorted, and clustered into operational taxonomic units (OTUs) at the default 97% similarity using UPARSE (version 7.0.1090)3 (Edgar, 2013), and dechimerized against the UCHIME reference dataset (Edgar et al., 2011). Taxonomic labels were assigned to the OTUs using the Silva 16S bacterial database and the UNITE fungal database (version 7.2) (Nilsson et al., 2019) using the RDP classifier (version 2.11)4 with a minimum confidence threshold of 0.7.
Determination of the Yield and Quality of Lettuce
In August 2019 and 2020, the lettuce crops were harvested after they had reached commercial standards. The yields of all plots were determined. The nitrate content was determined using the salicylic acid-sulfuric acid method (Cataldo et al., 1975). The 2,6-dichloroindophenol staining method was used to determine the vitamin C content (Arya et al., 2000), and the Coomassie Brilliant Blue method was used to estimate the levels of soluble protein (Sedmak and Grossberg, 1977). The soluble sugar content was measured using the anthrone-sulfuric acid assay (Wang et al., 2020).
Statistical Analysis
Microsoft Excel 2013 was used to calculate the mean and standard error (SE) of the crop yield and quality of lettuce and the chemical properties of soil in different treatment groups. One-way analysis of variance (ANOVA) was applied to evaluate the effects of different treatments on the lettuce yield and quality and soil chemical properties. Duncan’s multiple range test (significance level: p < 0.05) was used to compare the differences in mean values among different fertilization treatments. The SPSS software (version 23.0, IBM Corp., NY, United States) was used for one-way ANOVA.
Alpha diversity indices (coverage, ACE, Chao1, Shannon index, and Simpson index) at the OTU level were calculated using the Mothur software (version 1.30.2). For the analysis of beta diversity, we used unweighted UniFrac distances to generate principal coordinate analysis (PCoA) maps to assess the differences between the bacterial community members of soil samples with different fertilization treatments (Han et al., 2021). In addition, we used PCoA based on Bray–Curtis distances to evaluate differences between the fungal community members of soil samples with different fertilization treatments. Redundancy analysis (RDA) was applied to analyze the correlation between microbial community composition and soil chemical properties (Yang et al., 2019c). The Mantel test based on Bray–Curtis distances was used to determine the relationship between microbial community composition and soil properties at the OTU level. The PCoA, heatmap, and RDA were computed using the vegan package in R (Oksanen et al., 2015).
Results
Chemical Properties of Soils With Different Fertilization Treatments
Table 1 lists the chemical properties of soils with different fertilization treatments. Compared with the CK treatment, the four treatments (T1, T2, T3, and T4) with reduced chemical fertilizer combined with bio-organic fertilizer significantly increased the soil pH (Table 1). The SOM values of the T1, T2, T3, and T4 treatment groups were significantly higher compared to those in the CK treatment group by 61.94, 59.18, 65.60, and 71.12%, respectively. The EC in the T1, T2, T3, and T4 treatment groups was significantly lower than that in the CK group; however, the EC of T1 and T2 was significantly higher than that of T3 and T4. The T2 treatment group had the highest TN and TK levels, which were significantly higher (by 35.71 and 16.69%, respectively) than that in the CK group. There was no significant difference in TP among the T2, T3, T4, and CK treatments.
Alpha Diversity of Bacterial and Fungal Communities in Soils With Different Fertilization Treatments
Coverage refers to the coverage of each sample library. The higher the value, the higher the probability that the sequence in the sample will be detected. This index reflects whether the sequencing results represent the actual microorganisms in the sample. The coverage indices for both bacteria and fungi exceeded 0.96 (Table 2), indicating that the sequencing results reflected the actual microorganisms in the soil samples. The higher the values of community richness indices (Chao1 and ACE), the higher the community richness. In terms of the bacterial community, the Chao1 and ACE indices in the T1–T4 treatment groups were higher (2,950.93–3,058.16 and 2,942.68–3,060.67, respectively) than those in the CK treatment group (2,859.13 and 2,864.02, respectively). In particular, the Chao1 and ACE indices of the T3 treatment group were significantly higher than those of the CK group, although there was no significant difference between the other three treatment groups and the CK treatment group (Table 3, p < 0.05). For the fungal community, the Chao1 and ACE indices of the T1–T4 treatment groups were lower (582.66–690.62 and 585.28–695.51, respectively) than those in the CK treatment group (723.01 and 708.37, respectively). In particular, the Chao1 and ACE indices of the T2 treatment group were significantly lower than those of the CK treatment group, and there was no significant difference between the other three treatment groups and the CK treatment group (Table 3, p < 0.05). We also used the Shannon and Simpson diversity indices to assess the diversity of each sample. The higher the Shannon index, the higher the microbial diversity of the sample, and the higher the Simpson index, the lower the microbial diversity of the sample. The Shannon indices of bacteria were higher in the T2 and T3 treatment groups (6.393 and 6.408, respectively) than in the CK treatment group (6.388), but the difference was not significant (Table 2). The Simpson index of bacteria was the lowest in the T2 treatment group. The Shannon indices of fungi were significantly lower in the T1–T4 treatment groups than in the CK treatment group (p < 0.05), whereas the Simpson index showed the opposite trend.
Table 3. Mantel test results for correlations between soil chemical properties and bacterial and fungal community composition at the OTU level.
Beta Diversity of Bacterial and Fungal Communities in Soils With Different Fertilization Treatments
The PCoA analysis showed that the community composition of soil bacteria and fungi was altered by the reduction of chemical fertilizer combined with bio-organic fertilizer (Figure 1). For the bacterial community structure, the first and second principal coordinates explained 9.97 and 8.65% of the differences among the five treatments, respectively (Figure 1A). Moreover, the bacterial communities of the T1–T4 treatment groups were considerably separated from those of the CK treatment group (Figure 1A). The bacterial community structures of the T1 and T3 treatment groups appeared to be similar. For the fungal community structure, the first and second principal coordinates explained 36 and 17.21% of the variance, respectively (Figure 1B). The fungal communities of the T1–T4 treatment groups were clearly distinguished from those of the CK treatment group along the first principal component axis. Moreover, the fungal community structure was similar between the T3 and T4 treatment groups.
Figure 1. Principal coordinate analysis (PCoA) plots of bacterial (A) and fungal (B) community composition at the OTU level. CK, chemical fertilizer with conventional fertilization rate; T1, chemical fertilizer reduction by 30% + 6,000 kg ha–1 bio-organic fertilizer; T2, chemical fertilizer reduction by 30% + 9,000 kg ha–1 bio-organic fertilizer; T3, chemical fertilizer reduction by 40% + 6,000 kg ha–1 bio-organic fertilizer; T4, chemical fertilizer reduction by 40% + 9,000 kg ha–1 bio-organic fertilizer.
Composition and Relative Abundance of Bacterial and Fungal Communities in Soils With Different Fertilization Treatments
The phylum composition of the bacteria in the five treatment groups is shown in Figure 2A. The top five bacterial phyla were Actinobacteria (31.25–32.68%), Proteobacteria (17.22–19.88%), Chloroflexi (16.86–19.54%), Acidobacteria (12.65–16.01%), and Gemmatimonadota (4.31–4.61%). These five bacterial phyla accounted for more than 85% of the total sequence reads. Compared to the CK group, the relative abundance of Actinobacteria was higher in the T1, T2, T3, and T4 treatment groups. The relative abundance of Gemmatimonadota in the CK treatment was lower than those in the T1–T4 treatment groups. There was no significant difference in the relative abundance of other phyla among all treatments, except for Bacteroidetes (p < 0.05). The relative abundance of Bacteroidetes was significantly lower in the T1, T2, and T4 groups than in the CK treatment group (Figure 2C, p < 0.05).
Figure 2. Relative abundances of bacterial (A) and fungal (B) taxa at the phylum level. The relative abundances of bacteria (C) and fungi (D) were tested by the Kruskal–Wallis H-test. “Other” includes phyla with <1% average relative abundance. Values for individual treatments are the means of three replicate soil samples. The false discovery rate was used to control for multiple tests, and the Scheffe test was used to calculate the 95% confidence intervals. CK, T1, T2, T3, and T4 are as defined in the Figure 1 legend.
The dominant fungal phyla in all treatment groups were Ascomycota, Mortierellomycota, and Basidiomycota, which accounted for 65.10–72.87%, 18.39–27.28%, and 3.37–4.70% of all sequence reads, respectively (Figure 2B). Overall, these three fungal phyla accounted for more than 90% of the high-quality sequences. Compared with the CK treatment group, the T2 treatment decreased the relative abundance of Ascomycota and increased that of Basidiomycota. Furthermore, the Kruskal–Wallis H-test showed that the relative abundances of Mortierellomycota, Olpidiomycota, Glomeromycota, Blastocladiomycota, and Zoopagomycota were significantly different among the five treatment groups (Figure 2D, p < 0.05).
Relationship Between Microbial Community Structure and Soil Chemical Properties
The RDA analysis in Figure 3 shows the relationship between soil chemical properties and microbial community composition. The first two axes of the RDA analysis explained 40.46% of the total variance in the soil bacterial community (axis 1: 31.12%; axis 2: 9.34%) (Figure 3A). The TN showed the highest correlation with the soil bacterial community structure, followed by TP, SOM, TK, pH, and EC. The first two axes of the RDA analysis explained 42.50% of the total variance in the soil fungal community (axis 1: 38.95%; axis 2: 3.55%) (Figure 3B). Soil chemical properties influenced the soil fungal community composition in the following order: EC > TP > TN > SOM > pH > TK. The results of a Mantel test (Table 3) revealed that TN, SOM, pH, and TK were significantly (p < 0.05) or extremely significantly (p < 0.01) positively correlated with bacterial community composition, and that TN was significantly positively correlated with fungal community composition (p < 0.01).
Figure 3. Redundancy analysis (RDA) of bacterial (A) and fungal (B) communities with soil chemical properties. EC, electrical conductivity; SOM, soil organic matter; TN, total nitrogen; TP, total phosphorus; TK, total potassium. CK, T1, T2, T3, and T4 are as defined in the Figure 1 legend.
Lettuce Yield and Quality
Among the five fertilization treatments, the T2 treatment group had the highest yield (Figures 4A,B). The 2-year yields of the T1 and T2 treatment groups were significantly higher than that of the CK group (p < 0.05). Compared with the 2-year yield of the CK treatment group, the yield of the T2 treatment group increased by 19.43 and 16.63%, respectively. The nitrate content of lettuce was significantly lower in the T1–T4 treatment groups (by 25.25, 25.46, 32.88, and 32.44%, respectively) than in the CK group (Figure 4B, p < 0.05). However, there was no significant difference in vitamin C content among the five fertilization treatments (Figure 4C, p < 0.05). The soluble protein content of the T2 and T4 treatment groups were significantly higher (by 22.49 and 24.30%, respectively) compared to that in the CK group (Figure 4D, p < 0.05). The soluble sugar content of lettuce was higher in the T1–T4 treatment groups than in the CK group (Figure 4E). Compared to the soluble sugar content of the CK group, that of the T2 treatment increased significantly by 32.19%, and those of T1, T2, and T3 groups increased by 19.62, 22.52, and 12.29%, respectively (Figure 4F, p < 0.05).
Figure 4. Yield and quality parameters of lettuce grown in soils with different fertilization treatments. (A) Yield in 2019. (B) Yield in 2020. (C) Nitrate content of lettuce. (D) Vitamin C content of lettuce. (E) Soluble protein content of lettuce. (F) Soluble sugar content of lettuce. Values presented are the means ± SE (n ≥ 3). Vertical bars indicate the SE of the means. Different lowercase letters indicate significant differences among treatments at the same timepoint (Duncan’ s test, p < 0.05). CK, T1, T2, T3, and T4 are as defined in the Figure 1 legend.
Discussion
Effects of the Reduction of Chemical Fertilizer Combined With Bio-Organic Fertilizer Treatment on Soil Chemical Properties and Microbial Diversity
The results of the present study demonstrated that the soil chemical properties changed significantly after the application of chemical fertilizer combined with bio-organic fertilizer. Specifically, the pH and SOM of the four treatment groups (reduced chemical fertilizer combined with bio-organic fertilizer) were significantly improved compared with the values in the control (CK) group (Table 1, p < 0.05). This was consistent with the results of previous studies showing that the application of organic fertilizer improved the soil chemical properties in rice (Wang et al., 2014). The significant increase in soil pH may be ascribed to the mitigation of soil acidification through the decomposition of organic matter (Xiaojing et al., 2018). Organic fertilizers are an important source of SOM, and SOM with a higher C content is important for maintaining soil fertility and improving soil properties (Yang et al., 2019b; Zhu et al., 2019). Pinggera et al. (2015) found that adding an appropriate amount of nitrogen facilitates the microbial decomposition of organic matter and adjusts the C/N ratio to levels more suitable for the growth of soil microorganisms. This accelerates decomposition and the release of organic nutrients by microorganisms and increases the content of SOM, which may explain the increase in SOM in the T1–T4 treatment groups in this study. In addition, the TN content was increased by a 30% reduction in chemical fertilizers combined with 6,000 or 9,000 kg ha–1 bio-organic fertilizer (Table 1), which is consistent with the results of Yang et al. (2019b) and Wu et al. (2021b).
Previous studies on arable soil have shown that conventional fertilization with chemical fertilizers reduced bacterial richness and increased fungal richness (as measured using the Chao1 and ACE indices) (Shen et al., 2010; Zhong et al., 2010). In this study, the group with conventional chemical fertilization had the lowest Chao1 and ACE indices for bacteria and the highest indices for fungi (Table 2); thus, these results were consistent with those of the aforementioned studies. Ai et al. (2012) also reported that the application of mineral nitrogen fertilizers generally promotes the growth of fungi. In contrast, we found that the Chao1 and ACE indices were significantly lower in the T2 treatment group (30% fertilizer reduction + 9,000 kg ha–1 bio-organic fertilizer) than in the T3 and T4 treatment groups (40% fertilizer reduction + 6,000 or 9,000 kg ha–1 bio-organic fertilizer, p < 0.05). This may be because among the T1–T4 treatment groups, the T2 group had the highest pH, and studies have shown that within a pH range of 4–8.3, lower pH promotes the growth of fungi (Rousk et al., 2009). There was no significant change in the Shannon diversity index of bacteria in the treatment groups (p < 0.05), indicating that the bacteria showed adaptability to the soil under different fertilization systems. Luo et al. (2015) discovered that compared to inorganic fertilizers, organic fertilizers are more conducive to increasing fungal diversity. However, in this study, the Shannon diversity index of fungi was reduced by a treatment with reduced chemical fertilizers combined with bio-organic fertilizer, and the Simpson index showed the opposite trend. Studies using culture-dependent methods have shown that long-term organic fertilization—especially with manure—can reduce fungal diversity (He et al., 2008). These contrasting results are likely due to differences in experimental methods and soil types. Moreover, the PCoA analysis revealed that the bacterial and fungal community compositions were clearly separated between the control and treatment groups. Similar results have also been reported by Yu et al. (2013) and Ding et al. (2017).
Effects of Reduced Chemical Fertilizer Combined With Bio-Organic Fertilizer on Soil Microbial Community Composition
Inorganic or organic fertilization leads to the enrichment of specific bacteria or fungi that can effectively utilize these nutrients, thereby altering the composition of the microbial community. Elucidating soil microbial taxa at the phylum level can reveal the ecological consistency of a microbial community (Guo et al., 2018). In this study, high-throughput sequencing showed that the dominant bacterial phyla were Actinobacteria, Proteobacteria, Chloroflexi, Acidobacteria, and Gemmatimonadota (Figure 2A), which is consistent with previous studies on agricultural soils (Chen et al., 2016; Guo et al., 2018; Wu et al., 2021b). Actinobacteria are eutrophic bacteria that are considered an indicator taxon and play a key role in the soil carbon and nitrogen cycle (Zhu et al., 2019). The present study showed that the relative abundance of Actinobacteria was increased in the four treatment groups (T1–4). Wang et al. (2017a) reported that the relative abundance of Bacteroidetes increases after 10 years of inorganic fertilization. Similarly, in this study, the relative abundance of Bacteroidetes was significantly higher in the control group (conventional chemical fertilization) than in the combined fertilization groups (reduced chemical fertilizer and bio-organic fertilizer). Consistent with previous studies, Ascomycota was the most abundant fungal phylum in all fertilization treatments (Chen et al., 2016; Bei et al., 2018). Ascomycota and Basidiomycota are typical saprophytes that are environmentally friendly, decompose organic matter (cellulose, lignin, and pectin), and play an important role in degrading litter with high lignin content in the soil (Yelle et al., 2008; Guo et al., 2018; Liu et al., 2019). Compared with the CK treatment group, the T2 group (30% reduction of chemical fertilizer combined with 6,000 kg ha–1 bio-organic fertilizer) had a lower relative abundance of Ascomycota and increased relative abundance of Basidiomycota (Figure 2B). This may be due to a competitive relationship between the Ascomycota and Basidiomycota microbes that results in a negative correlation between the relative abundances of the two phyla (Ye et al., 2020a). Several studies have shown that compared with bacteria, fungi are more sensitive to invasive plants (Xiao et al., 2014), salt stress (Chowdhury et al., 2011), soil moisture (Kaisermann et al., 2015), fertilizer compost (Shen et al., 2019) and other factors. Yang et al. (2020b) demonstrated that compared to soil bacteria, soil fungi are more sensitive to the application of organic fertilizer (pig manure). Furthermore, the Kruskal–Wallis H- test also showed that the only bacterial phylum with significant differences among the five treatments was Bacteroidetes, whereas the other fungal phyla—including Mortierellomycota, Olpidiomycota, Glomeromycota, Blastocladiomycota, and Zoopagomycota—were significantly different between groups (Figure 2, p < 0.05).
Relationship Between the Microbial Community and Soil Environment Factors
TN is significantly associated with ammonia-oxidizing bacterial communities and is one of the main factors promoting changes in the soil microbial community (Qu et al., 2016; Yang et al., 2017). Actinobacteria is a dominant taxon in the bacterial community and is positively correlated with the soil nitrogen pool; therefore, TN has been proven to drive changes in the soil bacterial community (Nemergut et al., 2010). In this study, TN was found to be the main factor driving changes in the bacterial community (Figure 3A). pH is also considered the most significant factor determining the composition of the soil bacterial community (Zhao et al., 2015; Kim et al., 2016), and the results of the Mantel test showed that TN and pH were significantly correlated with soil bacterial community composition (Table 3, p < 0.05). Fungal community composition is affected by soil nutrient content (i.e., TN and TP) (Paungfoo-Lonhienne et al., 2015; He et al., 2016). Accordingly, the main driving factors of the fungal community were EC, TP, and TN in this study (Figure 3B). The results of the Mantel test showed that TN was significantly positively correlated with the fungal community composition (Table 3). Soil nutrients may indirectly affect fungal community composition by regulating changes in the aboveground portions of plants (Haiyan et al., 2016). Previous studies have demonstrated that SOM is the primary driver factor regulating soil fungal community composition (Liu et al., 2015). This is inconsistent with the results of this experiment and suggests that the mineralization of SOM by organic fertilizers is a long-term process. 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). Thus, the short-term effects of organic fertilizers on soil microbial communities appear to be more modest than those of chemical nitrogen fertilizers.
Effects of Reduced Chemical Fertilizer Combined With Bio-Organic Fertilizer on Lettuce Yield and Quality
Excessive chemical fertilization is related to several problems, including high levels of N–P–K in vegetable production and the safety of vegetable products in the market. Bio-organic fertilizers are a potential alternative to mineral fertilizers in crop production that can address these urgent problems (Naeem et al., 2006; Flores-Félix et al., 2013). The application of bio-organic fertilizers to improve plant nutrition can help reduce mineral fertilization and is becoming a popular strategy for sustainable agriculture. Bio-organic fertilizers are a new type of organic fertilizer that contains various beneficial microbial communities. These unique microbial communities can activate the soil, improve its physical and chemical properties, and increase soil microbial diversity (Bandyopadhyay et al., 2010; Xue et al., 2010). However, reducing the use of chemical fertilizers in agriculture without causing loss of productivity is a feasible—but difficult—challenge. Bio-organic fertilizers can replace 23–52% of nitrogen fertilizers without causing a loss of yield (Rose et al., 2014). In this study, we found that a 30% or 40% reduction in chemical fertilizers combined with 6,000 or 9,000 kg ha–1 bio-organic fertilizer did not cause a significant loss of yield, and the yields of the T1 and T2 treatment groups were higher than that of the control groups with conventional chemical fertilization (Figures 4A,B). Thus, our results reveal the potential role of bio-organic fertilizers in reducing the chemical fertilizer load in soil. Similar results have also been reported in cabbages (Sarkar and Rakshit, 2021), cucumbers (Zhao et al., 2018), tomatoes (Ye et al., 2020b), peppers (Wu et al., 2015), and other crops, and the use of bio-organic fertilizers has also been shown to increase crop yield.
Human diets containing vegetables with a high nitrate content can cause serious health hazards. Excessive application of mineral fertilizers—especially N fertilizers—is the main reason for high levels of nitrate accumulation in plant tissues (Wang et al., 2008). In the present study, the nitrate content of the T1–4 treatment groups was significantly lower (by 25.25–32.88%) than that of the CK group. This may be because the improvement in soil microorganisms reduces the nitrate content by increasing the activity of nitrate reductase, thus improving the safety and quality of vegetables (Chatterjee et al., 2012). Qi et al. (2021) found that the levels of soluble sugars and vitamin C increased in “yellow rose” cabbage treated with bio-organic fertilizer. In this study, we found that in lettuce, the levels of soluble sugars, soluble protein, and vitamin C increased to varying degrees in the T1–4 treatment groups (Figure 4). Similar results have also been reported in tomatoes (Brunetti et al., 2019) and cauliflowers (Yang et al., 2020a). Based on these findings, we speculate that the application of reduced chemical fertilizer combined with bio-organic fertilizer can promote the yield and quality of lettuce by improving soil chemical properties and microbial communities (Figure 5).
Figure 5. Potential mechanism of action of reduced chemical fertilizer combined with bio-organic fertilizer in regulating soil microecology to improve the yield and quality of lettuce. The red arrow indicates an increase, and the blue arrow indicates a decrease. The red box indicates bacterial or fungal phyla with significant differences between treatments.
Conclusion
In the present study, we revealed that the reduction of chemical fertilizer combined with bio-organic fertilizer prevented soil acidification and effectively improved the soil’s chemical properties. Treatments with reduced chemical fertilizer combined with bio-organic fertilizer increased the richness and diversity of the bacterial community and decreased those of the fungal community. The PCoA analysis results showed that the bacterial and fungal communities were separated between the treatment (reduced chemical fertilizer combined with bio-organic fertilizer) and CK groups. Different levels of reduced chemical fertilizers combined with bio-organic fertilizers also had an impact on the composition of the bacterial and fungal communities. The Kruskal–Wallis H-test showed that significant differences were more evident in the fungal phyla than in the bacterial phyla among the treatment groups, suggesting that the fungal community was more sensitive to reductions in chemical fertilizer combined with bio-organic fertilizer. The results of the RDA and Mantel test analysis indicated that the short-term impact of organic fertilizers on the soil microbial community was milder than that of chemical N fertilizer. Further, reduced chemical fertilizer combined with bio-organic fertilizer did not significantly affect the lettuce yield or increase it, and improved its safety and nutritional quality. Current agricultural practices present a contradiction between the rapid improvement of soil fertility and the reduction of pollution by chemical fertilizers. Against this background, the combined application of reduced chemical fertilizer and bio-organic fertilizer effectively improved soil fertility and microbial community structure and promoted the yield and quality of lettuce. This study provides a scientific basis for establishing an environmentally friendly fertilization technology and effectively promoting sustainable agriculture.
Data Availability Statement
The datasets (SRP359012) presented in this study can be found in the NCBI Sequence Read Archive (https://identifiers.org/ncbi/insdc.sra:SRP359012).
Author Contributions
NJ, JLy, and JY designed the study. NJ, LJ, SW, and JLi performed the work. NJ was involved with writing the manuscript. NJ, LJ, and FL analyzed the data. ZL, SL, YW, JLy, and JY revised the manuscript. All authors contributed to the article and approved the final version of this manuscript.
Funding
This research was funded by the Education Science and Technology Innovation Project of Gansu Province (GSSYLXM-02), the Special Project of Central Government Guiding Local Science and Technology Development (ZCYD-2021-07), Gansu People’s Livelihood Science and Technology Project (20CX9NA099), Fuxi Young Talents Fund of Gansu Agricultural University (GAUfx-04Y03), and Gansu Top Leading Talent Plan (GSBJLJ-2021-14).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We are grateful for the sequencing platform provided by Majorbio Biomedical Technology Co., Ltd. (Shanghai, China).
Footnotes
- ^ https://github.com/OpenGene/fastp
- ^ https://ccb.jhu.edu/software/FLASH/index.shtml
- ^ http://drive5.com/uparse/
- ^ http://sourceforge.net/projects/rdp-classifier/
References
Afreh, D., Zhang, J., Guan, D., Liu, K., Song, Z., Zheng, C., et al. (2018). Long-term fertilization on nitrogen use efficiency and greenhouse gas emissions in a double maize cropping system in subtropical China. Soil Till. Res. 180, 259–267. doi: 10.1016/j.still.2018.03.016
Ai, C., Liang, G., Sun, J., Wang, X., and Zhou, W. (2012). Responses of extracellular enzyme activities and microbial community in both the rhizosphere and bulk soil to long-term fertilization practices in a fluvo-aquic soil. Geoderma 173, 330–338. doi: 10.1016/j.geoderma.2011.07.020
Arya, S., Mahajan, M., and Jain, P. (2000). Non-spectrophotometric methods for the determination of Vitamin C. Anal. Chim. Acta 417, 1–14. doi: 10.1016/S0003-2670(00)00909-0
Bandyopadhyay, K., Misra, A., Ghosh, P., and Hati, K. (2010). Effect of integrated use of farmyard manure and chemical fertilizers on soil physical properties and productivity of soybean. Soil Till. Res. 110, 115–125. doi: 10.1016/j.still.2010.07.007
Bei, S., Zhang, Y., Li, T., Christie, P., Li, X., and Zhang, J. (2018). Response of the soil microbial community to different fertilizer inputs in a wheat-maize rotation on a calcareous soil. Agr. Ecosyst. Environ. 260, 58–69. doi: 10.1016/j.agee.2018.03.014
Belay, A., Claassens, A., and Wehner, F. (2002). Effect of direct nitrogen and potassium and residual phosphorus fertilizers on soil chemical properties, microbial components and maize yield under long-term crop rotation. Biol. Fert. Soils 35, 420–427. doi: 10.1007/s00374-002-0489-x
Blaalid, R., Kumar, S., Nilsson, R. H., Abarenkov, K., Kirk, P., and Kauserud, H. (2013). ITS 1 versus ITS 2 as DNA metabarcodes for fungi. Mol. Ecol. Resour. 13, 218–224. doi: 10.1111/1755-0998.12065
Brito, I., Goss, M. J., de Carvalho, M., Chatagnier, O., and van Tuinen, D. (2012). Impact of tillage system on arbuscular mycorrhiza fungal communities in the soil under Mediterranean conditions. Soil Till. Res. 121, 63–67. doi: 10.1016/j.still.2012.01.012
Brunetti, G., Traversa, A., De Mastro, F., and Cocozza, C. (2019). Short term effects of synergistic inorganic and organic fertilization on soil properties and yield and quality of plum tomato. Sci. Hortic. 252, 342–347. doi: 10.1016/j.scienta.2019.04.002
Caris-Veyrat, C., Amiot, M.-J., Tyssandier, V., Grasselly, D., Buret, M., Mikolajczak, M., et al. (2004). Influence of organic versus conventional agricultural practice on the antioxidant microconstituent content of tomatoes and derived purees; consequences on antioxidant plasma status in humans. J. Agr. Food Chem. 52, 6503–6509. doi: 10.1021/jf0346861
Cataldo, D., Maroon, M., Schrader, L. E., and Youngs, V. L. (1975). Rapid colorimetric determination of nitrate in plant tissue by nitration of salicylic acid. Commun. Soil Sci. Plan. 6, 71–80. doi: 10.1080/00103627509366547
Chatterjee, R., Jana, J., and Paul, P. (2012). Enhancement of head yield and quality of cabbage (Brassica oleracea) by combining different sources of nutrients. Indian J. Agr. Sci. 82, 324–328.
Chen, C., Zhang, J., Lu, M., Qin, C., Chen, Y., Yang, L., et al. (2016). Microbial communities of an arable soil treated for 8 years with organic and inorganic fertilizers. Biol. Fert. Soils 52, 455–467. doi: 10.1007/s00374-016-1089-5
Chen, S., Zhou, Y., Chen, Y., and Gu, J. (2018). FASTP: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884–i890. doi: 10.1101/274100
Chowdhury, N., Marschner, P., and Burns, R. (2011). Response of microbial activity and community structure to decreasing soil osmotic and matric potential. Plant Soil 344, 241–254. doi: 10.1007/s11104-011-0743-9
Ding, J., Jiang, X., Guan, D., Zhao, B., Ma, M., Zhou, B., et al. (2017). Influence of inorganic fertilizer and organic manure application on fungal communities in a long-term field experiment of Chinese Mollisols. Appl. Soil Ecol. 111, 114–122. doi: 10.1016/j.apsoil.2016.12.003
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200.
Flores-Félix, J. D., Menéndez, E., Rivera, L. P., Marcos-García, M., Martínez-Hidalgo, P., Mateos, P. F., et al. (2013). Use of Rhizobium leguminosarum as a potential biofertilizer for Lactuca sativa and Daucus carota crops. J. Plant Nutr. Soil Sci. 176, 876–882. doi: 10.1002/jpln.201300116
Geng, Y., Cao, G., Wang, L., and Wang, S. (2019). Effects of equal chemical fertilizer substitutions with organic manure on yield, dry matter, and nitrogen uptake of spring maize and soil nitrogen distribution. PLoS One 14:e0219512. doi: 10.1371/journal.pone.0219512
Gruzdeva, L., Matveeva, E., and Kovalenko, T. (2007). Changes in soil nematode communities under the impact of fertilizers. Eurasian Soil Sci. 40, 681–693.
Guo, J., Liu, W., Zhu, C., Luo, G., Kong, Y., Ling, N., et al. (2018). Bacterial rather than fungal community composition is associated with microbial activities and nutrient-use efficiencies in a paddy soil with short-term organic amendments. Plant Soil 424, 335–349. doi: 10.1007/s11104-017-3547-8
Haiyan, C., Xiang, X., Jian, Y., Adams, J. M., Zhang, K., Yuntao, L., et al. (2016). Effects of slope aspects on soil bacterial and arbuscular fungal communities in a boreal forest in China. Pedosphere 26, 226–234. doi: 10.1016/S1002-0160(15)60037-6
Han, J., Dong, Y., and Zhang, M. (2021). Chemical fertilizer reduction with organic fertilizer effectively improve soil fertility and microbial community from newly cultivated land in the Loess Plateau of China. Appl. Soil Ecol. 165:103966. doi: 10.1016/j.apsoil.2021.103966
Hartmann, M., Frey, B., Mayer, J., Mäder, P., and Widmer, F. (2015). Distinct soil microbial diversity under long-term organic and conventional farming. ISME J. 9, 1177–1194. doi: 10.1038/ismej.2014.210
He, D., Xiang, X., He, J.-S., Wang, C., Cao, G., Adams, J., et al. (2016). Composition of the soil fungal community is more sensitive to phosphorus than nitrogen addition in the alpine meadow on the Qinghai-Tibetan Plateau. Biol. Fert. Soils 52, 1059–1072. doi: 10.1007/s00374-016-1142-4
He, J.-Z., Zheng, Y., Chen, C.-R., He, Y.-Q., and Zhang, L.-M. (2008). Microbial composition and diversity of an upland red soil under long-term fertilization treatments as revealed by culture-dependent and culture-independent approaches. J. Soil Sediment 8, 349–358. doi: 10.1007/s11368-008-0025-1
Hermans, S. M., Buckley, H. L., Case, B. S., Curran-Cournane, F., Taylor, M., and Lear, G. (2017). Bacteria as emerging indicators of soil condition. Appl. Environ. Microb. 83, e2826–e2816. doi: 10.1128/AEM.02826-16
Kaisermann, A., Maron, P., Beaumelle, L., and Lata, J. (2015). Fungal communities are more sensitive indicators to non-extreme soil moisture variations than bacterial communities. Appl. Soil Ecol. 86, 158–164. doi: 10.1016/j.apsoil.2014.10.009
Kamaa, M., Mburu, H., Blanchart, E., Chibole, L., Chotte, J.-L., Kibunja, C., et al. (2011). Effects of organic and inorganic fertilization on soil bacterial and fungal microbial diversity in the Kabete long-term trial. Kenya. Biol. Fert. Soils 47, 315–321. doi: 10.1007/s00374-011-0539-3
Kim, J. M., Roh, A.-S., Choi, S.-C., Kim, E.-J., Choi, M.-T., Ahn, B.-K., et al. (2016). Soil pH and electrical conductivity are key edaphic factors shaping bacterial communities of greenhouse soils in Korea. J. Microbiol. 54, 838–845. doi: 10.1007/s12275-016-6526-5
Krashevska, V., Klarner, B., Widyastuti, R., Maraun, M., and Scheu, S. (2015). Impact of tropical lowland rainforest conversion into rubber and oil palm plantations on soil microbial communities. Biol. Fert. Soils 51, 697–705. doi: 10.1007/s00374-015-1021-4
Lazcano, C., Gómez-Brandón, M., Revilla, P., and Domínguez, J. (2013). Short-term effects of organic and inorganic fertilizers on soil microbial community structure and function. Biol. Fert. Soils 49, 723–733. doi: 10.1007/s00374-012-0761-7
Liang, Q., Chen, H., Gong, Y., Yang, H., Fan, M., and Kuzyakov, Y. (2014). Effects of 15 years of manure and mineral fertilizers on enzyme activities in particle-size fractions in a North China Plain soil. Eur. J. Soil Biol. 60, 112–119. doi: 10.1016/j.ejsobi.2013.11.009
Lin, Y., Ye, G., Kuzyakov, Y., Liu, D., Fan, J., and Ding, W. (2019). Long-term manure application increases soil organic matter and aggregation, and alters microbial community structure and keystone taxa. Soil Biol. Biochem. 134, 187–196. doi: 10.1016/j.soilbio.2019.03.030
Liu, G.y, Chen, L.l, Shi, X. r, Yuan, Z.y, Yuan, L. Y., Lock, T. R., et al. (2019). Changes in rhizosphere bacterial and fungal community composition with vegetation restoration in planted forests. Land Degrad. Dev. 30, 1147–1157. doi: 10.1002/ldr.3275
Liu, J., Sui, Y., Yu, Z., Shi, Y., Chu, H., Jin, J., et al. (2015). Soil carbon content drives the biogeographical distribution of fungal communities in the black soil zone of northeast China. Soil Biol. Biochem 83, 29–39. doi: 10.1016/j.soilbio.2015.01.009
Luo, P., Han, X., Wang, Y., Han, M., Shi, H., Liu, N., et al. (2015). Influence of long-term fertilization on soil microbial biomass, dehydrogenase activity, and bacterial and fungal community structure in a brown soil of northeast China. Ann. Microbiol. 65, 533–542. doi: 10.1007/s13213-014-0889-9
Luthria, D., Singh, A. P., Wilson, T., Vorsa, N., Banuelos, G. S., and Vinyard, B. T. (2010). Influence of conventional and organic agricultural practices on the phenolic content in eggplant pulp: plant-to-plant variation. Food Chem. 121, 406–411. doi: 10.1016/j.foodchem.2009.12.055
Lv, F., Song, J., Giltrap, D., Feng, Y., Yang, X., and Zhang, S. (2020). Crop yield and N2O emission affected by long-term organic manure substitution fertilizer under winter wheat-summer maize cropping system. Sci. Total Environ. 732:139321. doi: 10.1016/j.scitotenv.2020.139321
Magoč, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507
Miller, M., Zebarth, B., Dandie, C., Burton, D., Goyer, C., and Trevors, J. (2008). Crop residue influence on denitrification. Soil Biol. Biochem. 40, 2553–2562. doi: 10.1016/j.soilbio.2008.06.024
Morris, S. J., and Blackwood, C. B. (2015). “The ecology of the soil biota and their function,” in Soil Microbiology, Ecology and Biochemistry, ed. E. A. Paul (London: Elsevier), 273–309. doi: 10.1016/B978-0-12-415955-6.00010-4
Naeem, M., Iqbal, J., and Bakhsh, M. (2006). Comparative study of inorganic fertilizers and organic manures on yield and yield components of mungbean (Vigna radiat L.). J. Agric. Soc. Sci. 2, 227–229.
Nemergut, D. R., Cleveland, C. C., Wieder, W. R., Washenberger, C. L., and Townsend, A. R. (2010). Plot-scale manipulations of organic matter inputs to soils correlate with shifts in microbial community composition in a lowland tropical rain forest. Soil Biol. Biochem. 42, 2153–2160. doi: 10.1016/j.soilbio.2010.08.011
Nie, C., Zhu, P., Wang, M., Ma, S., and Wei, Z. (2017). Optimization of water-soluble polysaccharides from stem lettuce by response surface methodology and study on its characterization and bioactivities. Int. J. Biol. Macromol. 105, 912–923. doi: 10.1016/j.ijbiomac.2017.07.125
Nilsson, R. H., Larsson, K.-H., Taylor, A. F. S., Bengtsson-Palme, J., Jeppesen, T. S., Schigel, D., et al. (2019). The UNITE database for molecular identification of fungi: handling dark taxa and parallel taxonomic classifications. Nucleic Acids Res. 47, D259–D264. doi: 10.1093/nar/gky1022
Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P., O’Hara, R., et al. (2015). Vegan Community Ecology Package: Ordination Methods, Diversity Analysis and Other Functions for Community and Vegetation Ecologists. R Package ver, 2-3.
Oliveira, A. B., Moura, C. F., Gomes-Filho, E., Marco, C. A., Urban, L., and Miranda, M. R. A. (2013). The impact of organic farming on quality of tomatoes is associated to increased oxidative stress during fruit development. PLoS One 8:e56354. doi: 10.1371/journal.pone.0056354
Paungfoo-Lonhienne, C., Yeoh, Y. K., Kasinadhuni, N. R. P., Lonhienne, T. G., Robinson, N., Hugenholtz, P., et al. (2015). Nitrogen fertilizer dose alters fungal communities in sugarcane soil and rhizosphere. Sci. Rep. UK 5:8678. doi: 10.1038/srep08678
Pinggera, J., Geisseler, D., Piepho, H.-P., Joergensen, R. G., and Ludwig, B. (2015). Effect of substrate quality on the N uptake routes of soil microorganisms in different soil depths. Pedobiologia 58, 211–218. doi: 10.1016/j.pedobi.2015.10.002
Qi, Y., Jiang, F., Zhou, R., Wu, Y., Hou, X., Li, J., et al. (2021). Effects of reduced nitrogen with bio-organic fertilizer on soil properties, yield and quality of non-heading Chinese Cabbage. Agronomy 11:2196. doi: 10.3390/agronomy11112196
Qu, T.-b, Du, W.-c, Yuan, X., Yang, Z.-m, Liu, D.-b, Wang, D.-l, et al. (2016). Impacts of grazing intensity and plant community composition on soil bacterial community diversity in a steppe grassland. PLoS One 11:e0159680. doi: 10.1371/journal.pone.0159680
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
Rousk, J., Brookes, P. C., and Baath, E. (2009). Contrasting soil pH effects on fungal and bacterial growth suggest functional redundancy in carbon mineralization. Appl. Environ. Microb. 75, 1589–1596. doi: 10.1128/AEM.02775-08
Sarkar, D., and Rakshit, A. (2021). Bio-priming in combination with mineral fertilizer improves nutritional quality and yield of red cabbage under Middle Gangetic Plains India. Sci. Hortic. 283:110075. doi: 10.1016/j.scienta.2021.110075
Sedmak, J. J., and Grossberg, S. E. (1977). A rapid, sensitive, and versatile assay for protein using Coomassie brilliant blue G250. Anal. Biochem. 79, 544–552. doi: 10.1016/0003-2697(77)90428-6
Shen, J.-P., Zhang, L.-M., Guo, J.-F., Ray, J. L., and He, J.-Z. (2010). Impact of long-term fertilization practices on the abundance and composition of soil bacterial communities in Northeast China. Appl. Soil Ecol. 46, 119–124. doi: 10.1016/j.apsoil.2010.06.01
Shen, Z., Xue, C., Penton, C. R., Thomashow, L. S., Zhang, N., Wang, B., et al. (2019). Suppression of banana Panama disease induced by soil microbiome reconstruction through an integrated agricultural strategy. Soil Biol. Biochem. 128, 164–174. doi: 10.1016/j.soilbio.2018.10.016
Strickland, M. S., and Rousk, J. (2010). Considering fungal: bacterial dominance in soils–methods, controls, and ecosystem implications. Soil Biol. Biochem. 42, 1385–1395. doi: 10.1016/j.soilbio.2010.05.007
Sun, R., Zhang, X.-X., Guo, X., Wang, D., and Chu, H. (2015). Bacterial diversity in soils subjected to long-term chemical fertilization can be more stably maintained with the addition of livestock manure than wheat straw. Soil Biol. Biochem. 88, 9–18. doi: 10.1016/j.soilbio.2015.05.007
Tao, R., Liang, Y., Wakelin, S. A., and Chu, G. (2015). Supplementing chemical fertilizer with an organic component increases soil biological function and quality. Appl. Soil Ecol. 96, 42–51. doi: 10.1016/j.apsoil.2015.07.009
Tilman, D., Balzer, C., Hill, J., and Befort, B. L. (2011). Global food demand and the sustainable intensification of agriculture. P. Natl. Acad. Sci. U S A. 108, 20260–20264. doi: 10.1073/pnas.1116437108
Vallverdú-Queralt, A., Medina-Remón, A., Casals-Ribes, I., and Lamuela-Raventos, R. M. (2012). Is there any difference between the phenolic content of organic and conventional tomato juices? Food Chem. 130, 222–227. doi: 10.1016/j.foodchem.2011.07.017
Victoria, R., Banwart, S., Black, H., Ingram, J., Joosten, H., Milne, E., et al. (2012). The Benefits of Soil Carbon. F, Cheltenham: Edward Elgar Pub. Ltd. 2012, 19–33.
Wang, C., Lv, J., Coulter, J. A., Xie, J., Yu, J., Li, J., et al. (2020). Slow-release fertilizer improves the growth, quality, and nutrient utilization of wintering Chinese chives (Allium tuberosum Rottler ex Spreng.). Agronomy 10:381. doi: 10.3390/agronomy10030381
Wang, Y., Hu, N., Ge, T., Kuzyakov, Y., Wang, Z.-L., Li, Z., et al. (2017b). Soil aggregation regulates distributions of carbon, microbial community and enzyme activities after 23-year manure amendment. Appl. Soil Ecol. 111, 65–72. doi: 10.1016/j.apsoil.2016.11.015
Wang, J., Song, Y., Ma, T., Raza, W., Li, J., Howland, J. G., et al. (2017a). Impacts of inorganic and organic fertilization treatments on bacterial and fungal communities in a paddy soil. Appl. Soil Ecol. 112, 42–50. doi: 10.1016/j.apsoil.2017.01.005
Wang, J., Zhang, L., Lu, Q., Raza, W., Huang, Q., and Shen, Q. (2014). Ammonia oxidizer abundance in paddy soil profile with different fertilizer regimes. Appl. Soil Ecol. 84, 38–44. doi: 10.1016/j.apsoil.2014.06.009
Wang, Z. H., Li, S. X., and Malhi, S. (2008). Effects of fertilization and other agronomic measures on nutritional quality of crops. J. Sci. Food Agr. 88, 7–23. doi: 10.1002/jsfa.3084
Wolińska, A., Kuźniar, A., Zielenkiewicz, U., Izak, D., Szafranek-Nakonieczna, A., Banach, A., et al. (2017). Bacteroidetes as a sensitive biological indicator of agricultural soil usage revealed by a culture-independent approach. Appl. Soil Ecol. 119, 128–137.
Wolińska, A., Podlewski, J., Słomczewski, A., Grządziel, J., Gałązka, A., and Kuźniar, A. (2021). Fungal indicators of sensitivity and resistance to long-term maize monoculture: a culture-independent approach. Front. Microbiol. 12:799378. doi: 10.3389/fmicb.2021.799378
Wu, L., Jiang, Y., Zhao, F., He, X., Liu, H., and Yu, K. (2020). Increased organic fertilizer application and reduced chemical fertilizer application affect the soil properties and bacterial communities of grape rhizosphere soil. Sci. Rep. 10:9568. doi: 10.1038/s41598-020-66648-9
Wu, L., Li, Z., Zhao, F., Zhao, B., Phillip, F. O., Feng, J., et al. (2021a). Increased organic fertilizer and reduced chemical fertilizer increased fungal diversity and the abundance of beneficial fungi on the grape berry surface in arid areas. Front. Microbiol. 12:1066. doi: 10.3389/fmicb.2021.628503
Wu, X., Zhang, T., Zhao, J., Wang, L., Yang, D., Li, G., et al. (2021b). Variation of soil bacterial and fungal communities from fluvo-aquic soil under chemical fertilizer reduction combined with organic materials in North China Plain. J. Plant Nutr. Soil Sci. 21, 349–363. doi: 10.1007/s42729-020-00365-0
Wu, Y., Zhao, C., Farmer, J., and Sun, J. (2015). Effects of bio-organic fertilizer on pepper growth and Fusarium wilt biocontrol. Sci. Hortic. 193, 114–120. doi: 10.1016/j.scienta.2015.06.039
Xiao, H. F., Feng, Y. L., Schaefer, D. A., and Yang, X. D. (2014). Soil fungi rather than bacteria were modified by invasive plants, and that benefited invasive plant growth. Plant soil 378, 253–264. doi: 10.1007/s11104-014-2040-x
Xiao, L., Sun, Q., Yuan, H., and Lian, B. (2017). A practical soil management to improve soil quality by applying mineral organic fertilizer. Acta Geochimica 36, 198–204. doi: 10.1007/s11631-017-0139-5
Xiaojing, H., Junjie, L., Dan, W., Ping, Z., Xi’an, C., Baoku, Z., et al. (2018). Soil bacterial communities under different long-term fertilization regimes in three locations across the black soil region of Northeast China. Pedosphere 28, 751–763. doi: 10.1016/S1002-0160(18)60040-2
Xue, F., Yan, T., Yang, L., and Qiao, J. (2010). Influences of organic fertilizer application on soil biological properties. Zhongguo Shengtai Nongye Xuebao Chin. J. Eco Agric. 18, 1372–1377. doi: 10.3724/SP.J.1011.2010.01372
Yang, F., Tian, J., Fang, H., Gao, Y., Xu, M., Lou, Y., et al. (2019a). Functional soil organic matter fractions, microbial community, and enzyme activities in a mollisol under 35 years manure and mineral fertilization. J. Plant Nutr. Soil Sc. 19, 430–439. doi: 10.1007/s42729-019-00047-6
Yang, L., Bai, J., Zeng, N., Zhou, X., Liao, Y., Lu, Y., et al. (2019c). Diazotroph abundance and community structure are reshaped by straw return and mineral fertilizer in rice-rice-green manure rotation. Appl. Soil Ecol. 136, 11–20. doi: 10.1016/j.apsoil.2018.12.015
Yang, H., Ma, J., Rong, Z., Zeng, D., Wang, Y., Hu, S., et al. (2019b). Wheat straw return influences nitrogen-cycling and pathogen associated soil microbiota in a wheat–soybean rotation system. Front. Microbiol. 10:1811. doi: 10.3389/fmicb.2019.01811
Yang, Y., Li, X., Liu, J., Zhou, Z., Zhang, T., and Wang, X. (2020b). Fungal community structure in relation to manure rate in red soil in southern China. Appl. Soil Ecol. 147:103442. doi: 10.1016/j.apsoil.2019.103442
Yang, P., Lyu, J., Sohail, H., Yu, J., Xie, J., and Li, J. (2020a). Partial substitution of mineral fertilizer with biofertilizer enhances cauliflower nutritional quality, yield, and soil characteristics. Crop Sci. 60, 934–944. doi: 10.1002/csc2.20023
Yang, Y., Wang, Z., Hu, Y., and Zeng, Z. (2017). Irrigation frequency alters the abundance and community structure of ammonia-oxidizing archaea and bacteria in a northern Chinese upland soil. Eur. J. Soil Biol. 83, 34–42. doi: 10.1016/j.ejsobi.2017.09.005
Ye, G., Lin, Y., Luo, J., Di, H. J., Lindsey, S., Liu, D., et al. (2020a). Responses of soil fungal diversity and community composition to long-term fertilization: field experiment in an acidic Ultisol and literature synthesis. Appl. Soil Ecol. 145:103305. doi: 10.1016/j.apsoil.2019.06.008
Ye, L., Zhao, X., Bao, E., Li, J., Zou, Z., and Cao, K. (2020b). 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
Yelle, D. J., Ralph, J., Lu, F., and Hammel, K. E. (2008). Evidence for cleavage of lignin by a brown rot basidiomycete. Environ. Microbiol. 10, 1844–1849. doi: 10.1111/j.1462-2920.2008.01605.x
Yu, L., Nicolaisen, M., Larsen, J., and Ravnskov, S. (2013). Organic fertilization alters the community composition of root associated fungi in Pisum sativum. Soil Biol. Biochem. 58, 36–41. doi: 10.1016/j.soilbio.2012.11.004
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, C., Fu, S., Mathew, R. P., Lawrence, K. S., and Feng, Y. (2015). Soil microbial community structure and activity in a 100-year-old fertilization and crop rotation experiment. J. Plant Ecol. 8, 623–632. doi: 10.1093/jpe/rtv007
Zhao, J., Liu, J., Liang, H., Huang, J., Chen, Z., Nie, Y., et al. (2018). Manipulation of the rhizosphere microbial community through application of a new bio-organic fertilizer improves watermelon quality and health. PLoS One 13:e0192967. doi: 10.1371/journal.pone.0192967
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
Keywords: chemical fertilizer, bio-organic fertilizer, high-throughput sequencing, bacterial community, fungal community
Citation: Jin N, Jin L, Wang S, Li J, Liu F, Liu Z, Luo S, Wu Y, Lyu J and Yu J (2022) Reduced Chemical Fertilizer Combined With Bio-Organic Fertilizer Affects the Soil Microbial Community and Yield and Quality of Lettuce. Front. Microbiol. 13:863325. doi: 10.3389/fmicb.2022.863325
Received: 27 January 2022; Accepted: 07 March 2022;
Published: 21 April 2022.
Edited by:
Jesús Navas-Castillo, La Mayora Experimental Station (CSIC), SpainReviewed by:
Agnieszka Wolińska, The John Paul II Catholic University of Lublin, PolandJosé David Flores Félix, Universidade da Beira Interior, Portugal
Copyright © 2022 Jin, Jin, Wang, Li, Liu, Liu, Luo, Wu, Lyu and Yu. 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: Jian Lyu, bHZqaWFuZ3NAMTI2LmNvbQ==; Jihua Yu, eXVqaWh1YWdnQDE2My5jb20=