Three separate samples of CBM produced water were collected from a CBM well that is located within the Wall coal seam formation at a depth of 376 m below surface in the southern PRB, Wyoming, USA. The water sample was filtered through a series of 0.45 and 0.22 μm filters into gas tight Nalgene sample bottles, and acidified to pH <2 with 1% nitric acid (Trace metal grade, Fisher Scientific, Pittsburgh, PA, USA). Trace element concentrations in filtered water samples were quantified by inductively coupled plasma mass spectrometry (ICP-MS).

For the extraction of DNA 200 ml of water sample were drawn into a sterile 0.22 μm cellulose acetate filter (Millipore, Billerica, MA, USA) inside a presterilized 25-mm Swinnex filter holder, fast frozen on dry ice, and transported to the laboratory for storage at −80°C until further processing. Total DNA was extracted from the filters with the PowerWater DNA isolation kit according to manufacturer’s instructions (MoBio Laboratories, CA, USA).

Microcosms were initiated immediately after sampling from an actively gas-producing well stream by inoculating CBM produced water with pea-sized sub-bituminous coal sampled from the PRB. Ten replicates were set up with 30 ml of produced water and 5 g of coal in 120 ml anoxic serum vials capped with butyl rubber stoppers under a headspace of N₂/CO₂ [80:20 (v/v)] and incubated at 37°C without agitation. Yeast extract was not added, as coal is the only organic matter to be used for microbial growth. A separate set of microcosms was initiated from the same well stream by inoculating only CBM produced water in basal salt medium. All microcosms were set up anaerobically (Balch and Wolfe, 1976) in 160 ml of anoxic serum vials capped with butyl rubber stoppers. Upon arrival at the laboratory, various electron acceptors were tested and microbial growth and methane production were observed in microcosms incubated at 37°C without agitation. The electron donor H₂ (H₂/CO₂ (80/20 (v/v)) was replaced with 50 mM sodium acetate or 5 mM methanol with N₂/CO₂ as the headspace gas.

Long term methanogenic enrichment cultures were established from CBM produced water by transferring cultures grown on H₂/CO₂ [80/20 (v/v)] after 50–80 days of incubation during a period of 2.5 years. Enrichment cultures were established in 160 ml serum vials with H₂/CO₂ [80/20 (v/v)] and basal salt medium consisting of the following components; 15 ml of KH₂PO₄ (27.2 g/l), 15 ml Na₂HPO₄·2H₂O (35.6 g/l), 12.5 ml of a mixture of NH₄Cl (24 g/l), MgCl₂·6H₂O (8 g/l), and NaCl (24 g/l), and 40 ml of NaHCO₃ (100 g/l); 10 ml vitamin solution (Biotin 2.0 mg, Folic acid 2.0 mg, Pyridoxine-HCl 10.0 mg, Thiamine-HCl·2 H₂O 5.0 mg, Riboflavin 5.0 mg, Nicotinic acid 5.0 mg, d-Ca-pantothenate 5.0 mg, p-aminobenzoic acid 5.0 mg, Thioctic acid 5.0 mg/l). To this 1 ml of a well defined trace element solution (SL-10) was added, that contained (per liter) 1.5 g of FeCl₂·4H₂O, 70 mg of ZnCl₂, 100 mg of MnCl₂·4H₂O, 6 mg of H₃BO₃, 190 mg of CoCl₂·6H₂O, 2 mg of CuCl₂·2H₂O, 24 mg of NiCl₂·6H₂O, 36 mg of Na₂MoO₄·2H₂O, 15 mg of Na₂WO₄, 15 mg of Na₂SeO₃·5H₂O, and 10 ml of 25% HCl (Widdel et al., 1983). Resazurin (2.0 g/l) was added as a reducing indicator. The medium pH was adjusted to 7.2 ± 0.2 with dilute HCl. Cysteine at 0.5 mg ml⁻¹ HCl was added into medium as the sole reductant. The commonly used Na₂S·9H₂O was not preferred as a reductant to avoid a reaction with metals in the medium. It was observed that there is no significant effect on methane production without addition of Na₂S·9H₂O (data not shown).

In the first phase of this study, the effects of varying concentrations of trace elements on methane production were evaluated. Initially, 10% (v/v) actively growing log-phase enrichment was incubated in 160 ml of anoxic serum vials capped with butyl rubber stoppers with fresh medium described in the Section “Enrichments.” Yeast extract was not added, since it is a source of trace elements and its element content is unknown. The headspace of serum bottles was aseptically flushed with H₂/CO₂ [80/20 (v/v)]. No external electron acceptors were supplied and CO₂ was the sole electron acceptor. Increasing concentrations of trace element mix solutions (1×-TES, 2.5×-TES, and 5×-TES) were added into enrichment cultures prepared in triplicates. 1×-TES is the concentration of the SL-10 trace element solution for a mixture of eight trace elements (iron, nickel, cobalt, molybdenum, zinc, manganese, boron, and copper), which is given in detail above. Enrichments and abiotic controls were incubated for 6 weeks under strict anaerobic condition in the dark at 37 ± 1°C without agitation. Glassware utilized in the study was pre-cleaned and treated with 10% nitric acid over night and baked at 550°C for 6 h.

For the second and third phases of this study, 10% (v/v) actively growing log-phase enrichment was incubated into fresh medium free of trace elements three times over a period of four and a half months. The last transfer culture was used as seed for 1×-TES amended and unamended enrichment cultures. Medium preparation and incubation conditions were the same as in phase one. Enrichments were set up in 10 replicates and incubated over 5 weeks. The medium pH was adjusted to 7.2 ± 0.2 with dilute HCl, and checked every week. Every week two duplicate incubations were sacrificed for RNA isolation and total cell counting.

Total cell numbers were estimated with 4′, 6-diamidino-2-phenylindole (DAPI) nucleic acid staining (Becerra et al., 2009) in 1×-TES amended and unamended enrichment cultures and in CBM produced water samples. A sequence of 40 randomly chosen microscopic fields across each slide was examined under 1000× oil-immersion magnification using an epifluorescence microscope (Nikon Eclipse E400) equipped with a Hamamatsu digital camera to determine the average number of cells in the samples.

Methane accumulation in the headspace was monitored and quantified for all microcosms and enrichments by gas chromatography using a GC-17A (Shimadzu, Co., Kyoto) equipped with an Equity-1 column (30 m × 0.53 mm ID, 3.0 μm, Supelco, St. Louis, MO, USA) and a flame ionization detector with helium as the carrier gas. The injector and column temperatures were adjusted to 100°C and the detector was held at 125°C. Certificated standard CH₄ (Fisher Scientific, Pittsburgh, PA, USA) was used for calibration. Methane concentration in the headspace was calculated by comparing peak areas of sample and standard.

Total RNA was extracted from 1×-TES amended and unamended enrichment cultures at 2, 7, 14, 21, and 28 days using the method described by Chin et al. (2008) with some modifications. The complete procedure was performed at 2°C or on ice. A volume of 2.0 ml of enrichment culture was anaerobically removed from a serum bottle and transferred to 2.0 ml RNase-free screw-cap tubes. Tubes were centrifuged for 5 min and supernatant was decanted. Approximately 0.6 g of 0.1 mm diameter glass beads, 0.6 ml TPM buffer (50 mM Tris–HCl, 1.7% polyvinylpyrrolidon, 20 mM MgCl₂·6H₂O, in DEPC-treated water) and 2 μl yeast tRNA (mRNA carrier) were added to the cell pellet. The tubes were shaken upside down for 10 min, followed by ballistic cell destruction for 1 min at 4200 rpm using a mini-bead beater (BioSpec Products; Bartlesville, OK, USA). Cell debris and glass beads were concentrated by centrifugation for 3 min. The supernatant was transferred to a tube containing 1 μl RNase inhibitor (Ambion, Austin, TX, USA). The pellet was resuspended in 0.6 ml phenol-saturated lysis buffer [50 mM Tris/HCl (pH 7.0), 50 mM EDTA, 1% (w/v) sodium dodecyl sulfate (Ambion), and 6% water-saturated phenol. After an additional round of bead beating at the same conditions, the tubes were exposed to centrifugation for 5 min and the supernatant was pooled with the supernatant of the first round of bead-beating. A volume of 0.6 ml phenol (pH 4.3) was added to the pooled supernatant, the tubes were vortexed for 40 s and exposed to centrifugation for 5 min. The supernatant was extracted with phenol–chloroform–isoamylalcohol [pH 7.8–8.2; 25:24:1 (v/v)] and chloroform–isoamylalcohol [24:1 (v/v)]. The aqueous phase was transferred to a fresh 2 ml tube containing 1/10th volume of 3 M NaAc (pH 5.2) and 2 μl of linear acrylamide (5 mg/ml; Ambion), and filled with 100% cold ethanol up to a final volume of 2 ml. The tubes were invert mixed and contents concentrated by centrifugation for 30 min. The pellets were washed with 0.5 ml 70% cold ethanol, invert mixed, and again centrifuged for 10 min. The pellet was dried in a desiccator for 5 min and resuspended in water. RNA was treated with RNase free-DNase (Ambion) to remove DNA contamination. DNase treated RNA was further purified by precipitating with 1/10th volume 3 M NaAc (pH 5.2), and 4 volumes of cold ethanol followed by centrifugation for 20 min. The pellet was washed twice with 0.5 ml 70% cold ethanol, dried, and resuspended in water. Enrichment of mRNA was performed with DNase treated and purified RNA using the MICROB Express purification system (Ambion). DNA contamination was checked with agarose gel electrophoresis following reverse transcription-polymerase chain reaction (RT-PCR) by performing control experiments in which no reverse transcriptase was added to the isolated RNA before the PCR step. RNA concentration was determined by absorption at 260 nm with a Biophotometer (Eppendorf, Hamburg, Germany). Purified RNA was stored at -80°C.

The cDNA of mRNA for the gene mcrA was generated with the primer mcrA-rev, (Steinberg and Regan, 2008) with 0.5 μg mRNA as the template. Multiscribe Reverse Transcriptase™ (250 U; Applied Biosystems) was used according to the manufacturer’s instructions. The conditions used in the cDNA synthesis program were as follows: initial incubation at 25°C for 10 min, incubation at 37°C for 120 min, enzyme inactivation at 85°C for 5 s, followed by rapid cooling to 4°C. cDNA samples were stored at −20°C until use. Methyl-coenzyme M reductase (MCR) is the essential and unique enzyme for methane production. It catalyzes the final step of methanogenesis, which reduces a methyl group linked to coenzyme M to methane (Ellerman et al., 1988). Because this enzyme is found in all known methanogens and is absent from non-methanogenic Archaea and Bacteria (Chistoserdova et al., 1998; Thauer, 1998; Bapteste et al., 2005), the mcrA gene is used as a specific functional marker to detect methanogens (Weil et al., 1988; Hallam et al., 2003). Previous studies have suggested that the phylogeny of mcrA follows the 16S rRNA phylogeny within various environments, and thus has been used as an alternative phylogenetic tool to detect and identify methanogens (Springer et al., 1995; Lueders et al., 2001; Luton et al., 2002; Friedrich, 2005). The MCR operons are found in two forms; MCRI (mcrBDCGA) which encodes mcrA and MCRII (mrtBDGA or mrtBGA) which encodes mrtA. Although, MCRI exists in all methanogens, MCRII is only found in members of the orders Methanobacteriales and Methanococcales (Lehmacher and Klenk, 1994; Bult et al., 1996; Reeve et al., 1997).

Extracted genomic DNA was used as template for PCR amplification of 16S rRNA genes using archaea-specific primers ARC21F (DeLong, 1992) and 1492R (Lane, 1991). PCR cycles were as follows: 1 cycle at 94°C for 5 min, 30 cycles at 94°C for 1 min, 58°C for 1.5 min and 72°C for 2 min, and a final extension step at 72°C for 10 min. In addition, three mRNAs clone libraries were generated from mcrA-specific RT-PCR amplicons (491 bp length amplified with primers mlas and mcrA-rev; Steinberg and Regan, 2008) obtained from trace element amended (TES-T1) and unamended enrichment cultures (NO-TES) after 7 days of incubation, and from startup enrichment cultures designated as time zero (T0). The purified RT-PCR amplicons from enrichment cultures and Archaeal 16S rRNA genes from produced water were ligated into the pGEM^®-T Easy Vector (Promega) and transformed into Escherichia coli JM109. Positive clones were verified by PCR amplification with vector specific primers M13F and M13R according to the protocol supplied by manufacturer. Clones were randomly selected, and plasmid inserts were sequenced with the M13F primer by Sanger sequencing. Sequences were checked manually and edited for ambiguous bases using the software Chromas sequence viewer (Technelysium, Tewantin, Australia), and multiple sequence alignments were created using ClustalX v.1.83 (Thompson et al., 1997). Archaeal 16S rRNA gene sequences were assigned to an Operational Taxonomic Unit (OTU) using the software package DOTUR based on ≥97% similarity (Schloss and Handelsman, 2005). Phylogenetic and taxonomic analyses were conducted using a representative sequence from each OTU. Phylogenetic trees based on the mcrA amino acid sequence were constructed using the Neighbor-Joining method (Saitou and Nei, 1987) with the “percentage of acceptance mutations distance correction model (PAM)” provided by the software package MEGA4.1 (Kumar et al., 2008).

To estimate the confidence of tree topology, bootstrap resampling analysis for 1000 replicates was performed with the Neighbor-Joining method. After sequences were assigned to OTUs, the frequency data at species level were used to construct rarefaction curves and to calculate diversity and richness indices using DOTUR. Estimates of phylotype richness were calculated according to the Chao1 estimator that takes into account the number of different OTUs in a sample, and the abundance-based coverage estimate (ACE). Species richness and evenness were calculated using the Shannon diversity index (H_SHANNON) and the Simpson’s index (H_SIMPSON), which takes into account the level of phylotype dominance.

Dilution series of purified mcrA RT-PCR amplicons were used as calibration standards for the real-time PCR quantification. The purified RT-PCR products were quantified and prepared for serial dilution as described previously (Chin et al., 2004b). The detection limits were determined from two independent measurements with dilution series of purified RT-PCR products (standards) in real-time PCR (108–101 target molecules per reaction). The copy numbers of unknown samples were calculated after real-time amplification from the linear regression of the standard curve. The cDNA generated with mcrA gene-specific primers was quantified with real-time quantitative PCR using SYBR Green. The levels of mcrA transcripts (per μg mRNA) were determined in enrichment cultures for which methane production was measured in parallel. All reactions were carried out in 20 μl reaction volume containing 0.02 units of iProof High Fidelity polymerase (BioRad), buffer, 2.5 pmol of each dNTP, bovine serum albumin, and 20 pmol of each primer of the pair mlas-mcrArev (Steinberg and Regan, 2008). The temperature profile was composed of an initial incubation at 98°C for 1 min, followed by 35 cycles of 98°C for 10 s, 58°C annealing for 10 s, 65°C extension for 32 s, and a final extension at 65°C for 6 min. Dissociation curve analysis to detect the presence of primer dimers and non-specific amplification peaks was performed after the final extension. The size of the PCR products was checked with agarose gel electrophoresis, and the specificity of PCR products was verified by sequence analysis of the clone library. Quantitative analysis of the cDNA was carried out with the 7500 Real-Time PCR System (Applied Biosystem, Foster City, CA, USA) using 7500 Real-Time PCR System Sequence Detection software (v 1.3.1). The precision as well as the reproducibility of quantification were carefully optimized, and PCR products were checked for their correct lengths as described previously (Chin et al., 2004b).

The sequences for mcrA mRNA clones and for archaeal 16S rRNA gene clones generated in this study were deposited in the GenBank database under the following accession numbers: JQ917180 to JQ917194 for all mcrA sequences, and JQ917195 to JQ917203 for archaeal 16S rRNA sequences.

Statistical comparisons were performed using one-way ANOVA followed by the all pairwise multiple comparison procedures (Holm-Sidak method) on the probability level P < 0.001. Data analyses were performed using the software package SYSTAT 11 (Systat Software, 2004).

In situ conditions at well depth indicated a pH value of 6.7 and a temperature of 29.6°C. Physicochemical analysis of the CBM produced water included: dissolved inorganic carbon (DIC) 13.161 mm/kg; F, 0.68 mM; Cl, 0.31 mM; SO₄, 0.07 mM; Na, 6.24 mM; Mg, 0.12 mM; K, 0.13 mM; Ca, 0.26 mM; Ba, 0.365 mM; total cations at 7.11 mM, and total anions at 6.65 mM. Trace element concentrations (μg/l) were B, 32.05; Al, 0.100; V, 0.49; Mn, 6.84; Fe, 14.92; Co, 0.03; Ni, 0.76; Cu, 0.48; Zn, 21.61; As, 0.07; Mo, 0.62; Ag, 0.07; Ba, 329.100; Tl, 0.01; Pb, 0.02; and U, 0.01. The concentrations of trace elements in 1×-TES, which was the concentration of trace elements in the SL-10 solution, was far greater than in situ concentrations except for boron, which was at a lower concentration in trace element enrichment cultures (21.3 μg/l). The calculated concentrations of elements in 1×-TES enrichment cultures were 785 times more in Fe, 132 times more in Ni, four orders magnitude more in Co, 476 times more Mo, 73 times more Mn, 28 times more Zn and 31 times more Cu than in situ concentrations.

Total cell number in the CBM produced water sample was determined at 3.2 × 10⁵ by DAPI staining, which is within midrange of cell numbers observed in other deep aquifers (Pedersen, 1993). Significant amounts of methane were produced in microcosms inoculated with pie-size coals (9.2 ± 0.6 μmol/ml) after 6 months of incubation. This indicates that the methanogenic community had the ability to degrade coal into methane.

Since coal is potentially the source of trace elements (Jackson and Reddy, 2007), we did not use coal as a carbon source to test the effects of trace elements on methanogenesis. For controlled and reproducible conditions, the CBM produced water sample was supplemented with a methanogen specific medium and amended with electron acceptors such as acetate, methanol, or CO₂. After 5 months of incubation, methane production was detected only in microcosms with H₂/CO₂ (80/20 v/v) at the cumulative amount of 7.7 μmol/ml, and subsequent enrichments were continued based on this set of cultures.

An archaeal 16S RNA gene clone library was constructed to identify the methanogens in the CBM produced water sample. Nine OTUs were identified within the 67 clonal sequences based on a “cut-off” value of ≥97% nucleotide sequence identity. The archaeal community was dominated by the genus Methanobacterium sp. (60%) with 40% identified as uncultured archaea (Table 1). A total of 25 sequences were affiliated with Methanobacterium sp. GH (96–99% similarity) which was previously isolated from sediments of the Gahai Lake in China (Zhu et al., 2011). Fifteen sequences were affiliated with Methanobacterium sp. SA-12 (97%) found in a subsurface gold mine (unpublished data, GenBank accession number AY695843.1) and to an uncultured clone (97%) from a minerotrophic fen (Cadillo-Quiroz et al., 2008). Five clones were affiliated with an uncultured archaeon clone (93%) from a gas hydrate production research well, Canada (unpublished data); four clones were affiliated (92–97%) with clones obtained from the deep aquifer of gold mines, South Africa (Takai et al., 2001); 10 clones were affiliated (96%) with clones from deep subsurface sediments in the Pacific (Biddle et al., 2006; Zhang et al., 2010); and eight clones were affiliated with a clone (91%) from hydrothermal sediments in Japan (Nunoura et al., 2010).

The effect of increasing concentrations of trace elements (1×-TES, 2.5×-TES, and 5×-TES) on methane production in CBM produced water enrichments was tested weekly for a total duration of 42 days. The differences in the average values of cumulative methane production among the enrichment cultures amended with 1× and 2.5× trace element concentrations were statistically significant (P < 0.001), when compared with unamended enrichment cultures. At the end of the incubation period, the level of methane production was significantly lower (8.52 μmol/ml) in enrichment cultures without addition of trace elements than in enrichment cultures with addition of trace elements (11.7 μmol/ml; Figure 1A). This suggests that methanogenic activity became diminished under trace element limitation. It was also observed that there was a significant negative effect on methane production in enrichment cultures with higher concentration of trace elements (5×-TES) during 42 days of incubation (P = 0.822). Here we determined that the undiluted, 1×-TES concentration of trace element solution was the optimum concentration for methane production. Consequently, we used this concentration for further experiments to investigate the effect of trace element addition on mcrA transcript level and on the composition of methanogenic communities in the enrichment cultures.

Methane production rate during the first week in 1×-TES enrichment cultures, was the highest (2.57 μmol ml⁻¹ week⁻¹). While, methane production rate in unamended enrichment cultures reached its highest level at the end of the second week of incubation (2.64 μmol ml⁻¹ week⁻¹; Table 2). The cumulative methane production was significantly lower (P < 0.001) than 1×-TES enrichments over 6 weeks of incubation (Figure 1A). The methane production rate in 5×-TES enrichment cultures was low and remained stable during the first 3 weeks. However, after 6 weeks it reached the highest level (10.70 μmol ml⁻¹ week⁻¹; Table 2), which might be the result of adaptation over time of the methanogenic community to higher amounts of trace element concentrations.

Direct cell counts were obtained over 35 days to document cell growth in 1×-TES amended and unamended cultures (Figure 1B). Amended and unamended studies were started with the addition of 2.3 × 10⁷ and 2.5 × 10⁷ cells/ml, respectively. At the end of the incubation time, cell densities had increased to 9.61 × 10⁸ and to 4.26 × 10⁸, respectively. Total cell numbers were ∼125% higher in trace element amended than in unamended enrichment cultures after 5 weeks of incubation. Since these are mixed cultures, growth trends were not well resolved for lag, log, and stationary phases (Meyer-Dombard et al., 2012). Calculated generation times were ∼70 h during the first 168 h (day 7) for amended and unamended cultures. After 7–14 days, the generation times for amended cultures were ∼80 h and two times longer (∼159 h) for unamended cultures. Total cell numbers increased steadily during incubation time for unamended cultures with increasing doubling times ∼622, ∼791, and ∼1019 h for 14–21, 21–28, and 28–35 days, respectively. However, in amended cultures cell numbers fluctuated. Generation times were increased to ∼1281 h between 14 and 21 days, dropped back to ∼264 h between 21 and 28 days and were calculated as ∼1273 h between 28 and 35 days. Changes in pH values were observed in all incubations, but there was no significant difference between enrichment cultures during the entire incubation time (Figure 1B).

The levels of mcrA transcripts in 1×-TES amended and unamended enrichments were quantified in weekly time intervals using quantitative RT-PCR over 35 days of incubation. After 1 week, a quicker and stronger response was observed for mcrA transcript levels in amended enrichments compared to unamended enrichments. At this time point, the addition of trace elements enhanced mcrA transcript levels to its maximum (2.18 × 10⁷ copy numbers/μg mRNA), while the level of mcrA transcripts remained much lower in enrichment cultures without addition of trace elements (7.03 × 10⁶ copy numbers/μg mRNA). Unamended cultures reached their maximum levels at 14 days of incubation (1.3 × 10⁷), but it was still lower than amended enrichment cultures after 7 days of incubation (Figure 2). After 35 days of incubation, mcrA transcript numbers were relatively low, with 3.05 × 10⁴ and 1.05 × 10⁴ for amended and unamended enrichment cultures, respectively.

The increased levels of mcrA transcripts in response to trace element addition correlated positively with elevated rates of methane production (R² = 0.95), while a very weak correlation was determined in unamended cultures (R² = 0.69; Figure 3). Additionally, the rate of methane production in unamended enrichments was lower than in amended enrichments, which correlates with mcrA transcript levels over the first 7 days of incubation (Table 2; Figure 2).

Enrichments incubated for 1 week were chosen to construct clone libraries based on their highest methane production rates (2.57 μmol ml⁻¹ week⁻¹) and paralleled by highest mcrA transcript levels (2.18 × 10⁷ copy numbers/μg mRNA) as a proxy of highest methanogenic activity over 6 weeks of incubation (Table 2; Figure 2). Three clone libraries of RT-PCR amplicons of mcrA mRNA were constructed to compare the metabolically active methanogenic taxa in 1×-TES amended (designated T1:TES) and unamended (T1:No TES) enrichment cultures after 1 week of incubation, and in startup enrichment cultures (designated T0). Eighty-one unique sequences were clustered into individual phylotypes (OTUs) to calculate rarefaction and non-parametric estimators (Table 3). Methanobacterium formicicum- and Methanobacterium subterraneum-like sequences were detected in startup, amended, and unamended enrichment cultures.

Both richness and diversity of the methanogenic community changed during 7 days of incubation. Enrichment cultures with addition of trace elements had a slightly higher predicted methanogen richness and diversity than enrichment cultures without the addition of trace elements (Table 3).

The effects of trace elements on the relative abundance and distribution of methanogens were monitored. Addition of trace elements contributed specifically to the enhancement of M. subterraneum-like sequences over M. formicicum-like sequences (Figure 4). Amended and unamended enrichments were dominated by clones that grouped closely with the mcrA genes of M. subterraneum-like sequences (81.5% in amended enrichments and 46.2% in unamended enrichments) and M. formicicum-like sequences (18.5% in amended enrichments and 53.8% in unamended enrichments). The methanogenic community composition in the starting enrichment culture (T0) appeared to be unchanged during 1 week of incubation, and mcrA genes from M. subterraneum-like sequences comprised 46.2% and M. formicicum-like sequences 53.8% of the clone library, representing similar ratios between both taxa in enrichment cultures without addition of trace elements (Figure 4).

Rarefaction analysis based on 97% sequence similarity showed adequate sampling for the startup enrichment, whereas enrichments with and without addition of trace elements did not reach asymptotes (Figure 5). Non-parametric statistical analyses (Table 3) showed that the observed richness covered between 100, 86, and 100% of the estimated total methanogenic richness in startup enrichment, trace element amended and unamended enrichment cultures, respectively.

A phylogenetic tree based on amino acid sequences of the genes mcrA was constructed using a selection of sequences from cultured methanogens and environmental clones from the public database GenBank. Several clones from T1:TES, T1:No TES, and T0 clone libraries were found to be closely related (95–99%) with two sequences identified in subsurface shale environments (GenBank accession numbers EF117654, EF117630). Clone TES-3 was the most dominant clone in amended enrichments, representing 20 out of 27 clones. It exhibited 97% sequence identity with M. subterraneum which is the nearest isolated relative. Clone No-15 represented 11 of 26 clones and clustered most closely with M. subterraneum with 98% similarity. Only one clone from T1:TES and three clones from T1:No TES were closely related to unpublished clone rlm_673a (97%) which was identified in an anaerobic waste water digester (GenBank accession number JF460370.1). As many as 14 and 12 out of 26 clones in T0 clone libraries have 91 and 95–97% identities to the next closest pure culture, M. formicicum and M. subterraneum, respectively (Figure 6).

Microbial methane production in deep coal bed basins indicates the ability to sustain microbial life in these limiting environmental conditions (dark, anoxic, reduced, limited substrate replenishment, etc.). The associated physical, chemical, and biological restrictions may influence activity, function, and diversity of the microbial species in the anaerobic biodegradation of organic matter in coal. It is suggested that distribution and activity of the microbial communities are supported by water–coal interaction in subsurface environments (Green et al., 2008). The organic fraction (e.g., carbon, hydrogen, nitrogen, oxygen, and sulfur) and the inorganic constituents of coal (e.g., trace elements) support growth in this subsurface biosphere, and the geochemical potential determines the energy profile for microbial growth, maintenance, and survival under dormancy. Coal–water interactions are expected to provide trace elements for microorganisms (Finkelman, 1995). Yet, our water chemical analysis from the coal bed well in the PRB demonstrated that the amount of trace elements concentrations at this habitat was far lower than the concentrations required for optimum methanogenesis (Figure 1A).

Biodegradation of coal organic deposits into methane is limited to anoxic and reduced conditions, and by the interaction of at least four metabolically diverse groups of microorganisms, hydrolytic, fermentative, acetogenic, and methanogens (Wawrik et al., 2012). Although the rate-limiting step for coal biodegradation is the initial solubilization of a macromolecular, polycyclic, lignin-derived aromatic network into utilizable substrates by microorganisms (Fakoussa and Hofrichter, 1999), we argued that methanogenesis would be another limiting step in coal biodegradation to methane, if the amount of trace elements is limited. Methanogens are the terminal oxidizers in anaerobic biodegradation of organic matter and they are vital for the resultant production of methane (Strąpoć et al., 2008). This study provides the novel consideration that the availability of optimized trace element concentrations might have an influence on methanogenic community structures and on methanogenesis.

Trace elements in coal bed basins may be concentrated in three different ways: (1) concentration during the life of plants prior to coalification; (2) concentration during the decay of organic substances; and (3) concentration by aqueous solution after burial of plant materials (Goldschmidt, 1935). Trace element concentration should be unique in each coal bed deposit, since the composition of plant material and the water source for each coal bed basin can vary with the geological formation, surrounding geochemical conditions, and adjacent overlying mudstones, siltstones, and sandstones. This geological material might be the source of trace elements, which leach with the groundwater into coal bed basins. Furthermore, the source of CBM produced waters might be water accumulated during coal formations or derived from meteoric waters (Flores et al., 2008). Since coal has trace elements embedded within its structure at varying concentrations, we avoided to add coal to incubations to prevent uncontrolled trace element effects.

Optimal concentrations of trace elements exist for cellular activity, and concentrations beyond this optimum can have a toxic effect on cells (Swaine and Goodarzi, 1995). The results of this study demonstrated that an optimized concentration of trace elements added to the enrichment cultures from coal bed produced water is beneficial for methanogenic growth and activity. Some recent studies suggested that the addition of high amounts of trace elements leads to a decrease in methanogenic activity in environmental applications (Chakraborty et al., 2010; Pobeheim et al., 2010). Our results indicate that the lack of added trace elements (No-TES) to enrichment cultures limited methane production in a way similar to adding high concentrations of trace elements (2.5×-TES, and 5×-TES; Figure 1A).

In addition, the results from this study demonstrated that addition of trace elements to production water enrichments influenced methanogenic community structure (Table 3). This suggests that trace elements in coal bed basins may be an important factor influencing the methanogenic community structure in coal bed basins in situ. It also indicates that M. subterraneum-like group has a growth advantage over M. formicicum-like group in the applied concentration of trace elements after 7 days of incubation (Figure 4).

Trace elements, utilized by microorganisms, provide proteins with unique coordination, catalytic, and electron transfer properties. Shi et al. (2005) reported that almost one fourth of structurally characterized proteins required a coordinated metal atom. Cobalt (Co), nickel (Ni), iron (Fe), zinc (Zn), molybdenum (Mo), and/or tungsten (W) are known to be crucial for the activity of enzymes in methanogens (Takashima et al., 1990). The concentration and type of trace elements required for metabolism and growth may vary among different species of microorganisms (Basiliko and Yavitt, 2001). Burgess et al. (1999) reported that trace metals influence microbial waste degradation and species diversity within sewage sludge. Gough and Stahl (2010) have shown that individual microbial populations are selected by exposition to varying levels of metals in different regions of an anoxic freshwater lake. Chin et al. (2004) reported that structural changes in the methanogenic community lead to functional changes in methane production in rice fields.

The active microbial community can be represented by RNA rather than DNA analysis (Lloyd et al., 2010), and previous studies suggested a correlation between metabolic activity and the transcript level of key genes (Fleming et al., 1993; Nazaret et al., 1994; Chin et al., 2008). In fact, transcripts for mcrA can be employed to demonstrate that methanogens are metabolically active in the environment (Juottonen et al., 2008; Yuan et al., 2011). The results from this study suggest that addition of trace elements caused the methanogenic community shift in composition, and a range of mcrA transcript levels could be directly associated with a rate of methane production and thus with the activity of methanogens (Figure 3).

In this study, mcrA genes of M. formicicum and M. subterraneum, and the gene mrtA of M. formicicum were identified from both trace elements amended and unamended enrichment cultures (Figure 6). Since there is a high degree of DNA sequence conservation between mcrA and mrtA sequences (Luton et al., 2002), primers designed for detection of mcrA sequences also detect mrtA sequences. Since the gene mrtA has arisen by lateral transfer from the Methanococcales rather than gene duplication (Reeve et al., 1997), they form a distinct cluster on the phylogenetic tree closer to Methanococcales than Methanobacteriales mcrA sequences (Luton et al., 2002).

Molecular analyses and cultivation studies demonstrated that methanogens are one of the most abundant groups in subsurface environments such as petroleum reservoirs, faults, and marine sediments (Orphan et al., 2000; Shlimon et al., 2004; Moser et al., 2005). So far, very few methanogenic species were isolated from these environments (Kotelnikova et al., 1998; Lai and Wu, 2011). In this study, both results, the phylogenetic analysis of archaea and the methane production in microcosms with H₂/CO_2, supported that methane production dominantly occurs via the hydrogenotrophic pathway in the sampled coal bed well. Our result was supported Klein et al. (2008) who found Methanobacterium sp. in CBM produced water and also associated with coal samples from PRB, Wyoming. Furthermore, many other studies also reported findings of Methanobacteriales in a variety of deep subsurface coal basins such as the Jhari Gippsland basin, Australia (Midgley et al., 2010), Alberta coal beds of Canada (Penner et al., 2010), and coal fields located at Parbatpur in India (Singh et al., 2012). After ∼2.5 years of routine transfers of CBM produced water into methanogenic media, the enrichments were highly dominated by two species of Methanobacterium sp., which are M. subterraneum and M. formicicum. It is important to note that enrichment cultures successfully represented the majority of methanogens in the CBM produced water sample (Table 1). M. subterraneum was isolated from deep subterranean granitic aquifers, and its phenotypic and phylogenetic characteristics were described previously (Kotelnikova et al., 1998). M. formicicum was isolated from a domestic sewage sludge digest (Bryant and Boone, 1987) and in a municipal land field sample as an endosymbiont of a ciliate (Finlay and Fenchel, 1991). H₂/CO₂ and formate are the growth substrates for these two species but they cannot grow on methylamines, acetate, pyruvate, dimethyl sulfide, methanol, or other alcohols plus carbon dioxide. Growth of M. subterraneum strain A08 is inhibited by yeast extract (Kotelnikova et al., 1998). In the present study, H₂/CO₂, acetate, and methanol were each tested as growth substrate for microcosms, and H₂/CO₂ was found to be main carbon and energy source. Additionally, yeast extract, which was used in the initial enrichment culture, was not added to the trace element amended enrichments due to the unknown quantities of trace elements in it. Optimum temperature conditions for M. subterraneum strain A08 and M. formicicum are described as 20–40°C and 37–45°C, respectively. Optimum pH values for M. subterraneum strain A08 and M. formicicum are between pH 6.6–6.8 and pH 7.8–8.8, respectively (Kotelnikova et al., 1998). In this study, the pH value in enrichments was adjusted to 7.2 ± 0.2 and no differences were observed during 7 days of incubation at 37°C between trace element amended and non-amended enrichments. Among all growth factors in the enrichments such as temperature, pH, or added substrates, trace element addition was the only variable that could have an effect on methanogenesis.

There are still many unknowns regarding the effect of trace elements, due to their bioavailability and optimal concentration on facilitated methanogenesis, and their toxic concentrations for microbial growth and activity in coal bed basins in situ. Bioavailability and toxicity of trace elements in the environment is under the control of geochemical processes such as mineral dissolution, precipitation, and ion adsorption/desorption. In this study, highly enriched methanogenic cultures were used to investigate the effect of trace elements on methanogenesis. However, in the natural environment, abundance and concentration of trace elements might enhance the rate of carbon source degradation by methanogens either directly or in a more indirect way. For example, it might increase the specific metabolic activity of bacterial groups such as primary or secondary fermenters whose products are the substrates for methanogens. Further study should focus on effects of trace elements on bacterial activity in coal biodegradation pathways and their associated community compositions. Effects of individual trace elements on growth and activity of methanogens and on fermenters might provide further insights to enhance biogenic methane production in coal bed wells. Therefore, better understanding of trace elements as a limiting factor for methanogenic activity in a coal bed basin will have a broader impact on our knowledge about ecology and physiology of methanogens.