Next Article in Journal
Understanding the Dynamic of Rice Farming Systems in Southern Mozambique to Improve Production and Benefits to Smallholders
Next Article in Special Issue
Nitrogen Immobilisation and Microbial Biomass Build-Up Induced by Miscanthus x giganteus L. Based Fertilisers
Previous Article in Journal
Sulfonamides in Tomato from Commercial Greenhouses Irrigated with Reclaimed Wastewater: Uptake, Translocation and Food Safety
Previous Article in Special Issue
Short-Term Effect of Biochar on Microbial Biomass, Respiration and Enzymatic Activities in Wastewater Irrigated Soils in Urban Agroecosystems of the West African Savannah
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Manure and Chemical Fertilizer on Bacterial Community Structure and Soil Enzyme Activities in North China

1
College of Resources and Environment, Shanxi Agricultural University, Taiyuan 030031, China
2
Shanxi Institute of Organic Dryland Farming, Shanxi Agricultural University, Taiyuan 030031, China
*
Author to whom correspondence should be addressed.
Submission received: 24 March 2021 / Revised: 23 April 2021 / Accepted: 13 May 2021 / Published: 20 May 2021
(This article belongs to the Special Issue Soil Microbes and Nutrient Recycling in Agroecosystems)

Abstract

:
The application of organic fertilizer affects soil microbes and enzyme activities. In this study, we explored the effects of various long-term different fertilization treatments (manure, M; chemical fertilizer, NP; manure + chemical fertilizer, MNP; and no fertilizer, CK) on bacterial community structure and soil sucrase, urease, and alkaline phosphatase activities in Shaping, Hequ, China. High-throughput sequencing was used to amplify the third to the fourth hypervariable region of the 16S ribosomal RNA for analysis of the bacterial community structure. Enzyme activities were determined by colorimetry. Soil treated with MNP had the highest bacterial Abundance-based Coverage Estimator index and enzyme activities. The principal coordinates analysis results showed significant differences among the various fertilization treatments (p < 0.001). Proteobacteria, Actinobacteria, Acidobacteria, Gemmatimonadetes, and Chloroflexi were consistently dominant in all soil samples. The redundancy analysis and Monte Carlo permutation tests showed that the soil bacterial communities were significantly correlated with alkali-hydrolyzable nitrogen, organic matter, urease, and alkaline phosphatase. Our results reveal the fundamentally different effects that organic and inorganic fertilizers have on soil bacterial communities and their functions.

1. Introduction

Soil microorganisms and enzyme activities are important indicators in the characterization of soil fertility and these play a vital role in soil material transformation and energy flow [1,2]. Soil bacteria are one of the largest functional microbial subgroups in soil, participating in soil respiration, nutrient transformation, organic matter decomposition, and other processes [3,4]. Soil enzymes secreted by microorganisms play critical roles in nutrient cycling (such as nitrogen fixation, phosphorus adsorption and desorption, potassium release, etc.), as well as soil structure maintenance and crop production [5,6]. Organic matter in soil serves as a substrate for a diverse set of soil microorganisms, promoting enzymatic activity [2]. Previous studies have shown that soil microorganisms and enzyme activities are highly sensitive to changes in soil properties, including pH, the degree of organic matter, moisture content, the presence of mineral substances, etc. Furthermore, soil microorganisms are also regarded as vital factors driving the production and turnover of extracellular enzymes and their activities in soil [5,7]. Studies showed that the physical and chemical properties of soil could be largely influenced by fertilization management practices [8,9].
Fertilization is an important measure in agricultural production, which can improve soil nutrients and increase crop yield. Recently, a large number of chemical fertilizers have been imported into farmland to ensure food supply, and the result was an imbalance in soil nutrients, as well as a series of environmental problems [8,10]. Manure, however, could partly prevent many environmental problems caused by chemical fertilizers and also improve crop quality. It might be due to the fact that numerous beneficial microorganisms and organic materials enter soil through manure application, leading to changes in the soil environment that include differences in enzyme activity. Thus, it is necessary to investigate the soil’s bacterial community structure and soil enzymatic activity. One of the main challenges in microbiological research is to improve the identification of bacteria. Traditional methods have various defects. For instance, the plating method can only reflect culturable microorganisms that account for only about 1% of the total population [11,12]. Furthermore, neither single-strand conformation polymorphism (SSCP) nor terminal restriction fragment length polymorphism (T-RFLP) are suitable for bacterial analysis because they only capture the dominant community members as selected by their PCR primers [13]. High-throughput sequencing technology, however, has facilitated bacterial research, allowing for huge amounts of information to be obtained, and for the bacterial diversity in the environment to be fully reflected [14]. In this study, we explored the effects of the long-term application of manure and chemical fertilizer on bacterial community structure with high-throughput sequencing technology, and we studied the soil’s enzymatic activities by colorimetry in North China. The objectives of this study were as follows: (1) to elucidate the influences of different fertilization on bacterial community structure and soil enzyme activities, (2) and to search for the most suitable fertilization method.

2. Materials and Methods

2.1. Field Description and Experimental Design

The fertilization experiment was established in a maize field in Shaping, Hequ, China (39°12′18″ N, 111°15′41″ E). The altitude is 1089 m, with a frost-free period of about 140 days. The average temperature is 8.8 ℃ with an effective accumulative temperature of 3000–3360 ℃ per year. The experimental site has a light loam loess, classified as Hapli-Ustic Cambosols in the World Reference Base for Soil Resources (WRB) [15,16]. Analysis of soil samples taken from the experimental area showed that the basic physical and chemical properties of the surface 0–20 cm of soil before planting were as follows: organic matter (OM), 5.64 g/kg; total nitrogen (TN), 0.45 g/kg; total phosphorus (TP), 1.23 g/kg; alkali-hydrolyzable nitrogen (AN), 34.90 mg/kg; available phosphorus (AP), 2.69 mg/kg; cation exchange capacity (CEC), 6.97 cmol/kg; pH, 8.2.
The long-term fertilization trial began in 1988, and it lasted for 31 years. From 1988 to 2008, millet (Panicum miliaceum L.) and potatoes (Solanum tuberosum L.) were rotated; from 2009 to 2013, maize (Zea mays L.) was continuously cropped; from 2014 to 2017, potatoes and millet were rotated; and since 2018, potatoes and maize have been periodically rotated. The samples were obtained on 25 July 2019, when the crop was maize. The four treatments were arranged in a randomized block design with three replications, with each plot being 24 m2. The four treatments were (1) M—manure, 22,500 kg/ha; (2) NP—N: 120 kg/ha, P2O5: 75 kg/ha; (3) MNP—manure, 22,500 kg/ha, N: 120 kg/ha, P2O5: 75 kg/ha; and (4) CK—control treatment without fertilizer. Nitrogen was applied in the form of urea, while phosphorus was applied as superphosphate. The manure was made up of local cow dung with an average N content of 3.64 g/kg, P2O5 content of 2.46 g/kg, and K2O content of 7.87 g/kg. The fertilization details are shown in Table 1. All of the fertilizers were applied as basal fertilizers before planting.

