Hongjun Yang1,2
Jiaxin Ma1,2Zhenyang Rong1,2Dandan Zeng1,2
Yuanchao Wang1,2
Shuijin Hu3,4
Wenwu Ye1,2*Xiaobo Zheng1,2*- 1Department of Plant Pathology, Nanjing Agricultural University, Nanjing, China
- 2The Key Laboratory of Integrated Management of Crop Diseases and Pests, Ministry of Education, Nanjing, China
- 3College of Resources and Environmental Sciences, Nanjing Agricultural University, Nanjing, China
- 4Department of Entomology and Plant Pathology, North Carolina State University, Raleigh, NC, United States
Returning straw to soil is an effective way to sustain or improve soil quality and crop yields. However, a robust understanding of the impact of straw return on the composition of the soil microbial communities under field conditions has remained elusive. In this study, we characterized the effects of wheat straw return on soil bacterial and fungal communities in a wheat–soybean rotation system over a 3-year period, using Illumina-based 16S rRNA, and internal transcribed region (ITS) amplicon sequencing. Wheat straw return significantly affected the α-diversity of the soil bacterial, but not fungal, community. It enhanced the relative abundance of the bacterial phylum Proteobacteria and the fungal phylum Zygomycota, but reduced that of the bacterial phylum Acidobacteria, and the fungal phylum Ascomycota. Notably, it enriched the relative abundance of nitrogen-cycling bacterial genera such as Bradyrhizobium and Rhizobium. Preliminary analysis of soil chemical properties indicated that straw return soils had significantly higher total nitrogen (TN) contents than no straw return soils. In addition, the relative abundance of fungal genera containing pathogens was significantly lower in straw return soils relative to control soils, such as Fusarium, Alternaria, and Myrothecium. These results suggested a selection effect from the 3-year continuous straw return treatment and the soil bacterial and fungal communities were moderately changed.
Introduction
Burning of wheat straw after harvest is common in rural areas, which is a major seasonal source of air pollution and haze (Marschner et al., 2011; Qu et al., 2012; Xing et al., 2018) and contributes to the content of particulate matter of ≤2.5 μm diameter (PM2.5) in China (Qu et al., 2012; Xing et al., 2018). China produces ca. 0.9 billion tons of crop straw per year, which is rich in organic matter and nutrients (Valdez-Vazquez et al., 2016; Lu et al., 2018; Ma et al., 2018). Thus, straw burning also causes great loss of renewable resources (Yang et al., 2016; Lu et al., 2018). Consequently, straw return has been promoted as a means of improving soil fertility and crop yield (Marschner et al., 2011; Guo et al., 2016; Chen Y. et al., 2017; Li et al., 2018), while mitigating particulate air pollution.
The environmental impact on soils of straw return have been well studied, including soil water potential (Yang et al., 2016), temperature (Yang et al., 2016), enzyme activities (Ji et al., 2014; Zhao et al., 2016), soil organic matter fractions (Zhao et al., 2016; Chen Z. et al., 2017; Karlsson et al., 2017), soil quality and crop productivity (Hansen et al., 2017), soil greenhouse gas emissions (Bao et al., 2016; Zhou et al., 2017), soil chemical properties (Yu et al., 2018), and soil microbial communities (Marschner et al., 2011; Chen Y. et al., 2017; Li et al., 2017; Maarastawi et al., 2018). These results provide basic understanding in terms of how straw return may change soil carbon retention, soil quality, and soil ecosystem functions, and revealed a number of positive consequences, such as reducing soil water potential, increasing soil temperature and the activities of hydrolytic enzymes, and enhancing soil microbial functional diversity. However, there were limited studies on the potential effect of straw return to soil-borne pathogens. Straw return was usually deemed to the main causal agent for the soil-borne diseases, by providing good circumstances for pathogen growth, propagation, and accumulation, which then resulted in disease epidemic (Zhen et al., 2009; Li et al., 2012; Qiao et al., 2013). Whereas, some studies indicated that straw return could improve soil ecological environment (Lenka and Lal, 2013; Peng et al., 2013; Yuan et al., 2014), increase amounts of antagonistic microbes, disturb pathogen growth, increase plant resistance to pathogens, and then control soil-borne plant diseases (Osunlaja, 1990; Bailey and Lazarovits, 2003; Zhen et al., 2009).
Due to a new environmental regulation that prohibits the burning of straw, wheat straw return has increasingly been adopted in the Huang-Huai region of China. Huang-Huai region covers a huge area in Shandong, Anhui, Jiangsu, and Henan provinces between the Yellow River and the Haihe River, and is one of the main soybean-producing areas in China and under a rotation of summer soybean and winter wheat. In this study, we used Illumina-based 16S rRNA and internal transcribed region (ITS) amplicon sequencing to characterize the effects of straw return on soil bacterial and fungal communities, respectively, based on a three-season field experiment conducted from 2015 to 2017 in a wheat–soybean cropping system in three sites of the Huang-Huai region. We addressed the following questions: What are the effects of short-term wheat straw return on the soil bacterial and fungal communities under this system? How do microbes involved in the soil nitrogen cycle and wheat/soybean diseases respond to straw return?
Materials and Methods
Field Trial and Sampling
A three-season field experiment was conducted from 2015 to 2017 in a wheat–soybean cropping system in three sites of Huang-Huai region of China, including those at Jining, Shandong Province (35°27′N, 116°35′E, sandy loam soil); Xuzhou, Jiangsu Province (34°17′N, 117°17′E, yellow loam sand soil); and Suzhou, Anhui Province (33°38′N, 117°05′E, mortar black soil) (Supplementary Figure S1A). At each location, we sampled two production fields with different treatments: no wheat straw return (N), i.e., artificial harvest of wheat, wheat stubble less than 5 cm in height, and all wheat straw removed from the field; and wheat straw return (R), i.e., all wheat straw was crushed into pieces and mulched in the soil by no tillage after wheat harvest every year. The biomass of wheat straw was about 5000 kg.hm–2. Each treatment area was 30 × 8 m and encompassed five replicate subplots (Supplementary Figure S1B). The wheat–soybean cropping system dominates all three sites, with winter wheat seeded in early October and harvested in early June, and summer soybean seeded in mid-June and harvested in late September. The winter wheat cultivar is Jimai 22 (super-high yield cultivar) and the summer soybean cultivar is Zhonghuang13 (high-protein cultivar); both are popular in Huang-Huai region of China. The plots received basal fertilization of 50 kg N-P2O5-K2O ha–1 and 15 kg urea ha–1 at wheat seeding, and 10 kg N-P2O5-K2O ha–1 were applied at soybean seeding. The agronomic management and fertilization regimes at the three sites were similar.
Bulk soil samples were collected twice annually in June and August from 2015 to 2017: at 0 day (pre-planting) and 60 days after planting. For each site, 10 subplots (480 m2; i.e., 5 replicates per treatment) were chosen and 9 cores (20 cm depth) were collected from each subplot in an S-formation sampling method using shovels (3.8 cm diameter). The 9 cores from each subplot were mixed to form one composite sample, resulting in 20 samples per site per year. In total, 180 soil samples were analyzed (3 locations × 2 treatments × 5 subplots per treatment × 3 years × 2 time points per year). The soil samples were placed into separate sterile plastic bags and transported to the laboratory on ice. Each soil sample was sieved through a 2 mm mesh to remove roots and plant detritus. The samples were stored at –20°C until required.
Soil DNA Extraction
To minimize DNA extraction bias, DNA was extracted in quadruplicate from the composite soil samples. DNA was extracted from ca. 0.25 g bulk soil using the MoBio PowerSoil DNA Isolation Kit following the manufacturer’s protocol (MoBio Laboratories Inc., Carlsbad, CA, United States). The four DNA samples were pooled, and the DNA concentration and purity were quantified using a NanoDrop 1000 Spectrophotometer (Thermo Fisher Scientific).
PCR Amplification, Library Preparation, and Sequencing
Targeted metagenomic profiling of the samples was performed as described previously (Caporaso et al., 2011; McGuire et al., 2013). The V4 region of the bacterial 16S rRNA gene was amplified using the primers 515F (5′-GTGCCAGCMGCCGCGGTAA-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′) (Peiffer et al., 2013; Edwards et al., 2015), and the ITS1 region of the fungal ITS was amplified using the barcode primers ITS1F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and ITS2 (5′-GCTGCGTTCTTCATCGATGC-3′) (Xiong et al., 2017; Yang et al., 2018). PCR reaction mixtures (25 μL) contained Enzymatics Veraseq 2.0 Master Mix (NEB Phusion High-Fidelity PCR Master Mix), 13 μL PCR-grade water, 10 μL Master Mix, 0.5 μL forward primer (10 μM), 0.5 μL reverse primer (10 μM), and 1 μL (30 ng) sample DNA. The reaction mixtures were placed in a PCR Thermal Cycler Dice (TaKaRa, Japan) and thermal cycle was performed as follows: 94°C for 1 min, followed by 30 cycles of 45 s at 94°C, 30 s at 56°C for 16S or 52°C for ITS, and 90 s at 72°C. A final extension at 72°C for 10 min was followed by a hold at 4°C. The PCR products were run on a 1% agarose gel to verify successful amplification. Unsuccessful reactions were repeated but removed from the experiment if unsuccessful a second time. PCR products were cleaned using the Agencourt AMPure XP magnetic beads purification system (Beckman Coulter) (Stark et al., 2015). The final library was quantified by determining the average molecule length using the Agilent 2100 Bioanalyzer (Agilent DNA 1000 Reagents) and by real-time quantitative PCR (qPCR) (EvaGreenTM). Qualified libraries were paired-end sequenced using the MiSeq System with the PE250 strategy (PE251 + 8 + 8 + 251; MiSeq Reagent Kit).
Bioinformatics Analyses
Paired-end reads were generated on the Illumina MiSeq platform. The sequencing data were deposited in NCBI database (BioProject ID: PRJNA524011). After removing the adaptors and primer sequences, the raw sequences were assembled according to the unique barcode using QIIME (Caporaso et al., 2010). Split sequences for each sample were merged using FLASH (v. 1.2.11) (Magoc and Salzberg, 2011). The tags were clustered into operational taxonomic units (OTUs) with a 97% threshold using UPARSE (Edgar, 2013), and chimeras were filtered out using UCHIME (v. 4.2.40) (Edgar et al., 2011). Representative bacterial sequences were taxonomically classified using Ribosomal Database Project Classifier v. 2.2 against the Silva database (v. 128) (Quast et al., 2013), and representative fungal OTUs were classified using the UNITE database (v. 7.0) (Koljalg et al., 2013), with a cutoff confidence value of 0.6. OTUs not assigned to target species were removed from the data set.
Fungal and bacterial α-diversities were estimated by calculating the OTU richness and Shannon diversity indices in Mothur (v. 1.31.2)1. Principal coordinates analyses (PCoAs) on Bray–Curtis dissimilarities were performed using the “pcoa” function in the R package2. To test for significant differences between groups of samples, an analysis of similarity (ANOSIM) was performed in Vegan with 999 permutations based on Bray–Curtis distances between samples. For Mantel tests, the Bray-Curtis and Euclidean distances were used to construct dissimilarity matrices of communities and soil characteristics, respectively, using QIIME3. Mantel tests were used to calculate the correlations between the soil bacterial/fungal community composition and soil characteristics. The linear discriminant analysis (LDA) effect size (LEfSe) method was used to evaluate bacterial and fungal taxa significantly associated with wheat straw return (Segata et al., 2011). The LEfSe analyses were performed using the Galaxy web application and workflow framework4. The α-value threshold employed for the non-parametric factorial Kruskal–Wallis sum-rank test was 0.05, and the logarithmic LDA score threshold for feature discrimination was 2.0.
Analyses of Soil Chemical Properties
As a preliminary assay, the soils from 5 replicate subplots were pooled as one composite sample for each treatment, and only triplicate aliquots of each sample were subjected to elemental analyses. Total organic carbon (TOC) levels were measured using a TOC analyzer (TOC-L, Shimadzu, Japan). Total nitrogen (TN) levels were analyzed using a fully automatic azotometer (Kjeltec, 2300, Foss, Sweden).
Statistical Analyses
Soil chemical characteristics, the relative abundances of bacteria and fungi, and α-diversity indices between no wheat straw return soils and wheat straw return soils were compared using Wilcoxon rank-sum test (P < 0.05). The relative abundances of nitrogen-cycling bacteria and plant pathogen-associated microbes were compared also using the Wilcoxon rank-sum test (P < 0.05). Spearman’s rank-correlation coefficient was used to evaluate the relationships between α-diversity indices and soil characteristics. All statistical analyses were performed using SPSS v. 20.0 (SPSS Inc.). The compositions of the microbial communities in N and R soils were compared using the Wilcoxon rank-sum test in the R package Stats.
Results
Effects of Wheat Straw Return on α-Diversity of Soil Microbial Communities
We characterized the bacterial communities in the 180 bulk soil samples by sequencing the V4 region of the 16S rRNA gene. In all, 6,610,424 high-quality sequence reads were obtained with a median of 36,272 (range 30,312–44,663) per sample. The high-quality reads were clustered into 649,012 microbial OTUs based on a >97% sequence identity threshold. Measures of within-sample diversity (α-diversity) based on rarefied OTU tables (30,312 sequences per sample) revealed that the 60 days soils samples exhibited stronger responses to straw return in bacterial communities diversity from 2015 to 2017 (Figures 1A,B), although the responses differed among the three sites, and was rarely statistically significant (Supplementary Figure S2). Significantly higher bacterial diversity was observed in R soil than in N soil in Suzhou (P < 0.05, Wilcoxon rank-sum test; the first and third years of the 60 days periods). By contrast, R soil in Xuzhou contained lower bacterial diversity than N soil (P < 0.05, Wilcoxon rank-sum test; the second years of the 0 day period). The α-diversity in Jining did not differ significantly between the N and R treatments; however, the 0 d period had a lower bacterial α-diversity than the 60 days period (P < 0.05, Wilcoxon rank-sum test) in Jining (Figures 1A,B).
Figure 1. Richness (Sobs index) and α-diversity (Shannon’s diversity index) of the bacterial (A,B) and fungal (C,D) communities in N and R soils. Horizontal bars in boxes are medians; upper and lower box edges are the 75th and 25th quartiles, respectively; and whiskers are 1.5-fold the interquartile range. Wilcoxon rank-sum test results were displayed above the top whiskers; Significant differences between treatments and/or between years within a site and a period can be distinguished with different labels (P < 0.05, Wilcoxon rank-sum test).
We also evaluated the fungal communities in the 180 bulk soil samples by sequencing of the fungal ITS1 region. In all, 6,920,939 high-quality sequences were obtained with a median read count per sample of 36,272 (range 14,959–44,984), and 66,471 microbial OTUs were obtained. Measures of α-diversity revealed a fungal community-diversity gradient between N and R, calculated based on rarefied OTU tables (14,959 sequences per sample). The fungal α-diversity values did not differ significantly between N and R treatments in the three fields (P > 0.05, Wilcoxon rank-sum test; Figures 1C,D).
Effects of Wheat Straw Return on β-Diversity of Soil Microbial Communities
Differences in the composition of the microbial communities among samples and between groups of samples (β-diversity) were analyzed based on principal coordinate analysis (PCoA) and ANOSIM. Because of the strongest differences due to field location (Figures 2A,B), the samples from each site were analyzed separately to assess the impact of straw return in more detail. The response of bacterial and fungal community composition to wheat straw return was stronger in 60 days than in 0 day. Higher correlation to the factor of straw return were observed in fungal communities than in bacterial communities, especially in JN and XZ soils, and the correlation values were both increased year by year. Although the SZ soils exhibited stronger responses to straw return in both bacterial and fungal communities with a peak value at the second year, they exhibited decreased correlation values in bacterial communities year by year in 60 days (Figures 2C–H and Supplementary Table S1).
Figure 2. Principal coordinates analysis on bacterial and fungal community differences (Bray–Curtis dissimilarities) in the groups of samples. PCoAs showing the influence of field location and year on bacterial (A) and fungal (B) community composition. PCoAs of bacterial and fungal community composition in soils from different field locations (C–H). ANOSIM was applied to test for differences in community composition due to straw return and field location. R values are shown with *P < 0.05, ∗∗P < 0.01, and ∗∗∗P < 0.001.
Differential Abundance of Bacterial and Fungal Taxa Under Wheat Straw Return
In a phylum-level analysis, bacterial OTUs were classified into 45 phyla, predominantly Proteobacteria (31.6%), Acidobacteria (17.3%), Actinobacteria (10.9%), Bacteroidetes (9.0%), Chloroflexi (5.4%), and Planctomycetes (4.8%) (Figure 3A); these accounted for 86.5% of the bacterial sequences. Acidobacteria were more abundant in N soils, and Proteobacteria were more common in R soils. Fungal OTUs were predominantly associated with the phyla Ascomycota, Basidiomycota, and Zygomycota, which accounted for 90.3% of the total fungal sequences (Figure 3B). The relative abundance of Zygomycota was higher in R than in N soils, whereas that of Ascomycota showed the opposite trend (P < 0.05, Wilcoxon test).
Figure 3. Effects of wheat straw return on the composition of the soil microbial community. Relative abundances of the most abundant bacterial (A) and fungal (B) phyla in each stage, site, year, and straw treatment.
We further assessed the impact of straw return on the abundances of individual OTUs (relative abundance > 1%). For fungi, straw return led to an enrichment of OTUs belonging to classes Mortierellales (phylum Zygomycota) and Dothideomycetes (phylum Ascomycota) (Figure 4). The phylum with the most straw return-depleted OTUs across all site soils was Ascomycota. The majority of these OTUs were classified as Sordariomycetes, mostly represented by the family Nectriaceae and genus Fusarium (Figure 4). Thus, the straw return-responsive OTUs generally followed a consistent trend within certain phyla, classes, orders, families, and genera. For bacteria, among straw return-enriched OTUs, Proteobacteria was the most highly represented phylum in R soil communities (Figure 4). The relative abundances of several OTUs classified as Chloroflexi, mainly from the class Anaerolineae, increased under straw return. Within Proteobacteria, most straw return enriched OTUs were classified as Alphaproteobacteria and Gammaproteobacteria. In particular, Alphaproteobacteria was mainly represented by the orders Sphingomonadales and Rhizobiales (mostly in SZ and XZ soils). The phylum with the most straw return-depleted OTUs across all sites was Acidobacteria (Figure 4). The majority of these OTUs were classified as Acidobacteria, with Acidobacteria and Blastocatellales being the orders most broadly affected. OTUs from the class Gemmatimonadetes (phylum Gemmatimonadetes) were also depleted by straw return in SZ and XZ soils.
Figure 4. Straw return-responsive taxa (P < 0.05) of fungi and bacteria in each stage, site and year. The color of the cell indicates the log2 fold change in relative abundance with respect to the control treatment: an increase tends toward red, while a decrease tends toward green. The short dashes in the heatmap indicate the class (-), order (–), family (—), and genus (—-), the red fonts without short dash indicate the phylum level. Fungal genera marked in blue fonts. The plot displays all taxa detected as significantly affected by straw treatment in at least one site or year. “U” represents the relative abundance significantly increased in R soils, and “D” represents the relative abundance significantly decreased in R soils.
Based on the LDA effect size (LEfSe) method, we identified several bacterial, and fungal genera which were significantly associated with wheat straw return. The bacterial genera Rhizobium, Devosia, Altererythrobacter, and Chitinophaga were more abundant in R soils, and RB41, Microcoleus, Desulfurivibrio, and Elioraea were more common in N soils (Figure 5A). The fungal genera Mortierella, Pyrenochaetopsis, Pyrenophora, Aspergillus, Scedosporium, Phialocephala, Stropharia, Myrmecridium, Crocicreas, Dactylella, Arthrographis, and Pluteus were more abundant in R soils, and Fusarium, Crinipellis, Bjerkandera, Amauroascus, Myceliophthora, Lepidosphaeria, Myrothecium, and Magnaporthiopsis were more common in N soils (Figure 5B). Interestingly, the bacterial genus Rhizobium was related to nitrogen cycling (Figure 5A), and the fungal genera Fusarium and Myrothecium contain pathogens of soybean or wheat (Figure 5B).
Figure 5. Histogram of the LDA scores computed for differentially abundant bacterial (A) and fungal (B) genera between the N and R soils.
Microbes Involved in the Nitrogen Cycling
We further identified 19 genera associated with nitrogen cycling from the sequencing results, according to the description of previous reports (Nelson et al., 2016 and Supplementary Table S2). In most of the cases (72%), the total abundance of these nitrogen cycle-related bacteria was higher in R soils than in N soils, and among the three sites, the overall increasing rates in R soils were higher than in N soils during the 3 years, particular at 60 days (Figures 6A,B). Among these bacteria, the total abundance of nitrogen-fixing bacteria (Bradyrhizobium, Rhizobium, and Burkholderia) and denitrifying bacteria (Hyphomicrobium) was significantly increased by straw return. For example, Rhizobium was significantly enriched in R soils of Jining and Suzhou, and Bradyrhizobium was significantly enriched in R soils of Suzhou (Figures 6A,C,D).
Figure 6. OTUs involved in the nitrogen cycle in N and R soils from different sites. (A) Detailed relative abundances of 19 nitrogen cycle-related genera in indicated samples. Asterisk indicates a significant difference with P < 0.05 (Wilcoxon rank-sum test). (B) Accumulated relative abundance of 19 nitrogen cycle-related genera in N and R soils during the three-season experiment. (C,D) Relative abundance of the genera Rhizobium and Bradyrhizobium, respectively. Asterisk indicates a significant difference with P < 0.05 (Wilcoxon rank-sum test).
Microbes Related to Wheat and/or Soybean Pathogens
Comparing to the R soils, the N soils exhibited higher abundances of fungal genera that contain wheat and/or soybean pathogens, such as Fusarium, Alternaria, and Myrothecium (Bai et al., 2015; Zeng et al., 2016; Dewangan et al., 2017). In our amplicon sequencing annotation results, we identified 23 fungal genera and one bacterial genus containing plant pathogens. The 23 fungal genera were of the phyla Ascomycota and Basidiomycota. The relative abundances of these genera were low (<1%), with the exceptions of Fusarium, Alternaria, and Mycosphaerella (Figure 7A). Surprisingly, straw return decreased the total abundance of pathogenic fungi (Figures 7A,B) but increased that of the pathogenic bacterium Xanthomonadales (Figure 7C). In detail, the relative abundances of the pathogenic fungal genera Fusarium, Alternaria, and Myrothecium were significantly lower in R soils, while that of Mycosphaerella, which was significantly higher (P < 0.05, Wilcoxon rank-sum test) (Figures 7A,D–G).
Figure 7. OTUs involved in wheat or soybean pathogen-associated genera in N and R soils from different sites. (A) Detailed relative abundance of 23 fungal genera in indicated sample. Asterisk indicates a significant difference with P < 0.05 (Wilcoxon rank-sum test). (B) Accumulated relative abundance of 23 fungal genera in N and R soils during the three-season experiment. (C–G) Relative abundance of the bacteria order Xanthomonadales and the fungal genera Fusarium, Alternaria, Myrothecium, and Mycosphaerella, respectively. Asterisk indicates a significant difference with P < 0.05 (Wilcoxon rank-sum test).
Effects of Soil Total Organic Carbon and Total Nitrogen on Microbial Communities
Mantel tests revealed that contents of TOC and TN in soil were significantly correlated with the composition of the soil bacterial and fungal communities (R = 0.61, P < 0.01 for bacteria and R = 0.57, P < 0.01 for fungi). According to Spearman’s rank-correlation coefficients, the bacterial richness and diversity, and the fungal diversity but not richness, were significantly correlated with soil TOC and TN contents (P < 0.01). TOC and TN contents were slightly different between R than N soils; however, some results indicated that there might be tendency that TN content was increased after wheat straw return. For example, TN content was significantly higher in R than N soils of SZ (0 day) and XZ (60 days) collected in 2017 (P < 0.05, Wilcoxon rank-sum test; Figures 8A,B).
Figure 8. Total nitrogen (A) and total organic carbon (B) contents in the indicated soil samples, and soybean yields between no straw return and straw return fields (C). Significant differences between the treatments are indicated by arrows (P < 0.05, Wilcoxon rank-sum test).
Discussion
Despite its widespread applications to improve soil fertility and crop yields, the effects of wheat straw return on soil microbial communities in a wheat–soybean rotation system was unclear. Our results of investigation in the Huang-Huai region showed that wheat straw return altered the taxonomic and functional traits of soil microbial communities, including increases in the relative abundances of nitrogen cycle-related bacterial taxa and decreases in those of some pathogenic fungal taxa. These shifts are likely to impact belowground ecosystems.
To evaluate the consistent effects of straw return on soil microbial communities, we collected samples from three different fields. Comparing to the effects of straw return, the effects of field location was stronger. The bacterial and fungal α-diversity differed among the three sites, and the primary source of β-diversity was the field location (Figures 1, 2). This is expected and probably due to the differences of soil type and environmental factors (Brockett et al., 2012; Peiffer et al., 2013; Edwards et al., 2015; Almario et al., 2017). In addition, the response of the α-diversity of the soil bacterial community to wheat straw return differed among the three sites (Supplementary Figure S1), probably due to the influence of soil texture and agronomic management history (Cookson et al., 2008; Bach et al., 2010; Xiong et al., 2015; Zhang et al., 2017). Furthermore, planting year and sampling period affected the soil bacterial and fungal communities (Figures 1, 2). These results are in line with prior reports (Sugiyama et al., 2014; Zeng et al., 2014; Bai et al., 2015; Shen et al., 2015; Dang et al., 2017).
Based on the samples collected from three different fields across three seasons, we found that wheat straw return affected the relative abundance and diversity of the overall soil bacterial community (Figures 1A,B). Although there was no significant overall effect on fungal community (Figures 1C,D), stronger responses of fungal communities to wheat straw return were observed, especially in the second, and third years (Figures 2C–H and Supplementary Table S1). This may because that the life generation cycle of bacteria is more rapid than fungi, and bacteria dominate in the initial phases of crop straw decomposition and fungi dominate the latter stage (Poll et al., 2008; Marschner et al., 2011; Chen Z. et al., 2017). Furthermore, soil bacteria dominated the decomposition of crop residue in agricultural soils, and fungi were the main decomposers in prairie soils (Allison et al., 2005). The short-term responses of soil microbial communities to wheat straw return were relatively weak compared to the other factors, they became evident after excluding the variation caused by field location and year. For example, wheat straw return enhanced the relative abundance of the bacterial phylum Proteobacteria and the fungal phylum Zygomycota, but reduced that of the bacterial phylum Acidobacteria and the fungal phylum Ascomycota (Figures 3A,B). The Proteobacteria to Acidobacteria ratio has been suggested to be an indicator of the trophic level of soils (Smit et al., 2001; Castro et al., 2010; Gottel et al., 2011; Beckers et al., 2017), as Proteobacteria are favored in nutrient-rich soils and Acidobacteria in nutrient-poor soils. Alphaproteobacteria genera such as Rhizobium, Devosia, and Altererythrobacter were more abundant in R soils, and RB41 (Acidobacteria) were more common in N soils (Figure 5A). Ascomycota was more abundant in N soils, whereas Zygomycota was more common in R soils. Ascomycota includes various soil-borne plant-pathogenic fungi as well as non-pathogenic genera, and its relative abundance is higher in areas of low productivity (Chang et al., 2017). Among them, the relative abundances of Fusarium were significantly lower in R soils. Zygomycota is reportedly rare in agricultural soils (Huang et al., 2015; Wang et al., 2015). However, dominant Zygomycota species such as Mortierella may play a role in disease suppression (Tagawa et al., 2010; Wani et al., 2017; Xiong et al., 2017). In our results, Mortierella was the most abundant genus in R soils. These results indicated that wheat straw return could be beneficial to soil health.
Wheat straw return enriched bacteria genera related to the soil nitrogen cycle, including Nitrospira, Bradyrhizobium, Rhizobium, and Hyphomicrobium (Figure 6A). Nitrospira were thought to be capable only of nitrite oxidation but recent evidence suggests that they may oxidize ammonia (Daims et al., 2015; Longa et al., 2017; Kuypers et al., 2018). Bradyrhizobium and Rhizobium are related to nodulation, and their symbiosis with legumes enhances nitrogen fixation and increases crop productivity (Sugiyama et al., 2014; Chang et al., 2017; Kuypers et al., 2018). Hyphomicrobium has denitrifying activity and is widely used in wastewater treatment (Rissanen et al., 2016; Sun et al., 2017; Yasuda et al., 2017). The relative abundances of Bradyrhizobium and Hyphomicrobium have been reportedly to be higher in a straw-returned than control (no straw return) soil in the North China Plain (Chen Y. et al., 2017). In additional to microbial community traits, soil characteristics may also be affected by wheat straw return. In this study, TN contents were moderately higher in R than N soils (Figure 8A). This is in line with a previous report that wheat straw return significantly increased the TN content compared to the control treatment (Chen Y. et al., 2017). Although we found that wheat straw return did not significantly influence the soil TOC content (Figure 8B), we speculated that impacts of straw return on soil C and N contents are strongly dependent on initial soil C and N, the rates and duration of straw return, and local soil and climate conditions.
Wheat straw return significantly reduced the abundances of plant pathogen associated fungal genera such as Fusarium, Alternaria, and Myrothecium (Figures 7A,D,E,F). Some Fusarium species can cause Fusarium root rot of soybean and Fusarium head blight of wheat both diseases markedly decrease crop productivity (Bai et al., 2015; Cainong et al., 2015; Lenc, 2015; Lu et al., 2015). Similarly, some Alternaria species can cause leaf-spot disease of soybean, sunflower, Withania somnifera, and Pelargonium (Furukawa and Kishi, 2001; Oliveira et al., 2004; Pati et al., 2008; Kumar and Kumar, 2017). Plant-pathogenic Myrothecium species cause foliar spots in soybean, cotton, Anthurium, Nicandra physaloides, tomato, and cucumber (Quezado Duval et al., 2010; Kwon et al., 2014; Dewangan et al., 2017). However, the relative abundances of other microbial taxa that include some plant pathogens, such as Mycosphaerella, were enhanced in R soils (Figure 7G). Mycosphaerella species can cause leaf-spot disease in economically important crops including wheat, soybean, banana, citrus, and eucalyptus (Stukenbrock et al., 2011; Ohm et al., 2012; Zeng et al., 2017). In addition, the relative abundance of the bacterial order Xanthomonadales was increased in R soils (Figure 7C). Xanthomonadales species cause major losses of tomato, cabbage, pepper, banana, citrus, rice, grape, peach, plum, almond, coffee, and maple crops (Ryan et al., 2011; Naushad and Gupta, 2013; Garita-Cambronero et al., 2018).
Straw return is often considered to increase pathogen loads and promote diseases of crops (Li et al., 2012; Qiao et al., 2013). Yet, there are other examples where long-term straw return suppressed some plant diseases and promote healthy plant growth (Osunlaja, 1990; Bailey and Lazarovits, 2003). Tillage and straw incorporation into soil may reduce pathogen inocula, while facilitating nutrient cycling, probably because no-tillage, and straw-returning practices improving soil nutrition condition for soil microbial communities (Guo et al., 2016). In our 3-yr study, as an indication of plant healthy, crop yield in R soils did not differ significantly than that in N soils every year (Figure 8C). The lower pathogenic fungal diversity in our straw return soil may be to the result of the lower returning amount (ca.5000 kg.hm–2) of wheat straw before soybean production. Some previous studies observed that the occurrence of soil-borne diseases was significantly different under different straw amendment rates. For example, the indexes of wheat soil-borne diseases were reduced significantly at the maize straw amendment rates of 7500 and 3750 kg.hm–2), but increased dramatically when the amendment rate increased to 15000 kg/hm2 (Zhen et al., 2009). In our study, although the abundance of pathogenic fungi in R soil is lower than that in N, abundances of many fungal genera in R is still very high, and straw return could increase soil moisture, which may further promote the occurrence of diseases. Future studies using a series of root-associated compartments experiments would be helpful for disentangling the relationships between soil microbiota and straw return.
Conclusion
Our results from a short-term, multisite experiment showed that wheat straw return could alter the taxonomic, and functional traits of the soil bacterial and fungal communities. The bacterial phylum Proteobacteria and the fungal phylum Zygomycota were enriched in R soils, but the bacterial phylum Acidobacteria, and the fungal phylum Ascomycota were depleted. Wheat straw return enriched the population of bacteria related to the soil nitrogen cycle, such as Bradyrhizobium and Rhizobium. Wheat straw return significantly depleted plant pathogen associated fungal taxa, such as Fusarium and Alternaria. R soils had moderately higher TN contents than N soils. Long-term experiments are needed to examine whether our observed changes in the microbial community composition induced by wheat straw return persist over time.
Data Availability
The datasets generated for this study can be found in NCBI, BioProject: PRJNA524011.
Author Contributions
HY performed the experiments, analyzed the data, contributed reagents, materials, and analysis tools, prepared the figures and/or tables, and drafted the manuscript and approved its final version. JM performed the experiments, contributed reagents, materials, and analysis tools, and approved the final version of the manuscript. ZR and DZ contributed reagents, materials, and analysis tools, and approved the final version of the manuscript. YW conceived and designed the experiments, and reviewed the drafted manuscript and approved its final version. SH reviewed the drafted manuscript and approved its final version. WY conceived and designed the experiments, prepared the figures and/or tables, and reviewed the drafted manuscript and approved its final version. XZ conceived and designed the experiments and approved the final version of the manuscript.
Funding
This work was supported by grants to XZ, YW, and WY from the National Key R&D Program of China (2018YFD0201000), the Special Fund for Agro-scientific Research in the Public Interest (201303018), the China Agriculture Research System (CARS-004-PS14), and the National Natural Science Foundation of China (31721004 and 31772140).
Conflict of Interest Statement
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.
Acknowledgments
We thank Jining Academy of Agricultural Sciences, Suzhou Academy of Agricultural Sciences, and Xuzhou Academy of Agricultural Sciences for providing the site access and allowing us to harvest the samples from their field trial. We also thank Prof. Cunxiang Wu in the Chinese Academy of Agricultural Sciences, Prof. Haijian Zhi, Zhengguang Zhang, and Suomeng Dong in Nanjing Agricultural University, and Dr. Wei Zhao in the Anhui Academy of Agricultural Sciences for suggestion on experiment design.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.01811/full#supplementary-material
FIGURE S1 | Field experiment design. (A) The geographic coordinates of the three selected sites in the Huang-Huai region of China. (B) Schematic representation of field experiment layout. Each site had two treatments. Each treatment contained five replicated subplots.
FIGURE S2 | Bacterial and fungal alpha diversity estimated by Shannon’s diversity index. Significant differences between the treatments are distinguished by different letters (P < 0.05, Wilcoxon rank-sum test).
TABLE S1 | Influence of straw return on microbial community composition according to R-values derived from an ANOSIM. R values are shown with *P < 0.05, ∗∗P < 0.01, and ∗∗∗P < 0.001.
TABLE S2 | Annotation of the 19 nitrogen cycling related bacterial genera according to the description of previous reports.
Footnotes
- ^ https://www.mothur.org
- ^ https://www.r-project.org
- ^ http://qiime.org
- ^ http://huttenhower.sph.harvard.edu/galaxy
References
Allison, V. J., Miller, R. M., Jastrow, J. D., Matamala, R., and Zak, D. R. (2005). Changes in soil microbial community structure in a tallgrass prairie chronosequence. Soil Sci. Soc. Am. J. 69, 1412–1421. doi: 10.2136/sssaj2004.0252
Almario, J., Jeena, G., Wunder, J., Langen, G., Zuccaro, A., Coupland, G., et al. (2017). Root-associated fungal microbiota of nonmycorrhizal Arabis alpina and its contribution to plant phosphorus nutrition. Proc. Natl. Acad. Sci. U.S.A. 114, E9403–E9412. doi: 10.1073/pnas.1710455114
Bach, E. M., Baer, S. G., Meyer, C. K., and Six, J. (2010). Soil texture affects soil microbial and structural recovery during grassland restoration. Soil Biol. Biochem. 42, 2182–2191. doi: 10.1016/j.soilbio.2010.08.014
Bai, L., Cui, J., Jie, W., and Cai, B. (2015). Analysis of the community compositions of rhizosphere fungi in soybeans continuous cropping fields. Microbiol. Res. 180, 49–56. doi: 10.1016/j.micres.2015.07.007
Bailey, K. L., and Lazarovits, G. (2003). Suppressing soil-borne diseases with residue management and organic amendments. Soil Tillage Res. 72, 169–180. doi: 10.1016/s0167-1987(03)00086-2
Bao, Q., Huang, Y., Wang, F., Nie, S., Nicol, G. W., Yao, H., et al. (2016). Effect of nitrogen fertilizer and/or rice straw amendment on methanogenic archaeal communities and methane production from a rice paddy soil. Appl. Microbiol. Biotechnol. 100, 5989–5998. doi: 10.1007/s00253-016-7377-z
Beckers, B., Beeck, M. O. D., Weyens, N., Boerjan, W., and Vangronsveld, J. (2017). Structural variability and niche differentiation in the rhizosphere and endosphere bacterial microbiome of field-grown poplar trees. Microbiome 5:25. doi: 10.1186/s40168-017-0241-2
Brockett, B. F. T., Prescott, C. E., and Grayston, S. J. (2012). Soil moisture is the major factor influencing microbial community structure and enzyme activities across seven biogeoclimatic zones in western Canada. Soil Biol. Biochem. 44, 9–20. doi: 10.1016/j.soilbio.2011.09.003
Cainong, J. C., Bockus, W. W., Feng, Y., Chen, P., Qi, L., Sehgal, S. K., et al. (2015). Chromosome engineering, mapping, and transferring of resistance to Fusarium head blight disease from Elymus tsukushiensis into wheat. Theor. Appl. Genet. 128, 1019–1027. doi: 10.1007/s00122-015-2485-1
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Caporaso, J. G., Lauber, C. L., Walters, W. A., Berg-Lyons, D., Lozupone, C. A., Turnbaugh, P. J., et al. (2011). Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. U.S.A. 108, 4516–4522. doi: 10.1073/pnas.1000080107
Castro, H. F., Classen, A. T., Austin, E. E., Norby, R. J., and Schadt, C. W. (2010). Soil microbial community responses to multiple experimental climate change drivers. Appl. Environ. Microbiol. 76, 999–1007. doi: 10.1128/AEM.02874-9
Chang, H. X., Haudenshield, J. S., Bowen, C. R., and Hartman, G. L. (2017). Metagenome-wide association study and machine learning prediction of bulk soil microbiome and crop productivity. Front. Microbiol. 8:519. doi: 10.3389/fmicb.2017.00519
Chen, Y., Xin, L., Liu, J., Yuan, M., Liu, S., Jiang, W., et al. (2017). Changes in bacterial community of soil induced by long-term straw returning. Sci. Agric. 74, 349–356. doi: 10.1590/1678-992x-2016-2025
Chen, Z., Wang, H., Liu, X., Zhao, X., Lu, D., Zhou, J., et al. (2017). Changes in soil microbial community and organic carbon fractions under short-term straw return in a rice–wheat cropping system. Soil Tillage Res. 165, 121–127. doi: 10.1016/j.still.2016.07.018
Cookson, W. R., Murphy, D. V., and Roper, M. M. (2008). Characterizing the relationships between soil organic matter components and microbial function and composition along a tillage disturbance gradient. Soil Biol. Biochem. 40, 763–777. doi: 10.1016/j.soilbio.2007.10.011
Daims, H., Lebedeva, E. V., Pjevac, P., Han, P., Herbold, C., Albertsen, M., et al. (2015). Complete nitrification by Nitrospira bacteria. Nature 528, 504–509. doi: 10.1038/nature16461
Dang, P., Yu, X., Le, H., Liu, J., Shen, Z., and Zhao, Z. (2017). Effects of stand age and soil properties on soil bacterial and fungal community composition in Chinese pine plantations on the Loess Plateau. PLoS One 12:e0186501. doi: 10.1371/journal.pone.0186501
Dewangan, M., Kurre, A. K., Pandey, S. K., and Nag, U. K. (2017). Evaluation of fungicides against Myrothecium roridum causing leaf spot of soybean under in vitro conditions. Int. J. Curr. Microbiol. Appl. Sci. 6, 1953–1958. doi: 10.20546/ijcmas.2017.608.231
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381
Edwards, J., Johnson, C., Santos-Medellin, C., Lurie, E., Podishetty, N. K., Bhatnagar, S., et al. (2015). Structure, variation, and assembly of the root-associated microbiomes of rice. Proc. Natl. Acad. Sci. U.S.A. 112, E911–E920. doi: 10.1073/pnas.1414592112
Furukawa, T., and Kishi, K. (2001). Alternaria leaf spot on three species of Pelargoniurn caused by Alternaria alternata in Japan.pdf. J. Gen. Plant Pathol. 67, 268–272. doi: 10.1007/pl00013028
Garita-Cambronero, J., Palacio-Bielsa, A., and Cubero, J. (2018). Xanthomonas arboricola pv. Pruni, causal agent of bacterial spot of stone fruits and almond: its genomic and phenotypic characteristics in the X. Arboricola species context. Mol. Plant Pathol. 19, 2053–2065. doi: 10.1111/mpp.12679
Gottel, N. R., Castro, H. F., Kerley, M., Yang, Z., Pelletier, D. A., Podar, M., et al. (2011). Distinct microbial communities within the endosphere and rhizosphere of Populus deltoides roots across contrasting soil types. Appl. Environ. Microbiol. 77, 5934–5944. doi: 10.1128/AEM.05255-11
Guo, L., Zheng, S., Cao, C., and Li, C. (2016). Tillage practices and straw-returning methods affect topsoil bacterial community and organic C under a rice-wheat cropping system in central China. Sci. Rep. 6:33155. doi: 10.1038/srep33155
Hansen, V., Muller-Stover, D., Imparato, V., Krogh, P. H., Jensen, L. S., Dolmer, A., et al. (2017). The effects of straw or straw-derived gasification biochar applications on soil quality and crop productivity: a farm case study. J. Environ. Manage 186(Pt 1), 88–95. doi: 10.1016/j.jenvman.2016.10.041
Huang, X., Liu, L., Wen, T., Zhu, R., Zhang, J., and Cai, Z. (2015). Illumina MiSeq investigations on the changes of microbial community in the Fusarium oxysporum f.sp. cubense infected soil during and after reductive soil disinfestation. Microbiol. Res. 181, 33–42. doi: 10.1016/j.micres.2015.08.004
Ji, B., Hu, H., Zhao, Y., Mu, X., Liu, K., and Li, C. (2014). Effects of deep tillage and straw returning on soil microorganism and enzyme activities. Sci. World J. 2014:451493. doi: 10.1155/2014/451493
Karlsson, H., Ahlgren, S., Sandgren, M., Passoth, V., Wallberg, O., and Hansson, P. A. (2017). Greenhouse gas performance of biochemical biodiesel production from straw: soil organic carbon changes and time-dependent climate impact. Biotechnol. Biofuels 10:217. doi: 10.1186/s13068-017-0907-9
Koljalg, U., Nilsson, R. H., Abarenkov, K., Tedersoo, L., Taylor, A. F. S., Bahram, M., et al. (2013). Towards a unified paradigm for sequence-based identification of fungi. Mol. Ecol. 22, 5271–5277. doi: 10.1111/mec.12481
Kumar, S., and Kumar, P. (2017). Effect of date of sowing and varieties on severity of Alternaria leaf spot, rust and yield in soybean (Glycine max L. Merrill). Bull. Env. Pharmacol. Life Sci. 6, 448–451.
Kuypers, M. M. M., Marchant, H. K., and Kartal, B. (2018). The microbial nitrogen-cycling network. Nat. Rev. Microbiol. 16, 263–276. doi: 10.1038/nrmicro.2018.9
Kwon, H. W., Kim, J. Y., Choi, M. A., Son, S. Y., and Kim, S. H. (2014). Characterization of myrothecium roridum isolated from imported anthurium plant culture medium. Mycobiology 42, 82–85. doi: 10.5941/MYCO.2014.42.1.82
Lenc, L. (2015). Fusarium head blight (FHB) and Fusarium populations in grain of winter wheat grown in different cultivation systems. J. Plant Protect. Res. 55, 94–109. doi: 10.1515/jppr-2015-0013
Lenka, N. K., and Lal, R. (2013). Soil aggregation and greenhouse gas flux after 15 years of wheat straw and fertilizer management in a no-till system. Soil Tillage Res. 126, 78–89. doi: 10.1016/j.still.2012.08.011
Li, H., Dai, M., Dai, S., and Dong, X. (2018). Current status and environment impact of direct straw return in China’s cropland - A review. Ecotoxicol. Environ. Saf. 159, 293–300. doi: 10.1016/j.ecoenv.2018.05.014
Li, H., Liu, F., Gong, Z., Jing, X., and Huang, S. (2012). Effect of straw return back to paddy field on occurrence of rice major diseases. Crop Res. 26:4.
Li, P., Li, Y., Zheng, X., Ding, L., Ming, F., Pan, A., et al. (2017). Rice straw decomposition affects diversity and dynamics of soil fungal community, but not bacteria. J. Soils Sediments 18, 248–258. doi: 10.1007/s11368-017-1749-6
Longa, C. M. O., Nicola, L., Antonielli, L., Mescalchin, E., Zanzotti, R., Turco, E., et al. (2017). Soil microbiota respond to green manure in organic vineyards. J. Appl. Microbiol. 123, 1547–1560. doi: 10.1111/jam.13606
Lu, C., Zhang, H., Wang, Y., and Zheng, X. (2015). Rapid diagnosis of Fusarium root rot in soybean caused by Fusarium equiseti or Fusarium graminearum using loop-mediated isothermal amplification (LAMP) assays. Australas. Plant Pathol. 44, 437–443. doi: 10.1007/s13313-015-0361-8
Lu, W., Yu, S., Ma, Y., and Huang, H. (2018). Integrated economic and environmental analysis of agricultural straw reuse in edible fungi industry. PeerJ 6:e4624. doi: 10.7717/peerj.4624
Ma, S., Fang, C., Sun, X., Han, L., He, X., and Huang, G. (2018). Bacterial community succession during pig manure and wheat straw aerobic composting covered with a semi-permeable membrane under slight positive pressure. Bioresour. Technol. 259, 221–227. doi: 10.1016/j.biortech.2018.03.054
Maarastawi, S. A., Frindte, K., Linnartz, M., and Knief, C. (2018). Crop rotation and straw application impact microbial communities in Italian and Philippine soils and the rhizosphere of Zea mays. Front. Microbiol. 9:1295. doi: 10.3389/fmicb.2018.01295
Magoc, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507
Marschner, P., Umar, S., and Baumann, K. (2011). The microbial community composition changes rapidly in the early stages of decomposition of wheat residue. Soil Biol. Biochem. 43, 445–451. doi: 10.1016/j.soilbio.2010.11.015
McGuire, K. L., Payne, S. G., Palmer, M. I., Gillikin, C. M., Keefe, D., Kim, S. J., et al. (2013). Digging the New York city skyline: soil fungal communities in green roofs and city parks. PLoS One 8:e58020. doi: 10.1371/journal.pone.0058020
Naushad, H. S., and Gupta, R. S. (2013). Phylogenomics and molecular signatures for species from the plant pathogen-containing order xanthomonadales. PLoS One 8:e55216. doi: 10.1371/journal.pone.0055216
Nelson, M. B., Martiny, A. C., and Martiny, J. B. (2016). Global biogeography of microbial nitrogen-cycling traits in soil. Proc. Natl. Acad. Sci. U.S.A. 113, 8033–8040. doi: 10.1073/pnas.1601070113
Ohm, R. A., Feau, N., Henrissat, B., Schoch, C. L., Horwitz, B. A., Barry, K. W., et al. (2012). Diverse lifestyles and strategies of plant pathogenesis encoded in the genomes of eighteen Dothideomycetes fungi. PLoS Pathog. 8:e1003037. doi: 10.1371/journal.ppat.1003037
Oliveira, M. F. D., Neto, A. T., Leite, R. M. V. B. C., Castiglioni, V. B. R., and Arias, C. A. A. (2004). Mutation breeding in sunflower for resistance to alternaria leaf spot. Helia 27, 41–50. doi: 10.2298/HEL0441041d
Osunlaja, S. O. (1990). Effect of organic soil amendments on the incidence of stalk rot of maize. Plant Soil 127, 237–241. doi: 10.1007/bf00014431
Pati, P. K., Sharma, M., Salar, R. K., Sharma, A., Gupta, A. P., and Singh, B. (2008). Studies on leaf spot disease of Withania somnifera and its impact on secondary metabolites. Indian J. Microbiol. 48, 432–437. doi: 10.1007/s12088-008-0053-y
Peiffer, J. A., Spor, A., Koren, O., Jin, Z., Tringe, S. G., Dangl, J. L., et al. (2013). Diversity and heritability of the maize rhizosphere microbiome under field conditions. Proc. Natl. Acad. Sci. U.S.A. 110, 6548–6553. doi: 10.1073/pnas.1302837110
Peng, S., Guo, T., and Liu, G. (2013). The effects of arbuscular mycorrhizal hyphal networks on soil aggregations of purple soil in southwest China. Soil Biol. Biochem. 57, 411–417. doi: 10.1016/j.soilbio.2012.10.026
Poll, C., Marhan, S., Ingwersen, J., and Kandeler, E. (2008). Dynamics of litter carbon turnover and microbial abundance in a rye detritusphere. Soil Biol. Biochem. 40, 1306–1321. doi: 10.1016/j.soilbio.2007.04.002
Qiao, Y., Cao, C., Zhao, Z., Du, S., Zhang, Y., Liu, Y., et al. (2013). Effects of straw returning and N fertilizer application on yield, quality and occurrence of Fusarium head blight of Wheat. J. Triticeae Crops 33, 727–731. doi: 10.7606/j.issn.1009-1041.2013.04.018
Qu, C., Li, B., Wu, H., and Giesy, J. P. (2012). Controlling air pollution from straw burning in China calls for efficient recycling. Environ. Sci. Technol. 46, 7934–7936. doi: 10.1021/es302666s
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
Quezado Duval, A. M., Henz, G. P., Paz-Lima, M. L., Medeiros, A. R., Miranda, B. E. C., Pfenning, L. H., et al. (2010). New hosts of Myrothecium SPP. In brazil and a preliminary In Vitro assay of fungicides. Braz. J. Microbiol. 41, 246–252.
Rissanen, A. J., Ojala, A., Fred, T., Toivonen, J., and Tiirola, M. (2016). Methylophilaceae and Hyphomicrobium as target taxonomic groups in monitoring the function of methanol-fed denitrification biofilters in municipal wastewater treatment plants. J. Ind. Microbiol. Biotechnol. 44, 35–47. doi: 10.1007/s10295-016-1860-5
Ryan, R. P., Vorholter, F. J., Potnis, N., Jones, J. B., Van Sluys, M. A., Bogdanove, A. J., et al. (2011). Pathogenomics of Xanthomonas: understanding bacterium-plant interactions. Nat. Rev. Microbiol. 9, 344–355. doi: 10.1038/nrmicro2558
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60
Shen, Z., Ruan, Y., Wang, B., Zhong, S., Su, L., Li, R., et al. (2015). Effect of biofertilizer for suppressing Fusarium wilt disease of banana as well as enhancing microbial and chemical properties of soil under greenhouse trial. Appl. Soil Ecol. 93, 111–119. doi: 10.1016/j.apsoil.2015.04.013
Smit, E., Leeflang, P., Gommans, S., van den Broek, J., van Mil, S., and Wernars, K. (2001). Diversity and seasonal fluctuations of the dominant members of the bacterial soil community in a wheat field as determined by cultivation and molecular methods. Appl. Environ. Microbiol. 67, 2284–2291. doi: 10.1128/AEM.67.5.2284-2291.2001
Stark, S., Männistö, M. K., Ganzert, L., Tiirola, M., and Häggblom, M. M. (2015). Grazing intensity in subarctic tundra affects the temperature adaptation of soil microbial communities. Soil Biol. Biochem. 84, 147–157. doi: 10.1016/j.soilbio.2015.02.023
Stukenbrock, E. H., Bataillon, T., Dutheil, J. Y., Hansen, T. T., Li, R., Zala, M., et al. (2011). The making of a new pathogen: insights from comparative population genomics of the domesticated wheat pathogen Mycosphaerella graminicola and its wild sister species. Genome Res. 21, 2157–2166. doi: 10.1101/gr.118851.110
Sugiyama, A., Ueda, Y., Zushi, T., Takase, H., and Yazaki, K. (2014). Changes in the bacterial community of soybean rhizospheres during growth in the field. PLoS One 9:e100709. doi: 10.1371/journal.pone.0100709
Sun, H., Liu, F., Xu, S., Wu, S., Zhuang, G., Deng, Y., et al. (2017). Myriophyllum aquaticum constructed wetland effectively removes nitrogen in swine wastewater. Front. Microbiol. 8:1932. doi: 10.3389/fmicb.2017.01932
Tagawa, M., Tamaki, H., Manome, A., Koyama, O., and Kamagata, Y. (2010). Isolation and characterization of antagonistic fungi against potato scab pathogens from potato field soils. FEMS Microbiol. Lett. 305, 136–142. doi: 10.1111/j.1574-6968.2010.01928.x
Valdez-Vazquez, I., Torres-Aguirre, G. J., Molina, C., and Ruiz-Aguilar, G. M. L. (2016). Characterization of a lignocellulolytic consortium and methane production from untreated wheat straw: dependence on nitrogen and phosphorous content. BioResources 11, 4237–4251. doi: 10.15376/biores.11.2.4237-4251
Wang, B., Li, R., Ruan, Y., Ou, Y., Zhao, Y., and Shen, Q. (2015). Pineapple-banana rotation reduced the amount of Fusarium oxysporum more than maize-banana rotation mainly through modulating fungal communities. Soil Biol. Biochem. 86, 77–86. doi: 10.1016/j.soilbio.2015.02.021
Wani, Z. A., Kumar, A., Sultan, P., Bindu, K., Riyaz-Ul-Hassan, S., and Ashraf, N. (2017). Mortierella alpina CS10E4, an oleaginous fungal endophyte of Crocus sativus L. enhances apocarotenoid biosynthesis and stress tolerance in the host plant. Sci. Rep. 7:8598. doi: 10.1038/s41598-017-08974-z
Xing, X., Zhou, Y., Lang, J., Chen, D., Cheng, S., Han, L., et al. (2018). Spatiotemporal variation of domestic biomass burning emissions in rural China based on a new estimation of fuel consumption. Sci. Total Environ. 626, 274–286. doi: 10.1016/j.scitotenv.2018.01.048
Xiong, W., Li, R., Ren, Y., Liu, C., Zhao, Q., Wu, H., et al. (2017). Distinct roles for soil fungal and bacterial communities associated with the suppression of vanilla Fusarium wilt disease. Soil Biol. Biochem. 107, 198–207. doi: 10.1016/j.soilbio.2017.01.010
Xiong, W., Li, Z., Liu, H., Xue, C., Zhang, R., Wu, H., et al. (2015). The Effect of long-term continuous cropping of black pepper on soil bacterial communities as determined by 454 pyrosequencing. PLoS One 10:e0136946. doi: 10.1371/journal.pone.0136946
Yang, H., Feng, J., Zhai, S., Dai, Y., Xu, M., Wu, J., et al. (2016). Long-term ditch-buried straw return alters soil water potential, temperature, and microbial communities in a rice-wheat rotation system. Soil Tillage Res. 163, 21–31. doi: 10.1016/j.still.2016.05.003
Yang, H., Ye, W., Ma, J., Zeng, D., Rong, Z., Xu, M., et al. (2018). Endophytic fungal communities associated with field-grown soybean roots and seeds in the Huang-Huai region of China. PeerJ 6:e4713. doi: 10.7717/peerj.4713
Yasuda, T., Waki, M., Fukumoto, Y., Hanajima, D., Kuroda, K., Suzuki, K., et al. (2017). Community structure of denitrifying and total bacteria during nitrogen accumulation in an ammonia-loaded biofilter. J. Appl. Microbiol. 123, 1498–1511. doi: 10.1111/jam.13603
Yu, D., Wen, Z., Li, X., Song, X., Wu, H., and Yang, P. (2018). Effects of straw return on bacterial communities in a wheat-maize rotation system in the North China Plain. PLoS One 13:e0198087. doi: 10.1371/journal.pone.0198087
Yuan, L., Zhang, Z., Cao, X., Zhu, S., Zhang, X., and Wu, L. (2014). Responses of rice production, milled rice quality and soil properties to various nitrogen inputs and rice straw incorporation under continuous plastic film mulching cultivation. Field Crops Res. 155, 164–171. doi: 10.1016/j.fcr.2013.09.009
Zeng, D. D., Zhang, H. F., Tian, Q., Xu, M., Wang, Y. C., and Zheng, X. B. (2016). Detection of soybean seed-borne pathogens in huang-huai area using LAMP assays. J. Nanj. Agric. Univ. 39, 947–953.
Zeng, F., Lian, X., Zhang, G., Yu, X., Bradley, C. A., and Ming, R. (2017). A comparative genome analysis of Cercospora sojina with other members of the pathogen genus Mycosphaerella on different plant hosts. Genom Data 13, 54–63. doi: 10.1016/j.gdata.2017.07.007
Zeng, H., Tan, F., Zhang, Y., Feng, Y., Shu, Y., and Wang, J. (2014). Effects of cultivation and return of Bacillus thuringiensis (Bt) maize on the diversity of the arbuscular mycorrhizal community in soils and roots of subsequently cultivated conventional maize. Soil Biol. Biochem. 75, 254–263. doi: 10.1016/j.soilbio.2014.04.024
Zhang, H., Xu, W., Li, Y., Lyu, J., Cao, Y., and He, W. (2017). Changes of soil microbial communities during decomposition of straw residues under different land uses. J. Arid Land 9, 666–677. doi: 10.1007/s40333-017-0027-1
Zhao, S., Li, K., Zhou, W., Qiu, S., Huang, S., and He, P. (2016). Changes in soil microbial community, enzyme activities and organic matter fractions under long-term straw return in north-central China. Agric. Ecosyst. Environ. 216, 82–88. doi: 10.1016/j.agee.2015.09.028
Zhen, W., Wang, S., Zhang, C., and Ma, Z. (2009). Influence of maize straw amendment on soil-borne diseases of winter wheat. Front. Agric. China 3:7. doi: 10.1007/s11703-009-0003-4
Keywords: wheat straw return, wheat–soybean rotation, soil bacterial and fungal community, 16S rRNA and ITS amplicon sequencing, nitrogen-cycling microbe, plant pathogen
Citation: Yang H, Ma J, Rong Z, Zeng D, Wang Y, Hu S, Ye W and Zheng X (2019) Wheat Straw Return Influences Nitrogen-Cycling and Pathogen Associated Soil Microbiota in a Wheat–Soybean Rotation System. Front. Microbiol. 10:1811. doi: 10.3389/fmicb.2019.01811
Received: 26 February 2019; Accepted: 23 July 2019;
Published: 08 August 2019.
Edited by:
Manuel Delgado Baquerizo, Universidad Rey Juan Carlos, SpainReviewed by:
Felipe Bastida, Centro de Edafología y Biología Aplicada del Segura (CSIC), SpainHangwei Hu, The University of Melbourne, Australia
Copyright © 2019 Yang, Ma, Rong, Zeng, Wang, Hu, Ye and Zheng. 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: Wenwu Ye, yeww@njau.edu.cn