|Home | About | Journals | Submit | Contact Us | Français|
The transcription dynamics of subunit A of the key gene in methanogenesis (methyl coenzyme M reductase; mcrA) was studied to evaluate the relationship between process rate (methanogenesis) and gene transcription dynamics in a peat soil ecosystem. Soil methanogen process rates were determined during incubation of peat slurries at temperatures from 4 to 37°C, and real-time quantitative PCR was applied to quantify the abundances of mcrA genes and transcripts; corresponding transcriptional dynamics were calculated from mcrA transcript/gene ratios. Internal standards suggested unbiased recovery of mRNA abundances in comparison to DNA levels. In comparison to those in pure-culture studies, mcrA transcript/gene ratios indicated underestimation by 1 order of magnitude, possibly due to high proportions of inactive or dead methanogens. Methane production rates were temperature dependent, with maxima at 25°C, but changes in abundance and transcription of the mcrA gene showed no correlation with temperature. However, mcrA transcript/gene ratios correlated weakly (regression coefficient = 0.76) with rates of methanogenesis. Methanogen process rates increased over 3 orders of magnitude, while the corresponding maximum transcript/gene ratio increase was only 18-fold. mcrA transcript dynamics suggested steady-state expression in peat soil after incubation for 24 and 48 h, similar to that in stationary-phase cultures. mcrA transcript/gene ratios are therefore potential in situ indicators of methanogen process rate changes in complex soil systems.
In recent years, the analysis of bacterial gene expression in environmental samples has developed as an important approach for evaluating the in situ metabolic state of physiologically active microbial cells. Most studies of transcriptional activity have been restricted to environments with low concentrations of coextractable enzymatic inhibitors, organic radicals, or nucleases (21, 53, 56), possibly due to technical issues associated with the low efficiency of extraction and rapid degradation of environmental mRNA. More recently, attempts have been made to quantify the transcription of functional genes and correlated physiological processes in complex soil communities (40, 42). The properties of peat soils are particularly detrimental to high-quality nucleic acid extraction because of acidity (pH 3.8), high water content (80 to 90%), high concentrations of coextractable organic matter, high capacity for complexing nucleic acids (10), and high concentrations of stable free radicals (52) and phenolics (59). Not surprisingly, therefore, the transcriptional dynamics of the key microbial functional genes in peat ecosystems have not yet been investigated.
In peat soils, methane is produced by strictly anaerobic Euryarchaeota, mainly from CO2, H2, and acetate (13, 22, 27, 36, 37). The key functional enzyme in methane generation is methyl coenzyme M reductase (MR), which catalyzes the reduction of one methyl group bound to coenzyme M, releasing CH4 (57). MR is assumed to be specific for methanogens, and evolutionary phylogenies established from sequence analysis of the mcrA gene (which encodes subunit A of MR isoenzyme I) and of 16S rRNA genes are congruent, implying its conserved nature (30, 33).
Characterization of peat methanogen communities by analysis of the diversity of mcrA or 16S rRNA gene sequences (18, 23, 36, 55) is often carried out in conjunction with assessment of the methanogenic potential under laboratory conditions (12-14, 24) and/or field measurement of net methane emission rates (3, 19, 20, 22). Although molecular analysis of methanogen communities generates valuable information and possible correlation of sequence types with ecological niches (23, 35), this approach does not address the dynamics of methanogenesis, for example, during natural drying and rewetting cycles (26). Similarly, analysis of methanogen potential under laboratory conditions may be hampered by distortion of the natural processes leading to methanogenesis.
Recently, Cadillo-Quiroz and colleagues (6) quantified methanogens in peat soil by using real-time quantitative PCR (qPCR) amplification of 16S rRNA genes, and Colwell and colleagues (8) quantified the abundance of mcrA genes in marine sediments, extrapolating data to estimate in situ CH4 production rates. In situ methanogenic activity can potentially be assessed more directly by quantifying the expression of the functional genes encoding MR (47), while quantification of both genes and gene transcripts allows calculation of transcript/gene ratios (51, 61). Such ratios emulate transcript abundance per cell ratios that have been compared with growth and activity in methanogen pure-culture studies (31, 38, 43, 48), potentially allowing in situ analysis of growth state and activity.
The aim of this study was to assess the validity of this approach by determining the relationship of functional gene transcriptional activity and process rates in an unamended peat soil. Methane production rates were determined for peat soil slurries incubated at different temperatures and were compared with the transcriptional dynamics of the mcrA gene by reverse transcription-qPCR (RT-qPCR) and qPCR. A further aim was to assess DNA and mRNA recovery from peat and its influence on qPCR and RT-qPCR data, thereby enabling correction of biases with internal-standard recovery efficiencies.
Peat samples were collected from an ombrotrophic upland blanket peat at Lake Vyrnwy, North Wales, United Kingdom, in November 2007. Vegetation was dominated by immature mown Scottish heather (Calluna vulgaris) interspersed with cotton grass (Eriophorum spp.) on a peat moss (Sphagnum spp.) lawn. The water table was close to the peat surface. Soil samples were taken to a depth of 20 cm after removal of vegetation and roots. Peat was broken up into pieces of approximately 500 cm3 and stored in plastic bags at 4°C until further processing. pHH2O was 3.8 (standard error [SE], 0.05; n = 3). Water content, determined as weight loss after drying at 100°C for 48 h, was 82% (wt/vol) (SE, 0.30; n = 3). Samples were stored at 4°C for up to 6 months prior to analysis.
For preparation of peat slurries, large pieces of peat, stored for 2 to 6 months, were broken up, homogenized by blending (Philips HR1372/90 hand blender, United Kingdom), and gravimetrically adjusted to a water content of 90% (wt/vol). Homogenized peat was evacuated three times in an anaerobic jar at 900 Pa (Eppendorf 5301 diaphragm pump, Cambridge, United Kingdom), flushed with 99.9% N2 to establish anoxic conditions, and incubated at 4°C for 24 h to allow equilibration of residual traces of oxygen with the gas phase. Homogenized anoxic peat (15 g) was dispensed into 100-ml serum bottles in a glove box under constant N2 flow and sealed with butyl rubber septa. Serum bottles were again purged with N2 and incubated in quadruplicate at 4, 10, 15, 25, 32, and 37°C to assess temperature-dependent rates of CH4 production and mcrA transcription. To analyze mcrA transcription at constant CH4 production rates, triplicate peat slurries containing soil that had been stored for either 4 months or 6 months were incubated at 25°C for 16 days or 4 days, respectively, and were sampled with low or high frequency, respectively.
Methane production was assessed by direct injection of 0.5 ml of headspace gas, sampled using a 1-ml disposable syringe, into a Hewlett Packard 5890 series II gas chromatograph (Agilent Technologies, Stockport, United Kingdom) equipped with a flame ionization detector (FID-GC). Gases were separated with a 1.5-m Porapak Q 80:100 mesh column (Sigma-Aldrige, Poole, United Kingdom) with N2 as the carrier gas at a flow rate of 30 ml min−1. The GC injector, oven, and FID temperatures were 100°C, 50°C, and 350°C, respectively. Rates of methane production g−1 peat (dry weight) were calculated by linear regression of increases in CH4 concentration with time, using at least three consecutive measurements with a regression coefficient (r2) of >0.9. Temperature trial peat slurries were incubated for 6 days and analyzed every 2 days. Time course experiment results were analyzed after 0, 1, 2, 3, 5, 7, and 16 days or after 0, 6, 12, 24, 48, 72, and 96 h for peat stored for 4 or 6 months, respectively.
Peat slurry microcosms were sampled destructively immediately after GC analysis, and samples were frozen immediately in liquid nitrogen and stored at −80°C. Variation between batches during qPCR was avoided by restricting analysis to samples generating the maximum number of 96 qPCR reactions per single assay (including standards and controls). Thus, analysis was performed on slurries from peat stored for 2 months and incubated at 4, 15, 25, and 32°C, on slurries from peat stored for 4 months and sampled at 0, 1, 3, 5, and 16 days (for time course experiments), and on slurries from peat stored for 6 months and sampled at 0, 6, 12, 24, 48, and 72 h. Extractions were performed in duplicate per biological replicate as described previously (16, 41), with some modifications. Peat samples were semidefrosted on ice, and 0.5 g was transferred into prechilled, 2-ml screw-cap Blue Matrix Ribolyser tubes (Hybaid, Ashford, Middlesex, United Kingdom) preloaded with 0.5 ml cetyltrimethylammonium bromide extraction buffer for total nucleic acid extraction. Immediately before addition of cold Tris-buffered phenol and chloroform, 3.3 × 107 cells of Escherichia coli K-12 (NCIMB 10218) suspended in 10 μl of the RNA preservative RNAlater (Ambion, Applied Biosystems, Warrington, United Kingdom) was added to samples from the time course slurries constructed from peat stored for 4 months to provide an internal calibration standard for qPCR assays (see below). Disruption of cells was performed twice for 20 s each time, with intermittent cooling in a Hybaid Ribolyser (Thermo Scientific, Basingstoke, United Kingdom) at speed 4. After initial separation of the aqueous and organic phases, the aqueous phases were transferred to MaXtract high-density gel barrier tubes (Qiagen, Crawley, United Kingdom) and an additional phenol-chloroform extraction, followed by chloroform extraction, was carried out according to the manufacturer's instructions. Precipitation of nucleic acids was carried out for 2 h on ice with the addition of precipitation carrier linear acrylamide (final concentration, 15 μg ml−1; Ambion). Tubes were kept on ice at all times during extraction, and extracts were subdivided into two aliquots for preparation of DNA and RNA templates. DNA and RNA templates were further purified by using a NucleoSpin Extract II silica spin column purification system (Macherey Nagel, Düren, Germany) and an RNeasy minikit (Qiagen), respectively, in accordance with the manufacturers' protocols for high-yield purification. The latter included 10 μl β-mercaptoethanol (≥99.0%; Sigma-Aldrich) ml−1 RLT lysis/binding buffer. DNA residues were removed from RNA templates by digestion with 2 U Turbo DNase (Ambion) for 30 min at 37°C, with the addition of 40 U of the RNase inhibitor RNasin (Promega, Southampton, United Kingdom). DNase-digested RNA extracts were further purified by phenol-chloroform extraction followed by chloroform extraction in gel barrier tubes. Finally, RNA extracts were precipitated with twice the aqueous phase volume of 100% ethanol, 500 mmol ammonium acetate (Sigma-Aldrich), and 15 μg ml−1 linear acrylamide. Salts were removed by washing twice with 75% ethanol, and RNA was redissolved in 12.5 μl RNase-free H2O. DNA and RNA purity and concentrations were established by Nanodrop spectroscopy. DNA aliquots were adjusted to 2 ng μl−1 for qPCR, and RNA aliquots were adjusted to the highest possible common concentration.
cDNA was generated from 5 μl of RNA extract with random hexamer primers at a final concentration of 37.5 ng μl−1 with Superscript II reverse transcriptase (Invitrogen, Paisley, United Kingdom) in accordance with the manufacturer's instructions. The reverse transcription reaction mixture also included T4 gene 32 protein (25 ng μl−1; MP Biomedicals, Cambridge, United Kingdom), linear acrylamide (100 ng μl−1), bovine serum albumin (200 ng μl−1; Ambion), and 20 U of the RNase inhibitor RNasin. Generated cDNA was diluted 1 in 5 prior to use in RT-qPCR. For negative controls, aliquots of each RNA extract were pooled in batches of 10 extracts and diluted in proportion to the RT protocol. Performing RT-qPCR on negative controls confirmed the amplification of RNA templates free of DNA contamination. An additional control contained the RT mixture without the template to ensure that all reagents were free of contaminants.
mcrA DNA and mRNA gene fragment copy numbers were quantified by qPCR and RT-qPCR with the ML primer pair (24, 33). All qPCR reactions were carried out in 25 μl reactions with 5 μl of template DNA or cDNA added to 20 μl qPCR reaction mixture containing 2.5 U AccuSure polymerase (Bioline, London, United Kingdom), 1× AccuSure reaction buffer including 2 mmol MgSO4, 500 μmol of each dNTP and SYBR green I (10,000× concentration; Invitrogen) at a final concentration of 0.2×. To facilitate strand separation, dimethyl sulfoxide (5% final concentration; Sigma-Aldrich), betaine (0.65 M; Sigma-Aldrich), and 7-deaza-2′-deoxy-guanosine-5′-triphosphate (50 μmol; Roche Diagnostics, Burgess Hill, United Kingdom) were added to the qPCR master mix (39), and linear acrylamide (100 ng μl−1) and bovine serum albumin (200 ng μl−1) were added to increase enzymatic efficiency. The ML forward and reverse primers were each used at a 1.2 μM concentration. mcrA qPCR was carried out using a hot-start protocol, with an initial denaturation step of 15 min at 95°C and then 40 cycles of denaturation for 30 s at 95°C, annealing for 45 s at 55°C, and extension for 45 s at 68°C, followed by fluorescence quantification at the end of a 79°C step for 8 s to allow dissociation of possible primer dimers and unspecific amplification products. qPCR was finalized by an extension step at 68°C for 10 min and subsequent melt curve analysis for confirmation of amplicon specificity by verification of identical, clearly defined melting peaks and by amplicon size analysis using standard agarose gel electrophoresis. Standard curves for calibration of mcrA qPCR were created using triplicate 10-fold dilution series covering 7 orders of magnitude from 10 to 1 × 107 gene copies per qPCR reaction during each run. The standard consisted of a 782-bp PCR amplicon spanning the mcrA region of Methanosarcina barkeri Schnellen 1947 (DSM 804) created with primers M. barkeri-forward (CTCCATGGGTGAAATGCTTC) and M. barkeri-reverse (AAAGTTGAAGGGCAGGAGGT). M. barkeri primer sequences were aligned to the mcr operon sequence of Methanosarcina barkeri fusaro (NCBI accession no. NC_007355), using the eprimer3 primer prediction software program (49). qPCR reactions were performed in duplicate per DNA or cDNA template and were carried out using a MyIQ5 real-time PCR detection system (Bio-Rad, Hemel Hempstead, United Kingdom). DNA- and cDNA-associated qPCR-determined mcrA copy numbers were expressed as quantities ng−1 nucleic acid, g−1 soil (dry weight), or g−1 soil (dry weight) normalized to internal-standard-recovery efficiency (see below). A survey of available genomic sequences demonstrated the presence of one or two mcrA coding regions, but calculations assumed one mcrA gene copy per cell. Additionally, qPCR amplification efficiency of peat samples was assessed by subjecting serial dilutions of DNA and cDNA to mcrA qPCR with the ML primer pair as described above. Efficiency was calculated as the slope of linear regression of derived copy numbers against the dilution factor.
Slurry cDNA and DNA were also used to quantify recovery of an E. coli K-12 internal standard from slurries constructed from peat stored for 4 months. Coding sequences of E. coli K-12 substrain MG1655 (NCBI accession no. U00096) were subjected to BLASTn (1) searches against a database containing the genomic sequences of 12 different methanogens (NCBI accession no. NC_009515, NC_007955, NC_009635, NC_009135, NC_009634, NC_008942, NC_009051, NC_003551, NC_003552, NC_007355, NC_007681, and NC_000916) and three further environmental microbial genomes (NCBI accession no. NC_004757.1, NC_007614.1, and NC_002977). Ten E. coli K-12 coding sequences with no discernible sequence similarity were randomly selected as putative targets for qPCR internal standard assays. For each selected sequence, primers spanning 700 to 1,600 bp and 150 to 250 bp nested within were generated with eprimer3 for production of qPCR calibration standards and qPCR primer pairs, respectively. E. coli K-12 qPCR primers were screened for background PCR amplification with DNA and cDNA templates generated from unspiked peat, as described above, and for analysis of qPCR amplification efficiency on generated qPCR standards and on DNA and cDNA generated from mid-log-phase cultures of E. coli K-12. qPCR primers targeting 196 bp of E. coli coding sequence b0606, gene ahpF, demonstrated good product yield with E. coli K-12 DNA and cDNA templates, with no detectable PCR products on peat DNA and cDNA templates. To quantify internal-standard recovery, DNA and cDNA generated from the time course experiment and constructed from peat stored for 4 months were subjected to E. coli K-12 qPCR using E. coli K-12 qPCR primers b0606-qPCR-f (CGCATGACGTTGACTGAAAT) and b0606-qPCR-r (AGGATCTGACCACCAAAACG), each with the 1,554-bp E. coli K-12 b0606 amplicon as a calibration standard. E. coli K-12 b0606-f (TGCTCGACACAAATATGAAAAC) and b0606-r (TGGTGCGAATCAGGTAGTCA) were used to generate E. coli K-12 b0606 qPCR standards. To calculate E. coli K-12 RNA recovery efficiencies, E. coli K-12 cell suspensions were also subjected to nucleic acid extraction and reverse transcription as described above. E. coli K-12 cell suspensions were produced from mid-log-phase cultures grown for 6 h in Luria-Bertani medium. Pellets obtained by centrifugation of 2.5 ml culture at 10,000 × g for 15 min were resuspended in 1 ml of the RNA preservative RNAlater to generate internal standards or in 1 ml 4% formalin in phosphate-buffered saline for cell counts determined by DAPI (4′,6-diamidino-2-phenylindole) fluorescence-based counting (46) of triplicate samples, indicating 3.3 × 109 (SE, 1.9 × 108) cells ml−1. E. coli K-12 cell suspensions in RNAlater were combined and aliquoted to guarantee equal cell distribution and frozen at −80°C until further use. E. coli K-12 qPCRs were carried out as described above on peat slurry DNA and cDNA templates, E. coli K-12 cell suspensions, and qPCR standards and controls with a 0.4 μM final concentration of E. coli K-12 qPCR primers. The E. coli K-12 qPCR efficiencies were 84 to 90%, with linear regression coefficients (r2) of 0.993 to 0.997. To access the dynamic range of extraction and reverse transcription (42), serial dilutions of E. coli K-12 cell suspensions ranging from 3.3 × 108 to 3.3 × 105 cells per 10 μl were added in triplicate to 0.5 g of homogenized peat as described above, and nucleic acid extractions, qPCR, and RT-qPCR were carried out as described above.
CH4 was produced in all peat slurries, with rates spanning 3 orders of magnitude, from 1.1 to 595 nmol g−1 h−1 (Fig. (Fig.11 and and2).2). In the temperature experiment, CH4 production rates were highest at 25 and 30°C and there were no significant differences in rates between these two temperatures (P > 0.05), suggesting adaptation to temperate climate. Rates were lowest at 4°C, where increases in methane concentration were close to the detection limit of FID-GC. The temperature-dependent rate increase from 4 to 15°C was approximately 6.2 nmol g−1 h−1 °C−1, with substantial increases to 15 and 25/30°C. The methane production rates in both time course experiments at 25°C (Fig. 2A and B) were significantly lower (42.5 [SE, 4.6] nmol g−1 h−1 and 2.0 [SE 0.2] nmol g−1 h−1 for the slurries constructed from peat stored for 4 and 6 months, respectively), equivalent to those at 15°C and 4°C in the temperature experiment, with significant lag periods (24 and 48 h, respectively) (Fig. (Fig.2A2A).
Total DNA and RNA yields were high in all experiments, with DNA concentrations within the range of yields reported for soils with high organic content (5). Raw extracts after initial precipitation were faintly discolored, with low A260/A280 ratios, indicating residues of coextracted humic acids. Additional purification with silica spin columns increased ratios to >1.7 for total RNA and DNA extracts. Nucleic acid extraction yields and RNA/DNA ratios varied significantly (P < 0.05) between the three experiments (Table (Table1).1). Mean DNA yields from temperature experiment peat slurries were approximately 10-fold higher than those from slurries of peat stored for 6 months, but total RNA/DNA yield ratios were higher in peat stored for 4 months at 4°C and highest in peat stored for 6 months. Amplification of qPCR calibration standards was linear (r2 > 0.994) in all qPCR assays from 102 to 106 copies of the template per reaction, and qPCR efficiency in all assays was between 85% and 95%, allowing confident enumeration of >10 mcrA copies ng−1 template. The peat sample qPCR efficiencies calculated from the slope of linear regressions of template DNA and cDNA dilution series were within the same range, suggesting unbiased enumeration.
mcrA amplicons were generated during qPCR and RT-qPCR from all extracts, and qPCR specificity was confirmed by fragment size evaluation with conventional agarose gel electrophoresis. To assess possible qPCR inhibition, RT-qPCR and qPCR were carried out on serial dilutions (10−1 to 10−3) of total RNA and DNA extracts, and the proportionality of the resulting mcrA copy numbers to the dilution factors was confirmed (50). Comparison of SEs of technical and biological replicates suggested greater biological than technical variance for DNA extractions but greater technical variance for mRNA extractions. However, differences in mean SE were not significant. The mean copy numbers across all samples were 2.4 × 103 and 2.5 × 103 ng−1 nucleic acid, equivalent to 2.9 × 108 and 1.4 × 108 g−1 soil, for DNA and cDNA, respectively. mcrA gene and transcript copy numbers increased three- and fourfold, respectively, with increasing temperature and were highest at 25°C, corresponding to highest CH4 production rates. However, mcrA gene (2.3 × 102 ng−1 nucleic acid) and transcript (2 × 102 ng−1 nucleic acid) abundances were significant in slurries incubated at 4°C, where rates were close to the detection limit. During time course experiments, gene and transcript abundances decreased with incubation time, such that the highest mcrA gene and transcript copy numbers were detected in slurries with low or absent discernible CH4 production rates (Fig. (Fig.22 and and3).3). The initial mcrA transcript and gene abundances for the time course experiment conducted with peat stored for 4 months were 25- and 10-fold higher, respectively, than those observed after incubation for 16 days and 6- and 4-fold higher, respectively, than those observed after 96 h in the experiment conducted with peat stored for 6 months (Fig. (Fig.33).
mRNA transcript/gene ratios were highest at 25°C, increasing on average 2.5-fold from 4°C to 25°C (Fig. (Fig.1),1), with a maximum 4-fold increase. During the time course experiment with peat stored for 4 months, the initial ratio was relatively high (1.2× mean) and increased by day 2 (2× mean) and then decreased (0.7× mean). Subsequent decreases were not significant (P > 0.05), indicating steady-state transcriptional activity (Fig. (Fig.2B).2B). Variation in ratios was greatest during early stages of incubation with peat stored for 6 months, with the highest ratios at zero hour (1.5× mean). Again, ratios were more consistent at 48 to 96 h, with less variation (0.7 to 0.85× mean) (Fig. (Fig.2A).2A). Mean transcript/gene ratios were similar in temperature and time course experiments with peat stored for 4 months (1.6 and 1.5, respectively) but lower with peat stored for 6 months (0.15), which also showed the lowest CH4 production rate (Fig. (Fig.11 and 2A and B).
Data were analyzed further by assuming that temperature experiment mcrA transcript/gene abundance ratios and time course ratios during the later time points (based on the mean ratio from at least three consecutive data points with the least variation from the mean) reflected steady-state conditions. mcrA transcript/gene ratios from all experiments increased with corresponding CH4 production rates over an 18-fold range, implying dependence of process rate on transcription ratio (Fig. (Fig.4),4), with a slope of 0.26, an intercept of 1.9, and a regression coefficient of 0.79.
To assess potential bias in determining gene and transcript abundances, extraction efficiencies were determined by calculating recovery for E. coli K-12 cells added to peat stored for 4 months. qPCR of E. coli K-12 cell suspensions yielded 1.4 × 106 DNA copies μl−1 (SE, 5.1 × 105), 50% lower than DAPI fluorescence-determined cell counts (3.3 × 106 μl−1), and 1.1 × 104 cDNA copies μl−1. Recovery of E. coli K-12 cell suspensions added to peat soil in a dilution series over 4 orders of magnitude (3.3 × 105 to 3.3 × 108) was linear (Fig. (Fig.5)5) (r2 > 0.99) for both E. coli K-12 genes and transcripts from 1.4 × 105 to 1.4 × 108 and 1.1 × 103 to 1.1 × 106 for DNA and cDNA, respectively. Slopes of linear regressions indicated recovery efficiencies of 62% for DNA and 92% for cDNA, suggesting unbiased analysis of transcript dynamics. However, recoveries from time courses with peat stored for 4 months were lower and varied between individual samples. The mean extraction efficiencies were 25% (SE, 7.2) and 14% (SE, 1.6) for DNA and cDNA, respectively. Applying the individual standard recovery correction factors to the respective mcrA transcript and gene abundances for this experiment increased the means to 9.1 × 108 and 4.8 × 108 mcrA copies g−1 soil, respectively. This reduced the mean mcrA transcript/gene ratio to 0.7, while maintaining a similar range of ratios from minimum to maximum (3.4-fold). However, temporal changes in ratios obtained from internal-standard-corrected values gave results that were similar to, if somewhat higher than, those normalized to quantities ng−1 nucleic acid (Fig. 2B and C and 3B and C).
Recently, Nicolaisen and colleagues (42) evaluated mRNA extraction for assessing transcription dynamics in an agricultural soil. Our findings support their observations that rapid processing and maintenance of RNase-free and constantly cold conditions are necessary for extraction of high-quality environmental mRNA. Additional purification steps were required for peat, due to the high content of coextractable humic acids, but the proposed inclusion of antioxidants in the extraction buffer (7, 28) did not increase mcrA transcript or gene abundances (data not shown).
Variation between batches was minimized by analysis of sample numbers that could be processed simultaneously. The recovery efficiencies for extractions with added E. coli K-12 cell suspensions were constant over 4 orders of magnitude, and efficiency was high (>90% for mRNA), suggesting maximal recovery or possibly overestimation by qPCR (Fig. (Fig.5).5). All extractions were carried out using the same protocol and reagent lot by one person. Nevertheless, internal-standard-recovery efficiencies varied between batches of extractions and between individual extractions. In addition, the E. coli K-12 extraction efficiencies for DNA and mRNA obtained from long-term-incubation samples were significantly lower than the efficiencies for serial dilutions but were within the range of reported cell recoveries from soil (e.g., 31% DNA recovery from agricultural soil ), suggesting sensitivity to minor, unrecognized changes. Little is known about mRNA extraction efficiencies from environmental samples, but previous studies suggest that recoveries were significantly higher than reported previously (34, 53). Plasmids could be used as internal DNA standards, rather than cells (9), and T7 polymerase-amplified amplified RNA could also be used to calculate recovery efficiencies for mRNA. The addition of E. coli K-12 cell suspensions in RNA preservative RNAlater to peat prior to extraction allowed quantification of DNA and mRNA recovery efficiencies for each extraction. Efficiencies are thus less likely to be affected by differences in extraction of intracellular and extracellular nucleic acids. However, compared to that of transcript/gene ratios normalized to nucleic acid levels, the influence of recovery efficiencies was low because of similar DNA and mRNA recovery efficiencies. Confirmation of transcription dynamics by two independent normalization methods thus increases confidence in the validity of mRNA analysis in soil. Nevertheless, maximum ratios corrected for recovery efficiency were also significantly underestimated in comparison to what was found for stationary-phase Methanococcus vannielii cells (17). The mcrA gene abundances were similar to those obtained for an acidic peatland soil in which methanogen rRNA genes constituted approximately 55% of the total archaeal abundance (6), suggesting similar dominance of methanogens in the investigated peat soils.
The ratio of mRNA transcript/gene abundance has been suggested as a more direct measure of transcriptional activity than absolute abundances (42, 61). Relationships between mcrA transcription and methane production rates in three sets of peat soil slurries were studied. Although transcript abundance did not correlate well with process rate, there was evidence for linear correlation (regression coefficient r2 = 0.79) between methane production rate and transcript/gene ratio, whose use may overcome variation in reproducibility and differences in cell abundance (42, 61). This approach also allows some comparison with pure-culture studies of MR transcription and activity (15, 17, 43, 45). In batch cultures of the mesophile Methanococcus vannielii (17), 180, 200 to 400, and 50 mcrA mRNA transcripts cell−1 were observed during early lag phase, initial exponential growth, and stationary phase, respectively, suggesting increase in MR content during exponential growth and reduction in stationary-phase cells. These values are more than an order of magnitude greater than those obtained with the peat slurries toward the end of the incubations, when physiological conditions are likely to be comparable to those in stationary phase. Analysis of internal-standard recovery (Fig. (Fig.5)5) suggests that DNA and mRNA recovery efficiencies from peat are similar, as found for tfdA transcript and gene abundances in agricultural soil (42). The difference in ratio may therefore result from a higher proportion of dead or inactive methanogens in peat soil, with low levels of mcrA mRNA.
In vivo MR specific activities vary between strains and growth conditions. For example, Methanothermobacter thermautotrophicus, Methanimicrococcus blatticola, and Methanococcus maripaludis have specific activities of 3 to 5, 1.6, and 0.6 μmol CH4 min−1 mg−1 cell protein, respectively, (11, 32, 54). MR content will also vary with cell size (4, 60, 62). Temperature-dependent CH4 production rates were calculated per peat methanogen cell (μmol CH4 min−1 mcrA DNA−1), and temperature-dependent in vitro MR specific activities for the thermophile Methanothermobacter marburgensis (15) were used to estimate the MR cell content for peat methanogens required to achieve observed CH4 production rates (Fig. (Fig.6).6). MR was assumed to constitute approximately 6% of the cell dry weight for M. thermautotrophicus (2), equivalent to approximately 46,000 MR molecules per cell (assuming 300 kDa for MR, a cell size of 0.475 μm in diameter by 5.0 μm in length, and 50% protein content [dry weight] [2, 29, 62]). The recovery efficiency of methanogens was assumed to be 25% during qPCR, and 10% of the total population was assumed to be active, as indicated by internal-standard recovery and low transcript/gene ratios. Calculations using these assumptions predict that the number of MR enzymes required for observed CH4 production rates at 25°C exceeds the maximum number possible, even for larger cells (e.g., 250,000 for Methanosarcina balticum) or that the specific MR activity was ≥35 μmol CH4 mg−1 MR min−1 at 25°C. This suggests that MR specific activities at low temperatures in peat methanogens are significantly higher than those in M. marburgensis and/or that the number of active cells was underestimated by qPCR and extraction efficiency correction. The steep increase in extrapolated MR enzymes with increasing transcript/gene ratios (Fig. (Fig.6)6) further suggests that MR cell content was not maximal, probably due to substrate limitation, at the highest observed CH4 production rates in 25°C and 30°C peat slurries during the temperature experiment, even though the rate at 10°C was higher than those for a Finnish oligotrophic fen (14). These data indicate that the highest possible mcrA transcript/gene ratios were probably not reached during this experiment.
In peat slurries incubated at 4°C and in time course experiments with peat stored for 6 months, mcrA transcript abundance was significant at CH4 production rates close to the detection limit and were high at the beginning of the time course experiments. The half-life for mcrA at 37°C is 15 min (17). Temperate peat methanogens therefore maintain relatively high mcrA transcript abundances when activity is low or undetectable, as found in the mesophile M. vannielii, but not during early growth phase in the thermophile M. thermautotrophicus (43, 45). Initial time course samples were taken after exposure to ambient temperatures during preparation, which may have induced transcriptional activity in 4°C peat slurries, increasing initial ratios. Peat used for temperature and time course experiments had been stored for 2 and 6 months, respectively, and although gene and transcript abundances were similar, the latter exhibited significantly lower CH4 production rates and a longer lag period before detectable methane production. The low rates may have been due to other aspects of methanogen physiology following extended exposure to low temperatures, substrate limitation, and oxygen exposure but demonstrate a discrepancy between transcriptional and physiological activities. Nevertheless, while physiological activity in the separate experiments and slurries may have been affected by different parameters (temperature, substrate availability, physiological status, and pH change), the correlation of transcript/gene ratios with physiological rates suggests that ratios are indicative of in situ physiological activity, regardless of how it is induced.
In conclusion, the quantification of gene or transcript abundance alone allowed only limited insight into methanogen process dynamics in peat soil. The presence of dead, inactive, and dormant cells suggests that transcript abundance may more accurately reflect the functional activity of the peat methanogen community, and significant correlation was found between transcript/gene ratio and process rate. However, this approach for analysis of in situ function must be used with caution. MR specific activity is strongly temperature dependent (15), regulation of process rate may require transcription of several genes (25, 58), and ratios for natural ecosystems may not be directly comparable with those for pure cultures. Also, the approximation of in situ physiological MR activity and cell content requires detailed information, from pure-culture studies of representative acidophilic temperate peat methanogens (4), on mRNA transcript abundances and related enzyme activities under natural conditions. Additionally, the relationship of physiological rates and transcription per template ratio poses analytical problems for environmental applications. Relatively high “background” ratios equated to low physiological rates, while an 18-fold increase was associated with the highest observed process rates, and interpretation of transcript/gene ratios requires a better understanding of the mRNA regulatory system of the relevant organisms. Nevertheless, quantification of transcript/gene ratios provided a better indication of physiological activity than gene abundance. While this study identifies some potential difficulties in directly relating transcriptional and functional activities, the approach may have particular advantages in comparing in situ functional activities of different physiological groups within a soil community and in determining short-term responses to changes in environmental conditions.
This study was supported by funding from the NERC UK Population Biology Network (UKPopNet).
We thank the Royal Society for the Protection of Birds (RSPB) for allowing access to the Lake Vyrnwy special area of conservation.
Published ahead of print on 11 September 2009.