2.2. Soil Sampling and Analysis of Soil Properties

Soil were sampled from five points in each plot and mixed into one sample. There were 12 composite samples in total. After removing visible stone and debris, samples were collected in sterile plastic bags in an ice-box and transported to the laboratory. One part was used for analyzing physicochemical properties and enzyme activities after sieved through a 2 mm screen, and the other part was stored in −80℃ refrigerator for DNA extraction and sequencing.

2.3. Soil Physicochemical Analytical Procedures

Soil TN, TP, and total potassium (TK) were analyzed using a Vario Max element analyzer (Elementar Vario PYRO cube and Isoprime100, Hanau, Germany). Soil pH was determined at a soil/water ratio of 1:2.5 with a pH meter. OM was determined by K2Cr2O7 oxidation method. Soil AN was determined by the alkaline hydrolysis diffusion method. AP was determined by NaHCO3-extracted method. Available potassium (AK) was analyzed by CH3COONH4-extracted method. The details can be found from the reference [17].

2.4. Soil Enzyme Activities’ Determination

In this study, 3, 5-dinitrosalicylic acid colorimetry was used for sucrase (SUC) activity, and the weight of glucose (mg) in 1 g of soil after incubation for 24 h was used to represent SUC activity. The urease (URE) activity was indicated by sodium phenol-sodium hypochlorite colorimetric method, and the weight of NH3-N (mg) in 1 g soil after incubation for 24 h was used to represent the URE activity. The alkaline phosphatase (ALP) activity was represented by the phenyl disodium phosphate method. Detailed procedures for the three enzymes can be found from the reference [18].

2.5. Soil DNA Extraction and High-Throughput Sequencing

Total DNA of each sample was extracted form 0.5 g soil using the Fast DNA SPIN Isolation Kit (MP Biomedicals, Santa Ana, CA, USA) according to the manufacturer’s instructions. The bacterial third to the fourth hypervariable region of the 16S ribosomal RNA (16S rRNA) was amplified using primers 338F (5’-ACTCCTACGGGAGGCAGCA-3’) and 806R (5’-GGACTACHVGGGTWTCTAAT-3’). Specific barcodes of 7-bp were incorporated into the primers for multiplex sequencing. The first round PCR reactions were performed in a total volume of 50 μL mixture containing 10 μL of 5×Q5 Reaction Buffer, 10 μL of 5 × Q5 High-Fidelity GC Enhancer, 1 μL of dNTPs, 1 μL of each primer (10 μM), 0.2 μL of Q5 High-Fidelity DNA Polymerase (Transgen, Beijing, China), 2 μL of temple DNA and 24.8 μL of ddH2O. PCR process were as follows: 5 min of initial denaturation at 95 °C, followed by 15 cycles of denaturation at 95 °C for 1 min, annealing at 50 °C for 1 min and extension at 72 °C for 1 min, and a final extension of 7 min at 72 °C.
PCR products from the first PCR step were purified through VAHTSTM DNA Clean Beads (Vazyme, Nanjing, China). A second round PCR was then performed in a 40 μL reaction system which contained 20 μL 2 × Phusion High-Fidelity Master Mix (New England Biolabs, Boston, MA, USA), 8 μL ddH2O, 10 μM of each primer and 10 μL PCR products from the first step. PCR process was as follows: an initial denaturation at 98 °C for 30 s, followed by 10 cycles of denaturation at 98 °C for 10 s, annealing at 65 °C for 30 s and extension 72 °C for 30 s, and a final extension at 72 °C for 5 min. Finally, all PCR products were quantified by Quant-iT™ dsDNA HS Reagent (Thermo Fisher, Waltham, MA, USA) and pooled together. High-throughput sequencing analysis of bacterial 16S rRNA genes was performed with the purified and pooled sample using the Illumina Hiseq 2500 platform (2 × 250 paired ends) at Biomarker Technologies Corporation (Beijing, China).

2.6. Bioinformatics Analysis

The Quantitative Insights Into Microbial Ecology (QIIME, version 1.8.0) pipeline was employed to process the sequencing data [19]. Briefly, raw sequencing reads with exact matches to the barcodes were assigned to respective samples and identified as valid sequences. The low-quality sequences were removed. Paired-end reads were assembled using FLASH (version 1.2.7) [20]. After chimeras’ detection with UCHIME (version 4.2), the remaining high-quality sequences were clustered into operational taxonomic units (OTUs) at a 97% sequence identity threshold by UCLUST [21]. A representative sequence was selected from each OTU using default parameters. OTU taxonomic classification was conducted by BLAST searching the representative sequences set against the Silva Database (version 132) [22]. An OTU table was further generated to record the taxonomy of the OTUs and the abundance of each OTU in each sample. OTUs containing less than 0.001% of total sequences across all samples were discarded. To minimize sequencing depth differences across samples, an averaged, rounded, and rarefied OTU table was generated by averaging 100 evenly resampled OTU subsets under the 90% minimum sequencing depth for further analysis.
Sequence data analyses were mainly performed using QIIME and R (version 3.1.2, vegan package) [23,24]. The Abundance-based Coverage Estimator (ACE) and Shannon index were calculated based on the OTU table using QIIME. β-diversity analysis was performed to investigate the bacterial communities’ structural variation across samples using principal coordinate analysis (PCoA), which was based on the OTU level compositional profiles [25]. The significance of bacteria structure differentiation among groups was assessed by permutational multivariate analysis of variance (PERMANOVA) [26] with the distance legend of 0.01. The taxonomic compositions and abundances were visualized using MEGAN [27]. A Venn diagram was generated to visualize the shared and unique OTUs among groups based on the occurrence of OTUs across groups regardless of their relative abundances [28]. Redundancy analysis (RDA) and Monte Carlo permutation tests were carried out using Canoco for Windows 4.5 (Cabiy Information Techonogy Co., LTD., Shanghai, China). The species compositional similarity was analyzed by unweighted pair-group method with arithmetic mean (UPGMA) with R software.

