Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2002, p. 5537-5548, Vol. 68, No. 11
0099-2240/02/$04.00+0 DOI: 10.1128/AEM.68.11.5537-5548.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Fossil Fuels and Environmental Geochemistry and Centre for Molecular Ecology, University of Newcastle, Newcastle upon Tyne NE1 7RU,1 National Environment Technology Centre, AEA Technology, Abingdon OX14 3ED,3 National Environment Technology Centre, AEA Technology, Didcot, Oxfordshire OX11 OQJ, United Kingdom,4 Bedford Institute of Oceanography, Dartmouth, Nova Scotia B2Y 4A2, Canada2
Received 24 April 2002/ Accepted 20 August 2002
|
|
|---|
|
|
|---|
Bacteria are considered to represent the predominant agents of hydrocarbon degradation in the environment (27), and hydrocarbon-degrading bacteria are ubiquitous. More than 20 genera of marine hydrocarbon-degrading bacteria, distributed over several (sub)phyla (
-, ß-, and
-proteobacteria; Gram positives; Flexibacter-Cytophaga-Bacteroides) have been described so far (5, 13, 16-18, 22, 48). As a single species typically is capable of degrading only a limited number of the compounds found in crude oil, a consortium composed of many different bacterial species is usually involved in oil degradation.
Because of the high carbon content of oil and the low level of other nutrients essential for microbial growth, the rate and extent of degradation are, in general, limited by the low availability of nitrogen and phosphorus (2, 33). Consequently, growth of hydrocarbon-degrading bacteria and hydrocarbon degradation can be strongly enhanced by fertilization with inorganic N and P. This has proven an effective bioremediation treatment on several types of shorelines (4, 40, 41, 46).
Bioremediation studies have, in general, been dominated by an empirical approach, and optimum nutrient amendment levels are often informed by laboratory incubations. In the field, care must be taken in supplying optimum concentrations of inorganic nutrients. Too high concentrations may result in eutrophication, and too low concentrations may result in suboptimal biodegradation. A better understanding of the systematic effects of nutrient amendment on biodegradative microbial populations and the progress of bioremediation would assist the development of more rational bioremediation strategies (21, 37). Therefore, laboratory beach microcosm experiments were performed to determine the effects of different levels of N and P supply on oil degradation and bacterial community dynamics. Microbial communities were characterized by using cultivation-independent molecular techniques (16S rRNA gene [rDNA] PCR-based denaturing gradient gel electrophoresis [DGGE] [31] and clone libraries). Numerical analysis (34) was used to assess changes in microbial community structure and the effects of bioremediation treatments.
|
|
|---|
Microcosm experiments.
Sediment was homogenized with synthetic seawater (Instant Ocean; Aquarium Systems Inc.) in the laboratory, and 2.0 kg of wet sediment was placed into a microcosm (Fig. 1) and held at 20 ± 3°C. Fresh synthetic seawater provided each beach microcosm with two tidal cycles each day (12 h apart, 1 liter of seawater per cycle), with the ebb tide waters drained from each microcosm into 20-liter acid-washed plastic bottles. To allow equilibration of physical and chemical parameters, the microcosms were subjected to six tidal cycles before the addition of oil. Oil spilled at sea and washed onto a beach was simulated by weathering forties crude oil by distillation at 250°C in order to remove low-molecular-weight volatile hydrocarbons (e.g., <nC11 and BTEX [benzene, toluene, ethylbenzene, xylene]). Total resolvable hydrocarbons (TRH), nC11-to-nC35 alkanes, and aromatic hydrocarbons contributed, respectively, 31.2% ± 3.8%, 12.7% ± 0.6%, and 1.9% ± 0.2% of the total petroleum hydrocarbons (TPH) in the weathered oil. The weathered oil was mixed vigorously with synthetic seawater (25%, vol/vol), to form a stable water-in-oil emulsion. The emulsion was added to the microcosms at high tide at a level of 3.7 kg/m2 by pouring it onto the seawater. Inorganic nutrient solutions were added 24 h after oil addition and also after 7, 14, and 21 days. Each nutrient addition (50 ml in tap water) was calculated to provide a defined percentage of N (as sodium nitrate) and P (as potassium dihydrogen phosphate) relative to the total oil mass on a weight basis. The different treatments applied to the microcosms are shown in Table 1. Since the microcosms were labor intensive and resource demanding, only a few treatments were replicated (Table 1). Beach sediment treated with fertilizer alone was used as a control rather than untreated sediments because previous experiments had shown that nutrient-amended beach sediments produced DGGE profiles identical to those of untreated beach sediments and that the profiles obtained were stable over a 108-day sampling period (42).
![]() View larger version (22K): [in a new window] |
FIG. 1. Design of the microcosms used in this study.
|
|
View this table: [in a new window] |
TABLE 1. Overview of treatments of laboratory beach sand microcosms and molecular analyses of samples taken from the microcosms
|
Carbon dioxide measurements.
The evolution rate of carbon dioxide from each microcosm was determined before oil was added and daily thereafter at the same point in the tidal cycle. The headspace of a sealed microcosm was circulated through the cell of an infrared gas analyzer (Servomex) by a method detailed by Swannell et al. (39). The device was calibrated by using a CO2-free-air standard and a compressed-air cylinder containing 335 ppm (by volume) CO2. The stability of the analyzer over the duration of the measurements (<5 min) was determined by circulating the headspace of an empty glass bottle through the sample cell. The drift was always less than 4 ppm.
Nutrient analysis.
As an indirect measure of nutrient concentrations in pore water in the microcosms, inorganic nitrogen concentrations in the seawater collected after the simulated tidal cycles were determined by using a Technicon autoanalyzer system (28).
Oil chemistry.
Hydrocarbons in oiled sediments (10 g) spiked with squalane and 1,1-binaphthyl standards to determine extraction efficiency were extracted, analyzed by gas chromatography with flame ionization detection and mass spectrometric detection, and quantified as described previously (39). The efficiency of recovery of the added standards was, on average, 83%. Replicate analyses showed that the variability of measured values was always less than 10%. To distinguish between physical removal and biodegradation, TPH, total GC resolvable hydrocarbons, nC11 to nC35 alkanes, and polycyclic aromatic hydrocarbons (PAH) were expressed relative to 17
(72),25ß(72)-hopane, a degradation-resistant compound present in crude oil (4). The percentage biodegradation was calculated by dividing the concentration of individual compounds relative to that of 17
(72),25ß(72)-hopane at the end of the experiment by the concentrations relative to that of 17
(72),25ß(72)-hopane at the start of the experiment.
Statistical analysis of data on oil chemistry, inorganic nutrients, and carbon dioxide evolution.
Statistical analysis (parametric two-way analysis of variance, Pearson correlation) was performed by using Systat 7.0 (SPSS Inc.).
DNA extraction, PCR, and DGGE analysis.
DNA was extracted from 0.5-g samples (Table 1) by using the bead beating method described by Curtis and Craine with a Mikrodismembrator-U (B. Braun Biotech) (10). PCR was performed in a total volume of 50 µl containing 0.2 µM primer Vf-GC (corresponding to positions 341 to 358 of the Escherichia coli 16S rRNA [31]), 0.2 µM primer Vr (corresponding to positions 534 to 517 of the E. coli 16S rRNA [31]), 0.2 mM deoxynucleoside triphosphates, 1 U of BioTaq enzyme and the buffer supplied with the enzyme (Bioline, London, United Kingdom), and 1 µl of template DNA. Amplification was performed by using a Hybaid Omnigene Thermocycler as follows: 95°C for 3 min, followed by 30 cycles of 94°C for 0.5 min, 55°C for 1 min, and 72°C for 1 min, with a final elongation of 72°C for 10 min. DGGE was performed with the Bio-Rad DCode system. The PCR product was loaded onto 1-mm-thick 10% (wt/vol) polyacrylamide (37.5:1 acrylamide-bisacrylamide) gels containing a 30 to 60% linear denaturing gradient. One hundred percent denaturant is 7 M urea and 40% (vol/vol) deionized formamide. Gels were run in 1x TAE buffer (40 mM Tris-acetate, 1 mM Na-EDTA, pH 8.0) at 60°C and 200 V for 3.5 h. Gels were stained in 1x TAE buffer containing SYBR Green I (diluted 1:10,000; Sigma) and photographed under UV transillumination.
Statistical analysis of DGGE tracks.
Scanned negatives were analyzed by using Quantity One 4.1 software (Bio-Rad, Hercules, Calif.), and data were exported to Excel and analyzed in Systat 7.0 (SPSS Inc.). Similarities between tracks were calculated by using the Dice coefficient (band based) or the Pearson product-moment correlation coefficient (whole densitometric curve based) (34). Since cluster analysis of the resulting similarity matrices does not allow conclusions to be drawn regarding the statistical significance of differences between clusters or groups of samples (38), similarity coefficients from matrices were assigned to different groups and subsequently tested to determine whether their means were significantly different. Nonparametric analysis of variance (Mann-Whitney U test) was applied since similarity coefficients are not necessarily normally distributed (38).
Cloning, sequencing, and phylogenetic analysis of 16S rDNA.
Almost full-length 16S rRNA gene fragments were amplified by using primers pA and pH' (12). Except for the primers, the PCR conditions were similar to those described above. PCR products were cloned by using the AdvanTAge kit (Clontech, Palo Alto, Calif.), and the 16S rRNA gene libraries were screened by amplified ribosomal DNA restriction analysis (ARDRA). The number of clones screened from each library is given in Table 1. E. coli clones were categorized into different ARDRA types based upon the pattern obtained on simultaneous digestion (3 h, 37°C) with the restriction enzymes RsaI and HaeIII (5 U of each). Sequencing of several clones corresponding to dominantly occurring ARDRA types revealed that these ARDRA groups were internally homogeneous (data not shown). ARDRA types were named X-Yd-Z, in which X refers to the type of microcosm (Table 1), Y is the day of sampling (0, 6, or 26), and Z is the number of a representative clone in the clone library of the particular microcosm and sampling date. The distribution of ARDRA types present in different clone libraries was determined and used to calculate the Shannon-Weaver index {H = -
[ni · log(ni)], where ni is the relative contribution of clone type i to the whole library}. ARDRA types occurring more than once in a library were selected for sequence analysis. Sequence data were obtained with a single primer (pD', E. coli positions 534 to 518 [12]), generating 0.5 kb of sequence data. For a selection of clones, the nearly complete 16S rDNA was sequenced in both directions. Sequences were compared to sequences deposited in the GenBank DNA database by using the BLAST algorithm (1). Alignments were performed by using ClustalW and corrected manually. Distance analysis using the Jukes and Cantor correction (24) and bootstrap resampling (100 times) was done with the TREECON package (45), and the distance matrix was used to construct a tree via neighbor joining (35). Parsimony analysis was done by using DNAPARS from the PHYLIP package (14).
Nucleotide sequence accession numbers.
The nucleotide sequences determined in this study have been deposited in the GenBank database under accession numbers AF432269 to AF432344.
|
|
|---|
![]() View larger version (45K): [in a new window] |
FIG. 2. Chemical data from microcosms. (A) amounts of nitrogen in residing seawater of beach microcosms over time. Symbols: , 0% N; , 1% N; , 2% N; , 4% N; , 10% N. Arrows indicate addition of nutrients. (B) Daily carbon dioxide production, averaged over three independent runs, with bars indicating standard deviation, for microcosms. Symbols: , FO; , 0% N; , 0.75% N. For the FO and 0% N microcosms, standard deviation symbols were often smaller than the data symbols. Arrows indicate addition of nutrients. (C) Cumulative carbon dioxide production in microcosms. The open portion of the column indicates production during the first 7 days; carbon dioxide production during the last 21 days of the experiment is shown by the hatched portion of the column. Standard deviations are shown by error bars for the microcosms prepared in triplicate (FO, 0% N, and 0.75% N). (D) Percentage of biodegradation of TPH (open), TRH (black), n11-to-n33 alkanes (hatched), and PAH (grey) per treatment. Error bars indicate standard deviations.
|
Figure 2D shows the percentages of biodegradation of TPH, TRH, nC11-to-nC35 alkanes, and PAH. Significantly greater biodegradation of these compound classes was obtained with all of the bioremediation treatments than with the oil-only control (P < 0.001). TRH, nC11-to-nC35 alkanes, and aromatic hydrocarbons were almost completely removed (91% ±2%, 98% ± 1%, and 92% ± 4%, respectively). No significant differences in TRH, n-alkane, and aromatic hydrocarbon biodegradation were observed between the various nutrient amendments (P > 0.05). However, TPH degradation was significantly greater in the 6% N microcosm than in the 0.25% N, 4% N, and 10% N microcosms (P < 0.05).
Carbon dioxide production in the oiled, untreated controls (0% N) was similar (P > 0.05) to that in the unoiled microcosms treated with fertilizer (FO; Fig. 2B and C), in agreement with the absence of significant biodegradation of TPH, TRH, and alkanes (Fig. 2D).
Effect of nutrient amendment on bacterial community structure and dynamics.
Effect of bioremediation treatments on bacterial community structure was determined by using 16S rDNA-based PCR-DGGE for a number of treatments and sampling times (Table 1). In order to determine relationships between the community fingerprints, similarities were calculated on the basis of the absence or presence of bands (Dice coefficient) and on the basis of whole-track curve densitometric information (Pearson product-moment coefficient). Results based on the two coefficients were comparable; therefore, only the band-based similarities will be described. Analysis of triplicate samples taken from a single microcosm at a single sampling time showed that microbial communities were highly similar (95% ± 5%), indicating low spatial heterogeneity within the microcosms. However, obvious differences in community structure were evident both in time (Fig. 3 and 4) and between different treatments (Fig. 5 and 6). No obvious changes in the composition of the predominant bacterial communities occurred over time in the microcosm treated with fertilizer but no oil (FO). Although DGGE analysis of replicate microcosms for the FO treatment was not conducted in this study, data from a related field experiment that incorporates a fully replicated randomized block design indicate that DGGE profiles from the same beach sediments used in this study and treated with fertilizer only are highly similar over time and between replicate blocks (data not shown). The largest changes occurred in the microcosms remediated with 1 and 4% N (Fig. 4); the average similarity relative to day 0 was significantly lower for these two microcosms than in the other three oiled, nutrient-amended microcosms (P < 0.05). The most dramatic changes in community structure occurred during the first 6 days following oil addition, and profiles between days 6 and 26 were significantly more similar to each other than to profiles obtained from day 0 samples (Fig. 4; P < 0.05), suggesting that, following an initial significant change in community composition, the bacterial community structure remained relatively stable. Nevertheless, for all oiled microcosms, community changes after 6 days were significantly larger than in the microcosm that was not treated with oil (FO) over the same time period (P < 0.05).
![]() View larger version (100K): [in a new window] |
FIG. 3. Bacterial community dynamics during oil spill bioremediation in beach microcosms. Examples of changes in 16S rDNA PCR-DGGE profiles over time in microcosms; unoiled microcosms (FO), oiled, untreated microcosms (0% N); and oiled microcosms amended with 0.25 and 10% nitrogen. The values above the lanes indicate the numbers of days elapsed since nutrient addition.
|
![]() View larger version (39K): [in a new window] |
FIG. 4. Similarities (Dice coefficient) in DGGE profiles from microcosms receiving different treatments. Open columns show average similarities of profiles observed at days 6, 13, 20, and 26, relative to day 0, for each treatment. The average values were obtained from pairwise comparisons of DGGE profiles obtained at each time point relative to the time zero sample for each individual treatment. Thus, the data for the FO treatment, for example, represent the mean similarity obtained from four pairwise comparisons (the mean of the day 0 profile compared to the day 6 profile, the day 0 profile compared to the day 13 profile, the day 0 profile compared to the day 20 profile, and the day 0 profile compared to the day 26 profile). Error bars (n = 4) indicate standard deviations. Hatched columns show average similarity between profiles observed at days 6, 13, 20, and 26 for each treatment. In this case, the mean values were obtained by averaging the pairwise similarities of the day 6 profile compared to the day 13 profile, the day 6 profile compared to the day 20 profile, the day 6 profile compared to the day 26 profile, the day 13 profile compared to the day 20 profile, the day 13 profile compared to the day 26 profile, and the day 20 profile compared to the day 26 profile. Error bars (n = 6) indicate standard deviations. These data indicate that most of the change in the DGGE profiles occurred between days 0 and 6, since the mean similarities for all pairwise comparisons from day 6 onward show little variation.
|
![]() View larger version (124K): [in a new window] |
FIG. 5. Comparison of the effects of nutrient amendment on bacterial community structure at particular sampling times. Examples of 16S rDNA PCR-DGGE profiles from microcosms receiving different nutrient amendments, showing profiles at days 13 and 26. The values above the lanes indicate treatments (described in Table 1).
|
![]() View larger version (21K): [in a new window] |
FIG. 6. Similarities (Dice coefficient) in DGGE profiles between different treatments at each time point. Shown are average similarities of profiles observed for oiled, nutrient-amended microcosms (0.25, 0.5, 1, 4, and 10% N) relative to the oiled, unamended control ( ) and unoiled, nutrient-amended control ( ), respectively; average similarity between the different oiled, nutrient-amended microcosms ( ); and average similarity between DGGE profiles from replicate microcosms receiving the same nutrient amendment (0.75% N) ( ). Error bars indicate standard deviations. To aid in visualization, symbols representing different comparisons have been slightly offset. For each time point, the average values were obtained from pairwise comparisons of DGGE profiles obtained for each treatment relative to the oiled, unamended control. Thus, the data for the day 6 time point, for example, represent the mean similarity obtained from five pairwise comparisons (the mean of the 0.25% N profile compared to the 0% N profile, the 0.5% N profile compared to the 0% N profile, the 1.0% N profile compared to the 0% N profile, the 4.0% N profile compared to the 0% N profile, and the 10% N profile compared to the 0% N profile). Similar comparisons were made between the fertilizer-only control and the other treatments. For microcosms that received oil and different nutrient amendments, all possible pairwise similarities were averaged to obtain the mean value. These data indicate that, at each time point other than day 0, the variation in DGGE profiles between replicate microcosms ( ) is as great as the variation between different treatments ( ).
|
Phylogenetic analysis and diversity of bacterial communities.
To obtain a more specific picture of which bacterial taxa were stimulated by the oil and bioremediation treatments, phylogenetic analysis of cloned 16S rRNA genes was performed. The microcosms and sampling times indicated in Table 1 were selected because DGGE analysis revealed the greatest differences in the predominant populations. Clone libraries were screened by ARDRA, and representatives of ARDRA profiles occurring more than once in a library were partially sequenced (Table 2). It should be noted that the samples from the 0 and 4% N microcosms taken at day 0 were frozen immediately after oil addition and, although treated with oil, the bacterial populations would not have had time to respond to the oil addition and, in effect, represent bacterial populations not yet affected by oiling. For the five libraries from microcosms not affected by oiling (day 0 libraries of the FO, 0% N, and 4% N microcosms and day 6 and day 26 libraries of the FO microcosm), phylogenetic analysis of the clones occurring more than once in the libraries indicated that they belonged to a wide range of phylogenetic groups, with none dominating the clone libraries (Table 2 and Fig. 7). Almost half of the clones were unique, meaning that their ARDRA type was only encountered once in a clone library. The purpose of this study was not to catalogue the wider diversity of the bacterial populations present in the beach sediments, and these were not analyzed further. The high diversity was also indicated by the fact that, of the 49 clones sequenced, clones with similar ARDRA profiles and >99% sequence identity were encountered only twice in different clone libraries (indicated by the letters A and B in Table 2). In contrast, oiling and bioremediation resulted in strong dominance by a few ARDRA/sequence types (Table 2 and Fig. 7). This decrease in diversity was also illustrated by the calculation of the Shannon-Weaver index for the clone libraries (Fig. 8). Diversity remained high in the microcosm treated with fertilizer only. However, oiling strongly decreased the bacterial community diversity during the first 6 days of the experiment. The strongest effect was observed for the bioremediation treatment. Nevertheless, between days 6 and 26, diversity began to recover to preoiling levels (Fig. 8).
|
View this table: [in a new window] |
TABLE 2. Overview of ARDRA types occurring more than once in a clone library, their relative contributions to the clone library, and the closest relative in the GenBank database with similarity to the partially sequenced clone (0.4 kb)a
|
![]() View larger version (27K): [in a new window] |
FIG. 7. Relative contributions of different phylogenetic groups to microbial communities. (A) Microcosms not affected by oiling and remediation treatment (composed of five clone libraries [FO days 0, 6, and 26 and day 0 of the 0% N and 4.0% N microcosms]). (B) Oiled, untreated control (0% N) at days 6 and 26. (C) Oiled, 4% N-treated microcosm at days 6 and 26. CFB, Cytophaga-Flexibacter-Bacteroides; GNS, green non-sulfur bacteria.
|
![]() View larger version (14K): [in a new window] |
FIG. 8. Effect of oiling and bioremediation on changes in microbial diversity over time. The Shannon-Weaver index was calculated from the distribution of clones categorized on the basis of different ARDRA profiles. Symbols: , FO; , 0% N; , 4% N.
|
-proteobacteria (62.2% of the clones in the 0% N microcosm and 73.1% of the clones in the 4% N microcosm; Fig. 7).
-Proteobacteria (50.6% of the clones) continued to dominate the library from the microcosm treated with oil alone (0% N) after 26 days, but in the bioremediated microcosm (4% N),
-proteobacteria became dominant (63.3%). The nearly complete cloned 16S rDNA of representatives of the dominantly occurring ARDRA types was sequenced. The dominant
-proteobacterial clones (4%N-26d-5, -7, and -27) belonged to the aerobic anoxygenic phototrophic bacteria and were phylogenetically most closely related to Erythrobacter longus and E. citreus (Table 2 and Fig. 9A). ARDRA profiles indicative of these sequences were only found in the sample taken 26 days after 4% N fertilizer treatment and represented 50.0% of the clones. A different ARDRA type (0% N-26d-88) most closely related to Erythrobacter spp. was found in the day 26 library of the 0% N library, constituting 4.2% of the library (Table 2). Two dominant ARDRA types, found in several clone libraries (indicated by C and D in Table 2 and Fig. 9B) represented the alkane-degrading Alcanivorax/Fundibacter group (5, 48) of the
-proteobacteria. Sequences of ARDRA type C were distinct from Alcanivorax borkumensis and Fundibacter jadensis but formed a distinct group with these sequences (Fig. 9B). Type D sequences were almost identical to A. borkumensis (>99.8%). ARDRA profiles and sequences indicative of type C comprised 73 and 36% of the libraries constructed from the 4% N-amended microcosm at day 6 and the oil-only control 0% N microcosm at day 6, respectively. However, they were not detected in any other library. ARDRA profiles and sequences indicative of type D were only found in the libraries from microcosms treated with oil only (0% N) and comprised 4 and 32% of clones analyzed after 6 and 26 days, respectively. A small number of other clones most closely related to microorganisms known to degrade components of crude oil were also detected and are indicated in Table 2. All of these sequences constituted less than 6% of the clones in a library.
![]() View larger version (34K): [in a new window] |
FIG. 9. Phylogenetic tree based on almost complete 16S rDNA gene sequences of -proteobacterial (A) and -proteobacterial (B) clones occurring as dominant sequences in clone libraries of oiled microcosms. A neighbor-joining analysis with Jukes-and-Cantor correction was performed. Only bootstrap values of greater than 50% are shown.
|
|
|
|---|
N and P were supplied at different concentrations, which, according to resource ratio theory (44), should select different microbial communities dominated by the organisms most capable of utilizing the inorganic nutrients at the levels added to the polluted habitat. In principle, this could offer opportunities to direct oil degradation (21, 37). Clear-cut differences in bacterial community structure were found in microcosms treated with different levels of inorganic nutrients, which appears consistent with the predictions of resource ratio theory. However, statistical analysis of DGGE profiles revealed that the average similarities between DGGE profiles from three independently prepared microcosms receiving the same amount of inorganic nutrients (0.75% N) were not significantly different from the similarities between microcosms subjected to different nutrient amendments. Therefore, it cannot be concluded that the different communities selected resulted from differences in nutrient amendment alone. It is unlikely that heterogeneity within the microcosms, as the result of stratification or poor mixing of oil and sediment/water, contributed significantly to the observed highly variable composition in bacterial community structure between microcosms (only 65% ± 6% similarity). Replicate samples taken at the same time from a single microcosm, which would be subject to the same constraints on mixing and stratification as separate microcosms, were highly reproducible (95% ± 5%). Also, the community profiles between days 6 and 26, after major changes in bacterial communities had occurred, were similar for single microcosms (80% ± 4%).
Nevertheless, the results revealed that communities with highly different compositions (similarities of only 58 to 71%) showed similar rates and extents of oil degradation at different nutrient concentrations. This observation may be due largely to the strong selection for a few members of the Alcanivorax/Fundibacter group, which are capable of degrading alkanes (5, 48), a major oil component. Comparable observations have been made for a functionally stable methanogenic reactor fed with glucose (15). This reactor revealed considerable dynamics in microbial communities over a 605-day period, despite minimal differences in performance over time.
Bioremediation treatments are aimed at stimulating pollutant-degrading microorganisms to speed the recovery of contaminated ecosystem to a prepollution state in terms of biodiversity and ecosystem function. In this study, changes in the predominant bacterial populations occurred in all of the microcosms, except that which received inorganic nutrients but not oil. Oiling and especially bioremediation led to a strong decrease in bacterial community diversity at day 6, but a rapid recovery to near preoiling levels of diversity occurred subsequently. Still, despite having a similar level of biodiversity, the component organisms contributing to that diversity were somewhat different from the original community, as revealed by DGGE analysis and clone libraries. Following bioremediation,
-proteobacteria were dominant by day 26 of the experiment, whereas prior to oil contamination, they had only represented 8% of the clones in the libraries. Sequences from
-proteobacteria most closely related to Erythrobacter spp. were most commonly encountered. Erythrobacter spp. belong to the aerobic anoxygenic phototrophic bacteria (49). They metabolize organic substrates, with light enhancing their growth. Recently, it was revealed that aerobic anoxygenic phototrophic bacteria play a critical role in the carbon cycle in the ocean (26). Aerobic anoxygenic phototrophic bacteria grow well in high-nutrient media and are often isolated from environments with a high level of organic matter content (49), as well as from intertidal beach sediments (36). Recently, a PAH-degrading bacterium most closely related to Erythrobacter spp., Lutibacterium anuloederans LC8, was isolated from intertidal beach sediments (8). The rRNA gene of L. anuloederans shows 95.3 to 96.1% similarity to our Erythrobacter-like sequences. An
-proteobacterial sequence related to Erythrobacter spp. was also detected in beach sediments from an oil spill bioremediation field experiment conducted in Delaware (30).
Perhaps the most dramatic observation made in the present study was the rapid and strong selection for
-proteobacteria in oil-treated microcosms. The
-proteobacteria persisted in clone libraries from sediment samples treated with oil alone taken on day 26, while they were replaced by
-proteobacteria in microcosms treated with oil and inorganic nutrients. This may reflect slow ongoing alkane degradation in the microcosms not treated with inorganic nutrients and selection of organisms growing on residual PAH and secondary products of hydrocarbon degradation in the bioremediated microcosms. A field experiment on shoreline sediments in the Norwegian Arctic (19) revealed an increase in
-proteobacteria in oil-contaminated beach sediments, while the microbial community of beached oil paste after the Nakhodka oil spill accident in the Japan Sea was also dominated by
-proteobacteria (25). These results are consistent with our observation of a strong dominance by
-proteobacteria, and this may be a characteristic of bacterial communities associated with recently spilled oil. In our case, it was mainly members of the Alcanivorax/Fundibacter group that were selected. Alcanivorax borkumensis (48) and Fundibacter jadensis (5) were described recently. Both are capable of using only a few organic substrates, especially alkanes (5, 48) and the alkyl groups of n-alkylbenzenes and n-alkylcycloalkanes (11). Alkanes are major oil components, explaining why Alcanivorax/Fundibacter-like sequences made a significant contribution (up to 73%) to clone libraries early in the bioremediation process and for a longer period when the hydrocarbon degradation rate was lower. Because of biases associated with PCR and cloning (47), these clone percentages probably do not directly represent Alcanivorax cell numbers, but since all samples were analyzed in the same fashion and the predominance of Alcanivorax-like sequences in the clone libraries was so dramatic, it is likely that our results do reflect a genuine increase in the relative abundance of this group of bacteria. Furthermore, fluorescence in situ hybridization analyses have shown that Alcanivorax constituted more than 90% of the entire microbial community in laboratory incubations of seawater supplemented with oil (43).
Although in the contaminated microcosms without nutrient treatment (0% N), no significant oil degradation occurred over the 28-day incubation period, remarkably, the bacterial communities changed considerably and the predominance of Alcanivorax/Fundibacter-like sequences in the corresponding clone library was apparent. This may be explained by the release of nitrogen and phosphorus from biomass killed by the oil, supporting limited hydrocarbon degradation. Interestingly, hydrocarbon pollutants can induce prophages, resulting in lysis of a large proportion of the bacterial community (9, 23), and this may explain changes in community composition independent of extensive oil degradation. Even in the absence of significant hydrocarbon degradation, Alcanivorax/Fundibacter-like sequences constituted about 40% of the clone libraries constructed from samples taken at days 6 and 26 from the unamended, oil-contaminated control microcosm. To assess if this is a consequence of relatively very low rates of hydrocarbon degradation by Alcanivorax/Fundibacter-like bacteria under these conditions, information on their activity would be required, for example, measurement of mRNA that encodes enzymes involved in biodegradation.
The members of the Alcanivorax/Fundibacter group appear to have a cosmopolitan distribution, as their presence has been noted in coastal waters and on beaches of the United Kingdom (this study), Germany (5, 48), Italy (GenBank accession numbers AB0302701 to -4), Singapore (GenBank accession number AF062642), Japan (20, 25), and the United States (7, 32). They have been detected mainly in laboratory enrichments with oil components (5, 7, 20, 48), but more importantly, their occurrence in oil paste and seawater has been described following an oil spill in the Japan Sea (25). The worldwide distribution of Alcanivorax spp. indicates that they may be of considerable global significance in marine hydrocarbon degradation, and a project to sequence the genome of this important organism has recently been instigated (http://www.uni-bielefeld.de/presse/pm/engpm31.htm).
|
|
|---|
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»