- 1Ecosystems and Environment Research Programme, University of Helsinki, Helsinki, Finland
- 2Helsinki Institute of Sustainability Science (HELSUS), University of Helsinki, Helsinki, Finland
- 3Department of Microbiology, College of Resources, Sichuan Agricultural University, Chengdu, China
- 4Kemira, Kemira Oyj, Espoo R&D Center, Espoo, Finland
- 5Department of Agricultural Sciences and Viikki Plant Sciences Centre, Helsinki, Finland
Pastures are an important part of crop and food systems in cold climates. Understanding how fertilization and plant species affect soil bacterial community diversity and composition is the key for understanding the role of soil bacteria in sustainable agriculture. To study the response of soil bacteria to different fertilization and cropping managements, a 3-year (2013–2015) field study was established. In the split-plot design, fertilizer treatment (unfertilized control, organic fertilizer, and synthetic fertilizer) was the main plot factor, and plant treatment [clear fallow, red clover (Trifolium pratense), timothy (Phleum pratense), and a mixture of red clover and timothy] was the sub-plot factor. Soil bacterial community diversity and composition, soil properties, and crop growth were investigated through two growing seasons in 2014 and 2015, with different nitrogen input levels. The community diversity measures (richness, Shannon diversity, and Shannon evenness) and composition changed over time (P < 0.05) and at different time scales. The community diversity was lower in 2014 than in 2015. The temporal differences were greater than the differences between treatments. The overall correlations of Shannon diversity to soil pH, NO, NH, and surplus nitrogen were positive and that of bacterial richness to crop dry matter yield was negative (P < 0.05). The major differences in diversity and community composition were found between fallow and planted treatments and between organic and synthetic fertilizer treatments. The differences between the planted plots were restricted to individual operational taxonomic units (OTUs). Soil moisture, total carbon content, and total nitrogen content correlated consistently with the community composition (P < 0.05). Compared to the unfertilized control, the nitrogen fertilizer loading enhanced the temporal change of community composition in pure timothy and in the mixture more than that in red clover, which further emphasizes the complexity of interactions between fertilization and cropping treatments on soil bacteria.
Introduction
Grassy pastures are an important part of agricultural systems where ruminants are the main way for producing food under cold, dry, or other stressful conditions. Pastures for silage, hay, or grazing cover 26% of the arable region of Finland (Niemi and Väre, 2017), where the climate is defined as boreal and sub-boreal. To sustain or improve pasture production, farming systems have been widely intensified particularly through nitrogen (N) fertilizer inputs (Tilman et al., 2002; Stott and Gourley, 2016). However, reactive N as nitrate (NO) in either synthetic fertilizer or manure is highly mobile and easily leached into waterways or denitrified, resulting in substantial environmental pollution (Philippot et al., 2007; Ju et al., 2009). Soil microbes can directly or indirectly affect plant productivity by regulating the acquisition of limiting nutrients and the farming impacts on the environment by moderating the biogeochemical processes (Verwaal et al., 2006; Van et al., 2008; Strickland et al., 2009). Soil bacteria are critical for N flows through biological nitrogen fixation (BNF) that provides N for crop growth, nitrification by oxidizing NH into NO, denitrification by reducing NO into NO, N2O, and finally N2, and in dissimilatory NO reduction to ammonium (Philippot et al., 2007; Franche et al., 2009; Jia and Conrad, 2009; Yoshida et al., 2009; Robertson and Groffman, 2015).
Achieving or maintaining functional and diverse populations of soil organisms is one of the major criteria for sustainable land management (Dumanski et al., 1991). The effects of fertilizer and crop management on soil microbial communities may help us to understand the key roles of soil microbes on the biogeochemical processes of nutrient flows and crop production. Wang et al. (2017) reported a decrease in bacterial richness with synthetic fertilizer, and the simultaneous loss of soil biodiversity caused by intensification has raised concerns (Orgiazzi et al., 2016). In a comparison between grassland and managed fields, more intensive management resulted in greater differences in bacterial communities over time (Lauber et al., 2013). Generally, detecting the differences requires long-term studies (Shen et al., 2010; Hui et al., 2017). A meta-analysis by de Graaff et al. (2019) summarized that bacterial diversity was increased when synthetic N input rate was less than 150 kg ha–1 year–1, with organic N applied for more than 5 years, and the authors suggested that this was due to increased soil organic matter accumulation and retention. With elevated N input or across a N gradient, the community composition became progressively more distinct at higher N levels from those at no or lower levels of N fertilizer (Ramirez et al., 2010; Fierer et al., 2012). In addition, Soong et al. (2020) reported that soil organic carbon (C) content is the primary limitation for heterotrophic microbial communities.
The presence of N2-fixing legumes in legume–grass intercropping was considered as a component of sustainable intensification to enhance resource use efficiency and improve crop productivity (Tilman et al., 2002; Martin-Guay et al., 2018). Plant species could strongly affect the structure of rhizosphere microbial communities through root exudates (Haichar et al., 2008; Berg and Smalla, 2009). Nevertheless, little is known about the temporal response of soil bacterial community diversity and composition to the presence of legumes in field-grown crop mixtures with different fertilization treatments. Hence, we established a 3-year field study in the sub-boreal conditions of southern Finland. We have previously shown that the management and the season induced temporal changes in crop growth and soil properties, with the legume–grass mixture providing particular advantages over the pure stands of red clover (Trifolium pratense) and timothy (Phleum pretense) in terms of improved pasture yield and reduced potential N loss (Li et al., 2019). Given that the changes in soil organic C content, N availability, C/N ratio (Marschner et al., 2003; Ramirez et al., 2010), soil pH (Enwall et al., 2007; Shen et al., 2010), and soil moisture and temperature (Lauber et al., 2013) could drive the differences in microbial communities, we hypothesized that the response of community diversity measures (richness, Shannon diversity, and Shannon evenness) and composition would change over time and that the higher N input level in the third year would result in greater differences between treatments. Assuming that the advantages of the grass–legume mixture were linked to the differences in the biogeochemical processes of the nutrients caused by the bacterial community, we hypothesized that the soil bacterial communities in the mixture would be different from that in pure stands and that this would provide insight into the connections between soil bacteria and sustainable intensification of agriculture.
Materials and Methods
Field Establishment and Soil Sampling
The flat experimental field of 0.56 ha was at the Viikki Experimental Farm, University of Helsinki, Finland (60.227°N, 25.018°E, approximately 10 m above mean sea level). The soil texture was clay loam, with 32% clay, 36% silt, and 32% sand on average. The preceding crops were faba bean in 2010 with fertilizer N at 20 kg ha–1 year–1, followed by barley in 2011 and 2012 with fertilizer N at 60 kg ha–1 year–1. Prior to this experiment, soil pH was 6.41, total N content was 1.68 g kg –1, total C content was 25.3 g kg–1, NO3–N content was 5.42 mg kg–1, and NH4–N content was 4.51 mg kg–1.
The experiment started in May 2013 and ended in September 2015. The split-plot design included 48 plots (four replicates × three fertilizer treatments × four crop treatments). Each replicate block of 18 m × 8 m was divided into 6 m × 8 m main plots for the fertilizer treatments: no fertilizer (control), organic (urine in 2014 and manure slurry in 2015), and synthetic [Ca(NO3)2]. The main plot was divided into four 6 m × 2 m sub-plots for the crop treatments: fallow, pure red clover, pure timothy, and a mixture of red clover and timothy. The weeds in the fallow were removed with a rotary tiller. Barley was grown in the surrounding field in all 3 years, with N fertilization at 80 kg ha–1 year–1. Fertilizer was applied manually twice per year. In 2014, with reduced nitrogen input, organic fertilizer was spread at 35 and 20 kg N ha–1 and synthetic fertilizer was spread at 40 and 20 kg N ha–1 on 7 May and 8 July, respectively. In 2015, with normal nitrogen input, both the organic and the synthetic fertilizers were spread at 75 and 75 kg N ha–1 on 1 June and 16 July, respectively. Except for the small difference at the first time point, the N inputs were consistent in the organic and the synthetic treatments. The surplus N (kg ha–1) was calculated as fertilization N input – crop N yield.
The surface soil (0–20 cm) was sampled twice annually in 2014 and 2015 as 16 subsamples per sub-plot and mixed into one composite sample. The sampling dates were June 2014, September 2014, July 2015, and September 2015. In total, 191 soil samples were investigated. Field-moist soil samples were passed through a 5 mm sieve and preserved at −20°C until DNA isolation and dried before the biophysiochemical analysis. Soil NH and NO were extracted with 2 M KCl and measured with Lachat QuickChem 8000 (Lachat Instruments, Milwaukee, WI, United States), soil pH and EC were determined in a 1:2.5 (w/w) soil–water mixture, and total C and N were determined with TruSpec Elemental Determinator (LECO, United States). The details were reported in Li et al. (2019).
DNA Isolation and Amplicon Sequencing
Soil DNA was isolated from 0.25 g of soil with a Power Soil DNA Isolation Kit (MoBio, Carlsbad, CA, United States), following the manufacturer’s instructions. The quality of isolated DNA was checked with electrophoresis in 1% agarose gel, and the DNA was quantified fluorometrically using a PicoGreen dsDNA Quantification Reagent Kit (Molecular Probes, United States) in 96-well plates according to the manufacturer’s instructions. To amplify the V3 and V4 regions of the bacterial 16S rRNA gene, we used the primers 341F and 785R in duplicate 25 μl PCR reactions (Supplementary Table S1). After PCR, the amplicons from the duplicates were pooled and checked with electrophoresis in 1% agarose gel. The amplicons were sequenced using Illumina MiSeq sequencing (2 × 300 bp paired end) at the Institute of Biotechnology sequencing core facility, University of Helsinki.
Bioinformatics
The sequences were processed using Mothur v.1.39.5 (Schloss et al., 2009) with standard operating procedures SOP (Kozich et al., 2013) in the Taito supercluster of the Finnish IT Centre for Science (CSC, Espoo, Finland). In total, 11,020,235 sequences were aligned against the SILVA reference database v.132. After removing the chimeras, 2,568,584 sequences were classified into operational taxonomic units (OTUs) at 97% similarity level against the SILVA database v.132, with a bootstrap value of 80% for taxonomic assignment. Finally, 27,535 OTUs were obtained after removing the singletons and subsampled into the sequencing depth of 4,680 sequences per sample. The number of OTUs, Shannon diversity, and Shannon evenness was calculated in Mothur (v.1.39.5) (Kozich et al., 2013).
Statistical Analysis
All statistical analyses were performed in RStudio v. 1.1.383 (RStudio Team, 2016) with related packages based on R v. 3.5.0 (R Core Team, 2018). As the soil samples were collected from the 48 plots over four time points, the repeated-measures analysis of variance (rmANOVA) was performed to test the overall between-subjects and within-subject effects on soil bacterial community diversity and the relative abundances of individual taxa. Sampling time was treated as a repeated factor (a within-subject factor), block was a random factor, and fertilizer and crop treatments were fixed factors. As fertilizer was the main-plot factor, the overall main effect of fertilizer (df = 2) was tested against fertilizer × block (df = 6), and the effect of time × fertilizer was tested against time × fertilizer × block (df = 18). As crop was the sub-plot factor, the overall main effect of crop (df = 3) and the interaction of crop × fertilizer (df = 6) were tested against the main residual (df = 27). The effects of crop × time and crop × fertilizer × time were tested against the overall residual (df = 80).
To further test the fertilizer and the crop effects on soil bacterial community diversity and the relative abundance of individual taxon at each time point, the linear mixed-effects model was performed with lme function in nlme v.3.1-137 package (Pinheiro et al., 2018). In the model, according to the split-plot design, block was a random factor within block × fertilizer (frame), while fertilizer and crop factors were fixed factors with lme function (response variable ∼ fertilizer × crop, data = data, random = ∼ 1| block/frame). The normality of the distribution and the homogeneity of model residuals were tested by plotting the residuals against theoretical quantiles and fitted values. When the probability (P-value) of the effects was less than 0.05, the comparisons between two levels of the treatment and the Tukey tests were performed in package lsmeans (Lenth, 2016).
The repeated-measures correlations (rmcorr) were calculated to test the overall association of soil properties and crop growth characteristics with bacterial community diversity in package rmcorr v. 0.3.0 (Bakdash and Marusich, 2018). At each time point, the correlations of soil properties with bacterial community composition were tested by using Mantel test with the Spearman method.
Based on the Bray–Curtis dissimilarity between samples, a repeated-measures permutational multivariate analysis of variance (rm perMANOVA) was conducted using a permutation test with pseudo-F ratios (permutations = 999). This was performed using the adonis function in vegan v.2.5.2 (Oksanen et al., 2018) to test the community composition change over time and the overall treatment effects on the community composition. To know the community composition differences between treatments at each time point, perMANOVA was conducted at each time point. When P < 0.05, the comparisons between two levels of the treatment were performed with Pillai–Bartlett test (permutations = 999) in RVAideMemoire package (Hervé, 2018).
DESeq2 package v.1.20.0 (Love et al., 2014) was used for detecting the differential OTUs (log2 fold change in relative abundance) between two levels of the treatment based on the not subsampled OTU abundances. The differential OTUs were selected by the adjusted P < 0.01. All figures were prepared using the packages ggplot2 v.3.1.0 (Wickham, 2016) and randomcoloR v.1.1.0 (Ammar, 2017).
Sequence Accession Numbers
The sequence files have been deposited in the European Nucleotide Archive with project accession number PRJEB32412.
Results
The 27,535 16S rRNA gene OTUs were classified into 30 phyla (Figure 1), 82 classes (Supplementary Figure S1A), 198 orders (Supplementary Figure S1B), 355 families (Supplementary Figure S1C), and 698 genera (Supplementary Figure S1D). The most abundant phylum was Proteobacteria (relative abundance 0.25) that consisted of classes Alphaproteobacteria (0.163), Gammaproteobacteria (0.060), and Deltaproteobacteria (0.028). The other relatively abundant phyla were Actinobacteria (0.20), Acidobacteria (0.10), Bacteroidetes (0.11), Gemmatimonadetes (0.097), Chloroflexi (0.066), Planctomycetes (0.060), Verrucomicrobia (0.029), Patescibacteria (0.016), and Firmicutes (0.009) (Figure 1).
Figure 1. Soil bacterial phyla profile in a sub-boreal pasture. The abundances were defined based on 16S rRNA amplicon sequencing of soils from a split-plot field experiment with fertilizer and plant treatments. Fallow, no plant cover; mixture, red clover–timothy mixture; control, no fertilization; organic, organic fertilizer; synthetic, synthetic fertilizer.
Temporal Change of the Community Diversity and Composition
The bacterial richness [F(3,80) = 117.24, P < 0.001], Shannon diversity [F(3,80) = 72.54, P < 0.001], and Shannon evenness [F(3,80) = 27.53, P < 0.001] changed from year to year, but not from summer to autumn within a year (Table 1). The community diversity measures were generally higher in 2015 than in 2014 (Supplementary Table S2 and Figure 2). The Shannon diversity difference between crop treatments varied over time [F(9,80) = 2.22, P < 0.05] (Table 1). For example, the Shannon diversity in fallow was 0.15 lower than in planted plots in July 2015 [F(3,27) = 16.75, P < 0.001], whereas no difference was found in 2014 (Figure 2 and Supplementary Table S2). The overall repeated-measures correlation of the Shannon diversity with soil pH was 0.24 (P < 0.01), with NO it was 0.21 (P < 0.01), and with NH it was 0.38 (P < 0.001, Figure 3). The overall correlation of Shannon richness with crop dry matter yield was −0.20 (P < 0.05, Supplementary Figure S2) and with surplus nitrogen was 0.60 (P < 0.001, Supplementary Figure S3).
Table 1. Repeated-measures analysis of variance on soil bacterial community diversity in a sub-boreal pasture.
Figure 2. Boxplot (n = 12) of bacterial richness (A) and Shannon diversity (B) in a sub-boreal pasture over 2 years. The diversity indices were calculated based on 16S rRNA amplicon sequencing of soils from a split-plot field experiment with fertilizer (unfertilized control, organic fertilizer, and synthetic fertilizer) and plant (fallow, red clover, timothy, and a red clover–timothy mixture) treatments. Fallow: no plant cover; mixture, red clover–timothy mixture. The statistical analyses of differences between treatments at each time point and their temporal change are shown in Supplementary Table S2.
Figure 3. Repeated-measures correlation of soil bacterial diversity with soil pH, moisture, and NO3–N and NH4–N contents in a sub-boreal pasture. The diversity indices were calculated based on 16S rRNA amplicon sequencing of soils from a split-plot field experiment with fertilizer (unfertilized control, organic fertilizer, and synthetic fertilizer) and plant (fallow, red clover, timothy, and a red clover–timothy mixture) treatments.
The community composition changed over time [perMANOVA, F(3,80) = 7.42, P < 0.01], both year to year and from summer to autumn (Figure 4). The composition difference between crop treatments changed over time [perMANOVA, F(9,80) = 1.27, P < 0.01]. The changes over time exceeded the differences between crop treatments (Supplementary Figure S4A) and fertilizer treatments (Supplementary Figure S4B). Soil moisture (P < 0.05 in all time points) and soil total N and total C (P < 0.01 at the measured time points) contents consistently correlated with the community composition (Table 2 and Supplementary Table S3).
Figure 4. Non-metric multidimensional scaling of soil bacterial community change in a sub-boreal pasture over 2 years. The Bray–Curtis dissimilarities between samples were calculated based on 16S rRNA amplicon sequencing of soils from a split-plot field experiment with fertilizer (unfertilized control, organic fertilizer, and synthetic fertilizer) and plant (fallow, red clover, timothy, and a red clover–timothy mixture) treatments.
Table 2. Mantel test (Spearman) of the correlation between soil physicochemical factors and bacterial community composition in a sub-boreal pasture.
The relative abundances of Actinobacteria and Chloroflexi were higher in summer than in autumn, whereas those of Acidobacteria, Gammaproteobacteria, Bacteroidetes, Gemmatimonadetes, and Planctomecetes were lower (Figure 5).
Figure 5. Boxplot (n = 48) of temporal change of the most abundant phyla (class) in a sub-boreal pasture. The abundances are defined based on 16S rRNA amplicon sequencing of soils from a split-plot field experiment with fertilizer and plant treatments. Fallow, no plant cover; mixture, red clover–timothy mixture; control, no fertilization; organic, organic fertilizer; synthetic, synthetic fertilizer. The statistical analyses of differences between treatments and changes over time are shown in Supplementary Table S4.
Differences Between Fertilizer Treatments
At the end of the experiment, the Shannon evenness in the organic fertilizer treatment was higher than in the control and the synthetic fertilizer treatments [F(2,6) = 7.7, P < 0.05] (Supplementary Table S2). The community composition in the organic fertilizer treatment was different from those in the control and the synthetic fertilizer treatments [perMANOVA, F(2,6) = 1.54, P < 0.01] (Figure 6A). With time, the number of differential OTUs in the synthetic vs. the organic fertilizer treatments increased gradually from 2 to 45 (Table 3), but no differential OTUs were found between the control and the synthetic fertilizer plots. With the exception of one Chitinophagaceae OTU, the relative abundances of Bacteroidetes OTUs were higher in the organic fertilizer treatment than in the synthetic fertilizer treatment (Supplementary Figure S5A). At the end of the experiment, the relative abundances of Acidobacteria and Gemmatimonadetes OTUs were higher in the synthetic fertilizer treatment than in the organic fertilizer treatment (Supplementary Figure S5A).
Figure 6. Non-metric multidimensional scaling of soil bacterial communities in a sub-boreal pasture under different fertilizer (A) and plant (B) treatments. The abundances are defined based on 16S rRNA amplicon sequencing of soils from a split-plot field experiment. Fallow, no plant cover; mixture, red clover–timothy mixture; control, no fertilization; organic, organic fertilizer; synthetic, synthetic fertilizer.
Table 3. Number of differential operational taxonomic units (OTUs) between treatments in a sub-boreal pasture.
The relative abundances of Acidobacteria were generally higher in the control and the synthetic fertilizer treatments than in the organic fertilizer treatment (Supplementary Table S4). At the end of the experiment, the relative abundances of Proteobacteria and Fibrobacteres were higher in the organic fertilizer treatment than in the synthetic fertilizer treatment (P < 0.05), whereas those of Acidobacteria and Chloroflexi were higher in the control treatment than in the organic fertilizer treatment (P < 0.05) (Supplementary Table S4).
Differences Between Crop Treatments
In 2015, the community diversity measures in fallow plots were lower than those in planted plots (P < 0.05) (Supplementary Table S2). At each time point, the community composition in the fallow treatment was different from those in the planted treatments [perMANOVA, F(3,26) = 1.53 in June 2014, F(3,27) = 1.65 in September 2014, F(3,27) = 1.83 in July 2015, and F(3,27) = 1.64 in September 2015; P < 0.01] (Figure 6B and Supplementary Figure S6). The relative abundance of Fibrobacteres was higher and that of Nitrospirae was lower in the mixture and timothy treatments than in the fallow and red clover treatments (P < 0.05) (Figure 1 and Supplementary Table S4). The number of differential OTUs between the crop treatments changed over time (Table 3). At the end of the experiment, the number of differential OTUs was highest in timothy vs. fallow and lowest in red clover vs. timothy (Table 3). The differential OTUs in red clover vs. timothy treatments were mainly from the orders Flavobacteriales, Rhizobiales, and Saccharimonadales. Among the differential OTUs, the relative abundances of Gemmatimonadetes OTUs were higher in the fallow treatment than in the planted treatments (Supplementary Figure S5B). Compared to the unfertilized control, fertilization enhanced the temporal change in community composition in the mixture and the timothy plots (Supplementary Figure S7).
Discussion
We studied the dynamic response of soil bacteria to three fertilization and four crop management practices over four time points in a 3-year field study. As expected, the community diversity and the composition changed over time. The community diversity and the composition changed at different time scales, with the diversity changing between years, whereas the composition also changed between seasons within a year. Similarly, Hillebrand et al. (2018) reported that changes in species richness were often decoupled from the changes in composition.
The temporal change of the communities exceeded the differences between the fertilizer and the crop treatments, as also found in a very different climate in California, United States (Cruz-Martínez et al., 2009). Diversity increased slightly in all the treatments and correlated positively with soil NO and NH contents. The opposite and the negative correlations were found when the study was based on a single time point (Zeng et al., 2016), which emphasizes the importance of knowing the overall driving factors of diversity through time series analyses (Lauber et al., 2013). Microbial diversity is not a quantitative measurement with absolutes that require caution when comparing directly across different methods or studies (Mikkonen, 2012; Shade, 2017). Nevertheless, in accordance with our hypothesis, the higher N input in the third year resulted in greater differences between treatments: the number of OTUs differentiating organic from synthetic fertilizer treatments increased gradually from two in June 2014 to 45 in September 2015. The quality and the quantity of organic C influence the soil microbiome (Fierer, 2017), which may explain the OTU level differences between organic fertilization and other treatments. In contrast, no OTUs separated the control from the synthetic fertilizer treatment. The soil bacterial communities are primarily organic C limited, and only secondarily N limited (Soong et al., 2020), which may explain the similarity of the communities in the control and the synthetic fertilizer treatment. According to the report from Fierer et al. (2007), the continuous C input in the organic fertilizer treatment may have contributed to the higher abundances of the copiotrophic Bacteroidetes OTUs. The communities in the organic fertilizer plots were more even than those in the control and the synthetic fertilizer plots, as also found by Sun et al. (2004).
The community diversity and the composition differences were mainly between fallow and planted treatments. The differences between planted treatments were restricted to individual OTUs. The relative abundances of differential Flavobacteriales and Rhizobiales OTUs were higher in the red clover treatment than in the timothy treatment. Even though differences in absolute abundances cannot be concluded from differences in relative abundances (Gloor et al., 2017), it is likely that Rhizobiales abundances could have increased due to their capability to live in symbiosis with legumes (Lindström et al., 2019).
We found that the soil total C, total N, and moisture contents consistently correlated with the community differences between treatments. The differences between fallow and planted plots could be attributed to the higher soil C content under canopies and differences in soil moisture and temperature linked to crop shading (Trivedi et al., 2019). Soil moisture can directly affect microbial activity, decomposition rates, and nutrient diffusion rates (Manzoni et al., 2012). The higher relative abundance of Gemmatimonadetes in fallow may have been due to lower moisture since this phylum prefers drier soil (Fawaz, 2013). Other composition differences included the higher relative abundances of Planctomycetes and the lower relative abundances of Chloroflexi and Verrucomicrobia in the fallow than in the planted treatments. The differences in composition may indicate differences in the rates of biogeochemical processes (Schimel, 2001). However, a deeper understanding of the parallel changes between bacterial community compositions and the process rates requires further studies.
Generally, the soil had lower NO3–N content and higher moisture content in the autumn than in the summer (Li et al., 2019). From summer to autumn, the relative abundances of Acidobacteria and Bacteroidetes increased, while that of Actinobacteria decreased. In contrast, Lauber et al. (2013) found an opposite pattern across land-use types, indicating that the interaction between the local climate and the agricultural managements on soil microbiome is complex. In the present experiment, fertilization enhanced the temporal changes in community composition in the mixture and the timothy plots but not in the red clover plots, which further emphasizes the complexity of interactions among fertilizer and cropping managements in affecting soil bacterial communities.
Conclusion
In conclusion, we investigated the responses of soil bacterial diversity and community composition to pasture fertilizer and crop management procedures. Increased nitrogen input in the third year promoted diversity, which was associated with the overall increased pH as well as the increased availability of mineral N. Major differences in diversity and community composition were found between fallow and planted treatments, but not among the planted treatments, in spite of the higher botanical biodiversity in the mixture plots over the pure-stand plots. The relative abundances of individual OTUs differed between planted treatments and between organic and synthetic fertilizer treatments, such as the higher relative abundance of Rhizobiales OTUs in the clover-containing plots than in the timothy plots and the higher relative abundance of copiotrophic Bacteroidetes OTUs in the organic fertilizer plots than in the synthetic fertilizer plots. This suggests that the relative abundances of individual OTUs, in contrast to overall community diversity and composition, may help us to understand the mechanisms of bacterial response to agricultural practices and the possible connections of the response to the sustainability of agricultural practices.
Data Availability Statement
The datasets generated for this study can be found in the European Nucleotide Archive (ENA) with accession number PRJEB32412.
Author Contributions
HL did the soil DNA isolation, PCR, and sequences analysis under the supervision of AM and PP. HL did the statistical analysis and drafted the manuscript. All the authors agreed with the final version of the manuscript, contributed to interpreting the results, and revised the manuscript critically.
Funding
This work was supported by the Ministry of Agriculture and Forestry of Finland (KESTE project) and the Magnus Ehrnrooth Foundation.
Conflict of Interest
AM was employed by Kemira.
The remaining 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 Hannu Mikkola, Tapani Jokiniemi, and Markku Tykkyläinen for providing the field history information and Dr. P. J. Aphalo for the guidance on linear mixed-effects modeling. We acknowledge CSC-IT Center for Science, Finland, for computational resources. We appreciate the open-access publication fees support from Helsinki University Library. HL acknowledges the China Scholarship Council for a 4-year scholarship covering the stipend of her Ph.D. study at the University of Helsinki and Finnish Culture Foundation for a 6-month grant.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.01780/full#supplementary-material
References
Ammar, R. (2017). randomcoloR: Generate attractive random colors. R Package Version 1.1.0. Avaliable at: https://CRAN.R-project.org/package=randomcoloR (accessed November 24, 2019).
Bakdash, J. Z., and Marusich, L. R. (2018). rmcorr: Repeated Measures Correlation. R Package Version 0.3.0. Avaliable at: https://CRAN.R-project.org/package=rmcorr (accessed June 24, 2020).
Berg, G., and Smalla, K. (2009). Plant species and soil type cooperatively shape the structure and function of microbial communities in the rhizosphere. FEMS Microbiol. Ecol. 68, 1–13. doi: 10.1111/j.1574-6941.2009.00654.x
Cruz-Martínez, K., Suttle, K. B., Brodie, E. L., Power, M. E., Andersen, G. L., and Banfield, J. F. (2009). Despite strong seasonal responses, soil microbial consortia are more resilient to long-term changes in rainfall than overlying grassland. ISME J. 3, 738–744. doi: 10.1038/ismej.2009.16
de Graaff, M., Hornslein, N., Throop, H. L., Kardol, P., and van Diepen, L. T. (2019). “Effects of agricultural intensification on soil biodiversity and implications for ecosystem functioning: a meta-analysis,” in Advances in Agronomy, ed. D. Sparks (Amsterdam: Elsevier), 1–44. doi: 10.1016/bs.agron.2019.01.001
Dumanski, J., Eswaran, H., Pushparajah, E., and Smyth, A. J. (1991). Towards the Development of an International Framework. Bangkok: International Board for Soil Research and Management (IBSRAM). IBSRAM proceedings No 12.
Enwall, K., Nyberg, K., Bertilsson, S., Cederlund, H., Stenström, J., and Hallin, S. (2007). Long-term impact of fertilization on activity and composition of bacterial communities and metabolic guilds in agricultural soil. Soil Biol. Biochem. 39, 106–115. doi: 10.1016/j.soilbio.2006.06.015
Fawaz, M. N. (2013). Revealing the Ecological role of Gemmatimonadetes Through Cultivation and Molecular Analysis of Agricultural Soils. Master’s thesis, University of Tennessee, Knoxville, TN.
Fierer, N. (2017). Embracing the unknown: disentangling the complexities of the soil microbiome. Nat. Rev. Microbiol. 15, 579–590. doi: 10.1038/nrmicro.2017.87
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
Fierer, N., Lauber, C. L., Ramirez, K. S., Zaneveld, J., Bradford, M. A., and Knight, R. (2012). Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients. ISME J. 6, 1007–1017. doi: 10.1038/ismej.2011.159
Franche, C., Lindström, K., and Elmerich, C. (2009). Nitrogen-fixing bacteria associated with leguminous and non-leguminous plants. Plant Soil. 32, 35–59. doi: 10.1007/s11104-008-9833-8
Gloor, G. B., Macklaim, J. M., Pawlowsky-Glahn, V., and Egozcue, J. J. (2017). Microbiome datasets are compositional: and this is not optional. Front. Microbiol. 8:2224. doi: 10.3389/fmicb.2017.02224
Haichar, F. E. Z., Marol, C., Berge, O., Rangel-Castro, J., Prosser, J. I., Balesdent, J., et al. (2008). Plant host habitat and root exudates shape soil bacterial community structure. ISME J. 2, 1221–1230. doi: 10.1038/ismej.2008.80
Hillebrand, H., Blasius, B., Borer, E. T., Chase, J. M., Downing, J. A., Eriksson, B. K., et al. (2018). Biodiversity change is uncoupled from species richness trends: consequences for conservation and monitoring. J. Appl. Ecol. 55, 169–184. doi: 10.1111/1365-2664.12959
Hui, N., Jumpponen, A., Francini, G., Kotze, D. J., Liu, X., Romantschuk, M., et al. (2017). Soil microbial communities are shaped by vegetation type and park age in cities under cold climate. Environ. Microbiol. 19, 1281–1295. doi: 10.1111/1462-2920.13660
Jia, Z., and Conrad, R. (2009). Bacteria rather than archaea dominate microbial ammonia oxidation in an agricultural soil. Environ. Microbiol. 11, 1658–1671. doi: 10.1111/j.1462-2920.2009.01891.x
Ju, X., Xing, G., Chen, X., Zhang, S., Zhang, L., Liu, X., et al. (2009). Reducing environmental risk by improving N management in intensive Chinese agricultural systems. Proce. Natl. Acad. Sci. U.S.A. 106, 3041–3046. doi: 10.1073/pnas.0813417106
Kozich, J. J., Westcott, S. L., Baxter, N. T., Highlander, S. K., and Schloss, P. D. (2013). Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl. Environ. Microbiol. 79, 5112–5120. doi: 10.1128/AEM.01043-13
Lauber, C. L., Ramirez, K. S., Aanderud, Z., Lennon, J., and Fierer, N. (2013). Temporal variability in soil microbial communities across land-use types. ISME J. 7, 1641–1650. doi: 10.1038/ismej.2013.50
Lenth, R. V. (2016). Least-squares means: the R Package lsmeans. J. Stat. Softw. 69, 1–33. doi: 10.18637/jss.v069.i01
Li, H., Penttinen, P., Mikkola, H., and Lindström, K. (2019). Advantages of grass-legume mixture for improvement of crop growth and reducing potential nitrogen loss in a boreal climate. Agric. Food Sci. 28, 176–189. doi: 10.23986/afsci.75600
Lindström, K., and Mousavi, S. A. (2019). Effectiveness of nitrogen fixation in rhizobia. Microb. Biotechnol. 1–22. doi: 10.1111/1751-7915.13517
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8
Manzoni, S., Schimel, J. P., and Porporato, A. (2012). Responses of soil microbial communities to water stress: results from a meta-analysis. Ecology 93, 930–938. doi: 10.1890/11-0026.1
Marschner, P., Kandeler, E., and Marschner, B. (2003). Structure and function of the soil microbial community in a long-term fertilizer experiment. Soil Biol. Biochem. 35, 453–461. doi: 10.1016/S0038-0717(02)00297-3
Martin-Guay, M., Paquette, A., Dupras, J., and Rivest, D. (2018). The new green revolution: sustainable intensification of agriculture by intercropping. Sci. Total Environ. 615, 767–772. doi: 10.1016/j.scitotenv.2017.10.024
Mikkonen, A. (2012). The Potential of Microbial Ecological Indicators to Guide Ecosophisticated Management of Hydrocarbon-Contaminated Soils. Helsinki: University of Helsinki.
Niemi, J., and Väre, M. (2017). “Finnish agriculture and food sector 2016/17,” in Natural Resources and Bioeconomy Studies 49/2017, eds J. Niemi and M. Väre (Washington, DC: Natural Resources Institute), 69–79.
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2018). vegan: Community Ecology Package. R Package Version 2.5-2. Avaliable at: https://CRAN.R-project.org/package=vegan (accessed September 01, 2019).
Orgiazzi, A., Bardgett, R. D., and Barrios, E. (2016). Global Soil Biodiversity Atlas. Brussels: European Commission.
Philippot, L., Hallin, S., and Schloter, M. (2007). Ecology of denitrifying prokaryotes in agricultural soil. Adv. Agron. 96, 249–305. doi: 10.1016/S0065-2113(07)96003-4
Pinheiro, J., Bates, D., DebRoy, S., Sarkar, D., and R Core Team (2018). nlme: Linear and Nonlinear Mixed Effects Models 3.1-137. Avaliable at: https://CRAN.R-project.org/package=nlme (accessed May 24, 2020).
R Core Team (2018). R: A Language and Environment for Statistical Computing 2018. Vienna: R Foundation for Statistical Computing.
Ramirez, K. S., Lauber, C. L., Knight, R., Bradford, M. A., and Fierer, N. (2010). Consistent effects of nitrogen fertilization on soil bacterial communities in contrasting systems. Ecology 91, 3463–3470. doi: 10.1890/10-0426.1
Robertson, G. P., and Groffman, P. M. (2015). “Chapter 14 - nitrogen transformations,” in Soil Microbiology, Ecology and Biochemistry, 4th Edn, ed. E. Paul (Cambridge, MA: Academic Press), 421–446. doi: 10.1016/b978-0-12-415955-6.00014-1
Schimel, J. (2001). “1.13 - biogeochemical models: implicit versus explicit Microbiology,” in Global Biogeochemical Cycles in the Climate System, eds E. D. Schulze, D. S. Schimel, J. Lloyd, S. Harrison, I. C. Prentice, M. Heimann, et al. (Cambridge, MA: Academic Press), 177–183.
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Shade, A. (2017). Diversity is the question, not the answer. ISME J. 11, 1–6. doi: 10.1038/ismej.2016.118
Shen, J., Zhang, L., Guo, J., Ray, J. L., and He, J. (2010). Impact of long-term fertilization practices on the abundance and composition of soil bacterial communities in Northeast China. Appl. Soil Ecol. 46, 119–124. doi: 10.1016/j.apsoil.2010.06.015
Soong, J. L., Fuchslueger, L., Marañon-Jimenez, S., Torn, M. S., Janssens, I. A., Penuelas, J., et al. (2020). Microbial carbon limitation: the need for integrating microorganisms into our understanding of ecosystem carbon cycling. Global Change Biol. 26, 1953–1961. doi: 10.1111/gcb.14962
Stott, K. J., and Gourley, C. J. P. (2016). Intensification, nitrogen use and recovery in grazing-based dairy systems. Agric. Syst. 144, 101–112. doi: 10.1016/j.agsy.2016.01.003
Strickland, M. S., Lauber, C., Fierer, N., and Bradford, M. A. (2009). Testing the functional significance of microbial community composition. Ecology 90, 441–451. doi: 10.1890/08-0296.1
Sun, H. Y., Deng, S. P., and Raun, W. R. (2004). Bacterial community structure and diversity in a century-old manure-treated agroecosystem. Appl. Environ. Microbiol. 70, 58–68. doi: 10.1128/AEM.70.10.5868-5874.2004
Tilman, D., Cassman, K. G., Matson, P. A., Naylor, R., and Polasky, S. (2002). Agricultural sustainability and intensive production practices. Nature 418, 671–677. doi: 10.1038/nature01014
Trivedi, C., Reich, P. B., Maestre, F. T., Hu, H., Singh, B. K., and Delgado-Baquerizo, M. (2019). Plant-driven niche differentiation of ammonia-oxidizing bacteria and archaea in global drylands. ISME J. 13, 2727–2736. doi: 10.1038/s41396-019-0465-1
Van, D. H., 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
Verwaal, J., Rutten, M., Van Logtestijn, R., Bakker, R., Scheublin, T. R., Van, D. H., et al. (2006). Symbiotic bacteria as a determinant of plant community structure and plant productivity in dune grassland. FEMS Microbiol. 56, 178–187. doi: 10.1111/j.1574-6941.2006.00086.x
Wang, J., Song, Y., Ma, T., Raza, W., Li, J., Howland, J. G., et al. (2017). Impacts of inorganic and organic fertilization treatments on bacterial and fungal communities in a paddy soil. Appl. Soil Ecol. 112, 42–50. doi: 10.1016/j.apsoil.2017.01.005
Yoshida, M., Ishii, S., Otsuka, S., and Senoo, K. (2009). Temporal shifts in diversity and quantity of nirS and nirK in a rice paddy field soil. Soil Biol. Biochem. 41, 2044–2051. doi: 10.1016/j.soilbio.2009.07.012
Keywords: legume–grass mixture, bacterial community, seasonal change, pasture, sustainable agriculture
Citation: Li H, Penttinen P, Mikkonen A, Stoddard FL and Lindström K (2020) Response of Soil Bacterial Community Diversity and Composition to Time, Fertilization, and Plant Species in a Sub-Boreal Climate. Front. Microbiol. 11:1780. doi: 10.3389/fmicb.2020.01780
Received: 23 September 2019; Accepted: 07 July 2020;
Published: 05 August 2020.
Edited by:
Helena Freitas, University of Coimbra, PortugalReviewed by:
Weibing Xun, Nanjing Agricultural University, ChinaDiana Martín-Lammerding, Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria (INIA), Spain
Copyright © 2020 Li, Penttinen, Mikkonen, Stoddard and Lindström. 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: Honghong Li, honghong.li@helsinki.fi; Petri Penttinen, petri.penttinen@helsinki.fi