Previous Article | Next Article ![]()
Applied and Environmental Microbiology, April 2004, p. 1982-1989, Vol. 70, No. 4
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.4.1982-1989.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Department of Civil Engineering, University of Minnesota, Minneapolis, Minnesota 55455-0116
Received 31 October 2003/ Accepted 1 January 2004
|
|
|---|
|
|
|---|
MABs are unique compared to biofilms grown on inert substrata. Nutrient concentrations in MABs are highest at the biofilm-liquid boundary and decrease with depth into the biofilm; in contrast, oxygen concentrations are highest at the membrane and decrease toward the outside of the biofilm (Fig. 1). In thick MABs, oxygen partially penetrates the biofilm, generating an anoxic zone at the biofilm-liquid boundary (6). In such biofilms, the region of greatest microbial activity is variable and depends on the relative availability of oxygen and substrate (8, 12). Thick MABs also simultaneously provide favorable conditions for both nitrification (near the membrane) and denitrification (near the biofilm-liquid boundary) within a single biofilm.
![]() View larger version (79K): [in a new window] |
FIG. 1. Hypothetical structure of an MAB treating wastewater.
|
The purpose of this study, therefore, was to characterize the unique stratification of MABs with respect to bacterial community structure and function. We hypothesized that these biofilms would vary substantially in biomass density and respiratory activity as a function of biofilm depth. More importantly, AOB should be found near the membrane, where oxygen concentrations are high and organic concentrations are low. Denitrifying bacteria should preferentially grow in the outer anoxic region of the biofilm, where organic concentrations are high. MABs were grown under well-defined conditions of fluid flow and substrate concentrations. Biofilms were separately grown at different fluid velocities (2 and 14 cm s1) to compare the effect of fluid shear on MAB characteristics. The bacterial community structure was examined by PCR-denaturing gradient gel electrophoresis (PCR-DGGE) and PCR cloning of 16S rRNA genes and by competitive quantitative PCR of genes encoding enzymes required for ammonia oxidation (amoA) and for denitrification (nirS and nirK).
|
|
|---|
For each experiment, the reactor was inoculated with 1 ml of cryopreserved (15% [vol/vol] glycerol) activated sludge collected from the aeration tanks at the St. Paul, Minn., metropolitan wastewater treatment facility. The reactor was fed a mixture containing sodium acetate and ammonium chloride dissolved in dechlorinated tap water. Bulk reactor conditions were monitored daily and maintained at approximately 100 mg of chemical oxygen demand liter1 (mean = 102; standard deviation = 21) and 25 mg of ammonia nitrogen liter1 (mean = 25; standard deviation = 6). Each membrane was supplied with 30 ml of atmospheric air liter1 (pO2 = 0.21 atm) or an oxygen-air mixture (pO2 = 0.75 atm).
Biofilm samples (
1 cm by 1 cm) were removed from the reactor 6 to 8 weeks after inoculation. Samples were excised with a razor blade and immediately frozen at 15°C. Precise dimensions of each biofilm core were measured with a vernier caliper. Biofilms were then thin sliced (typical thickness, 100 to 300 µm) parallel to the membrane in a combined cryostat-microtome (Richard-Allan Scientific, Kalamazoo, Mich.). Biofilm slices were transferred to sterile microcentrifuge tubes and stored separately at 20°C until further analysis.
Analytical methods.
Dissolved oxygen profiles were obtained using a Clark-type microelectrode (tip diameter, 10 to 15 µm; response time, <5 s; spatial resolution, 20 to 30 µm; Unisense, Aarhus, Denmark). The microelectrode was lowered through an access window in the top of the reactor and into the biofilm using a computer-controlled micromanipulator (Oriel Inst., Stratford, Conn.) capable of maneuvering in 0.1-µm increments. Biofilm thickness was determined utilizing a 250x magnification microscope inspection system (Integrated Endoscopy, Irvine, Calif.). The microscope was focused alternately on the biofilm surface and a reference point at the bottom of the reactor. The distance between the membrane surface and the reference point was determined prior to biofilm growth, and the thickness was determined by the distance traveled by the micromanipulator (35).
Respiratory activity was measured by iodonitrotetrazolium chloride (INT) reduction to INT-formazan (39). Biofilm samples (
1 cm by 1 cm) were removed from the reactor and placed into 2 ml of INT solution (9.9 mM in 2% [vol/vol] N,N-dimethylformamide) and 1.5 ml of Tris buffer (1.0 M; pH 7.0) and incubated in the dark for 2 h at 40°C prior to freezing and slicing in the cryostat-microtome. INT-formazan was extracted from the biofilm slices using 1.0 ml of extraction solution (50% [vol/vol] N,N-dimethylformamide in ethanol) and incubated at room temperature in the dark for 1 h. INT-formazan concentrations were quantified using a UV-vis spectrophotometer (
= 464 nm; DU-530; Beckman, Fullerton, Calif.) and calibrated to a formazan standard curve.
Total protein concentrations were quantified by the Lowry method (19) using bovine serum albumin as a standard. Biofilm slices first were dried in a centrifugal drier (Labnet, Edison, N.J.) and resuspended in 1.0 N NaOH prior to analysis. Polysaccharide concentrations were determined by the anthrone method (15) using glucose as a standard. Biofilm slices first were dried in the centrifugal drier and resuspended in deionized water prior to analysis. Dry cell mass concentrations were determined by drying a known volume of biofilm at 103°C overnight and weighing the residue with an analytical balance.
DNA extraction.
Biofilm slices were suspended in 1 ml of lysis buffer (120 mM sodium phosphate buffer [pH 8.0], 5% [wt/vol] sodium dodecyl sulfate). Samples then underwent three successive freeze-thaw cycles followed by an incubation at 70°C for 90 min. Genomic DNA was extracted and purified from these samples using the FastDNA Spin kit per the manufacturer's instructions (QBioGene, Vista, Calif.).
PCR-DGGE.
Partial 16S rRNA genes were amplified from the extracted genomic DNA by PCR using a PTC 100 thermal cycler (MJ Research, Inc., Watertown, Mass.) using primers PRBA338F (18) and PRUN518R (23) with a GC clamp attached to the forward primer (23). The final 50-µl reaction mixture contained 1x PCR buffer (Promega, Madison, Wis.), 175 µmol of MgCl2, 4 nmol of deoxynucleoside triphosphates, 2% bovine serum albumin, 25 pmol of forward and reverse primers, 1.25 U of Taq polymerase (Promega), and
1 ng of template DNA. The PCR protocol included a 5-min initial denaturation at 94°C, 30 cycles of 92°C for 30 s, 55°C for 30 s, and 72°C for 30s, and a final extension for 10 min at 72°C.
DGGE was performed on a D-Code apparatus (Bio-Rad, Hercules, Calif.). Samples containing equal amounts of PCR amplicons were loaded onto 8% (wt/vol) polyacrylamide gels (37.5:1, acrylamide-bisacrylamide) in 0.5x Tris-acetate-EDTA (TAE) buffer (31) using a denaturing gradient ranging from 30 to 50% denaturant (100% denaturant contains 7 M urea, 40% [vol/vol] formamide in 0.5x TAE). Electrophoresis was performed at 60°C, initially at 20 V (15 min) and then at 200 V (200 min). Following electrophoresis, the gel was stained with SYBR Green I (Molecular Probes, Eugene, Oreg.; diluted 1:5,000 in 0.5x TAE). Gels were visualized on a UV transilluminator and photographed with a digital charge-coupled device camera (BioChemi system; UVP, Inc., Upland, Calif.). Photographs were enhanced for contrast and brightness using Adobe Photoshop version 6.0. Specific PCR-DGGE bands were manually excised from the gel, suspended in 30 µl of sterile water, and incubated overnight at room temperature. The PCR-DGGE process was then repeated using these samples as template until only a single band was detectable on a denaturing gel. A final PCR step was performed without the GC clamp attached to the forward primer. PCR products were then purified using the GeneClean II kit (Qbiogene) prior to nucleotide sequence determination.
PCR cloning.
Nearly complete 16S rRNA genes were amplified by PCR as described above except that primers PRBA8F (11) and PRBA1522R (11) were used. PCR products were purified, ligated into the pGEM-T Easy cloning vector (Promega), and transformed into competent Escherichia coli DH5
cells (31). Plasmids were purified by the alkaline lysis procedure (31). Plasmids were then digested with RsaI restriction enzyme, and their fragment lengths were compared by electrophoresis on a 2% agarose gel.
Nucleotide sequence determination.
Nucleotide sequences were determined fully in both directions from purified PCR-DGGE bands and plasmids at the Advanced Genetic Analysis Center at the University of Minnesota using an ABI 3100 genetic analyzer (Applied Biosystems, Foster City, Calif.). PCR-DGGE bands were sequenced using primers PRBA338F and PRUN518R. Nucleotide sequences were determined from plasmids using three forward primers (PRBA8F, PRUN518F, and PRUN907F) and three reverse primers (PRUN518R, PRUN907R, and PRBA1522R) (24). All reported sequences are the consensus of bidirectional sequence information. Reported nucleotide sequences do not include the original PCR primer sequences.
cPCR.
Nitrosomonas oligotropha-like populations were quantified in each biofilm slice by competitive PCR (cPCR) of ammonia monoxygenase (amoA) gene fragments (10). A competitor (10) of 100 nucleotides in length was synthesized by AlphaDNA (Montreal, Canada). PCR mixtures were as described above except that primers AMO598F (10) and AMO718R (10) and AmpliTaq polymerase (Roche Diagnostics, Indianapolis, Ind.) were used. The PCR protocol included a 5-min initial denaturation at 94°C, 40 cycles of 94°C for 1 min, 60°C for 1 min, and 72°C for 1 min, and a final extension for 7 min at 72°C. For each biofilm sample slice (target), cPCR amplifications were performed on four mixtures containing a constant target DNA concentration and decreasing competitor copy numbers. Products were resolved on a 3% (wt/vol) NuSieve 3:1 agarose gel (BioWhittaker, Rockland, Maine) in 1x Tris-borate-EDTA buffer (31) stained with ethidium bromide. Band intensities of target DNA and competitor DNA were quantified using LabWorks Image Acquisition software (UVP).
Denitrifying bacterial populations were quantified also by cPCR of two functionally redundant nitrite reductase gene fragments (nirS and nirK). To synthesize competitor templates, partial nirS and nirK genes were amplified by PCR from bacterial isolates known to contain each gene (nirS, Pseudomonas fluorescens ATCC 33512; nirK, Achromobacter cycloclastes ATCC 21921). Based on the nucleotide sequences of these PCR products (data not shown), internal PCR primers (NIRS-INT, 5'-GAC AGC AGC CAT CGC TA-3'; NIRK-INT, 5'-GCG AGC AGG ACT TCT AC-3') were designed for each gene fragment with the original forward primer attached to the 5' end. PCR with this primer resulted in a 102-nucleotide deletion from the nirS gene fragment (735 nucleotides) and a 59-nucleotide deletion from the nirK gene fragment (474 nucleotides). PCR amplicons of the synthesized competitor DNA were then ligated into the pGEM-T Easy cloning vector and transformed into competent E. coli DH5
cells (31). DNA concentrations in the plasmid extracts were quantified by staining with Hoechst 33258 dye and measured on a TD-700 fluorometer (Turner Designs, Sunnyvale, Calif.) using calf thymus as a DNA standard.
The PCR mixtures for cPCR of nirS and nirK gene fragments were the same as those described above except that different primer pairs were used (nirS, cd3 F and cd4 R; nirK, F1aCu and R3Cu) (13, 22). The PCR protocol for the amplification of nirS gene fragments included a 5-min initial denaturation at 94°C, 40 cycles of 94°C for 1 min, 50°C for 1 min, and 72°C for 1.5 min, and a final extension for 10 min at 72°C. The PCR protocol for the amplification of nirK gene fragments included a 3-min initial denaturation at 94°C, 35 cycles of 94°C for 0.5 min, 57°C for 1 min, and 72°C for 1 min, and a final extension for 10 min at 72°C. cPCR products were resolved on 2% agarose gels (Bio-Rad) in 1x TAE buffer (31) and stained with ethidium bromide.
Data analysis.
Total respiratory activity was computed by determining the area beneath the plot of INT-formazan concentration versus distance from the membrane by trapezoidal rule. An INT-formazan concentration of 1,000 µg ml1 was assumed as a baseline for these calculations.
The concentration of target gene was determined from the cPCR results by computing the y intercept of a linear regression of the log of the competitor copy number versus the log of the ratio of target to corrected competitor band intensity. Competitor band intensity was multiplied by a correction factor to account for the difference in the size of the target and the competitor (amoA, 121/100; nirS, 735/633; nirK, 474/415).
Nucleotide sequences were compared with sequences in the GenBank database (3) using the BLASTn program (2). Nucleotide sequences were checked for possible chimeric sequences using the CHECK_CHIMERA program at the Ribosomal Database Project website (20). Putative chimeric sequences were also manually split into subsections and resubmitted to the GenBank database to determine if the segments were from different phylogenetic groups.
Nucleotide sequence accession number.
Nucleotide sequences have been deposited in the GenBank database under accession numbers AY444967 to AY444996.
|
|
|---|
![]() View larger version (18K): [in a new window] |
FIG. 2. MAB thickness as a function of time for biofilms grown at a fluid velocity of 14 cm s1. , air-fed biofilm; , oxygen-fed biofilm. Error bars represent the standard deviations of the means (n = 9).
|
|
View this table: [in a new window] |
TABLE 1. Summary of pertinent MAB characteristicsa
|
![]() View larger version (17K): [in a new window] |
FIG. 3. Dissolved oxygen profiles of air-fed (a) and oxygen-fed (b) biofilms grown at a fluid velocity of 14 cm s1. Different profiles represent data collected on different days.
|
![]() View larger version (15K): [in a new window] |
FIG. 4. Protein concentration profiles within MABs grown at fluid velocities of 2 cm s1 (a) and 14 cm s1 (b). , air-fed biofilm; , oxygen-fed biofilm. Error bars represent the standard deviations of the means (n = 3).
|
Respiratory activity was also quantified in the biofilms (Fig. 5). For each biofilm, respiratory activity exhibited a hump-shaped pattern as a function of depth. Respiratory activity was low near the membrane and at the biofilm-liquid interface, with the highest activities occurring at an intermediate depth. The location of highest respiratory activity was closer to the membrane for the air-fed biofilms than the oxygen-fed biofilms at both fluid velocities (Table 1). The total respiratory activities were 40 to 80% higher in the biofilms grown at a fluid velocity of 14 cm s1 than in the MABs grown at 2 cm s1 (Table 1).
![]() View larger version (19K): [in a new window] |
FIG. 5. Profiles of respiratory activity as measured by the reduction of INT to INT-formazan within MABs grown at fluid velocities of 2 cm s1 (a) and 14 cm s1 (b). , air-fed biofilm; , oxygen-fed biofilm.
|
![]() View larger version (120K): [in a new window] |
FIG. 6. DGGE of PCR-amplified 16S rRNA gene fragments of an oxygen-fed biofilm grown at a fluid velocity of 2 cm s1. Individual gel lanes are identified by a number indicating the distance (in millimeters) from the membrane where the sample was collected. Letters and arrows identify specific bands that were excised from the gel and sequenced. The results of the nucleotide sequence analysis for these PCR-DGGE bands are shown in Table 2.
|
-Proteobacteria), ß-Proteobacteria, and
-Proteobacteria. |
View this table: [in a new window] |
TABLE 2. Sequence length of the PCR-DGGE bands and clones analyzed in this study and their closest phylogenetic match found in the GenBank databasea
|
-Proteobacteria and ß-Proteobacteria, respectively. In the biofilm slice furthest from the membrane, 29 clones were screened to reveal 11 different DNA sequences. Five of these clones were from the ß-Proteobacteria; other clones were from Bacteriodetes,
-Proteobacteria,
-Proteobacteria, and
-Proteobacteria. Although there were substantial differences between these two bacterial communities, there were two instances (clones M11 and L9 and clones M32 and L11) in which highly similar clones (>99.5% identical) were detected both near the membrane and at the biofilm-liquid interface.
Functional gene-based population profiles.
The bacterial community structure in the MABs was also analyzed by competitive quantitative PCR of genes coding for ammonia oxidation (amoA) or denitrification (nirK and nirS). In the air-fed biofilm grown at a fluid velocity of 2 cm s1, amoA gene fragments were detected only in the biofilm slice closest to the membrane, at a concentration near the detection limit (
400 gene copies, or 2 x 104 gene copies per ml of biofilm). Similarly, amoA genes were not detectable throughout the oxygen-fed biofilm. In both the air-fed and oxygen-fed biofilms grown at a fluid velocity of 2 cm s1, the concentrations of nirS gene fragments exhibited a hump-shaped pattern as a function of depth, with the highest concentration occurring in the middle of the biofilm (data not shown). The nirK gene fragment was difficult to quantify in either biofilm due to the coamplification of a similar-sized fragment.
In contrast, the population profiles in the biofilms grown at a fluid velocity of 14 cm s1 were substantially different (Fig. 7) than those from the MABs grown at 2 cm s1. In both the air-fed and oxygen-fed biofilms, the highest concentration of amoA gene fragments occurred near the membrane in the aerobic region of the biofilm. These concentrations gradually declined with distance from the membrane. In contrast, the concentrations of nirK gene fragments exhibited a hump-shaped pattern as a function of depth, with the highest concentration occurring in the middle of the biofilm. The nirS populations were not detectable near the membrane (air-fed biofilm, <1,000 µm; oxygen-fed biofilm, <1,900 µm) but increased rapidly as a function of distance from the membrane. The profile of nirS populations in the air-fed biofilm was such that the highest concentration of genes occurred near the biofilm-liquid interface. In contrast, the nirS population profile in the oxygen-fed biofilm was hump shaped, with the maximum concentration occurring in the middle of the biofilm. The concentrations of nirS populations were approximately 10-fold higher in the air-fed biofilm than in the oxygen-fed biofilm.
![]() View larger version (31K): [in a new window] |
FIG. 7. Population profiles of ammonia-oxidizing (amoA) and denitrifying (nirS and nirK) bacterial populations as determined by competitive quantitative PCR within air-fed (a) and oxygen-fed (b) biofilms grown at a fluid velocity of 14 cm s1. , nirS; , nirK; , amoA.
|
|
|
|---|
Our results also suggest a niche differentiation between the nirS and nirK denitrifying bacterial populations. In the biofilms grown at 14 cm s1, the nirK populations reached their maximum population density in the middle of the biofilm, whereas the nirS populations reached their maximum in the outer region of the biofilm. The conditions at these locations differed substantially in that the middle of the biofilm was nearer the aerobic zone and should have contained lower concentrations of organic nutrients, whereas the outer region of the biofilm was further from the aerobic zone but contained higher concentrations of organic nutrients. We therefore infer that these conditions are pertinent factors in the competition between nirS and nirK bacterial populations, although other factors that we did not quantify (e.g., nitrite and/or nitrate concentrations) could also be important.
PCR-DGGE and PCR cloning of 16S rRNA genes provided further evidence for the stratification of the bacterial community structure in MABs. These data also supported our cPCR results of amoA gene fragments in that none of the known AOB were detected in the biofilms grown at a fluid velocity of 2 cm s1. Previous researchers demonstrated that the phylogenies of the known AOB give similar results whether amoA or 16S rRNA gene sequences are analyzed (30). The results from the 16S rRNA approach, however, cannot be compared to the cPCR results for nirK and nirS gene fragments, because denitrifying bacteria are of broad phylogenetic distribution (4, 27).
Previous researchers have discussed the possible biases of cultivation-independent community analysis that can occur during nucleic acid extraction (40) and PCR (14, 40). In addition to these biases, our competitive quantitative PCR approach is limited by two other factors. First, gene copy numbers do not necessarily correlate to the number of cells comprising a particular population, because individual cells can have variable copies of a gene within their genome (17, 21). This complicates the comparison of gene copy number of the different population types (amoA, nirS, and nirK) detected herein. Second, all PCR-based assays for community analysis simultaneously must balance generality (i.e., are all of the specific gene targets amplified?) versus specificity (i.e., are all of the nontargeted genes excluded?). The cPCR technique used herein for AOB was selected because previous researchers have demonstrated that it had a high specificity for N. oligotropha-like populations (10). This PCR protocol, however, does not amplify all of the known AOB (e.g., Nitrosomonas europaea ATCC 19718 [data not shown]), and thus it is possible that our cPCR underestimated the number of physiologically relevant AOB in our biofilms. Similarly, the quantification of nirK and nirS genes by cPCR could have overestimated the abundance of these populations if genes other than those coding for nitrite reductase were amplified. The researchers who developed the PCR protocols used herein to amplify the nirK and nirS gene fragments emphasized generality (13) and specificity (22), respectively. Conversely, cPCR could have underestimated the nirK and nirS populations if the PCR amplifications failed to incorporate all of the diversity of nitrite reductase genes that existed in our biofilms. This concern over PCR generality and specificity is particularly pertinent because relatively few studies have performed PCR analysis of nirK and nirS populations.
In addition to the differences in bacterial community structure, substantial differences in bacterial physiology and activity were detected in MABs as functions of depth, flow velocity, and oxygen partial pressure. All of the biofilms in this study were sufficiently thick such that partial oxygen penetration was achieved, resulting in the formation of an anoxic region that ranged in size from 30 to 75% of the total MAB thickness. Although an anoxic region is necessary for anaerobic denitrification to occur, thick MABs can be detrimental to reactor performance (33). MAB thickness and oxygen penetration must therefore be optimized. Oxygen did not penetrate as far at the higher fluid velocity, which was linked to a shift in respiratory activity towards the membrane (i.e., the source of oxygen). This increase in oxygen penetration was possibly due to an increase in biomass density or to an increase in substrate penetration resulting from improved mass transfer efficiency. Our ongoing research will help elucidate the effects of fluid velocity and other reactor conditions on biofilm characteristics, with the goal of optimizing the performance of MAB reactors.
In conclusion, our research demonstrates that the bacterial community structure and activity in MABs is substantially stratified as a function of biofilm depth. Under the proper conditions, this stratification can lead to the simultaneous existence of both nitrifying and denitrifying bacterial populations within a single biofilm. MAB reactors, therefore, have potential utility for the biological treatment of wastewaters containing both carbonaceous and nitrogenous pollutants. Our research also has demonstrated that fluid velocity, which controls the mass transfer of substrate into biofilms, can be an important factor affecting the simultaneous development of both AOB and denitrifying bacterial populations in MABs.
|
|
|---|
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»