Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2007, p. 594-604, Vol. 73, No. 2
0099-2240/07/$08.00+0 doi:10.1128/AEM.01875-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
Environmental Biotechnology, Department of Biotechnology, Delft University of Technology, NL-2628 BC Delft, The Netherlands
Received 7 August 2006/ Accepted 1 November 2006
|
|
|---|
|
|
|---|
The major metabolic processes that take place in anaerobic bioreactors, such as methanogenesis, sulfidogenesis, and acetogenesis, are nowadays well understood. However, our knowledge of the diversity and dynamics of the microbial communities responsible for these processes is still limited. This is because microbial communities in large-scale biotechnological processes, such as wastewater treatment, are often treated as a black box (18). This is not due to the underestimation of the significance of the biological component but is due to our inability to isolate most of the microorganisms in pure cultures. Fortunately, molecular techniques have provided alternative approaches to overcome the problems associated with culture-dependent analysis of complex microbial communities (2).
Different molecular techniques, such as fluorescence in situ hybridization (FISH) (52), PCR-denaturing gradient gel electrophoresis (DGGE) (9), and DNA microarrays (25), have been used to study the diversity of SRB in natural and engineered ecosystems. Most of these studies have focused on the presence rather than on the activity of SRB in the samples. There have been only a few studies (64, 65) in which the metabolically active populations were monitored by targeting the mRNA of the [NiFe] hydrogenase gene fragments of SRB. However, the limited distribution of this gene among the SRB has restricted its use to the study of Desulfovibrio species only (62).
dsrAB
encodes the
and ß subunits of an enzyme that
catalyzes the six-electron reduction of sulfite to sulfide
(57). Due to a remarkably
high degree of conservation observed in dsrAB across
sulfate-reducing bacteria and archaea
(30), it is a potential
candidate for phylogenetic studies of these organisms. Previous studies
using partial sequences of dsrAB to evaluate the phylogeny of
different SRB lineages has revealed topology congruent with 16S rRNA
gene-based phylogenetic tree
(68). The dsrAB
gene-based molecular approach has been used to discriminate among SRB
in diverse environments
(4,
12,
15,
34,
43). Recently, Geets and
coworkers described DGGE of PCR-amplified dsrB gene fragments
to specifically follow population dynamics of SRB
(16).
This paper describes the results of a comparative study of SRB in lab- and full-scale sulfidogenic wastewater treatment reactors using 16S rRNA and dsrB gene fragments as molecular markers. The goals of our research were to obtain insight into the diversity of sulfate-reducing bacteria in different lab- and full-scale reactors treating sulfate-rich industrial wastewater and to identify the metabolically active community members by comparative DGGE analysis of PCR products obtained from both DNA and RNA of the 16S rRNA and the dsrAB genes. Apart from DGGE analysis, hybridization probes were designed to validate the presence and to determine the abundance of dominant SRB by using FISH.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Characteristics
of the reactors
|
Reverse transcription of RNA and PCR amplification.
Reverse
transcription of isolated RNA into cDNA was carried out using the
iScript cDNA synthesis kit (Bio-Rad, CA) according to the protocol
provided by the manufacturer, using 1 µl (80 to 100 ng) of the
RNA template. Amplification of 16S rRNA and dsrB gene
fragments was performed using the primer pairs 341F-GC (5' CCT
ACG GGA GGC AGC AG 3')/907R (5' CCG TCA ATT CMT TTG AGT
TT 3') (33) and
DSRp2060F-GC (5' CAA CAT CGT YCA YAC CCA GGG 3')
(16)/DSR4R (5'GTG
TAG CAG TTA CCG CA 3')
(63), respectively. We
used 1 µl of genomic DNA and 2 µl of cDNA as templates
for the amplification reactions. The protocol used for the
amplification of 16S rRNA gene fragments or cDNA was same as described
previously (33). However,
the protocol for amplification of the dsrB gene fragment as
described by Geets et al.
(16) was modified
slightly to increase the specificity of the amplification reaction. A
"touchdown" protocol was used, wherein the annealing
temperature was decreased from 65°C to 55°C in 20
cycles. Thermal cycling was carried out as follows: 5 min of initial
denaturation of DNA/cDNA at 95°C, followed by 20 cycles of
denaturation at 95°C for 40 s, a
"touchdown"-annealing step for 40 s, and
elongation at 72°C for 1 min. This was followed by another 30
cycles of denaturation at 95°C for 40 s, annealing at
55°C for 40 s, and elongation at 72°C for 1
min. Amplification was completed by a final elongation step at
72°C for 10 min. DNA from Desulfobulbus propionicus
was used as a positive control and water as a negative control in all
PCR amplifications. The quality of the PCR products was examined on 1%
(wt/vol) agarose gels, and the yield was quantified by absorption
spectrophotometry using the NanoDrop ND-1000 (NanoDrop Technologies,
DE).
DGGE of 16S rRNA and dsrB gene fragments.
DGGE was
performed as described by Schäfer and Muyzer
(56) using the D-Code
system (Bio-Rad Laboratories, CA). Electrophoresis was performed with
1-mm-thick 6% polyacrylamide gels (ratio of acrylamide to
bisacrylamide, 40:1) submerged in 1x TAE buffer (40 mM Tris, 40
mM acetic acid, 1 mM EDTA, pH 7.5) at a constant temperature of
60°C. PCR product in an amount ranging from 300 to 500 ng was
applied to the individual lanes on the gel. The electrophoresis
conditions for the 16S rRNA gene fragment were the same as
described previously
(56): 16 h at
100 V in a linear 20 to 80% denaturant gradient (100% denaturant is a
mixture of 7 M urea and 40% [vol/vol] formamide). However, the
conditions used for dsrB gene fragments were based on the
results of a perpendicular DGGE (see Fig.
2A) and a "time
travel" experiment (see Fig.
2B): 6 h at 150
V in a linear 30% to 65% denaturant gradient. After electrophoresis,
the gels were incubated for 30 min in Milli-Q water containing ethidium
bromide (0.5 µg/ml), rinsed for 20 min in Milli-Q water, and
photographed using a Bio-Rad GelDoc station (Bio-Rad, CA). Individual
bands were excised, resuspended in 20 µl of Milli-Q water, and
stored overnight at 4°C. A volume of 3 to 5 µl of the
supernatant was used for reamplification with the original primer sets.
The reamplified PCR products were run again on a denaturing gradient
gel to check their purity. Prior to sequencing, the PCR products were
purified using the Qiaquick PCR purification kit (QIAGEN GmbH, Hilden,
Germany).
![]() View larger version (47K): [in a new window] |
FIG. 2. (A)
Negative image of a perpendicular denaturing gradient gel of
PCR-amplified dsrB fragments from Desulfobulbus
propionicus obtained with primer pair DSRp2060F-GC and DSR4R. The
black dot indicates a urea-formamide concentration of 46%. The white
dots indicate the urea-formamide concentrations (30% and 65%) used for
the "time travel" experiment. (B)
"Time travel" experiment with PCR-amplified
dsrB fragments from Desulfobulbus propionicus (1) and
Desulfomicrobium escambience (2). A mixture of the two
fragments was loaded onto the gel every 15 min for a total of 240 min.
The electrophoresis conditions were 150 V for 4
h.
|
The partial dsrB gene sequences and the deduced amino acids were first analyzed and aligned using the BioEdit (version 7.0.5) sequence alignment editor (www.mbio.ncsu.edu/bioedit/bioedit.html). Subsequently, the sequences were imported into the ARB software program (26), in which the alignment was further corrected manually. The nucleotide sequences were aligned according to the alignment of deduced amino acid sequences. The alignment regions of insertions and deletions were omitted using a suitable alignment mask (indel filter). A full-length dsrAB consensus tree was constructed after comparing the topologies of phylogenetic trees calculated by maximum parsimony, neighbor-joining, and maximum-likelihood analyses. For tree reconstruction, only nearly full-length sequences were considered. Partial sequences were then inserted into the reconstructed tree by applying parsimony criteria, without allowing changes in the overall tree topology.
Design of oligonucleotide probes.
Specific probes for
the dominant sulfate-reducing bacteria were designed using the Probe
Design function in the ARB software
(26). The complete 16S
rRNA sequences of Desulfobacca acetoxidans and
Desulfovibrio sulfodismutans and the partial sequences
obtained from the excised DGGE bands were used for designing probes.
The probes were named with a number that indicates the position of the
first base in the target sequence (by Escherichia coli
numbering). The following fluorescently labeled oligonucleotide probes
were used in this study: (i) DSV827, a probe for members of the genus
Desulfovibrio (5' GGT CGC CCC CCG ACA CCT 3')
(this study); (ii) DSBA1017, a probe for Desulfobacca
acetoxidans (5' GTT GCC AGG CAC CCC CAT 3') (this
study); and (iii) DSR660, a probe for members of the genus
Desulfobulbus (5' GAA TTC CAC TTT CCC CTC TG
3') (11).
Desulfobacca acetoxidans DSM11109 and Desulfovibrio
sulfodismutans DSM3696 were used as reference strains to check the
specificity of the designed probes. Hybridization stringencies were
determined by performing hybridizations with increasing formamide
concentrations of 10% (vol/vol), 20% (vol/vol), 30% (vol/vol), and 40%
(vol/vol) using the reference
strains.
Whole-cell hybridization.
Samples from
the reactors were washed with 10 mM sodium phosphate buffer (pH 7.2)
containing 130 mM NaCl and were resuspended in the same buffer. The
cells were fixed with 4% (wt/vol) paraformaldehyde in potassium
phosphate buffer for at least 1 h. The cell suspension was
subsequently immobilized on Teflon-coated multiwell microscopic slides.
Hybridization was carried out according to the protocol described by
Pernthaler et al. (45),
using 30% (vol/vol) formamide for probes DSV827 and DSBA1017 and 40%
(vol/vol) formamide for probe DSR660. Following hybridization, the
slides were washed in washing buffer (5 mM EDTA, 20 mM Tris [pH 8.0],
112 mM NaCl, and 0.01% [wt/vol] sodium dodecyl sulfate) at 48°C
for 20 min and then rinsed with Milli-Q water. The slides
were embedded in Vectashield (Vector Laboratories, Burlingame, CA) and
observed with a Zeiss Axioplan 2 epifluorescence microscope. Images
were acquired with Leica FW4000 software. Cell counts of dominant SRB
were determined as described by Neef et al.
(36). The hybridized
cells were analyzed by two independent observers for determining the
fraction of positive signal from each probe relative to the signal
visualized with general probes for bacteria (i.e., probe EUB338
[5' GCT GCC TCC CGT AGG AGT 3']
[1]) and for SRB (i.e.,
probes SRB385 [5' CGG CGT CGC TGC GTC AGG 3'] and
SRB385Db [5' CGG CGT TGC TGC GTC AGG 3']
[1,
50]) or with the general
DNA stain DAPI (4',6'-diamidino-2-phenylindole). The
hybridization experiments were done in triplicate using different
fluorochromes for each probe; different microscopic fields on each
slide were analyzed to confirm the
results.
Nucleotide sequence accession numbers.
The sequences
determined in this study were submitted to GenBank under accession
numbers DQ514554 to
DQ514564 for the dsrB sequences and
DQ514565 to DQ514583 for
the 16S rRNA
sequences.
|
|
|---|
![]() View larger version (45K): [in a new window] |
FIG. 1. (A)
DGGE analysis of 16S rRNA gene fragments using DNA and RNA samples from
different sulfidogenic anaerobic bioreactors as templates. Lane 1, DNA
sample from reactor A; lane 2, RNA sample from reactor A; lane 3, DNA
sample from reactor B; lane 4, RNA sample from reactor B; lane 5, DNA
sample from reactor C; lane 6, RNA sample from reactor C; lane 7, DNA
sample from reactor D; lane 8, RNA sample from reactor D; lane 9, DNA
sample from reactor E; lane 10, RNA sample from reactor E; lane 11, DNA
sample from reactor F; lane 12, RNA sample from reactor F. Bands
indicated with numbers were excised from the gels and sequenced.
(B) Phylogenetic tree based on 16S rRNA gene sequences
obtained from the DGGE bands. Sequences determined in this study are in
boldface; the band number is preceded by 16S. The sequence accession
numbers are shown in parentheses. The sequence of Archaeoglobus
fulgidus was used
as an outgroup but was pruned from the tree. A black dot indicates a
bootstrap value of between 90 and 100%. The scale bar indicates 10%
sequence
difference.
|
Optimization of the dsr DGGE.
A PCR product of
approximately 390 bp was amplified using the primers
DSRp2060F-GC and DSR4R. The optimal gradient of the
denaturant concentration for the amplified product was determined by
performing a perpendicular DGGE. We loaded 600 ng of amplified PCR
product from a pure culture of Desulfobulbus propionicus into
a single large well on a polyacrylamide gel, in which the denaturant
gradient was perpendicular to the direction of electrophoresis
(65). After running
through the perpendicular gel, the dsrB gene fragment appeared
as a sigmoid curve (Fig.
2A). At low concentrations of denaturants (i.e., from 0% to approximately
40% urea-formamide), the fragment runs as a double-stranded molecule
with no denaturation of different melting domains. At a concentration
higher than 55% urea-formamide, the fragment undergoes direct melting
and is held together only by the GC clamp. At a denaturant
concentration of between 40 and 50%, the fragment displays a reduced
mobility with transitional denaturation of the melting domain, with a
melting temperature corresponding to 46% denaturants. From the
perpendicular gradient analysis, a denaturant concentration gradient of
30 to 65% was defined to resolve different dsrB sequence
variants in parallel denaturing gradient gels.
We also performed a "time travel" experiment with dsrB fragments obtained from Desulfobulbus propionicus and Desulfomicrobium escambiense to determine the optimal electrophoresis time (Fig. 2B). PCR fragments were loaded at 15-min intervals for up to 240 min (4 h) onto a polyacrylamide gel containing a 30% to 65% linear gradient of denaturants. After ca. 120 min the individual fragments from the two pure culture samples started to separate, and they were clearly distinguished after 4 h. Even after 4 h of electrophoresis the fragments were still migrating through the gel, although at a considerably reduced mobility. We thus decided on an electrophoresis time of 6 h at a constant voltage of 150 V to obtain a good separation between dsrB gene fragments from different sulfate reducers.
DGGE analysis of dsrB gene fragments.
The presence of dormant and
metabolically active members of SRB in the sulfidogenic bioreactors was
also studied using DGGE of dsrB gene fragments, amplified in
parallel from genomic DNA as well as mRNA. The DGGE analysis of the PCR
products is shown in Fig.
3A. A maximum of two dominant bands were observed in most of the reactor
samples, with an additional one to three less-intense bands appearing
in the DNA-derived profiles of reactors A and C. The DNA- and
mRNA-derived patterns in reactor B were similar and showed the presence
of one dominant band. The DGGE profile of the mRNA-derived PCR product
from a reactor C sample was markedly different from the corresponding
DNA-derived pattern. mRNA amplification products showed very weak bands
corresponding to bands at similar positions in the DNA-derived profile.
On the other hand, one intense band was also observed in the mRNA
profile, which after repeated attempts gave an ambiguous sequence and
hence could not be included in the phylogenetic studies. Reactors D and
E gave similar patterns for both the DNA- and mRNA-based PCR products;
the only difference is in the intensity of band 3. The intensity of
band 3 was higher in the RNA pattern than in the DNA-derived profile,
suggesting a highly active population of SRB represented by band 3. The
DNA- and RNA-derived patterns in reactor F were similar, although the
bands in the RNA-derived profile were less intense, suggesting
less-active SRB populations in reactor F. Bands 8 and 9
were found in the profiles of reactor F only, which suggests a
different SRB population dominating this reactor. The dominant bands
were excised and sequenced, and their phylogenetic affiliations were
analyzed and depicted in a phylogenetic tree (Fig.
3B). Bands with similar
melting behaviors had the same sequences. In general, the sequencing
results were consistent with the results obtained from the 16S rRNA
gene sequences. All lab- and full-scale reactors, except for reactor F,
showed the presence of Desulfobacca acetoxidans-like sequences
(i.e., bands 1, 2, and 3), as observed in the 16S rRNA DGGE profiles.
These bands were found in both the DNA- and the RNA-derived
dsrB profiles. The second most prominent group of bands (i.e.,
bands 4, 5, and 6) that occurred in all the reactors, except reactor B,
were closely related to Desulfobulbus rhabdoformis. One of the
sequences from reactor F (i.e., bands 8 and 9) was closely affiliated
with Desulfoarculus baarsii. Furthermore, sequences of
bacteria closely affiliated with the genera Desulfovibrio,
Desulfomicrobium, and Desulfococcus were found, but
the bands (bands 7, 10, and 11) were low in
intensity.
![]() View larger version (38K): [in a new window] |
FIG. 3. (A)
DGGE patterns of dsrB gene fragments using DNA and RNA samples
from different sulfidogenic anaerobic bioreactors as templates (see the
legend to Fig. 1A for
specifications of the samples). Bands indicated with numbers were
excised from the gels and sequenced. (B) Phylogenetic
consensus tree for dsrAB amino acid sequences deduced from
nearly full-length dsrAB sequences. Branching orders that were
not supported by all treeing methods are shown as multifurcations.
Partial sequences were individually added to the reconstructed
consensus tree by applying parsimony criteria without allowing changes
in the overall tree topology. Sequences determined in this study are in
boldface; the band number is preceded by DSR. The sequence accession
numbers are shown in parentheses. The sequences of
Thermodesulfovibrio islandicus and Thermodesulfovibrio
yellowstonii were used as an outgroup but were pruned from the
tree. A black dot indicates a bootstrap value of between 90 and 100%.
Bootstrap values were calculated only for nearly full-length
dsrAB sequences. The scale bar indicates 10% sequence
difference.
|
|
View this table: [in a new window] |
TABLE 2. Relative
abundances of SRB in different sulfidogenic reactors
|
![]() View larger version (52K): [in a new window] |
FIG. 4. Whole-cell
hybridization of reactor samples A to F with probe DSBA1017 (specific
for Desulfobacca acetoxidans and labeled with Fluos) (green),
probe DSR660 (specific for Desulfobulbus and labeled with Cy3)
(red), and probe DSV827 (specific for Desulfovibrio and
labeled with Cy5) (blue). Bar, 20
µm.
|
|
|
|---|
Microbial community analysis using 16S rRNA as a molecular marker.
DGGE analysis of the PCR-amplified 16S
rRNA gene was used to obtain an estimate of the microbial diversity
within the reactors. To obtain an overview of the metabolically active
members in the microbial communities, reverse transcription-PCR
amplification of 16S rRNA molecules was used. The analysis of rRNA as a
means of inferring the composition of the active members in the total
bacterial community has been proposed for many years
(46,
66). The rRNA content as
a function of growth rate in E. coli
(10) is well established.
Poulsen et al. (49) found
nearly a linear relation between rRNA content and growth rate of a
specific SRB population in anaerobic biofilms. Since then an increasing
number of studies have been conducted to characterize the metabolically
active bacterial and archaeal populations in diverse environments
(20,
28,
29,
31,
38). Comparative analysis
of PCR products obtained by amplification of rRNA genes and
reverse-transcribed rRNA resulted in different DGGE profiles,
indicating a marked difference in the microorganisms present and those
that are metabolically active. This difference in patterns may be
attributed in part to the difference in the ribosome content between
the metabolically active and the dormant populations, assuming that
actively growing cells contain increased levels of rRNA as has been
shown in previous studies
(10,
21,
49). A more complex and
diverse pattern was obtained from 16S rRNA gene-derived PCR amplicons
than from 16S rRNA-derived PCR products, indicating a lesser number of
metabolically active populations. DGGE bands 11, 12, 13, and 14 (Fig.
1A) were of low intensity
in the DNA-derived profiles compared to the RNA-derived profiles.
Assuming that the intensity of DGGE bands is an indirect measure of the
relative abundance (13,
32), bands 11, 12, 13,
and 14 could be attributed to relatively smaller but metabolically
active populations.
Considering the sulfidogenic nature of both the lab- and full-scale reactors, with high sulfate reduction and low electron donor-to-sulfate ratios, most of the 16S rRNA gene/rRNA sequences obtained from the reactor samples, as expected, showed high sequence identity to those of sulfate-reducing bacteria. Sequences were obtained from both completely and incompletely oxidizing sulfate reducers. Sequences of the bands 16S8, 16S12, 16S13, 16S14, and 16S18 grouped with two genera of incompletely oxidizing sulfate reducers, Desulfobulbus and Desulfovibrio (Fig. 1B). The sequence of band 16S8 was closely related to Desulfobulbus rhabdoformis (23), a gram-negative, mesophilic sulfate-reducing bacterium that can oxidize a broad range of substrates, such as propionate, lactate, pyruvate, and ethanol, to acetate and CO2.
The sequences of 16S12, 16S13, and 16S14 closely grouped with members of the genus Desulfovibrio, in particular with Desulfovibrio sulfodismutans, which was isolated from freshwater mud (3). Desulfovibrio sulfodismutans is characterized by the ability to obtain energy for growth by disproportionation of thiosulfate or sulfite to sulfate and sulfide. Growth by dissimilatory sulfate reduction has also been reported (3); however, it is slower than growth by disproportionation. The presence of Desulfovibrio sulfodismutans and Desulfobulbus species in wastewater treatment systems has often been reported in the literature (35, 39, 51, 53, 55, 64). The detection of Desulfovibrio and Desulfobulbus species suggests their fundamental role in the oxidation of carbon through sulfate reduction or through disproportionation of sulfite and thiosulfate in these reactors. Because, Desulfovibrio has a higher growth rate than Desulfobulbus (22), it may be assumed that Desulfovibrio species are mainly responsible for the oxidation of carbon in most of these reactors, although oxidation of some carbon compounds by Desulfobulbus species cannot be ruled out. However, the fermentation product propionate is degraded by Desulfobulbus spp. only, since members of genus Desulfovibrio cannot use propionate as an electron donor. In reactor F, the sequence of 16S18 also clustered with members of the genus Desulfovibrio, with Desulfovibrio gracilis (27) as the closest relative. Desulfovibrio gracilis has been described as moderately halophilic SRB, with optimal NaCl concentrations for growth being 50 to 60 g liter1. The higher salinity in reactor F (Table 1) may explain the presence of a Desulfovibrio gracilis-like SRB in this reactor. Sequence 16S17, obtained from the most dominant band in both the DNA- and RNA-derived DGGE patterns of reactor F, indicates the presence of an active population. Phylogenetic analysis revealed that this organism was closely related to an uncultured sulfate-reducing bacterium. The closest known sulfate-reducing bacterium was found to be Desulfosarcina variabilis. The availability of a variety of complex organic compounds in reactor F, rather than the input of one specific substrate as an energy source, might be responsible for the predominance of nutritionally versatile Desulfosarcina-like SRB (67).
Sequence 16S11, which could be detected in samples from all reactors except reactor F, showed close similarity to Desulfobacca acetoxidans, which oxidizes its carbon source completely in sulfate reduction (Fig. 1B). First isolated from a sulfidogenic bioreactor (41), Desulfobacca acetoxidans has been described as a mesophilic, gram-negative, completely oxidizing sulfate-reducing bacterium that can utilize acetate as the only source of organic carbon and electron donor. The specific growth rate of Desulfobacca acetoxidans (µmax = 0.013 to 0.017 h1) has been reported to be higher than those of acetate-degrading methanogenic archaea, such as Methanosaeta spp. (µmax = 0.003 to 0.012 h1), which dominate methanogenic bioreactors (41). Although archaea were not the focus of this study, it was expected that because of the high sulfate concentrations and limited electron donor quantity in these reactors and the better growth kinetics of Desulfobacca acetoxidans, it would outcompete the acetate-degrading methanogenic archaea.
The sequences 16S4, 16S5, 16S6, and 16S7 were detected in the lab-scale reactors only. The presence of the band was more prominent in reactor B, indicating the presence of a more abundant population. Phylogenetically the sequences grouped with the sequence of Chlorobium limicola (Fig. 1B), an anaerobic, photosynthetic green sulfur bacterium that can use sulfide or sulfur as an electron donor (14). The presence of these bacteria can be explained because the lab-scale reactors were made from glass and were exposed to light. The greenish color of the culture from the lab-scale reactor B at the time of sampling may be attributed to the presence of this bacterium. A likely syntrophic association of this bacterium with sulfate-reducing bacteria, oxidizing sulfide to sulfur and/or sulfate, may be assumed to be present in these reactors. Although present in high numbers in reactor B, the recycling of sulfur does not seem to significantly influence the sulfate reduction, as the percentages of sulfate removal in the two lab reactors were the same (Table 1).
Sequences 16S15 and 16S16 were closely related to the Fe(III)- and Mn(IV)-reducing deltaproteobacterium Geobacter bremensis (58). In the absence of ferric iron, Geobacter species can use alternate electron acceptors, such as nitrate, fumarate, and elemental sulfur, for the oxidation of variety of organic compounds. A coculture of Geobacter spp. in syntrophic association with sulfate-reducing bacteria in the presence of sulfate as an electron acceptor has been described before (7, 54). However, this culture grows very slowly, with a doubling time of more than 7 days.
Microbial community analysis using dsrB as a molecular marker.
We used the dissimilatory sulfite
reductase (dsrAB) gene as a molecular marker to elucidate the
community composition of sulfate-reducing bacteria in the reactors. For
this purpose, the primers originally designed as internal sequencing
primers by Pérez-Jiménez et al.
(44) and later modified
for DGGE analysis by Geets et al.
(16) were used. A 350-bp
dsrB gene fragment was amplified and used for DGGE analysis;
individual bands were sequenced and used for identification of the SRB.
Previous dsr gene-based studies to evaluate SRB communities
used primers DSR1F/DSR4R
(63), amplifying a 1.9-kb
dsr gene fragment. The subsequent phylogenetic studies were
then based either on sequencing the cloned 1.9-kb dsr gene
fragment (12) or on
restriction analysis of the cloned fragment
(24,
43). Although our study
was based on a much shorter sequence of 350 bp, previous studies have
shown that the general topology of the dsrAB trees based on
sequences of different lengths remains consistent
(44), and other studies
(5,
8) have indicated that
shorter sequences of the dsrAB gene can be successfully used
for the phylogenetic analysis of the SRB. We are well aware
that short sequences may not be suitable for detailed phylogenetic
inferences or for tracing the evolutionary history of sulfate-reducing
bacteria, but short sequences are suitable for identification of SRB in
environmental samples.
Detection of dsrB sequences indicated the presence of SRB in these reactors, but this does not provide evidence for metabolically active populations. Therefore, the expression of the dsrAB gene was studied by targeting mRNA, indicating the presence of active populations of sulfate reducers. Validation of this approach was given by Neretin and coworkers (37), who quantified dsr gene expression in Desulfobacterium autotrophicum by using real-time reverse transcription-PCR. They found a strong positive correlation between dsr mRNA concentration and cell-specific sulfate reduction rate. Our DGGE results with dsrB gene fragments indicated a relatively simple sulfate-reducing bacterial community, with similar populations present in all reactors except reactor F. This similarity is probably due to the similar environmental conditions (i.e., pH and temperature) or to the same carbon source being fed to most of the reactors. The complex organic compounds and higher salinity in reactor F (Table 1) might have contributed to the presence of a different SRB community in this reactor.
In most reactors, identification of SRB based on dsrB sequences showed a similarity to the SRB identified by comparative 16S rRNA gene sequence analysis. The prominent dsrB sequences (DSR1/DSR2/DSR3 and DSR4/DSR5/DSR6) were phylogenetically affiliated to Desulfobacca acetoxidans and Desulfobulbus rhabdoformis, respectively. The absence of Desulfobulbus rhabdoformis-like sequences in reactor B was consistent with the results obtained by 16S rRNA gene analysis. Desulfobacca acetoxidans-like sequences were not detected in reactor F, which again was in agreement with 16S rRNA gene analysis. The sequence of bands DSR8/DSR9 grouped with Desulfoarculus baarsii, which is capable of complete oxidation of organic compounds to CO2. The oxidation of acetate by Desulfoarculus baarsii is reported to take place at a low rate, with no substantial formation of biomass (6). This might explain the small amount of sulfate reduction (ca. 20%) in reactor F relative to the amount of sulfate reduction in other reactors.
Although a phylogenetically similar Desulfovibrio-like sequence (i.e., DSR10) was detected in reactor C, the result was in contrast to the 16S rRNA gene/rRNA DGGE, in which Desulfovibrio sulfodismutans-like sequences were found in all the reactors except reactor E. This may be attributed to PCR bias (48, 59), as two different primer sets amplifying two different genes, one targeting the whole bacterial population and the other specific to dissimilatory sulfate reducers, may result in different PCR biases. The dsrB DGGE resulted in the detection of two other SRB-like sequences (DSR7 and DSR11), which were closely related to Desulfomicrobium escambiense and Desulfococcus multivorans, respectively. These results may suggest an increased sensitivity in detection of sulfate-reducing bacteria by dsrB DGGE, which seems possible because the primers are specific for microorganism that are capable of dissimilatory sulfate reduction.
Whole-cell hybridization of SRB.
Since PCR-based
approaches for the analysis of microbial diversity in mixed populations
can be influenced by several constraints
(61), our results based
on PCR-DGGE do not necessarily reflect the abundance of target
sequences in the original sample. We therefore tried to confirm the
relevance of the sequence data by whole-cell hybridization using
fluorescently labeled oligonucleotide probes. Based on the comparative
analysis of the sequences retrieved and those of reference organisms,
three probes, DSBA1017, DSV827, and DSR660, were used in the
hybridization analysis to validate the results obtained by DGGE. The
hybridization results (Fig.
4) not only confirmed the
presence of dominant SRB sequences retrieved by DGGE but also gave an
estimate of the abundance of these bacteria relative to SRB or total
bacterial communities. FISH experiments indicated that
Desulfobacca-like sulfate-reducing bacteria made up a
significant part of the SRB community in the reactors, with the
exception of reactor F, in which no positive hybridization was observed
(Fig. 4 and Table
2). This is consistent
with the DGGE results. Again in accordance with the results obtained by
DGGE, the probe specific for members of the genus
Desulfobulbus showed that Desulfobulbus spp. were
more abundant in reactors D, E, and F than in reactors A and C (Fig.
4 and Table
2). In the latter
reactors, Desulfobulbus spp. do not seem to contribute
significantly to the overall SRB population. Cells detected with DSR660
were observed mainly as aggregates of a few hundred cells. Similar
clusters of cells of Desulfobulbus spp. have been observed
previously in wastewater biofilms
(40).
The detection of 30 to 40% of Desulfovibrio cells in most of the reactor samples with probe DSV827 (Table 2) implies that they are important in sulfate reduction in these reactors. The results are in general in agreement with the 16S rRNA gene DGGE, although dsrB DGGE failed to identify Desulfovibrio-like sequences. This might be due to PCR bias such as the preferential amplification of dsrB gene fragments from Desulfobacca- and Desulfobulbus-like species.
So far, most microbial ecology studies have focused on the diversity of microorganisms. However, more important for the cycling of chemical elements, such as sulfur, are the microorganisms that are active. In this study, we compared the structures and functions of sulfate-reducing bacterial communities in different lab- and full-scale wastewater treatment reactors by targeting both DNA and RNA of two different molecular markers, i.e., 16S rRNA and dsrB. Detection of the gene indicates the presence of microorganisms; however, detection of rRNA or the mRNA of the dsrB ensured that the corresponding SRB was metabolically active at the time of sampling. In general, congruent results were obtained with the two genes. Whole-cell hybridization with oligonucleotide probes targeting the 16S rRNAs of the dominant SRB populations confirmed the results obtained with PCR-DGGE and showed the relative abundances of these populations. The approach is an important step forward to gain insight into the niche differentiation of coexisting sulfate-reducing bacteria in different habitats.
We thank Paques B.V. (Balk, The Netherlands) for providing the samples for this study. We thank Miriam Foti for help in construction of the dsrAB phylogenetic tree and Marzia Miletto and Dimitry Sorokin for helpful discussions.
Published ahead of print on 10 November 2006. ![]()
|
|
|---|
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»