Previous Article | Next Article ![]()
Applied and Environmental Microbiology, September 2003, p. 5060-5069, Vol. 69, No. 9
0099-2240/03/$08.00+0 DOI: 10.1128/AEM.69.9.5060-5069.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
George E. Brown Jr. Salinity Laboratory, U.S. Department of AgricultureAgricultural Research Service, Riverside, California 92507,1 Orange County Water District, Fountain Valley, California 927282
Received 28 February 2003/ Accepted 1 July 2003
|
|
|---|
subunit of the ammonia monooxygenase gene (amoA) recovered from soil samples and DGGE bands. DGGE analysis of wetlands and manure samples revealed that the total bacterial community composition was dominated by bacteria from phylogenetic clusters related to Bacillus, Clostridium, Mycoplasma, Eubacterium, and Proteobacteria originally retrieved from the gastrointestinal tracts of mammals. The population of ammonia-oxidizing bacteria showed a higher percentage of Nitrosospira-like sequences from the wetland samples, while a higher percentage of Nitrosomonas-like sequences from manure, feces, raw washwater, and facultative pond was found. These results show that the wetland system is a natural process dependent upon the development of healthy microbial communities for optimal wastewater treatment. |
|
|---|
It has long been recognized that certain microbial groups in animal waste are responsible for breaking down various organic compounds and for the suppression of pathogens in waste. However, very little work has been done on the functional composition of the different bacterial groups in a subsurface flow wetland. Gersberg et al. (2) explored many of the biological parameters that affected the performance of horizontal subsurface flow wetlands. The authors discussed the important roles of aquatic plants, such as bulrushes, in translocating oxygen from the upper parts of the plant to the ends of the roots. They concluded that nitrifying bacterial communities may develop around this microaerobic zone and convert ammonium to nitrate. As the water flows into the anaerobic zone, the denitrifying bacteria convert nitrate to nitrogen gas. Therefore, an aerobic and anaerobic regime develops that contains complementary microbial communities. These communities can nitrify ammonium-laden water under aerobic conditions or convert the nitrate to nitrogen gas by denitrifying bacteria under anaerobic conditions. One advantage of this system over a surface flow system is that the water is being exposed to a far greater surface area of bacterial biofilms, resulting in a higher level of nutrient removal over a shorter time period.
The diversity of microorganisms in the wetland environment may be critical for the proper functioning and maintenance of the system. One such group of bacteria involved in different biological and chemical transformations of organic compounds in wastewater are the chemolithotrophic ammonia-oxidizing bacteria. These bacteria are responsible for the first, rate-limiting step in nitrification in which ammonium (NH4) is transformed to nitrite (NO2) and to nitrate (NO3) during natural nitrogen cycling (7, 15). Oved et al. (13) studied the effects of urban sewage effluent on the community composition and function of ammonia-oxidizing bacteria in soil by using denaturing gradient gel electrophoresis (DGGE) of PCR-amplified fragments of the amoA gene. They found a significant and consistent shift in the population composition of ammonia-oxidizing bacteria in soil irrigated with effluent dominated by Nitrosomonas-like populations, while noneffluent soil was dominated by Nitrosospira-like populations.
The use of the functional gene carrying the
subunit of amoA has been shown to be a good molecular marker for ammonia-oxidizing bacteria in soil (18, 22). The functional gene has been shown to detect only Nitrosospira in agricultural soils. More diverse species within this genus were determined by using the amoA gene rather than the 16S rRNA gene (4). Other investigators have also demonstrated the successful use of a specific set of PCR primers to amplify a fragment of amoA from a variety of pure cultures of ammonia-oxidizing bacteria from environmental samples (3, 9, 17, 18, 22).
This study examined the composition of the general bacterial population and ammonia-oxidizing bacteria in manure, feces, and dairy washwater prior to and after treatment in subsurface wetlands. The first objective was to characterize microbial communities in constructed wetland wastewater and to determine how the composition of the community may influence the final wastewater effluent quality. The second objective was to determine the effectiveness of the wetland treatment technology in reducing contaminants in dairy waste effluent. Present dairy washwater management practices involve storing washwater in ponds or spraying raw washwater onto crops and/or disposal land. Storage in lagoons affects the groundwater quality through the percolation of washwater, contributing significant amounts of nutrients to the groundwater. During a significant storm event, overflow from the storage lagoons also may occur, posing a threat to surface water quality from the high levels of organic material and pathogens present in washwater.
|
|
|---|
Wetland 1 was an end-loading design. The washwater entered through multiple inlets onto the coarse gravel, reeds, and bulrush and then drained into the collection box at the end of the basin. Wetland 2 was a side-loading basin. Washwater entered through multiple inlets along both sides of the wetland and passed through a narrow gravel bed containing reeds and bulrush before being collected through a perforated pipe along the center of the wetlands that drained into the collection box. The combination of open gravel and specific plants was used for the initial removal of suspended solids and subsequent nitrification of the wastewater in the microaerobic zone surrounding the roots of the bulrushes. The overall design of the wetlands was based on studies carried out by Gersberg et al. (2), with some modifications for higher biochemical oxygen demand (BOD), nutrient levels, and suspended solids.
Monthly samples were collected from six locations: wastewater lagoon, facultative pond, and wetlands 1 and 2 influents and effluents. Fresh cow and calf fecal samples and manure samples were collected for microbial analysis. Manure samples were taken from piles of fecal materials that had been deposited for about 2 weeks. Weekly water samples were collected for chemical analysis and determination of the quality of the final effluent water. All samples were collected between December 2000 and September 2001. Weekly samples were analyzed for total Kjeldahl nitrogen (TKN), ammonium nitrogen (NH4-N), nitrate-nitrogen (NO3-N), BOD, total suspended solids (TSS), orthophosphate (PO4-P), pH, potassium (K+), sodium (Na+), chloride (Cl-), total dissolved solids (TDS), chemical oxygen demand (COD), and fecal and total coliform bacteria. All analyses were done with protocols from the Standard Methods of the American Public Health Association (1).
DNA extraction and purification from soil and effluent.
Total bacterial DNA was extracted from 500 mg of fecal or manure samples and from 250 mg of pellets of concentrated effluent sample centrifuged at 3,000 x g for 10 min. DNA was extracted by using UltraClean fecal and water DNA kits (MO BIO, Inc., Solana Beach, Calif.) according to the manufacturer's protocol with slight modifications.
PCR primers and DGGE analysis of total bacterial community.
PCR was performed with about 50 ng of template DNA with the primers PRB 338f and PRUN 518r, located at the V3 region of the 16S rRNA genes of bacterioplankton (14) to assess bacterial community diversity. PCR mixtures for the bacterial 16S rRNA PCR DGGE sequence amplification contained 10 pmol of each primer, Ready-To-Go PCR beads from Amersham-Pharmacia Biotech (Piscataway, N.J.), and sterile distilled water in a final volume of 25 µl; PCR conditions used were those described by Ibekwe et al. (5).
DGGE was performed with 8% (wt/vol) acrylamide gels containing a linear chemical gradient ranging from 30 to 70% denaturant with 100% defined as 7 M urea and 40% formamide. Gels were run for 3 h at 200 V with the Dcode Universal Mutation System (Bio-Rad Laboratories, Hercules, Calif.). DNA was visualized after ethidium bromide staining by UV transillumination and was photographed with a Polaroid camera. Major bands were excised for identification of bacterial species. Bands were placed into sterilized vials with 20 µl of sterilized, distilled water and were stored overnight at 4°C to allow the DNA to passively diffuse out of the gel strips. Ten microliters of eluted DNA was used as the DNA template with the eubacteria primers. DNA was cloned by using the pGEM-T Easy vector system (Promega, Madison, Wis.) and was transformed into Escherichia coli JM109. Isolation of plasmids from E. coli was performed by using the Qiagen plasmid mini kit (Valencia, Calif.). The purified plasmids were sequenced with the ABI PRISM dye terminator cycle sequencing kit with AmpliTaq DNA polymerase, FS (Applied Biosystems, Foster City, Calif.).
DGGE and sequence analyses of ammonia oxidizers.
The diversity of the ammonia-oxidizing bacteria in the samples was performed with amoA primers targeting a partial stretch of the genes that encode the active-site polypeptide of ammonia monooxygenase (18). Products with the expected size of 491 bp were excised from the gel and were purified with a QiaexII gel extraction kit (Qiagen). Purified fragments were cloned as described above. A second amplification was done by a seminested PCR using the amoA primers for DGGE analysis. For the second assay, a GC clamp (5' CCGCCGCGCGGCGGGCGGGGCGGGGGCACGGGG-3') was attached to the forward primer to increase the DGGE gel separation (11). PCRs contained Ready-To-Go PCR beads, 50 pmol of each primer, template DNA, and sterile distilled water to a total volume of 50 µl. The reaction conditions were 94°C for 2 min, 30 cycles of 92°C for 60 s, 60°C for 60 s, and 68°C for 45 s, followed by a final single extension at 72°C for 5 min. DGGE was performed as described above, except a linear chemical gradient ranging from 30 to 60% denaturant was used. The DNA was visualized, the major bands were excised, and the DNA was diffused as described above. A total of 10 µl of eluted DNA was used as template with amoA primers without the GC clamp for reamplification. The products were purified with the QIAquick PCR purification kit and were cloned as described previously.
Statistical analysis and analysis of DGGE bands.
Data analyses of nutrient concentrations were performed by using SAS (20). Analysis of variances, means, and standard deviations for the individual nutrient in triplicate samples were determined (P = 0.05) to compare the concentration of each nutrient in each sample over the sampling period.
DNA fingerprints obtained from the 16S rRNA banding patterns on the DGGE gels were photographed and digitized by using an ImageMaster Labscan (Amersham-Pharmacia Biotech, Uppsala, Sweden). The lanes were normalized to contain the same amount of total signal after background subtraction. The gel images were straightened and aligned by using ImageMaster 1D Elite 3.01 (Amersham-Pharmacia Biotech) and were analyzed to give a densitometric curve for each gel. Band positions were converted to Rf values between 0 and 1, and profile similarity was calculated by determining Dice's coefficient for the total number of lane patterns. Cluster analyses of the lane patterns were constructed by using the Minitab statistical software (version 13; State College, Pa).
Phylogenetic analysis.
Sequence analyses were done by using the BLAST database (National Center for Biotechnology Information; www.ncbi.nlm.nih.gov). Partial amoA gene sequences were aligned with parts from the complete amoA gene sequences of ammonia-oxidizing bacteria obtained from the BLAST gene bank. Sequence alignments were performed by using the PILEUP program from the University of California-Riverside Genetics Computer Group (GCG programs). Matrices of evolutionary distances were computed by using the Phylip program with the Jukes-Cantor model (8). Phylogenetic trees were constructed and checked by bootstrap analysis (100 data sets) by using the program SEQBOOT. Bootstrap values represented the frequency of resampling that supported a specific branching pattern.
|
|
|---|
![]() View larger version (40K): [in a new window] |
FIG. 1. Percent removal of major chemical components in wetland. Analyses were done on weekly samples. Results are presented as an average removal rate for the 11 months that samples were collected. FP, facultative pond.
|
![]() View larger version (28K): [in a new window] |
FIG. 2. Changes in fecal and total coliform concentrations in and out of wetlands. Values shown here are based on the geometric mean from four samples averaged from each of the sampling months. Concentration of bacteria was determined from 100 ml of samples taken at different locations in the wetlands.
|
Values for the conservative ions determined for wetland influent and treated effluent were used in cross-correlation analyses to calculate the residence time in the wetlands (12). Based on the variation in concentrations of these constituents, the residence time was estimated at approximately 7 days. Cross-correlation analyses were possible based on the variations or spikes in the concentrations of ions that pulsed through the wetlands. The variations in ion concentrations were observable in both the wetlands influent and effluent, allowing for estimation of the residence time based on the time difference from when the spike occurred at the influent and at the effluent.
Bacterial community composition in manure and wetlands.
Analysis of the bacterial community composition by PCR-DGGE was performed on all manure, fecal, and wetland samples collected from December 2000 to September 2001. Banding patterns for the 16S rRNA DGGE-PCR amplicons are presented in Fig. 3A to D. Gels for other sampling times showed similar DGGE patterns (data not shown), while DGGE banding patterns from cow and calf feces were significantly different from each other. The gel images consist of Fig. 3A and part of 3B for fecal and manure samples. The number of bands per lane varied from 10 to 22. Some differences were noted in band position, intensity, and number of bands present in the fecal bacterial DGGE patterns throughout the collection period. Each calf or cow sample had its own unique profile, indicating variation within and between animals, as well as sampling time. However, the DGGE analyses from manure samples demonstrated relatively stable banding patterns throughout the collection period (Fig. 3A and B). Two bands (212 and 213) were common in all samples. The Image Master ID was used to create a database set for these samples. The set contained a total of 46 unique bands based on a synthetic band from the software as the standard band. The bands were found in different combinations among the DGGE banding patterns. Occurrence of different bands from the four gels showed that 24% of the bands occurred in more than 88% of the samples. Pattern comparisons were made among the four gels by using Dice's coefficient. From the database analysis, the Dice values indicated that all samples shared a large portion of the band set (data not shown), while a few had unique bands as indicated by the banding patterns shown in Fig. 3A to D.
![]() ![]() ![]() ![]() View larger version (202K): [in a new window] |
FIG. 3. DGGE analysis of 16S rRNA fragments of total bacterial population from fecal and manure samples. (A) Gel image from samples collected from manure in February (2mns), calf feces in April (4calfF), cow feces in April (4cowF), calf feces in July (7calfF), and cow feces in July (7cowF). (B) Gel image from samples collected from raw pond in February (2r), collected in December from wetland effluent (we), wetland influent (wi), facultative pond (p), and raw pond (r). Manure (9m1-9m3) and fecal samples were also collected in September from cow (9 cow1-3) and from calf (9 calf 1-3). (C) Numbers before letters represent months that samples were taken; W1 and W2, wetlands 1 and 2; I, influent; E, effluent. P and R represent facultative pond and raw pond, respectively. (D) Numbers before letters represent months that samples were taken. Amplified products were separated on a gradient gel of 30 to 70% denaturant. All labeled bands were excised from the gel, reamplified, and subjected to sequence analysis. These reamplification products were cloned and screened as described in the text.
|
![]() View larger version (34K): [in a new window] |
FIG. 4. Cluster analysis of microbial communities generated by the analysis of DGGE 16S rRNA PCR patterns representing the genetic similarity of the microbial community profiles obtained by PCR-DGGE. Data in panel A are derived from the gels shown in Fig. 3A and B (lanes 9m-9 calf1 3). Symbols are as shown in gels. (B) Data in panel B are derived from the gels shown in Fig. 3B (lanes 2r-12r), C, and D.
|
Analysis of predominant bacterial species by PCR-DGGE.
The analysis of predominant bacterial species was carried out with fecal, manure, and wetland samples collected during the study. Bands selected for analysis are shown in Fig. 3A to D. Table 1 shows the prominent bands recovered from the DGGE gel. In Fig. 3A and B all fecal and manure samples generated very complex and unique banding patterns, except for the manure samples that showed very stable banding patterns. In addition to the two prominent bands (212 and 213) derived from fecal and manure samples and present in all other samples, two other bands (4026 and 4228) derived from the raw waste and facultative ponds, respectively, were also found in all samples. The derived sequences from these bands confirmed 212 to be 100% similar to salt marsh clone LCP-89 and 213 to be 100% similar to uncultured Caulerpa texifolia bacterium. Bands 4026 and 4228 derived from the wetlands were 100% similar to uncultured delta proteobacterium and uncultured bacterium MK17, respectively. Most of the excised bands produced legible DNA sequences; the only exception was band 267. This rather diffuse band may have contained DNA from more than one bacteria species. The majority of the DGGE bands showed the highest levels of identity to clones recovered from the gastrointestinal tract. Most of the strains from this study had 99 to 100% similarity levels with sequences from the database (Table 1). Most of the prominent bands retrieved from the fecal and wetland samples (234, 255, 265, and 267) were 99 or 100% similar to bacterial sequences identified by Leser et al. (10) from pig gastrointestinal tract. Sequence identity of most of the bands in Table 1 showed that most of the clones were closely associated with gastrointestinal bacteria from different species or bacteria from sludge or wastewater environments.
|
View this table: [in a new window] |
TABLE 1. Sequence analysis of bands excised from DGGE gels derived from bacterial 16S rRNA extracted from feces, manure, and wastewater from wetland samples
|
subunit of ammonia monooxygenase (amoA). Results obtained from the DGGE analysis of the ammonia-oxidizing bacteria were very reproducible across the sampling points throughout the year (data not shown). DNA over the sampling time was pooled, and a total of nine samples and two positive controls were used for the DGGE analysis (Fig. 5). DGGE patterns from the sampling points showed four distinct dominant banding patterns from the nine samples and the controls. The two control samples migrated to the opposite end of the gel, with Nitrosomonas sp. strain ATCC 19178 migrating to the top portion of the gel and Nitrosospira sp. strain ATCC 25196 migrating to the lower portion of the gel. Sequence analysis also revealed that the bands in the lower third of the DGGE gel were phylogenetically related to Nitrosospira species, while the bands in the upper third of the gel were phylogenetically related to Nitrosomonas species (Fig. 5 and 6). However, most of the bands in the DGGE gel phylogenetically clustered with Nitrosomonas sp. strain ATCC 19178 (002, 004, and 007). These results are in agreement with previous results from DGGE analysis (13) using the same primers. Two out of three Nitrosomonas-like sequences obtained from the wetland samples were grouped into one cluster with Nitrosomonas europaea. The only Nitrosomonas-like sequence which did not cluster within this group was band 006, obtained from the wetland influent, which grouped with the N. mobilis cluster (Fig. 6). All sequences related to Nitrosomonas obtained from the wetland and manure samples belonged to the N. europaea-N. mobilis cluster (Fig. 6). None of the samples obtained from the wetlands and manure clustered with the Nitrosospira control sample (Fig. 6). Band 005 clustered with the soil isolate Nitrosospira sp. strain NAB-8-C11.
![]() View larger version (103K): [in a new window] |
FIG. 5. DGGE analysis of amoA fragments obtained from feces, manure, and wetland samples. Similar samples were pooled from all sampling dates (December 2000 to September 2001) for analysis. Nine samples were used for the final analysis plus the control samples for Nitrosomonas sp. strain ATCC (NS 19178) and Nitrosospira sp. strain ATCC (NTS 25196). F and R represent facultative pond and raw pond, respectively; WT1 and WT2, wetland 1 and 2, respectively; Inf, influent; Eff, effluent.
|
![]() View larger version (39K): [in a new window] |
FIG. 6. Phylogenetic tree constructed for partial amoA gene sequences and aligned by the GCG program from the University of California, Riverside, genetic group. The tree was produced by using a neighbor-joining algorithm. The tree shows the relationship between ammonia-oxidizing bacteria from wetland, manure, and fecal samples and sequences obtained from agricultural soils. The bar indicates an estimated 10% sequence divergence. The cluster nomenclature is that of Purkhold et al. (17).
|
The analyses of Nitrosospira-related sequences showed that the majority of the wetland samples were phylogenetically related to Nitrosospira. These samples clustered into two main groups, with most of them closely related to uncultured bacterium NAB-8-C11. The second group clustered with Nitrosospira sp. strain Nsp12 and Nitrosospira tenuis. To determine how closely related the wetland and agricultural soil clones were, five soil clones from eastern Washington state soils (4) were included in the analysis. The results showed that none of the wetland Nitrosospira sequences were in the same cluster as the sequences from agricultural soils (Fig. 6). All the sequences from agricultural soil were closely related to Nitrosospira briensis and Nitrosospira multiformis, which are known to be some of the most common soil ammonia-oxidizing bacteria (16). The Nitrosospira-like sequences identified in the wetlands, manure, and fecal samples in this study were clearly different from the Nitrosospira-like sequences identified in the agricultural soils under different management regimes.
|
|
|---|
The wetlands had greater ammonia-oxidizing bacteria diversity (especially Nitrosospira sp.) than the other components studied. It has been suggested that the level of ammonia-oxidizing bacteria diversity within a wastewater treatment facility has a major influence on process stability (19); the greater the diversity, the more stable the process. These authors concluded that a treatment plant with greater diversity would cope better with changing conditions. A reduction in the numbers of one organism may not mean process failure, as other organisms better adapted to the new conditions proliferate, resulting in a more functionally stable system. The results from this study support this notion, since a more diverse population of both Nitrosomonas and Nitrosospira was identified, and Nitrosospira was inherently more diversified in the wetlands. If diversity does play a major role in wastewater treatment in the constructed wetlands, the engineering of wetlands to have higher diversity may make processes such as nitrification and denitrification more efficient.
In this study, the wetland effluent is more suitable for on-site reuse and will reduce the amount of contaminants entering groundwater supplies as a result of percolation of washwater stored in ponds and sprayed on disposal lands. In addition, the removal of solids in the irrigation water prevents the formation of a hardened manure crust on the disposal surface, which reduces surface runoff. On-site reuse of the wetland-treated washwater through spray irrigation on disposal lands is optimized by the reduction of organic loads and solids in washwater, allowing more water of higher quality to be used for the irrigation of pastures as it has been adopted in this farm. Improved washwater quality also reduces the buildup of solids in irrigation lines and decreases the need for pump and line service by the dairies.
The removal of the main pollutants from the dairy washwater will have a beneficial effect on the surface and groundwater in the Chino Basin. In turn, this process will benefit the quality of water leading into the Santa Ana River and the Orange County groundwater basin. The wetland project will serve as an innovative model for waste management for the dairy industry and other confined-animal facilities. Providing a cost-effective, low-maintenance process that can be independently built and managed, wetland treatment systems throughout the Chino basin will have a significant long-term effect on the quality of the ground and surface water supply.
.
This research was supported by the 206 Manure and By-Product Utilization Project of the U.S. Department of AgricultureAgricultural Research Service.
|
|
|---|
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»