High Rates of Sulfate Reduction in a Low-Sulfate Hot Spring Microbial Mat Are Driven by a Low Level of Diversity of Sulfate-Respiring Microorganisms

ABSTRACT The importance of sulfate respiration in the microbial mat found in the low-sulfate thermal outflow of Mushroom Spring in Yellowstone National Park was evaluated using a combination of molecular, microelectrode, and radiotracer studies. Despite very low sulfate concentrations, this mat community was shown to sustain a highly active sulfur cycle. The highest rates of sulfate respiration were measured close to the surface of the mat late in the day when photosynthetic oxygen production ceased and were associated with a Thermodesulfovibrio-like population. Reduced activity at greater depths was correlated with novel populations of sulfate-reducing microorganisms, unrelated to characterized species, and most likely due to both sulfate and carbon limitation.

Cyanobacterial microbial mats are layered communities in which the combination of microbial activity and environmental gradients contributes to a fine vertical stratification of microbial populations and activities. They are self-sustaining communities that support all major biogeochemical cycles and are thought to be analogous to some of the earliest communities on Earth (9,11). These attributes have made microbial mats attractive models for better understanding both the evolution and the ecology of microbial systems supporting the processes that now sustain our biosphere.
Modern mat communities occur primarily in environments in which temperature or salinity suppresses grazing by higher eukaryotes (22,38,45). The mats that develop in low-sulfate, moderately thermophilic (50 to 72°C), moderately alkaline (pH ϳ8), siliceous hot spring effluents in Yellowstone National Park (YNP) are among the best characterized, and analyses of these systems have provided an increasingly complete view of their diversity and ecology (41)(42)(43). Most past studies of these mats have focused on the contribution of cyanobacteria (42), key biogeochemical players through both primary production and the generation of oxygen used by highly active populations of aerobes near the surface.
Although all major anaerobic processes are also recognized to be present in the deeper, anoxic regions of these mats, much less attention has been directed toward resolving their contribution to mat biogeochemistry. In particular, since these mats develop in low-sulfate springs, the role of sulfate-respiring prokaryotes (SRP) has received relatively little attention (13,33,48). However, a low concentration of sulfate can support a highly active population of SRP in communities that sustain populations mediating both the oxidizing and the reducing pathways of the sulfur cycle. This has been observed in certain freshwater lakes (21) and in hot springs supporting the growth of photosynthetic biofilms (15,16). The importance of sulfate respiration in low-sulfate environments is also of considerable significance for understanding the biogeochemistry of early Earth, in which oxidized sulfur species were much less abundant than today (6,20). The present study combined radiometric rate measurements, chemical profiling, and molecular phylogenetic analyses to investigate the role of sulfate respiration in the microbial mat found in Mushroom Spring, one of these low-sulfate thermal spring effluent environments. This study represents the first comprehensive investigation into the changing daily role of sulfate respiration and is the first to link activity measurements to phylogenetic identity in these low-sulfate microbial mats. It revealed a highly active sulfate-reducing community composed of only a few dsrAB genotypes, especially active near the surface during the switch from net photosynthetic to net respiratory conditions at the end of the day.

