Abstract
Cover crops can improve soil biological health and alter the composition of soil microbial communities in agricultural systems. However, the effects of diversified cover crops on soil microbial communities in continuous cropping systems are unclear. Here, using different soil biochemical analysis, quantitative PCR and 16S rRNA amplicon sequencing, we investigated the effects of cover crops, alone or in mixture, on soil physicochemical properties in 2019 and 2020, and soil bacterial communities in 2020 in a continuous pepper cropping system. A field trial was established before pepper planting and eight treatments were included: fallow (no cover crop; CK); three sole cover crop treatments: wheat (Triticum aestivum L.; W), faba bean (Vicia faba L.; B), and wild rocket (Diplotaxis tenuifolia; R); and four mixed treatments: wheat + wild rocket (WR), wheat + faba bean (WB), wild rocket + faba bean (RB), and wheat + wild rocket + faba bean (WRB). The pepper yield was increased in the WR and WB in 2019 and 2020, and in the WRB in 2020. Cover crops increased the soil pH, but decreased the concentrations of NH4+ and available phosphorus. Bacterial abundance was increased by cover crop treatments, and community structure was altered in the W, WB, and WRB treatments. Moreover, we found that pH was the key factor associated with the changes in the abundance and structure of the bacterial community. Cover crop treatments altered the bacterial community structure with shifts in the dominant genera, which have plant-growth-promoting and/or pathogen-antagonistic potentials, e.g., increased the abundances of Streptomyces, Arthrobacter, and Bacillus spp. in the W and WRB, and Gaiella spp. in the WB. Overall, we found that cover crops altered the soil physicochemical properties and bacterial community, and these changes varied with species composition of the cover crops, with wheat and its combination with legumes as most effective treatments. These results suggest that the diversification within cover crops could provide better crop yield stimulatory affects by altering soil biochemical environment.
Introduction
Modern agricultural crop production is usually focused on crop monocultures or minimum rotations, which usually negatively affect soil quality (e.g., imbalances in soil nutrition, decrease in soil microbial diversity, and the accumulation of autotoxic compounds) as well as crop yield (Gao et al., 2017). This is a negative plant-soil feedback phenomenon that has been described as “soil sickness.” Therefore, integrated soil management practices, e.g., increase in plant diversity in the field such as crop rotation, intercropping and cultivation of cover crop, have been used to alleviate soil sickness (Huang et al., 2013).
Cover crops are defined as the crops, which are grown primarily to cover and protect the soils from erosion and nutrient losses in between the period of main crop cultivation (Sarrantonio and Gallandt, 2003). Studies have shown that cover crops can be used as living mulch, green manure, and biofumigant (Snapp et al., 2005; Schipanski et al., 2014). Legumes, brassicas, cereals, and other broadleaf crops can be used as cover crops (Snapp et al., 2005; Schipanski et al., 2014). The cultivation of legumes as a cover crop is effective to improve soil N availability and moisture retention (Sullivan et al., 1991). Brassicaceae species, as cover crops, have been shown to control various soilborne diseases, and improve soil characteristics and crop yield (Larkin et al., 2010). Cereals as cover crops provide many benefits, for example, these crops increase soil fertility and organic carbon (Rorick and Kladivko, 2016), improve soil water retention (Villamil et al., 2006), and decrease bulk density (Moore et al., 2014). Overall, the use of cover crops during the fallow period in the field has the potential to improve the soil environment and increase crop productivity (Guo et al., 2008; Tian et al., 2011).
Soil microorganisms are critically important for maintaining soil functions and ecosystem sustainability because they are involved in the cycling of nutrients and the turnover of organic matter (Liu et al., 2020a). Studies have shown that cover crops can alter the dynamics of soil bacterial and fungal communities (Patkowska et al., 2016), stimulate beneficial microorganisms (Berendsen et al., 2012), and suppress soilborne pathogens (Manici et al., 2004). Cover crops vary in their characteristics and can have diverse effects on soil microbial community composition and functions (Thapa et al., 2021). For example, oats with dense rooting systems as cover crops can increase the abundance of arbuscular mycorrhizal fungi compared with the fallow field (Kabir and Koide, 2002). Legumes with high nitrogen (N) content can promote the proliferation of microbes as compared to non-leguminous cover crops (Murungu et al., 2011). Soil substrate availability and quality are related to plant diversity and can drive changes in soil microbial community structure (Thapa et al., 2021). Mixtures of cover crop species can vary in their effects on substrates, alter environmental conditions, and increase the diversity of microbial communities (Spedding et al., 2004). For example, grass cover crops mixed with legumes promote microbial activity because of their complementary effects on residue decomposition (Njeru et al., 2014). Diverse mixtures of cover crops increase the size of total microbial community and fungal community as compared to fallow (Thapa et al., 2021). However, the effect of multispecies cover crops on microbial community in continuous cropping systems remain poorly understood.
In this study, we evaluated the effects of three different cover crops (single as well as in mixtures) on main crop yield, soil physicochemical properties, and bacterial communities in the soils of continuous pepper cropping systems. Soil bacterial community composition and abundance were analyzed by 16S rRNA amplicon sequencing and quantitative PCR, respectively. The important carbon resources for soil microorganisms from plant root exudates and litter decomposition, and the chemical composition of these substances differ among plant species (Fierer et al., 2007; Eilers et al., 2010; Birouste et al., 2012). Therefore, we hypothesized that increased cover crop diversity would lead to higher soil bacterial community diversity and abundances. Moreover, since the soil bacterial communities are sensitive to soil environmental factors (Kuramae et al., 2011; Shen et al., 2013), we also hypothesized that cover crops may affect soil bacterial communities by altering the physicochemical properties of the soil.
Materials and Methods
Study Site and Treatments
The experiment was conducted from March 2019 to October 2020 at the Horticultural Experimental Station of Northeast Agricultural University (45°41′N, 126°37′E) Harbin, Heilongjiang Province, China. The soil was black soil (Mollisoil) and has been under continuous pepper cropping system for 7 years. The soil had the following features: organic matter 32.76 g kg−1, available phosphorus (AP) 105.90 mg kg−1, available potassium (AK) 222.21 mg kg−1, NO3−-N 131.88 mg kg−1, NH4+-N 16.72 mg kg−1, pH (1:2.5, w/v) 7.29, and EC (1:2.5, w/v) 0.64 mS cm−2.
The experiment was conducted using a pepper (Capsicum annuum L.)/cover crop rotation with a randomized block design consisting of eight treatments. There were three replicate plots for each treatment, and the size of each plot was 4.5 m2 (5 m × 0.9 m). Cover crops treatments consist on three sole cover crops treatments: wheat (Triticum aestivum L; W), faba bean (Vicia faba L.; B), and wild rocket (Diplotaxis tenuifolia; R); and four mixed cover crops treatments: wheat + wild rocket (WR), wheat + faba bean (WB), wild rocket + faba bean (RB), and wheat + wild rocket + faba bean (WRB), and a fallow treatment (no cover crops; CK).
Pepper and Cover Crop Management
Cover crops were planted in fallow periods. In cover crop plots, the first phase of rotation (cover crop-pepper) began with the planting of cover crops in March 2019 and ended in April 2019. Pepper was planted in May 2019 and harvested in August 2019. The second phase of rotation (crop cover-pepper) started in October 2019 and ended in April 2020, followed by pepper planting in May 2020 and harvested in October 2020 (Figure 1).
Figure 1
For planting sole cover crops, 34 g of wheat, 420 g of faba bean, and 3.96 g of wild rocket seeds were sown per plot. In two or three cover crop mixture treatments, the seeding rate decreased to 50 or 33% of the monoculture rate for each species, respectively. One month later, samples of aboveground biomass were hand clipped from a 1-m2 area (0.20-m2 quadrat area at five random spots) in each plot, and dry biomass was measured after being subjected to 65°C for 72 h (Supplementary Table S1). The rest of the harvested plant residues (above-ground) were naturally dried, cut into 10-cm pieces, and incorporated into the decomposed soil. Thirty days later, the pepper seedlings at the six-leaf stage were transplanted to the greenhouse. There were two rows in each plot, and 14 pepper plants were there in one row.
Soil Sampling
One hundred and twenty days after pepper transplantation in 2019 and 2020, bulk and rhizosphere soil samples were collected for soil physicochemical properties analysis and DNA extraction, respectively.
Bulk Soil Sampling and Soil Physicochemical Properties Analysis
Five soil cores (8 cm diameter, 20 cm deep) were randomly collected from the inner two rows of pepper in each plot and mixed to make a composite sample – one replicate. There were three replicates of bulk soil samples for each treatment. The soils were sieved through a 2-mm mesh to remove the stones, roots, and plant debris. The soil samples were stored at −4°C for soil physicochemical analysis.
Soil EC and pH were measured in a soil-water suspension at a water/soil ratio of 2.5:1 with a conductivity meter and glass electrode, respectively. Soil ammonia nitrogen (NH4+-N) and nitrate (NO3−-N), available P (AP) and K (AK) were extracted with 2 M KCl, 1 M NaHCO3, and 0.5 M CH3COONH4, respectively. Then, measurements were taken from the soil extract using a Continuous Flow Analyzer (SAN++, Skalar, Breda, Netherlands; Gao et al., 2017). The concentrations of soil total nitrogen (TN) and soil organic carbon (SOC) were determined by the micro-Kjeldahl method and wet oxidation redox titration method, respectively (Lin et al., 2018).
Rhizosphere Soil Sampling
The pepper roots were removed from the soil carefully. The soil loosely attached to pepper roots was removed by manual shaking. The soil tightly adhering to the roots was removed from the root surface by a sterile brush, sieved through a 2 mm mesh, and considered as rhizosphere soil. The rhizosphere soil samples from five plants in each treatment replicate were mixed to make a composite sample. The samples were stored at −80°C for DNA extraction; one part from 2019 and 2020 sampling each was used to determine the abundances of soil total bacteria, total fungi, Bacillus, and Pseudomonas spp., while other part was used for Illumina Miseq sequencing analysis for soil bacterial community composition.
DNA Extraction and Quantitative PCR
Soil total DNA was extracted from rhizosphere soil using a PowerSoil DNA Isolation Kit (MO BIO Laboratories, CA, United States) according to the manufacturer’s protocols. Agarose gel electrophoretic analysis and a NanoDrop 2000 Spectrophotometer (Thermo Scientific, United States) were used to assess the quantity and quality of DNA. The abundances of soil total bacteria, fungi, Bacillus, and Pseudomonas spp. in pepper rhizosphere were estimated by quantitative PCR assays conducted with an IQ5 real-time PCR system (Bio-Rad Lab, LA, United States). For total bacteria and fungi, the primer sets 338F/518R (Jin et al., 2019) and ITS1F/ITS4 (Bellemain et al., 2010) were used to amplify the V3-V4 hypervariable region of the bacterial 16S rRNA gene and ITS1 region of the fungal ITS gene, respectively. The primer sets PsF/PsR (Garbeva et al., 2004) and BacF/BacR (Garbeva et al., 2003) were used for Pseudomonas and Bacillus spp., respectively. Melting curves and agarose gel electrophoresis were used to verify the specificity of the amplicons. Standard curves were created using 10-fold dilution series of plasmid DNA with inserted target genes. The fluorescence signal was used to calculate Ct (cycle threshold) values, using the thermocycler software. The copy number of genes in the template DNA was calculated according to the linear equation and the Ct value conversion (Henry et al., 2006).
Illumina Miseq Sequencing and Data Processing
The compositions of soil bacterial communities in pepper rhizosphere in 2020 were analyzed with Illumina MiSeq. The V3-V4 hypervariable region of the bacterial 16S rRNA gene was amplified using the primers F338 (5′-ACTCCTACGGGAGGCAGCA-3′) and R806 (5′-ATGCAGGGACTACHVGGGTWTCTAAT-3′; Derakhshani et al., 2016). A 6-bp unique barcode was added to both the forward and reverse primers for each sample. Triplicate PCR amplifications were pooled and PCR products were purified with an Agarose Gel DNA Purification Kit (TaKaRa, China). The purified amplicons were then sequenced (2 × 300) on an Illumina Miseq platform at Majorbio Bio-Pharm Technology Co., Ltd., Shanghai, China. A total of 1,239,149 sequences were obtained from the 24 samples. The raw sequence reads were processed using the QIIME pipeline (Caporaso et al., 2010). Briefly, barcodes, adaptor sequences, and 30 low-quality bases at the end of each read were removed. Paired reads were joined (minimum overlapping read length of 20 base pairs) and quality filtered (Phred score of 20), and reads lower than 200 base pairs were removed. Chimaeras were screened and removed using USEARCH with the UCHIME algorithm (Edgar, 2013). The remaining high-quality sequences (97% similarity) were assigned to the same operational taxonomic units (OTUs) with UPARSE (Edgar, 2013). Taxonomies were assigned to the representative sequence of each OTU using the SILVA database release 128 (Quast et al., 2013). OTUs classified as chloroplasts and mitochondria, and singleton OTUs were eliminated. To avoid potential bias caused by the sequencing depth, all samples were randomly subsampled based on the lowest minimum number of sequence (20,008) per sample. The data set has been uploaded to the NCBI Sequence Read Archive (Accession Number SRP326928).
Statistical Analysis
Data of pepper yield, soil physicochemical properties, alpha diversity indices, and microbial abundances were analyzed by one-way ANOVA. The means of different treatments were compared with Tukey’s honestly significant difference (HSD) test at the 0.05 probability level. Relative weight analysis was performed using the “relweights” package in R (Kabacoff, 2015). A multivariate regression tree (MRT) was built to identify the most important soil environmental factors affecting bacterial community diversity using the “mvpart” package in R (Yang et al., 2019). Spearman correlation coefficients were calculated to test the relationships between pepper yield or the relative abundance of genera and soil physicochemical properties in SPSS software (Version 17.0). The Shannon index was calculated by QIIME (Caporaso et al., 2010). Principal coordinate analysis (PCoA) was used to analyze the composition of the bacterial community based on the Bray–Curtis distance dissimilarity. A redundancy analysis (RDA) was conducted to evaluate the relationships between the soil environmental factors and soil bacterial communities, and Mantel tests were used to determine the effects of environmental factors on the composition of the bacterial community using the “vegan” package in R. Species indicator analysis was conducted to identify the specific OTUs in different treatments using the “indicspecies” package in R (Cáceres and Legendre, 2009), and the results were displayed using iTOL tool (Letunic and Bork, 2007).
Results
Pepper Yield
Pepper yield was much higher in the WR and WB treatments than in CK in 2019, and no difference was observed among the other treatments. In 2020, pepper yield in WR, WB and WRB treatments were significantly higher than that in rest of the treatments, and significant differences were observed between the R and B treatments. The pepper yield in 2020 was higher than that in 2019 (Figure 2).
Figure 2
Soil Physicochemical Properties
The soil physicochemical properties of the treatments are shown in Table 1. The lowest soil pH was observed in CK. However, the EC and the concentrations of NH4+-N and AP were decreased in the cover crop treatments compared with CK. The concentration of NO3−-N was higher in W, WR, RB, and WRB treatments than in CK in both years. The concentration of AK was higher in the W and WR treatments in 2019 and the W and B treatments in 2020 than in CK, no differences were observed among the other treatments. The content of TN was higher in the cover crop treatments than in CK in both years, with the exception of the RB treatment in both years and the B treatment in 2019. The SOM was higher in the B, WR, and WB treatments in 2019 and the WB treatment in 2020 than in CK; SOM was highest in the WB treatment and lowest in the WRB treatment. C/N was higher in CK compared with the other treatments, with the exception of the B and WB treatments in 2019. In 2020, the C/N was higher in CK than in the WRB treatment (Table 1).
Table 1
| Treatment | pH | EC (uS.cm−1) | NO3− (mg kg−1) | NH4+ (mg kg−1) | TN (g kg−1) | AP (mg kg−1) | AK (mg kg−1) | SOM | C/N | |
|---|---|---|---|---|---|---|---|---|---|---|
| 2019 | CK | 6.96 ± 0.01 e | 408.33 ± 20.60 a | 36.68 ± 1.82 d | 25.24 ± 2.52 a | 3.36 ± 0.14 d | 81.44 ± 7.18 a | 201.67 ± 3.51 cd | 61.14 ± 2.60 cd | 18.17 ± 0.10 b |
| W | 7.06 ± 0.01 bc | 301.33 ± 22.50 b | 49.88 ± 2.18 abc | 20.32 ± 1.42 b | 3.76 ± 0.06 ab | 42.10 ± 2.23 b | 227.33 ± 3.21 a | 61.87 ± 1.03 bcd | 16.45 ± 0.04 e | |
| R | 7.06 ± 0.01 bc | 281.33 ± 18.45 b | 45.25 ± 2.29 c | 20.79 ± 1.42 b | 3.63 ± 0.04 bc | 47.63 ± 2.10 b | 211.67 ± 4.62 bc | 63.82 ± 0.69 abc | 17.58 ± 0.02 d | |
| B | 7.03 ± 0.01 d | 256.33 ± 11.68 bc | 44.14 ± 1.85 c | 19.68 ± 0.72 bc | 3.50 ± 0.02 cd | 46.28 ± 2.21 b | 201.67 ± 3.79 cd | 65.03 ± 0.34 ab | 18.60 ± 0.02 a | |
| WR | 7.08 ± 0.01 ab | 258.00 ± 10.54 bc | 53.28 ± 3.06 a | 17.95 ± 0.84 bc | 3.63 ± 0.02 bc | 46.33 ± 0.82 b | 216.00 ± 5.20 b | 65.28 ± 0.34 ab | 17.98 ± 0.02 c | |
| WB | 7.08 ± 0.01 ab | 269.33 ± 14.01 b | 30.62 ± 0.37 d | 18.37 ± 0.86 bc | 3.63 ± 0.02 bc | 43.58 ± 2.02 b | 202.00 ± 1.00 c | 67.47 ± 0.34 a | 18.59 ± 0.02 a | |
| RB | 7.05 ± 0.01 cd | 277.33 ± 14.05 b | 46.32 ± 2.42 bc | 16.60 ± 0.71 c | 3.43 ± 0.03 d | 46.92 ± 2.97 b | 198.00 ± 2.65 c | 62.60 ± 0.60 bcd | 18.25 ± 0.02 b | |
| WRB | 7.09 ± 0.01 a | 222.33 ± 14.19 c | 51.51 ± 2.37 ab | 19.39 ± 0.87 bc | 3.83 ± 0.05 a | 45.73 ± 0.23 b | 198.00 ± 2.66 c | 59.92 ± 0.69 d | 15.66 ± 0.02 f | |
| 2020 | CK | 6.91 ± 0.01 e | 448.67 ± 17.24 a | 36.84 ± 0.18 d | 21.04 ± 0.82 a | 3.27 ± 0.02 c | 59.45 ± 1.98 a | 221.67 ± 7.09 c | 60.30 ± 2.17 bc | 18.46 ± 0.65 ab |
| W | 7.08 ± 0.01 bc | 247.33 ± 15.31 b | 51.97 ± 4.65 a | 18.94 ± 0.60 b | 3.73 ± 0.04 a | 49.74 ± 0.46 bc | 262.33 ± 4.51 a | 61.78 ± 0.60 abc | 16.55 ± 0.16 bc | |
| R | 7.08 ± 0.01 bc | 230.67 ± 12.86 b | 43.90 ± 2.07 bcd | 18.75 ± 0.76 b | 3.54 ± 0.15 ab | 49.29 ± 0.86 c | 236.67 ± 9.61 bc | 63.50 ± 1.84 abc | 17.97 ± 0.52 ab | |
| B | 7.05 ± 0.01 d | 172.67 ± 14.01 d | 43.47 ± 1.48 cd | 17.91 ± 0.21 bc | 3.69 ± 0.07 ab | 45.65 ± 2.01 cd | 248.67 ± 7.02 ab | 66.70 ± 1.25 ab | 18.05 ± 0.34 ab | |
| WR | 7.12 ± 0.01 a | 174.67 ± 11.15 cd | 50.92 ± 2.06 ab | 18.05 ± 0.28 bc | 3.60 ± 0.01 ab | 48.64 ± 1.33 c | 224.33 ± 9.71 c | 64.24 ± 3.32 abc | 17.84 ± 0.92 ab | |
| WB | 7.11 ± 0.01 ab | 229.67 ± 17.56 b | 29.17 ± 1.28 e | 17.62 ± 0.57 bc | 3.60 ± 0.20 ab | 43.58 ± 2.02 d | 237.33 ± 7.51 bc | 67.43 ± 1.94 a | 18.74 ± 0.48 a | |
| RB | 7.07 ± 0.01 cd | 218.33 ± 17.79 bc | 45.65 ± 1.99 abc | 15.78 ± 0.57 d | 3.47 ± 0.16 bc | 49.14 ± 1.09 c | 221.33 ± 8.08 c | 62.02 ± 2.85 abc | 17.89 ± 0.82 ab | |
| WRB | 7.12 ± 0.01 a | 138.00 ± 18.73 d | 50.10 ± 4.03 abc | 16.48 ± 0.39 cd | 3.73 ± 0.04 a | 53.58 ± 1.22 b | 240.67 ± 7.64 abc | 59.07 ± 0.35 c | 15.82 ± 0.09 c |
Soil physicochemical properties in different treatments in 2019 and 2020.
Different letters in each column indicate significant differences among the samples (p < 0.05, HSD Tukey test). TN, soil total nitrogen; AP, available phosphorus; AK, available potassium; and SOM, soil organic matter. CK, no cover crop; W, wheat; R, wild rocket; B, faba bean; WR, wheat + wild rocket; WB, wheat + faba bean; RB, wild rocket + faba bean; and WRB, wheat + wild rocket + faba bean.
Abundances of Bacteria, Fungi, Pseudomonas, and Bacillus spp.
The abundances of bacteria were significantly higher in the cover crop treatments than in CK in 2019 and 2020, with the exception of the W treatment in 2019, and the abundance of bacteria was the highest in the WRB treatment (Figure 3A). The abundance of fungi was higher in all cover crop treatments in 2020 than in CK, while no difference among treatment was found in 2019 (Figure 3B). Cover crop treatments significantly increased the abundance of Pseudomonas in 2019 and 2020, with the exception of the B treatment in 2019 (Figure 3C). The abundance of Bacillus was significantly higher in cover crop treatments than in CK, with the exception of the W, R, and B treatments in 2019 (Figure 3D). Relative weight analysis showed that soil pH was the major factor affecting the abundance of bacteria, and it explained approximately 39.43% of the variation in the 16SrRNA gene copy number (Figure 4A).
Figure 3
Figure 4
Bacterial Community Diversity and Composition
The number of OTUs of the bacterial community was significantly increased in the B treatment compared with CK (Figure 5A). The Shannon index of the bacterial community was significantly higher in the WB treatment than in CK (Figure 5B). MRT analysis revealed that bacterial diversity was mainly shaped by soil available P (Figure 4B).
Figure 5
Actinobacteriota, Acidobacteriota, Proteobacteria, Chloroflexi, Gemmatimonadota, Bacteroidota, and Myxococcota were the dominant phyla in all samples, accounting for more than 90% of the total sequences. The relative abundances of Firmicutes, Patescibacteria, and Verrucomicrobiota were low (relative abundance <4%; Figure 6A). The relative abundances of Vicinamibacteria, Alphaproteobacteria, Actinobacteria, Gammaproteobacteria, Gemmatimonadetes, Thermoleophilia, Acidimicrobiia, and Bacteroidia were higher in all soil samples (relative abundance >3%), and these classes made up more than 66% of the total bacterial sequences. Eleven other bacterial classes, including Anaerolineae, Blastocatellia, Polyangia, KD4-96, and Bacilli were less abundant (relative abundance <3%) in all soil samples (Figure 6B).
Figure 6
The relative abundances of Streptomyces, Arthrobacter, and Bacillus were significantly increased in the W and WRB treatments compared with CK. The relative abundance of Streptomyces was significantly decreased in the B treatment relative to CK. The relative abundance of Lechevalieria was higher in the R and WR treatments than in CK. The relative abundances of Nitrospira and Gaiella were significantly increased in the WB and B treatments, respectively, compared with CK (Figure 7A). Correlation analysis indicated that soil physicochemical properties were closely correlated with the abundance of most of the genera (Figure 7B). The relative abundances of Arthrobacter, Bacillus, and Streptomyces were positively related to pH and NO3−, and negatively related to C/N. In contrast to Bacillus and Streptomyces, the relative abundance of Arthrobacter was positively related to AK and TN. The relative abundance of Lechevalieria was positively related to pH but negatively related to soil AK.
Figure 7
Structure of the Bacterial Community and Its Correlation With Soil Properties
Principal coordinate analysis of bacterial communities at the OTU level revealed variation in the bacterial community composition among samples; the first and second principal coordinates together explained 29.87% of the variation in community structure (Figure 8A). Among sole cover crop treatments, the community structures of the R and B treatments were most similar to that of CK, and the community structures of these treatments differed to that of the W treatment (Figure 8B). Among mixed species cover crop treatments, the community structures of the WR and RB treatments were most similar to that of CK, and the community structures of these treatments differed from the community structures of the WB and WRB treatments (Figure 8C). RDA (Figure 9) and the Mantel tests (Table 2) were conducted to determine the relationship between soil physicochemical properties and soil bacterial community composition. Among the environmental factors, pH and NH4+ were significantly correlated with the soil bacterial community composition.
Figure 8
Figure 9
Table 2
| Soil | r | p |
|---|---|---|
| pH | 0.343 | 0.007* |
| EC | 0.173 | 0.092 |
| NO-3 | 0.169 | 0.057 |
| NH + 4 | 0.254 | 0.013* |
| AP | 0.023 | 0.429 |
| AK | 0.005 | 0.459 |
| SOM | −0.026 | 0.620 |
| TN | 0.128 | 0.094 |
| C/N | 0.027 | 0.365 |
Mantel test correlations between bacteria community structure and soil physicochemical properties.
Only significant correlations (*p < 0.05) were shown in bold number. TN, soil total nitrogen; AP, available phosphorus; AK, available potassium; and SOM, soil organic matter.
Indicator Species Analysis
Indicator species analysis was carried out to identify the important OTUs in the different cover crop treatments. We identified indicator species, which are used as ecological indicators of habitat types or communities, based on the relative abundance and frequency of species. Two hundred and sixty four OTUs were identified to be sensitive to cover crop treatment. Most of these sensitive OTUs belonged to Proteobacteria, Actinobacteriota, and Patescibacteria phyla. Only eight OTUs were enriched solely in the CK. A total of 43, 51, and 26 OTUs were enriched solely in the W, WB, and WRB treatments, respectively (Figure 10).
Figure 10
Discussion
Effect of Cover Crops on Pepper Yield and Soil Characteristics
Cover crops usually have positive effects on main crop production and soil properties (e.g., moisture and nutrient status), but their effects greatly depend on the crop species (Marinari et al., 2015). In this study, pepper yield was increased in the WR and WB treatments in 2019 and 2020 and the WRB treatment in 2020 compared with CK (Figure 2). A positive correlation between pepper yield and soil available K was found regardless of year in our study (Supplementary Table S2), which was consistent with the finding reported by Tao et al. (2016), indicating that enhanced pepper yield might result from the increase of soil nutrient elements. Furthermore, we found that the pepper production in cover crop treatments in 2020 was higher than that in 2019 (Figure 2). Previous study showed that the net effects of biodiversity on ecosystem functions grow stronger with time because the magnitude of resource complementarity increases as experiment run longer (Cardinale et al., 2012). Soil nutrient status and pH are the important indices of soil properties (Raaijmakers et al., 2009; Thiele-Bruhn et al., 2012), and the deterioration of these properties is the main cause of “soil sickness” (Chang et al., 2017). In our study, we found that cover crops decreased the soil available P content in 2 years (Table 1). This finding is consistent with the results of Gao et al. (2017) showing that lower soil available P content might stem from the higher uptake of main crops or cover crops. In contrast, Wang et al. (2021) found that long-term cover crops could improve soil P availability. These results suggest that the duration of cover crops treatment might have an important effect on soil P nutrient. In addition, we found that cover crop treatments increased soil pH as compared to CK (Table 1). Similar results were found in the study reported by Balota and Chaves (2011) where pH was increased from 4 to 6 when cultivated with different summer legumes in the coffee crop. However, Wang et al. (2021) found that legumes as cover crop decreased soil pH and it was attributed to the secretion of acids and protons from their roots (Zhou et al., 2017). In our present study, despite the general increase in soil pH after being treated by cover crops, the pH was lowest in the B treatment, and this is in agreement with previous studies showing a higher soil pH in non-legumes than in legumes (Li et al., 2007; Maltais-Landry, 2015). The concentration of NH4+ was decreased in cover crop treatments after 2 years, which might be caused by the increase in nitrification activity stemming from the higher pH in cover crop treatments, which promoted the conversion of NH4+-N to NO3−-N (Islam et al., 2006; Lin et al., 2018). Soil microorganisms are involved in soil nutrient cycling (van der Heijden et al., 2008). Therefore, cover crops might indirectly affect soil properties through plant-mediated changes in the soil microbial community.
Effect of Cover Crops on Bacterial Abundance and Diversity
Increasing evidence indicates that cover crops can alter soil microbial communities (Nair and Ngouajio, 2012; Marinari et al., 2015; Detheridge et al., 2016). In this study, we found that cover crop treatments significantly increased the abundance of bacteria and fungi in 2020 (Figures 3A,B). This was consistent with the general observation that plant diversity has positive effects on soil microbial biomass due to the increases in soil carbon resources (Eisenhauer et al., 2010). Moreover, we also found that the abundance of bacteria and fungi in cover crop treatments in 2020 was higher than that in 2019 (Figures 3A,B). This might be partly explained by the time-lag response of soil biota to plant diversity changes (Eisenhauer et al., 2010, 2011). Soil microbiota have a significant effect on plant productivity in natural ecosystems (Wagg et al., 2014). Thus, we deduced that the difference of yield in different years might be related to the change of soil microbiota. Soil bacterial communities are sensitive to soil characteristics, especially to the soil pH (Kuramae et al., 2011; Shen et al., 2013). Similar result was found in our present study where the abundance of bacteria was mainly changed by soil pH (Figure 4A). Thus, our results suggest that the increased abundance of bacteria in the cover crop treatments compared with CK might be related to the higher pH.
High-throughput amplicon sequencing analysis in our study found that the Shannon index of the bacterial community only in the WB treatment was significantly higher than that of the CK in 2020 (Figure 5B), indicating our first hypothesis was not fully validated. Theoretically, higher plant diversity can be an efficient way to increase soil microbiota diversity due to the increase in soil habitat heterogeneity and carbon resource diversity (Eisenhauer et al., 2010, 2011), but some studies found that plant diversity have positive, neutral, or negative effects on soil microbe diversity (Wardle, 2006; Schlatter et al., 2015). It takes 4–6 years for soil biota to respond significantly to the changes in plant diversity (Eisenhauer et al., 2011), thus these inconsistent results may derive from the different duration of cover crops. Previous study found that multi-mixtures of grass and legumes improve microbial diversity when compared with single and control treatments (Vukicevich et al., 2016). Similarly, our results found that the bacterial diversity was higher in WB than in W. The possible reason could be that a mixture of high C/N cover crop (wheat) with legumes facilitates nutrient mineralization and provide a C source for diverse groups of soil organisms (Carrera et al., 2007; Dahal et al., 2020). MRT analysis indicated that bacterial diversity was mainly shaped by soil available P (Figure 4B). Long-term application of high phosphorus fertilizers to soil decreased the diversity of total bacteria (Liu et al., 2020b). In our research, the concentration of available P in CK treatment was higher than that in high P fertilizer treatment reported by Liu et al. (2020b). Thus, we also deduced that the higher Shannon index in the WB treatment was partly attributed to the lower soil available P.
Effect of Cover Crops on the Composition of Bacterial Communities
In our study, the main phyla in the cover crop system were Proteobacteria, Acidobacteria, Bacteroidetes, Actinobacteria, and Firmicutes (Figure 6A), which are related to soil C degradation (DeAngelis et al., 2013; Trivedi et al., 2015). This result is consistent with the finding of Zheng et al. (2018a). Moreover, we found that Firmicutes were more abundant in the W and WRB treatments as compared to CK, indicating that they were the main decomposing bacteria of cover crop residues in the W and WRB treatments. This was also supported by an analysis conducted at the genus levels (the relative abundance of Bacillus spp. was increased in the W and WRB treatments relative to CK (Figure 7A).
The positive interactions between plant and microbe can result from the increasing of abundance of plant-beneficial microorganisms (Bever et al., 2012). Our high-throughput amplicon sequencing analysis showed that cover crop treatments increased the relative abundances of several bacterial taxa, which have plant-growth-promoting and/or pathogen-antagonistic potentials in pepper rhizosphere (Manzanera et al., 2015; Zheng et al., 2018b; Mun et al., 2020). For example, the relative abundances of Streptomyces, Arthrobacter, and Bacillus spp. were significantly increased in both the W and WRB treatments as compared to CK (Figure 7A). The relative abundance of Pseudomonas spp. which was lower than 0.5% in all treatments, was significantly increased in the WRB treatment as compared to CK (Supplementary Figure S1). Similarly, quantitative PCR found that W and WRB treatments significantly increased the abundance of Bacillus and Pseudomonas spp. in 2020 (Figures 3C,D). Wheat can release malic acid and benzoxazinoids in root exudates, which attract Bacillus and Pseudomonas spp., respectively (Fang et al., 2008; Berendsen et al., 2012). Therefore, the increase in the relative abundance of potential plant-beneficial bacteria in the W and WRB treatments might be partly attributed to the release of secondary metabolites from wheat root exudates or litter decomposition (Berendsen et al., 2012). The sensitivity of bacterial genera to changes in the soil environment varies, thus, their responses to soil environmental variables might differ (Sun et al., 2020). In our study, the higher pH and NO3− content and lower C/N ratio in the W and WRB treatments might be another reason for the higher relative abundances of Streptomyces, Arthrobacter, and Bacillus spp., since the relative abundances of these three genera were positively related to the pH and soil NO3− content, but negatively related to the C/N ratio (Figure 7B). In addition, we found that the relative abundance of Lechevalieria was higher in the R and WR treatments than in CK (Figure 7A). The antibiotic produced by the genus Lechevalieria has been found highly effective against plant fungal pathogens (Lee et al., 2004). Thus, we deduced that the R and WR treatments might enhance the resistance of plants to pathogens. This hypothesis is consistent with a previous study showing that the use of wild rocket as green manure decreased the severity of cucumber Fusarium wilt disease (Jin et al., 2019). The positive correlation between the relative abundance of Lechevalieria and soil pH (Figure 7B) indicated that increasing soil pH could be beneficial for its growth (Aslam et al., 2020).
The change of soil microbial community structure has been shown to be affected by substrate availability and quality associated with plant diversity (Nautiyal et al., 2010; Mbuthia et al., 2015). Therefore, cover crop mixtures as increasing species diversity would result in a significant shift in microbial community structure by producing a variety of substrates and increasing habitat heterogeneity (Spedding et al., 2004; Carrera et al., 2007). A previous study showed that the use of oats and mixtures containing oats as cover crops increased the accumulation of total N and SOC compared with brassica-only and legume-only treatments (Ghimire et al., 2019). In our study, the use of wheat as a cover crop by itself or mixed with faba bean might have provided a continuous food supply for microbial communities and promoted the growth of soil microbe (Thapa et al., 2021). This may explain the shifts in bacterial community structure under wheat treatment and mixture treatments containing wheat (e.g., W, WB, and WRB), relative to CK (Figure 8). However, we found that the community structure of WR differed to that of the W (Figure 8A), which possibly owing to the difference in the soil environment. Wild rocket as a member of the Brassicaceae family contain glucosinolates in their tissue, which are the precursor for the production of bioactive compounds isothiocyanates (White and Weil, 2010; Koide and Peoples, 2012). Isothiocyanates could alter bacterial community structure (Cohen et al., 2005). However, our results found that the community structure of R and mixtures containing R as cover crops were similar to that of CK (Figure 8). Similarly, Klose et al. (2006) did not observe apparent changes of soil microbial community structure after treating with methyl isothiocyanate. We suggest, this may be due to the short half-life time of these compounds and can disappear within a few days after their release into the soil environment (Motisi et al., 2009; Hanschen et al., 2015). Several studies have shown that soil physicochemical properties, especially soil pH, are important factors affecting soil bacterial community structures (Shen et al., 2013; Liu et al., 2014). Our present study found that W, WB, and WRB treatments had higher pH and lower NH4+ concentration compared with CK, which stems from continued N fertilization (especially NH4+-N) of the continuous soil and/or the release of H+ into the soil by nitrification (Yu et al., 2019). RDA and Mantel tests were performed to further confirm that soil pH and NH4+ were significantly correlated with the soil bacterial community composition, which supported our second hypothesis. These results were in accordance with the results of Yu et al. (2019). Thus, cover crop treatments likely altered soil bacterial community structure through their effects on the chemical properties of the soil. This conclusion is consistent with the results of previous studies showing that plants type might indirectly affect the bacterial community by altering the pH and N status of soil (Nielsen et al., 2010; Shen et al., 2013).
Indicator species are considered as ecological indicators of environmental changes, community type, and habitat conditions because of their niche preferences (Cáceres and Legendre, 2009). PCoA analysis showed that the bacterial community structure of W, WB, and WRB treatments were different from CK (Figure 8). In this study, there were 264 OTUs that differed among the eight treatments (Figure 10). Only eight OTUs were enriched solely in the CK. However, a total of 43, 51, and 26 OTUs were enriched solely in the W, WB, and WRB treatments, respectively. The indicator species in W (OTU2333, OTU2959, and OTU755) and WRB (OTU5821, OTU3068, and OTU3228) treatments were most closely related to Pseudomonas, Bacillus, and Streptomyces spp., respectively. Root exudates from wheat play key role in recruiting several species of the Pseudomonas and Bacillus taxa (Berendsen et al., 2012). For WB treatment, the indicator species OTU4588 was most closely related to Gaiella (Figure 10), which has been shown to play a key role in the resistance of strawberry plants against F. oxysporum (Zhao et al., 2019). These results suggest that cover crop treatments altered the bacterial community structure with shifts in the dominant genus, which have plant-growth-promoting and/or pathogen-antagonistic potentials in W, WB, and WRB treatments. The root exudates from wheat increases in beneficial microorganisms may be responsible for the enhanced promoting effect of changes in soil biota induced by cover crops on pepper seedling growth; however, additional work is needed to test this hypothesis.
Conclusion
In this study, the effects of cover crops alone and their different combinations were evaluated on pepper yield, soil physicochemical properties, and soil bacterial communities in a continuous pepper cropping system. We found that the pepper yield was increased in the WR and WB in 2019 and 2020 and in the WRB in 2020. Cover crops increased the soil pH, decreased the concentrations of NH4+ and available phosphorus, and altered the bacterial abundance and community composition. A shift in certain bacterial taxa with mostly potential plant-beneficial bacteria specific to different cover crop combinations was observed. Moreover, soil pH was a strong indicator of these changes. Therefore, we conclude that mixing different cover crops could affect the soil environment differentially and regulate soil microbial community composition to benefit crop production. We also assumed that the enrichment of different taxa specific to different crop combination could be the result of differences in crop root exudates composition, which play critical role in rhizosphere assembly of microbial community. However, the assumption needs further study to underpin the detailed mechanism of how diversified cover crops regulate soil microbial community and increase crop yield.
Funding
This work was funded by National Key Research and Development Program of China (2019YFD1001900 and 2018YFD1000800), and Ministry of Finance, Ministry of Agriculture, and Rural Affairs, National Modern Industrial Technology System (CARS).
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.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number can be found below: https://www.ncbi.nlm.nih.gov/, SRP326928.
Author contributions
FW contributed to design this experiment. HG and GT performed the experiment. HG analyzed the data and wrote the manuscript. MR revised the manuscript. All authors contributed to the article and approved the submitted version.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.789034/full#supplementary-material
Supplementary Figure S1The relative abundance of Pseudomonas spp. in different treatments. Different lowercase letters in each column denote significant differences among the samples (p < 0.05, HSD Tukey test). CK, no cover crop; W, wheat; R, wild rocket; B, faba bean; WR, wheat + wild rocket; WB, wheat + faba bean; RB, wild rocket + faba bean; and WRB, wheat + wild rocket + faba bean.
References
1
AslamS.ShierW. T.SajidI. (2020). Taxonomic diversity, antimicrobial potential and metabolite profiling of aquatic actinobacteria from Kallar Kahar Lake. Pak. J. Zool.52:3. doi: 10.17582/journal.pjz/20190609230628
2
BalotaE. L.ChavesJ. C. D. (2011). Microbial activity in soil cultivated with different summer legumes in coffee crop. Braz. Arch. Biol. Technol.54, 35–44. doi: 10.1590/S1516-89132011000100005
3
BellemainE.CarlsenT.BrochmannC.CoissacE.TaberletP.KauserudH. (2010). ITS as an environmental DNA barcode for fungi: an in silico approach reveals potential PCR biases. BMC Microbiol.10:189. doi: 10.1186/1471-2180-10-189
4
BerendsenR. L.PieterseC.BakkerP. (2012). The rhizosphere microbiome and plant health. Trends Plant Sci.17, 478–486. doi: 10.1016/j.tplants.2012.04.001
5
BeverJ. D.PlattT. G.MortonE. R. (2012). Microbial population and community dynamics on plant roots and their feedbacks on plant communities. Annu. Rev. Microbiol.66, 265–283. doi: 10.1146/annurev-micro-092611-150107
6
BirousteM.KazakouE.BlanchardA.RoumetC. (2012). Plant traits and decomposition: are the relationships for roots comparable to those for leaves?Ann. Bot.109, 463–472. doi: 10.1093/aob/mcr297
7
CáceresM. D.LegendreP. (2009). Associations between species and groups of sites: indices and statistical inference. Ecology90, 3566–3574. doi: 10.1890/08-1823.1
8
CaporasoJ. G.KuczynskiJ.StombaughJ.BittingerK.BushmanF. D.CostelloK. E.et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods7, 335–336. doi: 10.1038/nmeth.f.303
9
CardinaleB. J.DuffyJ. E.GonzalezA.HooperD. U.PerringsC.VenailP.et al. (2012). Biodiversity loss and its impact on humanity. Nature486, 59–67. doi: 10.1038/nature11148
10
CarreraL. M.BuyerJ. S.VinyardB.Abdul-BakiA. A.SikoraL. J.Jr.TeasdaleJ. R. (2007). Effects of cover crops, compost, and manure amendments on soil microbial community structure in tomato production systems. Appl. Soil Ecol.37, 247–255. doi: 10.1016/j.apsoil.2007.08.003
11
ChangC. L.FuX. P.ZhouX. G.GuoM. Y.WuF. Z. (2017). Effects of seven different companion plants on cucumber productivity, soil chemical characteristics and pseudomonas community. J. Integr. Agric.16, 2206–2214. doi: 10.1016/S2095-3119(17)61698-8
12
CohenM. F.YamasakiH.MazzolaM. (2005). Brassica napus seed meal soil amendment modifies microbial community structure, nitric oxide production and incidence of Rhizoctonia root rot. Soil Biol. Biochem.37, 1215–1227. doi: 10.1016/j.soilbio.2004.11.027
13
DahalS.FranklinD.SubediA.CabreraM.HancockD.MahmudK.et al. (2020). Strategic grazing in beef-pastures for improved soil health and reduced runoff-nitrate-a step towards sustainability. Sustain. For.12:558. doi: 10.3390/su12020558
14
DeAngelisK. M.ChivianD.FortneyJ. L.ArkinA. P.SimmonsB.HazenT. C.et al. (2013). Changes in microbial dynamics during long-term decomposition in tropical forests. Soil Biol. Biochem.66, 60–68. doi: 10.1016/j.soilbio.2013.06.010
15
DerakhshaniH.TunH. M.KhafipourE. (2016). An extended single-index multiplexed 16S rRNA sequencing for microbial community analysis on MiSeq illumina platforms. J. Basic Microbiol.56, 321–326. doi: 10.1002/jobm.201500420
16
DetheridgeA. P.BrandG.FychanR.CrottyF. V.SandersonR.GriffithG. W.et al. (2016). The legacy effect of cover crops on soil fungal populations in a cereal rotation. Agric. Ecosyst. Environ.228, 49–61. doi: 10.1016/j.agee.2016.04.022
17
EdgarR. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods10, 996–998. doi: 10.1038/nmeth.2604
18
EilersK. G.LauberC. L.KnightR.FiererN. (2010). Shifts in bacterial community structure associated with inputs of low molecular weight carbon compounds to soil. Soil Biol. Biochem.42, 896–903. doi: 10.1016/j.soilbio.2010.02.003
19
EisenhauerN.BesslerH.EngelsC.GleixnerG.HabekostM.MilcuA.et al. (2010). Plant diversity effects on soil microorganisms support the singular hypothesis. Ecology91, 485–496. doi: 10.1890/08-2338.1
20
EisenhauerN.MilcuA.SabaisA. C. W.BesslerH.BrennerJ.EngelsC.et al. (2011). Plant diversity surpasses plant functional groups and plant productivity as driver of soil biota in the long term. PLoS One6:e1605. doi: 10.1371/journal.pone.0016055
21
FangZ. Y.AnZ. F.LiY. L. (2008). Dynamic change of organic acids secreted from wheat roots in Mn deficiency. Front. Agric. China2, 50–54. doi: 10.1007/s11703-008-0011-9
22
FiererN.BradfordM. A.JacksonR. B. (2007). Toward an ecological classification of soil bacteria. Ecology88, 1354–1364. doi: 10.1890/05-1839
23
GaoD. M.ZhouX. G.DuanY. D.FuX. P.WuF. Z. (2017). Wheat cover crop promoted cucumber seedling growth through regulating soil nutrient resources or soil microbial communities?Plant Soil418, 459–475. doi: 10.1007/s11104-017-3307-9
24
GarbevaP.van VeenJ.van ElsasJ. (2003). Predominant Bacillus spp. in agricultural soil under different management regimes detected via PCR-DGGE. Microb. Ecol.45, 302–316. doi: 10.1007/s00248-002-2034-8
25
GarbevaP.van VeenJ. A.van ElsasJ. D. (2004). Assessment of the diversity and antagonism to rhizoctonia solani AG3, of Pseudomonas species in soil from different agricultural regimes. FEMS Microbiol. Ecol.47, 51–64. doi: 10.1016/S0168-6496(03)00234-4
26
GhimireR.GhimireB.MesbahA. O.SainjuU. M.IdowuO. J. (2019). Soil health response of cover crops in winter wheat-fallow system. Agron. J.111, 2108–2115. doi: 10.2134/agronj2018.08.0492
27
GuoR. Y.LiX. L.ChristieP.ChenQ.JiangR. F.ZhangF. S. (2008). Influence of root zone nitrogen management and a summer catch crop on cucumber yield and soil mineral nitrogen dynamics in intensive production systems. Plant Soil313, 55–70. doi: 10.1007/s11104-008-9679-0
28
HanschenF. S.YimB.WinkelmannT.SmallaK.SchreinerM. (2015). Degradation of biofumigant isothiocyanates and allyl glucosinolate in soil and their effects on the microbial community composition. PLoS One10:e0132931. doi: 10.1371/journal.pone.0132931
29
HenryS.BruD.StresB.HalletS.PhilippotL. (2006). Quantitative detection of the nosz gene, encoding nitrous oxide reductase, and comparison of the abundances of 16s rRNA, narG, nirK, and nosZ genes in soils. Appl. Environ. Microbiol.72, 5181–5189. doi: 10.1128/AEM.00231-06
30
HuangL. F.SongL. X.XiaX. J.MaoW. H.ShiK.ZhouY. H.et al. (2013). Plant-soil feedbacks and soil sickness: From mechanisms to application in agriculture. J. Chem. Ecol.39, 232–242. doi: 10.1007/s10886-013-0244-9
31
IslamA.WhiteR. E.ChenD. (2006). Nitrification activity in acid soils of North-Eastern Victoria, Australia, as affected by liming and phosphorus fertilisation. Aust. J. Soil Res.44, 739–744. doi: 10.1071/SR06058
32
JinX.ZhangJ. H.ShiY. J.WuF. Z.ZhouX. G. (2019). Green manures of Indian mustard and wild rocket enhance cucumber resistance to Fusarium wilt through modulating rhizosphere bacterial community composition. Plant Soil441, 283–300. doi: 10.1007/s11104-019-04118-6
33
KabacoffR. (2015). R in Action: Data Analysis and Graphics with R.Shelter Island, NY, USA: Manning Publications Co.
34
KabirZ.KoideR. T. (2002). Effect of autumn and winter mycorrhizal cover crops on soil properties, nutrient uptake and yield of sweet corn in Pennsylvania, USA. Plant Soil238, 205–215. doi: 10.1023/A:1014408723664
35
KloseS.Acosta-MartinezV.AjwaH. A. (2006). Microbial community composition and enzyme activities in a sandy loam soil after fumigation with methyl bromide or alternative biocides. Soil Biol. Biochem.38, 1243–1254. doi: 10.1016/j.soilbio.2005.09.025
36
KoideR. T.PeoplesM. S. (2012). On the nature of temporary yield loss in maize following canola. Plant Soil360, 259–269. doi: 10.1007/s11104-012-1237-0
37
KuramaeE. E.EtienneY.WongL. C.PijlA. S.VanV. J. A.KowalchukG. A. (2011). Soil characteristics more strongly influence soil bacterial communities than land-use type. FEMS Microbiol. Ecol.79, 12–24. doi: 10.1111/j.1574-6941.2011.01192.x
38
LarkinR. P.GriffinT. S.HoneycuttC. W. (2010). Rotation and cover crop effects on soilborne potato diseases, tuber yield, and soil microbial communities. Plant Dis.94, 1491–1502. doi: 10.1094/PDIS-03-10-0172
39
LeeJ. Y.MoonS. S.YunB. S.YooI. D.HwangB. K. (2004). Thiobutacin, a novel antifungal and antioomycete antibiotic from lechevalieria aerocolonigenes. J. Nat. Prod.67, 2076–2078. doi: 10.1021/np049786v
40
LetunicI.BorkP. (2007). Interactive tree of life (iTOL): an online tool for phylogenetic tree display and annotation. Bioinformatics23, 127–128. doi: 10.1093/bioinformatics/btl529
41
LiL.LiS. M.SunJ. H.ZhouL. L.BaoX. G.ZhangH. G.et al. (2007). Diversity enhances agricultural productivity via rhizosphere phosphorus facilitation on phosphorus-deficient soils. Proc. Natl. Acad. Sci. U. S.A.104, 11192–11196. doi: 10.1073/pnas.0704591104
42
LinY. X.YeG. P.LiuD. Y.LedgardS.LuoJ. F.FanJ. B.et al. (2018). Long-term application of lime or pig manure rather than plant residues suppressed diazotroph abundance and diversity and altered community structure in an acidic Ultisol. Soil Biol. Biochem.123, 218–228. doi: 10.1016/j.soilbio.2018.05.018
43
LiuZ. X.LiuJ. J.YuZ. H.YaoQ.LiY. S.LiangA. Z.et al. (2020a). Long-term continuous cropping of soybean is comparable to crop rotation in mediating microbial abundance, diversity and community composition. Soil Tillage Res.197:104503. doi: 10.1016/j.still.2019.104503
44
LiuJ. S.MaQ.HuiX. L.RanJ. Y.MaQ. X.WangX. S.et al. (2020b). Long-term high-P fertilizer input decreased the total bacterial diversity but not phoD-harboring bacteria in wheat rhizosphere soil with available-P deficiency. Soil Biol. Biochem.149:107918. doi: 10.1016/j.soilbio.2020.107918
45
LiuJ. J.SuiY.YuZ.ShiY.ChuH.JinJ.et al. (2014). High throughput sequencing analysis of biogeographical distribution of bacterial communities in the black soils of Northeast China. Soil Biol. Biochem.70, 113–122. doi: 10.1016/j.soilbio.2013.12.014
46
Maltais-LandryG. (2015). Legumes have a greater effect on rhizosphere properties (pH, organic acids and enzyme activity) but a smaller impact on soil P compared to other cover crops. Plant Soil394, 139–154. doi: 10.1007/s11104-015-2518-1
47
ManiciL. M.CaputoF.BabiniV. (2004). Effect of green manure on Pythium spp. population and microbial communities in intensive cropping systems. Plant Soil263, 133–142. doi: 10.1023/B:PLSO.0000047720.40918.29
48
ManzaneraM.Narváez-ReinaldoJ. J.García-FontanaC.VílchezJ. I.González-LópezJ. (2015). Genome sequence of Arthrobacter koreensis 5J12A, a plant growth-promoting and desiccation-tolerant strain. Genome Announc.3, e00684–e00615. doi: 10.1128/genomeA.00648-15
49
MarinariS.RobertoM.BrunettiP.CampigliaE. (2015). Soil quality, microbial functions and tomato yield under cover crop mulching in the Mediterranean environment. Soil Tillage Res.145, 20–28. doi: 10.1016/j.still.2014.08.002
50
MbuthiaL. W.Acosta-MartínezV.DeBruynJ.SchaefferS.TylerD.OdoiE.et al. (2015). Long term tillage, cover crop, and fertilization effects on microbial community structure, activity: implications for soil quality. Soil Biol. Biochem.89, 24–34. doi: 10.1016/j.soilbio.2015.06.016
51
MooreE. B.WiedenhoeftM. H.KasparT. C.CambardellaC. A. (2014). Rye cover crop effects on soil quality in no-till corn silage-soybean cropping systems. Soil Sci. Soc. Am. J.78, 968–976. doi: 10.2136/sssaj2013.09.0401
52
MotisiN.MontfortF.DoréT.RomillacN.LucasP. (2009). Duration of control of two soilborne pathogens following incorporation of above- and below-ground residues of Brassica juncea into soil. Plant Pathol.58, 470–478. doi: 10.1111/j.1365-3059.2008.02017.x
53
MunB. G.LeeW. H.KangS. M.LeeS. U.LeeI. J. (2020). Streptomyces sp. lh 4 promotes plant growth and resistance against sclerotinia sclerotiorum in cucumber via modulation of enzymatic and defense pathways. Plant Soil448, 87–103. doi: 10.1007/s11104-019-04411-4
54
MurunguF. S.ChiduzaC.MuchaonyerwaP.MnkeniP. N. S. (2011). Decomposition, nitrogen and phosphorus mineralization from winter-grown cover crop residues and suitability for a small holder farming system in South Africa. Nutr. Cycl. Agroecosyst.89, 115–123. doi: 10.1007/s10705-010-9381-5
55
NairA.NgouajioM. (2012). Soil microbial biomass, functional microbial diversity, and nematode community structure as affected by cover crops and compost in an organic vegetable production system. Appl. Soil Ecol.58, 45–55. doi: 10.1016/j.apsoil.2012.03.008
56
NautiyalC. S.ChauhanP. S.BhatiaC. R. (2010). Changes in soil physicochemical properties and microbial functional diversity due to 14 years of conversion of grassland to organic agriculture in semi-arid agroecosystem. Soil Tillage Res.109, 55–60. doi: 10.1016/j.still.2010.04.008
57
NielsenU. N.OslerG. H. R.CampbellC. D.BurslemD. F. R. P.van der WalR. (2010). The influence of vegetation type, soil properties and precipitation on the composition of soil mite and microbial communities at the landscape scale. J. Biogeogr.37, 1317–1328. doi: 10.1111/j.1365-2699.2010.02281.x
58
NjeruE. M.AvioL.SbranaC.TurriniA.BocciG.BàrberiP.et al. (2014). First evidence for a major cover crop effect on arbuscular mycorrhizal fungi and organic maize growth. Agron. Sustain. Dev.34, 841–848. doi: 10.1007/s13593-013-0197-y
59
PatkowskaE.Baewicz-WoniakM.KonopińskiM.WachD. (2016). The effect of cover crops on the fungal and bacterial communities in the soil under carrot cultivation. Plant Soil Environ.62, 237–242. doi: 10.17221/117/2016-PSE
60
QuastC.PruesseE.YilmazP.GerkenJ.SchweerT.YarzaP.et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res.41, D590–D596. doi: 10.1093/nar/gks1219
61
RaaijmakersJ. M.PaulitzT. C.SteinbergC.AlabouvetteC.Moënne-LoccozY. (2009). The rhizosphere: A playground and battlefield for soilborne pathogens and beneficial microorganisms. Plant Soil321, 341–361. doi: 10.1007/s11104-008-9568-6
62
RorickJ. D.KladivkoE. J. (2016). Cereal rye cover crop effects on soil carbon and physical properties in southeastern Indiana. J. Soil Water Conserv.72, 260–265. doi: 10.2489/jswc.72.3.260
63
SarrantonioM.GallandtE. (2003). The role of cover crops in north American cropping systems. J. Crop. Prod.8, 53–74. doi: 10.1300/J144v08n01_04
64
SchipanskiM. E.BarbercheckM.DouglasM. R.FinneyD. M.HaiderK.KayeJ. P.et al. (2014). A framework for evaluating ecosystem services provided by cover crops in agroecosystems. Agric. Syst.125, 12–22. doi: 10.1016/j.agsy.2013.11.004
65
SchlatterD. C.BakkerM. G.BradeenJ. M.KinkelL. L. (2015). Plant community richness andmicrobial interactions structure bacterial communities in soil. Ecology96, 134–142. doi: 10.1890/13-1648.1
66
ShenC. C.XiongJ. B.ZhangH. Y.FengY. Z.LinX. G.LiX. Y.et al. (2013). Soil pH drives the spatial distribution of bacterial communities along elevation on Changbai Mountain. Soil Biol. Biochem.57, 204–211. doi: 10.1016/j.soilbio.2012.07.013
67
SnappS. S.SwintonS. M.LabartaR.MutchD.BlackJ. R.LeepR.et al. (2005). Evaluating cover crops for benefits, costs and performance within cropping system niches. Agron. J.97, 322–332. doi: 10.2134/agronj2005.0322a
68
SpeddingT. A.HamelC.MehuysG. R.MadramootooC. A. (2004). Soil microbial dynamics in maize-growing soil under different tillage and residue management systems. Soil Biol. Biochem.36, 499–512. doi: 10.1016/j.soilbio.2003.10.026
69
SullivanP. G.ParrishD. J.LunaJ. M. (1991). Cover crop contributions to N supply and water conservation in corn production. Am. J. Alternative Agr.6, 106–113. doi: 10.1017/S0889189300003982
70
SunL.LiJ.WangQ.ZhangY. H.XuZ. G.WangR.et al. (2020). Effects of eight years of conservation tillage on the soil physicochemical properties and bacterial communities in a rain-fed agroecosystem of the loess plateau, China. Land Degrad. Dev.31, 2475–2489. doi: 10.1002/ldr.3619
71
TaoJ. M.LiuX. D.LiangY. L.NiuJ. J.XiaoY. H.GuY. B.et al. (2016). Maize growth responses to soil microbes and soil properties after fertilization with different green manures. Appl. Microbiol. Biotechnol.101, 1289–1299. doi: 10.1007/s00253-016-7938-1
72
ThapaV. R.GhimireR.Acosta-MartínezV.MarsalisM. A.SchipanskiM. E. (2021). Cover crop biomass and species composition affect soil microbial community structure and enzyme activities in semiarid cropping systems. Appl. Soil Ecol.157:103735. doi: 10.1016/j.apsoil.2020.103735
73
Thiele-BruhnS.BloemJ.de VriesF. T.KalbitzK.WaggC. (2012). Linking soil biodiversity and agricultural soil management. Curr. Opin. Environ. Sustain.4, 523–528. doi: 10.1016/j.cosust.2012.06.004
74
TianY.ZhangX.LiuJ.GaoL. H. (2011). Effects of summer cover crop and residue management on cucumber growth in intensive Chinese production systems: soil nutrients, microbial properties and nematodes. Plant Soil339, 299–315. doi: 10.1007/s11104-010-0579-8
75
TrivediP.RochesterI. J.TrivediC.Van NostrandJ. D.ZhouJ.KarunaratneS.et al. (2015). Soil aggregate size mediates the impacts of cropping regimes on soil carbon and microbial communities. Soil Biol. Biochem.91, 169–181. doi: 10.1016/j.soilbio.2015.08.034
76
van der HeijdenM. G.BardgettR. D.van StraalenN. M. (2008). The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol. Lett.11, 296–310. doi: 10.1111/j.1461-0248.2007.01139.x
77
VillamilM. B.BolleroG. A.DarmodyR. G.SimmonsF. W.BullockD. G. (2006). No-till corn/soybean systems including winter cover crops: effects on soil properties. Soil Sci. Soc. Am. J.70, 1936–1944. doi: 10.2136/sssaj2005.0350
78
VukicevichE.LoweryT.BowenP.Úbez-TorresJ. R.HartM. (2016). Cover crops to increase soil microbial diversity and mitigate decline in perennial agriculture. A review. Agron. Sustain. Dev.36:48. doi: 10.1007/s13593-016-0385-7
79
WaggC.BenderS. F.WidmerF.van der HeijdenM. G. A. (2014). Soil biodiversity and soil community composition determine ecosystem multifunctionality. Proc. Natl. Acad. Sci. U. S. A.111, 5266–5270. doi: 10.1073/pnas.1320054111
80
WangY.HuangQ.GaoH.ZhangR.LiG. (2021). Long-term cover crops improved soil phosphorus availability in a rain-fed apple orchard. Chemosphere275:130093. doi: 10.1016/j.chemosphere.2021.130093
81
WardleD. A. (2006). The influence of biotic interactions on soil biodiversity. Ecol. Lett.9, 870–886. doi: 10.1111/j.1461-0248.2006.00931.x
82
WhiteC. M.WeilR. R. (2010). Forage radish and cereal rye cover crop effects on mycorrhizal fungus colonization of maize roots. Plant Soil328, 507–521. doi: 10.1007/s11104-009-0131-x
83
YangL.BaiJ. S.ZengN. H.ZhouX.LiaoY. L.LuY. H.et al. (2019). 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
84
YuL. L.TangY. L.WangZ. G.GouY. G.WangJ. W. (2019). Nitrogen-cycling genes and rhizosphere microbial community with reduced nitrogen application in maize/soybean strip intercropping. Nutr. Cycl. Agroecosyst.113, 35–49. doi: 10.1007/s10705-018-9960-4
85
ZhaoF. Y.ZhangY. Y.DongW. G.ZhangY. Q.ZhangG. X.SunZ. P.et al. (2019). Vermicompost can suppress fusarium oxysporum f. sp. lycopersici via generation of beneficial bacteria in a long-term tomato monoculture soil. Plant Soil440, 491–505. doi: 10.1007/s11104-019-04104-y
86
ZhengW.GongQ. L.ZhaoZ. Y.LiuJ.ZhaiB. N.WangZ. H.et al. (2018a). Changes in the soil bacterial community structure and enzyme activities after intercrop mulch with cover crop for eight years in an orchard. Eur. J. Soil Biol.86, 34–41. doi: 10.1016/j.ejsobi.2018.01.009
87
ZhengW. J.ZengS. Q.BaisH.LamannaJ. M.HusseyD. S.JacobsonD. L.et al. (2018b). Plant growth-promoting rhizobacteria (PGPR) reduce evaporation and increase soil water retention. Water Resour. Res.54, 3673–3687. doi: 10.1029/2018WR022656
88
ZhouY.ZhuH.FuS.YaoQ. (2017). Variation in soil microbial community structure associated with different legume species is greater than that associated with different grass species. Front. Microbiol.8:1007. doi: 10.3389/fmicb.2017.01007
Summary
Keywords
cover crop, bacterial community, continuous cropping, pepper, bacterial diversity
Citation
Gao H, Tian G, Khashi u Rahman M and Wu F (2022) Cover Crop Species Composition Alters the Soil Bacterial Community in a Continuous Pepper Cropping System. Front. Microbiol. 12:789034. doi: 10.3389/fmicb.2021.789034
Received
04 October 2021
Accepted
22 November 2021
Published
03 January 2022
Volume
12 - 2021
Edited by
Tofazzal Islam, Bangabandhu Sheikh Mujibur Rahman Agricultural University, Bangladesh
Reviewed by
Esther Menendez, University of Evora, Portugal; Tengxiang Lian, South China Agricultural University, China
Updates
Copyright
© 2022 Gao, Tian, Khashi u Rahman and Wu.
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: Fengzhi Wu, fzwu2006@aliyun.com
This article was submitted to Microbe and Virus Interactions With Plants, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.