Previous Article | Next Article ![]()
Applied and Environmental Microbiology, December 2005, p. 8773-8783, Vol. 71, No. 12
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.12.8773-8783.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Fred J. Brockman,1
James K. Fredrickson,1
David L. Balkwill,2
Michael E. Dollhopf,2,
Barbara Sherwood Lollar,3
Lisa M. Pratt,4
Erik Boice,4
Gordon Southam,5
Greg Wanger,5
Brett J. Baker,6
Susan M. Pfiffner,7
Li-Hung Lin,8,
and
T. C. Onstott8
Environmental Microbiology Group, Pacific Northwest National Laboratory, Richland, Washington 99352,1 Department of Biomedical Sciences Biology, The Florida State University, Tallahassee, Florida 32306,2 Department of Geosciences, University of Toronto, Toronto, Ontario M5S 3B1, Canada,3 Department of Geological Sciences, Biogeochemical Laboratories, Indiana University, Bloomington, Indiana 47405,4 Department of Earth Sciences, University of Western Ontario, London, Ontario N6A 5B7, Canada,5 Earth and Planetary Sciences, University of California, Berkeley, California 94720,6 Center of Biomarker Analysis, University of Tennessee, Knoxville, Tennessee 37932,7 Department of Geological and Geophysical Sciences, Princeton University, Princeton, New Jersey 085448
Received 6 April 2005/ Accepted 15 August 2005
|
|
|---|
|
|
|---|
Although a reasonable understanding of the energetic foundations for marine and terrestrial subsurface sedimentary ecosystems to a depth of 500 m has been achieved (20, 55), mechanisms for the energetic support of igneous and metamorphic rock systems remain incompletely addressed. Best known are igneous rock-hosted aquifers at depths of <450 m, where chemolithoautotrophic subsurface ecosystems, founded upon biological methanogenesis or acetogenesis and supported by geological H2 production and inorganic C, are thought to be limited in photosynthetically derived organic C (12, 56, 68). Relatively little is known, however, about the distribution of such communities at much greater depths or whether microbial ecosystems of other deep rock types (e.g., sedimentary and metamorphic) operate with similar energy and C sources.
Microbiological sampling opportunities in the deep subsurface are extremely limited, and a persistent challenge remains the differentiation of indigenous from introduced microorganisms (52). Several studies, however, have shown that in mine boreholes, introduced microbes fail to persist (47, 57), especially when the unidirectional flushing of hot, anaerobic fracture water reestablishes a redox environment similar to that of the predrilling conditions (47). South African mines are unusual in that they routinely attain depths of >3 kilometers below land surface (kmbls), thereby accessing thermophilic, anaerobic environments distinct from more frequently studied, shallower subsurface aquifers (56, 68). Mine drilling programs occasionally intersect pressurized fluid-filled fracture systems, some of which flush for years at rates of liters min1. Once purged of drilling-associated microbial and chemical contamination, flowing boreholes can be viewed as artificial artesian fractures and functional conduits into native deep subsurface habitats.
Here we describe such an artesian borehole, which being situated >1.5 km from the nearest current or historical mining activity, enabled the sampling of minimally impacted strata. Originating at 2.716 kmbls and extending downwards to 3.360 kmbls, the borehole intersected three faults. The borehole was uncased from 10 m below the outlet and had continuously flowed over the 12 months following drilling, producing >2.8 x 106 liters (
2,400 borehole volumes) of 54°C, anaerobic water. Using a bailer, a profile of the fluid chemistry and microbial community structure was obtained from source faults deep within the borehole to the outlet.
|
|
|---|
![]() View larger version (43K): [in a new window] |
FIG. 1. (A) Plan view showing location of borehole D8A. The Driefontein property is outlined, and major geological features are labeled. Mined areas are hatched. The locations of the comparative boreholes mentioned in the text, W6-38-bh1 (4) and E5-46-bh1 (47), are shown. (B) Geological cross section intersected by borehole. The depth, in meters below the surface, is depicted along the vertical axis, and the angular downhole distance is also shown. Faults reported are denoted with arrowheads and dashed lines. Sampling intervals are denoted by circles.
|
Sampling.
The outlet sample was collected with an expandable steel packer (Cementation Engineers, Johannesburg, South Africa) customized with a nylon insert to eliminate metal-water contact. Water and gas samples were collected anaerobically using a Delrin plastic (DuPont) valved manifold (46). To obtain samples from depth, a chlorine bleach-disinfected bailer was deployed with a wire winch. The bailer was tested to 100 m using sterile rhodamine WT (1,000 mg liter1) as a tracer. Upon retrieval, gas accumulated below the upper ball valve, which was subsampled via a dedicated sampling port with a gas-tight syringe. Water samples for microbiological and geochemical analyses were collected anaerobically from the bottom of the bailer.
Filtered (0.22-µm nylon Acrodisk; Gelman) samples included anions, cations, and NH3, preserved as described previously (47). Bottles for unfiltered samples were preloaded with preservatives and capped without headspace. These included dissolved inorganic carbon (DIC) and [
13C]DIC (140-ml serum vials with 200 µl of saturated HgCl2), NH3 (140-ml serum vials with 50 µl concentrated H2SO4), HS (20-ml glass serum vials with 500 µl 2 M Zn acetate),
34S (60-ml plastic syringes with 500 µl 2 N Zn acetate),
2H and
18O (140-ml and 20-ml serum vials), dissolved gases (140-ml evacuated serum vials with 200 µl of saturated HgCl2), organic acids (45-ml serum vials), dissolved organic carbon (DOC) and total organic carbon (TOC) (45-ml glass vials with Teflon septa containing 2 ml concentrated HCl), and HPO42 (acid-washed 250-ml Nalgene bottles). The anion, HPO42, and organic acid sample bottles were partially filled and frozen. The
34S samples were refrigerated. Due to volume limitations, some samples from the bailer (anions, cations, DIC, NH3, and HPO42) were present in reduced amounts.
On-site measurements.
Conductivity, pH, temperature, and Eh were measured at the borehole outlet with hand-held probes (HI98130 and HI98201; Hanna Instruments, Woonsocket, RI). Samples were verified as anaerobic with Chemets dissolved O2 kits (0.1 to 1 ppm [0.03 to 0.3 mM] K-7501; Chemetrics, Inc., Calverton, VA). Downhole temperatures were obtained by attaching eight-point temperature labels (range, 37 to 65°C; EW-08068-20; Cole Parmer) to the bailer. These strips were certified accurate to ±3% by the manufacturer and agreed to within ±1°C with meter and glass thermometer readings at the outlet.
Laboratory procedures.
Anions and organic acids were measured by ion chromatography (DX-320; Dionex, Sunnyvale, CA). The HPO42 concentrations were determined using the MAGIC method (11). Cations were measured by inductively coupled plasma-atomic emission spectroscopy (ICP-AES) (ICP 4300 DV Optima; Perkin-Elmer, Wellesley, MA). NH4+ concentrations were determined by nesslerization (22). The HS contents were determined by two methods. For some samples (Table 1), HS was measured by filtering the ZnS precipitate from the water, dissolving it with 0.1 M HCl, and measuring the Zn content by ICP-AES. Alternatively, HS was determined gravimetrically from the ZnS precipitate in the samples for which the
-34S content was determined. The nonpurgeable TOC and DOC were measured from the acidified samples as CO2 generated by catalytic combustion (EPA method 415.1), using an infrared detector (Tekmar Dohrman DC-190). For DOC, the sample was centrifuged at 10,000 x g for 10 min, and the supernatant was analyzed. DIC was determined by adding phosphoric acid to the sealed serum vials, purging the samples with high-purity N2, and measuring the CO2 gas concentration by infrared absorption (LI-6252; LiCOR Biosciences, Lincoln, NE).
|
View this table: [in a new window] |
TABLE 1. Major anions and cations in borehole D8A
|
The dissolved gases were analyzed with a Kappa-5 reduced gas analyzer (RGA 5; Trace Analytical, Sparks, MD) equipped with a flame ionization detector (FID) and a Hg vapor detector (for H2, CO, and CH4). Gas separation occurred on individual columns for the RGA and FID, and the high-purity N2 carrier gas flow rate was 40 ml min1. Compositional analyses of gas samples were performed at the Stable Isotope Laboratory at the University of Toronto as described by Ward et al. (74). Dissolved gas concentrations were derived from the gas volume abundance, the ratio of water to gas flow rates, and Henry's law constants following the procedure of Andrews and Wilson (2). Losses during retrieval caused dissolved gas concentrations to be lower in the bailer samples than in the packed outlet sample (Table 2). For this reason, the packed sample values were used for thermodynamic calculations. Because the measured gas concentrations may be underestimates due to phase separation at depth within the partially depressurized fault zone (41), they are considered minimum estimates. The
2H and
18O for three samples were also determined, using the procedures of Ward et al. (74).
|
View this table: [in a new window] |
TABLE 2. Gas concentrations and stable isotope analyses
|
G) of about 20 kJ mol1 to capture energy from a given reaction (63), although microbial activity has been shown to proceed in the lab at 5 to 15 kJ mol1 with syntrophic microbial consortia (28, 33, 64). Finally, the maximum rate at which this energy could be accrued was set by the rate of diffusion of the rate-limiting reactant to the microorganism (50). This rate (mol cell1 s1) is approximated by 4
DrC, where C is the concentration of the rate-limiting reactant (mol kg1), D is the diffusivity of the reactant (cm2 s1), and r is the radius of the microorganism (determined by flow cytometry). The reactant diffusivity increases with temperature according to the Stokes-Einstein relationship, and the values used were from the work of Cussler (17). This assumption presumes that deep subsurface microorganisms are nonmotile, which is likely given their extraordinarily low rates of growth (36, 58). The free energy flux (J cell1 s1) for a specific microbial redox reaction is equal to 4
DrC
G. The rate needs to be at least equivalent to the microorganism's maintenance energy demand in order for the pathway to be viable. In the case of a mesophilic nitrifying bacterium, for example, this maintenance demand is on the order of 1016 J cell1 s1 (71).
Flow cytometry.
Samples were preserved in the field with 3.7% formaldehyde. The concentration and distribution of fluorescing particles and forward scatter and side scatter versus fluorescence intensity plots were determined on a FACScan flow cytometer (Becton Dickinson Inc., San Jose, CA). An internal standard consisting of 10 µl of 1.0-µm carboxylate-modified TransFluorSpheres (Molecular Probes, Eugene, OR) at 1.8 x 105 spheres µl1 was added to 1 ml of sample. Samples were treated with the fluorescent nucleic acid dye SYTO13 (Molecular Probes) at 5 nM. Untreated aliquots of each sample were run to correct for autofluorescence. The average particle sizes in different areas of the plots were estimated for both treated and untreated parameters using Cell Quest software, and cell numbers were determined by the differences between the treated and untreated analyses. The average cell size was estimated by comparison with singlet and doublet clusters of the fluorescent spheres.
Microscopy.
Samples for electron microscopy were collected in 2-ml cryovials and fixed in the field with 2% (vol/vol) glutaraldehyde. One hundred microliters of each sample was suspended in 5 ml of filtered (0.2 µm) deionized water (Barnstead E-pure) and passed through a 13-mm, 0.4-µm-pore-size Isopore membrane filter (Millipore). Samples were processed through an ethanol dehydration series (25%, 50%, 75%, 100%, and 100% [30 min each]), followed by critical point drying using a Samdri PVT 3B CPD apparatus (Tousimis Research Corp.). Gold-coated samples were examined using a Hitachi S-4500 field emission scanning electron microscope (SEM) equipped with an EDAX system at an accelerating voltage of 10 kV with a 15-mm working distance. A total of 771 and 1,606 cells were counted from the 125-m and 250-m samples, using fluorescence in situ hybridization (FISH) analysis according to the procedure of Baker et al. (3), with EUB338 and ARC915 probes.
PLFA analysis.
For phospholipid fatty acid (PLFA) analysis, biomasses from 8 liters of the outlet sample and 140 ml of the bailer samples were collected on Anodisc filters (0.2 µm by 47 or 25 mm; Laguna Niguel) (47). Compositional analyses of PLFAs were performed using procedures described by White and Ringelberg (75). Approximate bacterial cell densities were calculated using a conversion factor of 2.5 x 104 cells pmol1 PLFA (5).
Molecular analyses.
Biomasses from 7 liters of sample D8A-Outlet-110701 and 140-ml bailer samples from 390, 550, and 648 m downhole were collected and preserved as previously described (46, 47). Filters were frozen at 20°C in the field, transported on dry ice, and stored at 80°C. DNA was extracted using an UltraClean Soil DNA mega prep kit (MoBio Laboratories, Solana Beach, CA), with shaking at 70°C and vortexing with beads. Final extracts (6 ml) were concentrated by precipitation in isopropanol-sodium acetate (9:1), followed by a 70% ethanol wash and suspension in 25 µl 10 mM Tris, 0.1 mM EDTA (47).
16S rRNA genes were amplified using La Taq polymerase (Takara Mirus Bio, Madison, WI) and primers at 0.2 µM. To remove possible contamination prior to template addition, digestion with HaeIII (0.1 units) was performed (incubation at 37°C for 20 min, followed by inactivation at 80°C for 25 min). One microliter of environmental DNA extract was added, and the reaction was incubated at 95°C for 5 min (1 cycle); at 95°C for 30 s, 50°C for 1 min, and 72°C for 1 min 30 s (35 cycles); and at 72°C for 20 min (1 cycle). Products were purified using UltraClean PCR clean-up kits (MoBio). Bacterial 16S rRNA genes were amplified using the primers 27f (10) and 1492R (39). Possibly due to the poor DNA extraction efficiency resulting from the gentle lysis procedure, reliable amplification of archaeal rRNA genes required a half-nested reaction in which the product generated using primers 21F (19) and 1492R (39) was used as the template for a second round of PCR with primers 21F and 958R (19). 16S rRNA gene cloning was done with TA cloning kits (Invitrogen). For each clone library, inserts were amplified from 96 clones, and the resulting products were screened by restriction fragment length polymorphism (RFLP) analysis. Representatives of each identified ribotype were purified using UltraClean PCR clean-up kits and then sequenced directly with vector primers (Big Dye reaction kit and ABI 3100 DNA sequencer; PE Applied Biosystems).
Chimeric sequences were identified using Bellerophon (29) and RDP ChimeraCheck (13). The closest available matches to nonchimeric sequences were identified using the NCBI BLAST search tool (1). Sequences were managed within ARB (http://www.arb-home.de) and aligned with an ARB database built upon the 2002 Hugenholtz RDP alignment, using version 3.01 within the ARB_EDIT_4.1 tool followed by manual adjustments. Clone types (99% similarity) were designated after the alignment of sequences with ARB and comparisons of unambiguously aligned positions, as determined using the Lane mask (39). Maximum likelihood trees were generated using the Phylip FastDnaMl algorithm. Parsimony analysis and the calculation of bootstrap probabilities were performed with PAUP (Sinauer and Associates, Sunderland, MA). Phylogenetic trees based on Bayesian inferences were generated using MRBAYES (30), with 100,000 generations using a gamma rate correction, a sampling frequency of 100, and a burn-in value set to omit the first 350 trees from the consensus.
Terminal restriction fragment length polymorphism (T-RFLP) analysis of community 16S rRNA sequences was performed as previously described (47). To relate individual T-RFs to library clones, the 16S rRNA clone sequences were digested in silico, and the predicted T-RFs were compared to those produced by T-RFLP using an empirically derived accuracy criterion (47).
Functional genes were amplified using the same reagents as those used for 16S rRNA genes. The gene encoding the
subunit of methyl coenzyme M reductase (mcrA) was amplified from the 648-m sample with the primers (ME1F and ME2R) and PCR conditions described by Hales et al. (26). The dissimilatory sulfite reductase gene (dsrAB) was amplified from the 550-m sample as described by Baker et al. (3). PCR products were gel purified and cloned using TOPO-TA kits (Invitrogen). From each of the libraries, 48 random clones were typed by RFLP analysis. Clones were aligned with published representative sequences using Clustal X and imported into ARB, where neighbor-joining trees (10,000 bootstraps) were constructed. dsrAB sequences were imported into ARB, and phylogeny was generated using maximum likelihood. Bootstrap values and Bayesian inferences were also used to statistically support the nodes.
|
|
|---|
Water age and source.
Isotopically, the water in borehole D8A is meteoric in origin but distinct from that of the overlying dolomite (Table 2; Fig. 2). The
2H versus
18O measurements agreed with those for fracture water from nearby 4 Shaft of Kloof mine F3, as reported by Takai et al. (70), for which an age of 3 to 30 Ma was determined (41). Because the average chlorinity of D8A water also matched that of the Kloof sample (27.8 mM) (41) but the He concentration in D8A was four times greater that in the Kloof borehole water, at 1.3 mM (not shown) versus 0.3 mM (41), the age of D8A water is greater than that of water from the Kloof borehole. Using the models of Lippmann et al. (41), the estimated age is 4 (U = 240 ppm) to 53 Ma (U = 18 ppm), based on the flux and in situ He models, respectively.
![]() View larger version (28K): [in a new window] |
FIG. 2. 2H versus 18O plot showing positions of three samples from borehole D8A (crosses) in relation to other samples from nearby mines (diamonds) (41), Vaal River water (open circle) (21), and Pretoria annual average precipitation (closed circle) (31). The local mean water line (LMWL) is depicted as a solid line (44), and the global mean water line (GMWL) is shown as a dashed line (16). The samples from D8A occupy the same positions as those from the Kloof Mine sample, for which an age of 3 to 30 Ma was determined (K4-41-FW1) (41). A sample from the Witwatersrand supergroup quartzite (W6-38-FW3) (41) also plots very close to the samples from borehole D8A, whereas samples from the overlying dolomite aquifer, mine service water (51), and a Ventersdorp supergroup metavolcanic sample (E5-46-bh1) (47) occupy different positions. The figure is based on the work of Lippmann et al. (41), with permission.
|
The
2H value for H2 (685
VSMOW [Vienna Standard Mean Ocean Water]), relative to that for H2O (26.9
VSMOW) (Table 2), yielded an isotopic equilibration temperature of 60.5°C (8, 43), which is significantly greater than the 54°C measured downhole (Table 1). Based upon heat flow measurements obtained from West Witwatersrand line quartzite (49 ± 3 mW m2) (34) and a one-dimensional thermal conductivity model (49), this temperature constrains the water source to between 3.8 and 5.0 kmbls. The >2.8 x 106 liters expelled prior to sampling would fill the volume of a hypothetical fracture with an average width of 0.5 cm and an area approximated by a plane of 250 by 2,035 m, which is compatible with the estimated source depth of the fault water. The fact that neither surface nor mine water breakthrough had occurred indicates that D8A is linked to a voluminous, deep fracture system with minimal connection to shallower sources. A similar inference was reached for the 3.1-kmbls Kloof Mine shaft 4 borehole water (70), whose 60°C temperature, exceeding the predicted in situ rock temperature by
14°C, was attributed to upward migration along interconnected fractures (49). Collectively, the aqueous geochemistry and isotopic data are best explained by a deep source of ancient meteoric water. Since the water bears no evidence of interaction with the Ventersdorp supergroup metabasalt, the flow path must be confined to the underlying Witwatersrand supergroup quartzite.
Microbial abundance and community structure.
PLFA analysis indicated a bacterial biomass ranging from
6 x 103 to
2 x 104 cells ml1, with a possible increase (about three times) at shallower depths (Table 3). Total microbial cell densities given by flow cytometry were
1 x 104 to 4 x 104 cells ml1 for samples from the 250-m depth to the outlet but were
50 times higher for the two deepest samples (
1 x 106 cells ml1) (Table 3). Since PLFA is a measure of bacterial biomass, the discrepancy between these methods probably reflects an enrichment of archaea in the deeper samples. FISH analyses of the 125-m and 250-m samples revealed that 9% and 10%, respectively, of the cells were filamentous archaea, with the remainder being rod-shaped bacteria (Table 3). This observation is consistent with the PLFA cell densities being slightly less than the flow cytometry cell densities. Flow cytometry also indicated that most of the observed increase in total cell counts with increased depths is attributable to weakly stained, large-diameter (2 to 11 µm) particles (data not shown). SEM observations indicated that these particles were comprised of Si and Fe mineral colloids with attached filamentous microbial cells (data not shown). This colloidal material was visually apparent as a darker cast in the deeper samples. It was impossible to obtain accurate FISH counts in the deeper samples due to interference by mineral colloids, but it was qualitatively apparent that the filaments far outnumbered the rod-shaped bacterial cells that dominated the shallower samples.
|
View this table: [in a new window] |
TABLE 3. Cell number estimates by flow cytometry and phospholipid fatty acid analysis
|
|
View this table: [in a new window] |
TABLE 4. Clone library summary
|
![]() View larger version (37K): [in a new window] |
FIG. 3. Phylogenetic analyses. Clones from Driefontein Mine boreholes are depicted with bold lettering (47). Those specific for this study are underlined. (A) Phylogenetic dendrogram based on a maximum likelihood analysis (Phylip FastDnaMl) of 862 unambiguously aligned positions of 16S rRNA gene sequences. (B) Maximum likelihood tree of dissimilatory sulfite reductase (dsrAB) genes. (C) Phylogenetic tree based on neighbor-joining analyses of methyl-coenzyme M (mcrA) reductase gene sequences. Maximum likelihood tree topologies were compared against trees generated by Bayesian inferences and parsimony. Closed circles indicate nodes conserved in all three treatments, with Bayesian probability and parsimony bootstrap values of >75%. Nodes conserved in the majority of tree topologies and with 50 to 74% probabilities and bootstraps are designated by open circles. Branch points supported by fewer than two trees or having <50% probability and bootstrap values have no designation. In the neighbor-joining tree, branch points supported by >79% bootstrap values are designated by solid circles, those supported by 60 to 79% bootstrap values are designated with open circles, and those having bootstrap values of <60% are unlabeled. Bars, changes per nucleotide.
|
Key functional genes presumed to be involved in the energy metabolism of Desulfotomaculum (the dissimilatory sulfite reductase gene, dsrAB, for SO42 reduction) (73) and of Methanobacteriaceae (the methyl-coenzyme M reductase gene, mcrA, for methanogenesis) (26) were also detected. A single dsrAB clone type (Table 4; Fig. 3B), which was most closely affiliated with Desulfotomaculum spp., and two mcrA clone types (Table 4; Fig. 3C) from within the Methanobacteriaceae were identified, confirming the genomic potential for SO42 reduction and methanogenesis.
The microbial community structure of D8A fault water is therefore populated by a very simple microbial community dominated by a single bacterium (probably a sulfate reducer) and several closely related methanogens. Whether the methanogens or sulfate reducers dominate the microbial community within the fault is impossible to ascertain, because the filamentous methanogens may have preferentially adhered to the mineral colloids that formed and then settled to the bottom of the borehole as the hot fracture water cooled upon its ascent to the borehole outlet. This phenomenon would have gone undetected if not for the depth profile results. Colloid formation in deep artesian boreholes may be a common occurrence that needs to be considered if the emanating fluids are to be used to quantify the microbial community structure, as they have been in previous reports (12, 15, 25, 47, 68).
Predicted in situ microbial reactions.
No significant difference was observed between DOC and TOC values, nor did either exhibit any trend with the borehole depth (Table 1). Acetate was the most abundant organic acid, ranging from 0.013 to 0.040 mM. The
13C value for DIC ranged from 13.7 to 18.3
Vienna Pee Dee Belemnite, averaging 15.9
, with the deepest samples showing the least isotopic depletion (Table 2). Dividing cells were not observed by SEM or FISH, suggesting that the microorganisms were not growing at the time of sampling. The mM HS versus µM SO42 concentrations were uniform along the length of the borehole (Table 1), indicating that SO42 reduction had occurred within the fault system rather than the borehole water column.
These observations, combined with the short residence time in the borehole (ca. 3.5 h), argue for a steady-state balance between the observed water chemistry and the microbial population structure within the source fault water. To determine if the calculated Gibbs free energy of reaction (
G) and the steady-state free energy flux favored SO42 reduction and methanogenesis over other possible physiologies, calculations were performed using the borehole water chemistry for 83 possible microbial reactions at 43°C and 54°C and the measured outlet pH of 9.0. Calculations were also performed assuming a temperature of 61°C and a pH of 7.8, the predicted conditions for the fault water deduced from the borehole outlet pH and DIC after adjusting for the degassing of CO2 within the borehole. The microbial reduction of Fe3+ was excluded due its reactivity at high levels of HS. The free energy flux calculations predicted that the three most thermodynamically favorable reactions were methanogenic (Table 5, reactions 1 to 3), with two reactions utilizing CO. Reactions 1 and 2 in Table 5 are known to be catalyzed by Methanobacterium spp. (69).
|
View this table: [in a new window] |
TABLE 5. Free energy and steady-state free energy flux for possible microbially mediated reactionsa
|
Desulfotomaculum and Methanobacterium spp. are known to intimately interact during the transfer of H2 and/or formate in thermophilic (55°C) anaerobic fermentation reactor biogranules utilizing low-molecular-weight, reduced organic compounds (32, 59). The high H2 concentrations in the D8A water would seem to preclude this process, and cell assemblages, where low H2 can be maintained in microsites by consumption rates that exceed diffusive flux, were not observed microscopically. Similarly, anaerobic CH4 oxidation coupled to SO42 reduction (Table 5, reaction 13) was less favorable than the individual methanogenic and SO42-reducing reactions. The
13C value for the
121 µM DIC was greater than that for the mM CH4, which would not be the case for anaerobic CH4 oxidation by metabolic coupling between a SO42 reducer and a methanogen. Furthermore, molecular evidence for predicted anaerobic CH4 oxidizers (e.g., ANME-1) was not detected in either the 16S rRNA or mcrA clone library.
Microbial activity.
The
34Ssulfide value was enriched in 34S (average, 12.1
Cañon Diablo troilite [CDT]) (Table 2) relative to the sulfide value and more consistent with sulfate values that have been reported for the open borehole E5-46-Bh1 (i.e.,
34Ssulfide = 1.4
CDT and
34Ssulfate = 19
CDT for E5-46-Bh1) (47). The sulfide isotopic enrichment and low SO42 concentrations are consistent with microbially driven fractionation under SO42 limitation, with minimal oxidation of the product sulfide. This is consistent with a hydrologically isolated environment and suggests that the microorganism represented by the D8A-1 clone performs sulfate reduction within the fault zone. The isotopically depleted
13C value for the DIC indicates that some of this sulfate reduction is tied to the oxidation of organic C, most likely acetate.
The
13C values for hydrocarbon gases trended towards greater depletion with increasing carbon chain lengths (C1 versus C2 + C3), whereas the
2H became enriched from C1 to C2 (Table 2), both of which are indicative of abiogenic production (66, 74). If biological methanogenesis were occurring within the faults, then the concentrations of biogenic CH4 would have to be >100 µM to isotopically shift the >17 mM abiogenic CH4. The
13C value for the
121 µM DIC is not as isotopically enriched as that reported by Stevens and McKinley for the Columbia River basalt aquifer (68), which suggests that autotrophically produced CH4 would have to be <100 µM.
This ecosystem appears be controlled in a manner altogether different from that for subseafloor sediments (20, 55) and shallow terrestrial sedimentary environments, which are generally electron donor limited and represent the terminal portions of complex food webs which begin with detrital organic carbon and electron acceptors that are ultimately derived from the photosynthetically oxidizing surface world. In sediments, the H2 concentration is controlled by the major terminal electron-accepting process, with methanogenesis and SO42 reduction depleting H2 to 7 to 10 and 1 to 1.5 nM, respectively (42). The D8A H2 concentration is at least 3 orders of magnitude (Table 2) greater than would be predicted for either process. The potential cellular electron flux from H2 is
16 times that of SO42, but only about one-half that of HCO3, taking into account the H+ demands of the methanogenic reaction (Table 5). If the cellular HCO3 flux truly exceeds the H2 flux, then this may explain why the
13C value for the DIC is not greatly isotopically enriched, since it is not limiting. The answer to why the autotrophic methanogens do not deplete the pool of H2 may reside in the free energy for this reaction, which hovers between 63 and 70 kJ/mol (Table 5), close to the energy required for 1 mol of ATP. Significant reduction of the H2 concentration in the presence of elevated CH4 concentrations would reduce the free energy to below the minimum required for ATP production. The resulting H2 concentration greatly exceeds the requirements imposed by the relatively miniscule SO42 flux. The HCO3 concentration is controlled by the high Ca2+ concentration, high pH, high temperature, and solubility of calcite. This limitation on the electron acceptor flux and the minimal free energy requirement for methanogenesis may also account for the substantial excesses in H2 that have been noted for other terrestrial subsurface environments (38, 68). Thus, given the H2 abundance, both methanogenic and SO42 reduction metabolisms can coexist.
The Columbia River basaltic aquifers are similar to the D8A fault water in that they are alkaline (7.9 to 9.9) and anaerobic, but they are also cooler (18 to 20°C), are much less saline (0.1 to 4 mM), contain much higher DIC concentrations (1 to 3 mM), are much less sulfidogenic (0.2 to 30 µM), and have much lower concentrations of CH4 (2 to 209 µM). Rock/water microcosms incubated with a suite of well water samples under excess H2 and CO2 indicated that the balance between methanogens and SO42 reducers was controlled by SO42 concentrations, with methanogens dominating under conditions of SO42 limitation (68). A subsequent quantitative rRNA analysis of the microbial community structure of some of the same wells (methanogenic [DB-11] and sulfidogenic [DC-06]) revealed that the active microbial community in the shallower methanogenic well was dominated by bacteria (64% bacteria versus 3% archaea) and that in the deeper sulfidogenic well was even more so (92% bacteria versus 1.8% archaea) (24). Gram-positive bacteria and
-proteobacteria comprised a small portion of the bacterial component. Overall, the Columbia River basaltic aquifers represent a fresher, shallower, cooler, and microbially more diverse ecosystem than the deep faults intersected by D8A.
The D8A fault water ecosystem is also different in many respects from the 59°C, 70-m-deep groundwater of Lidy Hot Spring (12). The Lidy Spring water has a lower pH (6.9), lower salinity (0.2 mM), and much higher DIC (55 mM) and SO42 (1.5 mM) concentrations than the D8A water. Its CH4 concentration (0.1 mM), H2 concentration (0.013 µM), and DOC concentration (20 µM) are much lower, making it an electron donor-limited environment, unlike D8A. The 16S rRNA gene library was almost entirely dominated (
99%) by methanogens (12), and sulfate reducers were not detected, despite the abundance of SO42.
The D8A results do bear a striking resemblance to those for the Lost City seafloor vent field, which has been interpreted as a microbial ecosystem driven by abiotic CH4 and H2 generation (7). The 60 to 75°C water venting from the Lost City carbonate towers (35) is similar to D8A borehole water in that it has a high pH (9 to 11), low DIC concentrations, high concentrations of abiogenic CH4 (1 to 2 mM), and elevated H2 concentrations (<1 to 15 mM) but differs by having a greater salinity and higher SO42 concentrations (1 to 4 mM). The high-temperature microbial community is dominated by a methanogen, Methanosarcinales, and by thermophilic, sulfate-reducing Firmicutes (35), similar to D8A. Unlike the case for Lost City, however, in D8A the H2 and possibly the SO42 are generated by radiolysis of water in pyrite-rich quartzite (40), and the formation of abiogenic hydrocarbons occurs in the absence of ultramafic rock strata (67).
The identification of ocean vent microorganisms in the boreholes of deep South African Au mines has been reported previously, e.g., Pyrococcus abysii was reported by Takai et al. (70). In this study, the microbial ecosystem of D8A is most similar to that associated with off-axis ocean crust fluids and vents that form when seawater circulates through the ocean crust a few million years after it has moved from the axis and cooled. These thermophilic microorganisms could be inhabitants of the deep ocean selectively colonizing the ocean crust during the infiltration of SO42-rich seawater as part of hydrothermal circulation, or they may be permanent residents of the off-axis ocean crust with its unique thermal alkaline environment that are migrating horizontally as the ocean crust slowly moves away from the central axis.
The quartzite of the Witwatersrand Basin experienced sterilizing temperatures long after the last tectonic/metamorphic episode at 2 Ga. At 90 Ma,
2 kilometers of the overlying Karoo supergroup sediments and volcanics was removed by erosion, and the underlying strata present at 3.2 kmbls began to cool from 120°C, reaching the present-day temperature at about 35 Ma (49). The 4- to 53-Ma He model ages for the fault water overlap with this post-Mesozoic uplift and cooling interval and are consistent with two origins for the microbial inhabitants of the D8A faults. They either existed in the overlying aquifers and selectively colonized the fault zone during meteoric water infiltration, or they have been residents of the thermophilic, alkaline environment of the Witwatersrand subsurface and have slowly migrated downwards through the crust with uplift and cooling. Regardless of the origin of the D8A microbial community, oceanic and continental crusts with distinct thermal histories appear to host geochemically and microbially similar ecosystems. The Lost City hydrothermal field has been advanced as an exciting natural laboratory for the study of ecosystems driven by abiotic CH4 and H2 (7). The functionally analogous system described here significantly expands this concept to include an alternative mechanism for abiotic CH4 and H2 production within metamorphic and sedimentary strata of the continental deep subsurface.
We gratefully acknowledge the support of Gold Fields Ltd., Driefontein Consolidated Mine, the Universities of Witwatersrand and Orange Free State, and Rob Wilson and colleagues of Turgis Technology (Pty.) Ltd. Special thanks to Dawie Nel, the Driefontein senior geologist, and Ben Kotze, the 9 Shaft master shaft sinker. Finally, we acknowledge the insightful comments of the reviewers, which greatly improved our interpretations of the data.
Present address: Florida State University, Department of Oceanography, Tallahassee, FL 32306. ![]()
Present address: Department of Civil and Environmental Engineering, Marquette University, P.O. Box 1881, Milwaukee, WI 53201. ![]()
Present address: Institute of Zoology, National Taiwan University, No. 1, Sec. 4, Roosevelt Road, Taipei 106, Taiwan. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»