2.7. Statistical Analysis

Statistical analyses were carried out using a one-way analysis of variance (ANOVA) procedure in the SPSS software (IBM SPSS Statistics version 22) to check the normal distribution and homoscedasticity, and to detect differences in soil parameters, soil enzyme activities, the bacterial relative abundance, α-diversity indices, etc. Significant differences between data were determined with the least significant difference (LSD) test at the p = 0.05 level.

3. Results

3.1. Soil Physicochemical Characteristics

As shown in Table 2, long-term fertilization significantly promoted soil TP, AN and AP (p < 0.05) compared to the control. The TN level was the highest in treatment MNP, followed by treatment M, and then treatment NP, while M was not significantly different from MNP. The TP value was the highest in treatment MNP, followed by treatment NP and treatment M. Again, however, there was no significant difference between NP and MNP. Soil pH (H2O) values decreased with the application of fertilizers, although no significant differences were observed. F- and p-values of the factors can be found in Table S1.

3.2. Bacterial α-Diversity

After the removal of chimeras, a total number of 748,247 sequences were obtained in this study, which are available through the NCBI Sequence Read Archive (Accession: PRJNA706481). The OTU numbers in the samples did not significantly increase with the number of sequences sampled and the curve gradually flattened out (Supplementary Figure S1), indicating that the sample size was large enough and that the library could satisfactorily characterize the bacterial communities. The ACE and Shannon indices under different treatments were estimated. As shown in Table 3, the ACE was the highest in treatment MNP with a combination of manure and chemical fertilizer, followed by treatments M, NP, and CK, respectively. The Shannon index in the M treatment was significantly higher than in the CK treatment. However, no significant differences were observed between the other treatments.

3.3. Soil Bacterial Community Composition

As shown in Figure 1, 1701 OTUs were shared by different treatments, accounting for 23.4% of the total OTU of all the samples. Numbers of specific OTUs of treatments M, NP, MNP, and CK were 1, 0, 25, and 0, respectively. The combined application of manure and chemical fertilizer notably increased the number of unique OTUs in treatment MNP compared with the other treatments.

3.4. Effect of Different Treatments on Soil Bacterial Communities

PCoA was used to compare the similarity of soil bacterial communities among different treatments at the OTU level. As shown in Figure 2, the plot identified two principal component factors related to the percentage abundance of groups, explaining 63.91% (PC1) and 7.99% (PC3) of the total variation. PCoA showed that the four treatments could be separated clearly. The three repeats of each treatment were clustered together, showing good repeatability. Samples of treatment NP were parted from the other treatments by the first axis, while samples of treatments NP and CK were separated from samples of treatments M and MNP by the second axis. Samples of treatments M and MNP tended to cluster together, which could also be seen in the clustering tree (Supplementary Figure S2). In Figure 3, PERMANOVA showed significant differences in the bacterial community structures among treatments M, NP, MNP, and CK at the OTU level (p < 0.001).

3.5. Taxonomic Composition Analysis at the Phylum and Genus Level

The taxonomic distributions of bacterial communities were evaluated at different levels of classification. As shown in Figure 4, Proteobacteria, Actinobacteria, Acidobacteria, Gemmatimonadetes, and Chloroflexi were the most abundant phyla, accounting for approximately 90% of the bacteria detected in all 12 soil samples. Compared with treatment CK, the percentage of Proteobacteria increased by 2.49% in treatment M, 0.08% in treatment NP, and 1.94% in treatment MNP. The percentage of Gemmatimonadetes increased in treatments M, NP, and MNP by 0.29%, 1.58%, and 0.45%, respectively. However, compared with CK, the relative abundance of Chloroflexi decreased in treatments M, NP, and MNP by 0.97%, 0.53%, and 1.11%, respectively. The bacterial difference significance analysis can be found in Table S2.
As shown in Figure 5, at genus level, the top ten genera accounted for about 32% of the total genera and belonged to four phyla: Gemmatimonadetes, Acidobacteria, Proteobacteria, and Actinobacteria. Compared with CK, the percentage of MND1 decreased by 1.42% in treatment M, 1.56% in treatment NP, and 2.26% in treatment MNP, respectively. The relative abundance of RB41 decreased by 0.23%, 0.86%, and 0.85% in treatments M, NP, and MNP, respectively. On the contrast, the percentage of Arthrobacter increased by 0.03% in treatment M, 0.88% in treatment NP, and 0.29% in treatment MNP, respectively. The relative abundance of Steroidobacter increased by 0.60%, 0.21%, and 0.74% in treatments M, NP, and MNP, respectively. The bacterial difference significance analysis can be found in Table S3.

3.6. Effect of Different Treatments on Soil Enzyme Activities

Fertilization significantly increased soil enzyme activities compared with CK (Table 4). Three enzyme activities were the highest in MNP treatment, in which the nutrients were likewise the most abundant, followed by the M and NP treatments. There was no significant difference in the activities of the three soil enzymes between treatments M and NP, but they were slightly higher in the M treatment than in the NP treatment. Compared with CK, SUC activity increased in treatments M, NP, and MNP by 124.1%, 80.9%, and 145.6%, respectively. The increasing range of URE and ALP were 63.9–109.7% and 4.11–67.9%.

3.7. Correlations

