Huan Gao
Gangming Tian1,2
Muhammad Khashi u Rahman
Fengzhi Wu- 1Department of Horticulture and Landscape Architecture, Northeast Agricultural University, Harbin, China
- 2Key Laboratory of Biology and Genetic Improvement of Horticultural Crops (Northeast Region), Ministry of Agriculture and Rural Affairs, Northeast Agricultural University, Harbin, China
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. Diagram of cover crop-pepper rotation and timeline of soil sampling for soil microbial community and physicochemical properties analysis.
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. Pepper yield after cultivation of different cover crops in 2019 and 2020. CK, no cover crop; cover crop treatments: 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. Different lowercase letters in each column denote significant differences among the samples [p < 0.05, one-way ANOVA; honestly significant difference (HSD) Tukey test].
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. Soil physicochemical properties in different treatments in 2019 and 2020.
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. The abundances of total bacteria (A), total fungi (B), Pseudomonas spp. (C), and Bacillus spp. (D) in 2019 and 2020 in the pepper rhizosphere cultivated after different cover crops. 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. Different lowercase letters in each column denote significant differences among the samples (p < 0.05; HSD Tukey test).
Figure 4. Relative influence of soil physicochemical properties on abundance of bacteria was evaluated using relative weight analysis (A). Multivariate regression tree (MRT) analysis of alpha diversity (observed species, Chao 1, Shannon and Simpson indexes) of bacterial communities and soil physicochemical properties (B). Treatments and the number of samples included in the analysis are shown under bar plots. AK, available potassium; AP, available phosphorus; C/N, carbon to nitrogen ratio; SOM, soil organic matter; and TN, total nitrogen.
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. The number of operational taxonomic units (OTUs) (A) and Shannon index (B) of pepper rhizosphere bacterial communities in different treatments. 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. Different lowercase letters in each column denote significant differences among the samples (p < 0.05; HSD Tukey test).
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. Relative abundance of dominant bacterial phyla (A) and classes (B) in pepper rhizosphere. Bacterial phyla and classes with average relative abundance >1% in at least one treatment were shown. 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.
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. The bacterial significantly altered genera with relative abundances >1% in at least one treatment were shown. Different lowercase letters in each column denote significant differences among the samples (p < 0.05, HSD Tukey test; (A). Spearman correlation analysis between the relative abundances of dominant bacterial genera and soil physicochemical properties. r indicates the correlation coefficient; *p < 0.05, **p < 0.01. AK, available potassium; AP, available phosphorus; C/N, carbon to nitrogen ratio; SOM, soil organic matter; and TN, total nitrogen (B). 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.
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. Principal Coordinates Analysis (PCoA) plot of bacterial communities based on OTUs from 24 soil samples (A). The subplans of (B) and (C) indicated the separated PCoA 3D plots, which showed the response of the bacterial composition in sole cover crop treatments and in mixed species cover crop treatments, respectively. 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.
Figure 9. Redundancy analysis (RDA) showing correlations between soil bacterial community and soil physicochemical properties. 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.
Table 2. Mantel test correlations between bacteria community structure and soil physicochemical properties.
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. The maximum likelihood tree of the cover crops sensitive OTUs (CK, no cover crop; cover crop treatments: 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). The background of the OTU names is colored according to their taxonomic affiliations at the phylum level. Taxa that could be assigned to the genus level were shown as genus, otherwise were shown as OTU id. The inner strip indicates the associations between sensitive OTUs and different cover crop treatments. The size of each circle represents the relative abundances of each OTU in each cover crop treatment.
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.
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.
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).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.789034/full#supplementary-material
References
Aslam, S., Shier, W. T., and Sajid, I. (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
Balota, E. L., and Chaves, J. 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
Bellemain, E., Carlsen, T., Brochmann, C., Coissac, E., Taberlet, P., and Kauserud, H. (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
Berendsen, R. L., Pieterse, C., and Bakker, P. (2012). The rhizosphere microbiome and plant health. Trends Plant Sci. 17, 478–486. doi: 10.1016/j.tplants.2012.04.001
Bever, J. D., Platt, T. G., and Morton, E. 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
Birouste, M., Kazakou, E., Blanchard, A., and Roumet, C. (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
Cáceres, M. D., and Legendre, P. (2009). Associations between species and groups of sites: indices and statistical inference. Ecology 90, 3566–3574. doi: 10.1890/08-1823.1
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, K. E., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Cardinale, B. J., Duffy, J. E., Gonzalez, A., Hooper, D. U., Perrings, C., Venail, P., et al. (2012). Biodiversity loss and its impact on humanity. Nature 486, 59–67. doi: 10.1038/nature11148
Carrera, L. M., Buyer, J. S., Vinyard, B., Abdul-Baki, A. A., Sikora, L. J. Jr., and Teasdale, J. 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
Chang, C. L., Fu, X. P., Zhou, X. G., Guo, M. Y., and Wu, F. 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
Cohen, M. F., Yamasaki, H., and Mazzola, M. (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
Dahal, S., Franklin, D., Subedi, A., Cabrera, M., Hancock, D., Mahmud, K., 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
DeAngelis, K. M., Chivian, D., Fortney, J. L., Arkin, A. P., Simmons, B., Hazen, T. 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
Derakhshani, H., Tun, H. M., and Khafipour, E. (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
Detheridge, A. P., Brand, G., Fychan, R., Crotty, F. V., Sanderson, R., Griffith, G. 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
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604
Eilers, K. G., Lauber, C. L., Knight, R., and Fierer, N. (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
Eisenhauer, N., Bessler, H., Engels, C., Gleixner, G., Habekost, M., Milcu, A., et al. (2010). Plant diversity effects on soil microorganisms support the singular hypothesis. Ecology 91, 485–496. doi: 10.1890/08-2338.1
Eisenhauer, N., Milcu, A., Sabais, A. C. W., Bessler, H., Brenner, J., Engels, C., et al. (2011). Plant diversity surpasses plant functional groups and plant productivity as driver of soil biota in the long term. PLoS One 6:e1605. doi: 10.1371/journal.pone.0016055
Fang, Z. Y., An, Z. F., and Li, Y. L. (2008). Dynamic change of organic acids secreted from wheat roots in Mn deficiency. Front. Agric. China 2, 50–54. doi: 10.1007/s11703-008-0011-9
Fierer, N., Bradford, M. A., and Jackson, R. B. (2007). Toward an ecological classification of soil bacteria. Ecology 88, 1354–1364. doi: 10.1890/05-1839
Gao, D. M., Zhou, X. G., Duan, Y. D., Fu, X. P., and Wu, F. Z. (2017). Wheat cover crop promoted cucumber seedling growth through regulating soil nutrient resources or soil microbial communities? Plant Soil 418, 459–475. doi: 10.1007/s11104-017-3307-9
Garbeva, P., van Veen, J., and van Elsas, J. (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
Garbeva, P., van Veen, J. A., and van Elsas, J. 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
Ghimire, R., Ghimire, B., Mesbah, A. O., Sainju, U. M., and Idowu, O. J. (2019). Soil health response of cover crops in winter wheat-fallow system. Agron. J. 111, 2108–2115. doi: 10.2134/agronj2018.08.0492
Guo, R. Y., Li, X. L., Christie, P., Chen, Q., Jiang, R. F., and Zhang, F. 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 Soil 313, 55–70. doi: 10.1007/s11104-008-9679-0
Hanschen, F. S., Yim, B., Winkelmann, T., Smalla, K., and Schreiner, M. (2015). Degradation of biofumigant isothiocyanates and allyl glucosinolate in soil and their effects on the microbial community composition. PLoS One 10:e0132931. doi: 10.1371/journal.pone.0132931
Henry, S., Bru, D., Stres, B., Hallet, S., and Philippot, L. (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
Huang, L. F., Song, L. X., Xia, X. J., Mao, W. H., Shi, K., Zhou, Y. 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
Islam, A., White, R. E., and Chen, D. (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
Jin, X., Zhang, J. H., Shi, Y. J., Wu, F. Z., and Zhou, X. G. (2019). Green manures of Indian mustard and wild rocket enhance cucumber resistance to Fusarium wilt through modulating rhizosphere bacterial community composition. Plant Soil 441, 283–300. doi: 10.1007/s11104-019-04118-6
Kabacoff, R. (2015). R in Action: Data Analysis and Graphics with R. Shelter Island, NY, USA: Manning Publications Co.
Kabir, Z., and Koide, R. T. (2002). Effect of autumn and winter mycorrhizal cover crops on soil properties, nutrient uptake and yield of sweet corn in Pennsylvania, USA. Plant Soil 238, 205–215. doi: 10.1023/A:1014408723664
Klose, S., Acosta-Martinez, V., and Ajwa, H. 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
Koide, R. T., and Peoples, M. S. (2012). On the nature of temporary yield loss in maize following canola. Plant Soil 360, 259–269. doi: 10.1007/s11104-012-1237-0
Kuramae, E. E., Etienne, Y., Wong, L. C., Pijl, A. S., Van, V. J. A., and Kowalchuk, G. 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
Larkin, R. P., Griffin, T. S., and Honeycutt, C. 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
Lee, J. Y., Moon, S. S., Yun, B. S., Yoo, I. D., and Hwang, B. K. (2004). Thiobutacin, a novel antifungal and antioomycete antibiotic from lechevalieria aerocolonigenes. J. Nat. Prod. 67, 2076–2078. doi: 10.1021/np049786v
Letunic, I., and Bork, P. (2007). Interactive tree of life (iTOL): an online tool for phylogenetic tree display and annotation. Bioinformatics 23, 127–128. doi: 10.1093/bioinformatics/btl529
Li, L., Li, S. M., Sun, J. H., Zhou, L. L., Bao, X. G., Zhang, H. 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
Lin, Y. X., Ye, G. P., Liu, D. Y., Ledgard, S., Luo, J. F., Fan, J. 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
Liu, Z. X., Liu, J. J., Yu, Z. H., Yao, Q., Li, Y. S., Liang, A. 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
Liu, J. S., Ma, Q., Hui, X. L., Ran, J. Y., Ma, Q. X., Wang, X. 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
Liu, J. J., Sui, Y., Yu, Z., Shi, Y., Chu, H., Jin, J., 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
Maltais-Landry, G. (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 Soil 394, 139–154. doi: 10.1007/s11104-015-2518-1
Manici, L. M., Caputo, F., and Babini, V. (2004). Effect of green manure on Pythium spp. population and microbial communities in intensive cropping systems. Plant Soil 263, 133–142. doi: 10.1023/B:PLSO.0000047720.40918.29
Manzanera, M., Narváez-Reinaldo, J. J., García-Fontana, C., Vílchez, J. I., and González-López, J. (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
Marinari, S., Roberto, M., Brunetti, P., and Campiglia, E. (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
Mbuthia, L. W., Acosta-Martínez, V., DeBruyn, J., Schaeffer, S., Tyler, D., Odoi, E., 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
Moore, E. B., Wiedenhoeft, M. H., Kaspar, T. C., and Cambardella, C. 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
Motisi, N., Montfort, F., Doré, T., Romillac, N., and Lucas, P. (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
Mun, B. G., Lee, W. H., Kang, S. M., Lee, S. U., and Lee, I. J. (2020). Streptomyces sp. lh 4 promotes plant growth and resistance against sclerotinia sclerotiorum in cucumber via modulation of enzymatic and defense pathways. Plant Soil 448, 87–103. doi: 10.1007/s11104-019-04411-4
Murungu, F. S., Chiduza, C., Muchaonyerwa, P., and Mnkeni, P. 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
Nair, A., and Ngouajio, M. (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
Nautiyal, C. S., Chauhan, P. S., and Bhatia, C. 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
Nielsen, U. N., Osler, G. H. R., Campbell, C. D., Burslem, D. F. R. P., and van der Wal, R. (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
Njeru, E. M., Avio, L., Sbrana, C., Turrini, A., Bocci, G., Bàrberi, P., 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
Patkowska, E., Baewicz-Woniak, M., Konopiński, M., and Wach, D. (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
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., 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
Raaijmakers, J. M., Paulitz, T. C., Steinberg, C., Alabouvette, C., and Moënne-Loccoz, Y. (2009). The rhizosphere: A playground and battlefield for soilborne pathogens and beneficial microorganisms. Plant Soil 321, 341–361. doi: 10.1007/s11104-008-9568-6
Rorick, J. D., and Kladivko, E. 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
Sarrantonio, M., and Gallandt, E. (2003). The role of cover crops in north American cropping systems. J. Crop. Prod. 8, 53–74. doi: 10.1300/J144v08n01_04
Schipanski, M. E., Barbercheck, M., Douglas, M. R., Finney, D. M., Haider, K., Kaye, J. 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
Schlatter, D. C., Bakker, M. G., Bradeen, J. M., and Kinkel, L. L. (2015). Plant community richness andmicrobial interactions structure bacterial communities in soil. Ecology 96, 134–142. doi: 10.1890/13-1648.1
Shen, C. C., Xiong, J. B., Zhang, H. Y., Feng, Y. Z., Lin, X. G., Li, X. 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
Snapp, S. S., Swinton, S. M., Labarta, R., Mutch, D., Black, J. R., Leep, R., 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
Spedding, T. A., Hamel, C., Mehuys, G. R., and Madramootoo, C. 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
Sullivan, P. G., Parrish, D. J., and Luna, J. 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
Sun, L., Li, J., Wang, Q., Zhang, Y. H., Xu, Z. G., Wang, R., 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
Tao, J. M., Liu, X. D., Liang, Y. L., Niu, J. J., Xiao, Y. H., Gu, Y. 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
Thapa, V. R., Ghimire, R., Acosta-Martínez, V., Marsalis, M. A., and Schipanski, M. 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
Thiele-Bruhn, S., Bloem, J., de Vries, F. T., Kalbitz, K., and Wagg, C. (2012). Linking soil biodiversity and agricultural soil management. Curr. Opin. Environ. Sustain. 4, 523–528. doi: 10.1016/j.cosust.2012.06.004
Tian, Y., Zhang, X., Liu, J., and Gao, L. 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 Soil 339, 299–315. doi: 10.1007/s11104-010-0579-8
Trivedi, P., Rochester, I. J., Trivedi, C., Van Nostrand, J. D., Zhou, J., Karunaratne, S., 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
van der Heijden, M. G., Bardgett, R. D., and van Straalen, N. 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
Villamil, M. B., Bollero, G. A., Darmody, R. G., Simmons, F. W., and Bullock, D. 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
Vukicevich, E., Lowery, T., Bowen, P., Úbez-Torres, J. R., and Hart, M. (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
Wagg, C., Bender, S. F., Widmer, F., and van der Heijden, M. 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
Wang, Y., Huang, Q., Gao, H., Zhang, R., and Li, G. (2021). Long-term cover crops improved soil phosphorus availability in a rain-fed apple orchard. Chemosphere 275:130093. doi: 10.1016/j.chemosphere.2021.130093
Wardle, D. A. (2006). The influence of biotic interactions on soil biodiversity. Ecol. Lett. 9, 870–886. doi: 10.1111/j.1461-0248.2006.00931.x
White, C. M., and Weil, R. R. (2010). Forage radish and cereal rye cover crop effects on mycorrhizal fungus colonization of maize roots. Plant Soil 328, 507–521. doi: 10.1007/s11104-009-0131-x
Yang, L., Bai, J. S., Zeng, N. H., Zhou, X., Liao, Y. L., Lu, Y. 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
Yu, L. L., Tang, Y. L., Wang, Z. G., Gou, Y. G., and Wang, J. 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
Zhao, F. Y., Zhang, Y. Y., Dong, W. G., Zhang, Y. Q., Zhang, G. X., Sun, Z. 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 Soil 440, 491–505. doi: 10.1007/s11104-019-04104-y
Zheng, W., Gong, Q. L., Zhao, Z. Y., Liu, J., Zhai, B. N., Wang, Z. 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
Zheng, W. J., Zeng, S. Q., Bais, H., Lamanna, J. M., Hussey, D. S., Jacobson, D. 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
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
Edited by:
Tofazzal Islam, Bangabandhu Sheikh Mujibur Rahman Agricultural University, BangladeshReviewed by:
Esther Menendez, University of Evora, PortugalTengxiang Lian, South China Agricultural University, China
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, Znp3dTIwMDZAYWxpeXVuLmNvbQ==