Original Research ARTICLE
Soil Conditions Rather Than Long-Term Exposure to Elevated CO2 Affect Soil Microbial Communities Associated with N-Cycling
- 1Department of Biogeochemistry, Max Planck Institute for Terrestrial Microbiology, Marburg, Germany
- 2Department of Plant Ecology, University of Giessen, Giessen, Germany
- 3Water Resources and Environmental Engineering Research Unit, University of Oulu, Oulu, Finland
- 4Department of Ecological Microbiology, University of Bayreuth, Bayreuth, Germany
- 5Institute of Microbiology, Leibniz Universität Hannover, Hannover, Germany
- 6Climate Change Research for Special Crops, Department of Soil Science and Plant Nutrition, Geisenheim University, Geisenheim, Germany
- 7School of Biology and Environmental Science, University College Dublin, Dublin, Ireland
- 8University of Kiel, Kiel, Germany
Continuously rising atmospheric CO2 concentrations may lead to an increased transfer of organic C from plants to the soil through rhizodeposition and may affect the interaction between the C- and N-cycle. For instance, fumigation of soils with elevated CO2 (eCO2) concentrations (20% higher compared to current atmospheric concentrations) at the Giessen Free-Air Carbon Dioxide Enrichment (GiFACE) sites resulted in a more than 2-fold increase of long-term N2O emissions and an increase in dissimilatory reduction of nitrate compared to ambient CO2 (aCO2). We hypothesized that the observed differences in soil functioning were based on differences in the abundance and composition of microbial communities in general and especially of those which are responsible for N-transformations in soil. We also expected eCO2 effects on soil parameters, such as on nitrate as previously reported. To explore the impact of long-term eCO2 on soil microbial communities, we applied a molecular approach (qPCR, T-RFLP, and 454 pyrosequencing). Microbial groups were analyzed in soil of three sets of two FACE plots (three replicate samples from each plot), which were fumigated with eCO2 and aCO2, respectively. N-fixers, denitrifiers, archaeal and bacterial ammonia oxidizers, and dissimilatory nitrate reducers producing ammonia were targeted by analysis of functional marker genes, and the overall archaeal community by 16S rRNA genes. Remarkably, soil parameters as well as the abundance and composition of microbial communities in the top soil under eCO2 differed only slightly from soil under aCO2. Wherever differences in microbial community abundance and composition were detected, they were not linked to CO2 level but rather determined by differences in soil parameters (e.g., soil moisture content) due to the localization of the GiFACE sets in the experimental field. We concluded that +20% eCO2 had little to no effect on the overall microbial community involved in N-cycling in the soil but that spatial heterogeneity over extended periods had shaped microbial communities at particular sites in the field. Hence, microbial community composition and abundance alone cannot explain the functional differences leading to higher N2O emissions under eCO2 and future studies should aim at exploring the active members of the soil microbial community.
Due to anthropogenic emissions, atmospheric CO2 concentrations are rising by about 1% per year and are expected to double in this century (IPCC, 2014) causing well known-climatic effects. Observations from the world-wide longest lasting CO2 enrichment study, the Giessen Free Air Carbon Dioxide Enrichment (GiFACE since 1998, ongoing), showed that elevated atmospheric CO2 (eCO2) concentrations also exert several impacts on soil communities. For instance, plant biomass was stimulated by ~12–15% (Kammann et al., 2005; Andresen et al., 2017). It is thus hypothesized that an increased transfer of organic C from plants to the soil through rhizodeposition occurs which affects soil microbial communities with implications for the interaction between C- and N-cycling (Freeman et al., 2004; Denef et al., 2007). A meta-analysis of greenhouse gas emission data from CO2 enrichment experiments demonstrated that increased CO2 generally stimulated emissions of nitrous oxide (N2O), another potent greenhouse gas, from terrestrial ecosystems (van Groenigen et al., 2011). At GiFACE for instance, long-term N2O emissions under eCO2 increased more than 2-fold compared to ambient CO2 (aCO2), but the underlying mechanisms are not fully resolved yet (Kammann et al., 2008). In two 15N tracing laboratory experiments with soils from FACE sites in Giessen and New Zealand gross N-transformations under eCO2 shifted toward a higher importance of heterotrophic processes (Müller et al., 2009; Rütting et al., 2010). In addition, turnover of ammonia (heterotrophic nitrification) and the rates of dissimilatory reduction of nitrate to ammonia (DNRA) increased, while turnover of nitrate was reduced. At GiFACE, ammonia concentrations under eCO2 were on average 17% higher while nitrate concentrations were significantly lower than at ambient CO2 (Müller et al., 2009).
Changes in gross N-transformations and gaseous N emissions are dependent on the dynamics and activity of microbial communities. In soils, N2O is mainly produced by denitrifiers and nitrifiers (Conrad, 1996; Butterbach-Bahl et al., 2013) and alterations in the functioning of denitrifiers and ammonia oxidizers in soils exposed to eCO2 were clearly discernable (e.g., Barnard et al., 2005, 2006). However, little information is available to date on how these functional shifts may be related to shifts in the underlying microbial communities and the understanding of potential feedback effects resulting in higher N2O emissions is still limited. Several studies found profound differences in abundance and composition between the overall microbial communities in soils exposed to elevated and ambient CO2 (Denef et al., 2007; Drigo et al., 2008, 2009; He et al., 2010; Deng et al., 2012; Xu et al., 2013; Dunbar et al., 2014; Xiong et al., 2015; Xia et al., 2017). Elevated levels of CO2 were also reported to influence microbial communities associated with N-cycling. Two early cultivation based studies showed an enhanced abundance of nitrate dissimilating Pseudomonas in the rhizosphere of grasses at eCO2 (Fromin et al., 2005; Roussel-Delif et al., 2005). Lesaulnier et al. (2008) found a significant decrease of nitrate reducers and crenarchaeal ammonia oxidizers with eCO2 and field exposure of a grassland ecosystem to eCO2 for 10 years significantly increased the abundance of N-fixers and nirS-type denitrifiers (He et al., 2010). Interestingly, in two out of three replicate FACE plots studied at GiFACE, the ratio of N2O reducers to nitrite reducers was lower under eCO2 (Regan et al., 2011) and may thus explain higher N2O fluxes from the soil (Philippot et al., 2011). In a California grassland, the structure and abundance of the ammonia oxidizing bacterial community was altered by eCO2, strongly interacting with the factor precipitation (Horz et al., 2004). Horz et al. (2004) also showed that multifactorial global change (eCO2, temperature, precipitation, N-deposition) fed back into the enrichment of a specific clade of ammonia oxidizers related to Nitrosospira spp. with higher potential for nitrification.
Alteration of microbial communities associated with soil functioning such as N-cycling suggests concomitant alterations of potential functional activity and hence of ecosystem functioning (He et al., 2010). We hypothesized that the increased N2O emissions in response to eCO2 (Kammann et al., 2008) and shifts in N-transformations observed at Giessen FACE during long-term exposure to elevated CO2 levels (Müller et al., 2009; Rütting et al., 2010) lead to differences in soil parameters (e.g., soil nitrate) and can be explained by changes of microbial communities in general and particularly of communities associated with N-cycling. To explore the microbial communities, we applied a molecular approach to study the abundance and composition based on functional marker genes for denitrification (nirK/nirS, nosZ), ammonia oxidation (bacterial and archaeal amoA), nitrogen fixation (nifH), dissimilatory nitrate reduction to ammonia (DNRA, nrfA) as well as archaeal and bacterial communities (16S rRNA genes). We used different molecular techniques (qPCR, T-RFLP, and 454 pyrosequencing), of which each satisfies specific demands in microbial community analyses. In addition to bias introduced by the use of gene specific primers, each technique, however, has its limitations but their complementary results have the potential to provide more comprehensive insights. Screening techniques (qPCR and T-RFLP) allow for a comparative assessment of microbial community abundance and composition in high numbers of replicate samples. An appropriate level of replication is a prerequisite for downstream statistical analyses of the data sets. While these approaches do not provide information on the identity of organisms, thousands of sequences generated through 454 pyrosequencing on the other hand, allow separating genotypes of functional marker genes and thus provide an in depth analysis of community composition.
Materials and Methods
Site Description and Sampling
Soil samples were taken from the GiFACE experiment site (50°32′N and 8°43.3′E; 172 m a.s.l.) near Giessen, Germany. Within the GiFACE experiment CO2 fumigation on an old grassland site (>100 years) was started in May 1998 to study the response of a semi-natural grassland to long-term, moderate atmospheric CO2 enrichment of 20% above ambient (Jäger et al., 2003). The whole facility consists of six circular plots, each with 8 m internal diameter. Two plots build one set each (numbered 1, 2, and 3) with an ambient (aCO2, labeled with A) and an elevated (eCO2, labeled with E) CO2 plot. In 1998, the ambient CO2 concentration in the aCO2 plots was 364 and 399 ppm in eCO2. In 2012, when this study was conducted, the CO2 concentrations were 390 and 422 ppm, respectively. The three sets of plots are located on a gradual terrain slope (2°) in the direction of the rivulet Lückebach, which causes a gradient in soil moisture during spring and summer, and therefore drier conditions in set 1 compared to the others. The evapotranspiration during the plant growth period differed between aCO2 and eCO2 plots and resulted in soil moisture differences (Figure S2). The soil in the GiFACE plots was classified as a Fluvic Gleysol and has a sandy clay loam texture on top of a clay layer that varies in depth between the three sets (Jäger et al., 2003). The soil was characterized by a mean C and N content of 4.5 and 0.45%, respectively, and had a pH of ~6.2. Mean annual precipitation was 550 mm and mean annual air temperature was 9.6°C during the observation period from 1996 to 2003. Vegetation is the same in all plots and is dominated by 12 grass species, 2 legumes, and 15 non-leguminous herbs, and is characterized for all six plots as an Arrhenatheretum elatioris Br. Bl. Filipendula ulmaria subcommunity. The grassland has not been plowed for at least 100 years. It has been managed for several decades as a hay meadow with two cuts per year, and fertilized in mid-April with granular mineral calcium-ammonium-nitrate fertilizer at the rate of 40 kg N ha−1 yr−1 since 1996; before 1996, it was fertilized at 50–100 kg N ha−1 yr−1 (Kammann et al., 2008). The aboveground plant biomass harvest in eCO2 plots was in most years significantly higher (8–16%) than in aCO2 plots, this CO2 fertilization effect was most prominent in years with average temperatures and soil moisture and disappeared in years with extreme climatic conditions (Obermeier et al., 2017).
In July 2012, three replicate soil core samples were taken inside each of the six plots at a depth of 0–7.5 cm at random locations east, south and west of the center. July was chosen for sampling because usually the most pronounced differences in soil functioning between ambient and elevated plots were found during the summer period. The samples (18 in total) were homogenized and divided into two equal portions. One portion of each sample was stored at −20°C until further molecular analyses in the laboratory and the other portion was stored at 4°C until soil analytics.
Measurement of Soil Parameters
Soil parameters were regularly recorded since the start of the GiFACE facility in 1997. N2O flux, soil moisture content, and precipitation at the field site was measured as described by Kammann et al. (2008) and Regan et al. (2011). In brief, N2O fluxes were determined in triplicate in each GiFACE plot using chambers of 30 cm height (frustum shape; 0.184 m3 volume). Chambers were sealed for 60–90 min to permanently installed soil frames, and sampled four times in 20–30-min intervals (longer in winter-time where fluxes were lower) with 60-ml PE syringes, and N2O fluxes were calculated by linear regression. Samples were analyzed within 24 h after collection on a gas chromatograph (HP6890) equipped with an ECD. Part of the dataset (1997–2006, 2008) used in this study was published previously (Kammann et al., 2008; Regan et al., 2011), data for 2007 and from 2009 to 2013 were additionally included. Flux data as well as soil moisture content and precipitation for all plots were then outlined for these dates.
Soil pH, water content, nitrate ()-, nitrite ()-, and ammonia ()-concentrations as well as total carbon (C) and nitrogen (N) content were determined from each soil core. Soil pH was determined after shaking a soil sample (10 g) in 25 mL CaCl2 solution (0.01 M CaCl2 × 2H2O; Merck, Germany) for 20 min followed by settling for 1 h in the dark at room temperature (Schinner et al., 1996). Suspension pH was measured with an InLab® semi-micro electrode (Mettler-Toledo GmbH, Giessen, Germany). Soil moisture content (%) was determined gravimetrically by drying 1 g of homogenized soil for 3 days at 65°C in a drying oven (Memmert GmbH & Co. KG, Schwabach, Germany). The moisture content was calculated from the sample weight before and after drying. Afterwards, dried samples were ground after addition of liquid nitrogen and aliquots were analyzed at the Chemical Department of the Phillips-University Marburg (Germany) with a CHN-elemental analyzer to determine the total C/H/N percentage concentration of the soil.
To measure , , and concentrations, 1 g of soil sample was suspended in 1 mL of Nuclease-free H2O and subsequently sterile-filtered with a disposable Filter Unit (0.2 μm; Whatman, MAGV, Germany). Concentrations of and were analyzed by ion chromatography (IC; Skyam GmbH Eresing, Germany; 70°C oven temperature) equipped with a LCA A14 column (Skyam GmbH, Eresing, Germany) using a 50 μL injection volume and Na2CO3 as eluent (flow of 1.5 mL min−1). The concentration of in the soil samples was measured fluorometrically in triplicates by microscale analysis (Murase et al., 2006).
Nucleic Acid Extraction
DNA was extracted from 0.35 g soil using the NucleoSpin® Soil Kit (Machery-Nagel GmbH & Co. KG, Düren, Germany) following the manufacturer's protocol. Afterwards, the amount and purity of extracted DNA was determined with a NanDrop1000 Spectrophotometer (Thermo Scientific, Langenselbold, Germany). The concentration of DNA ranged from 100 to 120 ng μL−1. The ratios of A260/280 and A260/230 were 1.6–1.9 and 1.8–2.1, respectively, which indicated a high purity of the extracted DNA with minimum contamination of e.g., proteins and phenol (A260/280) and carbohydrates (A260/230).
Quantification of Functional Marker and 16S rRNA Genes
Copy numbers of genes encoding the denitrification associated enzymes nitrite reductase (nirK/nirS) and nitrous oxide reductase (nosZ), nitrogen fixation associated dinitrogenase (nifH), nitrification associated archaeal and bacterial ammonia monooxygenase (amoA), nitrite reductase associated with the dissimilatory reduction of nitrate to ammonia (nrfA) as well as archaeal and bacterial 16S rRNA were quantified by qPCR. Primers and PCR conditions used are given in Table S1. A typical reaction mixture contained 12.5 μL of SybrGreen Jump-Start ReadyMix (Sigma-Aldrich, Taufkirchen, Germany), 0.5 μM of each primer, 3–4.0 mM MgCl2, 2 μL of soil DNA except for amplification of nosZ, for which 3 μL of DNA were used. For the amplification of functional marker genes involved in nitrogen cycling 200 ng BSA mL−1 were added. All assays were performed in an iCycler (Applied Biosystems, Darmstadt, Germany). Standard curves were obtained using serial 10-fold dilutions of a known amount of plasmid DNA (108 to 101 gene copies) containing the respective gene fragment. Negative controls were always run with water instead of template DNA. PCR reactions were done with 1:50 and 1:100 diluted DNA extracts. Efficiencies for all assays were between 80 and 97% with r2-values between 0.971 and 0.996.
Analysis of the Composition of Functional Marker and 16S rRNA Genes
The composition of microbial communities containing nirK/nirS, nosZ, archaeal and bacterial amoA, nifH, nrfA, and archaeal 16S rRNA genes was explored by terminal restriction length polymorphism (T-RFLP) and barcode labeled 454 pyrosequencing analyses of PCR amplified gene fragments. Overall bacterial community composition based on bacterial 16S rRNA genes was assessed by T-RFLP, pyrosequencing data are available from de Menezes et al. (2016). Primers and PCR conditions used are given in Table S2. The quantity and quality of PCR amplicons were analyzed by gel electrophoresis (1.5% w/v agarose) and staining gels with 3 × GelRed Nucleic Acid Stain (Biotium, Köln, Deutschland). PCR products of the expected size were excised from the gel and purified using the DNA Wizard® SV Gel-and-PCR-Clean-up system (Promega, Mannheim, Germany).
For T-RFLP, forward or reverse primers were 5′-6-carboxyfluorescein labeled and amplicons were hydrolyzed by the restriction enzymes (FastDigest, Fermentas, St. Leon-Rot, Germany) HaeIII (nirK/nirS), HhaI (nosZ, nifH, nrfA, and amoA) and MspI and Taq1 (archaeal and bacterial 16SrRNA, respectively). Afterwards, reaction products were purified using the SigmaSpin™ Sequencing Reaction Clean-up Columns (Sigma-Aldrich) according to the manufacturer's instructions. Fluorescently labeled restriction fragments were separated on an ABI PRISM 3100 Genetic Analyzer (Applera Deutschland GmbH, Darmstadt, Germany). The lengths of fluorescently labeled terminal restriction fragments (T-RFs) were determined by comparison with an internal DNA fragment length standard (X-Rhodamine MapMarker® 30–1,000 bp; BioVentures, Murfreesboro, TN) using GeneMapper software (Applied Biosystems). Peaks with >1% of the total fluorescence of a sample and >30 bp length were analyzed by aligning fragments to the internal standard. Reproducibility of patterns was confirmed for repeated terminal restriction fragment length polymorphism (T-RFLP) analysis using the same DNA extracts of selected samples. A difference of <2 base pairs in estimated length between different profiles was the basis for considering fragments identical in size. Peak heights from different samples were normalized to identical total fluorescence units by an iterative normalization procedure (Dunbar et al., 2001).
For pyrosequencing, DNA extracts from the three replicate samples of each plot were pooled and PCR amplified using the primers used for T-RFLP but with barcode labels (6 bp) which were designed to differentiate between GiFACE plots (E1: ACACAC; E2: ATGTAT; E3: AGCAGC; A1: ATCATC; A2: AGACTA; A3: AGTCAT) and with annealing temperatures increased by 2°C due to barcode tagging (Table S2). DNA concentration was 90–180 ng as determined by a Qubit® 2.0 Fluorometer using the Quant-iT TM dsDNA BR Assay Kit (Invitrogen Darmstadt, Germany). Libraries were built by pooling amplicons (200 ng each) of each gene from soil of six GiFACE plots and subjected to barcode labeled 454 pyrosequencing (GATC, Köln, Germany).
Sequence processing and analysis was done in Qiime 1.3 (qiime.org). Pyrosequencing and PCR errors of the reads were corrected using the AmpliconNoise pipeline (Quince et al., 2011). Sequences of functional marker genes (nirK/nirS, nosZ, archaeal and bacterial amoA, nifH, and nrfA) were clustered as described previously (Caporaso et al., 2010; Palmer and Horn, 2012) using threshold similarities of 92%, because this reflects the threshold value beyond which the number of OTUs stays stable (Palmer and Horn, 2012; Palmer et al., 2012). Archaeal 16S rRNA gene sequences were clustered at 97% threshold similarities. Representative sequences were determined for each OTU. For statistical comparison of gene diversity in the plots, alpha-diversity measures were calculated in Qiime from rarefied OTU tables as described elsewhere (Hughes and Hellmann, 2005; Palmer and Horn, 2012). Rarefied OTU tables were generated by randomly subsampling original OTU tables 100 times. A sampling depth of 400 sequences was chosen for archaeal and bacterial amoA, nifH, nirK, nirS, nosZ, and nrfA to allow comparison of diversity between the different functional marker genes, as the number of sequences obtained exceeded 400 for all genes and soils. Rarefied OTU tables of archaeal 16S rRNA gene sequences were generated at a sampling depth of 150 sequences. The low number of sequences, however, clustered in up to 23 OTUs and was therefore considered sufficient to represent the most abundant taxa. The raw sequences were deposited at NCBI under the biosample accession numbers SAMN07212133 (archaeal 16S rRNA), SAMN07212134 (archaeal amoA), SAMN07212135 (bacterial amoA), SAMN07212136 (nifH), SAMN07212137 (nirK), SAMN07212138 (nirS), SAMN07212139 (nosZ), and SAMN07212140 (nrfA). The sequences representing the OTUs for each functional marker gene were deposited at NCBI under the biosample accession numbers SAMN07276914 (archaeal 16S rRNA), SAMN07276915 (archaeal amoA), SAMN07276916 (bacterial amoA), SAMN07276917 (nifH), SAMN07276918 (nirK), SAMN07276919 (nirS), SAMN07276920 (nosZ), and SAMN07276921 (nrfA).
Statistical Analyses of Collected Data
All statistical analyses were done using the statistical software R (version 3.0.1, R Development Core Team, 2013). Significant differences in copy numbers of archaeal/bacterial 16S rRNA genes, archaeal and bacterial amoA, nirK, nirS, nosZ, nifH, and nrfA were assessed using ANOVA (P < 0.05). All quantitative data were log-transformed prior to analysis to satisfy the assumptions of homoscedasticity and normally distributed residuals.
The effect of soil parameters on T-RFLP community profiles was explored by canonical correspondence analysis (CCA). Statistical significance of the CCA was assessed using permutation test (1,000 iterations). All community composition data were log-transformed before analysis, in order to reach normal distribution.
N2O Flux and Soil Moisture Content over a Period of 12 Years
N2O flux data collected since the start of the GiFACE facility in 1998 until June 2013 revealed an average flux of 21 μg N2O-N m−2 h−1 (1.84 kg N2O-N ha−1 yr−1) for the soil fumigated with eCO2 (E) and of 10 μg N2O-N m−2 h−1 (0.88 kg N2O-N ha−1 yr−1) for soil under aCO2 (A). The largest differences between fluxes from soil fumigated with eCO2 and aCO2 occurred during the first 3 years after the start of the experiment and the highest frequency of events occurred in years 2005–2006. Only one event occurred where the N2O flux from soil under aCO2 was higher (P < 0.05; Figure 1). Mostly but not generally, a rain event, which resulted in increased soil moisture content, preceded higher N2O fluxes (Figures S1, S2). For example, on day 866 (September 7, 1999) and on day 4055 (May 31, 2008) N2O fluxes from soil under eCO2 exceeded fluxes from soil under aCO2 by 73 and 104 μg N2O-N m−2 h−1, respectively. Both flux events were preceded by precipitation of 25.6 mm (day 865) and of 56.8 mm (day 4,054) which accounted for more than half of the cumulative precipitation in September 1999 (30.8 mm) and May 2008 (72.8 mm; Figure S1) which in consequence increased soil moisture concentration levels by more than 10% (Figure S2).
Figure 1. N2O flux measurements at GiFACE from 1997 to 2013 shown as the differences in fluxes between eCO2 and aCO2 plots. (A) Difference of mean N2O fluxes from soil at elevated (E) and ambient (A) CO2. Triangles mark occasions where N2O fluxes from eCO2 plots were significantly larger than from aCO2 plots (black triangles) or vice versa (white triangles) tested by ANOVA (P < 0.05). (B–D) Differences in N2O fluxes in the three sets (B, E1/A1; C, E2/A2; D, E3/A3).
Soil characteristics differed only marginally and non-significantly between plots fumigated with eCO2 and aCO2 (Table 1). Significant differences occurred only between GiFACE sets (Table S3) but not between eCO2 and aCO2 plots. The soil was moderately acidic with pH ranging from 5.45 to 6.10. Differences existed mostly between soil of the first set and the two other sets. In set E1/A1 pH (5.55) was lower than in set E2/A2 (6.03) and concentration and C-content were lower than in E2/A2 and E3/A3, respectively. N-content was lower at E1/A1, while -concentration and C:N ratio was higher than in the two other sets. Water-content of the soil samples was similar in all plots.
Abundance of Microbial Groups Associated with Soil Nitrogen Cycling
Total bacterial 16S rRNA gene copy numbers were in the order of 1 × 109 g−1 dw soil, while archaeal 16S rRNA gene copy numbers ranged between 5 × 107 and 1 × 108 g−1 dw soil for all plots (Figure 2). The abundance of bacteria and archaea did not differ significantly between plots or sets (Figure 2). The absolute copy numbers of the functional marker genes nirK, nosZ, nrfA, nifH, and archaeal amoA (Figure 2) and their numbers relative to total (bacterial + archaeal) 16S rRNA gene copies were similarly high in all plots (Table 2, Figure S2). The genes nirS and bacterial amoA were 5- to 10-fold and 100-fold less abundant than the other functional marker genes, respectively (Figure 2). Absolute and relative copy numbers of nirS in plot A2 were similar to numbers in plot E2 but significantly higher than in each individual plot of the sets E1/A1 and E3/A3 (Figure 2, Table S4). Relative numbers were also significantly higher in set E2/A2 than in the other two sets (Table 2). Absolute and relative numbers of bacterial amoA in plot A3 were significantly higher than in plots A1 and E1 (Figure 2, Table S4) and relative numbers of bacterial amoA in set E3/A3 were significantly higher than in set E1/A1 (Table 2). Comparison of archaeal amoA and 16S rRNA gene copy numbers indicated that a large fraction of archaea harbored a copy of the amoA gene (ratios close to one, data not shown). For denitrification genes, copy numbers of nirK exceeded nirS which is also reflected in higher ratios of nosZ/nirS of 2–13 compared to 0.36–0.58 for nosZ/nirK which was low in all plot (Figure S3).
Figure 2. Abundance of dissimilatory nitrate reducers, denitrifiers, nitrogen fixers, ammonia oxidizers and total bacteria and archaea based on quantitative PCR analysis of the functional marker genes (nrfA, nosZ, nirS, nirK, nifH, archaeal and bacterial amoA) as well as of 16S rRNA genes, respectively. Bars indicate the total gene copy numbers. (Mean ± SD, n = 3). Different letters indicate significant differences in the abundance of a functional group between plots.
Table 2. Abundance of functional marker genes (archaeal and bacterial amoA, nirK, nirS, nosZ, nrfA, and nifH) relative to total 16S rRNA gene abundance (archaeal + bacterial) in soil of GiFACE sets E1/A1, E2/A3 and E3/A3.
Composition of Microbial Communities Involved in Soil Nitrogen Cycling
Applying a threshold similarity of 92% (97% for archaeal 16S rRNA genes) to sequences obtained from pyrosequencing, coverage of the libraries was in the range of 77.2–100% (Table 3). For amoA (archaeal and bacterial) the number of operational taxonomic units (OTUs) was low with only 3–7 OTUs observed and 3-8 OTUs estimated. Hence diversity calculated from rarefied tables (Shannon Diversity index Harch. amoA = 0.90–1.59; Hbact. amoA = 0.64–1.37) was also low. Evenness of the archaeal ammonia oxidizer community and of the overall archaeal communities ranged from 0.47 to 0.64 and from 0.65–0.73, respectively, as the communities were dominated by only two amoA (OTU 1, 51.5–71.1%; OTU 2, 26.6–35.4%) as well as by one 16S rRNA (OTU 1, 34.7–42.0%) genotypes in all GiFACE plots (Table S5). These genotypes were closely related to amoA and the 16S rRNA gene from Candidatus Nitrosphaera gargensis and Cand. Nitrosphaera viennensis, respectively. Evenness (E = 0.25–0.50) was even lower for bacterial ammonia oxidizers which were dominated by two amoA genotypes (OTU 1, 75.4–84.2%; OTU 2, 14.2–22.5%) closely related to amoA of Nitrospira spp.
Table 3. Analysis of representative, Qiime-clustered sequences of PCR amplified gene fragments from GiFACE soil.
Numbers of observed OTUs for marker genes for N-fixation, denitrification and DNRA were at least one order of magnitude higher than for archaeal and bacterial amoA genes (Table 3). Communities of N-fixers, denitrifiers, and dissimilatory nitrate reducers were also more diverse than ammonia oxidizers (HnifH = 4.53–5.25; HnirK = 4.11; HnirS = 2.80–3.75; HnosZ = 2.90–4.64; HnrfA = 1.72–6.53) and also more even (EnifH = 0.76–0.82; EnirK = 0.62–0.71; EnirS = 0.51–0.65; EnosZ = 0.53–0.76; EnrfA = 0.41–0.87). The lowest evenness levels of E = 0.41 and 0.47 were found for DNRA communities of plot E3 and A3.
Generally, OTUs were most closely related to genes originating from as yet uncultured species but sequence identities of >71% to genes from cultivated species known to be involved in N-cycling confirmed that these genes were indeed derived from organisms of the respective target group. OTUs representing species of Bradyrhizobium were most abundant among N-fixers and nitrite reducers, while an OTU representing Rhodopseudomonas palustris dominated the nosZ-containing denitrifier communities. Communities of organisms capable of DNRA in sets E1/A1 and E2/A2 were not dominated by single OTUs and sequences were most closely related to nrfA from Bacteroides spp., Anaromyxobacter spp., Sorangium spp., and Geobacter spp.
Influence of Soil Characteristics on the Composition of Microbial Communities Involved in Soil Nitrogen Cycling
CCA based on T-RFLP data clustered N-fixer communities (nifH) into three distinct groups which corresponded to set E1/A1, E2/A2, and E3/A3 (Figure 3). Communities of denitrifiers (nirK/nirS, nosZ) and archaeal ammonia oxidizers (archaeal amoA) as well as the overall archaeal community (16S rRNA genes) in set E1/A1 were distinct from those in the other two sets. Lower pH and soil nitrate concentration as well as higher ammonia concentration separated communities of set E1/A1 from those of the other two sets (Table 1, Figure 3). The DNRA (nrfA), bacterial ammonia oxidizer (bacterial amoA) and overall bacterial communities (16S rRNA) showed no clustering according to GiFACE sets (Figure 3).
Figure 3. Canonical correspondence analysis (CCA) biplots based on T-RFLP community analyses of nirK (A), nirS (B), nosZ (C), nifH (D), archaeal amoA (E), bacterial amoA (F), nrfA (G), bacterial 16S rRNA genes (H), and archaeal 16S rRNA genes (I). Arrows indicate the direction and relative importance (arrow lengths) of soil parameters associated with the clustering of the communities. For each gene the most important environmental variables are displayed and highlighted by an asterisk if significant in the model (ANOVA: P < 0.05). Square, triangle, and circle symbols represent sets E1/A1, E2/A2, and E3/A3, respectively. Closed symbols represent fumigation with eCO2 and open symbols the control plot at aCO2. (n = 3).
CCA identified pH (16.4–29.6% of the variance; Table 4) and ammonia concentration (12.9–30.7% of the variance; Table 4) as the most important soil parameters to shape the soil microbial communities. Both exerted a significant impact on microbial community composition independent of the gene considered. Nitrate concentration (5.2–20.3% of the variance; Table 4) also determined the composition of the microbial communities except for the communities of nosZ-containing denitrifiers, dissimilatory nitrate reducers, and archaeal ammonia oxidizers. The level of CO2, whether elevated or ambient, influenced the composition of dissimilatory nitrate reducer communities but the other communities were unaffected (Table 4).
Table 4. Proportion of variance in soil microbial communities in soil of GiFACE plots explained by environmental variables (percentage of total variation).
Exploring whether CO2 exerted an influence on community composition in single sets showed that different microbial communities were affected. The level of CO2 determined the composition of the archaeal community in set E1/A1, the bacterial community and nirS-type denitrifiers in set E2/A2, and of dissimilatory nitrate reducers in set E3/A3 (Table 5).
Table 5. Influence of elevated atmospheric CO2 on the composition of microbial communities associated with nitrogen cycling in soil of GiFACE sets E1/A1, E2/A2, and E3/A3.
Given the functional differences in N-cycling observed by FACE experiments in Giessen and around the world (Kammann et al., 2008; Müller et al., 2009; Rütting et al., 2010; van Groenigen et al., 2011), we hypothesized adaptation of the soil microbial communities, i.e., differences in abundance and composition, to long-term exposure to elevated CO2. For the study period (1998–2013) we found 2-fold higher average N2O fluxes at eCO2. Thus, compared to the initial experimental period lasting from 1998 to 2006, the difference between fluxes at eCO2 and aCO2 (0.90 under aCO2 vs. 2.07 kg N2O–N ha−1 y−1 under eCO2) remained stable (Kammann et al., 2008). The continuous difference in functionality, altered N transformation rates as well as lower nitrate and higher ammonia concentrations observed at the GiFACE under elevated CO2 (Müller et al., 2009) suggest differences in the underlying microbial communities and particularly of communities associated with N-cycling. Our study, however, shows that soil microbial communities were surprisingly unaffected by elevated levels of CO2. In contrast to our hypothesis, the abundance and composition of the soil microbial communities associated with N-cycling in a given GiFACE set (E1 vs. A1, A2 vs. E2 and E3 vs. A3) were largely unaffected by CO2 level. A lack of response of microbial communities against long-term exposure to eCO2 was reported previously (Haase et al., 2008; Nelson et al., 2010; Marhan et al., 2011; Regan et al., 2011; Pujol Pereira et al., 2013; Dunbar et al., 2014). In addition, de Menezes et al. (2016) showed for GiFACE that eCO2 exerted no substantial effects on the composition of the overall soil bacterial community. Likewise, microbial communities involved in denitrification, ammonia oxidation, and DNRA remained unaffected by eCO2 in other studies (Deiglmayr et al., 2004; Haase et al., 2008; Marhan et al., 2011; Pujol Pereira et al., 2013).
These results, however, disagree with other studies that showed effects of eCO2 on soil microbial communities in general (Drigo et al., 2008, 2009; He et al., 2010; Deng et al., 2012; Xu et al., 2013; Dunbar et al., 2014; Xiong et al., 2015; Xia et al., 2017) or on communities associated with N-cycling (Horz et al., 2004; Fromin et al., 2005; Roussel-Delif et al., 2005; Lesaulnier et al., 2008; He et al., 2010; Regan et al., 2011). Yet, a comparison of results from different FACE facilities is ambiguous as the CO2 concentration applied (up to an overall concentration of +50%) varied. Denef et al. (2007) for instance, demonstrated even for GiFACE differences in PFLA patterns between soil fumigated with eCO2 and aCO2 albeit at concentrations of +50%. Changes in community composition and abundance were also observed in rice root samples, but again in response to higher elevation of CO2 (Okubo et al., 2015). Higher eCO2 levels may have resulted in a higher C-input into the soil by plants in these studies even though the total amount of nitrogen is almost equal between the study sites (Feng et al., 2015). Depending on the plant population the amount of C provided by the plants differs between 20 and 50% of total CO2 uptake (Kuzyakov and Domanski, 2000) and only a small fraction is available for microbial biomass production (van Veen et al., 1991). Thus, it remains unclear whether an elevation of CO2 by +20% suffices to increase C-inputs into the rhizosphere and raises the question whether a threshold in CO2 concentration may exist where not only the activity of the soil microbial communities is affected but also a community response will be detectable.
Haase et al. (2008) attributed the lack of response in microbial community abundance to unaltered C-flux from the whole root system of Phaseolus vulgaris into the soil. At GiFACE, the additional C assimilation also did not result in increased soil C sequestration. Instead, a loss of soil C, together with the breakup of large macroaggregates, was detected and caused enhanced ecosystem respiration under eCO2 (Lenhart, 2008). Influence of higher labile C input by the plant-root system may occur only directly at the root-soil interface and would then be rapidly consumed by microorganisms attached or located around the roots (Haase et al., 2008). It was also reported that fungal biomass was more strongly influenced by elevated CO2 than bacterial biomass (Jones et al., 1998; Drigo et al., 2009), but other studies found a negligible effect on fungal communities by eCO2 (Guenet et al., 2012; Dunbar et al., 2014; Lee et al., 2015). In our study, the level of CO2 had a general effect on the composition of dissimilatory nitrate reducer communities and affected the composition of additional but distinct communities in one but not all GiFACE sets. The effect on dissimilatory nitrate reducers agrees well with the findings by Müller et al. (2009) that DNRA rates increased by ~150% while soil nitrate content decreased under eCO2.
The persisting differences in N2O fluxes (Kammann et al., 2008) suggested differences in soil parameters at elevated and ambient CO2 levels but neither lower soil nitrate levels under eCO2 occurred nor was any other soil parameter affected by eCO2. Presumably, the impact of cultivation as permanent grassland for over 100 years had a more profound effect on the soil characteristics than 14 years of moderate exposure to eCO2. In addition, the increase in atmospheric CO2 concentration from ~300 to 400 ppm in the last 100 years is larger than the experimental exposure at GiFACE. Hence, we assume that long-term cultivation and the increase in CO2 prior to the experimental period led to the development of microbial communities which are adapted to the prevailing soil conditions but seem unresponsive to moderately increased CO2 levels. The grassland soil has been under long-term stable management as a hay meadow and has not been plowed for at least 100 years (Kammann et al., 2008). Plowing, fertilization regime, and cropping with annual vs. perennial plants were shown to influence N2O fluxes and ammonia oxidizer and denitrifier soil bacterial communities (Thompson et al., 2016). Thompson et al. found decreased abundance and diversity of denitrifiers after plowing and ammonia oxidizer (archaeal and bacterial) communities which differed between soil cropped with annual and perennial plants.
In all plots a large fraction of sequences belonged to only a few OTUs which were most closely related to genes from a limited range of species. These groups may thus represent the well-adapted key players of N-cycling in the soil and occurred in almost identical relative abundance under eCO2 and aCO2. N-fixers and denitrifiers (nirK, nirS, and nosZ) were primarily represented by OTUs similar to functional genes of Bradyrhizobia, a group with known N-fixation and denitrification capability (Bedmar et al., 2005). Ammonia oxidizers were dominated by Nitrosopira spp., the prevalent bacterial ammonia oxidizer species in soils (Kowalchuk et al., 2000) and Nitrososphaera spp., a typical archaeal species in soils related to Nitrosphaera viennensis which was isolated from soil (Tourna et al., 2011; Stieglmeier et al., 2014). Nitrososphaera spp. was representative for both the archaeal ammonia oxidizer as well as the total archaeal community. Although bacterial and archaeal ammonia oxidizers clustered only into few OTUs, our results agree with previous findings on the relevance of these species in soils. OTUs representing Bacteriodetes and Anaeromyxobacter dominated the microbial community involved in DNRA. These species are known to also harbor a clade II nosZ gene (Sanford et al., 2012; Jones et al., 2013) which we did not study here but which may contribute to N2O metabolism at GiFACE.
Communities varied between different sets at GiFACE and were related to differences in soil parameters determined by the localization of the sets in the experimental field. A previous study found differences in the water level in the deeper soil layers at GiFACE in the order: E1/A1 < E3/A3 < E2/A2 (Lenhart, 2008). de Menezes et al. (2016) attributed differences in bacterial community composition at GiFACE to the soil moisture gradient. Long-term differences in the water level of the deeper soil layers probably determined the differences in pH and the concentration of C- and N-compounds observed in the top soil layer (0–7.5 cm) of the sets. Microbial community composition at GiFACE was influenced by pH, nitrate, and ammonia. Soil parameters are known as the predominant drivers determining the distribution of microorganisms and shaping their communities (e.g., Zhou et al., 2008). Variation in e.g., soil denitrifier (Enwall et al., 2010) and ammonia oxidizer communities (Wessén et al., 2010) was previously found to occur with spatial heterogeneity at scales similar to those at the GiFACE experimental site. Regan et al. (2011) also found stronger influence of the location of the GiFACE sets or the soil depth on the abundance of amoA, nirK, nirS, and nosZ than of eCO2. Likewise, Marhan et al. (2011) observed a similar trend and that temporal variation and soil depth had a greater effect on the abundance of nitrate reducers and bacteria than eCO2. Microbial communities in a soybean agroecosystem were also not affected by eCO2 but by lower soil moisture compared to ambient conditions (Pujol Pereira et al., 2013). Experimentally reduced precipitation increased the mass of microaggregates concomitant with a higher abundance of bacteria and nosZ-containing denitrifiers in microhabitats protected from reduced soil moisture (Pujol Pereira et al., 2013).
The long-term cultivation as grassland as well as the increase in atmospheric CO2 in the past 100 years seems to have shaped soil microbial communities which density and composition remained unaltered by moderately eCO2. Despite a lack of response of microbial communities associated with N-cycling, N-transformations were altered and average N2O fluxes doubled under eCO2. Large fluxes were related to a high N-status of the soil after fertilization and plant growth in spring. This correlation suggests differential response to altered environmental conditions against the background of a community which remained unaltered. Palmer et al. (2016), for instance, observed differential transcriptional activation of denitrifiers in response to drying-rewetting and flooding. Hence, future studies should explore in more detail how elevated CO2 in conjunction with massive N inputs during fertilization impact microbial communities in the soil and whether this leads to a short-term activation of microbial groups involved in N-cycling and hence higher production of N2O.
KB designed the study, performed the field sampling, performed all lab work (nucleic-acid extractions, T-RFLP analysis, qPCR analysis, preparation of 454 pyrosequencing), performed statistical analysis, evaluated the data and wrote the manuscript. KK performed the analysis of 454 pyrosequencing data and revised the manuscript. MH performed the analysis of 454 pyrosequencing data and revised the manuscript. GM designed the study, provided climate data and revised the manuscript. CK provided climate and N2O data and revised the manuscript. CG designed the gas analyses study and performed sampling, provided N2O data and revised the manuscript. CM developed the research idea, designed the study and revised the manuscript. GB developed the research idea, designed the research, and wrote the manuscript.
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.
This work is part of the dissertation of KB which can be accessed online (Brenzinger, 2015). Its publication is in line with the author's university policy. The project has been funded by a grant of the German Science Foundation (DFG)-Project MU1302/6-1 to CM and GB, by the Max-Planck Society, by FACE2FACE within the LOEWE research fund for the development of scientific and economic excellence of the federal state Hessian, and by the Hessian Agency for Nature Conservation, Environment and Geology (HLNUG).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2017.01976/full#supplementary-material
Andresen, L. C., Yuan, N., Seibert, R., Moser, G., Kammann, C., Luterbacher, J., et al. (2017). Biomass responses in a temperate European grassland through 17 years of elevated CO2. Glob. Change Biol. 2017, 1–11. doi: 10.1111/gcb.13705
Barnard, R., Le Roux, X., Hungate, B. A., Cleland, E. E., Blankinship, J. C., Barthes, L., et al. (2006). Several components of global change alter nitrifying and denitrifying activities in an annual grassland. Funct. Ecol. 20, 557–564. doi: 10.1111/j.1365-2435.2006.01146.x
Barnard, R., Leadley, P. W., Lensi, R., and Barthes, L. (2005). Plant, soil microbial and soil inorganic nitrogen responses to elevated CO2: a study in microcosms of Holcus lanatus. Acta Oecol. 27, 171–178. doi: 10.1016/j.actao.2004.11.005
Bedmar, E. J., Robles, E. F., and Delgado, M. J. (2005). The complete denitrification pathway of the symbiotic, nitrogen-fixing bacterium Bradyrhizobium japonicum. Biochem. Soc. Trans. 33, 141–144. doi: 10.1042/BST0330141
Butterbach-Bahl, K., Baggs, E. M., Dannenmann, M., Kiese, R., and Zechmeister-Boltenstern, S. (2013). Nitrous oxide emissions from soils: how well do we understand the processes and their controls? Philos. Trans. R. Soc. Lond. B. Biol. Sci. 368:20130122. doi: 10.1098/rstb.2013.0122
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
de Menezes, A. B., Müller, C., Clipson, N., and Doyle, E. (2016). The soil microbiome at the Gi-FACE experiment responds to a moisture gradient but not to CO2 enrichment. Microbiology 162, 1572–1582. doi: 10.1099/mic.0.000341
Deiglmayr, K., Philippot, L., Hartwig, U. A., and Kandeler, E. (2004). Structure and activity of the nitrate-reducing community in the rhizosphere of Lolium perenne and Trifolium repens under long-term elevated atmospheric pCO2. FEMS Microbiol. Ecol. 49, 445–454. doi: 10.1016/j.femsec.2004.04.017
Denef, K., Bubenheim, H., Lenhart, K., Vermeulen, J., van Cleemput, O., Boeckx, P., et al. (2007). Community shifts and carbon translocation within metabolically-active rhizosphere microorganisms in grasslands under elevated CO2. Biogeosciences 4, 769–779. doi: 10.5194/bg-4-769-2007
Deng, Y., He, Z., Xu, M., Qin, Y., van Nostrand, J. D., Wu, L., et al. (2012). Elevated carbon dioxide alters the structure of soil microbial communities. Appl. Environ. Microbiol. 78, 2991–2995. doi: 10.1128/AEM.06924-11
Drigo, B., Kowalchuk, G. A., and van Veen, J. A. (2008). Climate change goes underground: effects of elevated atmospheric CO2 on microbial community structure and activities in the rhizosphere. Biol. Fertil. Soils 44, 667–679. doi: 10.1007/s00374-008-0277-3
Drigo, B., van Veen, J. A., and Kowalchuk, G. A. (2009). Specific rhizosphere bacterial and fungal groups respond differently to elevated atmospheric CO2. ISME J. 3, 1204–1217. doi: 10.1038/ismej.2009.65
Dunbar, J., Gallegos-Graves, L. V., Steven, B., Mueller, R., Hesse, C., Zak, D. R., et al. (2014). Surface soil fungal and bacterial communities in aspen stands are resilient to eleven years of elevated CO2 and O3. Soil Biol. Biochem. 76, 227–234. doi: 10.1016/j.soilbio.2014.05.027
Dunbar, J., Ticknor, L. O., and Kuske, C. R. (2001). Phylogenetic specificity and reproducibility and new method for analysis of terminal restriction fragment profiles of 16S rRNA genes from bacterial communities. Appl. Environ. Microbiol. 67, 190–197. doi: 10.1128/AEM.67.1.190-197.2001
Enwall, K., Throbäck, I. N., Stenberg, M., Söderström, M., and Hallin, S. (2010). Soil resources influence spatial patterns of denitrifying communities at scales compatible with land management. Appl. Environ. Microbiol. 76, 2243–2250. doi: 10.1128/AEM.02197-09
Feng, Z., Rütting, T., Pleijel, H., Wallin, G., Reich, P. B., Kammann, C. I., et al. (2015). Constraints to nitrogen acquisition of terrestrial plants under elevated CO2. Glob. Change Biol. 21, 3152–3168. doi: 10.1111/gcb.12938
Fromin, N., Tarnawski, S., Roussel-Delif, L., Hemelin, J., Baggs, E. M., and Aragno, M. (2005). Nitrogen fertiliser rate affects the frequency of nitrate-dissimilating Pseudomonas spp. in the rhizosphere of Lolium perenne grown under elevated CO2 (Swiss FACE). Soil Biol. Biochem. 37, 1962–1965. doi: 10.1016/j.soilbio.2005.02.003
Guenet, B., Lenhart, K., Leloup, J., Giusti-Miller, S., Pouteau, V., et al. (2012). The impact of long-term CO2 enrichment and moisture levels on soil microbial community structure and enzyme activities. Geoderma 170, 331–336. doi: 10.1016/j.geoderma.2011.12.002
Haase, S., Philippot, L., Neumann, G., Marhan, S., and Kandeler, E. (2008). Local response of bacterial densities and enzyme activities to elevated atmospheric CO2 and different N supply in the rhizosphere of Phaseolus vulgaris L. Soil Biol. Biochem. 40, 1225–1234. doi: 10.1016/j.soilbio.2007.12.025
He, Z., Xu, M., Deng, Y., Kang, S., Kellogg, L., Wu, L., et al. (2010). Metagenomic analysis reveals a marked divergence in the structure of belowground microbial communities at elevated CO2. Ecol. Lett. 13, 564–575. doi: 10.1111/j.1461-0248.2010.01453.x
Horz, H. P., Barbrook, A., Field, C. B., and Bohannan, B. J. M. (2004). Ammonia-oxidizing bacteria respond to multifactorial global change. Proc. Natl. Acad. Sci. U.S.A. 101, 15136–15141. doi: 10.1073/pnas.0406616101
IPCC (2014) “Climate change 2014: synthesis report,” in Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change eds R. K. Pachauri and L. A. Meyer (Geneva: IPCC) 151.
Jäger, H. J., Schmidt, S. W., Kammann, C., Grünhage, L., Müller, C., and Hanewald, K. (2003). The University of Giessen free-air carbon dioxide enrichment study: description of the experimental site and of a new enrichment system. J. Appl. Bot. 77, 117–127.
Jones, C. M., Graf, D. R. H., Bru, D., Philippot, L., and Hallin, S. (2013). The unaccounted yet abundant nitrous oxide-reducing microbial community: a potential nitrous oxide sink. ISME J. 7, 417–426. doi: 10.1038/ismej.2012.125
Jones, T. H., Thompson, L. J., Lawton, J. H., Bezemer, T. M., Bardgett, R. D., Blackburn, T. M., et al. (1998). Impacts of rising atmospheric carbon dioxide on model terrestrial ecosystems. Science 280, 441–443. doi: 10.1126/science.280.5362.441
Kammann, C., Grünhage, L., Grüters, U., Janze, S., and Jäger, H.-J. (2005). Response of aboveground grassland biomass and soil moisture to moderate long-term CO2 enrichment. Basic Appl. Ecol. 6, 351–365. doi: 10.1016/j.baae.2005.01.011
Kowalchuk, G. A., Stienstra, A. W., Heilig, G. H. J., Stephen, J. R., and Woldendorp, J. W. (2000). Molecular analysis of ammonia-oxidising bacteria in soil of successional grasslands of the Drentsche A (The Netherlands). FEMS Microbiol. Ecol. 31, 207–214. doi: 10.1111/j.1574-6941.2000.tb00685.x
Lee, S. H., Kim, S. Y., Ding, W. X., and Kang, H. (2015). Impact of elevated CO2 and N addition on bacteria, fungi, and archaea in a marsh ecosystem with various types of plants. Appl. Microbiol. Biotechnol. 99, 5295–5305. doi: 10.1007/s00253-015-6385-8
Lenhart, K. (2008). The effects of Long-Term Free Air CO2 Enrichment (FACE) on Soil Aggregation, Soil Carbon Input, and Ecosystem CO2 Dynamics in a Temperate Grassland Ecosystem. dissertation, Justus-Liebig University, Giessen.
Lesaulnier, C., Papamichail, D., McCorkle, S., Ollivier, B., Skiena, S., Taghavi, S., et al. (2008). Elevated atmospheric CO2 affects soil microbial diversity associated with trembling aspen. Environ. Microbiol. 10, 926–941. doi: 10.1111/j.1462-2920.2007.01512.x
Marhan, S., Philippot, L., Bru, D., Rudolph, S., Franzaring, J., Hogy, P., et al. (2011). Abundance and activity of nitrate reducers in an arable soil are more affected by temporal variation and soil depth than by elevated atmospheric CO2. FEMS Microbiol. Ecol. 76, 209–219. doi: 10.1111/j.1574-6941.2011.01048.x
Müller, C., Rütting, T., Abbasi, M. K., Laughlin, R. J., Kammann, C., Clough, T. J., et al. (2009). Effect of elevated CO2 on soil N dynamics in a temperate grassland soil. Soil Biol. Biochem. 41, 1996–2001. doi: 10.1016/j.soilbio.2009.07.003
Murase, J., Noll, M., and Frenzel, P. (2006). Impact of protists on the activity and structure of the bacterial community in a rice field soil. Appl. Environ. Microbiol. 72, 5436–5444. doi: 10.1128/AEM.00207-06
Nelson, D. M., Cann, I. K. O., and Mackie, R. I. (2010). Response of archaeal communities in the rhizosphere of maize and soybean to elevated atmospheric CO2 concentrations. PLoS ONE 5:e15897. doi: 10.1371/journal.pone.0015897
Obermeier, W. A., Lehnert, L. W., Kammann, C. I., Müller, C., Grünhage, L., Luterbacher, J., et al. (2017). Reduced CO2 fertilization effect in temperate C3 grasslands under more extreme weather conditions. Nat. Clim. Change 7, 137–141. doi: 10.1038/nclimate3191
Okubo, T., Liu, D., Tsurumaru, H., Ikeda, S., Asakawa, S., Tokida, T., et al. (2015). Elevated atmospheric CO2 levels affect community structure of rice root-associated bacteria. Front. Microbiol. 6:136. doi: 10.3389/fmicb.2015.00136
Palmer, K., and Horn, M. A. (2012). Actinobacterial nitrate reducers and proteobacterial denitrifiers are abundant in N2O-metabolizing palsa peat. Appl. Environ. Microbiol. 78, 5584–5596. doi: 10.1128/AEM.00810-12
Palmer, K., Biasi, C., and Horn, M. A. (2012). Contrasting denitrifier communities relate to contrasting N2O emission patterns from acidic peat soils in arctic tundra. ISME J. 6, 1058–1077. doi: 10.1038/ismej.2011.172
Palmer, K., Köpp, J., Gebauer, G., and Horn, M. A. (2016). Drying-rewetting and flooding impact denitrifier activity rather than community structure in a moderately acidic fen. Front. Microbiol. 7:727. doi: 10.3389/fmicb.2016.00727
Philippot, L., Andert, J., Jones, C. M., Bru, D., and Hallin, S. (2011). Importance of denitrifiers lacking the genes encoding the nitrous oxide reductase for N2O emissions from soil. Glob. Change Biol. 17, 1497–1504. doi: 10.1111/j.1365-2486.2010.02334.x
Pujol Pereira, E. I. P., Chung, H., Scow, K., and Six, J. (2013). Microbial communities and soil structure are affected by reduced precipitation but not by elevated carbon dioxide. Soil Sci. Soc. Am. J. 77, 482–488. doi: 10.2136/sssaj2012.0218
R Development Core Team (2013). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing. Available online at: http://www.R-project.org/
Regan, K., Kammann, C., Hartung, K., Lenhart, K., Müller, C., Philippot, L., et al. (2011). Can differences in microbial abundances help explain enhanced N2O emissions in a permanent grassland under elevated atmospheric CO2? Glob. Change Biol. 17, 3176–3186. doi: 10.1111/j.1365-2486.2011.02470.x
Roussel-Delif, L., Tarnawski, S., Hamelin, J., Philippot, L., Aragno, M., and Fromin, N. (2005). Frequency and diversity of nitrate reductase genes among nitrate-dissimilating Pseudomonas in the rhizosphere of perennial grasses grown in field conditions. Microb. Ecol. 49, 63–72. doi: 10.1007/s00248-003-0228-3
Rütting, T., Clough, T. J., Müller, C., Lieffering, M., and Newton, P. C. (2010). Ten years of elevated atmospheric carbon dioxide alters soil nitrogen transformations in a sheep-grazed pasture. Glob. Change Biol. 16, 2530–2542. doi: 10.1111/j.1365-2486.2009.02089.x
Sanford, R. A., Wagner, D. D., Wu, Q., Chee-Sanford, J. C., Thomas, S. H., Cruz-García, C., et al. (2012). Unexpected nondenitrifier nitrous oxide reductase gene diversity and abundance in soils. Proc. Natl. Acad. Sci. U.S.A. 109, 19709–19714. doi: 10.1073/pnas.1211238109
Stieglmeier, M., Klingl, A., Alves, R. J., Simon, K. M. R., Melcher, M., Leisch, N., et al. (2014). Nitrososphaera viennensis sp. nov., an aerobic and mesophilic ammonia-oxidizing archaeon from soil and member of the archaeal phylum Thaumarchaeota. Int. J. Syst. Evol. Microbiol. 64, 2738–2752. doi: 10.1099/ijs.0.063172-0
Thompson, K. A., Bent, E., Abalos, D., Wagner-Riddle, C., and Dunfield, K. E. (2016). Soil microbial communities as potential regulators of NO fluxes in annual and perennial cropping systems. Soil Biol. Biochem. 103, 262–273. doi: 10.1016/j.soilbio.2016.08.030
Tourna, M., Stieglmeier, M., Spang, A., Könneke, M., Schintlmeister, A., Urich, T., et al. (2011). Nitrososphaera viennensis, an ammonia oxidizing archaeon from soil. Proc. Natl. Acad. Sci. U.S.A. 108, 8420–8425. doi: 10.1073/pnas.1013488108
van Veen, J. A., Liljeroth, E., Lekkerkerk, L. J. A., and Van de Geijn, S. C. (1991). Carbon fluxes in plant-soil systems at elevated atmospheric CO2 levels. Ecol. Appl. 1, 175–181. doi: 10.2307/1941810
Wessén, E., Hallin, S., and Philippot, L. (2010). Differential responses of bacterial and archaeal groups at high taxonomical ranks to soil management. Soil Biol. Biochem. 42, 1759–1765. doi: 10.1016/j.soilbio.2010.06.013
Xia, W., Jia, Z., Bowatte, S., and Newton, P. C. D. (2017). Impact of elevated atmospheric CO2 on soil bacteria community in a grazed pasture after 12-year enrichment. Geoderma 285, 19–26. doi: 10.1016/j.geoderma.2016.09.015
Xiong, J., He, Z., Shi, S., Kent, A., Deng, Y., Wu, L., et al. (2015). Elevated CO2 shifts the functional structure and metabolic potentials of soil microbial communities in a C4 agroecosystem. Sci. Rep. 5:9316. doi: 10.1038/srep09316
Xu, M., He, Z., Deng, Y., Wu, L., van Nostrand, J. D., Hobbie, S. E., et al. (2013). Elevated CO2 influences microbial carbon and nitrogen cycling. BMC Microbiol. 13:124. doi: 10.1186/1471-2180-13-124
Keywords: elevated CO2, N2O, denitrifiers, ammonia oxidizers, N-fixers, DNRA, FACE
Citation: Brenzinger K, Kujala K, Horn MA, Moser G, Guillet C, Kammann C, Müller C and Braker G (2017) Soil Conditions Rather Than Long-Term Exposure to Elevated CO2 Affect Soil Microbial Communities Associated with N-Cycling. Front. Microbiol. 8:1976. doi: 10.3389/fmicb.2017.01976
Received: 02 April 2017; Accepted: 25 September 2017;
Published: 18 October 2017.
Edited by:Jeanette M. Norton, Utah State University, United States
Reviewed by:Zhili He, University of Oklahoma, United States
Thomas Pommier, Institut National de la Recherche Agronomique (INRA), France
Copyright © 2017 Brenzinger, Kujala, Horn, Moser, Guillet, Kammann, Müller and Braker. 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: Gesche Braker, email@example.com