We performed RDA to understand the correlations between phylum-level bacterial community composition and soil properties and that of enzyme activities. These eight environmental variables together explained 82.9% of the bacterial community variances. In addition, the first axis of the RDA accounted for a higher variation of 62.9%, whereas the second axis accounted for a lower of 20.0% (Figure 6a). The bacterial community composition had significant correlation with AN (F = 10.892, p = 0.002) and OM (F = 10.847, p = 0.002).
Figure 6b illustrated the relationship between enzyme activities and soil characteristics. The first and second axes explained 91.2% and 3.0% of the total change of the soil enzyme activities. The enzyme activities showed significant correlation with AP (F = 44.136, p = 0.002) and AN (F = 34.551, p = 0.002).
RDA between phylum-level bacterial community composition and soil enzyme activities was illustrated in Figure 6c and soil enzyme activities were taken as explanatory variables. The first and second axes explained 64.9% and 8.9% of the total changes of the bacterial community. The bacterial community composition showed significant correlation with URE (F = 15.559, p = 0.002) and ALP (F = 14.798, p = 0.002).

4. Discussion

4.1. Effects of Treatments on the Soil Physicochemical Properties

In the present study, no fertilizer had been applied in treatment CK for 31 years, leading to the lowest chemical properties. Long-term fertilization improved soil nutrients including AN, AP, AK, OM, TN, and TP. In this study, fertilization lowered soil pH, which was consistent with the findings of Wang, who reported a decline in soil pH with the application of fertilizer [29]. pH was regarded as one of the main factors related to fertilization management. A large number of studies have shown that soil pH decreased with excessive application of chemical N fertilizer, resulting in soil acidification [30,31]. In this study, no significant differences were observed in these treatments, which might be due to the fact that the calcareous soils containing higher carbonates to buffer soil pH [32].

4.2. Effects of Treatments on the Bacterial Diversity and Community Composition

Alpha diversity analysis showed significant difference (p < 0.05) of ACE index under different fertilization measures, and these were consistent with the previous results [8]. The combination of manure and chemical fertilizer significantly increased the bacterial ACE index compared with the other fertilization measures. The main reason was that organic materials themselves contain a variety of nutrient elements, which could increase the exogenous carbon, promote the formation of aggregate structure, and provide good conditions for the growth and reproduction of bacteria [33].
PCoA and PERMANOVA showed significant differences (p < 0.001) in bacterial community structure under different fertilization treatments, consistent with the results of Liang, who documented a clear separation of bacteria under different fertilization modes [34]. Fertilization directly affected the soil’s bacterial community structure by changing nutrients in the soil and by further impacting the biological activity of bacteria in soil. Additionally, the change of soil pH caused by fertilization had an indirect effect on the bacterial community structure. This change might be due to the fact that the optimal pH range for bacterial growth and reproduction was narrow [35]. The application of manure could cause significant differences on bacterial community structures compared with chemical fertilizer [36].
The availability of soil nutrients was usually linked with the transformation from oligotrophic microorganisms to copiotrophic microorganisms. For example, a high nutrient availability in soil promoted the growth and reproduction of copiotrophic microorganisms, while in nutrient-restricted soil environments, the number of oligotrophic microorganisms that grow slowly would increase [37,38]. Proteobacteria, Actinobacteria, Acidobacteria, Gemmatimonadetes, and Chloroflexi were the most abundant phyla, accounting for approximately 90% of the total bacteria, consistent with the previous results reported in different systems [39,40]. Among these phyla, Proteobacteria is considered to be a kind of copiotrophic microorganism widely existing in soil and whose abundance is relatively higher in treatments with manure. Most of the taxa within Acidobacteria are thought to be slow growing, oligotrophic microorganisms [37]. In the present study, the relative abundance of Acidobacteria was higher in CK without fertilizer than in NP and MNP by 2.16% and 2.34%, respectively, which indicated that Acidobacteria was more suitable to nutrient-restricted soil environments [38]. Gemmatimonadetes is dominant in soils with low water content [41]. The NP treatment in this study seemed to be more suitable for its survival. Chloroflexi is considered an oligotrophic bacteria in the soil and prefers to decompose some refractory organic compounds [37]. In the present study, the relative abundance of Chloroflexi decreased in treatments with fertilizer compared with CK. However, the abundance was relatively higher in treatments M and MNP, compared with treatment NP.
Most of dominant genera are widely distributed in the environment. For example, Sphingomonas has been found in both aqueous (both fresh- and seawater) and terrestrial habitats, plant root systems, and others. MND1 can be detected more frequently in soil, aquatic and subsurface ecosystems [40]. Arthrobacter is commonly found in soils, the aerial surfaces of plants, and wastewater sediments [41]. Some of the genera have the ability to use organic matter to grow and reproduce. For example, Sphingomonas distributes widely in the environment due to its ability to use a wide range of organic compounds and to grow and survive in low-nutrient conditions. Arthrobacter can degrade unusual and polymeric compounds and plays an important role in biodegrading agrochemicals and pollutants. In the present study, the relative abundance of these genera tended to be higher in treatments with manure than the others, which was consistent with these properties. In the previous study, Steroidobacter was reported to be a complex organic compound degrading bacteria [42], whose abundance was the highest in MNP treatment, followed by M treatment, and then by NP treatment in this study. On the contrary, the abundance of RB41 was the highest in CK treatment, consistent with Yan who observed the highest relative abundance of RB41 in abandoned farmland compared to other treatments with fertilizer [43].

4.3. Effect of Different Treatments on Soil Enzyme Activities

Enzymatic activity in soil is largely related to the microorganisms present, to available nutrients and to root exudates, as well [44,45]. Compared with CK, fertilization significantly increased the activities of SUC, URE, and ALP (p < 0.05). With the increase of nutrients input, soil enzyme activities gradually increased. Enzyme activities were the highest in treatment MNP. The results indicated that soil enzyme activity was directly affected by the amount of fertilizer applicated. Previous studies have shown that manure application was of paramount importance to improve soil enzyme activities, which significantly increased enzyme activities more than that of chemical fertilizer application [46].

4.4. Relationship between Soil Bacteria Community Composition and Environmental Factors and Soil Enzyme Activities

