Previous Article | Next Article ![]()
Applied and Environmental Microbiology, October 2006, p. 6800-6807, Vol. 72, No. 10
0099-2240/06/$08.00+0 doi:10.1128/AEM.00971-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Södertörn University College, School of Life Sciences, SE-141 89 Huddinge, Sweden,1 Swedish University of Agricultural Sciences, Department of Microbiology, SE-750 07 Uppsala, Sweden2
Received 25 April 2006/ Accepted 11 August 2006
|
|
|---|
|
|
|---|
The classical approach for identification of viable microorganisms in environmental samples, including sediment, is plate counting on agar medium. In the Baltic water column, a large number of bacteria have been successfully isolated (20, 25, 37). However, the number of microbially diverse populations in sediments is predicted to be much higher than that in the water column (45), and the composition of viable microorganisms in Baltic Sea sediments is still largely unexplored.
Approaches for studying activity of microbial communities in sediments include measurements of net rates of biochemical processes (e.g., O2 uptake, CO2 production, ATP production, etc.). Another approach for assessing bacterial community productivity in marine environments is the measurement of microbial incorporation of radioactively labeled thymidine into newly synthesized DNA (14, 15). However, this method cannot discriminate the relative contributions of different microbial taxa to community productivity.
In recent years, there have been developments in molecular approaches for linking function or activity with identity. For example, Mar-FISH is based on combination of uptake of radioactive substrates with fluorescent in situ hybridization (FISH) (24), and stable isotope probing can be combined with molecular fingerprinting approaches to link microbial identity (biomarker) and activity (isotope assimilation) (38). Another promising molecular approach that has recently been used for identification of growing cells in environmental samples is based on incorporation of the thymidine analogue bromodeoxyuridine (BrdU) into the DNA of cells during DNA replication (3, 46). The DNA with incorporated BrdU can be selectively extracted by immunocapture and analyzed by molecular fingerprinting techniques, such as T-RFLP, to determine the composition of the growing members of the community (3, 4, 21).
In this study, we used a combination of BrdU immunocapture and T-RFLP fingerprinting to assess the composition of growing bacteria in Baltic Sea sediments. We sampled sediments outside a paper mill located south of Västervik, on the southeast coast of Sweden. This site was the focus of a large Swedish remediation initiative from 2001 to 2003. Sediment samples were collected from the inner parts of the bay in the former industrial area where the Västervik paper mill operated between 1915 and 1980. The same site was sampled twice. The first sampling was in 2002, from an area of the bay where the sediments were highly polluted with polycyclic aromatic hydrocarbons (PAHs), polychlorinated biphenyls (PCBs), and mercury (Hg), which also were reported to leak to surrounding water, air, and soil environments (34). The second sampling was performed exactly 1 year later after removal of the top highly polluted sediment layers by dredging. When dredging was completed, only low levels of pollutants could be detected in the sediments.
This study was performed to monitor the community structures and to identify actively growing bacteria in surface sediments before and after the Västervik remediation project. The same region of the Baltic Sea was studied on two consecutive years, but with two drastically different environmental conditions due to removal of the surface polluted sediment layers. A combination of molecular approaches was used to examine the actively growing bacteria in these sediments.
|
|
|---|
Chemical analysis of pollutants.
The Swedish Environmental Institute (IVL; Gothenburg, Sweden) collected samples and performed all chemical analyses prior to dredging from 1999 to 2000 and after dredging in July 2003. To quantify PCBs, sediments were pretreated with concentrated sulfuric acid and purified using an alumina column (150 mesh, chromatographic grade, combusted at 400°C for 4 h and stored at 120°C before use) and eluted with pentane for removal of matrix interference. The purified extracts were concentrated to a final volume of 1.0 ml using a Kuderna Danish evaporator. Analysis was performed by gas chromatography on a Varian model 3800 (Palo Alto, CA), with an electron capture detector, which was equipped with a 50-m-long capillary column (CP-Sil 8 CB; Chrompack, Holland). PCB components were identified and quantified using internal and external standards (certified standard solution from Larodan-Sweden, NO813). For PAH quantification, sediment samples were purified in a silicate column and eluted with pentane as described above and then reextracted to an organic state, hexane, by additions of water. The purified extracts were then concentrated to 1 ml. PAH concentrations were analyzed using a mass spectrometer, Varian Saturn 2000 Ion-Trap (Palo Alto), equipped with a 25-m-long capillary column (CP-Sil 8 CB; Chrompack, Holland). PAH components were identified and quantified using internal and external standards (certified standard solution from Larodan-Sweden, SMR 1491). The total mercury concentration was determined by gas chromatography (6, 8). Before gas chromatography, samples were incubated with nitric acid during autoclaving. After reduction of the sediment samples with tin chloride, the mercury concentration was analyzed using a Perkin-Elmer 2100 flameless atomic absorption spectrophotometer (Perkin-Elmer, Norwalk, CT). The level of detection was 0.01 mg mercury/kg dry weight (DW) sediment.
DNA extraction and BrdU immunocapture.
DNA extractions were performed from each of the three separate, replicate sediment cores taken from the four sampling stations on the two sampling occasions. DNA was extracted from 1-ml portions of the upper sediment layers in each core that were either preincubated with BrdU as described above, or without BrdU, using a FastDNASPIN kit (for soil) following the manufacturer's instructions (Q-BIOgene, Carlsbad, CA). Immunochemical purification of BrdU-labeled DNA was performed by a modification of the methods described by Urbach et al. (46) and Jansson and coworkers (3, 4, 21). Monoclonal anti-BrdU antibodies (3 µl; Roche, Basel, Switzerland) were mixed at a 1:9 ratio with herring sperm DNA (1.25 mg ml1 in phosphate-buffered saline [PBS]; Promega, Madison, WI) which had been denatured at 100°C for 5 min, quickly transferred to ice, and kept on ice for 5 min. The mixture was incubated for 1 h at room temperature. Magnetic beads (Dynabeads) coated with goat anti-mouse immunoglobulin G (DYNAL, Oslo, Norway) were washed three times with 1 mg ml1 acetylated bovine serum albumin (BSA) in PBS buffer (PBS-BSA) using a magnetic particle concentrator (DYNAL) and resuspended in PBS-BSA to the initial concentration. Sediment DNA samples (1 µg total in 25 µl of PBS) were denatured (heated for 5 min at 100°C, quickly transferred to ice, and kept on ice for 5 min), mixed with 30 µl of the herring sperm-DNA antibody mixture, and incubated for 1 h in the dark at room temperature with constant agitation. The samples were mixed with 10-µl portions of Dynabeads, and the incubation was continued for an additional 1 h. To elute the BrdU-containing DNA fraction, 100 µl of 1.7 mM BrdU (in PBS-BSA) was added, and the samples were incubated for 1 h in the dark at room temperature with constant agitation. The DNA was separated from the beads with the magnetic particle concentrator, isolated by two rounds of ethanol precipitation, and dissolved in 20 µl of sterile distilled water.
The efficacy of BrdU immunocapture was confirmed by comparing visible bands on agarose gels after DNA extraction, immunocapture, and PCR amplification with control sediment to which no BrdU was added and no PCR-amplified product was visible (data not shown).
T-RFLP.
16S rRNA genes (BrdU-labeled or nonlabeled DNA) were PCR amplified from all DNA extracts from sediment samples (see above) using the following bacterial primers, as previously described (12): fD1-Hex (AGAGTTTGATCMTGGCTCAG) (48), labeled at the 5' end with hexachlorofluorescein; and 926r (CCGTCAATTCCTTTRAGTTT) (31). The PCR products (12 µl, approximately 20 ng DNA) were digested for 2 h (50 µl total volume) using HaeIII restriction enzyme (Promega), and the fragments were separated by electrophoresis on an ABI PRISM 377 DNA sequencer (Applied Biosystems, Foster City, CA). The lengths of the fluorescent terminal restriction fragments (TRFs) were determined using GeneScan, version 3.1, analysis software (Applied Biosystems), and their relative peak areas were determined by dividing the individual TRF area by the total area of peaks within the thresholds of 60 to 490 bp. Only peaks with relative fluorescence intensities over 50 U and that contributed to the total community abundance at 1.5% or more were included in the later analyses.
Cloning.
Cloning and sequencing of 16S rRNA genes PCR amplified from BrdU-extracted DNA were performed to enable possible species identification of TRFs of interest in the T-RFLP profiles. The same PCR products from DNA extracts used for T-RFLP were used for cloning so that the results could be directly compared. The PCR products from the three replicates and four sampling sites at each sampling time were pooled since the T-RFLP results showed that there was minimal variation between these samples, resulting in two samples corresponding to DNA from actively growing community members before and after dredging, respectively. The amplified, approximately 900-bp BrdU-incubated 16S rRNA gene fragments (total of 5 ng) were cloned in separate reactions using a TOPO TA cloning kit (Invitrogen, Carlsbad, CA) following the manufacturer's instructions. Cloned inserts were screened by RFLP using the r926 and fD1 primers and amplification conditions described. The PCR products were digested with HaeIII and RsaI (Promega) in individual reactions and then pooled together. The resulting restriction fragments were separated by gel electrophoresis in 2.2% agarose gels. Gels were stained with ethidium bromide (1 mg liter1). In total, 220 clones before dredging and 210 clones after dredging were screened. The relative abundance of each individual phylotype, which was represented by a unique restriction pattern in RFLP analysis, was calculated by dividing the numbers of individual clone types by the total number of clones that were analyzed.
Isolation and identification of PAH-degrading bacteria.
One hundred sixty milliliters of frozen (80°C) polluted surface sediments that were collected in four 50-ml Falcon tubes was thawed at 4°C overnight. These sediment samples were pooled and homogenized by gentle stirring and spiked with 10 ppm phenanthrene as described by Brinch and colleagues (7). Three replicates, each containing 50 ml of the spiked sediment slurry samples, were poured into 200-ml Pyrex flasks with a tightened screw-cap top and incubated in darkness at 5°C for 7 days. After 7 days of incubation, samples (1 g) were homogenized in 10 ml PBS buffer and serially diluted in PBS before plating onto 10% tryptic soy agar media (Difco Laboratories, Detroit, MI). After inoculation, the plates were covered with a thin layer of phenanthrene by the sublimation method as previously described (1). After 12 h of incubation at 5°C in the dark, colonies and surrounding agar were assessed for color changes and clearing zones under UV light illumination indicative of phenanthrene degradation (26). Potentially positive colonies were also tested for catechol dioxygenase activity by spraying 0.5 M catechol solution over the plates and assessing eventual color change of colonies or/and agar medium, as previously described (50).
Sequencing and phylogenetic analysis.
16S rRNA genes of isolates cultivated from the sediment samples and from one clone of each RFLP group in each sample were partially sequenced on an Applied Biosystems International PRISM 377 DNA sequencer using the fD1 and r926 primers described above. Sequences were compared to sequences in the GenBank database by using the blastn software (2) and also the BLAST software from NCBI's Map Viewer to find the closest relatives. Multiple sequence alignments were performed with the MUSCLE software, version 3.5, available at http://www.drive5.com/muscle (11). Maximum likelihood trees (13) were calculated using the software PHYML version 2.4 (19). One thousand bootstrap resamplings were analyzed. Trees were edited in BBEdit, version 6.1 (Bare Bones Software, Inc., Bedford, MA), and drawn in NJplot (36). Reference 16S rRNA gene sequences derived from NCBI or from bacteria with completely sequenced genomes were downloaded from NCBI's Map Viewer. All 16S rRNA gene sequences were also tested for chimeras in Chimera Check, version 2.7, Ribosomal Database Project II (9). Chimeric 16S rRNA gene sequences that were detected within the data set were discarded. Sequences with 97% sequence similarity or more were considered to belong to the same species as the currently deposited sequences in the database. This takes into account microvariations, which may be induced by PCR and cloning (42), as well as variations between 16S rRNA gene copies, which usually differ by 1% or less (28).
Statistical analyses.
Correspondence analysis (CA) (30) was performed on the T-RFLP data resulting from three DNA extractions from each of three sediment cores taken at each sampling site. T-RFLP data from BrdU-incubated and -nonincubated samples were prepared in the same matrix and analyzed using ADE-4 software (44). TRF data were entered according to TRF size (bp) and relative abundance. TRF peaks were clustered together in the range of ±1.0 bp. The statistical significance (P < 0.001) of clustering patterns in CA plots was determined by one-way analysis of variance (ANOVA) using Prism 4 software (GraphPad Software, Inc., San Diego, CA).
Completeness of clone libraries was calculated by using coverage (C) values from the equation C = 1 (n1/N) x 100 (18), where n is the number of unique clones and N is the total number of clones examined.
Nucleotide sequence accession number.
The nucleotide sequences have been deposited in the EMBL nucleotide sequence database under accession no. AM167975 to AM167991, AM167995, and AM167997 to AM168001 for clones and AM260538 to AM260540 and AM286417 to AM286420 for isolates.
|
|
|---|
![]() View larger version (9K): [in a new window] |
FIG. 1. Pollutant concentration in DW sediment from Örserumsviken, Västervik, Sweden, before (gray bars) and after (black bars) dredging was commenced.
|
![]() View larger version (9K): [in a new window] |
FIG. 2. Ordination diagram from CA of T-RFLP data derived from surface sediment (0 to 2 cm) of triplicate sediment cores collected from Örserumsviken before and after dredging. Numbers correspond to sampling sites. Color labeling: red, actively growing bacterial community before dredging; blue, actively growing bacterial community after dredging; green, total bacterial community before dredging; black, total bacterial community after dredging. Distances between the profiles that clustered together were significant (P < 0.001). Eigenvalues for first and second ordination axes were 37% and 32%, respectively.
|
![]() View larger version (41K): [in a new window] |
FIG. 4. Phylogenetic classification of sequenced BrdU-incorporated16S rRNA genes. Clones sequenced in this study have the abbreviation VAS printed in plain and bold text, representing samples taken before and after dredging, respectively. Alignments of approximately 900 characters were used to generate the tree. The phenogram was drawn by the maximum likelihood method, and bootstrap values are presented as percentage of 1,000 resamplings. "U" indicates the bootstrap value is lower than 50%. The number of clones of individual sequences in each clone library is indicated inside parentheses. The predicted length of HaeIII-generated TRFs, based on in silico restriction of clone sequences, is indicated to the right of the brackets. Uncl., uncultured.
|
![]() View larger version (28K): [in a new window] |
FIG. 3. Clone frequency of actively growing bacterial community members from sediment collected from Öreserumsviken before and after dredging. Color labeling of bacterial genus or class: red, Spirochaeta; black, Actinobacteria; green, Cyanobacteria; gray, Alphaproteobacteria; yellow, Deltaproteobacteria; blue, Gammaproteobacteria; orange, Flavobacteria.
|
|
View this table: [in a new window] |
TABLE 1. Relative abundances of 16S rRNA gene sequences of actively growing bacteria represented in T-RFLP profiles and clone libraries
|
-proteobacterium, accession no. AY539822 (99% sequence similarity [5/70 colonies]); Exiguobacterium oxidotolerans, accession no. AB105164 (100% sequence similarity [14/70 colonies]), and Arthrobacter sp., accession no. AJ967024 (99% sequence similarity [12/70 colonies]). The pseudomonad and the gammaproteobacterium were also identified in the clone library of the active community in the polluted sediments (clones VAS1 and VAS7, each contributing with 5/70 and 7/70 colonies, respectively), indicating in these cases good matches between the cultivation-based and molecular approaches for assessing growing cells. However, no representatives of Exiguobacterium oxidotolerans or Arthrobacter sp. were present in the clone library. The remaining five isolates that were isolated from phenanthrene-incubated agar plates showed no catechol dioxygenase activity. Their 16S rRNA gene sequences had the following matches: Serratia sp., accession no. AY689059 (98% sequence similarity [4% of relative colony abundance]); a
-proteobacterium, accession no. DQ226202 (96% sequence similarity [10% of relative colony abundance]); Pseudomonas putida, accession no. DQ178233 (98% sequence similarity [7% of relative colony abundance]); Stenotrophomonas maltophilia, accession no. AY581129 (96% sequence similarity [7% of relative colony abundance]), and Arctic Sea Ice bacterium, accession no. AF468434 (97% sequence similarity [7% of relative colony abundance]).
Community analyses.
The numbers of TRFs detected in T-RFLP electropherograms were used to gauge the relative richness of the bacterial communities, although one should keep in mind that only dominant members of the community are detected using this approach and that there is some bias in the methodology used. Before dredging, the numbers of TRFs were 15 and 60 for actively growing and for total community members, respectively. After dredging, the numbers of TRFs were 14 and 64 for actively growing and for total community members, respectively. There were therefore considerably fewer TRFs represented in the growing community compared to the total community at both sampling periods. This low number of growing bacteria was also reflected by comparison to the low numbers of unique sequences in clone libraries created from BrdU-labeled 16S rRNA genes, with only 13 and 10 unique sequences obtained at the first and second sampling years, respectively. Hence, approximately the same numbers of actively growing bacteria were detected using the two molecular approaches, T-RFLP and cloning and sequencing.
|
|
|---|
Based on results from clone libraries and T-RFLP analysis of 16S rRNA gene fragments, we found that the number of actively growing bacterial communities was approximately 75% lower than for the dominant members of the total communities. Also, none of the TRFs in the BrdU-extracted DNA were present in electropherograms of directly extracted DNA. The low number of 16S rRNA genes in the clone libraries was striking, since screening of approximately 200 clones per library (Fig. 4) was sufficient to reach complete coverage for each library (i.e., no new 16S rRNA gene sequences could be detected). One possible explanation for the low numbers could be insufficient BrdU incorporation by all growing cells in the community. In previous studies, it was found that some anaerobic bacteria did not incorporate the similar molecule, thymidine, when it was added exogenously during DNA replication (17, 22, 49). However, in our study, representatives of the Deltaproteobacteria with close phylogenetic matches to anaerobes were capable of BrdU uptake (Fig. 4). Also, in this study, we probably selected for the fast-growing bacterial community members since sediments were incubated with BrdU for only 10 h. However, this was done to limit the ability of the bacteria to use BrdU as a carbon source, as previously observed to be problematic for radioactively labeled thymidine during production and growth rate measurements (40). In a separate study, we confirmed that BrdU was not being used as a carbon source in these sediments (Edlund and Jansson, unpublished). Interestingly, we recently found that there were very few growing bacteria detected in soil using the BrdU immunocapture approach (21), reinforcing the hypothesis that most bacteria in the environment persist in a dormant, nongrowing physiological state.
We observed significant differences in the compositions of both the actively growing communities and the total/dominant members of the communities at the two different sampling periods, and there are several possible reasons for these observations. For example, the difference could be due to removal of pollutants that repressed some bacterial populations (5, 16) or that conversely served as carbon sources for specific populations (10, 29), or the difference could be due to the selection of resistant (e.g., mercury resistant) populations before dredging. Another explanation for the significant differences in microbial community structures before and after dredging could be due to succession since the upper sediment layers were completely removed. Therefore, the community profiles representing sediment samples after dredging may reflect a new niche that was recolonized by immigration and recruitment (39) of Baltic bacteria from the water column and the surrounding environment.
Sequences corresponding to Deltaproteobacteria were only obtained in the clone library from the polluted sediments, supporting previous findings that this group is ubiquitous in organic-rich marine sediments (23). Two sequences, VAS2 and VAS5, belonging to the purple sulfur bacteria (Chromatiaceae), known to require anoxic conditions, hydrogen sulfide, and illumination, were obtained from sediments before dredging, suggesting that that these photosynthetic sulfide oxidizers were also actively growing in the polluted sediments. However, these species were probably metabolizing sulfur and carbohydrate resources that were stored inside their cells, as previously described by Stal and Moezelaar (43), during the BrdU incubation period since the samples were incubated in darkness.
Members of the Gammaproteobacteria were growing in the sediments on both sampling occasions (Fig. 3). However, only one operational taxonomic unit (OTU), with closest identity to a 16S rRNA gene sequence with NCBI accession no. AY539822, dominated the community after dredging while four OTUs were present before dredging (Fig. 4). In previous studies concerning PAH-polluted environments, members of the Gammaproteobacteria were also shown to be metabolically active and capable of degrading low-molecular-weight PAHs (27, 41, 47). In line with these findings, we identified a Pseudomonas sp. (accession no. AF521651) and an unnamed species belonging to Gammaproteobacteria (accession no. AY539822) by cloning DNA from cells growing in the polluted sediments. We also cultivated a matching Pseudomonas sp. and an isolate of Gammaproteobacteria that were capable of phenanthrene transformation. Another isolate in this study, Exiguobacterium oxidotolerans (accession no. AB105164), was able to transform phenanthrene but was not represented in the clone library. However, we previously isolated a similar strain that had 99% 16S rRNA sequence identity to this bacterium in sediments from a different region of the Baltic Sea along a eutrophication gradient (12), indicating that this species may be common in coastal Baltic Sea surface sediments.
Several of the 16S rRNA gene sequences from the clone libraries that could be matched to specific TRFs had similar relative abundance values in T-RFLP profiles and in clone libraries (Table 1), demonstrating the utility of this combination of culture-independent approaches. The T-RFLP data enabled us to determine which of the phylogroups were most dominant, whereas the sequence information enabled us to ascertain their identities. Thus, a combination of approaches is useful to obtain more information about the identities of key microbial species in complex environmental samples, rather than relying on one technique alone. This study may serve as a good base for future studies of key bacterial species involved in recycling nutrients and in degradation of pollutants in marine sediments.
This project was financed by The Baltic Sea Foundation.
Published ahead of print on 1 September 2006. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»