Previous Article | Next Article ![]()
Applied and Environmental Microbiology, March 2005, p. 1405-1416, Vol. 71, No. 3
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.3.1405-1416.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Department of Marine Sciences, University of Georgia, Athens, Georgia,1 Savannah River Ecology Laboratory, University of Georgia, Aiken, South Carolina,2 Departamento de Microbiologia y Biologia Celular, Facultad de Farmacia, Universidad de La Laguna, La Laguna, Tenerife, Spain3
Received 2 August 2004/ Accepted 3 October 2004
|
|
|---|
-Proteobacteria); Caulobacter (
-Proteobacteria); and Brachymonas and Xenophilus (ß-Proteobacteria) in the May-02 sample and to Ketogulonicigenium (Roseobacter lineage,
-Proteobacteria) and novel
-Proteobacteria in the Sept-02 sample. Our study suggests that diverse bacterioplankton participate in the metabolism of DMSP in coastal marine systems and that their relative importance varies temporally. |
|
|---|
Previous studies have made significant progress in identifying which bacterial taxa mediate DMSP degradation in natural marine communities. Culturing approaches have identified DMSP-utilizing bacteria in the Roseobacter lineage (
-Proteobacteria) and in the genera Alcaligenes (ß-Proteobacteria) and Pseudomonas (
-Proteobacteria) (1, 5, 6, 12, 21, 22, 38). Culture-independent studies of surface waters from estuaries and open ocean sites with high DMSP turnover rates have revealed DMSP-utilizing bacteria in the Bacteroidetes, Roseobacter, ß-Proteobacteria, and
-Proteobacteria lineages (1, 13, 16, 24, 35, 39, 40). Most of these field studies, however, have targeted a broad taxonomic grouping (at the phylum or class level), and few have resolution at the level of genus and below.
Flow cytometry (FCM) discriminates between populations within complex bacterial communities based on fluorescence and size differences among the cells. In conjunction with a sorting unit, defined FCM populations can be physically separated and then subjected to further taxonomic analysis (2, 3, 23, 28, 32, 36). However, the combined methods are challenging to use, primarily because of the low bacterioplankton densities typical of seawater and the large number of sorted cells needed for subsequent DNA-based analysis.
Coastal salt marshes have one of the highest natural DMS emission rates and represent an active DMSP-degrading ecosystem (4). In this study, a culture-independent approach was employed to identify bacteria involved in DMSP degradation at the species level. We amended water from a salt marsh tidal creek with DMSP and identified cells responding to the addition by FCM sorting of cells with elevated nucleic acid (NA) content. Methods were optimized to obtain DNA template from as few as 25,000 FCM sorted preserved bacterial cells. Terminal restriction fragment length polymorphism (T-RFLP) analysis and sequencing of 16S rRNA genes in the sorted populations were used to identify bacterial taxa involved in DMSP degradation.
|
|
|---|
FCM.
Flow-cytometric analysis was performed with a MoFlo flow cytometer (DakoCytomation, Glostrup, Denmark). Prior to running on the instrument, preserved bacterioplankton were stained with Sybr Green II (1:5,000 dilution of the commercial stock; Molecular Probes Inc.) in the dark at room temperature for 20 min (18) and then mixed with an internal bead standard (1-µm-diameter yellow-green fluorescent beads; Fluoresbrite YG Microspheres; Polysciences, Warrington, Pa.). A sterilized phosphate-buffered saline (PBS) solution (Puraflo; Dakocytomation, Fort Collins, Colo.) served as the sheath fluid. Data acquisition was triggered by green fluorescence (FL1). All signals were collected with logarithmic amplification.
All FCM populations were gated initially on the FCM cytogram of one DMSP replicate after 24 h of incubation, and the gating frames were kept fixed through all the measurements within a single sample date. Gate notation was based on the value of FL1, which is related to the cell NA content. The FCM populations that exhibited high FL1 values were designated the "high-NA-content groups" (HNAs), and the others were designated the "low-NA-content groups" (LNAs). The FCM populations observed for each sample are shown in Fig. 1. Because the FCM population gates chosen were specific to each date, the gate notation refers to different ranges of fluorescence and side scatter for each sample date.
![]() View larger version (38K): [in a new window] |
FIG. 1. Flow-cytometric analysis of control (left panel) and DMSP-amended (right panel) samples collected after 24 h of incubation in May 2002, September 2002, January 2003, and May 2003. Beads and HNA and LNA cells are delimited by gates on each cytogram. Boldface letters are used to label subpopulations that were analyzed with T-RFLP. Asterisks are used to label FCM populations used for clone library construction.
|
DNA extraction and PCR amplifications.
Sorted cells from each FCM population were filtered onto a 25-mm-diameter, 0.2-µm-pore-size polycarbonate filter (Poretics Products) and washed three times with 2 ml of PBS followed by 2 ml of deionized water. The filter was air dried, and its edge was trimmed off. One quarter of the trimmed filter was stored at 20°C, while the other three quarters were cut into pieces and transferred into an 0.5-ml tube. Ten microliters of Lyse-N-Go PCR reagent (Pierce Inc., Rockford, Ill.) was added to the tube to lyse the bacterial cells on the filter according to the manufacturer's protocol. A second lysis step was conducted by adding 82 µl of autoclaved water to the same tube and subjecting the tube to three freeze-thaw cycles. The resulting lysate served as the DNA template for PCR amplification.
16S rRNA gene amplification was carried out with Ready-To-Go PCR beads (Amersham Pharmacia, Piscataway, N.J.) with 0.4 µM concentrations of both forward (8F, 5'-AGAGTTTGATCCTGGCTCAG-3') and reverse (1492R, 5'-TACGGYTACCTTGTTACGACTT-3') primers. A touchdown PCR program was performed with annealing temperature sequentially decreasing from 62 to 52°C by 1°C per cycle, followed by 15 cycles at 52°C. In each cycle, denaturing (at 95°C), annealing (at 62 to 52°C), and extension (at 72°C) steps were of 1-min duration. An initial 3-min denaturation and final 10-min extension step were also included in the PCR program. PCR amplification was confirmed by electrophoresis on ethidium bromide-stained 1% agarose gels. The amplicons were excised from the gel and cleaned with the QIAquick gel extraction kit (Qiagen, Valencia, Calif.) before subsequent T-RFLP and sequencing analysis.
T-RFLP analysis.
Direct PCR amplification of 16S rRNA genes from sorted cells with a 6-carboxyfluorescein (FAM)-labeled forward primer was inconsistent, probably due to the difficulty in binding the FAM-labeled primer to the low-copy-number templates. Instead, 16S rRNA genes of sorted cells were amplified first with nonlabeled primers 8F and 1492R as described above and then reamplified with 0.4 µM FAM-labeled 8F and nonlabeled 785R (5'-CTACCAGGGTATCTAATCC-3') primers with two Ready-To-Go PCR beads in 50-µl-volume reaction mixtures with a single annealing temperature of 59°C. PCR amplification was confirmed by electrophoresis on ethidium bromide-stained 1% agarose gels. The reamplified amplicons were excised from the gel, cleaned with the QIAquick gel extraction kit, and digested with CfoI restriction enzyme (Roche, Indianapolis, Ind.) at 37°C for 3 h, after which an ethanol precipitation was performed. The restricted amplicons were resuspended in 12 µl of deionized formamide plus 0.7 µl of DNA fragment length standard (GeneScan 2500 TAMRA; Applied Biosystems, Warrington, United Kingdom). The terminal restriction fragment (T-RF) lengths were determined on an ABI Prism 310 genetic analyzer (Applied Biosystems, Foster City, Calif.) in GeneScan mode.
T-RFLP data were obtained successfully for most sorted bacterial FCM populations. The exceptions were as follows: (i) one set of the May-02 triplicates was used up in the process of optimizing the conditions for DNA extraction and PCR amplification, (ii) some FCM populations were too small to provide enough template for PCR amplification of 16S rRNA genes (Fig. 1), and (iii) one FCM population (Sept-02, LNA1) failed to amplify even though cell numbers were relatively high.
Identification of bacteria associated with responsive T-RFs.
16S rRNA genes directly amplified from sorted cells were used to construct clone libraries of 16S rRNA genes for selected samples (May-02 and Sept-02). Libraries were constructed using a TA cloning kit (Invitrogen Corporation, Carlsbad, Calif.) with the pCR 2.1 vector according to the manufacturer's instructions.
Sequencing was carried out on an ABI PRISM 310 genetic analyzer in sequencing mode, with the BigDye terminator cycle-sequencing kit (PE Biosystems, Foster City, Calif.). 8F served as the sequencing primer. The average length of usable sequence was 400 bp. Bacterial identities were determined using the GenBank blastn search (http://www.ncbi.nlm.nih.gov/BLAST/) and the RDP-II sequence match program (http://rdp.cme.msu.edu/html/analyses.html).
The bacteria associated with T-RFs showing significant differences between control and DMSP-amended microcosms (responsive T-RFs) were putatively assigned using the following process. Initially, 20 clones were randomly selected from each library for sequencing, and the length of the T-RF for each clone was calculated based on sequence data with TRFLPtools, a Visual Basic program for Microsoft Excel (13, 34). If any responsive T-RFs remained unidentified, a clone pool approach was used to rapidly screen the 16S rRNA gene library. Briefly, five random clones were pooled prior to plasmid extraction with the QIAprep Spin Miniprep kit (Qiagen Inc.). T-RFLP analysis was carried out by the method described previously with 500 ng of extracted plasmid DNAs as templates for the FAM-labeled PCR amplification. If any matching T-RFs were found, each of the five clones was analyzed individually to identify the specific clone(s) of interest in the pooled sample. This process continued until the identities of all T-RFs were assigned or more than 200 clones were screened. Then, clones putatively responsible for responsive T-RFs were coinjected with community amplicons in a single T-RFLP analysis. Finally, the inserts of the clones of interest were sequenced to determine taxonomic identity.
T-RFLP analysis of small populations.
We determined whether apparent bacterial community structures were significantly affected by the number of cells used in the T-RFLP analysis, since sorted FCM populations differed in cell number by as much as 25-fold (Table 1). The total number of cells in the May-02 water sample was determined by epifluorescence microscopy counting after 4',6'-diamidino-2-phenylindole (DAPI) staining. Aliquots containing various numbers of cells (2 x 109, 2 x 107, 2 x 106, 2 x 105, and 2 x 103) were filtered onto 25-mm-diameter, 0.2-µm-pore-size polycarbonate filters. PCR and T-RFLP analysis of the 16S rRNA genes were then performed according to the procedure described above.
|
View this table: [in a new window] |
TABLE 1. Cell number and percentage of the total cell count for HNA and LNA subpopulations in four salt marsh bacterioplankton samples after 24 h of incubationa
|
A hierarchical cluster analysis was performed using Primer v5 (Primer-E Ltd, Plymouth, United Kingdom) to make quantitative comparisons of T-RFLP profiles among bacterial FCM populations. Before analysis, T-RFLP output data were standardized as described previously (34), except that relative peak area instead of peak height was used as a proxy for the relative abundance of bacterial taxa associated with each T-RF. T-RFs with lengths greater than 600 bp or relative peak areas less than 2% of total area were excluded from the analysis.
Nucleotide sequence accession numbers.
The GenBank accession numbers for 16S rRNA gene sequences determined in this study are AF547399, AF547405 to AF547409, AF547413 to AF547415, AF547417, AF547418, AF547421, AF547422, AF547424, AF547425, AF547427, AY476738 to AY476741, AY476749 to AY476751, AY476763, AY476764, AY476767, AY476769 to AY476771, AY476774 to AY476776, AY476778 to AY476781, AY476784, AY476785, AY476794, AY476795, AY476797, AY476799 to AY476802, AY476804, and AY476805.
|
|
|---|
![]() View larger version (18K): [in a new window] |
FIG. 2. Total bacterial cell numbers (means ± standard deviations) counted by FCM in control and DMSP treatments over the 48-h incubation for each sample date.
|
T-RFLP and sequence analysis of sorted populations.
T-RFLP analyses of the replicate subsamples for a given sample date were highly similar as revealed by cluster analysis (Fig. 3). Thus, there was low variability among the replicate incubations within microcosms and consistency in results of molecular analysis. The cluster analysis of T-RFLP data also grouped all samples based on collection date, regardless of the sorting gates used. For the May-02 and Sept-02 samples (those in which a bacterial cell NA-size response to DMSP amendments was evident), FCM populations of control and DMSP-amended samples clustered separately with high dissimilarity values. However, for Jan-03 and May-03 samples (those in which little or no bacterial cell NA-size response to DMSP additions was evident), populations from control and DMSP-amended microcosms were intermingled.
![]() View larger version (33K): [in a new window] |
FIG. 3. Cluster analysis of T-RFLP data after 24 h of incubation. Similarity calculation was based on the Bray-Curtis dissimilarity on square-root-transformed relative peak area. FCM populations from triplicate incubations are shown, except for the May-02 sample, for which only duplicates were available. Sample codes indicate treatment (D, DMSP amendments; C, controls), replicates (1, 2, or 3), and sorted subpopulation (LNA2, LNA3, or HNA).
|
-Proteobacteria (mainly Roseobacter lineage), ß-Proteobacteria,
-Proteobacteria, Bacteroidetes, and Actinobacteria (Table 2). The seventh T-RF (length, 115 bp; T-115) was unidentified. Bacteria associated with T-115 (unidentified), T-202 (ß-Proteobacteria), and T-350
360 (
-Proteobacteria) were stimulated by DMSP amendments as evidenced by greater relative abundance in the elevated-NA-larger-size populations. Bacteria associated with T-32 (Actinobacteria), T-57 (
-Proteobacteria, Roseobacter lineage), T-98 (Bacteroidetes), and T-565 (
-Proteobacteria) decreased in relative abundance in the higher bacterial cell NA-size class after DMSP amendments. In the Sept-02 experiment (Fig. 4b), there were six responsive T-RFs. Bacteria represented by T-231 (
-Proteobacteria) and T-518 (Roseobacter lineage,
-Proteobacteria) were relatively enriched in FCM populations with elevated NA content after DMSP amendments, while T-57 (Roseobacter lineage,
-Proteobacteria), T-98 (Bacteroidetes), T-202 (ß-Proteobacteria), and T-310 (
-Proteobacteria) reacted in the opposite manner.
![]() ![]() View larger version (69K): [in a new window] |
FIG. 4. T-RFLP chromatogram of May-02 (a) and Sept-02 (b) samples. Empirically determined taxonomic identities of T-RFs and their lengths are indicated beneath the chromatogram. The average contribution of each identified T-RF is indicated by percentages. Alpha, -Proteobacteria; beta, ß-Proteobacteria; gamma, -Proteobacteria; unid, unidentified. Sample codes are as in Fig. 3.
|
|
View this table: [in a new window] |
TABLE 2. Taxonomic identification of bacterial groups associated with peaks in T-RFLP profiles based on 16S rRNA gene clone library analysisa
|
T-RFLP analysis of small populations.
Following optimization of the DNA extraction method, PCR amplification and T-RFLP analysis were carried out on bacterial population sizes as low as 2,000 cells from unpreserved samples and 20,000 cells from preserved samples. Bacterial populations containing 2 x 109, 2 x 106, and 2 x 103 unpreserved cells from the May-02 water sample were analyzed by T-RFLP. Regardless of the initial cell number, T-RFLP chromatograms typically contained the same major peaks in approximately the same relative abundances (Fig. 5). Based on this finding, FCM populations differing in size by as much as 500,000 cells were compared in the T-RFLP analyses.
![]() View larger version (19K): [in a new window] |
FIG. 5. Community structures as revealed by T-RFLP chromatograms from different-sized bacterioplankton populations of 2 x 109 (A), 2 x 106 (B), and 2 x 103 (C) cells collected in May 2002.
|
|
|
|---|
The taxa identified as responsive to DMSP amendments probably include those assimilating DMSP directly as well as those utilizing products of DMSP degradation. Taxonomic groups assimilating DMSP or its degradation products but maintaining the same NA content and size would be overlooked with our approach (as was the case for cells in the Jan-03 and May-03 samples). The high concentrations of DMSP in amended microcosms (20 µM) relative to natural concentrations (20 to 100 nM) (17) and the preincubation used to establish carbon-limited conditions might have introduced bias. Nonetheless, this approach allowed an initial identification of bacteria that respond to DMSP and that may be important under natural conditions.
T-RFLP analysis was used in this study because it provides an efficient way to view the community composition of many sorted FCM populations. However, like fluorescent in situ hybridization and denaturing gradient gel electrophoresis, this quick examination of community-level 16S rRNA heterogeneity includes a tradeoff with taxonomic resolution (7). To partially overcome this limitation, we also constructed 16S rRNA gene clone libraries for selected samples to access the within-T-RF variability (Table 3). The sequence similarity of clones within dominant T-RFs was highly variable and ranged from as low as 74.1% (for T-98, representing the Bacteroidetes group in May-02) to as high as 99.3% (for T-57, representing the Roseobacter lineage in Sept-02). High variability within Bacteroidetes T-RFs generated with commonly used restriction enzymes has been found previously (7).
|
View this table: [in a new window] |
TABLE 3. Contribution of major bacterial taxa to the total bacterioplankton assemblage in control microcosms after 24 h of incubation based on percentage of total area under the T-RFLP chromatogram
|
-Proteobacteria) is ubiquitous in marine environments and is a major taxonomic group in many marine bacterioplankton communities (10). Recently, evidence from both lab and field studies has suggested that members of the Roseobacter lineage carry out DMSP transformations in marine systems (12, 22, 37). Roseobacter group isolates can degrade DMSP by two distinct pathways (cleavage and demethylation-demethiolation [1, 12, 15]) and may be important in regulating DMS formation in the surface ocean (31). During blooms of DMSP-producing algae in the open ocean, members of the Roseobacter lineage have been found in high abundance (13, 39). More recently, they have been directly linked to DMSP utilization by microautoradiography-fluorescent in situ hybridization methodology (24, 35). Several members of Roseobacter lineages identified in the present study as responding positively to DMSP additions had high 16S rRNA gene sequence similarity to previously identified DMSP-utilizing bacteria, including Ruegeria atlantica and Sulfitobacter pontiacus (1, 11-13) (Table 1). However, the response of other members of the Roseobacter lineage indicated that they either were uninvolved with DMSP metabolism or were outcompeted by other taxa (Table 2; Fig. 6). This intragroup difference demonstrates the challenge of assigning biogeochemical functions to large and diverse marine bacterioplankton groups and clearly would be missed by the use of analysis methods with lower taxonomic resolution.
![]() View larger version (50K): [in a new window] |
FIG. 6. Phylogenetic tree based on partial sequences of 16S rRNA genes showing the relationship among clones from FCM populations. Clones from the HNA populations formed after DMSP addition are labeled with boldface. Clones from LNA populations in the DSMP and control treatments are shown in lightface. The tree was constructed using Jukes-Cantor distance and neighbor-joining methods excluding positions with <50% sequence conservation. Bacillus subtilis was used as the outgroup. Bootstrap values higher than 50% are indicated at the branch nodes. The scale bar indicates the amount of genetic change in terms of the number of nucleotide substitutions per site. The GenBank accession numbers of the reference sequences are shown in parentheses. Empirically determined T-RFLP fragment lengths are shown in brackets. Only one representative is shown if more than one clone had identical sequences in one sample.
|
-Proteobacteria were readily isolated from estuarine water and surface sediment in a coastal marsh (1). Also,
-Proteobacteria have been found to be abundant in the bacterioplankton community associated with algal blooms in the ocean (13, 40) and can assimilate DMSP in natural bacterial assemblages (35). In this study, bacteria related to Brachymonas and Xenophilus (ß-Proteobacteria) in the May-02 sample and a novel
-proteobacterium in the Sept-02 sample were enriched in the FCM populations with elevated NA content after DMSP addition, indicating active linkage to DMSP utilization. In contrast, bacteria related to Actinobacteria, Flexibacter and Flavobacterium (Bacteroidetes), and Neptunomonas (
-Proteobacteria) in the May-02 sample and Cytophaga (Bacteroidetes), Brachymonas, Rhodoferax (ß-Proteobacteria), and Pseudoalteromonas and Colwellia (
-Proteobacteria) in the Sept-02 sample were concentrated in the FCM populations with unchanged or low NA content (LNA1 and LNA2) (Table 2; Fig. 4b). However, because not all DMSP-utilizing bacteria responded with an increase in bacterial cell NA content or size (Fig. 2), we cannot assume that these taxa play no role in DMSP degradation.
Cluster analysis grouped bacterial FCM populations into four clusters coincident with sampling dates, indicating that temporal factors had a major influence on bacterial community composition. The bacterial assemblages in the control microcosms for the two samples in which an NA-size response to DMSP additions was evident (May-02 and Sept-02) were similar, with relatively low contributions from ß-Proteobacteria (Table 3). The bacterial assemblages in the control microcosms of the other two samples (with no NA-size response to DMSP; Jan-03 and May-03) were likewise similar to one another but distinct from the May-02 and Sept-02 samples, with a higher abundance of ß-Proteobacteria and few
-Proteobacteria or Bacteroidetes (Table 3). These differences in bacterial community composition may underlie the disparity in cell parameter response to DMSP amendments observed among samples, and they underscore the influence of temporal variability in the composition and physiology of the bacterioplankton that transform DMSP in coastal environments.
This study was supported by NSF grant 0084164 to the Sapelo Island Microbial Observatory (SIMO).
|
|
|---|
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»