SUC, URE, and ALP were mainly involved in soil’s carbon, nitrogen, and phosphorus cycles, which were not only limited by soil nutrients, but also affected by other soil properties. In this study, AP and AN were two significant factors related to the changes of enzyme activities, which were mainly caused by fertilization regimes. Part of the soil enzymes are derived from the decomposed materials of bacterial cell exudates and residues, and in turn, soil bacterial diversity is greatly affected by soil enzyme activities [47]. Moreover, the composition of the bacterial community was influenced by its environment, including soil properties. Furthermore, through the RDA analysis, we found that AN, OM, URE, and ALP were main factors attributed to bacterial community structure in this region. This might be due to OM-provided substrates and energy for bacterial reproduction and metabolism [2]. Alkali-hydrolyzable nitrogen and carbon were coupled to participate in many bacterial metabolic processes. Previous studies also reported that microbial community composition had a significant correlation with URE [47].
In the present study, MNP treatment had the best effect on the improvement of bacterial diversity and enzyme activity, not only because of the maximum input of nitrogen, phosphorus, and potassium, but also because of the key advantages of combining of manure and chemical fertilizer application. Admittedly, manure contains a variety of nutrient elements, which can increase the input of exogenous carbon, promote the formation of aggregate structures, reduce soil bulk density, enhance the activity of microorganisms, promote the growth of crop roots and above-ground parts, and even increase the activity of soil enzymes [48,49,50]. However, a single application of manure alone is not conducive to the growth of crops, possibly because the slow and delayed release of nutrients in organic fertilizer cannot meet the needs of crops for nutrients in time [51]. Chemical fertilizers can provide quick nutrients, but studies have shown that a single application of chemical fertilizer can easily lead to soil degradation and is not suitable for the growth and reproduction of microorganisms. By contrast, manure has a good buffering ability to protect crops from adverse fluctuations and enhance the anti-interference ability of the soil [52,53]. Thus, manure and chemical fertilizer cooperation can foster strengths and circumvent weaknesses of either technique, and this can provide a more balanced supply of nutrients, just as it can also promote the growth of crops more effectively.

5. Conclusions

Various degrees of difference in bacterial community structures, soil enzyme activities and soil properties were observed with manure, chemical fertilizers, or both. As a whole, the combination of manure and chemical fertilizer showed great advantages, which could be a promising method in agricultural production. The proportion of the combination, however, has yet to be explored for different crops in the future.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/agronomy11051017/s1, Figure S1: Rarefraction Curve. M: manure; NP: chemical fertilizer; MNP: manure + chemical fertilizer; CK: no fertilizer. Figure S2: UPGMA clustering tree. M: manure; NP: chemical fertilizer; MNP: manure + chemical fertilizer; CK: no fertilizer. The diatance equals to 0.01. Table S1: F and p-value of soil properties between groups from the ANOVA. Table S2: Difference analysis of bacteria at phylum level. Table S3: Difference analysis of bacteria at genus level.

Author Contributions

Conceptualization: H.Z.; Writing-Original draft preparation: Z.L. and W.X.; Resources: Z.Y. Supervision: X.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Key R&D plan of Shanxi Province (201703D211002-3), Key R&D plan of Shanxi Province (201703D211002-4), Applied and Basic Research Program of Shanxi Province (201901D211557), Applied Basic Research Program of Shanxi Academy of Agricultural Sciences (YBSJJ2012).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original sequence data are available through the NCBI Sequence Read Archive (Accession: PRJNA706481).

Acknowledgments

