Previous Article | Next Article ![]()
Applied and Environmental Microbiology, December 2008, p. 7321-7328, Vol. 74, No. 23
0099-2240/08/$08.00+0 doi:10.1128/AEM.01266-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Rich Boden,
Hélène Moussard,
Hendrik Schäfer,* and
J. Colin Murrell
Department of Biological Sciences, University of Warwick, Coventry CV4 7AL, United Kingdom
Received 7 June 2008/ Accepted 1 October 2008
|
|
|---|
|
|
|---|
DNA SIP recently identified Methylophaga-like organisms as active methylotrophs that assimilated methanol and methylamine in surface waters of the English Channel (36). That study also demonstrated that 16S rRNA gene sequences representing clades of uncultivated Gammaproteobacteria were also retrieved from the heavy DNA for each of these compounds that clustered close to Methylophaga. A DNA SIP experiment with methanol substrate dilution to concentrations anticipated to reflect those in situ (34) confirmed the involvement of Methylophaga spp. in methanol consumption and retrieved functional genes involved in methanol metabolism from these active methylotrophs using metagenomic libraries.
The goal of the current study was to extend our previous observations that were made under nonbloom conditions by studying methylotrophic populations in the context of a phytoplankton bloom dominated by Emiliania huxleyi and Karenia mikimotoi (formerly Gyrodinium aureolum). Both the coccolithophores (e.g., Emiliania) and the small dinoflagellates (e.g., Karenia) are associated with dimethylsulfoniopropionate production (22, 29), and phytoplankton blooms are known to produce relevant C1 compounds or their precursors, including methanol (13), methylated sulfur compounds (24), and methyl halides (42). As with our previous marine SIP studies (34, 36), seawater samples were incubated with methanol and methylamine, and in this investigation, SIP incubations were also carried out with 13C-labeled methane, dimethylamine, methyl bromide, and DMS to identify microbial populations that are actively involved in the cycling of these C1 compounds during phytoplankton blooms in situ.
|
|
|---|
Incubation with 13C-labeled substrates.
Samples taken from the edge of the bloom were chosen to set up SIP incubations with several 13C-labeled C1 substrates. Seawater sample aliquots of 750 ml were added to 1-liter serum bottles with the addition of 0.1% (750 µl) marine ammonium mineral salts medium (modified from that described in reference 12) and substrate. A total of 75 µmol of 13C-labeled substrate was added to bottles for methanol, monomethylamine, dimethylamine, methyl bromide, and methane (final concentration of 100 µM, assuming complete dissolution). For DMS, 187.5 µmol of substrate was added to make up a final concentration of 250 µM. All serum bottles were crimp sealed with butyl rubber bungs to prevent the loss of volatile substrates. All 13C-labeled compounds had a purity of 99% or higher and were obtained from Cambridge Isotope Laboratories (Hook Hampshire, United Kingdom), except methylated amines (Sigma, Gillingham, United Kingdom). [13C2]DMS was prepared by a method adapted from that described for labeled dimethyl sulfoxide synthesis (5). Sodium sulfide nonahydrate (6.5 g) was dissolved in 6.5 ml of sterile deionized water in a glass test tube and cooled to 0°C in an ice-water bath, with vigorous stirring. Subsequently, 5 g of [13C]methyl iodide (Cambridge Isotope Laboratories LTD, Andover, MA) was added dropwise over a period of 30 min prior to incubating the reaction mixture at 0°C for 5 h with stirring. Five milliliters each of 2 M sodium hydroxide solution and 1 M sodium thiosulfate solution was added, the reaction mixture vessel was then connected to a receiving tube held at –170°C in liquid nitrogen, and the reaction mixture was allowed to warm to 40°C in a water bath. The [13C2]DMS was distilled from the reaction mixture for 90 min and then redistilled into a sterile receiving vessel for 1 h. Sterile deionized water was added to the receiving vessel to dissolve the [13C2]DMS, and the resulting solution was transferred, with washings, to a sterile 1-liter serum vial which was then sealed with a butyl rubber bung. The concentration and purity of the [13C2]DMS solution were assessed by gas chromatography with a flame ionization detector. A total of 250 ml of a 7 mM solution of pure [13C2]DMS was obtained.
For all substrates, parallel incubations were set up as 12C-unlabeled controls and 13C- and 12C-labeled substrate incubations were harvested at a single time point. With the exception of monomethylamine and dimethylamine, substrate utilization was monitored by gas chromatography using a flame ionization detector. Measurement of DMS and methyl bromide concentrations in sterile seawater controls confirmed that the degradation observed for SIP incubations was due to biological processes and not to chemical decomposition. The concentrations of the methylamines were assumed to mirror those of methanol; recovery of [13C]DNA from methylamine and dimethylamine incubations confirmed that methylated amines had been assimilated. Following substrate depletion, SIP incubations were filtered through 0.22-µm Sterivex filters and frozen at –80°C until processed for nucleic acid extraction.
DNA extraction, SIP gradient centrifugation, and fractionation.
Total nucleic acids were extracted directly from Sterivex filters according to a previously described protocol (36). Briefly, lysozyme, proteinase K, and sodium dodecyl sulfate were used to lyse cells, and lysates were transferred to 15-ml phase lock tubes (Qiagen, West Sussex, United Kingdom) for phenol-chloroform and chloroform extractions. Purified DNA was quantified on a 1% (wt/vol) agarose gel. Aliquots (1 to 5 µg) of DNA extracts from each of the SIP incubations were added to cesium chloride (CsCl) solution (average density of
1.725 g ml–1) and transferred to an ultracentrifuge gradient tube for centrifugation and fractionation as previously described (37). Briefly, tubes were added to a Vti 65.2 rotor (Beckman Coulter, Fullerton, CA) and centrifuged at 44,100 rpm (177,000 x gav [average]) for 40 h at 20°C. Gradients were fractionated from bottom (fraction 1, highest density) to top (fraction 12, lowest density) into 425-µl fractions. DNA was purified from CsCl and quantitatively recovered by precipitation with glycogen (20 µg) and polyethylene glycol (30% PEG 6000 and 1.6 M NaCl). Purified DNA was suspended in 30 µl of sterile LoTE buffer (3 mM Tris [pH 8], 0.2 mM EDTA), and 5-µl aliquots were run on a 1% (wt/vol) agarose gel for quantification and to identify the distribution of [13C]DNA relative to background, unlabeled [12C]DNA (see Fig. S1 in the supplemental material). These data indicated that the 13C labeling of DNA was very high for methanol, monomethylamine, and DMS incubations; most of the DNA for these 13C-incubated samples eluted in heavy fractions (fractions 7 and 8,
1.725 to 1.730 g ml–1). The detection of 13C-labeled DNA confirmed that the substrate was incorporated into microbial biomass. For dimethylamine and methyl bromide incubations, the extent of DNA labeling was less pronounced. For methyl bromide incubations, there was almost no difference between the smears of DNA across gradients associated with 12C- and 13C-methyl bromide SIP incubations (see Fig. S1 in the supplemental material).
DGGE.
For each fractionated gradient, two fractions were selected for analysis of the "heavy" [13C]DNA (fractions 7 and 8,
1.725 to 1.730 g ml–1), and one fraction was selected for the characterization of "light" [12C]DNA (either fraction 11 or 12,
1.710 to 1.705 g ml–1). One-microliter aliquots of gradient fractions were used as the template for PCR to obtain 16S rRNA gene fragments suitable for denaturing gradient gel electrophoresis (DGGE) analysis. Each 50-µl reaction mixture consisted of 25 pmol of primers 341f-GC and 534r (32), 1x (NH4)2SO4 buffer (Fermentas, York, United Kingdom), 1.5 mM MgCl, 33.6 µg nonacetylated bovine serum albumin (Sigma, Gillingham, United Kingdom), 40 nmol of each deoxynucleoside triphosphate, 1.25 U Taq polymerase (Fermentas). The reaction mixture tubes were loaded directly into the block at 95°C (simplified hot start), followed by an initial denaturation at 95°C for 5 min and 30 cycles of 94°C for 1 min (denaturing), 55°C for 1 min (annealing), and 72°C for 1 min (extension). A final extension at 72°C for 7 min was followed by a holding step at 10°C. Five-microliter aliquots were quantified on a 1% (wt/vol) agarose gel.
For DGGE, 5-µl aliquots (100 to 300 ng) were run on a 10% polyacrylamide gel with a 30 to 70% denaturing gradient (100% denaturant is 7.0 M urea and 40% deionized formamide) according to the D-Code system instructions (Bio-Rad, Hercules, CA). Gels were run overnight (14 h) at 85 V and then stained for 1 h in SYBR Green I (Invitrogen, Paisley, United Kingdom). Gel images were captured with a FLA-5000 imaging system (Fujifilm, Tokyo, Japan). Bands selected for sequence analyses were sampled from the gel by means of sterile pipette tips and amplified from the gel, using the PCR conditions described above for DGGE. Sequencing was done with the 341f primer and a BigDye terminator version 3.1 kit (Applied Biosystems, Foster City), and the sequencing products were run on an ABI Prism 3130 genetic analyzer (Applied Biosystems) by the Molecular Biology Service, University of Warwick, Warwick, United Kingdom. DGGE band sequences were approximately 150 bases in length.
For determining the relatedness of the DGGE fingerprints, gels were imported into Gelcompar II (Applied Maths, Sint-Martens-Latem, Belgium) and normalized to ladder bands and additional internal standard bands. A dendrogram produced with the unweighted-pair group method using average linkages was generated by performing a Pearson correlation with background-subtracted densitometric curves, which takes band intensities into account. The output of the clustering analysis was independent of the input order of DGGE fingerprints.
16S rRNA gene libraries.
Clone libraries of bacterial 16S rRNA genes were generated from the original seawater samples (outside of bloom, edge of bloom, inside of bloom; 36 clones sequenced from each) and for the heavy DNA associated with the five substrates (methanol, monomethylamine, dimethylamine, methyl bromide, and DMS) that yielded 13C-labeled DNA (24 clones sequenced from each). The PCR to amplify the 16S rRNA gene used primers 27f and 1492r (25) and the same amplification reaction as that used for DGGE, except with an extension time of 1.5 min. Products were cloned into the TOPO-TA vector according to the manufacturer's protocol (Invitrogen). Screening was done as described previously (38), and cloned 16S rRNA gene inserts were sequenced at the Edinburgh node of the NERC Molecular Genetics Facility, using the 27f primer. Pintail (3) software was used to identify suspected chimeras and identified one heavy-band sequence which was likely chimeric in origin and several water library sequences that were likely chimeric; these sequences were excluded from further analyses. For seawater samples, classification of 16S rRNA gene sequences was done using an RDP-II classifier (53) after manually verifying base calls. For the 16S rRNA gene libraries constructed using [13C]DNA from SIP experiments, manually verified 16S rRNA gene sequences were compared to those in GenBank (6) to retrieve the three closest matches for each library sequence. Sequences were aligned within ARB (30) software, and an alignment was exported to MEGA4 (49). Evolutionary distances were computed using the maximum-composite likelihood method (50) and are presented in the units of the number of base substitutions per site. All positions containing gaps and missing data were eliminated from the data set (complete deletion option). There was a total of 466 nucleotide positions in the final data set. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1,000 replicates) are shown next to the branches (9). The tree was drawn to scale, with branch lengths given in the same units as those of the evolutionary distances used to infer the phylogenetic tree. In the absence of cultivated methylotrophic organisms that fell within the groups of 16S rRNA genes derived from SIP experiments, clades were defined based on the consistent association with particular substrates, but a specific cutoff value was not used.
Nucleotide sequence accession numbers.
All sequences were deposited in GenBank for the marine samples taken from the edge (accession no. EU399242 to EU399272), the inside (accession no. EU399273 to EU399306), and the outside (accession no. EU399307to EU399340) of the phytoplankton bloom. The 16S rRNA gene clone library sequences from heavy DNA were deposited with the following accession numbers for SIP incubations with dimethylamine (accession no. EU399341 to EU399364), DMS (accession no. EU399365 to EU399386), methyl bromide (accession no. EU399387 to EU399407), monomethylamine (accession no. EU399408 to EU399428), and methanol (accession no. EU399429 to EU399451). DGGE band sequences from heavy DNA were deposited with the following accession numbers for SIP incubations with methanol (accession no. EU399452 to EU399457), monomethylamine (accession no. EU399458 to EU399464), dimethylamine (accession no. EU399465 to EU399469), DMS (accession no. EU399470 to EU399474), and methyl bromide (accession no. EU399475 to EU399477).
|
|
|---|
![]() View larger version (75K): [in a new window] |
FIG. 1. Physical location and bacterial composition of samples. (A) A grayscale Aqua-MODIS satellite radiance image generated from a mixed-wavelength composite image (551-nm, 488-nm, and 443-nm wavebands) on 26 July 2006. Emiliania huxleyi coccolithophore reflectance is readily visible in this image. A chlorophyll a satellite image from the same day shows chlorophyll abundance corresponding to the Karenia bloom adjacent to the bloom on the north side of the English Channel (image not shown). Satellite data were received and processed by the NERC Earth Observation Data Acquisition and Analysis Service (NEODAAS) at Dundee University and Plymouth Marine Laboratory (www.neodaas.ac.uk). (B) Bacterial DGGE fingerprints of samples taken from the locations indicated in panel A. (C) Frequency of 16S rRNA gene clones belonging to major phylogenetic groups across the different gene libraries analyzed.
|
0.63% in the headspace) for seawater incubations with methane were unchanged for several months (data not shown), and these incubations were not analyzed further.
16S rRNA gene fingerprinting of DNA from SIP experiments.
DGGE was used to profile the bacterial communities associated ["heavy"] and "light" fractions for both the 12C-labeled control samples and the 13C-labeled incubated samples. The 12C-labeled incubated samples showed no profile differences between the "heavy" and the "light" fractions (data not shown), whereas unique fingerprints were evident for all C1 substrate incubations (Fig. 2). As expected, the fingerprints for the light fractions of all tubes clustered together. The "background" bacterial communities in each SIP incubation were more similar to one another than to the "heavy" 13C-labeled fraction fingerprints of the same incubation. However, for the methyl bromide SIP incubation, the fingerprint of the "heavy" DNA was less clearly unique for the light DNA than for the other substrate incubations, reflecting the possibility that only a small amount of DNA was labeled and was just detectable above background [12C]DNA. All other DGGE fingerprints from heavy fractions (fractions 7 and 8) of 13C-labeled substrate incubations clustered in distinct clades apart from the "light" DNA, with monomethylamine and methanol fingerprints clustering closely, with some similarity to the dimethylamine fingerprints. Dimethylsulfide [13C]DNA fingerprints were distinct from all other patterns in this study, reflecting the unique composition of active methylotrophs enriched in these SIP incubations. Individual bands from fingerprints representing "heavy" DNA from 13C1 incubations were selected for PCR reamplification and sequencing. These sequences were used to assign band sequences to specific SIP-related 16S rRNA gene clades derived from this study and from a previous SIP study that was carried out under nonbloom conditions in the English Channel (36). The results indicate that Methylophaga spp. were associated with methyl bromide, methanol, and methylamine SIP incubations, whereas additional clades were affiliated with dimethylamine and DMS, likely contributing to their more distinct fingerprint profiles (Fig. 2).
![]() View larger version (86K): [in a new window] |
FIG. 2. DGGE fingerprint comparison of "light" and "heavy" DNA associated with DNA SIP incubations with different carbon sources (MMA, monomethylamine; DMA, dimethylamine; MOH, methanol; MBr, methyl bromide; and DMS). The dendrogram scale bar refers to the percent similarity of Pearson correlations between fingerprint densitometric curves. Shading of the triangle pointers indicates the phylogenetic affiliations of sequenced bands, most of which were associated with clades in Fig. 3. Numbers at the bottom right of the fingerprints correspond to sequenced bands submitted to GenBank. For example, the open triangle for fraction 7 of the [13C]methylbromide SIP (MBr_7) is labeled MBr_7_3 for the GenBank submission. Several bands were not associated with clades but were affiliated with sequences in Fig. 3, as follows: MBr_7_3 is identical to MBr_587_7; DMA_7_6 is closest to MBr_587_24; DMS_7_2 is identical to DMS_584_3.
|
Methanol-assimilating phylotypes.
Phylogenetic analysis demonstrated that sequences associated with SIP incubations with methanol clustered in the Methylophaga species clade of the Gammaproteobacteria together with sequences from a previous SIP incubation (36) and with several characterized Methylophaga isolates (Fig. 3 and see Fig. S2 in the supplemental material). In addition, one cloned 16S rRNA gene sequence obtained from the methanol SIP was a member of a clade of unknown phylogenetic affiliation, which also contained one DMS and DMA-SIP-derived cloned 16S rRNA gene. The methanol clade identified in a previous methanol SIP experiment (Fig. 3; see Fig. S2 in the supplemental material) has the closest affiliation to Gammaproteobacteria sequences in GenBank that were retrieved from multiple marine Arctic surface sediments or detected on the surface of submerged artificial substrates incubated in marine water near China. This clade, however, was not detected in the current study.
![]() View larger version (59K): [in a new window] |
FIG. 3. Phylogenetic affiliations of 16S rRNA gene sequences obtained by 13C-labeled SIP incubations with methanol (MOH), monomethylamine (MMA), dimethylamine (DMA), methyl bromide (MBr), and DMS SIP incubations. Selected GenBank sequences from uncultivated clones and reference strains are included for comparison. Bootstrap values are included for all branch points on this neighbor-joining tree. GenBank accession numbers are shown in parentheses. The scale bar (in the tree) represents 5% sequence divergence. The collapsed clades are expanded in Fig. S2, S3, and S4 in the supplemental material. The division-level affiliation of sequences is indicated in the boxes at the right. CFB, Cytophaga-Flavobacterium-Bacteroides.
|
Phylotypes assimilating DMS.
Almost all 16S rRNA gene sequences derived from the DMS SIP "heavy" DNA were nearly identical and formed an additional clade with low relative diversity (Fig. 3; see Fig. S4 in the supplemental material). The DMS clade was most closely related to the methanol SIP clade associated with the class Gammaproteobacteria and identified in a previous study (35); it shared close similarity (96%) with sequences retrieved from clone libraries associated with DMS-enriched seawater samples from the Sargasso Sea (52) and was approximately 91% similar to the Methylophaga sp. clade, based on the percent similarity between sequences of DMS_584_22 and Methylophaga marina (GenBank accession no. X95459) over 722 bases. Another sequence from the "heavy" DNA of the DMS SIP was affiliated with a clade of unknown phylogeny.
|
|
|---|
This study represents a comprehensive survey of active methylotrophs in a marine surface water sample during a bloom of phytoplankton associated with the production of dimethylsulfoniopropionate. The methylotrophs detected in this survey are consistent with the results of our pilot study with only methanol and monomethylamine under nonbloom conditions obtained a year prior to the current sampling event (36); however, the use of a wider range of C1 substrates allowed the identification of a larger diversity of methylotrophs than that found previously, including populations that assimilate dimethylamine, DMS, and methyl bromide. The DMS SIP clones obtained were most closely related to clones obtained from DMS enrichments from Pensacola, FL, and the Sargasso Sea by Vila-Costa and colleagues (52), suggesting that the latter had similar metabolic activities and indeed represented DMS-degrading populations. Those sequences were classified as "uncultivated Methylophaga"; however, given the relatively low similarity of the 16S rRNA gene sequences of these cloned 16S rRNA gene sequences to those of Methylophaga isolates (around 92%) and their distinct clustering, supported by bootstrap analysis (see Fig. S2 and S4 in the supplemental material), it is also possible that these represent DMS-degrading populations belonging to a different genus. Conversely, none of the DMS SIP clones was closely related to previously isolated DMS-degrading Methylophaga isolates (44), which belonged to the Methylophaga clade detected with methanol, monomethylamine, and methyl bromide. This strongly suggests that populations closely related to the isolated strains may have a preference for other C1 substrates and/or are outcompeted by those represented by the DMS clade under the specific incubation conditions. The methyl bromide SIP sequences indicate that methyl bromide may be used by members of the genus Methylophaga and an organism with a 16S rRNA gene sequence most similar to that of Phaeobacter gallaeciensis (formerly Roseobacter) within the family Rhodobacteraceae. The notion that the Phaeobacter-related population degraded methyl halides would be supported by previous cultivation-based identification of marine methyl halide-degrading organisms which were closely related (43, 45, 46); however, screening of several Methylophaga isolates has failed to show their ability to degrade methyl halides (H. Schäfer, unpublished). The observation of Methylophaga-like sequences in the [13C]methyl bromide incubation could therefore be due to the slow hydrolytic conversion of methyl bromide to methanol (1) and subsequent utilization of the resulting methanol by these organisms. If Methylophaga populations in the methyl bromide incubations became labeled with 13C due to the uptake of methanol produced by conversion of methyl halides to methanol, this would further underline their ability to take up methanol at ambient concentrations and support using SIP incubations with elevated substrate concentrations to investigate substrate-responsive populations in seawater. Together, these data suggest that marine waters harbor a diverse suite of active methylotrophs that, apart from Methylophaga spp., have been unnoticed by previous cultivation studies (8, 17, 23, 43-45) and are almost completely without representation in marine clone libraries. The sequences represented here are important targets for directed cultivation and focused activity-based studies of marine methylotrophy.
Given the focus of past marine metagenomic studies on abundant community members, it is perhaps not surprising that few genes (phylogenetic or "functional") have reflected the predominance of methylotrophic bacteria. Although formaldehyde oxidation genes were identified in the Sargasso Sea metagenomic libraries (51), genes for methane, methylamine, and methanol oxidation were not detected (18). Furthermore, the only presumed methylotroph 16S rRNA gene sequences identified in a marine metagenomic library was from Methylophilus spp., and these sequences occurred at
0.4% of the total 16S rRNA gene data set from the global ocean survey (41). The contribution of Methylophilus to marine C1 cycling remains unclear, and Methylophilus spp. have not been detected in [13C]DNA from the incubations carried out in this study. One possibility is that Methylophilus spp. represent K-selected organisms, which adapted to concentrations of carbon and nutrients that are lower than those used in this study. Cultivation-based approaches (17), enrichment cultures (44, 52), and SIP incubations (current study; 34, 36) have all demonstrated that Methylophaga spp. and related Gammaproteobacteria from multiple disparate marine samples (including estuary sediment; unpublished data) are present in the seawater samples and rapidly respond to the presence of C1 substrates. It is possible that these organisms may represent low-abundance and r-selected bacteria capable of opportunistic growth in the presence of relatively high concentrations of growth substrates during phytoplankton blooms, for example.
This study represents a comprehensive cultivation-independent survey of active marine methylotrophs and demonstrates that previously unrecognized bacterial groups are present in seawater, which are capable of responding to the presence of added C1 substrates. The presence of numerous clades of presumed substrate-specific methylotrophs presents a challenge to microbiologists to focus cultivation and quantitative molecular approaches to better understand the metabolism and distribution dynamics of these organisms with potentially enormous biogeochemical significance.
This work was supported by funding from the Natural Environment Research Council (United Kingdom) Aquatic Microbial Metagenomics and Biogeochemical Cycles grant, reference NE/C001 923/1. H.S. was supported by a NERC postdoctoral fellowship (NE/B501404/1). J.D.N. acknowledges support from the Natural Sciences and Engineering Research Council (Canada).
Published ahead of print on 10 October 2008. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
Present address: Department of Biology, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»