Previous Article | Next Article ![]()
Applied and Environmental Microbiology, August 2006, p. 5445-5452, Vol. 72, No. 8
0099-2240/06/$08.00+0 doi:10.1128/AEM.00809-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Department of Microbiology, Swedish University of Agricultural Sciences, Box 7025, SE 750 07 Uppsala, Sweden
Received 6 April 2006/ Accepted 1 June 2006
|
|
|---|
|
|
|---|
Although much is known about process performance after the addition of external carbon sources, little effort has been made to explain the specific effects on the denitrifying communities in activated sludge. The selective pressure from the carbon compound added might affect the composition of the denitrifying community, which would influence metabolic diversity and functional properties of the denitrifiers. For example, it was shown that acetate as the major carbon source in a laboratory-scale reactor selected for denitrifiers with nonflocculating properties, which would adversely affect sedimentation if they became abundant (8). Methanol addition to activated sludge often requires an adaptation period of up to several months before its effect on denitrification rates is fully recognized, but the response to ethanol or acetate seems to be faster (12, 13, 19, 22, 25, 29). Recent studies have shown that ethanol and acetate select for the Azoarcus, Dechloromonas, Thauera, and Acidovorax-like denitrifiers commonly found in treatment processes, which could explain the rapid adaptation to these external carbon sources (8, 17). The lag period observed before the denitrification rate increases after methanol addition has been assumed to be an effect of population shifts within the denitrifying community, and this supposition is strongly supported by the results of Ginige et al. (7). By employing stable isotope probing, they identified members of the family Methylophilales that were enriched in a laboratory-scale reactor operated under denitrifying conditions with methanol as the sole carbon source. Older studies using culture-dependent techniques reported that methanol selects for Hyphomicrobium spp. and Paracoccus spp. in postdenitrification processes (1, 3, 24, 26, 37), i.e., with anoxic zones for denitrification placed before the aerobic ones. Predenitrification processes, with denitrification as the first step, are more complex, and a large variety of carbon compounds is present in the influent wastewater. These systems should therefore allow for more diverse denitrifying communities to thrive than those in the cases mentioned above, where fewer easily degradable carbon sources other than the organic compound added were available.
The objective of this work was to compare the effects of methanol and ethanol as external carbon sources on the genetic and metabolic diversity of denitrifying communities in predenitrifying activated sludge. The genes nirK and nirS, which encode the two structurally different but functionally equal nitrite reductases in the denitrification pathway, were exploited as markers to target the dominant nir genotypes in two parallel pilot-scale plants fed with methanol or ethanol and a full-scale plant without external carbon sources. To compare all samples, we analyzed the nirK communities by denaturing gradient gel electrophoresis (DGGE), whereas clone libraries of nirK and nirS genes were compared in samples from processes completely adapted to methanol or ethanol and from full-scale and pilot-scale processes without external carbon addition. In addition, the community function and metabolic diversity were assessed by non-growth-dependent electron donor use profiling.
|
|
|---|
From the methanol line in the pilot plant, samples were taken 6 days prior to methanol addition (RP2) and after 21 (MP1) and 65 days (MP2) of continuous methanol dosage. To compare the effects of methanol and ethanol addition, the parallel ethanol line was also sampled 6 days prior to ethanol addition (RP1) and on day 28 (EP1) and day 65 (EP2). The full-scale plant without an external carbon source was used as a reference and was sampled twice during the experimental period. The first sample (RF1) was taken on the same day as the RP1 and RP2 samples to compare processes without external carbon addition, and the second (RF2) was taken on day 30 to determine the stability and reproducibility of community patterns and metabolic profiles in the reference plant. On each sampling occasion, a 1-liter sample was withdrawn from the return-sludge flow from the last anaerobic zone in each line and transferred to three 250-ml flasks under anaerobic conditions, and the activity was measured within 1 h of sampling. In addition, 10 1-ml aliquots of each sludge sample were centrifuged for 10 min at 14,000 rpm, and the pellets were frozen at 20°C prior to DNA extraction.
Potential denitrification activity assay for metabolic profiles.
The potential denitrification activities with 10 different electron donors (formate, acetate, propionate, butyrate, ethanol, methanol, glycerol, glucose, glutamine, and benzoate) were measured in triplicate using the acetylene inhibition technique (42) to generate electron donor use profiles for samples from the processes without external carbon addition (RP1, RP2, RF1, and RF2), with methanol addition (MP2), and with ethanol addition (EP2). The electron donors were selected to cover a diversity of metabolic pathways. Rates obtained with only nitrate added (NA) in the assay were the controls. Activated sludge (10 ml) was added to 39 ml 10 mM sodium phosphate buffer (pH 7.2) in glass flasks (118 ml) with an anaerobic atmosphere but was otherwise treated as described by Hallin and Pell (11). The substrate was added by injecting 1 ml of a stock solution to reach an initial concentration of 2 mM KNO3-N and 100 mg C liter1. Samples were incubated at 15°C for 1 h, and 0.5-ml gas samples were withdrawn every 12 min.
Nitrous oxide was analyzed on a gas chromatograph (Chrompack 9000; Chrompack, Rotterdam, The Netherlands) equipped with a 63Ni electron capture detector. Denitrification rates were calculated from linear regression of the N2O produced during incubation and expressed on a mixed-liquor volatile suspended solids (MLVSS) basis (11). MLVSS were measured according to Swedish standards (SS-EN 872, Swedish Standards Institute [http://www.sis.se]). In short, the amount of suspended solids was first determined by vacuum filtration of the sludge sample through a preignited 0.7-µm fiberglass filter and drying at 105°C. The amount of MLVSS was then calculated from the weight loss after incineration at 550°C.
Comparative study.
The data from the electron donor use profiles were merged with those from previously published profiles obtained from other experiments at the same treatment plant (12). Variables from the data in the previous study that had no equivalent in the new data were excluded and vice versa. The data were normalized by dividing each individual variable by the total sum of rates for each sample profile. The total data set was then subjected to principal component analysis in order to elucidate patterns of electron donor use profiles in the differently operated processes and over time, using the Unscrambler software package (Camos AS, Oslo, Norway).
DNA extraction and PCR amplification of nirK and nirS.
Genomic DNAs were extracted in triplicate from all samples, except for the RP1 sample, which was lost, by using a FastDNA Spin kit for soil (Qbiogene Inc., Carlsbad, CA) according to the manufacturer's protocol, with a few modifications. The sludge pellets were first dissolved in 978 µl phosphate buffer supplied with the kit and then transferred to multimix tubes to which 122 µl MT buffer from the kit was added. The isolated DNAs were stored at 20°C.
Amplification of a 473-bp nirK fragment for DGGE analysis and cloning and of a 425-bp nirS fragment for cloning was performed in triplicate reactions for each DNA extract, using the F1aCu/R3Cu primer set, with and without a GC clamp, and the cd3aF/R3cd primer set (36), respectively. The oligonucleotides were purchased from Invitrogen (Carlsbad, CA). The gene fragments were amplified with initial denaturation of the DNA at 94°C for 2 min, followed by 28 cycles of 30 s at 94°C, 1 min at 57°C, and 1 min at 72°C. The reactions were completed after 10 min at 72°C. PCRs were carried out in a minicycler (MJ Research, Waltham, MA), and each 25-µl reaction mix contained 2.5 µl of 10x PCR buffer (500 mM KCl, 15 mM MgCl2, and 100 mM Tris-HCl, pH 9.0, at room temperature), a 200 µM concentration of each deoxynucleoside triphosphate, 1.25 U of Taq polymerase (GE Healthcare, Giles, United Kingdom), 1.0 µM of each primer, and approximately 20 ng DNA. For nirK and nirS reactions, bovine serum albumin (GE Healthcare) was added to reach final concentrations of 400 and 1,000 ng µl1, respectively.
DGGE of nirK.
The nirK amplicons with a GC clamp were resolved by DGGE analysis. Prior to DGGE, the three PCR products from the same DNA extract were pooled, but the triplicate DNA extracts from each sample were treated separately throughout the subsequent analysis. The nirK PCR products were separated using a Dcode system (Bio-Rad Laboratories Inc., Hercules, CA) according to the method of Throbäck et al. (36). Gels composed of 37.5:1 acrylamide:bisacrylamide (7% [vol/vol]) with a gradient of 50 to 70% denaturant (formamide and urea) were cast with a gradient maker (Bio-Rad Laboratories Inc.). After electrophoresis at 130 V for 13 h at 60°C in TAE buffer (40 mM Tris-acetate, 1 mM EDTA, pH 7.5), the migration patterns were visualized by staining with 1:10,000 SYBR green (Molecular Probes, Eugene, Oreg.) for 30 min followed by UV transillumination. Images were documented with a video camera (Ikegami Tsushinki Co., Tokyo, Japan), and the pictures were scanned to create digital images.
Cloning of nirK and nirS and RFLP analysis.
Four clone libraries for each of the amplified nirK and nirS genes were constructed from samples taken from the pilot-scale plant before methanol addition (RP2) and after 65 days of continuous methanol dosage (MP2), from the parallel plant fed with ethanol on day 65 (EP2), and from the full-scale plant at the start of the experimental period (RF1). Prior to cloning, the nirK or nirS PCR products from the same DNA extract were pooled, separated by agarose gel (1%) electrophoresis, and purified with a JETquick gel extraction spin kit (Genomed, Löhne, Germany) to avoid insertions of incorrect amplicons. The PCR products were then cloned by using a TOPO TA cloning kit for sequencing (Invitrogen, Carlsbad, CA). Small amounts of cell material from randomly picked clones were collected with toothpicks and resuspended in 25-µl preprepared PCR mixtures, and the inserts were amplified as described above. In total, 70 nirK and 60 nirS PCR products of the correct size from each library were digested overnight at 37°C, using HaeIII for nirK and HhaI for nirS, and then screened by restriction fragment length polymorphism (RFLP). The digested products were separated by 2% DNA agarose gel electrophoresis at 80 V for 2.0 h.
The RFLP patterns were compared by using Quantity One 1-D analysis software (Bio-Rad Laboratories Inc.) and were grouped into pattern types. Library coverage (C) was estimated as follows: C = 1 nN1, where n is the number of different RFLP pattern types from a clone library that were encountered only once and N is the total number of clones analyzed. The diversity of RFLP pattern types for nirK and nirS was analyzed using Analytic Rarefaction software (http://www.uga.edu/
strata/software/Software.html).
Sequencing and phylogenetic analysis.
A selection of nirK and nirS clones from all RFLP pattern types represented by at least two clones were grown at 37°C overnight in LB medium, and the plasmids were isolated for sequencing using the QIAprep Spin miniprep kit protocol (QIAGEN, Valencia, CA). The nirK and nirS inserts were sequenced on one strand by Macrogen Inc. (Seoul, Korea), using an ABI3730 XL automatic DNA sequencer and the vector primer M13F. For sequencing of nirK bands from DGGE gels, the middle sections of the visible bands were excised from the denaturing gels, placed in 160 µl ultrapure water, and stored at 70°C until used for sequencing. To elute the DNA from the polyacrylamide gels, the samples were thawed for 1 h at room temperature, frozen at 70°C for 1 h, and then thawed again at 8°C for 12 h. The eluted fragments were PCR amplified with the nirK primers F1aCu and R3Cu without a GC clamp, using 4 µl template DNA in a total reaction volume of 50 µl, but otherwise following the method described above. Prior to sequencing, 45 µl of each PCR product was purified with a MicroSpin S-400 HR column (GE Healthcare). Both strands were sequenced in 10-µl reaction mixtures with a DYEnamic ET Terminator cycle sequencing kit (GE Healthcare), using an ABI PRISM 377 (Perkin-Elmer, Wellesley, MA) automated DNA sequencer with F1aCu and R3Cu as sequencing primers. The DGGE bands harboring more than one nirK sequence were further cloned by the procedure described above. From each mixed DGGE band, 10 clones with the correct insert were chosen for sequencing. Both strands of the insert were sequenced as described above, using the vector primers M13F and M13R.
The nirK and nirS sequences obtained in this study were aligned with nirK and nirS gene sequences from cultivated strains and a selection of sequences from environmental nir clones obtained in other studies available from the NCBI database, using CLUSTAL W software (http://www.ebi.ac.uk/clustalw/). Tree analysis was performed with TREECON software (39). Distance matrix analyses were performed with the Jukes and Cantor correction (20). The trees were reconstructed by using the neighbor-joining method of Saitou and Nei (31), and tree topology was evaluated by bootstrap analysis using 100 replicates. The resulting phylograms are given in the supplemental material.
Nucleotide sequence accession numbers.
The partial nirK and nirS gene sequences from this study have been deposited in the GenBank database under accession numbers DQ182155 to DQ182227 and DQ182111 to DQ182154, respectively.
|
|
|---|
![]() View larger version (18K): [in a new window] |
FIG. 1. Potential denitrification rates with 10 organic compounds supplied in excess (means ± standard deviations; n = 3) to generate metabolic profiles for samples from (a) the pilot-scale plants prior to ethanol (RP1) and methanol (RP2) addition and the full-scale plant without an external carbon source on two occasions (RF1 and RF2) and (b) the pilot-scale plants with ethanol (EP2) and methanol (MP2) addition on day 65 in comparison to the pilot-scale plants prior to the addition of ethanol (RP1) and methanol (RP2).
|
![]() View larger version (13K): [in a new window] |
FIG. 2. Score plot (a) and loading plot (b) obtained by principal component analysis of potential denitrification rates with different electron donors to compare the metabolic profiles reported by Hallin and Pell (12) to those from this study, including those for the pilot-scale plants prior to ethanol (RP1) and methanol (RP2) addition, the full-scale plant without an external carbon source on two occasions (RF1 and RF2), and the pilot-scale plants with methanol (MP2) and ethanol (EP2) addition on day 65. Samples in score plot: P, pilot-scale plants; F, full-scale plants; E, addition of ethanol; M, addition of methanol; R, no external carbon source. Variables in loading plot: ME, methanol; ET, ethanol; AC, acetate; PR, propionate; BU, butyrate; GU, glucose; GY, glycerol; BE, benzoate; NA, no additional electron donor. Ellipses indicate processes with similar modes of operation.
|
![]() View larger version (163K): [in a new window] |
FIG. 3. DGGE analysis with SYBR green staining of PCR-amplified partial nirK genes from the pilot-scale plant prior to methanol addition (RP2), the pilot-scale plant with methanol addition on day 21 (MP1) and day 65 (MP2), the pilot-scale plant with ethanol addition on day 28 (EP1) and day 65 (EP2), and the full-scale plant without external carbon addition at the start of the experiment (RF1) and on day 30 (RF2).
|
nirK and nirS clone libraries.
A total of 280 nirK- and 240 nirS-positive clones were screened with RFLP. The coverage for the eight clone libraries, based on RFLP pattern types, varied from 69 to 91% (Table 1). The rarefaction curves for each library did not reach a plateau, indicating that both nirK and nirS genotype diversity was high in all samples (Fig. 4). Nevertheless, the diversity of nirS genotypes was higher than the nirK diversity within each plant operation mode. The lowest level of nir diversity was found in the clone libraries from the full-scale plant without an external carbon source. The rarefaction curves indicated less nirK diversity in the processes with methanol or ethanol addition than in the pilot-scale plant without supplemental carbon (RP2), but there was more nirS diversity in the methanol sample. Of the 280 nirK clones screened, 79 different RFLP pattern types were found, and 41 of these patterns were represented by more than one clone. For nirS, 96 different RFLP patterns were detected, and 36 of these were represented by at least two clones among the 240 screened.
|
View this table: [in a new window] |
TABLE 1. Coverage of nirK and nirS clone libraries based on RFLP pattern typesa
|
![]() View larger version (16K): [in a new window] |
FIG. 4. Rarefaction curves showing the numbers of RFLP pattern types obtained from the numbers of nirK (white) and nirS (black) clones screened from the pilot-scale plant sludge prior to methanol addition (RP2; squares), after methanol adaptation on day 65 (MP2; circles), and after ethanol adaptation on day 65 (EP2; triangles) and the sludge from the full-scale plant without external carbon addition at the start of the experiment (RF1; diamonds). The solid line shows the maximum possible level of diversity.
|
![]() View larger version (46K): [in a new window] |
FIG. 5. Distributions of RFLP pattern types represented by at least two nirK clones in (a) the pilot-scale plant prior to methanol addition (RP2), (b) the methanol-fed pilot-scale plant on day 65 (MP2), (c) the ethanol-fed plant on day 65 (EP2), and (d) the full-scale plant without an external carbon source at the start of the experiment (RF1). White pie slices represent RFLP pattern types specific for the sample, and the numbers refer to the RFLP pattern types.
|
![]() View larger version (35K): [in a new window] |
FIG. 6. Distributions of RFLP pattern types represented by at least two nirS clones in (a) the pilot-scale plant prior to methanol addition (RP2), (b) the methanol-fed pilot-scale plant on day 65 (MP2), (c) the ethanol-fed plant on day 65 (EP2), and (d) the full-scale plant without an external carbon source at the start of the experiment (RF1). White pie slices represent RFLP pattern types specific for the sample, and the numbers refer to the RFLP pattern types.
|
|
|
|---|
The comparative analysis showed that treatment plants within a certain operational mode did not differ substantially from each other. The metabolic pattern stability observed was most likely due to stable operation of the treatment plant and similarity of the influent wastewater over time at the treatment facility. The pattern from principal component analysis was related to that obtained in the study by Hallin and Pell (12), even though the analysis in the present study included not only electron donor use profiles with fewer variables but also profiles from both full-scale and pilot-scale plant experiments conducted in different years. This strengthens the proposed preference for certain electron donors in activated sludge from plants operated either without an external carbon source or with methanol or ethanol addition, as discussed above. The profiles obtained from the pilot-scale plants treated with ethanol in the two experiments run in different years were more similar to each other than were those from full-scale and pilot-scale plants with methanol addition. This is likely explained by more controlled operation of the pilot-scale plants giving more reproducible results but also shows that there are differences that might depend on the size of the treatment plant. Nevertheless, the similarity in the metabolic profiles among the methanol-amended plants indicates that the pilot-scale plant was sufficiently large to mimic a full-scale plant.
Not only was activity affected differently by the external carbon sources, but the diversity in terms of richness in the denitrifying communities was also affected. In comparison to that in the sample obtained prior to commencement of methanol addition, the diversity of nirK genotypes decreased, while the nirS diversity increased when methanol was added to the process. The amounts of nirK diversity were similar in the methanol and ethanol plants, but the amount of nirS diversity was higher in the methanol plant than in the plant with ethanol addition. This indicates that methanol caused the growth of additional denitrifier populations that mainly use methanol even in predenitrification processes, as hypothesized from activity measurements in a previous study by Hallin et al. (10). The level of diversity of nirS genotypes was always higher than the level of nirK diversity within each operational mode. This is in agreement with the nirK and nirS distributions in cultured denitrifiers from a municipal wastewater treatment plant (43) and in clone libraries from a saline metallurgic wastewater treatment plant (41). Results from a biofilm reactor treating wastewater suggest that the concentration of organic compounds is a major factor in the competition between nirK and nirS genotypes, with lower concentrations favoring nirK (4). Yan et al. (40) speculated that nitrate levels might alter the proportions of nirK and nirS denitrifiers in groundwater and that the community dynamics of nirK genotypes affect the nirS community and vice versa.
The composition of nirK and nirS genotypes in the different samples implies that the denitrifier community structure was affected by external carbon sources. The effect was most pronounced on the nirS denitrifiers, while the nirK communities appeared to be more stable. All four nirS libraries differed from each other, showing distinctly different RFLP pattern types. The sample from the ethanol-fed plant differed the most, and the clone library was dominated by treatment-specific nirS RFLP types. The similarity among the nirK libraries fits with the DGGE analysis of nirK genes, which demonstrated that the same dominant denitrifiers of the nirK type were present on all sampling occasions. The resolution is lower for DGGE analysis than for clone library analysis, and therefore minor differences among the treatments could be detected by DGGE. The addition of methanol to the plant resulted in the appearance of an additional band, suggesting that a new nirK population was on the rise. This band was also detected in samples from the ethanol-fed plant. The sequence from this band clustered with other nirK sequences from this study, but nirK sequences from cultivated strains were not found within this cluster. The sequence was identical to sequences of RFLP type 10 in the nirK libraries, but this RFLP type did not dominate any of the clone libraries. Surprisingly, it was found in all libraries except that from the methanol plants, although the coverage of the nirK clone library from this sample was 84%.
As reported in other environmental studies of the functional genes in the denitrification pathway (e.g., see references 5 and 33), most of the dominant nir types in our study clustered with other environmental clones. This could indicate a large diverse population of unknown denitrifying bacteria in activated sludge. Since methanol added to the process increased the capacity of the denitrifiers in the sludge to utilize both methanol and formate for denitrification, methylotrophic denitrifiers would certainly be present. We speculate that some of these unknown nir genotypes are affiliated with nir genes from as yet uncultivated denitrifying methylotrophs. Recent results have shown that many bacteria involved in C1 cycling remain uncultivated and uncharacterized (27), and methanol was shown to select for methylotrophs not previously known as denitrifiers in activated sludge (7). Nevertheless, some of the dominant nirS types in our study were related to nirS genes from Paracoccus, Thauera, Azoarcus, and Alcaligenes faecalis, although none of them could be assigned to any specific operational mode. Accordingly, both cultivation and 16S rRNA-based studies have suggested that the families Comamonadaceae and Rhodocyclaceae, dominated by Azoarcus, Paracoccus, Hydrogenophaga, and Acidovorax species, are key players in denitrification in municipal and industrial activated sludge processes (6, 16, 21, 32). However, the identification of denitrifiers from the retrieved 16S rRNA sequences is difficult. Surprisingly, only a few functional gene approaches to target denitrifiers in municipal activated sludge have been published. Sequencing of the dominant nirK, nirS, and nosZ genes in activated sludge samples from two different full-scale treatment plants revealed that some genotypes were distantly related to those from Rhizobium spp., Azospirillum spp., Rhodobacter sphaeroides, and Paracoccus denitrificans, although the majority clustered with other environmental clones (36). In contrast, most nosZ clones clustered with nosZ sequences from well-known denitrifiers within the
(Azospirillum)- and
-proteobacteria (Pseudomonas) in a biofilm pilot-scale reactor treating urban wastewater (9). Nearly 70% of the nirS clones analyzed in a sequencing batch reactor process with combined nitrogen and phosphorus removal were similar to nirS genes in Azoarcus and Thauera, within the Rhodocyclus group (38).
To conclude, we used a community-level approach to study the effects of methanol or ethanol as an external carbon source on both the structure and function of denitrifying bacterial communities in municipal activated sludge. Our results support an earlier hypothesis that the denitrifier community composition changes not only in postdenitrification processes, but also in predenitrification systems, when external carbon sources are added to speed up denitrification rates. Methanol had a greater impact than ethanol on the metabolic profiles of the communities, and this was linked to an increased level of diversity of denitrifiers of the nirS type. This indicates that denitrifying bacteria favored by methanol were on the rise in the community. The composition of nirS genotypes in the communities varied greatly among the samples, whereas the nirK communities were more stable, suggesting that there might be a niche differentiation between nirS and nirK genotypes in activated sludge processes. The results do not indicate any severe negative effect of either methanol or ethanol on denitrification and the denitrifying communities. Based on this study, we cannot conclude which carbon source to prefer. In the operation of sewage treatment plants, other aspects, such as sludge production, sludge settling properties, economy, etc., are equally important to consider when making a choice.
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
|
|
|---|
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»