INTRODUCTION
High-latitude peatlands store approximately one-third of global soil carbon and may pose a climatic threat if rising global temperatures accelerate the release of this stored carbon in gaseous forms, as either carbon dioxide or methane (
1–3). Mineral-poor (ombrotrophic) peatlands receive most of their nutrient inputs from atmospheric deposition and contain
Sphagnum moss as their primary plant cover (
2,
4). The peat moss
Sphagnum is a keystone genus in these ecosystems and is responsible for much of the primary production and recalcitrant dead organic matter (
5,
6).
Sphagnum mosses also host complex microbiomes (
7–10), including N
2 fixers (diazotrophs) that are significant nitrogen sources for peatland ecosystems (
11).
Despite decades of research, there is still much debate on the identity of the dominant diazotrophs in ombrotrophic peatlands. Early work implicated
Cyanobacteria (
12–14) or heterotrophic bacteria (
15) based primarily on microscopic studies, while more recent molecular analyses argue for the importance of methanotrophic
Beijerinckiaceae (
16) as major diazotrophs in
Sphagnum peat bogs (
17–20). Possible contributions from other potential diazotrophs, such as strictly anaerobic methanogenic
Euryarchaeota, remain unknown. However, it is quite possible that diverse diazotrophs exist within defined niches of peatland environments (
21).
Diazotrophy is catalyzed by the nitrogenase metalloenzyme, a complex of thee subunits (H, D, and K) that contains abundant iron as Fe-S clusters. This enzyme is extremely O
2 sensitive (
22) and must be protected from exposure to O
2 for diazotrophy to occur (
23). The most common form of nitrogenase, encoded by
nif genes, contains molybdenum (Mo) as its cofactor. When Mo is scarce, some species of
Bacteria and
Archaea express nitrogenases containing vanadium ([V]
vnf genes) or iron ([Fe]
anf genes) in place of Mo, but these “alternative” nitrogenases are less efficient than the Mo form (
24,
25). The most conserved nitrogenase gene,
nifH (
26), has become the marker gene of choice for environmental diazotrophy (
27–29). Phylogenetic studies show five
nifH clusters: aerobic bacteria (cluster I), alternative nitrogenases (cluster II), anaerobic bacteria and archaea (cluster III), uncharacterized sequences (cluster IV), and paralogs related to chlorophyll biosynthesis (cluster V) (
30). Because
vnfH and
anfH genes in cluster II cannot be differentiated by sequence alone, the D subunit (
nifD-vnfD-anfD) has become the preferred marker gene for studies of alternative nitrogenases (
31). Consistent with higher concentrations of V than Mo in most rocks (
32), microbes from diverse soils have been shown to contain
vnfD genes (
31,
33–37). Given that oligotrophic conditions dominate in peatlands, trace metals may limit diazotrophy. However, little is known about trace metal availability and the role of alternative nitrogenase pathways in ombrotrophic peatlands.
Similarly, methane monooxygenase ([MMO] the enzyme that catalyzes the first step of methane oxidation) occurs in particulate copper (Cu)-containing (pMMO) and soluble Fe-containing (sMMO) forms. While pMMO has more specific substrate requirements, pathways that employ sMMO can use a wider range of compounds (
38). Both forms of MMO are inhibited by acetylene (C
2H
2) (
39,
40). In organisms with both sets of genes, pMMO is expressed when Cu is abundant, whereas Cu limitation induces sMMO expression (
41). The dominant peatland methanotrophs in
Alphaproteobacteria and
Gammaproteobacteria tend to possess both MMOs (
42–46), although
Methylocella and
Methyloferula species containing solely sMMO have been isolated from peat bogs (
47–49). While most studies have primarily targeted the
pmoA gene (
43,
45),
mmoX genes and transcripts have also been reported from peatlands (
46,
50,
51), raising questions about the relative importance of each form for peatland methane oxidation.
The acetylene reduction assay (ARA) is commonly used as a proxy for diazotroph activity (
52,
53). This assay is effective for capturing the potential activity of diazotrophs that are not inhibited by C
2H
2, such as
Cyanobacteria and nonmethanotrophic
Proteobacteria (e.g.,
Bradyrhizobiaceae) (
54). However, a number of functional guilds of microorganisms, including methanotrophs, methanogens, sulfate reducers, and nitrifiers, are inhibited by C
2H
2 (
55–60). If these or other C
2H
2-sensitive microbes perform diazotrophy and/or provide substrates to other diazotrophs (see Fulweiler et al. [
61]), ARA may underestimate diazotrophy in that system. Thus, recent studies have shifted to tracking diazotrophy by incorporation of the stable isotope tracer,
15N
2 (
20,
21,
62,
63).
In this study of the S1 peat bog at the Marcell Experimental Forest in northern Minnesota, USA, dissolved macronutrients (NH
4+, NO
3−, and PO
43−) and micronutrients (Fe, Cu, V, and Mo) were profiled along with the community composition and abundance of diazotrophic microorganisms. We also performed separate laboratory incubation experiments to measure potential rates of ARA and
15N
2 incorporation to (i) assess environmental controls (light, O
2, and CH
4) on diazotrophy; (ii) quantify the effect of C
2H
2 on rates of diazotrophy and methanotrophy; and (iii) search for diagnostic markers for alternative nitrogenase activity, such as a low conversion factor of ARA to
15N
2 incorporation (
31) and C
2H
2 reduction to ethane (
64). Finally, we make recommendations on universal
nifH primers for amplicon sequencing and quantitative PCR based on our findings.
MATERIALS AND METHODS
Site description and sample collection.
Samples were collected from the S1 (black spruce
Sphagnum spp.) peat bog at Marcell Experimental Forest (MEF; 47°30.476′N, 93°27.162′W), the site of the DOE SPRUCE (
Spruce and
Peatland
Responses
Under
Climatic and
Environmental Change) experiment in northern Minnesota, USA (
95). The S1 bog is ombrotrophic and acidic (average pH, 3.5 to 4 [
66,
96]). Over the summer months, the water table is ±5 cm from the hollow surface (
66,
97). Dissolved O
2 levels decrease to below detection (∼20 ppb) within the top 5 cm of the bog. Three locations were sampled along S1 bog transect 3 (T3) at near, middle, and far sites (see Lin et al. [
98] for further details). Surface (0- to 10-cm depth) peat was collected from hollows dominated by a mixture of
Sphagnum fallax and
S. angustifolium and from hummocks dominated by
S. magellanicum. Peat depth cores (0 to 200 cm) were sampled from hollows where the water level reached the surface of the
Sphagnum layer.
Macronutrients.
Peat porewater was collected using piezometers from depths of 0, 10, 25, 50, 75, 150, and 200 cm. Piezometers were recharged the same day as collection, and porewater was pumped to the surface, filtered through sterile 0.2-μm polyethersulfone membrane filters, and stored frozen until analysis. Nitrate (NO
3−) and nitrite (NO
2−) were analyzed using the spectrophotometric assay described by García-Robledo et al. (
99). Ammonium (NH
4+) concentrations were determined with the indophenol blue assay (
100). Phosphate concentrations were measured with the molybdate-antimony ascorbic acid colorimetric assay (
101).
Micronutrients.
Peat porewater was collected from two locations in the S1 bog from cores at depths of 0 to 30 cm, 30 to 50 cm, and 100 to 150 cm by filtration through 0.15-μm Rhizon soil samplers (Rhizosphere Research Products). All plastics were washed with HCl prior to sampling; Rhizon soil samplers were cleaned by pumping 10 ml of 1 N HCl through them, followed by rinsing with ultrapure water until the pH returned to neutral (∼100 ml/filter). After collection, samples were acidified with 0.32 M HNO3 (Fisher Optima) and analyzed using a Thermo Element 2 high-resolution inductively coupled plasma mass spectrometer (HR-ICP-MS; National High Magnetic Field Laboratory, Florida State University). Initial analyses resulted in the frequent clogging of the nebulizer, likely due to the abundance of dissolved organic carbon. Therefore, samples were diluted 1:10 to minimize interruptions from nebulizer clogs. Concentrations were quantified with a 7-point external calibration using standards prepared in 0.32 N HNO3 from a multi-element standard mix (High-Purity Standards).
To generate an organic-free sample matrix suitable for ICP-MS analysis without contaminating or diluting the sample, subsequent samples were digested as follows: 1-ml aliquots of the porewater samples were heated in 15-ml Teflon beakers (Savillex) with 1 ml of 16 N HNO3 (Ultrex II, JT Baker) and 100 μl of 30% H2O2 (Ultrex II, JT Baker) for 36 h at 230°C in a trace-metal-clean polypropylene exhaust hood. The HNO3-H2O2 mixture oxidizes any dissolved organic matter (DOM) to CO2, but the resulting matrix is too acidic for direct ICP-MS introduction. Therefore, samples were evaporated to near dryness, resuspended in a 0.32 N HNO3 matrix suitable for ICP-MS analysis, and analyzed using the Element2 ICP-MS along with parallel blank solutions.
Quantification and sequencing of gene and transcript amplicons.
Peat was frozen on dry ice at the field site in July 2013 or in liquid N
2 after 7-day incubations at 25°C in the light under a degassed (80% N
2 plus 20% CO
2) headspace with 1% C
2H
2, with or without 1% CH
4, for June 2014 incubations (see “Acetylene reduction and methane consumption rates” below). DNA and RNA were extracted with Mo Bio PowerSoil DNA and total RNA extraction kits, respectively, as described by Lin et al. (
46). RNA was cleaned with a Turbo DNA-free kit (Ambion). Nucleic acid purity was analyzed for the 260/280 absorbance ratio (of 1.8 to 2.0) on a NanoDrop spectrophotometer. cDNA was synthesized using the GoScript reverse transcription system (Promega) according to the manufacturer's protocol.
Plasmid standards for qPCR were constructed according to the method of Lin et al. (
102). Primer pairs are given in
Table 1. The gene fragments of
nifH,
pmoA, and
mcrA for constructing plasmid standards for qPCR were amplified from genomic DNA of
Rhodobacter sphaeroides,
Methylococcus capsulatus strain Bath, and S1 peat bog peat soil, respectively. To prepare cDNA standards, plasmid DNA with a positive gene insert was linearized with NcoI restriction enzyme following the manufacturer's protocol (Promega) and purified by using a MinElute PCR purification kit (Qiagen). RNA was synthesized from the linearized plasmid DNA by using the Riboprobe
in vitro transcription system (Promega) followed by cDNA synthesis using the GoScript reverse transcription system (Promega) according to the manufacturer's protocols.
The abundance of functional gene transcripts was quantified in samples run in duplicates on a StepOnePlus real-time PCR system (ABI) using Power SYBR green PCR master mix. Reaction mixtures of 20 μl comprised 2 μl of template cDNA (10 to 100 ng/μl) added to 10 μl of SYBR green master mix, 0.5 to 1.6 μl of each forward and reverse primer (0.3 to 0.8 μM final concentration;
Table 1), and 4.8 to 6.5 μl of PCR-grade water. Samples were run against a cDNA standard curve (10
1 to 10
7 copies of plasmid gene fragment) on a StepOnePlus qPCR instrument with 96 wells with an initial denaturation step of 2 to 5 min at 95°C and 40 cycles of denaturation at 95°C for 15 to 30 s, annealing at 55 to 64°C for 30 to 45 s, extension at 72°C for 30 to 45 s, and data acquisition at 83 to 86.5°C for 16 to 30 s. To minimize the effects of inhibitors in assays, peat DNA was diluted to 1/40 of the original concentrations, and duplicate 20-μl reaction mixtures, each containing 2 μl of diluted DNA, were run for each sample. Functional gene and transcript copy numbers were normalized to the dry weight of peat or 16S rRNA transcript copies for incubation samples. Amplicons were sent to the University of Illinois at Chicago for DNA sequencing using a 454 platform. Raw sequences were demultiplexed, trimmed, and quality filtered in CLC bio software. The phylogenies of
vnfD-anfD sequences were inferred using the maximum likelihood method based on the Kimura 2-parameter model in MEGA5 (
103).
Acetylene reduction and methane consumption rates.
Samples of bulk peat (Sphagnum spp. and surrounding soil) were collected from 0- to 10- and 10- to 30-cm depths in September 2013, April 2014, June 2014, September 2014, and August 2015 and stored at 4°C until the start of laboratory incubations. Samples from the 0- to 10-cm depth were gently homogenized so as not to rupture Sphagnum sp. tissues, while peat samples from the 10- to 30-cm depth were fully homogenized. For each sample, 5 g of bulk peat was placed in 70-ml glass serum bottles, stoppered with black bromobutyl stoppers (Geo-Microbial Technologies) (pretreated by boiling 3 times in 0.1 M NaOH), and sealed with an aluminum crimp seal. Headspaces were oxic (room air, 80% N2 plus 20% O2) or degassed (100% N2 or 80% N2 plus 20% CO2) with or without 1% C2H2 or 1% CH4. Treatments were incubated for 1 week at 25°C in the light or dark. A gas chromatograph with a flame ionization detector (SRI Instruments) equipped with a HayeSep N column was used to quantify CH4, C2H2, and C2H4. Samples were measured for C2H4 production daily until C2H4 production was linear (∼7 days). Controls not amended with C2H2 did not produce ethylene (C2H4). Incubations of hollow peat from June 2014 incubated under an oxic headspace with and without 1% C2H2 were also monitored for the consumption of 1% CH4. Statistical analysis was performed with JMP Pro (v. 12.1.0) using the Tukey-Kramer honestly significant difference (HSD) comparison of all means.
15N2 incorporation rates.
In September 2014, samples were quantified for N2 fixation rates by 15N2 incorporation in parallel with ARA measurements. Incubations were set up as described above and supplemented with 7 ml of 98% 15N2 (Cambridge Isotope Laboratories, Tewksbury, MA, USA). After 7 days, samples were dried at 80°C, homogenized into a fine powder, and analyzed for N content and δ15N by isotope ratio mass spectrometry (IRMS) with a MICRO cube elemental analyzer and IsoPrime100 IRMS (Elementar) at the University of California, Berkeley, corrected relative to National Institute of Standards and Technology (Gaithersburg, MD, USA) standards.
Metagenomic analyses.
Metagenomes were generated in a previous study (
104). Diazotrophic and methanotrophic pathways were investigated using the following bioinformatics approaches. Briefly, Illumina reads were filtered by quality (Phred33 score threshold of Q25) using Trim Galore (Babraham Bioinformatics) and a minimum sequence length cutoff of 100 bp. The sequences were then queried using RAPSearch2 (
105) against the NCBI nr database of nonredundant protein sequences as of November 2013. Sequences with bit scores of 50 and higher were retained to determine the total number of functional genes for normalization across the different samples. The taxonomic composition of protein-coding sequences was determined based on the taxonomic annotation of each gene according to the NCBI nr taxonomy in MEGAN5 (
106; minimum score, 50; maximum expected, 0.01; top percent, 10; minimum complexity, 0.3).
To classify sequences by nitrogenase cluster type, genes were analyzed using BLASTX (E value, 0.1; bit score, 50) versus a custom
nifH database that includes a phylogenetic tree to distinguish the principal clusters (I, V, and III) in the
nifH phylogeny, as well as paralogous cluster IV
nifH-like sequences (
27). Abundances of
nifH genes from the four clusters were normalized to those of total protein-coding genes from RAPSearch2 output sequences. The relative abundance of particulate (
pmoA) versus soluble (
mmoX) methane monooxygenase was based on previous analyses reported by Lin et al. (
46).
Accession number(s).
Metagenomes were reported in a previous study (
104) and deposited in BioProject
PRJNA382698 (SAMN06712535-06712540).
pmoA cDNA amplicons were reported in a previous study (
43) and deposited in BioProject
PRJNA311735 .
nifH,
mcrA,
nifD, and
vnfD-anfD cDNA amplicons were deposited in BioProjects
PRJNA382268,
PRJNA382282,
PRJNA382288, and
PRJNA382295, respectively.