We thank Biomarker Technologies Corporation for assistance in bioinformatics analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chaparro, J.M.; Sheflin, A.M.; Manter, D.K.; Vivanco, J.M. Manipulating the soil microbiome to increase soil health and plant fertility. Biol. Fert. Soils 2012, 48, 489–499. [Google Scholar] [CrossRef]
  2. Chen, H.; Liang, Q.; Gong, Y.; Kuzyakov, Y.; Fan, M.; Plante, A.F. Reduced tillage and increased residue retention increase enzyme activity and carbon and nitrogen concentrations in soil particle size fractions in a long-term field experiment on Loess Plateau in China. Soil Tillage Res. 2019, 194, 104–296. [Google Scholar] [CrossRef]
  3. Postma-Blaauw, M.B.; Goede, R.G.M.D.; Bloem, J.; Faber, J.H.; Brussaard, L. Soil biota community structure and abundance under agricultural intensification and extensification. Ecology 2010, 91, 460–473. [Google Scholar] [CrossRef]
  4. Ollivier, J.; Töwe, S.; Bannert, A.; Hai, B.; Kastl, E.-M.; Meyer, A.; Su, M.X.; Kleineidam, K.; Schloter, M. Nitrogen turnover in soil and global change. FEMS Microbiol. Ecol. 2011, 78, 3–16. [Google Scholar] [CrossRef] [Green Version]
  5. Burns, R.; Deforest, J.; Marxsen, J.; Sinsabaugh, R.; Stromberger, M.; Wallenstein, M.; Weintraub, M.; Zoppini, A. Soil enzymes in a changing environment: Current knowledge and future directions. Soil Biol. Biochem. 2013, 58, 216–234. [Google Scholar] [CrossRef]
  6. Zhao, S.; Li, K.; Zhou, W.; Qiu, S.; Huang, S.; He, P. Changes in soil microbial community, enzyme activities and organic matter fractions under long-term straw return in north-central China. Agr. Ecosyst. Environ. 2016, 216, 82–88. [Google Scholar] [CrossRef]
  7. Schimel, J.P.; Weintraub, M.N. The implications of exoenzyme activity on microbial carbon and nitrogen limitation in soil: A theoretical model. Soil Biol. Biochem. 2003, 35. [Google Scholar] [CrossRef]
  8. Hartmann, M.; Frey, B.; Mayer, J.; Mder, P.; Widmer, F. Distinct soil microbial diversity under long-term organic and conventional farming. ISME J. 2014, 9, 1177–1194. [Google Scholar] [CrossRef] [Green Version]
  9. Luo, P.; Han, X.; Wang, Y.; Han, M.; Shi, H.; Liu, N.; Bai, H. Influence of long-term fertilization on soil microbial biomass, dehydrogenase activity, and bacterial and fungal community structure in a brown soil of northeast China. Annu. Microbiol. 2015, 65, 533–542. [Google Scholar] [CrossRef] [Green Version]
  10. Yang, Y.; Wang, P.; Zeng, Z. Dynamics of Bacterial Communities in a 30-Year Fertilized Paddy Field under Different Organic–Inorganic Fertilization Strategies. Agronomy 2019, 9, 14. [Google Scholar] [CrossRef] [Green Version]
  11. Hahn, M.W.; Stadler, P.; Wu, Q.L.; Pöckl, M. The filtration-acclimatization method for isolation of an important fraction of the not readily cultivable bacteria. J. Microbiol. Meth. 2004, 57, 379–390. [Google Scholar] [CrossRef]
  12. Guan, L.; Cho, K.H.; Lee, J.H. Analysis of the cultivable bacterial community in jeotgal, a Korean salted andfermented seafood, and identification of its dominant bacteria. Food Microbiol. 2011, 28, 101–113. [Google Scholar] [CrossRef]
  13. Tebbe, C.C.; Dohrmann, A.B.; Hemkemeyer, M.; Nther, A. Microbial Community Profiling: SSCP and T-RFLP Techniques. Hydrocarbon and Lipid Microbiology Protocols, Springer Protocols Handbooks; Springer Nature: Basingstoke, UK, 2015. [Google Scholar]
  14. MacLean, D.; Jones, J.D.; Studholme, D.J. Application of next-generation sequencing technologies to microbial genetics. Nat. Rev. Microbiol. 2009, 7, 287–296. [Google Scholar] [CrossRef] [PubMed]
  15. ISS-CAS, Institute of Soil Science, Chinese Academy of Science Scientific Data Base. China Soil Data Base. Available online: http://www.soil.csdb.cn/ (accessed on 10 March 2021).
  16. IUSS Working Group WRB. World Reference Base for Soil Resources 2006. A Framework for International Classification, Correlation and Communication, 2nd ed.; World Soil Resources Reports No. 103; FAO: Rome, Italy, 2006. [Google Scholar]
  17. Bao, S.D. Analysis Method of Soil and Agricultural Chemistry, 3rd ed.; China Agricultural Press: Beijing, China, 2000; pp. 25–108. [Google Scholar]
  18. Guan, S.Y. Soil Enzyme and its Study Method; China AgriculturePress: Beijing, China, 1986; pp. 274–280. [Google Scholar]
  19. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Pena, A.G.; Goodrich, J.K.; Gordon, J.I.; et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 2010, 7, 335–336. [Google Scholar] [CrossRef] [Green Version]
  20. Magoc, T.; Salzberg, S.L. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 2011, 27, 2957–2963. [Google Scholar] [CrossRef]
  21. Edgar, R.C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010, 26, 2460–2461. [Google Scholar] [CrossRef] [Green Version]
  22. Quast, C.; Pruesse, E.; Yilmaz, P.; Gerken, J.; Schweer, T.; Yarza, P.; Peplies, J.; Gloeckner, F.O. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nuclceic Acids Res. 2012, 41, 590–596. [Google Scholar] [CrossRef]
  23. Varsos, C.; Patkos, T.; Oulas, A.; Pavloudi, C.; Gougousis, A.; Zeeshan Ijaz, U.; Filiopoulou, I.; Pattakos, N.; Vanden Berghe, E.; Fernández-Guerra, A.; et al. Optimized r functions for analysis of ecological community data using the r virtual laboratory (rvlab). Biodivers. Data J. 2016, 4, 8357–8383. [Google Scholar] [CrossRef] [Green Version]
  24. Jari Oksanen, F.; Guillaume, B.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. Package “vegan”: Community ecology package. Time Int. 2012, 1997, 15–17. [Google Scholar] [CrossRef]
  25. Ramette, A. Multivariate analyses in microbial ecology. FEMS Microbiol. Ecol. 2007, 62, 142–160. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. McArdle, B.H.; Anderson, M.J. Fitting multivariate models to community data: A comment on distance-based redundancy analysis. Ecology 2001, 82, 290–297. [Google Scholar] [CrossRef]
  27. Huson, D.H.; Mitra, S.; Ruscheweyh, H.J.; Weber, N.; Schuster, S.C. Integrative analysis of environmental sequences using MEGAN4. Genome Res. 2011, 21, 1552–1560. [Google Scholar] [CrossRef] [Green Version]
  28. Zaura, E.; Keijser, B.J.F.; Huse, S.M.; Crielaard, W. Defining the healthy “core microbiome” of oral microbial communities. BMC Microbiol. 2009, 9, 1–12. [Google Scholar] [CrossRef] [Green Version]
  29. Wang, W.; Luo, X.; Chen, Y.; Ye, X.; Wang, H.; Cao, Z.; Ran, W.; Cui, Z. Succession of Composition and Function of Soil Bacterial Communities During Key Rice Growth Stages. Front. Microbiol. 2019, 10, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Tian, D.; Niu, S. A global analysis of soil acidification caused by nitrogen addition. Environ. Res. Lett. 2015, 10, 1714–1721. [Google Scholar] [CrossRef]
  31. Schroder, J.L.; Zhang, H.; Girma, K.; Raun, W.R.; Penn, C.J.; Payton, M.E. Soil acidification from long-term use of nitrogen fertilizers on winter wheat. Soil Sci. Soc. Am. J. 2011, 75, 957–964. [Google Scholar] [CrossRef]
  32. Bolan, N.; Hedley, M.; White, R. Processes of soil acidification during nitrogen cycling with emphasis on legume based pastures. Plant Soil 1991, 134, 53–63. [Google Scholar] [CrossRef]
  33. Powlson, D.S.; Brookes, P.C.; Christensen, B.T. Measurement of soil microbial biomass provides an early indication of changes in total soil organic matter due to straw incorporation. Soil Biol. Biochem. 1987, 19, 159–164. [Google Scholar] [CrossRef]
  34. Liang, R.; Hou, R.; Li, J.; Lyu, Y.; Ouyang, Z. Effects of different fertilizers on rhizosphere bacterial communities of winter wheat in the north china plain. Agronomy 2020, 10, 93. [Google Scholar] [CrossRef] [Green Version]
  35. Lauber, C.L.; Hamady, M.; Knight, R.; Fierer, N. Pyrosequencing-based assessment of soil p H as a predictor of soil bacterial community structure at the continental scale. Appl. Environ. Microbiol. 2009, 75, 5111–5120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Sun, R.; Zhang, X.X.; Guo, X.; Wang, D.; Chu, H. 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. 2015, 88, 9–18. [Google Scholar] [CrossRef]
  37. Fierer, N.; Bradford, M.A.; Jackson, R.B. Toward an ecological classification of soil bacteria. Ecology 2007, 88, 1354–1364. [Google Scholar] [CrossRef]
  38. Fierer, N.; Lauber, C.L.; Ramirez, K.S.; Zaneveld, J.; Bradford, M.A.; Knight, R. Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients. ISME J. 2012, 6, 1007–1017. [Google Scholar] [CrossRef] [Green Version]
  39. DeBruyn, J.M.; Nixon, L.T.; Fawaz, M.N.; Johnson, A.M.; Radosevich, M. Global biogeography and quantitative seasonal dynamics of Gemmatimonadetes in soil. Appl. Environ. Microbiol. 2011, 77, 6295–6300. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Spain, A.M.; Krumholz, L.R.; Elshahed, M.S. Abundance, composition, diversity and novelty of soil Proteobacteria. ISME J. 2009, 3, 992–1000. [Google Scholar] [CrossRef] [PubMed]
  41. Funke, G.; Pagano-Niederer, M.; Berit, S.; Falsen, E. Characteristics of Arthrobacter cumminsii, the most frequently encountered Arthrobacter species in human clinical specimens. J. Clin. Microbiol. 1998, 36, 1539–1543. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Alcaraz, L.D.; Peimbert, M.; Barajas, H.R.; Dorantes-Acosta, A.E.; Bowman, J.L.; Arteaga-Vázquez, M.A. Marchantia liverworts as a proxy to plants’ basal microbiomes. Sci. Rep. 2018, 8, 12712–12723. [Google Scholar] [CrossRef]
  43. Yan, S.; Song, J.; Fan, J.; Yan, C.; Dong, S.; Ma, C.; Gong, Z. Changes in soil organic carbon fractions and microbial community under rice straw return in Northeast China. Glob. Ecol. Conserv. 2020, 22, 962–973. [Google Scholar] [CrossRef]
  44. Moreno-Espíndola, I.P.; Ferrara-Guerrero, M.J.; Luna-Guido, M.L.; Ramírez-Villanueva, D.A.; de León-Lorenzana, A.S.; Gómez-Acata, S.; González-Terreros, E.; Ramírez-Barajas, B.; Navarro-Noya, Y.E.; Sánchez-Rodríguez, L.M.; et al. The Bacterial community structure and microbial activity in a traditional organic milpa farming system under different soil moisture conditions. Front. Microbiol. 2018, 9, 1–19. [Google Scholar] [CrossRef]
  45. Caldwell, B.A. Enzyme activities as a component of soil biodiversity: A review. Pedobiologia 2005, 49, 637–644. [Google Scholar] [CrossRef]
  46. Amadou, A.; Alin Song, A.; Tang, Z.; Li, Y.; Wang, E.; Lu, Y.; Liu, X.; Yi, K.; Zhang, B.; Fan, F. The effects of organic and mineral fertilization on soil enzyme activities and bacterial community in the below- and above-ground parts of wheat. Agronomy 2020, 10, 1452. [Google Scholar] [CrossRef]
  47. Xu, L.; Yi, M.; Yi, H.; Guo, E.; Zhang, A. Manure and mineral fertilization change enzyme activity and bacterial community in millet rhizosphere soils. World J. Microb. Biot. 2018, 34, 1–13. [Google Scholar] [CrossRef]
  48. Ros, M.; Pascual, J.A.; Garcia, C.; Hernandez, M.T.; Insam, H. Hydrolase activities, microbial biomass and bacterial community in a soil after long-term amendment with different composts. Soil Biol. Biochem. 2006, 38, 3443–3452. [Google Scholar] [CrossRef]
  49. Dorado, J.; Zancada, M.C.; Almendros, G.; López-Fando, C. Changes in soil properties and humic substances after long-term amendments with manure and crop residues in dryland farming system. J. Plant Nutr. Soil Sci. 2003, 166, 31–38. [Google Scholar] [CrossRef]
  50. Bastida, F.; Kandeler, E.; Hernández, T.; García, C. Long-term effect of municipal solid waste amendment on microbial abundance and humus-associated enzyme activities under semiarid conditions. Microb. Ecol. 2008, 55, 651–661. [Google Scholar] [CrossRef]
  51. Ginting, D.; Kessavalou, A.; Eghball, B.; Doran, J.W. Greenhouse gas emissions and soil indicators four years after manure and compost applications. J. Environ. Qual. 2003, 32, 23–32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Diacono, M.; Montemurro, F. Long-term effects of organic amendments on soil fertility. A review. Agron. Sustain. Dev. 2010, 30, 401–422. [Google Scholar] [CrossRef] [Green Version]
  53. Watson, C.A.; Atkinson, D.; Gosling, P.; Jackson, L.R.; Rayns, F.W. Managing soil fertility in organic farming systems. Soil Use Manag. 2002, 18, 239–247. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Venn Diagram showing unique and overlapped OTUs between different fertilization treatments. M: manure; NP: NP chemical fertilizer; MNP: manure + NP chemical fertilizer; CK: no fertilizer.
