Department of Environmental Microbiology, UFZ Centre for Environmental Research Leipzig-Halle, Leipzig, Germany
Received 3 November 2005/ Accepted 16 February 2006
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The degradation of hydrocarbons as major components of mineral oil is generally brought about by microorganisms. At the salinity of seawater (
3% total salts, mainly NaCl) or below, most low-molecular-weight components of mineral oil are easily degraded by a variety of microorganisms. However, the range of degrading organisms decreases with increasing salinity. There are only a few reports on the degradation of hydrocarbons in hypersaline environments. Enrichment cultures from the Great Solar Lake grew on mineral oil at salinities up to 17.2% (70), and a bacterial consortium and an extremely halophilic archaeon grew on hexadecane and crude oil at salinities of 15% and 32%, respectively (5). Degradation of hydrocarbons (mainly n-alkanes) by archaea (Halobacteria) at 15 to 32% NaCl was also reported (40). An extremely halophilic archaeon from a salt marsh degraded various saturated and aromatic hydrocarbons at salinities of up to 31% and grew optimally at 22% NaCl (10). In our own laboratory, microbial communities from saline soils in Patagonia degraded diesel fuel up to a salinity of 17.5% (59). Although some members of the communities were identified, their contributions to the diesel fuel biodegradation remained obscure.
Information on the phylogenetic diversity of microbial communities can be obtained by molecular methods, like fingerprinting or cloning and sequencing of PCR-amplified rRNA genes. These methods have been applied to isolates and consortia enriched from natural hypersaline environments, such as salt lakes, solar salterns, and soda lakes (1, 8, 24, 52). The 16S rRNA approach was also used to assess the microbial diversity in hypersaline environments and to follow transitions in community structure upon seasonal fluctuations or along spatial gradients in salinity, temperature, and availability of substrates and minerals (3, 6, 7, 8, 9, 11, 20, 29, 43, 46, 52).
Whereas the metabolic capabilities of oil-degrading communities grown under hypersaline conditions have been investigated by community profiling with BIOLOG plates (43, 44), information on the contributions of community members to oil biodegradation is more difficult to obtain. PCR screening for functional genes and expression profiling based on mRNA using reverse transcriptase PCR or microarrays are options when the degradation pathways and the genes involved are known (for reviews, see references 71, 73, 77, and 78). Other culture-independent methods linking degradation activities with the identities of community members are based on the uptake of isotope-labeled substrates. Examples are a combination of fluorescence in situ hybridization and microautoradiography that has been applied to halophiles (60) and stable isotope probing of DNA or RNA (for a review, see reference 57). Multiparametric flow cytometry has the potential to identify the active community members without requiring knowledge of metabolic genes. It provides information on the physiological states and activities of individuals and populations within communities, e.g., by indicating proliferation activity, while raising the possibility of identifying individuals by specific fluorescence labeling (50) or by rRNA-based methods after cell sorting (49, 53).
In the present study, the structural dynamics of a bacterial consortium isolated from saline soil were characterized during growth on diesel fuel at salinities up to 20%. Biodegradation and growth were followed by measurement of oxygen consumption and by analyses of alterations in community composition by single-strand conformation polymorphism (SSCP) fingerprinting and community-level substrate utilization profiles using the BIOLOG system. Multiparametric flow cytometry and high-speed cell sorting were applied to separate emerging proliferating subcommunities and to identify the dominant phylotypes therein by 16S rRNA gene sequencing. The aim of this work was to establish a high-speed detection tool capable of determining the degradation potential within natural environments via the investigation of the diversity and dynamics of a bacterial community and the identification of those consortium members that were actively involved in the degradation process. The information can be used to support bioremediation processes in soils of contaminated saline, arid regions.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Isolation and propagation of the bacterial consortium.
The microorganisms were extracted after mixing and sonicating the soil in 0.2% Na4P2O7 solution (adjusted to pH 8.5) as described previously (58), except that 7.5% NaCl was added. The consortium was enriched in shaken flasks containing 250 ml mineral medium amended with diesel fuel (800 µl liter1), paraffin (400 µl liter1), and yeast extract (120 mg liter1). Clay particles were removed by repeated transfer into fresh medium. The composition of the mineral medium was (liter1) 870.9 mg K2HPO4, 680.5 mg KH2PO4, 760 mg NH4Cl, 71.2 mg MgSO4 · 7H2O, 5.47 mg CaCl2, 0.44 mg ZnSO4 · 7H2O, 0.812 mg MnSO4 · 4H2O, 0.785 mg CuSO4 · 5H2O, 0.252 mg Na2MoO4 · 2H2O, 4.48 mg FeSO4 · 7H2O. Finally, 7.5% NaCl was added to the medium.
Degradation experiments.
Biodegradation of the diesel fuel was continuously monitored by following oxygen consumption at 25°C in Sapromat D12 respirometers (Voith, Heidenheim, Germany). The 500-ml flasks contained 240 ml mineral medium with 0, 7.5, 15, or 20% NaCl; 250 µl (215 mg) diesel fuel as a sole carbon source; and 10 ml inoculum corresponding to 30 mg dry mass. Degradation experiments took place in 24 flasks that were sacrificed for chemical and biological analyses. One flask served as a substrate-free control. A second flask was sacrificed right after inoculation to obtain initial values for the protein content, CFU, cell numbers, and BIOLOG profiles. Six flasks each were cultivated at 15 and 20% NaCl, and five flasks each were cultivated at 0% and 7.5% NaCl. Entire flasks were harvested to obtain biomass for the determination of protein contents, CFU, cell numbers, BIOLOG profiles, SSCP fingerprints (after 42 and 84 days of cultivation at the different salt concentrations), and flow cytometry, as well as cell sorting (after 21, 42, 63, and 84 days of cultivation at the different salt concentrations). The oxygen consumption rate was continuously measured in two parallel flasks at 0 and 7.5% NaCl each and in three parallel flasks at 15 and 20% NaCl each. The residual diesel fuel content at the end of the experiment (after 84 days) was determined by gas chromatography (34).
Cell numbers and protein content.
CFU counts were obtained in triplicate on R2A agar (Difco) supplemented with 2% NaCl (for samples from the culture without NaCl addition) or 7.5% NaCl (for samples from all higher salinities), each amended with 70 mg liter1 cycloheximide. Total cell counts were determined by fluorescence microscopy after dispersion, dilution, staining with DAPI (4',6'-diamidino-2-phenylindole), and filtering on 0.2-µm Nuclepore polycarbonate filters in triplicate. The protein content was determined in duplicate according to the method of Bradford (12).
Substrate utilization profiles.
Community level substrate utilization profiles were obtained by the BIOLOG method (BIOLOG Inc., Hayward, CA), which indicates substrate utilization, since respiring organisms reduce tetrazolium violet to colorimetrically distinguishable formazan. Cell suspensions from degradation experiments were diluted 1:10 with 7.5% NaCl solution (with the exception that 3% NaCl solution was used for samples from the degradation experiment without NaCl addition) and used to inoculate BIOLOG GN and GP plates. The number of different substrates was 128, since both types of plates contained 95 substrates with an overlap of 62 substrates. The plates were incubated at 25°C and analyzed after 24, 48, 72, and 96 h.
Multiparametric flow cytometry.
Flow cytometric measurements were carried out using a MoFlo cell sorter (DakoCytomation, Fort Collins, CO) equipped with two water-cooled argon ion lasers (Innova 90C and Innova 70C from Coherent, Santa Clara, CA). Excitation of 580 mW at 488 nm was used to analyze the forward scatter (FSC) and side scatter (SSC) as a trigger signal at the first observation point, using a neutral-density filter with an optical density of 2.3. DAPI was excited by 180 mW of multiline UV (333 to 365 nm) at the second observation point. The orthogonal signal was first reflected by a beam splitter and then recorded after reflection by a 555-nm long-pass dichroic mirror and passage by a 505-nm short-pass dichroic mirror and a band pass of 488/10. DAPI fluorescence was passed through a 450/65 band pass filter. Photomultiplier tubes (models R 928 and R 3896) were obtained from Hamamatsu Photonics (Hamamatsu City, Japan). Amplification was carried out at linear or logarithmic scales, depending on the application. Data were acquired and analyzed using Summit software (DakoCytomation, Fort Collins, CO). The dot plots were gated in such a way that visibly distinct subcommunities were included in the gates and cell counts therein were calculated. Fluorescent beads (polybead microspheres [diameter, 0.483 µm; Flow Check BB/Green Compensation Kit, blue alignment grade, reference 23520], Polyscience) were used to align the MoFlo (coefficient of variation, about 2%). Also, an internal DAPI-stained bacterial cell standard was introduced for tuning the device up to a coefficient of variation not higher than 6%.
Cell sorting was performed using the four-way-sort option at high speed (12 ms1). The most accurate sort mode (single and one-drop mode; highest purity, 99%) was chosen for separating 5,000 cells per second. The cells were sorted into nucleic-acid-free glass flasks. The cells were separated from the mixed culture using DNA-DAPI fluorescence intensity and forward scatter signals in several independent experiments with different gate settings. Dominant and apparently growing subcommunities were separated in order to facilitate their molecular identification. Samples of incubations at all four NaCl concentrations were taken after 42 days and 84 days. Up to four subcommunities per sample could be separated simultaneously due to instrumental restrictions. To separate additional subcommunities, the sample taken after 84 days from the 20% NaCl incubation was run a second time. To obtain sufficient DNA for the generation of 16S rRNA gene clone libraries, at least 106 cells per subcommunity were sorted.
DNA content.
Cells were harvested from the cultures, centrifuged (3,200 x g; 6 min), immediately resuspended in sodium azide (10 µl of a 10% solution in bidistilled water), and stored at 4°C. This procedure was found to preserve the cells for at least 3 months. Just prior to being stained, the stored cells were washed in sodium chloride-phosphate buffer (400 mM Na2HPO4-NaH2PO4, 150 mM NaCl, pH 7.2) and diluted to an optical density at 700 nm of 0.035 in the same buffer. For DNA determination, 2 ml of diluted cell suspension was treated with 1 ml solution A (2.1 g citric acid-0.5 g Tween 20 in 100 ml bidistilled water) for 10 min, washed, and resuspended in 2 ml solution B (0.24 µM DAPI [Sigma], 400 mM Na2HPO4, pH 7.0) for at least 10 min in the dark at room temperature using a modification of a standard procedure (47). The staining method is quantitative for individuals in a population, since each cell binds an amount of the AT-specific DNA dye DAPI that is proportional to its AT content. Based on knowledge of the cell cycle of a population, it allows actively proliferating cells to be distinguished, because the staining intensity correlates with the number of chromosomes of the individuals within a population. Proliferation activity was thus inferred from high chromosome numbers in populations and from signs of a split-up of populations into subpopulations of different chromosome contents, visible as stacks or towers of subcommunities along the DNA axis in cytometry dot plots (51). The microscopically proven absence of cell aggregation facilitated cell sorting into clearly separate subcommunities.
Scatter behavior.
FSC is related to cell size, and SSC is related to cell granularity. Data were obtained by examining the light-scattering behavior of individual cells, mediated by the 488-nm line of the argon ion laser. Usually, shifts in the FSC/SSC mean values were recorded in parallel with the fluorescence measurements. At least 100,000 cells per sample were analyzed.
Microscopy.
To ensure staining of all cells within the consortium, blue DAPI fluorescence was observed with an epifluorescence microscope equipped with a UV light from a 100-W mercury arc lamp (Axioscope, Germany) using a UG1 filter cube (excitation filter 20 UG1; IF 17L420) and recorded with a DXC-9100P camera (Sony, Japan). Micrographs were subjected to image analysis with the Openlab 3.1.4 software (Improvision).
SSCP fingerprinting and reamplification of SSCP bands.
SSCP fingerprints were obtained from the inoculum and after 42 and 84 days of incubation with diesel fuel at various salinities. The SSCP analyses were performed according to the method of Schwieger and Tebbe (62), with some modifications. Cells were harvested by centrifugation, and total genomic DNA was extracted with a FastDNA SPIN Kit for soil (QBIOgene, Germany), according to the manufacturer's instructions. The V4-V5 region of the bacterial 16S rRNA gene was amplified using the primers UniBac515f (5'-GTG CCA GCA GCC GCG-3') and UniBac927r-Ph (5'-CCC GTC AAT TYM TTT GAG TT-3'); the reverse primer was phosphorylated at the 5' end. PCR was performed in 100-µl reaction mixtures containing 50 µl Taq PCR Master Mix (QIAGEN, Germany), 3 mM MgCl2, 10 pmol of each primer (supplied by MWG Biotech, Germany), and 1 µl template DNA (in dilutions of 1, 1:10, and 1:100) in a PTC-200 Thermal Cycler (MJ Research). The cycle parameters were as follows: 3 min at 94°C and 30 cycles of 20 s at 94°C, 30 s at 54°C, and 1 min at 72°C, followed by a 10-min extension step at 72°C. After the results were checked by electrophoresis on a 1.5% agarose gel, the PCR products were purified with an E.Z.N.A. Cycle-Pure Kit (peqLab Biotechnologie GmbH, Germany), eluted in 30 µl bidistilled water, and digested with 6 units Lambda exonuclease (New England Biolabs, Germany) for 1 h at 37°C. After addition of 1 volume of 2x SSCP sample buffer (95% formamide, 10 mM NaOH, 0.25% bromophenol blue, 0.25% xylene cyanol FF), the samples were denatured for 2 min at 95°C, immediately chilled on ice for 5 min, and stored at 20°C until electrophoresis. A 0.6% MDE gel (Biozym, Germany) was prepared on a Polybond film (Biometra, Germany) according to the supplier's protocol. Electrophoresis was run in a TGGE Maxi chamber (Biometra, Germany) in 1x Tris-borate-EDTA buffer (61) at 400 V and 26°C for 16 to 20 h. The gel was silver stained according to the method of Bassam et al. (4) with the modification that 3% NaOH was used as a developing bath. Bands of interest were cut from the wet gel, and the DNA was extracted as described previously (26). The extraction procedure was repeated, the supernatants were pooled, and DNA was precipitated with ethanol (61) and dissolved in 10 µl bidistilled water. For reamplification, 1 to 4 µl of the DNA was applied as a template in a 25-µl reaction mixture as described above, with the exception that the reverse primer, UniBac927r, was not phosphorylated. Reamplification was performed with 35 cycles according to the protocol described above. After agarose gel electrophoresis and purification, the PCR products were used for direct sequencing or cloning.
16S rRNA gene cloning and ARDRA screening of clone libraries.
Bacterial 16S rRNA gene fragments were amplified by PCR using the primers 27F and 1492R (41). PCR was performed in 25-µl reaction mixtures containing 12.5 µl Taq PCR Master Mix (QIAGEN, Germany), 5 pmol of each primer (supplied by MWG Biotech, Germany), and 1 to 4 µl template DNA (extracted from sorted cells after 42 and 84 days of cultivation, as described above) with a PTC-200 Thermal Cycler (MJ Research). The cycle parameters were as follows: 4 min at 94°C and 30 cycles of 45 s at 94°C, 1 min at 58°C, and 2 min at 72°C, followed by a 10-min extension step at 72°C. After the results were checked by electrophoresis on a 1.5% agarose gel, the PCR products were purified with an E.Z.N.A. Cycle-Pure Kit (peqLab Biotechnologie GmbH, Germany) and cloned using a QIAGEN PCR Cloning Kit. For each sample of sorted cells, about 200 clones were picked and screened for the appropriate insert size by PCR using the vector-specific primers M13uni(21) (5'-TGT AAA ACG ACG GCC AGT-3') and M13rev(29) (5'-CAG GAA ACA GCT ATG ACC-3'). Positive clones were screened by amplified ribosomal DNA restriction analysis (ARDRA) as follows: the PCR products were digested with HaeIII (New England Biolabs, Germany) and separated electrophoretically on 2% MetaPhor agarose (Cambrex; distributed by Biozym, Germany). ARDRA patterns were analyzed using the Phoretix 1D software (Nonlinear Dynamics, United Kingdom). A hierarchical cluster analysis was performed applying the Jaccard similarity index and the Complete Linkage algorithm for calculating a similarity dendrogram. Representative clones from clusters comprising at least six operational taxonomic units were selected for partial sequencing. From each cluster, at least two clones were sequenced. The numbers of clones affiliated with the identified phylotypes (i.e., displaying identical ARDRA patterns) were recorded using the Phoretix database, which had been set up from the ARDRA patterns. Reamplified DNA extracted from SSCP bands was cloned as described above. For each band, 24 clones were picked and screened for the appropriate insert size and different ARDRA patterns. Representatives of the different patterns were sequenced.
16S rRNA gene sequence analysis.
The feasibility of 16S rRNA gene sequencing of flow cytometrically sorted, fixed cells was checked with pure cultures of Escherichia coli HB101 and Cupriavidus necator JMP134. Both cultures were subjected to sodium azide fixation, DAPI staining, sorting, and DNA extraction. Successful amplification and partial sequencing of the 16S rRNA genes revealed that DAPI staining of DNA did not affect the PCR by, e.g., blocking the Taq polymerase or cause mutations to an extent that affected the sequence analysis (data not shown).
Partial DNA sequencing of cloned 16S rRNA gene amplicons was performed using a BigDye RR Terminator AmpliTaq FS Kit version 3.1 (Applied Biosystems, Germany) and the sequencing primers 27F and 519R (41). Cloned amplicons from SSCP bands were sequenced using the vector-specific primers M13 uni (21) and M13 rev (29), respectively, as sequencing primers. Direct sequencing of reamplified SSCP bands was done with the primers 530f (41) and UniBac927r. Electrophoresis and data collection were carried out on an ABI PRISM 3100 Genetic Analyzer (Applied Biosystems, Germany). Data were analyzed by the ABI PRISM DNA Sequencing Analysis software, and sequences of both complementary strands were assembled by the ABI PRISM Autoassembler software. The BLASTN program (http://www.ncbi.nlm.nih.gov/BLAST; 2) was used to search for similar sequences in the GenBank, EMBL, and DDBJ nucleotide sequence databases, and the Sequence Match tool was used to search for similar sequences compiled by the Ribosomal Database Project II Release 9 (http://rdp.cme.msu.edu; 19). In order to verify the taxonomic classification resulting from the database retrieval, with particular sequences, a phylogenetic analysis was performed using the ARB package version Linux Beta 030822 (http://www.arb-home.de; 45).
The partial 16S rRNA sequences were aligned to the ARB database (release of January 2003) by the ARB editor, with manual correction. The aligned sequences were added to the ARB tree using the quick-add parsimony tool, applying phylum-specific filters.
Statistical analysis of community dynamics.
In order to obtain a more integral picture of how subcommunities changed over time, and depending on the salinity, statistical analyses were performed on the data. To estimate relative abundances for each taxon, the total cell numbers per gate as derived by flow cytometry and taxon identities obtained from clone libraries were combined. Taxon determinations using PCR-generated clone libraries are limited to the most abundant species in the libraries. Therefore, this method might underestimate species numbers in natural samples by overlooking rare species. However, within this study, general shifts in composition and frequencies determined by the detectable (and most abundant) taxa were assumed to correspond to the key degraders of diesel fuel in the examined saline soil. Their changes therefore represent the general community dynamics within the samples as reactions to the different salinities. Thus, the abundance of each taxon in a subcommunity was calculated by multiplying the proportion of the taxon within the gate by the total cell number of that gate. The total abundance of a taxon in a sample was the sum of its abundances in all different gates where it was detected, ignoring its possible occurrence outside the gates. The total cell numbers of the gates were correlated with the total cell counts of a sample, again ignoring cells outside the gates.
Species richness (number of taxa) and abundances (cell counts per taxon) were estimated. For univariate description of the samples, the Shannon-Wiener diversity index (on a log10 basis) was calculated (63). Additionally, an iterative multiregression approach was used to identify those taxa that best explained changes between the samples. For this purpose, taxa were randomly chosen, and their abundances were compared to the pooled data from all samples. As a result, those taxa that best explained the overall changes between the samples were estimated. Taxa selected by this method may not be the most abundant in all samples but those that best explained the differences between the samples.
To get a more detailed insight into community processes, a multidimensional-scaling (MDS) analysis based on taxon abundance data of the subcommunities was performed. Prior to the analysis, the data were square root transformed in order to level statistical weights between common and rare species. The appropriate transformation method was estimated by plotting the log-transformed standard deviation against the log-mean-transformed subcommunity taxon abundance data and then determining the slope of the trend line (for further explanation, see reference 18). Finally, for MDS analysis, a similarity matrix of the taxon abundance database on the Bray-Curtis similarity matrix was computed. The Bray-Curtis similarity is a widely used index in community analysis and expresses how similar two samples are to each other. It integrates both the presence and absence of species between samples and considers the proportional abundances of the species within samples (25). MDS attempts to place the samples on a (usually two-dimensional) sheet so that the rank order of the distances between samples on the sheet exactly agrees with the rank order of the matching similarities. The distance on the plot (i.e., the similarity) of two or more subcommunities indicates to what extent they share the same species and abundances (16). For community analyses, the ecological computation package PRIMER was used (17).
Nucleotide sequence accession numbers.
The partial 16S rRNA gene sequences that were determined have been deposited in the GenBank, EMBL, and DDBJ nucleotide sequence databases under accession no. AY918097 to -918118 and DQ153878 to -153948.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
|
|
Phylotypes related to an unidentified alphaproteobacterium that were also found by SSCP analysis showed similar signs of proliferation. They dominated subcommunities 6 and 9 (15% NaCl; 60.6% of all cells) and subcommunities 10 and 12 (20% NaCl; 81.2% of all cells) after 42 days of cultivation (Fig. 4). Phylotypes affiliated with Ralstonia and Halomonas were associated with these subpopulations. At the end of the cultivation, Ralstonia became dominant in the corresponding subcommunities 22 and 25 at 15% NaCl and 26 and 27 at 20% NaCl. Multiplication was also inferred from similarly shaped subcommunities in gates 7/8 and 11/13 after 42 days and 23/24 and ungated (Fig. 4)/29 after 84 days of incubation. These groups generally contributed less than 20% to the respective samples and contained mainly large cells of Ralstonia sp. The DNA contents given by channels 65 and 115 for gates 8 and 7, respectively (Fig. 4), are indicative of doubling. Members of the genus Ralstonia were abundant at all salinities and over the whole time course (present in 27 out of 30 subcommunities and dominating 20 of them) (Fig. 4). Halomonas, a genus also identified by SSCP fingerprinting (Fig. 2), contributed significantly to many subcommunities during growth under saline conditions. Also, Dietzia was present at significant abundances in dominant subcommunities (gates 5, 6, and 9), mainly after 42 days of cultivation, although SSCP analysis (Fig. 2) indicated the general presence of the genus.
The prevalence of the other phylotypes displaying prominent SSCP bands, like Naxibacter, Marinobacter, and Idiomarina, was not reflected in the 16S rRNA gene clone libraries generated from sorted subcommunities. Conversely, numerous phylotypes not detected in SSCP bands contributed significantly to the sorted subcommunities. Examples are Acidovorax, Janibacter, and phylotypes affiliated with the Ectothiorhodospiraceae and the Methylophilaceae (Fig. 4). Surprisingly, even phylotypes related to Desulfobacula (Desulfobacteraceae) and Desulfosporosinus (Peptococcaceae) were detected, which may indicate sulfate reduction.
Community analysis.
The highest taxon richness (10 taxa), according to the data obtained by analyzing the sorted subcommunities, was found in the community grown at 20% NaCl after 84 days of incubation, while the lowest number was seen in the 0% NaCl sample after 42 days (3 taxa). In all other samples, five to seven taxa were detected (Table 4). The total abundances of sorted subcommunities ranged from 0.76 x 108 (at 0% NaCl after 84 days) to 4.11 x 108 (at 15% NaCl after 42 days) cells per ml. Diversity index values revealed that the samples with the highest abundances and/or taxon numbers did not result in high diversity index scores. The samples with the highest taxon numbers (e.g., the 20% NaCl sample after 84 days) had only low abundances, while those with the highest abundances (e.g., the 15% NaCl sample after 42 days) had only low to moderate taxon numbers but were characterized by the same dominant taxa (Fig. 5). Diversity values (Table 4) ranged from a minimum of 0.32 (in the 0% NaCl sample after 42 days; six species at 1.22 x 108 cells per ml) to 0.75 (in the 7.5% NaCl sample after 84 days; seven species at 1.49 x 108 cells per ml). Low diversity values, which usually characterize samples with low species numbers and unevenly distributed, high abundances, are often an indication of stressed and poorly adapted communities, while high diversity values hint at a better adaptation of the microbial consortium. Multiregressional analysis revealed that 92.6% of the variation between the samples could be explained by the occurrence of the characteristic phylotypes affiliated with Ralstonia, Sphingomonas, the unidentified alphaproteobacterium, Desulfobacteraceae, and Peptococcaceae.
|
| DISCUSSION |
|---|
|
|
|---|
Information on the taxonomic composition of the community was obtained by SSCP fingerprinting of the PCR-amplified V4-V5 region of 16S rRNA genes performed at different cultivation times for all cultures. Flow cytometry was used to get quantitative information on the community dynamics and for the identification of active populations. Flow cytometry is a well-established method for depicting bacterial consortia in natural environments (15, 22, 55, 65, 68) and was recently combined with cell sorting as a prerequisite for phylogenetic analyses of subcommunities (49, 53). We used narrow ranges of DNA contents and cell sizes as the criteria for cell sorting. We assumed that distinct groupings in DNA-versus-cell-size plots exhibiting high DNA contents or showing signs of a cooccurrence of different cell cycle stages represented actively proliferating bacteria. This assumption was generally confirmed by the phylogenetic identification of clearly dominant populations in these groupings. It appears justified to conclude that these proliferating populations contributed to the degradation of the diesel fuel, since other carbon sources were absent. Similar proliferation patterns obtained by flow cytometry had already been successfully employed to study the dynamics of binary bacterial communities under various cultivation conditions (50, 69). Such an analysis based on DNA contents can be regarded as highly reliable and quantitative, since there is no risk of active excretion of the dye by fixed cells (56).
The overall community behavior was in agreement with commonly accepted concepts about ecological community shifts after disturbance (28, 31, 32), in that originally dominant specialists were replaced by the copiotrophic and fast-growing generalists Ralstonia and Halomonas as the leading taxa in all subcommunities during entire cultivations. This held true especially in those incubations that were disturbed by more extreme changes in salinity (i.e., 0% and 20% NaCl). The degradation efficiencies inferred from the oxygen consumption were low in these incubations. A number of additional phylotypes possibly representing specialists appeared toward the ends of these cultivations, e.g., Sphingomonas, Janibacter, members of the Ectothiorhodospiraceae and the Methylophilaceae, or putative sulfate reducers, such as species related to Desulfobacula (Desulfobacteraceae) and Desulfosporosinus (Peptococcaceae). The increase in species numbers also agrees well with common ecological models of late successional stages of well-adapted communities (28, 31, 32). Whereas species numbers in the presumably better-preadapted communities incubated at 7.5% and 15% NaCl (i.e., near natural salinities) remained stable over time, the increase in species numbers in cultures grown at 0% and 20% indicates considerable community adaptation, which also explains the delayed diesel fuel degradation. A Sphingomonas sp. appearing toward the end of the incubation at 0% NaCl and dominating the culture showed the clearest signs of proliferation in cytometry plots. It is thus likely that it was mainly the need of the population to adapt to the changed salinity or, more probably, to the utilization of diesel fuel, followed by an increase in cell number, that was responsible for the delayed but finally efficient diesel fuel degradation. Sphingomonas species are common soil bacteria known as versatile degraders of hydrophobic organic compounds, such as polycyclic aromatic hydrocarbons (35, 37, 66).
Proliferation was also observed for a dominant subcommunity mainly consisting of an alphaproteobacterium and species of Ralstonia, Dietzia, and Halomonas. The dominance of these taxa at 7.5% and 15% NaCl throughout the incubation was confirmed by SSCP analysis. The proportions of the strains changed during the incubation. Ralstonia and Dietzia species are known to degrade various hydrocarbons (13, 14, 54, 64, 72, 76), and a halotolerant Dietzia sp. has been isolated from a similar diesel fuel-degrading consortium (59). A highly uniform, proliferating subcommunity of minor abundance present in gates 7 and 8 (Fig. 4), but also visible in other samples until the end of the incubation, contained a different species of Ralstonia, as seen from its larger cells.
In addition to the identification of proliferation patterns, the development of new subcommunities could be quantitatively followed by flow cytometry and MDS. The narrow DNA distribution of the subcommunity confined by gate 28 strongly indicates that the comprised Halomonas and alphaproteobacteria belong to different species than those in gates 6 and 9, which exhibited a very different, uniform DNA pattern (Fig. 4). This subcommunity also included Janibacter, a genus of actinomycetes known to degrade several xenobiotics (33, 42, 75), and a member of the Bacillaceae related to Virgibacillus marismortui, a halophilic bacterium.
Another conspicuous subcommunity (Fig. 4) arose at the highest salinity to quantities between 1 and 2% of the whole community. The cells within this subcommunity were characterized by a very high, uniform DNA content and variable but generally large cell sizes. Unfortunately, insufficient DNA for PCR amplification was retrieved from this subcommunity. It may be that any of the marine genera Idiomarina, Marinobacter, and Alcanivorax that formed strong SSCP bands were comprised in this subcommunity. Strains of Halomonas, Marinobacter, and Alcanivorax are important degraders of hydrocarbons in marine environments (30, 38) and have been suspected to be involved in the degradation of diesel fuel components at high salinities (27). Alternatively, these genera may have been comprised in other subcommunities but were missed because the clone libraries were too small to cover all major phylotypes. It should be pointed out that our method of community identification counterselects pleomorphic species and those with cell cycle behaviors that do not favor the subpopulations with uniform DNA contents.
Janibacter, Acidovorax, and phylotypes affiliated with the Ectothiorhodospiraceae and the Methylophilaceae, as well as putative sulfate reducers, were detected only in sorted subcommunities, but not in SSCP fingerprints. Low biomass, and therefore insufficient DNA, might be the reason for the failure to detect these rare phylotypes in community fingerprints, whereas cell sorting reduces the dominance of other organisms. Oxygen depletion during rapid diesel fuel degradation might explain the occurrence of anaerobic bacteria related to Desulfobacula (Desulfobacteraceae) or Desulfosporosinus (Peptococcaceae). However, it appears too daring to conclude that anaerobic degradation of diesel fuel components (66, 74) took place in our incubations. Higher than average stability of these organisms in overall declining communities may also have favored their detection toward the end of the cultivation. Although sulfate reducers are considered obligate anaerobes, there is evidence that some species can tolerate oxygen (39, 48) or may even respire aerobically (21, 23, 36, 67).
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
| REFERENCES |
|---|
|
|
|---|