MATERIALS AND METHODS
Study site and sample collection. During the summers of 2000, 2001, and 2003, in situ spring water pH and temperature in the effluent channel of Mushroom Spring, a thermal spring located in the lower geyser basin of YNP, WY, were measured using a pH 330 combined pH and temperature probe (WTW, Ft. Myers, FL). Samples between 56 and 60°C, a temperature range over which a 2to 4-cm-thick microbial mat was growing on top of the siliceous substrate, were collected. Replicate 1-cm-diameter glass cores were used to sample the mat for nucleic acid extraction, sulfate determinations, and sulfate reduction rate (SRR) measurements. Cores for nucleic acid extraction were immediately frozen on dry ice and returned to the laboratory, where they were stored at Ϫ80°C until extraction.
Sulfate measurements. In 2000, bulk porewater sulfate concentrations ([SO 4 ] values) were determined for 1-cm cores by capillary electrophoresis using a Quanta 2000 capillary ion analyzer (Waters Corp., Milford, MA) equipped with a negative power supply. In 2001 and 2003, the sulfate was analyzed by ion chromatography. In June 2001, a single set of core samples (n ϭ 4) was collected at 09:00 with a 30-ml syringe with the tip removed. Cores were sliced with a sterile razor blade in 2-mm slices for the first 4 mm and into 4-mm slices below that until fully sampled (typically 2 to 3 cm). On 6 and 7 September 2003, diel samples were taken at 20:00, 23:00, 08:35, and 11:41, in parallel with the first four SRR sampling points (see below). Replicate cores (n ϭ 3 to 5) were sampled with a 20-ml syringe with the tip removed and sliced using very fine monofilament fishing line (4 lb) at the following intervals from the surface of the mat: 0 to 2 mm, 2 to 4 mm, 4 to 6 mm, 6 to 8 mm, 8 to 10 mm, and 10 to 15 mm. In both 2001 and 2003, core slices were placed into a preweighed 15-ml centrifuge tube preloaded with 1 ml of 0.1% ZnCl 2 and frozen on dry ice. The samples were returned to the University of Southern Denmark, where the sulfate was analyzed using an ICS-2000 anion exchange chromatograph (Dionex). The eluents were 3.5 mM Na 2 CO 3 and 1 mM NaHCO 3 , and the flow rate was 2 ml min Ϫ1 . The detection limit was 1 M sulfate.
Diel SRR assays. During July 2000, two sets of triplicate core samples were taken during daytime and incubated in the hot spring effluent following vertical injection of 10 l (0.37 MBq, 10 Ci) of Na 2 35 SO 4 (specific activities, 38.8 to 59.2 TBq mmol Ϫ1 ) through the center of the core with a Hamilton syringe. One set was left under ambient-light conditions, while another set was incubated in the dark to control for light effects. In July 2001 and September 2003, triplicate core samples were taken at multiple time points over the diel period. To minimize the effect of reoxidation of sulfide generated (19), 5-min incubation times were used in all experiments. Following incubation, the cores were immediately removed from the spring and sliced, and further sulfate reduction was stopped and the sulfide bound by placing each slice into a 2-ml screw cap tube containing 0.5 ml of 20% zinc acetate solution and freezing it on dry ice. Slicing was performed at an ambient-air temperature by using the monofilament line over the following intervals: 1 mm for the first 4 mm, 2 mm for the next 6 mm, 0.5 cm for the next cm, and 1 cm below that until the core was fully sampled (usually 2 to 3 cm). The entire slicing/killing procedure typically took approximately 5 min; however, the surface slices were taken first and reactions were stopped immediately. SRRs were calculated according to the equations of Fossing and Jørgensen (17), using an incubation time of 5 min. This may result in our calculated SRRs being slight overestimates of specific rates for some of the deeper-sliced samples; however, consistent sampling practices across replicate diel sampling time points and across years allowed us to assess relative diel changes in rates with confidence. Sulfate reduction was confirmed to be due to bacterial activity by using triplicate cores that were incubated with 10 l of 66 mM sodium molybdate injected immediately after the radiolabeled substrate was added. These rates were never higher than 6% of the midday rates for each layer. Additionally, a nonincubation control core was stopped by addition of 0.5 ml of 20% zinc acetate and freezing on dry ice immediately after injection with Na 2 35 SO 4 . These nonincubation control values were subtracted from rates for each layer and were never greater than 10% of the lowest rates measured. In 2000, samples were analyzed for sulfate reduction by using the single-step chromium reduction method (17).  (37) was used for analyzing the SRRs of these samples. Similar rates in replicate cores processed using both methods were found (unpublished observation).
Microelectrode profile measurements. Oxygen microelectrode profile measurements were made in parallel with SRR incubations in 2000 and 2001, using a Clark-style oxygen microsensor with a 10-m tip manufactured by Diamond General (Ann Arbor, MI) connected to a Keithley 485 picoammeter. Vertical profiles were made by sequentially moving the electrode, using a micromanipulator, down into the mat and recording the electrical current. In 2003, both oxygen and sulfide profiles were made using Clark-style microsensors (Unisense, Aarhus, Denmark).
Hydrogen sulfide profile data from the late evening, after oxygen concentrations in the mat had dropped and sulfide was detectable in the mat, were used to calculate specific SRRs by using diffusion coefficients and parabolic function fitting as described in Kühl and Jørgensen for biofilms (26). Rates were calculated as the product of the diffusion coefficient and the coefficient of X 2 for the parabolic function calculated from the H 2 S profile. As reported in the Kühl and Jørgensen study, we multiplied the diffusion coefficient by 80% to account for the reduced diffusion within the mat compared to that in the water environment in which the coefficient was originally measured (26). However, we used a higher sulfide diffusion coefficient of 3.87 ϫ 10 Ϫ5 (35) to account for the higher temperature in the hot spring mat environment.
Light measurements. Irradiance was measured in 2003 using an LI-1000 data logger with an LI-190 quantum sensor (Li-Cor, Lincoln, NE).
Nucleic acid isolation. Total microbial mat DNA was extracted in both 2000 and 2001 from the 0-to 1-mm layer, dominated by cyanobacterial pigments, and from the 1-to 2-mm layer, rich in carotenoid pigments of single-core samples. In 2001, an additional slice below the oxycline between 15 and 18 mm was sampled from the same core. All layers were excised with a sterile razor blade, and extractions were performed using an UltraClean Soil DNA kit (MoBio, Solana Beach, CA), following the manufacturer's instructions.
PCR amplification, cloning, and sequencing. PCR amplifications were performed using modified versions (Table 1) of the DSR1f and DSR4r primers, targeting a large fragment of the dissimilatory sulfite reductase gene (dsrAB) (39). PCR mixtures contained 2 mM MgCl 2 , 0.2 mM of each deoxynucleoside triphosphate, 2 l of 10ϫ PCR buffer, 1.5 U Platinum Taq DNA polymerase (Gibco/Life technologies, Rockville, MD), a DNA template at 0.04 to 0.4 ng/ml, and 10 pmol each of both primers in a total volume of 20 l. Because the optimal annealing temperatures of the degenerate dsrAB primers can vary for different templates, amplifications were carried out with a gradient thermocycler (Ther-moHybaid, Franklin, MA) programmed as follows: an initial denaturation for 5 min at 95°C, followed by 30 cycles at 95°C for 30 s, exposure to a gradient of annealing temperatures from 50°C to 65°C for 30 s, and elongation at 72°C for 1 min. The reaction was completed with a final elongation step at 72°C for 10 min. The best amplified dsrAB gene products were ligated into the pCR4-TOPO cloning vector and transformed into ONE SHOT competent Escherichia coli cells, following the manufacturer's directions (TOPO TA cloning kit for sequencing; Invitrogen/Life Technologies, Carlsbad, CA). Cloned inserts were first screened for proper size (ca. 1,800 bp) by amplification using M13 vector primers and then sequenced with M13F/R vector primers. After initial screening of clones for sequence similarity, full-length sequences from select clones from each major clade were obtained using sequence-specific internal primers designed for this study (Table 1). Sequencing reactions were performed using a DYEnamic ET dye terminator cycle sequencing kit (Amersham Biosciences, Piscataway, NJ) and run on a Megabace 1000 system (Amersham Biosciences) at the University of Washington Center for Marine Biotechnology.
Phylogenetic and statistical analyses. Phylogenetic analyses were performed on alignments of the inferred amino acid sequences of amplified dsrAB DNA sequences. Alignments of dsrAB were performed using the alignment tool of the ARB software package (http://www.arb-home.de) (28) and comparing the sequences to those in our own dsrAB database. Alignments were then refined manually by visual inspection. Regions of uncertain alignment or containing missing data were excluded from analysis. Phylogenetic trees were constructed using neighbor-joining and parsimony methods implemented in PAUP* v. 4.0b10 software (D. L. Swofford, Sinauer Associates, Sunderland, MA). Bootstrap resampling of the neighbor-joining and parsimony trees was performed using 1,000 replicates for each.
Regression analyses and one-and two-way analyses of variance (ANOVA) were performed using StatView v. 5.0 (SAS Institute, Cary, NC). Planned pair-

RESULTS
Hot spring conditions. Physical and chemical conditions in Mushroom Spring and associated effluent channels were consistent among the 2000, 2001, and 2003 sampling dates and consistent with reports dating back to at least 1988 (2). The pH levels in the outflow channels were 8.1 to 8.3, and thick (2-to 4-cm) mats were consistently found between 56 and 60°C. Sulfate levels in the Mushroom Spring microbial mat were very low. During the 2000 sampling period, the bulk porewater [SO 4 ] in the mat was 193 M. In 2001, the value was 184 M in the overlying water, but the mean porewater concentration within specific mat slices was 102 M in the upper 1 mm and decreased rapidly to ϳ10 to 15 M within 1 cm (Fig. 1) (Fig. 2). Pairwise comparisons among depth intervals confirmed that mean porewater [SO 4 ] values were significantly higher in the upper 0-to 2-mm and 2-to 4-mm layers than for any lower-depth interval (Fisher's PLSD; P Ͻ 0.0001 and P Ͻ 0.005, respectively). Mean sulfate values changed significantly over the diel cycle, with the highest levels at midday and the lowest levels at nighttime (P ϭ 0.0104) (Fig. 2). A significant depth-by-time interaction term (P ϭ 0.0038) reflects a less steep decrease in sulfate between the 0-to 2-and 2-to 4-mm intervals during the midday sampling time than at other sampling times. These diel changes may reflect the changes in sulfate consumption due to sulfate respiration (see below).
Diel cycles of sulfate respiration. Significant daytime SRRs were measured in the mat during the initial 2000 sampling trip, revealing low but greater-than-nonincubation-control rates near the surface, increasing below the oxycline (ca. 3 mm) to a peak between 15 and 20 mm and declining below that (Fig.  3A). Incubations under dark conditions resulted in a significant increase in mean specific SRR (7.3 versus 143 nmol cm Ϫ3 day Ϫ1 ; one-way ANOVA; P ϭ 0.0192) for the upper (0-to 4-mm) zone of the mat, whereas at the permanently anoxic depths below 4 mm, mean SRRs were nearly identical in light and dark (125 versus 128 nmol cm Ϫ3 day Ϫ1 ; P ϭ 0.9409). Similar daytime results showing low sulfate respiration activity in the presence of photosynthetic oxygen were detected in subsequent sampling trips in 2001 and 2003, although the depths with peak SRRs in those years were shallower at ca. 6 to 8 mm, with the SRRs decreasing below that depth to near zero (Fig. 3).
In both 2001 and 2003, rate measurements at different depths were performed periodically over the diel period, allowing a more comprehensive view of the dynamics of SRR in this system ( Fig. 1 and 2, respectively). In both years, highly significant changes were detected with two-way ANOVA of SRR with depth (P Ͻ 0.001) and time of day (P Ͻ 0.001). Significant depth-by-time interaction terms (P ϭ 0.0068 and P Ͻ 0.001, respectively) indicate that the depth profiles changed over the course of the day. In general, SRRs were most variable in the upper 4 to 6 mm of the mat and lowest during daytime hours when photosynthetically generated oxygen was at or above saturation, peaking during early evening hours after oxygen was consumed in most of the mat and remaining high a few hours later in complete darkness ( Fig. 1  and 2).
Evening-and nighttime-specific SRR values were as much as 27 to 50 times higher than midday values for a given millimetric interval. In contrast to the midday SRR peak found deep in the mat as described above, rates during the dusk and nighttime points were highest at or just below the surface, decreasing monotonically below that in parallel with declining [SO 4 ] values ( Fig. 1 and 2). Evening and nighttime SRRs averaged over the upper 6 mm of the mat were significantly higher than dawn SRRs for layers below 6 mm did not significantly change over the diel period (P Ͼ 0.05 for both years). A similar daily pattern was found when depth-integrated SRRs (mol cm Ϫ2 day Ϫ1 ) for entire cores were calculated (see inset values in Fig. 1 and 2). In 2001, integrated rates for dusk and nighttime samples were significantly higher than those for the midday samples (Fisher's PLSD; P Ͻ 0.01 for both). In 2003, integrated rates for the first evening samples were significantly higher than those for both the dawn and the midday sampling periods (Fisher's PLSD; P Ͻ 0.05 for both). The nighttime and the second evening samples had higher average integrated rates, but these differences lacked statistical significance. Despite the significant depression of SRRs above the oxycline (ϳ6 mm), 14 to 45% of the integrated sulfate respiration activity occurred in this portion of the mat at midday. Therefore, significant sulfate respiration occurred even in the presence of oxygen.
In 2003, the relationship between SRR and [SO 4 ] was assessed. During the evening transition period, in layers without any oxygen (below 1 mm), there was a significant positive correlation between [SO 4 ] and SRR (R 2 ϭ 0.479, P ϭ 0.005) (panel 1 in Fig. 2), whereas at nighttime, characterized by diminished levels of sulfate at the surface (panel 2 in Fig. 2), no correlation between SRR and [SO 4 ] was found (R 2 ϭ 0.007, P ϭ 0.7124).
Diel changes in microelectrode profiles. Oxygen levels changed predictably throughout the diel period, increasing rapidly at sunrise, reaching supersaturation during the daytime, when high irradiance drives cyanobacterial photosynthesis, and dropping to zero in all but the upper millimeter by evening, when heterotrophic consumption dominates (Fig. 2). The peak oxygen levels measured in 2000 and 2001 during July were ca. 35 to 55% higher than the September 2003 levels, but the depth of oxygen penetration was somewhat deeper in 2003 (Fig. 3).
Hydrogen sulfide levels were measured during the transition period during the second evening period in 2003. H 2 S profiles were very consistent with the corresponding evening patterns of SRR described above; sulfide levels increased to a peak at a few millimeters of depth and then declined below that (Fig. 4). There were Ͼ2-fold increases in sulfide concentrations calculated from profiles made between 19:39 and 20:20. Microelectrode depth profiles showed peak H 2 S concentrations in the mat just below 2 mm of depth at 19:39, closely mirroring the peak SRRs measured during evening sample periods. Peak H 2 S levels shifted downward in the mat during subsequent measurements until they peaked between 3.25 and 3.5 mm at 20:20. Specific SRRs calculated from the sulfide profiles by using the parabolic fitting approach increased from 2.97 to 4.86 mol cm Ϫ3 day Ϫ1 as H 2 S levels increased during the evening. These values represent total rates calculated for the depths over which sulfide increased in the upper portion of the parabolic curve (see the darker shaded area in Fig. 4). The calculated SRR values were 86 to 140% of the highest radioisotopically determined specific SRRs for evening periods for 2003.
Diversity of sulfate-reducing populations in Mushroom Spring. Despite screening of over 50 clones, molecular analyses revealed only a single dsrAB phylotype in the upper two millimetric layers (Fig. 5, clade 1) in direct association with the Synechococcus-Chloroflexus layers that were routinely exposed to high levels of oxygen. Surface sequences were 99% identical at the amino acid level and most closely related to the Thermodesulfovibrio clade, with amino acid identities of ca. 85% with T. yellowstonii and T. islandicus. Additionally, closely related members of this phylogenetic grouping (98% amino acid similarity) were recovered from the anoxic 15-to 18-mm layer (Fig. 5, clade 1). This clade had very strong bootstrap support (100% for each) in both neighbor joining and parsimony analyses (Fig. 5).
A greater diversity of dsrAB sequence types was detected in the anoxic layer, as three other clades were detected in the 15to 18-mm sample (Fig. 5, clades 2 to 4). None of these three clades was closely related to any cultured SRP. One group of Mushroom Spring dsrAB sequences (Fig. 5, clade 3) had wellsupported phylogenetic affiliation (100% of bootstrap replicates with both phylogenetic approaches) with sequences recovered by members of our group during a previous survey of hot springs in YNP where sulfate reduction was detected (16). The other two groups of sequences were not significantly similar to any sequences in the GenBank database (BLAST E values of ϾeϪ39 for both).

DISCUSSION
Three decades of study of laminated mats in low-sulfate thermal springs have provided a better understanding of the microbial diversity of these communities and highlighted the importance of a number of metabolic processes, including ox-ygenic photosynthesis by the unicellular cyanobacterium Synechococcus lividus (31), cross-feeding between cyanobacteria and photoheterotrophic Chloroflexus-type filaments (2), nitrogen fixation (34), methanogenesis (40,47), and acetogenesis, and the importance of anaerobic processes in decomposition and mineralization (1,12,34). Although the contribution of sulfate-reducing bacteria to processes in photosynthetic mats in high-sulfate hypersaline environments is well documented (11,24,25,29,30,32), their contributions to processes in mats that develop in low-sulfate environments are less well studied. However, the study of sulfate respiration in microbial mats in low-sulfate environments is especially relevant to understanding the evolution of these communities, since isotopic data indicate that sulfate levels, even in marine environments, were likely to have been Ͻ200 M on the Archean Earth (20).
The presence of sulfate-respiring microorganisms in mats that develop in low-sulfate geothermal outflows in YNP was first shown by chemical analyses (13), demonstrating that over a 48-h period sulfide formation was associated with declining sulfate concentrations with depth. Two additional reports, using radiotracers to survey the distribution of sulfate reduction, confirmed activity in Octopus Spring (48) and Mushroom Spring (33) mats, showing rates comparable to our daytime values (91 to 483 nmol ml Ϫ1 day Ϫ1 ). Our analyses, examining depth-specific and depth-integrated changes in SRR over a diel period, now provide a more complete accounting of this process and its association with the populations likely contributing to the observed activities.
Although sulfate reduction was measurable at all depth intervals examined, the highest rates near the surface were observed shortly after the collapse of the oxygen chemocline at dusk or in immediate response to shading. This pattern differs from that observed in mats that develop in hypersaline habitats, where elevated daytime SRRs (7,18,36) have been attributed to increased temperature (7) and photosynthetic productivity (5,8,18). The SRR patterns in Mushroom Spring were more similar to the nightly increase in nitrogen fixation, another oxygen-sensitive process, observed in coastal microbial mat populations (3). Despite these differences from hypersaline mats, the Mushroom Spring data may also reflect an influence of irradiance-driven photosynthetic productivity on sulfate reduction. Previous studies of hypersaline mats have shown that supplementation with photoexcretion products such as glycolate stimulate sulfate respiration activity (18). Chloroflexus filaments in Mushroom Spring have been reported to consume glycolate excreted by Synechococcus (2). The dramatic increase in specific SRRs at the Mushroom Spring mat surface may simply be related to the decrease of oxygen that was inhibiting sulfate respiration or reoxidizing sulfide (see below) during the day; however, the subsequent decline of depth-integrated rates from the evening to a minimum at dawn may be due to the gradual depletion of photosynthetically derived carbon in the mat. The connection between SRR and photosynthetic productivity may also have contributed to the up-to-fourfold decrease in depth-integrated SRRs between the July 2001 and September 2003 sampling dates, perhaps reflecting seasonal changes in day length and peak irradiance between those two study periods. This is supported by microelectrode profile data that showed peak photosynthetic oxygen levels near 900 M in September compared to levels around 1,400 M in 2001. Although daytime SRRs in both 2001 and 2003 were comparable, it may be that the size of the evening SRR peak is related instead to the length of the day and the potentially larger accumulated pool of carbon formed over the whole day.
Our studies also point to the possible presence of an active sulfur cycle near the surface of the mat during daytime. Although appreciable sulfate reduction at the surface was not observed following short 5-minute incubations in the light, the almost instantaneous 20-fold increases in rates in response to shading observed in the 2000 dark control experiment is suggestive of an active sulfur cycle in this region. Reoxidation of sulfide, via either aerobic oxidation or anoxic photosynthesis, is inhibited when both oxygenic and anoxygenic photosynthesis is attenuated by shading (31). These data point to the existence of a resident population of SRP active near the surface, possibly respiring sulfate in the presence of oxygen or, alternatively, rapidly initiating sulfate respiration when the oxygen level drops following attenuation of light. Thus, our measurements may underestimate true daytime SRRs at the surface.
Overall, the SRRs that we detected in the current study are comparable to rates detected in a cyanobacterial mat in a high-sulfate hot spring (Bath Lake; 6 to 8 mM) in Mammoth Basin, YNP (44). Notably, the high evening-specific SRRs that we measured in 2001 (ca. 14 mol cm Ϫ3 day Ϫ1 ) are comparable to the highest rates of sulfate reduction reported by Habicht and Canfield (19). In that study, SRRs of 23 to 37 mol cm Ϫ3 day Ϫ1 were detected in microbial mat and sulfuretum sediment samples incubated in bags under artificially elevated temperatures (30°C). The high in situ rates in Mushroom Spring may be partly due to enhanced reaction rates at high temperatures (56 to 58°C), although past reports have shown similar optimal specific SRRs for cultures adapted to mesophilic or thermophilic temperatures (6).
The radiometric SRR measurements confirmed the sulfide profile data set. Past authors have attempted to measure sulfide in high-temperature springs but have encountered difficulty in calibrating voltametric electrodes (31). This report, using Clark-style amperometric sensors, represents the first measurement of sulfide in hot springs. Our results were quite typical compared with profiles collected from sediments and an aerobic biofilm (26,27), except that the decrease in sulfide with depth that we saw does not typically occur in higher-sulfate environments. Calculated specific SRRs based on the measured sulfide profiles were very consistent with radiometric determinations of specific rates, indicating that sulfide profiles may appropriately be used to estimate SRRs in high-temper- The microelectrode profiles and radiometric data both showed similar patterns of declining sulfide production in the lower portion of the mat. This vertical pattern of declining SRRs has been repeatedly observed in marine sedimentary habitats (14,23,46). Those studies have typically suggested that carbon limitation is responsible for the decline of activity with depth. Because of the vanishingly low [SO 4 ] in Mushroom Spring, it may be a combined effect of carbon and/or sulfate limitation that is responsible for the decline in this ecosystem.
Analysis of dsrAB sequence types revealed four major clades of sulfate-reducing phylotypes that were differentially distributed with depth in the mat. Three were recovered only from the exclusively anoxic zone of the mat, consistent with a tradi-tional understanding of the habitat preference of SRP. One of these groups (clade 3) showed significant homology with sequences isolated from other thermal habitats in YNP. Interestingly, temperatures among these sites vary and neither Obsidian Pool (89°C; pH 6.9; [SO 4 ], 0.54 mM) nor Black Sediment Pool (69°C; pH 6.6; [SO 4 ], 0.6 mM), where those related sequences were obtained, contain photosynthetic microbial mats; however, like Mushroom Spring, they both have relatively low sulfate concentrations. This suggests that this dsrAB clade may be broadly distributed among low-sulfate thermal habitats.
In contrast, Thermodesulfovibrio-like dsrAB genotypes were recovered from both oxic and anoxic depths. Although SRP in Octopus Spring have been reported (43), including a Deltaproteobacterium-and a Thermodesulfobacterium-like sequence, FIG. 5. Unrooted dendrogram depicting relationships among dsrAB gene sequences for clones obtained from mat sections from within (0 to 1 mm and 1 to 2 mm) (in bold) and below (15 to 18 mm) (in italics) the midday oxycline and closely related sequences. The tree was constructed using an alignment of 481 amino acid positions (gaps and ambiguous residues were excluded using a filter in ARB). Symbols at branches represent nodes where bootstrap confidence values of Ͼ50% (Ⅺ), Ͼ75% (E), and Ͼ90% (F) were observed for bootstrapped neighbor-joining and parsimony trees and expressed as percentages of 1,000 replicates for both. Nodes without values were unsupported by the majority of tree reconstructions in one or both methods. The total number of replicate clones in each clade or subclade in each of the clone libraries is shown in parentheses.
this is the first observation of a Thermodesulfovibrio-like population in a hot spring microbial mat. A significant contribution of these Thermodesulfovibrio-like species to sulfate reduction was also supported by the recovery of this sequence type in total bacterial 16S rRNA gene clone libraries (2 to 4% of clones) generated from DNA extracted from the surface of the mat (unpublished observations). Thus, this possibly more oxygen-tolerant population is likely responsible for the immediate production of sulfide observed near the surface following shading of the mat. Variable depth distributions of different SRP phylotypes have also been observed in other stratified systems, including microbial mats (29,32). A trend observed in both our study and one of those hypersaline mat studies was that some dsrAB populations were uniformly distributed with depth, while others were more localized (29). However, analyses of dsrAB phylotypes in the hypersaline mats of Solar Lake (Sinai, Egypt) and Exportadora de Sal Saltworks (Guerrero Negro, Mexico) (30, 32; J. G. Dillon, S. Miller, B. Bebout, N. Pinel, and D. A. Stahl., unpublished data) have revealed a diversity of sequence types significantly greater than that in the samples recovered from Mushroom Spring, consistent with the view that the conditions in hot springs restrict the overall diversity of biological communities (4,43). In particular, the low sulfate availability may specifically restrict the diversity of SRP in this spring.
Our findings demonstrate the importance of sulfate respiration as a terminal process in thermophilic microbial mats, despite fluctuating oxygen levels and very low [SO 4 ] values. Dramatic changes in sulfate reduction between day and night, especially during the evening transition, emphasize the dynamic nature of photosynthetic microbial mats and the limitations of basing conclusions on single sampling points. Other studies of sulfate-rich saline habitats have detected significant SRRs even within the hyperoxygenated surface layers of photosynthetic mats during daytime (7). Although we observed a depression of rates during daytime, as was true in that study, we have detected some of the highest in situ sulfate respiration rates ever reported. These results, along with those of recent studies of pure cultures (10), have challenged the traditional view that all SRP are obligate anaerobes and suggest that further work for improving the understanding of the complex biogeochemical interplay of photosynthesis and sulfate respiration within the context of microbial mats is needed. This should include work investigating the potential transfer of photosynthetic carbon to SRP over the diel period. Future studies, including targeted efforts (cultivation based and cultivation independent) to study potentially oxygen-tolerant SRP, such as the Thermodesulfovibrio-like group identified in this study, are also needed to confirm the role of this group in the Mushroom Spring mat. More generally, an expanded study of the role of sulfate respiration in other low-sulfate thermal habitats, both with and without microbial mats, is warranted to determine the overall importance of the sulfur cycle in these environments.