Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2005, p. 6885-6899, Vol. 71, No. 11
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.11.6885-6899.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Emma Noyes,2,
Marina G. Kalyuzhnaya,1
Mary E. Lidstrom,1,3 and
Ludmila Chistoserdova1*
Department of Chemical Engineering,1 Department of Microbiology, University of Washington, Seattle, Washington 98195,3 Omak High School, Omak, Washington 988412
Received 14 April 2005/ Accepted 13 June 2005
|
|
|---|
|
|
|---|
The recently developed cultivation-independent approaches allow identification of microbial communities active in particular metabolic processes. One of these approaches involves detection and analysis of 16S rRNA or mRNA transcribed from functional genes (i.e., nifH, pmoA, merR, and denitrification genes) via reverse transcription followed by PCR (RT-PCR), cloning, and sequencing (10, 32, 35, 45). Another approach is the use of stable-isotope-labeled substrates (stable-isotope probing [SIP]) (37). One of the most commonly used SIP techniques relies on incorporation of 13C originating from a specific substrate into nucleic acids (DNA and RNA) of microorganisms actively involved in its utilization. The 13C-labeled (heavy) nucleic acids (DNA or RNA) are separated from the 12C-nucleic acids by isopycnic centrifugation and used as templates in PCR amplification to determine the identities of the respective microbes. This approach has been successfully applied to identify members of microbial communities active in utilization of methane, methanol, propionate, and aromatic compounds (11, 17, 27, 28, 37, 38). While it is understood that SIP experiments should be conducted in conditions as closely resembling the conditions in situ as possible in order to obtain a true picture of the activities and processes native to a given environment, the success of SIP inherently depends on altering the in situ conditions to a degree. A higher-than-in-situ concentration of the substrate in question is required for successful labeling (37, 38), and prolonged incubations often are required (37, 38), which may result in enrichment of DNAs of fast-growing microbes. While possessing a potential of uncovering unexpected and unsuspected participants in certain environments (11, 37, 38), SIP results need critical interpretation.
The top layer of the sediment of Lake Washington is an environment in which steep gradients of methane and oxygen are present, and aerobic methane oxidation occurs at high rates (25). The methanotroph population in Lake Washington has been characterized in previous work (2, 6, 7). However, much less is known about metabolism of C1 substrates that are less reduced than methane. While recent tests employing primers targeting genes of the H4MPT-linked formaldehyde oxidation pathway revealed the presence of a diverse population of microbes that do not fall into known methanotrophic groups (18, 19), nothing is known about the respective roles and activities of these organisms in cycling of C1 compounds. The aim of this study was to identify active members of the bacterial community in the sediment of Lake Washington, with emphasis on C1 utilizers, via two complementary culture-independent approaches, analysis of mRNA and rRNA by using RT-PCR-based techniques and DNA-SIP with a variety of C1-labeled substrates.
|
|
|---|
Stable-isotope probing.
Microcosms were
set up, consisting of 15 ml of sediment slurry mixed with an equal
volume of 0.2 µm-pore-filtered Lake Washington water placed
into sterile 50-ml screw-cap plastic tubes and supplemented with one of
following (final concentrations): [13C]methanol (0.05%,
vol/vol), [13C]methylamine (20 mM),
[13C]formaldehyde (1 mM), or [13C]formate
(10 mM). All substrates were 99 atom% 13C
and were purchased from Sigma-Aldrich, with the exception of
[13C]methanol, which was provided by the National Stable
Isotope Resource at Los Alamos National Laboratory. Microcosms were
incubated at room temperature (20 to 25°C) with agitation (150
rpm) for up to 10 days. The incubation temperature was chosen based on
preliminary experiments in which [13C]DNA fraction
complexities for microcosms incubated at in situ temperature
(8°C) and microcosms incubated at higher temperatures were
compared by restriction fragment length polymorphism (RFLP) analysis.
The DNA complexities of the samples were similar, with microcosms
incubated at higher temperatures accumulating the [13C]DNA
fraction more rapidly (data not
shown).
Nucleic acid extraction.
Total nucleic
acids (RNA and DNA) were isolated from the sediment samples and from
each microcosm according to the protocol described by Griffiths et al.
(14), with some
modifications. Briefly, 0.5 g (wet weight) of sample and
0.5 g of 0.1-mm zirconia-silica beads (Biospec Products) were
suspended in 750 µl of extraction buffer (a mixture of equal
volumes of 10% CTAB [cetyltrimethylammonium bromide] in 1.6 M NaCl and
0.2 M phosphate buffer, pH 8.0), to which 75 µl each of 10%
sodium dodecyl sulfate and 10% lauroyl sarcosine were added. After
addition of 750 µl of phenol-chloroform-isoamyl alcohol
(25:24:1), the mixtures were homogenized in a minibeater
(Biospec Products) for 30 to 60 s at 4°C
(75% of the maximum power) and thencentrifuged at
16,000 x g for 5 min at 4°C. The aqueous phase
was mixed withan equal volume of chloroform-isoamyl
alcohol (24:1) and centrifuged at 16,000 x
g for 5 min at 4°C. Nucleic acids were precipitated
for 1 to 2 h at room temperature by adding either
MgCl2 (final concentration, 2 mM), 0.1 volume of 3 M sodium
acetate and 0.7 volume of isopropanol, or 2 volumes of 30% PEG
6000-1.6 M NaCl. Nucleic acids were recovered by centrifugation
at 18,000 x g for 10 min at 4°C, washed with
75% ethanol, and resuspended in 50 µl of sterile nuclease-free
water. Nucleic acids isolated from four to five replicates were pooled.
Nucleic acids were checked for quality and quantity by electrophoresis
in agarose gels containing ethidium bromide. DNAs isolated from the
microcosm experiments were purified from the coextracted RNA by
incubating preparations for 1 h at 37°C with 1 U of
DNase-free RNase (QIAGEN). RNAs isolated from sediment samples were
purified from the coextracted DNA by incubating twice for 30 min at
37°C with 4 U of DNase I (Ambion). Total RNA was diluted 50
times for use in
RT-PCRs.
Ultracentrifugation and DNA recovery.
DNA extracted
from the microcosms was prepared for CsCl-ethidium bromide density
gradient ultracentrifugation as described previously
(37) and centrifuged at
265,000 x g (Beckman VTi 65 rotor) for 16 h
at 20°C. [13C]DNA and [12C]DNA fractions
from each microcosm were collected using 19-gauge needles. DNA from
each fraction was purified following standard procedures
(40) and used in a second
CsCl-ethidium bromide density gradient ultracentrifugation, as
described above. After centrifugation, discrete fractions of the
gradients (<500 µl) were collected and processed as
described above. DNA samples were stored at
20°C.
PCR amplifications.
Aliquots (1
µl) of the [13C]DNA or [12C]DNA isolated
from each microcosm were used as templates in PCRs employing either 16S
rRNA gene-specific (8F/1407R
[34,
40]) or
fae-specific
(18) primers. PCR
mixtures (final volume, 20 µl) contained: 1x buffer,
1.5 mM MgCl2, 0.2 µM of each deoxynucleoside
triphosphate, 0.2 µM of the primers, and 0.5 U of DNA
polymerase (Invitrogen). Cycling was performed as described previously
(18,
33). Thirty cycles were
used to obtain 16S rRNA gene fragments (approximately
1.4 kb). fae fragments were amplified in two
steps, as described by Kalyuzhnaya et al.
(18), to increase
amplification specificity. After 30 cycles as described above, 1
µl of the amplification mixture was used as a template in a new
PCR mixture, as described above, and 25 additional cycles were
performed, to amplify a product of approximately 300
bp.
RT-PCR amplifications.
Aliquots (2 µl) of RNA
isolated from the sediment were used as templates in reverse
transcription reactions followed by amplification by PCR, according to
the manufacturer's instructions (QIAGEN One-Step RT-PCR kit). The
primer sets 8F/1407R, pmo189/mb661
(7), fae1F/fae1R, and
fae2F/fae2R were used to detect, respectively, eubacterial 16S rRNA,
pmoA, and fae gene transcripts. In addition, a
group-specific pair of primers was used to detect transcripts from a
divergent fae clone previously detected in Lake Washington,
clone L1N9 (18): primer
L1N9fae-f
(5'-ATGGAGGCACCCATGGCAG-3') and
primer L1N9fae-r
(5'-CTCAGACGCCTTCGACACC-3').
Twenty-five cycles were used to amplify 16S rRNA gene products, and 35
cycles were used to amplify pmoA products. fae
products were obtained using a two-step amplification protocol
essentially as described above, except that 35 cycles were used in the
first step, followed by 20 additional
cycles.
Clone library construction.
Amplicons of
the expected size were purified using the QIAGEN gel extraction kit and
cloned into the pCR2.1 Topo TA cloning vector, according to the
manufacturer's instructions.
RFLP analysis and sequencing.
PCR products obtained from
single-colony PCR or using purified plasmids as templates were
categorized into operational taxonomic units (OTUs), based on their
restriction patterns obtained by digestion with the following
restriction enzymes: RsaI for 16S rRNA and pmoA genes and
either BclI/SacII/NcoI/HincII/BstUI or HaeIII for fae genes.
DNA fragments were resolved in 2 to 2.5% agarose gels. The sampling
effort in each library was evaluated by calculating the coverage
(C) (13)
according to the equation C= 1
(n/N), where n is the number of
OTUs containing unique sequences and N is the number of clones
analyzed in the library. For 16S rRNA gene libraries, OTUs represented
by at least two clones and at least 10% of the OTUs represented by a
single clone were analyzed by sequencing using the primer 515-537R
(5'-CCGTMTTACCGCGCTGCTGGCA-3') or
the M13F primer
(5'-GTAAAACGACGGCCAG-3')
and the Big Dye v3.1 sequencing kit (Applied Biosystems). Reaction
analysis was performed using an ABI3790XL high-throughput capillary DNA
analyzer (Applied Biosystems) by the Department of Biochemistry DNA
sequencing facility at the University of Washington. To sequence
representatives of OTUs from pmoA and fae libraries,
the M13F primer was used. Sequences of the 16S rRNA gene (
400
bp), fae, and pmoA were aligned using the FastAligner
v3.0 program (ARB software package
[http://www.arb-home.de])
against their closest relatives in the GenBank databases as determined
using BLASTN and TBLASTX searches
(http://www.ncbi.nih.nlm.edu/BLAST).
Similarity matrixes were constructed, and threshold values of 97% for
16S rRNA gene sequences
(41) and 94% for
pmoA and fae sequences
(23,
42) were used to group
the sequences into phylotypes. One representative of each phylotype was
completely sequenced using primers M13F and M13R
(5'-CAGGAAACAGCTATGAC-3'), with
the exception of a few phylotypes related to
Beta-/Gammaproteobacteria that were partial
sequences. Affiliations with phylogenetic groups mentioned were defined
according to the Ribosomal Database Project (RDP) Classifier
(http://rdp.cme.msu.edu/classifier.jsp).
Phylogenetic analyses.
Nearly-full-length 16S rRNA gene
sequences were submitted to CHECK-CHIMERA, available on the Ribosomal
Database Project release 8.1
(http://rdp8.cme.msu.edu/html/analyses.html),
in order to identify chimeras. Phylogenetic analyses were performed
using the ARB software package
(http://www.arb-home.de).
The 16S rRNA gene phylogenetic analyses were performed by the
maximum-likelihood method
(39), using 1,285 to
1,392 nucleotide positions. The functional genes were translated into
amino acid sequences, and these were included in phylogenetic analyses
using the neighbor-joining method
(39) (Dayhoff
PAM model). Ninety-two positions for Fae and 130 positions
for PmoA were used. The lists of sequences included in the phylogenetic
analyses are available upon
request.
Nucleotide sequence accession numbers.
The 16S rRNA gene,
fae, and pmoA sequences obtained in this study were
deposited in GenBank under accession numbers
DQ066940
to
DQ067043,
DQ067044
to
DQ067063,
and
DQ067064to
DQ067087,
respectively.
|
|
|---|
1.4 kb,
300 bp, and
500 bp, respectively) were obtained and
cloned to generate specific gene libraries. Libraries of 151 clones
containing 16S rRNA gene fragments, 98 clones containing fae
gene fragments, and 85 clones containing pmoA gene fragments
were generated, and these were each screened by RFLP. Based on RFLP
patterns, the 151 clones of the 16S rRNA library were assigned to 91
OTUs, the 98 clones of the fae library were assigned to 8
OTUs, and the 85 clones of the pmoA library were assigned to
19 OTUs. Partial sequences were determined for each OTU identified in
the 16S rRNA gene library, and these were submitted to the RDP
Classifier using default parameters (confidence threshold, 80%) to
determine phylogenetic affiliations of the respective bacteria. Of the
91 OTUs in the 16S rRNA gene library, 52 (83 clones) were assigned to
the phylum Proteobacteria, 16 (27 clones) were assigned to the
phylum Bacteroidetes, 2 (4 clones) were assigned to the phylum
Nitrospira, 1 (3 clones) was assigned to the phylum
Actinobacteria, and 1 (3 clones) was assigned to the candidate
division WS3 (9). Nineteen
OTUs (31 clones) were categorized as "unclassified
bacteria" by the RDP Classifier. Nearly complete 16S rRNA gene
sequences (
1.3 to 1.4 kb) were determined for OTUs that
contained at least two clones and for at least 10% of the OTUs that
contained a single clone (a total of 50 sequences), and these were
included into a similarity matrix to sort the sequences into
phylotypes. A total of 40 phylotypes were identified (Table
1; Fig. 1). Of these, 22 were assigned to the phylum Proteobacteria (Fig. 1A, B, and C). These were
further classified as follows: 9 phylotypes (29 clones) were identified
as Betaproteobacteria and were assigned to uncultured
organisms within the families Commamonadaceae,
Oxalobacteraceae, Nitrosomonadaceae, and
Methylophilaceae; 10 phylotypes (25 clones) were identified as
Gammaproteobacteria, of which 6 phylotypes (16 clones) were
related to methanotrophic bacteria of the genera
Methylobacter and Methylomonas, while the
remaining phylotypes were assigned to organisms belonging to the genera
Achromatium and Thioploca/Beggiatoa and to
uncultured Gammaproteobacteria; and 3 phylotypes (5 clones)
were identified as Deltaproteobacteria and were related to
uncultured organisms of the order Myxococcales and
Syntrophus. Seven phylotypes (22 clones) were assigned to the
phylum Bacteroidetes, within which some were related to
uncultured organisms of the order Bacteroidales of the
families of Sphingobacteriaceae, Flexibacteraceae,
and Saprospiraceae (Table
1; Fig.
1D). The remaining
phylotypes were assigned to uncultured representatives of the
phyla Actinobacteria, Gemmatimonadetes,
Acidobacteria, candidate division WS3,
Verrumicrobium, and Spirochaetes (Table
1; Fig.
1D and
E). |
View this table: [in a new window] |
TABLE 1. 16S rRNA gene phylotype distribution in libraries obtained from total RNA isolated from Lake Washington sediment and from [13C]DNA and [12C]DNA fractions obtained from the microcosm experimentsa
|
![]() ![]() ![]() ![]() ![]() View larger version (142K): [in a new window] |
FIG.1. Phylogenetic trees obtained by maximum-likelihood analysis, reflecting the relationships of 16S rRNA gene sequences amplified from RNA isolated from Lake Washington sediment (*) and from [13C]DNA fractions (). Sequences determined in this study are in boldface. The scale bars indicate the number of expected nucleic acid substitutions per site per unit of branch length. (A) Beta- and Gammaproteobacteria phylogenetic tree. Alphaproteobacteria were used to root the tree; 1,389 characters were used to infer the tree. Partial sequences (<600 bp) were added to the tree by using a maximum-parsimony option within ARB; 433 characters were used. (B) Alphaproteobacteria phylogenetic tree. Deltaproteobacteria were used to root the tree; 1,365 characters were used to infer the tree. (C) Deltaproteobacteria phylogenetic tree; 1,392 characters were used to infer the tree. (D) Acidobacteria, Nitrospirae, Nitrospina, Actinobacteria, Gemmatimonadetes, and Bacteroidetes phylogenetic tree. Deionoccus and Thermus were used to root the tree; 1,287 characters were used to infer the tree. (E) Verrumicrobium, Planctomycetes, Chlorobi, Spirochaetes, and Chloroflexi phylogenetic tree. Deionoccus and Thermus were used to root the tree; 1,285 characters were used to infer the tree.
|
|
View this table: [in a new window] |
TABLE 2. Fae
phylotype distribution in libraries obtained from total RNA isolated
from Lake Washington sediment and from [13C]DNA
fractions obtained from the microcosm experiments
|
![]() View larger version (25K): [in a new window] |
FIG. 2. Phylogenetic tree reflecting the relationships of fae gene sequences amplified from RNA isolated from Lake Washington sediment (*) and from [13C]DNA fractions (). The tree topology was obtained from inferred amino acid sequences (92 positions) by using neighbor-joining analysis. Sequences determined in this study are in boldface. The scale bar indicates the number of expected amino acid substitutions per site per unit of branch length.
|
![]() View larger version (26K): [in a new window] |
FIG. 3. Phylogenetic tree reflecting the relationships of pmoA gene sequences amplified from RNA isolated from Lake Washington sediment. The tree topology was obtained from inferred amino acid sequences (130 positions) by using neighbor-joining analysis. Sequences determined in this study are in boldface. The scale bar indicates the number of expected amino acid substitutions per site per unit of branch length.
|
Diversity of bacterial 16S rRNA in [13C]DNA versus [12C]DNA fractions isolated from microcosms.
Two 16S rRNA gene
libraries were constructed from each microcosm, using, respectively,
the [13C]DNA fraction and the [12C]DNA fraction
as templates for PCR. The libraries were analyzed by RFLP and
categorized into OTUs as described above. The 45 clones of the library
originating from the [13C]DNA fraction of the methanol
microcosm were assigned to 10 OTUs, 3 of which were represented by 32
clones. The 43 clones of the library originating from the
[13C]DNA fraction of the methylamine microcosm were assigned
to four OTUs, of which two were represented by 41 clones. The 88 clones
of the library originating from the [13C]DNA fraction of the
formaldehyde microcosm were assigned to 31 OTUs, of which 4 were
represented by 52 clones. The 41 clones of the library originating from
the [13C]DNA fraction of the formate microcosm were assigned
to 29 OTUs. Coverage values were calculated as described in Materials
and Methods. For methanol and methylamine microcosm libraries the
values were above 90%, indicating that the sampling effort was likely
sufficient to describe the 16S rRNA gene diversity represented in these
libraries. Lower coverage values were obtained for formaldehyde and
formate microcosms (50% and 75%, respectively), indicating that these
libraries were only partially sampled.
The 16S rRNA gene diversity in the libraries originating from [13C]DNA was compared to the diversity in libraries originating from the [12C]DNA microcosm fractions. The 48 clones of the library originating from the [12C]DNA fraction of the methanol microcosm were assigned to 31 OTUs, the 45 clones of the library originating from the [12C]DNA fraction of the methylamine microcosm were assigned to 38 OTUs, the 94 clones of the library originating from the [12C]DNA fraction of the formaldehyde microcosm were assigned to 43 OTUs, and the 46 clones of the library originating from the [12C]DNA fraction of the formate microcosm were assigned to 30 OTUs. The coverage values calculated for these libraries (29% to 69%) were lower than the values for the [13C]DNA-based libraries, reflecting a higher diversity of the 16S rRNA gene in these libraries than in the [13C]DNA-derived libraries.
Partial sequences were obtained for OTUs in each library containing at least two clones and for 10% of the OTUs containing a single clone, and these were included into a similarity matrix in order to categorize the sequences into phylotypes. Of the 90 OTUs sequenced, a total of 68 phylotypes were identified (Table 1). The phylogenetic analysis of the sequences obtained in [13C]DNA fractions revealed that these phylotypes were associated with diverse bacterial lineages, mainly in the phyla Proteobacteria, Bacteroidetes, Gemmatimonadetes, Acidobacteria, and Actinobacteria (Table 1; Fig. 1A, B, C, D, and E).
Phylotypes that were more represented in the libraries originating from [13C]DNA fractions than in the libraries originating from [12C]DNA fractions were of special interest for this study, as enrichment in [13C]DNA in the respective organisms could be indicative of active metabolism of C1 compounds. Two phylotypes following this pattern were closely related to the sequences of known methylotrophic species belonging to the family Methylophilaceae, and these were enriched in the methanol, methylamine, and formaldehyde microcosms (Table 1). However, other phylotypes following a similar pattern were not related to any known methylotrophs (Table 1). Phylotypes related to Novosphingobium were enriched in the methanol microcosm, phylotypes related to the family Gemmatimonadaceae were enriched in the formaldehyde and formate microcosms, and a phylotype related to Xanthomonanadaceae was enriched in the formate microcosm.
Diversity of fae in [13C]DNA fractions isolated from microcosms.
A single fae clone library was
constructed for each microcosm, using the [13C]DNA fraction
as a template. The libraries were analyzed by RFLP and categorized into
OTUs as described above (Table
2). The 31 clones of the
library originating from the methanol microcosm were assigned to seven
OTUs, of which three contained 27 clones. The 47 clones of the library
originating from the methylamine microcosm were assigned to six OTUs,
of which three contained 44 clones. The 41 clones of the library
originating from the formaldehyde microcosm were assigned to seven
OTUs, of which two contained 30 clones. The 43 clones of the library
originating from the formate microcosm were assigned to five OTUs, of
which two contained 32 clones. Representatives of all the OTUs
identified were sequenced and grouped into a total of 15 phylotypes
(Table 2; Fig.
3). Of the three
phylotypes identified in the methanol microcosm, one was associated
with Methylosinus sequences, the second was closely related to
the Rubrivivax sequences, and the third clustered with the
sequences with no cultivated representatives
(18). The methylamine
microcosm also revealed the presence of fae sequences
belonging to Methylosinus spp. and sequences related to the
Rubrivivax-like sequences, with the rest of the phylotypes
affiliated with the unknown organisms belonging to the Alpha-
and Beta-/Gammaproteobacteria. Four (27 clones) of
the five phylotypes identified in the formaldehyde microcosm were
related to the sequences clustering with those of
Beta-/Gammaproteobacteria (related to both
Rubrivivax and Burkholderia), while the fifth
phylotype (14 clones) clustered with sequences previously identified in
Lake Washington that are not related to the sequences of
Proteobacteria, Planctomycetes, or archaeal homologs
of fae (18). The
phylotypes identified in the formate microcosm were related to those of
unidentified Alphaproteobacteria and
Beta-/Gammaproteobacteria (Rubrivivax
relatives) and to unaffiliated
sequences.
|
|
|---|
Much less is known about bacterial C1 metabolism in this site that involves substrates less reduced than methane. However, recent tests employing primers targeting one of the most ubiquitous C1 oxidation pathways, the H4MPT-linked formaldehyde oxidation pathway, revealed the presence of a diverse population of microbes capable of C1 oxidation (18, 19). Further experiments involving incubations of Lake Washington sediment communities in the presence of various C1 substrates have resulted in the enrichment of a number of novel phylotypes, pointing toward the presence of yet-unidentified methylotrophs in Lake Washington (19). In this work, we continued to characterize the active methylotroph population in the Lake Washington sediment via two complementary culture-independent approaches, analysis of rRNA and mRNA by using the RT-PCR-based technique and SIP with a variety of C1 substrates, with emphasis on non-methane utilizers. RT-PCR analysis of both rRNA and mRNA has been previously employed to characterize active microbial populations in a number of environments, including methanotroph populations (10, 15, 32, 35, 45). SIP has been previously successfully applied to detect active microbial populations involved in aerobic methane and methanol oxidation, aerobic aromatic compound oxidation, or anaerobic propionate oxidation in environments such as soil (36, 38), rice field (27), bacterial mat (16), contaminated aquifer (17) and soda lake (26). Here we extended both the range of environments and the range of the C1 substrates used in SIP experiments.
Analysis of 16S rRNA isolated directly from the Lake Washington sediment revealed a high diversity of microorganisms present. The majority of the sequences recovered were affiliated with the Beta-, Delta-, and Gammaproteobacteria, phylogenetic groups that include organisms of very diverse phenotypes, lifestyles, and trophic capabilities (22). In addition, a number of the sequences were associated with phylogenetic groups lacking or with very few cultured representatives (e.g., Bacteroidetes, Actinobacteria, Acidobacteria, Gemmatimonadetes, and candidate division WS3). The presence of such a variety of microorganisms in the sample may reflect the wide range of environmental conditions and metabolic activities that take place in the top layer of Lake Washington sediment, one of these being methane oxidation. The data obtained in this work support the previous observations on the dominant role of the gammaproteobacterial methanotrophs, based on rRNA and functional gene mRNA (pmoA) detection. Another functional gene employed in this study, fae, was used to identify active populations capable of C1 oxidation downstream of formaldehyde. Most of the fae sequences uncovered were affiliated with those of Methylosinus spp. or belonged to uncultivated organisms of the Beta- and Gammaproteobacteria, indicating that the latter may also contribute to the cycling of C1 compounds in Lake Washington. The lack of fae sequences related to known methanotrophs of the Gammaproteobacteria (Methylomonas spp. and Methylobacter spp.) from PCR-amplified clone libraries has been reported before and has been attributed to the toxic effect of this gene on Escherichia coli (18), underscoring the importance of using multiple culture-independent methods for assessing functional diversity in natural habitats.
In a previously conducted inventory of fae genes in Lake Washington, based on PCR amplification from environmental DNA, many of the sequences fell within a novel fae cluster and diverged significantly from known fae or fae homologs characterized for Proteobacteria, Planctomycetes, or Archaea (18). The phylogenetic identities of the respective microbes, their lifestyles, and their roles in elemental cycling in the site remain unknown. In this study, transcripts of the sequences related to this novel fae cluster were identified (phylotypes pLWFae-18 and pLWFae-19), suggesting that these divergent fae genes are expressed in Lake Washington and suggesting a role for this unknown group in C1 cycling.
While the analysis of fae transcripts has revealed the presence of active organisms other than methanotrophs that are capable of C1 metabolism, the range of C1 substrates potentially used by these organisms remained largely unknown. We used the stable-isotope probing technique in order to link the identity of these organisms to their substrate specificity and thus to their potential function. The analysis of 16S rRNA gene libraries implied that Methylophylaceae may be involved not only in utilization of methanol and methylamine but also in the utilization of formaldehyde in the environment. This approach also allowed us to identify novel phyla potentially involved in the use of C1 compounds. Our data suggest that Alphaproteobacteria of the genus Novosphingobium are candidates for methanol utilization, Gemmatimonadaceae are candidates for utilization of both formaldehyde and formate, and Xanthomonadaceae are candidates for utilization of formate. These results now suggest experiments focused on isolation and taxonomic characterization of the candidate organisms to test these predictions. The fae libraries also provided evidence suggesting that Methylosinus may be involved not only in the oxidation of methane but also in methanol and formaldehyde oxidation. In addition, many of the sequences in the fae libraries belonged to unidentified species, suggesting that novel groups of bacteria are involved in C1 cycling in this habitat. As the SIP technique has certain inherent biases, for example, toward fast-growing organisms, these results should be used only as hints for future experiments uncovering the potential roles of these phyla in C1 metabolism in Lake Washington sediment.
In conclusion, we here describe the use of a combination of two culture-independent molecular techniques (RT-PCR and SIP of DNA) to identify active bacterial populations involved in utilization of C1 compounds in the sediment of Lake Washington, a freshwater lake. The diversity of 16S rRNA in our library reflected the high level of complexity of the active microbial community in the sediment. Using functional genes that are signatures of methylotrophy, pmoA and fae, we identified genes from known methanotrophs and methylotrophs such as Methylobacter, Methylomonas, Methylosinus, Methylobacterium and also novel sequences affiliated with uncultivated organisms of the Alpha-, Beta-, and Gammaproteobacteria, as well as divergent lineages possibly representing novel phyla. Using SIP, we were able to identify a few known methylotroph species as being active in utilization of methanol, methyl- amine, formaldehyde, and formate. However, most of the sequences obtained were not affiliated with any cultivated organisms. The combined results of these two approaches suggest that many bacteria involved in C1 cycling in the site remain uncultivated and uncharacterized. These data highlight the existing gaps in the understanding of C1 cycling, especially downstream of methane, in Lake Washington and likely globally and form the basis for culture-dependent investigations to identify these new groups of bacteria.
We are grateful to the crew of R/V Clifford Barnes and S. Stolyar for help with sample acquisition, to D. Stahl for sharing a UNIX server, and to G. Jacobson for technical assistance. [13C]methanol was provided by the National Stable Isotope Resource at Los Alamos National Laboratory.
Present address: CEA/Cadarache, DSV/DEVM/LEMiR, B
t 161, 13108 St-Paul lez Durance, France. ![]()
Present address: Department of Microbiology, University of Washington, Seattle, WA 98195. ![]()
|
|
|---|
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»