Figure 1. Venn Diagram showing unique and overlapped OTUs between different fertilization treatments. M: manure; NP: NP chemical fertilizer; MNP: manure + NP chemical fertilizer; CK: no fertilizer.
Agronomy 11 01017 g001
Figure 2. PCoA of the bacterial community compositions in soil under different treatments based on OTUs. M: manure; NP: chemical fertilizer; MNP: manure + chemical fertilizer; CK: no fertilizer.
Figure 2. PCoA of the bacterial community compositions in soil under different treatments based on OTUs. M: manure; NP: chemical fertilizer; MNP: manure + chemical fertilizer; CK: no fertilizer.
Agronomy 11 01017 g002
Figure 3. PERMANOVA analysis under different treatments. M: manure; NP: chemical fertilizer; MNP: manure + chemical fertilizer; CK: no fertilizer.
Figure 3. PERMANOVA analysis under different treatments. M: manure; NP: chemical fertilizer; MNP: manure + chemical fertilizer; CK: no fertilizer.
Agronomy 11 01017 g003
Figure 4. Taxonomic composition and abundance distribution of bacteria at the phylum level under different treatments. M: manure; NP: chemical fertilizer; MNP: manure + chemical fertilizer; CK: no fertilizer.
Figure 4. Taxonomic composition and abundance distribution of bacteria at the phylum level under different treatments. M: manure; NP: chemical fertilizer; MNP: manure + chemical fertilizer; CK: no fertilizer.
Agronomy 11 01017 g004
Figure 5. Taxonomic composition and abundance distribution of bacteria at the genus level under different treatments. M: manure; NP: chemical fertilizer; MNP: manure + chemical fertilizer; CK: no fertilizer.
Figure 5. Taxonomic composition and abundance distribution of bacteria at the genus level under different treatments. M: manure; NP: chemical fertilizer; MNP: manure + chemical fertilizer; CK: no fertilizer.
Agronomy 11 01017 g005
Figure 6. RDA depicting the relationship between the soil properties, the phylum-level soil bacteria and soil enzyme activities. (a) The relationship between soil bacterial community composition and properties. (b) The relationship between soil enzyme activities and properties. (c) The relationship between soil bacterial community composition and enzyme activities.
Figure 6. RDA depicting the relationship between the soil properties, the phylum-level soil bacteria and soil enzyme activities. (a) The relationship between soil bacterial community composition and properties. (b) The relationship between soil enzyme activities and properties. (c) The relationship between soil bacterial community composition and enzyme activities.
Agronomy 11 01017 g006
Table 1. Experimental treatments and fertilization.
Table 1. Experimental treatments and fertilization.
TreatmentsChemical Fertilizer (kg/hm2)Manure (kg/hm2)Total Nutrient (kg/hm2)
NP2O5K2ONP2O5K2ONP2O5K2O
Manure (M)00080551808055180
Chemical fertilizer (NP)120750000120750
Manure + Chemical fertilizer (MNP)1207508055180200130180
No fertilizer (CK)000000000
Table 2. Soil physicochemical properties under different treatments.
Table 2. Soil physicochemical properties under different treatments.
Treatments pHAN
(mg/kg)
AP
(mg/kg)
AK
(mg/kg)
OM
(g/kg)
TN
(g/kg)
TP
(g/kg)
TK
(g/kg)
M8.40 ± 0.23 a63.09 ± 3.06 ab22.53 ± 1.73 a244.78 ± 16.67 a11.76 ± 2.26 a0.81 ± 0.10 a0.80 ± 0.03 b12.79 ± 0.65 a
NP8.51 ± 0.13a59.15 ± 7.35 b14.40 ± 5.52 b71.27 ± 4.25 c7.82 ± 0.57 b0.74 ± 0.14 ab0.84 ± 0.04 ab12.58 ± 0.34 a
MNP8.40 ± 0.18 a74.55 ± 4.55 a26.90 ± 2.80 a223.90 ± 3.04 b12.27 ± 1.14 a0.85 ± 0.12 a0.85 ± 0.12 a12.36 ± 0.74 a
CK8.59 ± 0.10 a37.98 ± 1.23 c5.10 ± 0.70 c60.01 ± 8.30 c6.52 ± 1.48 b0.56 ± 0.06 b0.66 ± 0.03 c12.78 ± 0.17 a
Note: M: manure; NP: NP chemical fertilizer; MNP: manure + NP chemical fertilizer; CK: no fertilizer. AN: alkali-hydrolyzable Nitrogen; AP: available phosphorus; AK: available potassium; OM: organic matter; TN: total nitrogen; TP: total phosphorus; TK: total potassium; Different letters in the same column indicate significant differences (ANOVA followed by LSD post hoc test, n = 3, p < 0.05, average value, SD standard deviation).
Table 3. α-diversity indices of the bacteria under different treatments.
Table 3. α-diversity indices of the bacteria under different treatments.
TreatmentsACEShannon
M1828.12 ± 6.06 b6.57 ± 0.01 a
NP1769.66 ± 7.89 c6.53 ± 0.02 ab
MNP1867.82 ± 8.81 a6.55 ± 0.05 ab
CK1665.51 ± 5.54 d6.48 ± 0.02 b
Note: M: manure; NP: NP chemical fertilizer; MNP: manure + NP chemical fertilizer; CK: no fertilizer. Different letters in the same column indicate significant differences (ANOVA followed by LSD post hoc test, n = 3, p < 0.05, average value, SD standard deviation).
Table 4. Effects of different treatments on soil enzyme activities.
Table 4. Effects of different treatments on soil enzyme activities.
TreatmentsSUC Activity
(mg·Glucose·g−1·24 h−1)
URE Activity
(mg·NH3-N·g−1·24 h−1)
ALP Activity
(mg Phenol·g−1·24 h−1)
M19.55 ± 3.227 ab1.983 ± 0.021 b1.820 ± 0.080 b
NP15.78 ± 1.520 b1.900 ± 0.290 b1.404 ± 0.153 bc
MNP21.42 ± 0.708 a2.537 ± 0.218 a2.264 ± 0.278 a
CK8.723 ± 0.083 c1.210 ± 0.128 c1.348 ± 0.016 c
Note: M: manure; NP: NP chemical fertilizer; MNP: manure + NP chemical fertilizer; CK: no fertilizer. Different letters in the same column indicate significant differences (ANOVA followed by LSD post hoc test, n = 3, p < 0.05, average value, SD standard deviation).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, Z.; Xie, W.; Yang, Z.; Huang, X.; Zhou, H. Effects of Manure and Chemical Fertilizer on Bacterial Community Structure and Soil Enzyme Activities in North China. Agronomy 2021, 11, 1017. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11051017

AMA Style

Liu Z, Xie W, Yang Z, Huang X, Zhou H. Effects of Manure and Chemical Fertilizer on Bacterial Community Structure and Soil Enzyme Activities in North China. Agronomy. 2021; 11(5):1017. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11051017

Chicago/Turabian Style

Liu, Zhiping, Wenyan Xie, Zhenxing Yang, Xuefang Huang, and Huaiping Zhou. 2021. "Effects of Manure and Chemical Fertilizer on Bacterial Community Structure and Soil Enzyme Activities in North China" Agronomy 11, no. 5: 1017. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11051017

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop