| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

Department of Earth and Planetary Sciences, Harvard University, 20 Oxford St., Cambridge, Massachusetts 02138,1 Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, California 91125,2 Department of Biological Sciences, University of Southern California, Los Angeles, California 900893
Received 6 May 2007/ Accepted 4 December 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The bacterial genera most commonly identified as hydrocarbon degraders include Pseudomonas, Alcanivorax, Marinobacter, and Cycloclasticus (gammaproteobacteria) (5, 14, 19, 74); Sphingomonas (alphaproteobacteria) (18); gram-positive bacterial genera such as Staphylococcus and Geobacillus (43, 76); and deltaproteobacterial genera such as Desulfobacter, Desulfococcus, Desulfosarcina, and Desulfatibacillum (1, 12, 24, 51, 63). However, because many of the methods for measuring bacterial degradation of petroleum rely on enrichment cultures or isolates (58, 68, 72), the rates at which hydrocarbons are consumed in situ can be difficult to determine. Studies of the degradation of model hydrocarbons by species in pure or enrichment cultures cannot reflect the rates of conversion of oil into biomass by complex microbial communities in their native environments. Since petroleum is a mixture of many compounds with different degradation potentials, estimates based on the uptake of radiolabeled substrates (usually 14C-labeled n-alkanes) (69) likely do not represent the utilization rate for the total oil. This is especially evident in gas chromatograms of contaminated sediments, where typically an "unresolved complex mixture" (UCM) is a dominant feature in addition to n-alkanes (54). Newer methods of stable-isotope labeling (8) or probing (52) combined with molecular analyses of community diversity by PCR-denaturing gradient gel electrophoresis (35, 56) are promising, but they can require long incubation times and the results are difficult to convert to estimates of oil utilization rates. These methods also are subject to the challenges described above, namely, selection of a single molecule as a proxy for a complex, heterogeneous mixture.
To quantify the microbial incorporation of oil in coastal salt marsh habitats, isotope labeling experiments may not be necessary. It is possible to use whole oil and thus avoid many of the concerns raised above. Variations in the natural abundance of the stable isotopes of carbon, 12C and 13C, can be used to trace the carbon sources used by heterotrophic organisms (7, 13). These ratios are expressed as values of
13C:
![]() |
13C value of approximately –13
(48). The values of
13C measured for heterotrophic bacteria living on salt marsh detritus are similar to the values for Spartina, as the marsh grass is the primary producer in these ecosystems. In contrast, fossil fuels, which are derived mainly from marine products, have a lower natural abundance of 13C, resulting in more negative values of
13C. Although petroleum mixtures contain a wide variety of hydrocarbon structures, the
13C value of most oils is near –27
(55). In this work, our goal was to utilize the difference between the natural isotopic compositions of these substrates (13C-enriched Spartina and 13C-depleted petroleum) to monitor uptake of oil-derived carbon into salt marsh bacteria. We studied whole-community uptake of hydrocarbons using natural-abundance 13C while simultaneously investigating the microbes likely to be responsible for the uptake. The 16S rRNA molecule (rRNA) is a logical "biomarker" that can link isotopic composition and taxonomy (39, 40, 46). Specifically, we combined (i) phylogenetic approaches to obtain sequences of 16S rRNA genes with (ii) carbon isotopic analyses of native rRNA molecules. Since rRNA in actively growing bacteria can account for 20% of the cellular mass (dry weight), it is possible to purify sufficient quantities of total bacterial or group-specific rRNA molecules to subject the material directly to isotopic analysis.
To this end, we used a moving-wire interface that allows 13C analysis of complex, nonvolatile organic samples, such as nucleic acids (9). Recent improvements in the sensitivity of the device, now termed "spooling-wire microcombustion" (SWiM) (15, 60), made this work feasible by providing sample requirements compatible with RNA purification protocols. Typical elemental analyzers require several orders of magnitude more material. The rRNA-SWiM approach is likely to have many future applications in investigations of carbon metabolism in heterogeneous microbial communities.
In September 1969, the oil barge Florida spilled 700,000 liters of diesel fuel (no. 2 fuel oil) in Buzzards Bay near West Falmouth, MA (22). This event has become one of the most well-studied oil spills in history. Thirty years after the spill, sediments in nearby Wild Harbor contain remains of degraded oil. The abundance and composition of the residue, known as a UCM, changed only slightly between 1973 and 2000, and the residue remains confined to a narrow 2- to 3-cm horizon now located approximately 10 cm deep in the marsh sediments (54). In the current work, Wild Harbor served as a control. Specifically, we expected to detect no petroleum-related signal in the
13C values of rRNA, since active biodegradation is now undetectable in the contaminated sediment (62). We then used natural-abundance 13C-rRNA-SWiM measurement to quantify the microbial assimilation of hydrocarbons from a modern, simulated spill in a similar environment. We selected nearby Woodneck Marsh, an uncontaminated site, for the simulated spill. This location is a tidally flushed salt marsh similar to Wild Harbor and is also adjacent to Sippewissett Marsh, where community microbial diversity and metabolism have been studied extensively (21, 25, 61). Here we describe both the microbial diversity enriched by addition of no. 2 fuel oil and the use of carbon isotopic signatures of rRNA to quantify incorporation of oil-derived carbon into microbial biomass.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Woodneck Marsh bulk samples.
Bulk samples of Spartina and wild rose were collected from Woodneck Beach Salt Marsh, Falmouth, MA. In addition, bulk biomass was skimmed from the surface of a shallow tidal pool, and a mussel was collected from the same location. These samples were frozen in N2(l) and homogenized with a mortar and pestle. A portion of each bulk sample was prepared for analysis of
13C of total organic carbon using standard sealed-tube combustion procedures (44). Values of
13C were determined using purified CO2 by dual-inlet isotope ratio mass spectrometry (VG Prism).
Woodneck Marsh incubation.
Twelve shallow cores (diameter, 14 cm; depth, 5 cm) were harvested in midsummer and immediately incubated in sterile (baked in air at 450°C for 8 h) Pyrex crystallizing dishes for 0, 1, or 2 weeks under ambient temperature and sunlight conditions; 8 of these cores were incubated with no. 2 fuel oil (Mobil diesel), and 4 were incubated without fuel oil. An excess of no. 2 fuel oil was supplied as an irregular layer of oil-coated sand on top of the eight inoculated samples. The sand was used to keep the oil from forming a continuous film and thus to allow gas exchange. Oil-free zones were observed throughout the experiment, and some areas of the sediment surface remained exposed. To simulate tidal flushing and prevent desiccation,
20 ml of sterile-filtered marsh water was added to the preparations daily. Four of the samples were harvested after 1 week, and the other four were harvested after 2 weeks. The samples were frozen at –80°C for extraction of hydrocarbons and nucleic acids.
Purification of bacterial rRNA. (i) RNA extraction.
Total RNA was extracted using methods described previously (46). Some samples were lysed in 50 mM sodium acetate (pH 5.2) containing 0.5% sodium dodecyl sulfate and 0.1 mg/ml proteinase K (Life Technologies), followed by extraction with phenol-chloroform-isoamyl alcohol (25:24:1) and 100% chloroform, and some samples were extracted with RNAWiz (Ambion) used according to the manufacturer's instructions. Nucleic acids were precipitated with isopropanol and washed three times with 70% ethanol, and the pellets were resuspended in Nanopure H2O. Visual inspection of electrophoresis gels confirmed that the total extracts were enriched in RNA but that some samples also contained small amounts of DNA. Humic acids were removed effectively by the RNAWiz method in particular, and significant A320 generally was not observed for these samples.
(ii) Magnetic bead capture of small-subunit rRNA.
Small-subunit rRNA (16S and 23S rRNA) of bacteria was captured from the total RNA extracts using an Ambion MICROBExpress bacterial mRNA enrichment kit, with minor modifications to the manufacturer's protocol. The MICROBExpress kit contains a probe cocktail designed to bind, capture, and remove bacterial 16S and 23S rRNA molecules to enrich the mRNA in the remaining solution (compatible genera are listed at http://www.ambion.com/techlib/misc/microbe.html). Here, the kit was employed to recover the bacterial rRNA fraction from mixed environmental RNA, and the uncaptured (mRNA) fraction was discarded. The total RNA, capture probe cocktail, and binding solution were hybridized and combined with the included paramagnetic beads, as instructed. After capture of the bead-probe-rRNA complex using a magnetic particle separator as described previously (46), the complex was washed and eluted in Nanopure H2O at 90°C for 60 s. Recovered rRNA was precipitated in 500 µl of isopropanol plus 10 µl of 5 M NaCl, centrifuged (10,000 x g, 15 min), and washed with 70% ethanol.
To remove residual carbon from buffers used in the capture procedure, samples were cleaned using an Ambion MEGAclear RNA purification kit. RNA recovered from the MEGAclear elution step was precipitated in 500 µl of isopropanol plus 10 µl of 5 M NaCl, centrifuged, and washed as described above. Recovered pellets were redissolved in 400 ml of Nanopure H2O and reprecipitated using high-performance liquid chromatography-grade isopropanol (Fisher). The final pellets were washed three times in 70% ethanol, dried at 60C, and redissolved in 10 µl of Nanopure H2O for isotopic analysis. Method blanks were prepared by processing aliquots of Nanopure water that were the same size using the entire hybridization, capture, precipitation, and clean-up procedure.
Isotopic analysis using a SWiM interface.
The SWiM system was described by Sessions et al. (60) and is a modification of the original device built by Brand and Dobberstein (9). Briefly, samples in
1-µl droplets of water are applied to a wire (0.25-mm nickel wire moving at a rate of 0.8 cm/s), which passes sequentially through a cleaning oven (950°C), a drying oven (120°C), and then a combustion oven (750°C) containing wireform CuO. The atmosphere is excluded from the combustion furnace by positive pressure of the helium carrier gas. A portion of the combustion gases flows through a countercurrent Nafion membrane to remove H2O, through an open split, and then to a ThermoFinnigan 252 isotope ratio mass spectrometer. CO2 reference gas for calibration of isotopic ratios was supplied via a custom-built interface similar to the Conflo III device, and samples were analyzed as sets of five to seven injections spaced at 45-s intervals. The carbon content and
13C values were calculated from the peak areas of m/z 44 and 45 ion chromatograms. The precision of
13C values attainable with this system is roughly 1
for samples containing as little as 10 ng C and 0.2
for samples containing 100 ng C (60).
Analytical blanks were estimated by repeatedly measuring 1-µl aliquots of sample-free extracts (described above). They typically contained <10% of the typical carbon concentration in processed RNA samples and had
13C values of –23 to –24
. The contributions to the overall analytical blank arising from the SWiM system are estimated to be <0.2 ng C (55), with the balance of C presumably derived from RNA extraction procedures. Blank contributions were subtracted from all reported sample
13C values, and the error was propagated to the stated uncertainties using equations given by Hayes (26).
Analysis of hydrocarbons.
Aliquots of homogenized material were extracted by vortexing wet sediment in a 1:1 mixture of CH2Cl2 and H2O in 50-ml Teflon tubes. The tubes were centrifuged at 5,000 x g, and the CH2Cl2 layer was moved to a separatory funnel, where it was washed with H2O. The total organic extract was passed over anhydrous Na2SO4, reacted with HCl-activated Cu to remove elemental sulfur, and fractionated over SiO2 gel in hexane-CH2Cl2 (95:5). Hydrocarbons were analyzed by gas chromatography-mass spectrometry using an Agilent 6890N gas chromatograph interfaced to a 5973 MSD equipped with an HP-5MS column (30 m by 0.25 mm [inside diameter]; film thickness, 0.25 µm).The gas chromatograph temperature program was 65°C for 2 min, increasing to 130°C (0 min) at a rate of 20°C/min, to 280°C (0 min) at a rate of 6°C/min, and finally to 320°C (25 min) at a rate of 3°C/min.
Phylogenetic characterization.
Total DNA was extracted using an UltraClean Mega soil DNA kit from MoBio Labs and was concentrated and purified using a Qiagen QIAquick kit. Partial 16S rRNA genes were amplified from Woodneck Marsh samples using universal bacterial primers (0.4 µM each) 27F (5'-AGAGTTTGATCCTGGCTCAG-3') and 1100R (5'-AGGGTTGCGCTCGTTG-3'), a total DNA template (20 ng), and Sigma JumpStart RedTaq Ready Mix for PCR. The thermal cycling conditions were 2 min at 94°C, followed by 30 cycles of 1 min at 94°C, 1 min at 55°C, and 1.5 min at 72°C and then a final extension step of 10 min at 72°C. Duplicates were run for each sample, and 2.5-µl portions of the duplicates were combined in one reaction mixture and amplified for an additional five cycles using a modified "low bias and reconditioning" protocol (67). PCR products were visualized on a 1.0% agarose gel using SYBR green. All reactions were performed using 25-µl (final volume) mixtures. Negative controls (reagents, lab water) were included with every PCR.
PCR products from all samples were cloned using an Invitrogen TOPO TA cloning kit for sequencing with E. coli One Shot TOP10 chemically competent cells. Clones were incubated overnight at 37°C on plates containing LB medium with ampicillin (50 µg/ml), inoculated into 5 ml of liquid medium, and incubated for 24 h at 37°C. Plasmids were purified using a Qiagen QIAprep spin miniprep kit. Sequencing was done at the Dana-Farber/Harvard DNA Resource Core (http://dnaseq.med.harvard.edu/).
Preliminary 16S rRNA sequence alignments were obtained using BioEdit (http://www.mbio.ncsu.edu/BioEdit/bioedit.html) and then were manually inspected for realignment and trimmed to
700 bp. A check for chimeric sequences suggested that there were no chimeras (Bellerophon server) (30). Maximum likelihood trees were created from the clone sequences using PHYML (20; http://atgc.lirmm.fr/phyml/). The substitution model was HKY, and 100 bootstrap replicates were calculated. Sequences from related species and from similar environmental clones were downloaded from the National Center for Biotechnology Information website (http://www.ncbi.nlm.nih.gov/) and included in the alignments and trees. Aquifex aeolicus served as the outgroup.
Rarefaction curves were calculated using distance-based operational taxonomic unit (OTU) and richness determination (DOTUR) (59). Distance matrices for the DNA sequences of the clones were calculated using default parameters of dnadist in PHYLIP (http://evolution.genetics.washington.edu/phylip.html). Sequences were then assigned to OTUs based on the unweighted-pair group method using average linkages (average neighbor) clustering algorithm implemented in DOTUR with default parameters for precision (0.01) and bootstrapping (1,000).
| RESULTS |
|---|
|
|
|---|
|
The values of
13C for the total RNA extracts ranged from –15.5 to –14.3
(Table 1). The differences between the values for extracts obtained in the region of the UCM and the values for the horizons above and below the residue of the spill were insignificant (Fig. 1c). The values of
13C for the captured fraction of bacterial rRNA were systematically more negative by
1
than the values for the total RNA. Three of the four values were in the range from –16.1 to –15.7
, while one sample from within the oil spill horizon (9 to 10 cm) had a slightly more positive
13C value than the other samples (–14.9
). The difference is not significant given the analytical uncertainties, and the direction of the change is opposite what would be expected if the bacteria in this horizon were using oil as a substrate (the
13C for oil is approximately –27
).
|
(i) Composition and degradation of oil.
The control samples contained no detectable hydrocarbon-rich UCM and therefore had an insignificant petroleum component. All of the oil-inoculated samples still contained readily detectable petroleum (Fig. 2a). The chromatograms included a UCM similar to that found in the deeper, Wild Harbor sediments. As a result of the quantity of fresh oil added, n-alkanes were also visible in the chromatograms (Fig. 2a). We examined the composition of this oil before and after incubation by using gas chromatography-mass spectrometry. The net uptake of oil represented a very small fraction of the total amount added, and the composition remained unchanged; we thus presumed that the isotopic composition of the residual oil did not differ from the original composition.
|
The values of
13C for the total RNA extracts ranged from –11.5 to –13.0
(Table 1). These values also decreased moderately over the course of the experiment (Fig. 2c), and the difference,
1.5
, was greater than the isotopic heterogeneity among the replicate samples (the sample standard deviations were ±0.2 to 0.4
). All treatment samples displayed isotopic depletion relative to the control (no-oil) samples.
(iii) Bacterial rRNA.
The values of
13C for the captured fraction of bacterial rRNA ranged from –10.3 to –14.7
. These values reflect a 4.4
difference between the control samples and the oil incubation samples at 2 weeks (Fig. 2c). The results are consistent with the incorporation of oil-derived carbon into bacterial biomass. The 4.4
decrease in the value of
13C of captured bacterial rRNA could be used to estimate the fraction of new bacterial biomass (fnew) in the samples that was oil derived, either directly (heterotrophically) or indirectly (autotrophic incorporation of respired CO2).
![]() |
± 0.6
, an assumed composition for oil of –27.0
± 0.5
, and the final value of –14.7
± 0.5
, the isotopic mass balance suggested that at 2 weeks, 26% ± 4% of the carbon in the captured bacterial rRNA originated from the oil.
Woodneck Marsh microbial community.
Partial gene sequences of 16S rRNA genes were obtained from three of the 2-week incubation samples and from one of the control samples (also after 2 weeks of incubation). Forty-six clones were sequenced from control sample SM-10, and 138 clones were sequenced from samples SM-3, SM-6, and SM-9 (treatments). The clones sequenced from the control sample included Alpha-, Delta-, and Gammaproteobacteria at the highest frequencies (Fig. 3a). In contrast, the treatment samples contained primarily Gammaproteobacteria, Firmicutes, Spirochaetes, and other or unknown organisms, while the relative number of Deltaproteobacteria clones decreased (Fig. 3a). The number of unique OTUs (the OTU threshold was set at 3% divergence) detected at similar sampling levels remained approximately the same, regardless of the presence of oil. Calculations using DOTUR (58) yielded 94 unique OTUs among the 138 sequences in the treatment samples (Fig. 3b) and 33 OTUs among the 46 sequences in the control sample (Fig. 3b). These numbers should be interpreted with caution, however, as the numbers of sequences examined were not the same for the two samples (46 sequences versus 138 sequences) and an estimate of community diversity obtained using the set of only 46 control sequences, in particular, was not likely to be very robust. The individual treatment samples had a diversity similar to that of the control sample; SM-3 contained 36 OTUs, SM-6 contained 43 OTUs, and SM-9 contained 35 OTUs. These numbers are consistent with similar levels of total diversity in all samples.
|
|
No Firmicutes were identified in the control sample. However, 18% of the clones detected in the treatment samples were related to Firmicutes. The species detected were putative anaerobes closely related to Clostridium and Alkaliphilus (Fig. 4b) and were similar to other clones from gas hydrates, oil wells, and anaerobic hydrocarbon enrichments (53, 75).
Many of the new Deltaproteobacteria in the treatment library were similar to clones from hydrocarbon-contaminated systems. These Deltaproteobacteria sequences, however, were more evenly distributed among both the control and treatment samples (Fig. 4c). In both types of samples, clones were found that were related to alkane-degrading species, such as "Candidatus Desulfococcus oleovorans" Hxd3 (1, 63). Clones whose sequences were 100% identical to the sequence of a clone from an oil-degrading hydrocarbon mat from Etang de Berre (29) were also found in both the treatment and control samples; this group also clustered more closely with Desulfococcus spp. than with other sulfate-reducing Deltaproteobacteria (Fig. 4c) The preponderance of the remaining clones appeared to be related to sulfate-reducing species of Deltaproteobacteria and included Desulfovibrio and Desulfotalea spp. No clones closely related to Geobacter or Pelobacter spp. were detected.
Many of the clones in the other major categories (Fig. 3a) also were closely related to clones from hydrocarbon-rich environments, consistent with the overall observation that the experimental clone library reflected many species that have previously been shown to be involved, directly or indirectly, in the metabolism of hydrocarbons (29, 64, 66, 70).
Additional measurements.
The SWiM device was used to obtain
13C data for RNA extracts of pure cultures for method development (Table 2). Concurrent analyses of nucleotide standards, bulk RNA standards, and total nucleic acid extracts produced
13C results that had an analytical precision (based on three to five injections) of around ± 0.2
and that were accurate to within ±1.0
for samples containing 100 nmol of carbon (60). Values of
13C measured for RNA using the SWiM device were consistently 1 to 2
enriched in 13C relative to the biomass from which they were derived (shown for E. coli in Table 2; also measured for Ferroplasma acidarmanus [data not shown]). The results agree with previous data which established that the 13C content of RNA is influenced by the ribosugar component (7), because sugars tend to be enriched in 13C relative to total biomass.
|
| DISCUSSION |
|---|
|
|
|---|
13C value for total RNA in Wild Harbor, which is 3
less isotopically (approximately –14.5
) than the value for untreated total RNA from the surface community in Woodneck Marsh (–11.5
). This effect could be the result of differences in the carbon substrates and assimilation pathways used by the microbial communities in the two locations. In the Woodneck Marsh sediments, both the phylogenetic distribution and the 13C content of bacterial rRNA were affected by the addition of fresh fuel oil. The data demonstrate that the salt marsh microbial community responded rapidly to the introduction of hydrocarbon contaminants. The observation that this carbon entered the active metabolic pool—as detected using the labile 16S and 23S rRNA molecules—is consistent with the incorporation of carbon from petroleum into microbial biomass. Biodegradation of hydrocarbons proceeds first through the more labile n-alkanes to more recalcitrant molecules, such as aromatic hydrocarbons and branched alkanes (54, 58, 69). The duration of the experiment captured only the initial degradation of the labile fraction of oil, and the experiment was not sensitive enough to detect differential utilization of structurally or isotopically distinct components; a longer experiment is needed to determine changes in the microbial community after the depletion of n-alkanes.
However, in this short experiment we estimated that after 2 weeks, up to
26% of the bacterial carbon was derived from the oil. Stable isotope labeling and analysis of captured rRNA have been used previously to assess microbial uptake of a variety of "tracer" carbon compounds (40). Other studies have examined the catabolic conversion of hydrocarbon substrates into CO2 (2), while more recent approaches have included isotope labeling methods (DNA stable isotope probing or lipids) to detect the flow of oil-derived carbon into active microbes (8, 41, 47). Many of these studies utilized labeled substrates and/or examined only a single compound. To our knowledge, the study presented here is the first study to use the natural abundance of 13C in rRNA to trace the microbial incorporation of carbon from total oil as a substrate. Conversely, we could not detect if some components of the oil were preferentially consumed relative to other components, other than to note that the total composition of the mixture did not appear to change with time. Regardless, in general, this experiment expanded the potential for using non-isotope-labeled nucleic acids as tracers for natural biogeochemical processes. In the case of petroleum biodegradation, this approach may be particularly well suited.
Oil spills affect sediment-associated microbial consortia by stimulating the growth of species involved directly and indirectly in metabolism of the complex hydrocarbon mixtures (28, 32, 57, 75). In the Woodneck Marsh incubations we observed differences in the community composition between the controls and the treatment experiments, both within groups (Gammaproteobacteria) and with the appearance of new groups (Firmicutes). The frequency of our clone phylogenies shifted from primarily Alpha-, Gamma-, and Deltaproteobacteria in the control samples to mainly Gammaproteobacteria, Firmicutes, and Spirochaetes in the oil-treated samples. Sulfate-reducing Deltaproteobacteria were present in both types of samples, but the relative number of clones decreased and the diversities were comparable in the control and treatment samples. Also, we did not find examples of frequently isolated oil-degrading genera, such as Azoarcus, Cycloclasticus, or Alcanivorax. Alcanivorax and Cycloclasticus are two of the best-known genera of marine hydrocarbon degraders (14, 33, 34, 56). The failure to detect clones related to these genera may have been due to environmental differences between the Woodneck Marsh salt marsh sediments and other, more open, oceanic sites or could have been due to the relatively short incubation time of the experiment (14 days). The overall biodiversity did not decrease, consistent with some previous reports for hydrocarbon-contaminated marine systems (29, 56).
Many genera of Gammaproteobacteria degrade hydrocarbons in pure or enrichment cultures, and clones associated with Gammaproteobacteria commonly are detected in environmental incubations and in hydrocarbon-contaminated sediment samples (16, 29, 68, 74, 75). Our enrichment of a set of clones with high similarity to Idiomarina spp. (Fig. 4a) also is consistent with the previous detection of this genus in association with oil-rich microbial mats and sediments, (29, 36, 57), although the role that Idiomarina plays in the degradation of hydrocarbons is not known.
The sediment microbial community changed visibly when the oil was added, and gradually white filaments appeared to replace pink and purple filaments. This observation is consistent with the results of our clone library sampling analysis; decreased numbers of Alphaproteobacteria clone representatives and increased numbers of clones that are related to sulfur oxidizers (among the total clones of Gammaproteobacteria) were detected in the library generated for the treatment samples. Many of the new clones observed in the oil-treated samples were associated with Thiomicrospira (73) and other sulfur-oxidizing environmental clones (49, 65). Local anoxia also was indicated by the results for the Woodneck Marsh treatment samples. Microbial consortia associated with hydrocarbon degradation in marine sediments frequently contain sulfate reducers (Deltaproteobacteria) (49, 57, 66), and sulfate-reducing Deltaproteobacteria clones were detected throughout the experiment. Microaerophilic or anoxic conditions also were suggested by the presence of clones affiliated with the genus Clostridium, which were absent from the controls.
While we found that the phylogenetic distributions of clones were different for the treatment and control samples, our results do not quantitatively reflect the population changes induced by the mock oil spill. Biases can be introduced by differential copy numbers of 16S rRNA genes or through artifacts of extraction or PCR (10, 17, 38, 50). The results described here thus provide an overview of differences between the oil-treated and control samples. Future work should combine quantitative PCR with the greater phylogenetic specificity of the 13C-rRNA-SWiM approach in order to obtain a more detailed understanding of the quantitative role that specific microbial groups play in oil degradation.
This work provides a foundation for studies of the metabolic pathways of carbon assimilation in heterogeneous microbial consortia. Here, measurements obtained using the natural isotopic contrast between 13C-depleted oil and a 13C-enriched salt marsh ecosystem showed that hydrocarbons serve as a carbon source at the same time that their oxidation serves as a source of energy. The 13C content of purified bacterial rRNA showed the formation of 13C-depleted biomass over a period of 2 weeks. Most of the isotopically depleted signal was presumed to reflect carbon assimilated from the hydrocarbon amendment, and clonal enrichment of organotrophs is consistent with this observation.
Further study should allow us to investigate hydrocarbon uptake by environmental species using rRNA capture probes having greater specificity. Genus-specific 13C-rRNA-SWiM should be able to quantify the clones detected in the treatment samples, including Idiomarina spp., that reflect carbon assimilated directly from oil. We also noted that while the suite of RNA capture probes used here recovers most genera of bacteria, it is not universally compatible. Some excluded organisms of potential environmental interest are Pirellula (and presumably other Planctomycetales), Chloroflexus (and presumably other Chloroflexales), and Thermatogales (presumably including any of their low-temperature relatives). These exceptions highlight the need for group-targeted probes in future work.
The 13C-rRNA-SWiM approach offers several advantages over traditional substrate incubations. Such approaches provide a valuable new method for studying hydrocarbon assimilation by specific organisms and complement other methods, such as DNA isotope probing (52), labeled-substrate RNA capture (40), and analysis of labeled phospholipid fatty acids (8, 47). Additionally, the investment in technology and labor is arguably less resource-intensive than ion probe (secondary ion mass spectrometry) (45) techniques, although to date our laboratories house the only SWiM interfaces available. The next steps toward broader application of this approach, therefore, include development of new capture probes and broader collaborations between microbiologists and isotope geochemists at the interface of these analytical approaches. Routine
13C analyses of rRNA now require <50 ng C, which is equivalent to the RNA content of 104 to 106 cells, depending on the growth rate. This should open new windows for understanding carbon flows and pathways in complex microbial consortia.
| ACKNOWLEDGMENTS |
|---|
13C by closed-tube combustion. This work was supported by a grant from the Reinhart Coastal Research Center of Woods Hole Oceanographic Institution (to A.P. and K.J.E.) and by grant NSF-EAR-0311937 (to A.P.).
| FOOTNOTES |
|---|
Published ahead of print on 14 December 2007. ![]()
| REFERENCES |
|---|
|
|
|---|
| ||||||||