ORIGINAL RESEARCH article
Sec. Terrestrial Microbiology
Long-term forest soil warming alters microbial communities in temperate forest soils
- 1Department of Microbiology, University of Massachusetts, Amherst, MA, USA
- 2Department of Natural Resources and the Environment, University of New Hampshire, Durham, NH, USA
- 3Department of Biology, University of Massachusetts, Amherst, MA, USA
- 4Marine Biological Labs, Woods Hole, MA, USA
Soil microbes are major drivers of soil carbon cycling, yet we lack an understanding of how climate warming will affect microbial communities. Three ongoing field studies at the Harvard Forest Long-term Ecological Research (LTER) site (Petersham, MA) have warmed soils 5°C above ambient temperatures for 5, 8, and 20 years. We used this chronosequence to test the hypothesis that soil microbial communities have changed in response to chronic warming. Bacterial community composition was studied using Illumina sequencing of the 16S ribosomal RNA gene, and bacterial and fungal abundance were assessed using quantitative PCR. Only the 20-year warmed site exhibited significant change in bacterial community structure in the organic soil horizon, with no significant changes in the mineral soil. The dominant taxa, abundant at 0.1% or greater, represented 0.3% of the richness but nearly 50% of the observations (sequences). Individual members of the Actinobacteria, Alphaproteobacteria and Acidobacteria showed strong warming responses, with one Actinomycete decreasing from 4.5 to 1% relative abundance with warming. Ribosomal RNA copy number can obfuscate community profiles, but is also correlated with maximum growth rate or trophic strategy among bacteria. Ribosomal RNA copy number correction did not affect community profiles, but rRNA copy number was significantly decreased in warming plots compared to controls. Increased bacterial evenness, shifting beta diversity, decreased fungal abundance and increased abundance of bacteria with low rRNA operon copy number, including Alphaproteobacteria and Acidobacteria, together suggest that more or alternative niche space is being created over the course of long-term warming.
Earth's climate is warming, and this is exacerbated by both biophysical (e.g., albedo) and biogeochemical [e.g., carbon (C) cycle] feedbacks (IPCC, 2013). Microbes are key players in every biogeochemical cycle, regulating greenhouse gas fluxes between soils and the atmosphere (Falkowski et al., 2008). Despite their pivotal role, we know little about how microbes respond to environmental change, and microbial dynamics are only beginning to be represented in ecosystem models (Reid, 2011; Treseder et al., 2012). Based on the improved predictive capacity of soil carbon cycling models when microbial physiology is considered, it is clear that microbes are important for understanding and predicting ecosystem processes (Allison et al., 2010; Li et al., 2014). New genomic approaches hint at what the most abundant organisms are, though their ecological roles are still unclear. A better understanding of microbial dynamics is critical for projecting the rate and magnitude of climate change (Melillo et al., 1990, 2002; Heimann and Reichstein, 2008).
At the Harvard Forest Long-Term Ecological Research (LTER) site, warming-induced disruptions to ecosystem C cycles have been observed as part of three ongoing warming studies. In the longest-running site (20 years, Table 1), soil respiration rates in the warming treatment were initially higher than controls, then slowed to become equal to controls after 12 years (Melillo et al., 2002). Most warming studies, including our own, have observed a short-term increase in CO2 emissions with warming (Rustad et al., 2001; Contosta et al., 2011; Wu et al., 2011; Davidson et al., 2012) as well as slowed CO2 emissions following this initial increase. Slowed CO2-C loss is concomitant with a depletion of labile C (Bradford et al., 2008; Frey et al., 2008; Melillo et al., 2011), as well as a decline in microbial biomass (Frey et al., 2008), thermal adaptation of soil microbes (Bradford et al., 2008) and a shift in microbial carbon use efficiency (Frey et al., 2013). More recently in our longest running (20-year) warming experiment, soil respiration has begun to increase again in warmed soils compared to controls. Partitioning the two components of respiration—microbial and root respiration—showed that on an annual basis, the majority of CO2 (70–80%) from forest soil plots was microbial regardless of heat treatment (Melillo et al., 2002, 2011). Climate warming and other cascading effects (such as drying) have been suggested to destabilize stored soil C in a manner that is dependent on organic matter age, protection, and other conditions that may alter microbial access and decay (Fierer et al., 2009; Schimel and Schaeffer, 2012). The non-linear response of soil respiration to warming, combined with the prior observations of microbial adaptation and that most soil respiration is microbial has lead us to the hypothesis that long-term warming has caused changes in the soil microbial community, which may be exacerbating C cycle feedbacks to the climate system.
Table 1. Experimental sites that make up the warming chronosequence at Harvard Forest Long-Term Ecological Research (LTER) site.
Studies that have examined both microbial community characteristics and activity as affected by warming, whether as lab experiments (Zogg et al., 1997; Frey et al., 2008) or field studies (Bardgett et al., 1999; Belay-Tedla et al., 2009; Wu et al., 2011; Zhou et al., 2012; Pold and DeAngelis, 2013) reveal that although the physiological response (e.g., CO2 flux) across experiments and biomes tends to be consistent in that warming increases C mineralization, the microbial community response is not, with microbial communities in different ecosystems responding differently to warming. Over the short term the net effect of warming on soil microbes tends to be increased microbial activity, including increased soil respiration (Rustad et al., 2001; Melillo et al., 2011; Wu et al., 2011). This response to warming is likely partly attributable to changes in the active fraction of the biomass, rather than changing community's constituents (Zogg et al., 1997; Andrews et al., 2000; Pietikåinen et al., 2005). A metaanalysis of 75 manipulative climate change experiments showed that not all soil microbial communities respond similarly to warming (Blankinship et al., 2011). Long-term experiments at the Kessler Farm Field Laboratory (KFFL) in the plains of central Oklahoma found increased diversity under warming and drought, suggesting that warming may have somehow “primed” the community to be more resilient and resistant to further disturbance (Zhang et al., 2005). Furthermore, Zhang and colleagues found that warming was associated with decreased net N mineralization and a significant shift in the substrate utilization profiles, indicating a change in the substrate availability for the community. Metagenomic sequencing on grassland soils from experimental plots that had been warmed continuously (+2°C) for 8 years showed that gene signatures for degradation of more labile C compounds—starch, chitin, hemicellulose, and cellulose—were activated, whereas genes involved in recalcitrant C degradation, such as lignin, were not stimulated by warming (Zhou et al., 2012; Luo et al., 2014). They concluded that warming in grasslands could result in a weakened positive feedback between the terrestrial carbon cycle and climate. While results from the early years of the Harvard Forest warming experiment are consistent with the grasslands study, the long-term effects point to a very different conclusion about climate feedbacks.
The Harvard Forest soil warming experiments offer a unique opportunity to understand how climate change affects soil microbial community composition over the course of long-term warming. This is important because bacterial activities are possibly involved in a positive feedback to climate (Melillo et al., 2002; Frey et al., 2008; Brzostek et al., 2012). Three forest plots experienced 5°C above ambient soil temperatures for 5 (SWaN Plots), 8 (Barre Woods) and 20 (Prospect Hill) years at the time of sampling, forming an experimental warming “chronosequence” (Table 1). Fatty acid methyl ester (FAME) analysis of the Prospect Hill warming experiment after 12 years of heating 5°C above ambient showed a significant decrease in fungal abundance and an increase in Gram-positive bacteria and Actinobacteria in warmed soils compared to controls (Frey et al., 2008). Previous results indicate observed changes in C source utilization (Frey et al., 2008), plant-microbial interactions (Butler et al., 2012), mass-specific microbial respiration rates (Bradford et al., 2008), and soil chemistry (Melillo et al., 2011, 2002). Together, these results suggest a substantial change in microbial substrate availability which would be consistent with changing bacterial community profiles and potential altered feedbacks to climate.
In this study, we used high-throughput sequencing of the V4 variable region of the 16S ribosomal RNA gene to test the hypothesis that bacterial communities have changed over two decades of simulated climate change in organic and mineral soil horizons. However, utilization of the ribosomal RNA gene for phylogenetic analysis carries the caveat that the rRNA operon often exists at multiple copies per genome, which both confounds estimates of genomic abundance by over-estimating taxa with multiple rRNA operon copies (Crosby and Criddle, 2003; Větrovský and Baldrian, 2013) but also suggests life strategy, where taxa with multiple rRNA operon copies have faster growth rates, while taxa with single copies tend to be associated with oligotrophic environments (Klappenbach et al., 2000; Stevenson and Schmidt, 2004). We estimated copy number by phylogenetic inference using the available complete genome sequences in GenBank (Kembel et al., 2012), and used this to examine the effect of copy number correction on community profiles. We also used these estimates to test the additional hypothesis that average bacterial copy number was decreased by long-term warming, which would be consistent with observations of microbial adaptation (Bradford et al., 2008) and decreased soil carbon with long-term warming (Melillo et al., 2011). Finally, we performed quantitative PCR (Q-PCR) of bacterial and fungal communities, as well as the phyla Actinobacteria, Acidobacteria, and class Alphaproteobacteria, to measure changes in absolute abundance of these dominant and dynamic groups.
Materials and Methods
All three experimental sites at the Harvard Forest Long-Term Ecological Research (LTER) site (Petersham, MA) were established in mixed hardwood forest stands with dominant tree species being paper and black birch (Betula papyrifera and lenta), red maple (Acer rubrum), black and red oak (Quercus velutina and rubra), and American beech (Fagus grandifolia). The soils are coarse-loamy inceptisols. The warmed plots at all three sites have been heated continuously by the use of resistance heating cables buried at 10 cm depth in the mineral horizon (Melillo et al., 2002). Controls in the longest running experiment are disturbance controls (which had cables buried but never activated); effects of disturbance on microbial community composition (Frey et al., 2008), soil inorganic nitrogen, and carbon dioxide flux (Peterjohn et al., 1994) have been minimal. Heating was to 5°C above ambient, a temperature increase which falls within the range of worst-case-scenario model projections for increased global air temperature by the year 2100 as projected by the Intergovernmental Panel on Climate Change (IPCC, 2013). Mean monthly temperatures at Harvard Forest range from 19°C in July to −5°C in January and average annual precipitation is 112 cm, distributed relatively evenly throughout the year.
Soil samples were collected October 25–27, 2011 from each of the three sites in the chronosequence with four replicates from each treatment at Barre Woods, Prospect Hill and SWaN (Table 1). Soil was collected from the organic horizon (which varies in depth from 1–5 cm) and the upper 10 cm of mineral soil in the heated and control plots of all three studies, with four replicates totaling 48 samples (3 sites × 2 warming treatments × 2 soil horizons × 4 replicates). The organic horizon was collected as intact 20 × 20 cm blocks to the depth of the mineral soil, and the underlying mineral soil was sampled to 10 cm depth in the same locations using a custom-made stainless steel auger (9 cm diameter). Subsamples of both horizons were taken and immediately flash frozen in the field and stored at -80°C until extraction.
Nucleic Acid Extraction
Soils were extracted twice based on previously published methods (DeAngelis et al., 2010) with a few modifications. Frozen soils were extracted for DNA and RNA simultaneously (Griffiths et al., 2000) using modified CTAB extraction buffer (0.25 M phosphate buffer (pH 8), 5% hexadecyltrimethylammonium bromide (CTAB) in 1M NaCl) and 50 μl of 0.1 M ammonium aluminum sulfate (Braid et al., 2003). Three replicate extractions were performed for each sample, then pooled and put through the Qiagen All DNA/ RNA Mini kit (Qiagen, Valencia, CA).
Amplification and Sequencing
Library generation and sequencing of the V4 region of the 16S ribosomal RNA gene were performed as per a recently published protocol (Caporaso et al., 2012). The V4 region of the 16S ribosomal RNA gene was amplified on an Eppendorf ProS thermal cycler using the primers 515F and 806R, where the forward primers contained a subset of 48 of the 8 bp barcoded primers. Reactions were performed in a final volume of 25 μl using Takara ExTaq with 200 pM of each primer, 25 μg of BSA and 2 units of DNA polymerase (Takara Mirus Bio, Madison, WI) with 10 ng template per reaction. PCR amplifications were performed at 50°C annealing temperatures (Tm), with an initial denaturation (5 min) followed by 30 cycles of 95°C (30 s), Tm (25 s) and 72°C (120 s), and a final extension of 72°C (10 min). Triplicate amplification reactions were verified by agarose gel electrophoresis, pooled, then cleaned using Qiagen MinElute kit (Qiagen Sciences, Valencia, CA). Cleaned amplicon pools were quantified using picogreen and quality assessed by nanodrop. Sequencing was performed by the Molecular Biology Core Facility at the Dana-Farber Cancer Institute using standard operating procedures. Due to the low diversity of the library 50% by DNA mass PhiX spike was added to the pooled, barcoded sample just before running. Samples were sequenced using the MiSeq platform to generate 2 × 150 bp paired-end reads.
Sequencing Data Analysis
Paired-end sequences were assessed for quality using FastQC (Andrews, 2010), and initial quality filtering and assembly was performed using FLASH using default parameters (Magoč and Salzberg, 2011). Sequences were then binned into operational taxonomic units (OTUs) and taxonomies assigned using the subsampled open-referenced OTU picking method in QIIME (Caporaso et al., 2010) at 99% sequence identity based on the October 2012 greengenes taxonomy (DeSantis et al., 2006), where the fasta reference file was truncated to include just the V4 region (Werner et al., 2012). Chimeric sequences were detected using UCHIME (Edgar et al., 2011; Wright et al., 2012). We removed sequences observed only once or twice (singletons or doubletons), as well as erroneous sequences that were probable chimeras (DeSantis et al., 2006). A second community matrix of dominant taxa was generated with only taxa present in the data set at 0.1% relative abundance or greater. Community matrices were rarefied to the number of taxa in the sample with the lowest number of observations for most analyses excluding diversity analyses, since there tends to be a positive correlation between richness and depth of sequence sampling. Sequences are available in GenBank under accession numbers SRP040706, BioProject ID PRJNA242868.
Total bacterial and fungal abundances were measured for all 48 samples by quantitative PCR (Q-PCR) using 341F (5′-CCT ACG GGA GGC AGC AG-3′) and 534R (5′-ATT ACC GCG GCT GCT GGC-3′) (Muyzer et al., 1993) primers for total bacteria and ITS1f (5′-TCC GTA GGT GAA CCT GCG G-3′) and 5.8 s (5′-CGC TGC GTT CTT CAT CG-3′) (Fierer et al., 2005) for total fungi. Standard curves were based on linear PCR product from Klebsiella sp. str. BRL6-2 for bacteria and Saccharomyces cerevisiea for fungi, after purification, quantification and dilution. All samples were run as technical duplicates on an Eppendorf Realplex2, using the QuantiFast SYBR Green PCR Kit (Qiagen) with 10ng of DNA per 25ul reaction and primers at a final concentration of 1uM. PCR program was 5 min 95°C initially, followed by 40 × (10 s 95°C, 10 s 58°C, 20 s 60°C), and melt curve (60°C–95°C). Efficiencies were 101.6% for fungal, with an R2 of 0.9923, and 107.8% for bacteria, with an R2 of 0.9972. Total Acidobacteria, Actinobacteria and Alphaproteobacteria abundances were measured for all eight Prospect Hill organic horizon samples by quantitative PCR (Q-PCR). Primer sets for Acidobacteria were Acid31 (5′-GAT CCT GGC TCA GAA TC-3′) and Eub518 (5′-ATT ACC GCG GCT GCT GG-3′) (Fierer et al., 2005); for Actinobacteria were Act920F3 (5′-TAC GGC CGC AAG GCT A-3′) and Act1200R (5′-TCR TCC CCA CCT TCC TCC G-3′); and for Alphaproteobacteria were α682F (5′-CIA GTG TAG AGG TGA AAT T-3′) and 908αR (5′-CCC CGT CAA TTC CTT TGA GTT-3′) (Bacchetti De Gregoris et al., 2011). To generate standard curves, we used Micrococcus luteus (ATCC 381), Caulobacter crescentus str. CB15N (Peter Chien, pers. comm.), and Acidobacterium capsulatum (DSMZ 11244) for Actinobacteria, Alphaproteobacteria, and Acidobacteria, respectively. These standard curves were based on PCR product except for Acidobacteria, which was based on plasmid. The identity of all organisms and specificity of amplified products were confirmed by Sanger sequencing. Actinobacteria had a two-step amplification cycle with 10 s at 95°C followed by 30 s at 61.5°C. Acidobacteria and Alphaproteobacteria had a three-step cycle with 10 s at 95°C, 20 s at 55°C, and 10 s at 60°C. Soil dry weights were obtained by drying 10 g fresh soil in a 105°C oven for 3 days, and abundances are reported as counts per dry weight of soil. Counts represent genome equivalents not considering the confounding factor of multiple small subunit ribosomal RNA operon copies per genome.
The experimental design was a fully replicated design with three warming experiments: Prospect Hill (20 years warming), Barre Woods (8 years warming), and SWaN (5 years warming); one treatment of warming or no-heated control; two soil horizons (organic and mineral); and four field replicates, resulting in 48 samples total. For the copy number correction, we used a recently published method (Kembel et al., 2012), where the copy number corrected OTU table was generated using copy numbers estimated in R based on ancestral state reconstruction (Matsen et al., 2010). Weighted average copy number is based on multiplying relative abundance of OTUs by their estimated copy number, then calculating average copy number per sample. The UniFrac distance matrices were generated using FastUnifrac (Hamady et al., 2009). Beta diversity was estimated based on weighted Unifrac distance matrices used to generate principal coordinates plots, and Procrustes analysis was performed to compare principal component plots of copy number corrected environmental file and uncorrected environmental file. Evenness was measured by Pielou's J, and richness measured as total number of taxa. We used seven methods to identify indicator OTUs of warming in the Prospect Hill organic horizons: (1) indicator value analysis (Dufrêne and Legendre, 1997); (2) volcano plots using two as the minimal fold change, 0.05 as the p value threshold in a t-test adjusted for multiple comparisons completed in METAGENassist after normalizing the data using Pareto Scaling (Arndt et al., 2012); (3) Nearest Shrunken Centroid classification (NSC) (Tibshirani et al., 2003); (4) partial least squares discriminatory analysis; (5) Bayesian groupings; (6) rank mobility based on mean abundance and reporting the top 5% of OTUs which showed the greatest change in rank; (7) and paired Student's T-test using the Benjamini-Hochberg correction (Benjamini and Hochberg, 1995). To compare whether Q-PCR copies differed between warmed and control plots, we used Bayesian inference (Kruschke, 2013), because this method yields richer inference considering the relatively small sample size. All statistical analyses were performed in R using the RStudio interface (RStudio, 2012; R Core Team, 2014), including packages reshape (Wickham, 2012), vegan (Oksanen et al., 2011), ggplot2 (Wickham, 2009), limma (Smyth, 2004), pplacer (Matsen et al., 2010), indicspecies (De Cáceres and Legendre, 2009), pamr (Tibshirani et al., 2003), caret (Kuhn, 2008), phyloseq (McMurdie and Holmes, 2013), pls (Mevik and Wehrens, 2007), and BEST (Kruschke, 2013).
Sequencing produced 3,487,689 high-quality sequenced observations of the V4 16S ribosomal RNA gene region (Table S1, Figure S1). Sequencing depths before rarefaction ranged from 33,545 to 112,502 sequences per sample, with a mean of 63,293 reads and median of 62,838 reads. For most analyses, communities were rarefied to 33,545 sequences. Sequence reads were binned into operational taxonomic units (OTUs) based on 99% sequence identity (Werner et al., 2012), both because bacteria with nearly-identical 16S rDNA sequences may represent variable genotypes and different species (Suau et al., 1999; Hold et al., 2002; Konstantinidis and Tiedje, 2005) and because many functional traits are phylogenetically conserved up to 0.01% ribosomal RNA gene sequence dissimilarity (Martiny et al., 2013). After removal of singletons, doubletons and chimeras 2,938,751 observations remained (Table S1) that were binned into 45,875 OTUs, a scale of diversity on par with current estimates of soil bacterial diversity (Torsvik et al., 2002; Gans et al., 2005; Roesch et al., 2007).
Microbial Community Responds to Soil Warming After 20 Years
Soil warming had a statistically significant impact on bacterial community structure, but only after 20 years of warming. After 20 years, organic warmed soils were significantly different from those of control plots (Figure 1A). The sites warmed for 5 or 8 years showed no significant treatment effect on beta diversity, but bacterial communities in organic horizons warmed for 20 years began to resemble mineral soils. The strongest effect on bacterial community structure in our analysis was soil horizon (organic versus mineral), followed by site and then warming treatment (Table S2), with no significant interaction between soil horizon and treatment.
Figure 1. Measures of diversity for the warming chronosequence, based on the dominant subset community (N = 155). (A) One PCoA ordination was performed on all sites, both treatments and both soil types, but the three sites are shown separately for clarity: SWaN Plots, Barre Woods, and Prospect Hill. (B) Diversity as measured by Shannon's H index is shown for the three sites as box (first and third quartiles) and whisker (95% CI) plots where the solid bar is median. (C) Rank abundance curve of the 155 dominant species in the community, with inset showing the three most abundant taxa and their relative abundance in heated versus control treatments averaged for the three sites. Statistical significance is indicated as *P < 0.05, **P < 0.01.
While soil warming affected community structure only in the organic horizon, diversity as measured by Shannon's H index was higher at the 20-year warmed Prospect Hill site (P < 0.01, Figure 1B, Table S3A), and was on average lower at the SWaN site compared to the other two sites (P < 0.001, Table S3B). In the Prospect Hill site, changes in diversity were driven more by increasing evenness (T-test P < 0.01; Bayesian effect size 1.60, 95% range −3.1 to −0.16) than changes in richness (T-test P = 0.06; Bayesian effect size −0.983, 95% range −2.05–0.20). The observed increase in bacterial diversity was driven strongly by decreased abundance of a few dominant taxa.
Dominant Subset Reflects Total Community Response to Warming
The bacterial community was observed to be highly uneven: only 155 taxa (0.3% of total) present at 0.1% relative abundance or higher accounted for over half of the 3 million observations (Figure 1C). Comparison of the full community to the dominant subset community showed that the two were highly correlated (Procrustes R = 0.987, P < 0.001; Mantel test R = 0.980, P < 0.001). The same trends emerged for the dominant community as were observed for the full community, where soil horizon, then site, were major drivers of community structure, and that warming only affected the organic horizon of the 20-year long-term warmed soils (Table S2). This is perhaps not surprising considering that these dominant taxa, though comprising less than 1% of the richness, represented over half of all observations. The treatment effect on organic soil communities was stronger when looking at the dominant compared to the full community (Table S2), suggesting that dominant taxa were more strongly affected by warming.
The composition of the dominant subset community was similar to that of the total community, though the richness and phylogenetic range were much depleted. Compared to the 43,909 OTUs found in 33 phyla and 133 families in the total rarefied community, the dominant subset community was comprised of 155 OTUs found in 6 phyla and only 19 families. These included the phyla (and subphyla) Acidobacteria (57 OTUs), Alphaproteobacteria (40), Actinobacteria (25), Verrucomicrobia (11), Gammaproteobacteria (11), Planctomycetes (6), Firmicutes (2), Betaproteobacteria (2), and Deltaproteobacteria (1).
Most bacterial taxa in the dominant subset responded positively to long-term warming (Figure 2), though we observed no change in absolute abundance of total bacteria with warming treatment by Q-PCR (Figure 3A, Table S4). In the organic soil horizon, bacterial abundance was unaffected by warming, while fungal abundance was marginally decreased with warming (Bayesian effect size 0.744; t-test P = 0.07). Power analysis showed that 6.5 biological (field) replicates per group would distinguish the treatment effect on fungi in the organic horizon (for power = 0.90, significance level 0.95). Fungal ITS Q-PCR counts in the mineral horizon were unaffected by warming. The Q-PCR data showed that there were more fungal rRNA operon copies in the organic horizon compared to the mineral soil (P < 0.05), with the same amount of bacterial rRNA operon copies in the two soils. Bacteria 16S ribosomal RNA gene counts outnumbered fungal ITS counts by an order of magnitude in the organic layer and mineral layers (P < 0.001 for both soils).
Figure 2. Phylogenetic tree of the dominant subset community (n = 155) in the Prospect Hill 20-year warmed sites, where (from left to right) the color strips denote phylum-level classification; the circles denote a significant indicator species test (P < 0.05) with blue circles indicative of control treatments and red circles indicative of heated treatments; bars show fold change in OTU abundance with warming treatment. The seven indicator species analyses are (1) Dufrêne and Legendre's IndVal, (2) fold change (volcano-plots), (3) nearest shrunken centroid, (4) PLS-DA loadings test, (5) Bayesian group comparison, (6) Rank abundance test, (7) Student's t-test. Taxa with two or more significant indicator species tests are marked (*).
Figure 3. Quantitative PCR of bacteria and fungi, and the phyla Actinobacteria, Acidobacteria and Alphaproteobacteria showing abundance of dominant microbial phylogenetic groups for both treatments in organic horizon (averaged by site, which was not a significant factor). Fungi were less abundant, and bacteria unchanged in heated (+5°C) compared to control (Ctl) organic horizon soils across the three sites (A). In the organic horizons of the Prospect Hill site, Actinobacteria were unchanged, while Alpha-proteobacteria and Acidobacteria were enriched by long-term warming (B). Means are shown as as box (first and third quartiles) and whisker (95% CI) plots where the solid bar is median.
Select Taxa were Especially Sensitive to 20 Years of Warming
To identify the taxa responsible for the observed changes in the organic horizon bacterial community after 20 years of warming, indicator species analysis was performed using seven different methods (Table S5). Taxa that are more likely indicators of warmed or control conditions for the 20-year warmed site were defined as significant by multiple indicator species analyses (Figure 2). These analyses revealed that 19 of the 155 dominant taxa were differentially abundant in the organic soil horizon of the heated compared to control plots, with 15 increased in the warmed plots and four decreased. All six Acidobacteria detected as indicator species were indicators of heating treatment, with increased abundance in warmed compared to control soils. The four Actinobacteria were divided evenly between being indicators for heated or control treatments, and all belonged to the class Actinobacteria and order Actinomycetales. The last nine species were all Proteobacteria, and the two Proteobacteria that were indicators for control treatments were both Alphaproteobacteria in the order Rhodospirilales. The seven Proteobacteria that were indicators for warming treatments consisted of five Rhizobiales (Alphaproteobacteria), Syntrophobacterales (Deltaproteobacteria), and Xanthomonadales (Gammaproteobacteria). Measures of absolute abundance (based on Q-PCR) of Actinobacteria, Alphaproteobacteria and Acidobacteria were performed because taxa in these groups were dominant in our MiSeq observations: 15% Actinobacteria (25 of 155 OTUs), 26% Alphaproteobacteria (40 of 155), and 37% Acidobacteria (57 of 155). Absolute copies of Actinobacteria did not differ between warmed and control plots, though heated plots had more Alphaproteobacteria (P < 0.05) and trended toward having increased abundances of Acidobacteria (P = 0.08) compared to control plots in the Prospect Hill organic horizon (Figure 3B, Table S4).
Mean Ribosomal RNA Copy Number is Significantly Depleted by Warming
The organic horizon in the 20-year warmed site (Prospect Hill) had lower average ribosomal RNA copy number of the bacteria present in heated compared to control soils (2.40 Control, 2.28 Heated, P < 0.05; Figure 4). Community structure did not change substantially when copy number was corrected based on Procrustes analysis of community profiles (Monte Carlo P < 0.05), with a high degree of similarity (R = 0.9511) between the copy number corrected and uncorrected PCoA community models (Figure S2). The most abundant taxa tended to have fewer than four copies of the 16S ribosomal RNA operon, while lower abundance taxa had as many as 14 estimated copies (Figure S3). Of the top 15 most dominant taxa, 13 had only one or two copies of the 16S ribosomal RNA operon. When copy number correction was applied, only 28 taxa changed rank such that they were no longer in the dominant subset community (comprised of 155 taxa).
Figure 4. Mean ribosomal RNA copy number was calculated for Prospect Hill (20 years warmed) soil communities, and these calculations were based on phylogenetic inference (see methods for details). Means are shown as as box (first and third quartiles) and whisker (95% CI) plots where the solid bar is median; statistical significance was determined based on ANOVA (P < 0.05).
Mounting evidence suggests that soil microbes play a role in elevated CO2 emissions and soil organic matter loss that is symptomatic of long-term warming in temperate forest ecosystems (Bardgett et al., 1999; Schimel and Gulledge, 2004; Frey et al., 2008, 2013). We set out to test the primary hypothesis that the soil bacterial community is altered by warming, and statistically significant differences were only apparent after 20 years. There are many factors in the environment that correlate with warming effects on microbial feedbacks to the climate system, including associated soil moisture and drought, changes in plant communities, and N deposition. For example, Blankinship and colleagues' meta-analysis of 75 manipulative climate change experiments found that warming was more likely to have a negative effect on microbial abundance (density) in cool, dry locations (Blankinship et al., 2011). Observed changes in beta diversity may be due to the loss of labile C (Frey et al., 2008) that represents a significant change in substrate available to resident microbes. Changes in substrate availability in soil are well known to affect changes in microbial community structure (Cleveland et al., 2007; Fierer et al., 2007; España et al., 2011). The extent to which changes in microbial substrate utilization will result in net changes in carbon cycle feedbacks to the atmosphere remain to be examined, and will necessitate an understanding of why communities are changing, as well as how populations' carbon use efficiencies are changing and the extent to which new soil carbon pools (essentially, new niche space) are being degraded by the changing microbial populations.
Our secondary hypothesis, that average bacterial copy number was decreased by long-term warming, was also supported by the ribosomal RNA copy number estimation evidence, showing that that long-term warming favors bacteria with an oligotrophic lifestyle (Klappenbach et al., 2000; Stevenson and Schmidt, 2004). Within a bacterial genome, the number of rRNA gene operon copies tends to correlate with maximum growth rate (Stevenson and Schmidt, 2004), the ability to change growth rates quickly (Klappenbach et al., 2000), and other traits including limited mobility and fewer types of more high-affinity transporters (Lauro et al., 2009), though there are exceptions to these generalizations (Blazewicz et al., 2013). Organisms with many copies of the rRNA gene operon are broadly considered to be copiotrophs, adapted for exploitation of varying and high-quality substrates, while those with single or few copies are considered to be oligotrophs, adapted to extract maximum resources out of a limited supply (Klappenbach et al., 2000; Stevenson and Schmidt, 2004). However, oligotrophy can also occur under conditions where privatization of resources is possible, e.g., conditions of high spatial structure or high heterogeneity (Pfeiffer et al., 2001; Stevenson and Schmidt, 2004; Lennon et al., 2012; Bachmann et al., 2013). It is possible that long-term warming has caused a change in soil structure that has increased the spatial heterogeneity, porosity, or other physical structure of the soils including physical or chemical protection of soil carbon that may contribute to the observed decrease in ribosomal operon copy number, though further study would be required to test this hypothesis.
We assume that the dominant subset of the community likely includes taxa that have important functions in the soil due to their success, though this is based on measures of relative abundance. Our measures of relative abundance based on QPCR suggest that of all bacteria, Actinobacteria comprise 8.8% in control and 9.4% in warmed, that Alphaproteobacteria comprise 25% in control and 44% in warmed soils, and that Acidobacteria comprise 3.4% in control and 7.4% in warmed soils (Table S4). However, the fraction of active populations at any one time may be as low as 10% for soil (Lennon and Jones, 2011), which confounds hypothesized links between function and observed community profiles by soil DNA evidence. Our original hypotheses were that long-term warming would induce a shift in the soil microbial communities, and that observed decreased soil C represented a decrease in microbial substrate availability that would increase the incidence of oligotrophy. These hypotheses are supported by our data, which include increased evenness with long-term warming, decreased ribosomal RNA copy number, increased community evenness and increased relative and absolute abundance of known or suspected oligotrophic taxa.
Ultimately, we are interested in changes in community structure insofar as they can reveal indicators of microbial feedbacks to climate, and because of this, we turned to indicator species and QPCR of key dynamic groups: Alphaproteobacteria, Acidobacteria, and Actinobacteria. Of the Alphaproteobacteria that changed with warming, the Rhizobiales mostly increased with warming while the Rhodospirillales mostly decreased with warming. The Rhizobiales include the most dominant Alphaproteobacteria, an unknown Hyphomicrobium spp. (Figure 1C), as well as members that are purple sulfur and non-sulfur bacteria, a group known for being able to grow under a wide range of conditions (Larimer et al., 2004). These taxa contain well-known plant root-associated microbes, though these experimental plots are too small to take into account changes in plant physiology or community, which also affect the observed microbial community in a warmer world (Bardgett et al., 1999). Acidobacteria can grow on complex polymers, including plant hemicellulose or cellulose and fungal chitin (Eichorst et al., 2011). The phyla Alphaproteobacteria, Acidobacteria and Actinobacteria contain many representative taxa known to degrade recalcitrant C and/or that have plant-specific associations (Barret et al., 2011). Further studies may elucidate the genetic or functional differences between the groups thus far represented only by the V4 region sequence and observed changes in relative abundance in a warmer world, though observation of genetic evolution in response to long-term warming or evolved functional changes, such as extracellular enzyme temperature optima, increased tolerance to low water potential conditions, or increased capacity for uptake and degradation of lower quality and quantity carbon will all require measures of physiology of isolated organisms in the lab.
Though absolute copies of Actinobacteria did not differ between warmed and control plots (Figure 3B), increased absolute abundance of Actinobacteria was observed at the Prospect Hill sites by FAME analysis after 12 years of warming (Frey et al., 2008). The most dominant taxon in the control plots, an Actinomycetales (class Actinobacteria, NewOTU5624), decreased 78% in response to warming, declining from 4.5 to 1.02% relative abundance (P < 0.05); the second most abundant taxon was also an Actinomycetales (NewOTU189491) and was unaffected by by warming (Figure 1C). In a separate study at Harvard Forest, Actinomycetes were shown to increase in relative abundance with the addition of labile, but not recalcitrant C source (Goldfarb et al., 2011), though these isolates were only incubated with lignin as recalcitrant C for 48 h. There are other studies that suggest Actinobacteria may be a rich reservoir of extracellular peroxidases including lignin peroxidases (Godden et al., 1992; Kirby, 2006), though further studies will determine how far lignin activity among the phylum Actinobacteria extends beyond the streptomycetes (Le Roes-Hill et al., 2011). We hypothesize that the first and second most dominant taxa may be representative of copiotrophic and oligotrophic groups, respectively: their respective estimated copy numbers support this (4.07 for NewOTU5624, and 2.75 for NewOTU189491). This would also explain the net zero change in Actinomycetes by Q-PCR. The reduced amount of labile C in the warmed soils (Bradford et al., 2008) would drive opposing responses of two phylogenetically similar but functionally divergent groups to warming–decreased relative abundance of the potential copiotroph and increase of the potential oligotroph.
Observed changes in fungal biomass were observed by FAME analysis after 12 years of warming at the 20-year warmed site, which showed that warming caused decreased fungal abundance in both the mineral and organic soil horizons (Frey et al., 2008). Reduced fungal biomass has also been observed with warming in other sites (Waldrop et al., 2004; Rinnan et al., 2007; Frey et al., 2008) and though this current study focused mainly on bacteria, the importance of fungi in this systems is being separately studied. Fungi generally dominate primary decomposition in the organic horizon of temperate soils (Berg et al., 1998; Thevenot et al., 2010), and fungal laccase, phenol oxidase and peroxidase activities and genes encoding these enzymes have been found in greater abundance in upper layer, high organic-matter content soils relative to deeper, mineral soils (Luis et al., 2005; Sinsabaugh, 2010). The significant loss in organic horizon soil carbon with long-term warming is likely a combination of increased activity of primary decomposers (generally fungi) in the organic horizon, and increased demand (either through activity or abundance) of secondary decomposers (generally bacteria) in the mineral horizon. Like fungi, Actinomycetes are usually filamentous, and their dominance could suggest changing niches in soil and different contributions to soil C cycling due to long-term warming (Six et al., 2006; De Boer et al., 2008, 2005). Functional analyses of organisms from the organic and mineral horizons separately will suggest mechanisms for the observed substantially depleted soil organic layer mass with warming.
While a community shift was only observed after 20 years, functional changes and thermal acclimation were observed much earlier (Melillo et al., 2002; Bradford et al., 2008; Frey et al., 2013). Shorter term studies have also noticed no change in microbial community by rRNA gene profiles with warming despite increases in respiration, such as temperate mountain forest soil warmed 4°C for 4 years (Kuffner et al., 2012) and mature spruce forest soils warmed 4°C during snow-free seasons for 4 years (Schindlbacher et al., 2011), which is consistent with our observations of the 5-year warmed SWaN plots. Genomic studies such as this permit us to determine whether a change in the community accompanies such changes in function, where no change in community would suggest a high degree of functional redundancy and a capacity for soil community resilience despite the warming treatment (Shade et al., 2012). For example, the KFFL field grassland soils are also dominated by Actinobacteria, Alphaproteobacteria, and Acidobacteria, as are the old-field soils studied by Castro and colleagues, and in both cases, warming alone had a smaller effect than when it was studied in conjunction with elevated CO2 (Castro et al., 2010) or with drought (Sheik et al., 2011). This study is valuable in that it permits examination of the long-term effects of warming without confounding environmental and climate factors. The change in community structure that we observed after 20 years suggests that some other aspect of the soil niche space caused pressure on the community, resulting in the change (Schimel et al., 2007). While in some instances a strong positive correlation has been observed between taxonomic and functional richness of bacteria (Konstantinidis and Tiedje, 2007; Richter and Rosselló-Móra, 2009), we expect that expanded functional diversity that has been observed in these soils (Frey et al., 2013) more likely has its origins in changing active populations or genetic adaptation.
The long-term effects of warming, including altered community structure, decreased fungal biomass, increased evenness and decreased ribosomal copy number as an indicator of oligotrophy, all suggest positive warming-induced climate feedbacks. Indications of increased oligotrophy are perhaps the most alarming, because these may suggest that long-term warming is causing decreased physical protection of older or more recalcitrant soil C pools, assuming that physically protected C pools are more recalcitrant and that accessing these pools requires strategies (different enzymes, exploratory growth) consistent with oligotrophy (Schmidt et al., 2011; Bödeker et al., 2014). Testing hypotheses of mechanisms of how warming changes microbial communities, by distinguishing increased enzyme activity from altered enzyme production (Conant et al., 2011), as well as through decreased soil moisture (Toberman et al., 2008; Peñuelas et al., 2012), would benefit from accompanying direct observations on microbial activity, including comparative microbial physiology and genomics of isolated dominant strains. Genomic studies will remain valuable for understanding genomic context for changes in function, and should ultimately enable incorporation of microbial parameters into modeling efforts for prediction of microbial feedbacks to changing climate in a warmer world.
Our data support the main hypothesis that long-term warming induces changes in microbial community composition, changes that are not seen in intermediate lengths of time. The strong community unevenness and similarity between dominant and whole community beta diversity suggests that a few keystone species may be responsible for a large proportion of the soil C cycling activity. The decreased abundance of dominant bacterial taxa as well as of total fungi (Frey et al., 2008) supports our second hypothesis that there is shifting niche space that may be evidence of changing C availability. The reduced ability of fungi and some Actinobacteria to survive may be due to dwindling resources, shifts in C quality, or a reduction in fine root biomass due to long-term warming, and may have created an opportunity for other oligotrophic bacteria (Butler et al., 2012; Koranda et al., 2013). Understanding the specific contributions of the dominant taxa to soil C cycling would be a powerful tool for modeling the relationship between microbial diversity and changes in climate.
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.
The authors wish to gratefully acknowledge Melissa Knorr, Brian Godbois, and Chris Cook at the University of New Hampshire for soil sample collection; Zach Herbert at the DFCI Molecular Biology Core Facility for providing technical assistance with sequencing; and Erin Nuccio at Lawrence Livermore National Laboratory for providing the Greengenes reference files truncated to the V4 region for OTU calling. This work was supported by funding from the University of Massachusetts Amherst to DeAngelis and the National Science Foundation Long-term Ecological Research (LTER) Program.
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fmicb.2015.00104/abstract
Andrews, J. A., Matamala, R., Westover, K. M., and Schlesinger, W. H. (2000). Temperature effects on the diversity of soil heterotrophs and the δ 13C of soil-respired CO2. Soil Biol. Biochem. 32, 699–706. doi: 10.1016/S0038-0717(99)00206-0
Andrews, S. (2010). FASTQC. A Quality Control Tool for High Throughput Sequence Data. Available online at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc
Arndt, D., Xia, J., Liu, Y., Zhou, Y., Guo, A. C., Cruz, J. A., et al. (2012). METAGENassist: a comprehensive web server for comparative metagenomics. Nucleic Acids Res. 40, W88–W95. doi: 10.1093/nar/gks497
Bacchetti De Gregoris, T., Aldred, N., Clare, A. S., and Burgess, J. G. (2011). Improvement of phylum- and class-specific primers for real-time PCR quantification of bacterial taxa. J. Microbiol. Methods 86, 351–356. doi: 10.1016/j.mimet.2011.06.010
Bachmann, H., Fischlechner, M., Rabbers, I., Barfa, N., Santos, F. B., dos, Molenaar, D., et al. (2013). Availability of public goods shapes the evolution of competing metabolic strategies. Proc. Natl. Acad. Sci. U.S.A. 110, 14302–14307. doi: 10.1073/pnas.1308523110
Bardgett, R. D., Kandeler, E., Tscherko, D., Hobbs, P. J., Bezemer, T. M., Jones, T. H., et al. (1999). Below-ground microbial community development in a high temperature world. Oikos 85, 193. doi: 10.2307/3546486
Barret, M., Morrissey, J., and O'Gara, F. (2011). Functional genomics analysis of plant growth-promoting rhizobacterial traits involved in rhizosphere competence. Biol. Fertil. Soils 47, 729–743. doi: 10.1007/s00374-011-0605-x
Belay-Tedla, A., Zhou, X., Su, B., Wan, S., and Luo, Y. (2009). Labile, recalcitrant, and microbial carbon and nitrogen pools of a tallgrass prairie soil in the US Great Plains subjected to experimental warming and clipping. Soil Biol. Biochem. 41, 110–116. doi: 10.1016/j.soilbio.2008.10.003
Berg, M. P., Kniese, J. P., and Verhoef, H. A. (1998). Dynamics and stratification of bacteria and fungi in the organic layers of a scots pine forest soil. Biol. Fertil. Soils 26, 313–322. doi: 10.1007/s003740050382
Blazewicz, S. J., Barnard, R. L., Daly, R. A., and Firestone, M. K. (2013). Evaluating rRNA as an indicator of microbial activity in environmental communities: limitations and uses. ISME J. 7, 2061–2068. doi: 10.1038/ismej.2013.102
Bödeker, I. T. M., Clemmensen, K. E., de Boer, W., Martin, F., Olson, Å., and Lindahl, B. D. (2014). Ectomycorrhizal cortinarius species participate in enzymatic oxidation of humus in northern forest ecosystems. New Phytol. 203, 245–256. doi: 10.1111/nph.12791
Bradford, M. A., Davies, C. A., Frey, S. D., Maddox, T. R., Melillo, J. M., Mohan, J. E., et al. (2008). Thermal adaptation of soil microbial respiration to elevated temperature. Ecol. Lett. 11, 1316–1327. doi: 10.1111/j.1461-0248.2008.01251.x
Brzostek, E. R., Blair, J. M., Dukes, J. S., Frey, S. D., Hobbie, S. E., Melillo, J. M., et al. (2012). The effect of experimental warming and precipitation change on proteolytic enzyme activity: positive feedbacks to nitrogen availability are not universal. Glob. Change Biol. 18, 2617–2625. doi: 10.1111/j.1365-2486.2012.02685.x
Butler, S. M., Melillo, J. M., Johnson, J., Mohan, J., Steudler, P. A., Lux, H., et al. (2012). Soil warming changes the nitrogen cycle: nitrogen-driven changes in species composition. Oecologia 168, 819–828. doi: 10.1007/s00442-011-2133-7
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. Meth. 7, 335–336. doi: 10.1038/nmeth.f.303
Caporaso, J. G., Lauber, C. L., Walters, W. A., Berg-Lyons, D., Huntley, J., Fierer, N., et al. (2012). Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 6, 1621–1624. doi: 10.1038/ismej.2012.8
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-09
Cleveland, C., Nemergut, D., Schmidt, S., and Townsend, A. (2007). Increases in soil respiration following labile carbon additions linked to rapid shifts in soil microbial community composition. Biogeochemistry 82, 229–240. doi: 10.1007/s10533-006-9065-z
Conant, R. T., Ryan, M. G., Ågren, G. I., Birge, H. E., Davidson, E. A., Eliasson, P. E., et al. (2011). Temperature and soil organic matter decomposition rates – synthesis of current knowledge and a way forward. Glob. Change Biol. 17, 3392–3404. doi: 10.1111/j.1365-2486.2011.02496.x
Contosta, A. R., Frey, S. D., and Cooper, A. B. (2011). Seasonal dynamics of soil respiration and N mineralization in chronically warmed and fertilized soils. Ecosphere 2:art36. doi: 10.1890/ES10-00133.1
Davidson, E. A., Samanta, S., Caramori, S. S., and Savage, K. (2012). The Dual Arrhenius and Michaelis–Menten kinetics model for decomposition of soil organic matter at hourly to seasonal time scales. Glob. Change Biol. 18, 371–384. doi: 10.1111/j.1365-2486.2011.02546.x
DeAngelis, K. M., Silver, W. L., Thompson, A. W., and Firestone, M. K. (2010). Microbial communities acclimate to recurring changes in soil redox potential status. Environ. Microbiol. 12, 3137–3149. doi: 10.1111/j.1462-2920.2010.02286.x
De Boer, W., de Ridder-Duine, A. S., Klein Gunnewiek, P. J. A., Smant, W., and Van Veen, J. A. (2008). Rhizosphere bacteria from sites with higher fungal densities exhibit greater levels of potential antifungal properties. Soil Biol. Biochem. 40, 1542–1544. doi: 10.1016/j.soilbio.2007.12.030
De Boer, W., Folman, L. B., Summerbell, R. C., and Boddy, L. (2005). Living in a fungal world: impact of fungi on soil bacterial niche development. FEMS Microbiol. Rev. 29, 795–811. doi: 10.1016/j.femsre.2004.11.005
DeSantis, T. Z., Hugenholtz, P., Larsen, N., Rojas, M., Brodie, E. L., Keller, K., et al. (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 72, 5069–5072. doi: 10.1128/AEM.03006-05
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
Eichorst, S. A., Kuske, C. R., and Schmidt, T. M. (2011). Influence of plant polymers on the distribution and cultivation of bacteria in the phylum Acidobacteria. Appl. Environ. Microbiol. 77, 586–596. doi: 10.1128/AEM.01080-10
España, M., Rasche, F., Kandeler, E., Brune, T., Rodriguez, B., Bending, G. D., et al. (2011). Identification of active bacteria involved in decomposition of complex maize and soybean residues in a tropical Vertisol using 15N-DNA stable isotope probing. Pedobiologia 54, 187–193. doi: 10.1016/j.pedobi.2011.03.001
Fierer, N., Jackson, J. A., Vilgalys, R., and Jackson, R. B. (2005). Assessment of soil microbial community structure by use of taxon-specific quantitative PCR assays. Appl. Environ. Microbiol. 71, 4117–4120. doi: 10.1128/AEM.71.7.4117-4120.2005
Frey, S. D., Drijber, R., Smith, H., and Melillo, J. (2008). Microbial biomass, functional capacity, and community structure after 12 years of soil warming. Soil Biol. Biochem. 40, 2904–2907. doi: 10.1016/j.soilbio.2008.07.020
Godden, B., Ball, A. S., Helvenstein, P., Mccarthy, A. J., and Penninckx, M. J. (1992). Towards elucidation of the lignin degradation pathway in actinomycetes. J. Gen. Microbiol. 138, 2441–2448. doi: 10.1099/00221287-138-11-2441
Goldfarb, K. C., Karaoz, U., Hanson, C. A., Santee, C. A., Bradford, M. A., Treseder, K. K., et al. (2011). Differential growth responses of soil bacterial taxa to carbon substrates of varying chemical recalcitrance. Front. Microbiol. 2:94. doi: 10.3389/fmicb.2011.00094
Griffiths, R. I., Whiteley, A. S., O'Donnell, A. G., and Bailey, M. J. (2000). Rapid method for coextraction of DNA and RNA from natural environments for analysis of ribosomal DNA- and rRNA-based microbial community composition. Appl. Environ. Microbiol. 66, 5488–5491. doi: 10.1128/AEM.66.12.5488-5491.2000
Hamady, M., Lozupone, C., and Knight, R. (2009). Fast UniFrac: facilitating high-throughput phylogenetic analyses of microbial communities including analysis of pyrosequencing and PhyloChip data. ISME J. 4, 17–27. doi: 10.1038/ismej.2009.97
Hold, G. L., Pryde, S. E., Russell, V. J., Furrie, E., and Flint, H. J. (2002). Assessment of microbial diversity in human colonic samples by 16S rDNA sequence analysis. FEMS Microbiol. Ecol. 39, 33–39. doi: 10.1111/j.1574-6941.2002.tb00904.x
IPCC. (2013). “Climate Change 2013: The Physical Science Basis,” in Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, eds T. F. Stocker, D. Qin, G.-K. Plattner, M. Tignor, S. K. Allen, J. Boschung, et al. (Cambridge;New York, NY: Cambridge University Press), 1535.
Kembel, S. W., Wu, M., Eisen, J. A., and Green, J. L. (2012). Incorporating 16S Gene Copy Number Information Improves Estimates of Microbial Diversity and Abundance. PLoS Comput. Biol. 8:e1002743. doi: 10.1371/journal.pcbi.1002743
Klappenbach, J. A., Dunbar, J. M., and Schmidt, T. M. (2000). rRNA operon copy number reflects ecological strategies of bacteria. Appl. Environ. Microbiol. 66, 1328–1333. doi: 10.1128/AEM.66.4.1328-1333.2000
Konstantinidis, K. T., and Tiedje, J. M. (2007). Prokaryotic taxonomy and phylogeny in the genomic era: advancements and challenges ahead. Curr. Opin. Microbiol. 10, 504–509. doi: 10.1016/j.mib.2007.08.006
Koranda, M., Kaiser, C., Fuchslueger, L., Kitzler, B., Sessitsch, A., Zechmeister-Boltenstern, S., et al. (2013). Seasonal variation in functional properties of microbial communities in beech forest soil. Soil Biol. Biochem. 60, 95–104. doi: 10.1016/j.soilbio.2013.01.025
Kuffner, M., Hai, B., Rattei, T., Melodelima, C., Schloter, M., Zechmeister-Boltenstern, S., et al. (2012). Effects of season and experimental warming on the bacterial community in a temperate mountain forest soil assessed by 16S rRNA gene pyrosequencing. FEMS Microbiol. Ecol. 82, 551–562. doi: 10.1111/j.1574-6941.2012.01420.x
Kuhn, M. (2008). Building Predictive Models in R Using the caret Package. J. Stat. Softw. 28. Avaliable online at: http://econpapers.repec.org/article/jssjstsof/28_3ai05.htm
Larimer, F. W., Chain, P., Hauser, L., Lamerdin, J., Malfatti, S., Do, L., et al. (2004). Complete genome sequence of the metabolically versatile photosynthetic bacterium Rhodopseudomonas palustris. Nat. Biotech. 22, 55–61. doi: 10.1038/nbt923
Lauro, F. M., McDougald, D., Thomas, T., Williams, T. J., Egan, S., Rice, S., et al. (2009). The genomic basis of trophic strategy in marine bacteria. Proc. Natl. Acad. Sci. U.S.A. 106, 15527–15533. doi: 10.1073/pnas.0903507106
Lennon, J. T., Aanderud, Z. T., Lehmkuhl, B. K., and Schoolmaster, D. R. (2012). Mapping the niche space of soil microorganisms using taxonomy and traits. Ecology 93, 1867–1879. doi: 10.1890/11-1745.1
Li, J., Wang, G., Allison, S. D., Mayes, M. A., and Luo, Y. (2014). Soil carbon sensitivity to temperature and carbon use efficiency compared across microbial-ecosystem models of varying complexity. Biogeochemistry 119, 67–84. doi: 10.1007/s10533-013-9948-8
Luis, P., Kellner, H., Zimdars, B., Langer, U., Martin, F., and Buscot, F. (2005). Patchiness and spatial distribution of laccase genes of ectomycorrhizal, saprotrophic, and unknown basidiomycetes in the upper horizons of a mixed forest cambisol. Microb. Ecol. 50, 570–579. doi: 10.1007/s00248-005-5047-2
Luo, C., Rodriguez, -R. L. M., Johnston, E. R., Wu, L., Cheng, L., Xue, K., et al. (2014). Soil Microbial Community Responses to a Decade of Warming as Revealed by Comparative Metagenomics. Appl. Environ. Microbiol. 80, 1777–1786. doi: 10.1128/AEM.03712-13
Matsen, F. A., Kodner, R. B., and Armbrust, E. V. (2010). pplacer: linear time maximum-likelihood and Bayesian phylogenetic placement of sequences onto a fixed reference tree. BMC Bioinformatics 11:538. doi: 10.1186/1471-2105-11-538
Melillo, J. M., Butler, S., Johnson, J., Mohan, J., Steudler, P., Lux, H., et al. (2011). Soil warming, carbon–nitrogen interactions, and forest carbon budgets. Proc. Natl. Acad. Sci. U.S.A. 108, 9508–9512. doi: 10.1073/pnas.1018189108
Melillo, J. M., Callaghan, T. V., Woodward, F. I., Salati, E., and Sinha, S. K. (1990). Effects on the ecosystems. Climate Change—The IPCC Scientific Assessment Cambridge (UK);Cambridge University Press. Avaliable online at: http://agris.fao.org/agris-search/search/display.do?f=1993/ES/ES93008.xml;ES9202067 (Accessed April 4, 2012).
Melillo, J. M., Steudler, P. A., Aber, J. D., Newkirk, K., Lux, H., Bowles, F. P., et al. (2002). Soil warming and carbon-cycle feedbacks to the climate system. Science 298, 2173–2176. doi: 10.1126/science.1074153
Mevik, B.-H., and Wehrens, R. (2007). The pls package: principal component and partial least squares regression in R. J. Stat. Softw. 18. Avaliable online at: http://econpapers.repec.org/article/jssjstsof/18_3Ai02.htm
Muyzer, G., de Waal, E. C., and Uitterlinden, A. G. (1993). Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA. Appl. Environ. Microbiol. 59, 695–700.
Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O'Hara, R. B., et al. (2011). Vegan: Comunity Ecology Package. Avaliable online at: http://CRAN.R-project.org/package=vegan
Peñuelas, J., Rico, L., Ogaya, R., Jump, A. S., and Terradas, J. (2012). Summer season and long-term drought increase the richness of bacteria and fungi in the foliar phyllosphere of Quercus ilex in a mixed Mediterranean forest. Plant Biol. 14, 565–575. doi: 10.1111/j.1438-8677.2011.00532.x
Peterjohn, W. T., Melillo, J. M., Steudler, P. A., Newkirk, K. M., Bowles, F. P., and Aber, J. D. (1994). Responses of trace gas fluxes and N availability to experimentally elevated soil temperatures. Ecol. Appl. 4, 617–625. doi: 10.2307/1941962
Pietikåinen, J., Pettersson, M., and Bååth, E. (2005). Comparison of temperature effects on soil respiration and bacterial and fungal growth rates. FEMS Microbiol. Ecol. 52, 49–58. doi: 10.1016/j.femsec.2004.10.002
Pold, G., and DeAngelis, K. (2013). Up against the wall: the effects of climate warming on soil microbial diversity and the potential for feedbacks to the carbon cycle. Diversity 5, 409–425. doi: 10.3390/d5020409
R Core Team. (2014). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing. Available online at: http://www.R-project.org/.
Rinnan, R., Michelsen, A., Bååth, E., and Jonasson, S. (2007). Fifteen years of climate change manipulations alter soil microbial communities in a subarctic heath ecosystem. Glob. Change Biol. 13, 28–39. doi: 10.1111/j.1365-2486.2006.01263.x
Roesch, L. F. W., Fulthorpe, R. R., Riva, A., Casella, G., Hadwin, A. K. M., Kent, A. D., et al. (2007). Pyrosequencing enumerates and contrasts soil microbial diversity. ISME J. 1, 283–290. doi: 10.1038/ismej.2007.53
RStudio. (2012). RStudio: Integrated development environment for R (Version 0.96.122). Boston, MA Avaliable online at: http://www.rstudio.org/.
Rustad, L. R., Campbell, J. C., Marion, G. M., Norby, R. N., Mitchell, M. M., Hartley, A. H., et al. (2001). A meta-analysis of the response of soil respiration, net nitrogen mineralization, and aboveground plant growth to experimental ecosystem warming. Oecologia 126, 543–562. doi: 10.1007/s004420000544
Schindlbacher, A., Rodler, A., Kuffner, M., Kitzler, B., Sessitsch, A., and Zechmeister-Boltenstern, S. (2011). Experimental warming effects on the microbial community of a temperate mountain forest soil. Soil Biol. Biochem. 43, 1417–1425. doi: 10.1016/j.soilbio.2011.03.005
Schmidt, M. W. I., Torn, M. S., Abiven, S., Dittmar, T., Guggenberger, G., Janssens, I. A., et al. (2011). Persistence of soil organic matter as an ecosystem property. Nature 478, 49–56. doi: 10.1038/nature10386
Shade, A., Peter, H., Allison, S. D., Baho, D. L., Berga, M., Burgmann, H., et al. (2012). Fundamentals of Microbial Community Resistance and Resilience. Front. Microbiol. 3:417 doi: 10.3389/fmicb.2012.00417
Sheik, C. S., Beasley, W. H., Elshahed, M. S., Zhou, X., Luo, Y., and Krumholz, L. R. (2011). Effect of warming and drought on grassland microbial communities. ISME J. 5, 1692–1700. doi: 10.1038/ismej.2011.32
Smyth, G. (2004). “limma: linear models for microarray Data,” in Bioinformatics and Computational Biology Solutions Using R and Bioconductor, (New York, NY: Springer), 397–420. doi: 10.1007/0-387-29362-0_23
Suau, A., Bonnet, R., Sutren, M., Godon, J.-J., Gibson, G. R., Collins, M. D., et al. (1999). Direct analysis of genes encoding 16s rRNA from complex communities reveals many novel molecular species within the human gut. Appl. Environ. Microbiol. 65, 4799–4807.
Tibshirani, R., Hastie, T., Narasimhan, B., and Chu, G. (2003). Class prediction by nearest shrunken centroids, with applications to DNA microarrays. Statist. Sci. 18, 104–117. doi: 10.1214/ss/1056397488
Toberman, H., Freeman, C., Evans, C., Fenner, N., and Artz, R. R. E. (2008). Summer drought decreases soil fungal diversity and associated phenol oxidase activity in upland Calluna heathland soil. FEMS Microbiol. Ecol. 66, 426–436. doi: 10.1111/j.1574-6941.2008.00560.x
Treseder, K., Balser, T., Bradford, M., Brodie, E., Dubinsky, E., Eviner, V., et al. (2012). Integrating microbial ecology into ecosystem models: challenges and priorities. Biogeochemistry 109, 7–18. doi: 10.1007/s10533-011-9636-5
Větrovský, T., and Baldrian, P. (2013). The variability of the 16s rRNA gene in bacterial genomes and its consequences for bacterial community analyses. PLoS ONE 8:e57923. doi: 10.1371/journal.pone.0057923
Waldrop, M. P., Zak, D. R., and Sinsabaugh, R. L. (2004). Microbial community response to nitrogen deposition in northern forest ecosystems. Soil Biol. Biochem. 36, 1443–1451. doi: 10.1016/j.soilbio.2004.04.023
Werner, J. J., Koren, O., Hugenholtz, P., DeSantis, T. Z., Walters, W. A., Caporaso, J. G., et al. (2012). Impact of training sets on classification of high-throughput bacterial 16s rRNA gene surveys. ISME J. 6, 94–103. doi: 10.1038/ismej.2011.82
Wickham, H. (2012). Reshaping data with the reshape package. J. Stat. Softw. 21. Avaliable online at: http://econpapers.repec.org/article/jssjstsof/21_3ai12.htm
Wright, E. S., Yilmaz, L. S., and Noguera, D. R. (2012). DECIPHER, a search-based approach to chimera identification for 16S rRNA sequences. Appl. Environ. Microbiol. 78, 717–725. doi: 10.1128/AEM.06516-11
Wu, Z., Dijkstra, P., Koch, G. W., Peñuelas, J., and Hungate, B. A. (2011). Responses of terrestrial ecosystems to temperature and precipitation change: a meta-analysis of experimental manipulation. Glob. Change Biol. 17, 927–942. doi: 10.1111/j.1365-2486.2010.02302.x
Zhang, W., Parker, K. M., Luo, Y., Wan, S., Wallace, L. L., and Hu, S. (2005). Soil microbial responses to experimental warming and clipping in a tallgrass prairie. Glob. Change Biol. 11, 266–277. doi: 10.1111/j.1365-2486.2005.00902.x
Zogg, G. P., Zak, D. R., Ringelberg, D. B., White, D. C., MacDonald, N. W., and Pregitzer, K. S. (1997). Compositional and Functional Shifts in Microbial Communities Due to Soil Warming. Soil Sci. Soc. Am. J. 61, 475–481. doi: 10.2136/sssaj1997.03615995006100020015x
Keywords: climate change, microbial ecology, ribosomal RNA, rrn operon copy number, trophic strategy
Citation: DeAngelis KM, Pold G, Topçuoğlu BD, van Diepen LTA, Varney RM, Blanchard JL, Melillo J and Frey SD (2015) Long-term forest soil warming alters microbial communities in temperate forest soils. Front. Microbiol. 6:104. doi: 10.3389/fmicb.2015.00104
Received: 04 October 2014; Accepted: 27 January 2015;
Published online: 13 February 2015.
Edited by:Stuart Findlay, Cary Institute of Ecosystem Studies, USA
Reviewed by:Kirsten Hofmockel, Iowa State University, USA
Hongchen Jiang, China University of Geosciences, Wuhan, China
Copyright © 2015 DeAngelis, Pold, Topçuoğlu, van Diepen, Varney, Blanchard, Melillo and Frey. 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) or licensor 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: Kristen M. DeAngelis, Department of Microbiology, University of Massachusetts, 639 North Pleasant St., 203 Morrill IVN, Amherst, MA 01002, USA e-mail: firstname.lastname@example.org