Previous Article | Next Article ![]()
Applied and Environmental Microbiology, August 2004, p. 4950-4960, Vol. 70, No. 8
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.8.4950-4960.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
and James M. Tiedje1,2,4*
Center for Microbial Ecology,1 Department of Microbiology and Molecular Genetics,2 Department of Civil and Environmental Engineering,3 and Department of Crop and Soil Sciences, Michigan State University, East Lansing, Michigan 48824-13254
Received 25 December 2003/ Accepted 31 March 2004
|
|
|---|
|
|
|---|
Cometabolism has been successfully used for degradation of TCE in bioreactors (12, 35) and in groundwater communities (21, 28) provided the appropriate growth substrate. However, the majority of the studies have focused on the kinetics and efficiency of TCE degradation without considering the microbial communities involved. In only one study did the workers investigate the dynamics and diversity of microbial communities undergoing TCE cometabolism in a test aquifer (18, 19). Understanding the structure and stability of microbial communities is important for the success of bioremediation strategies since particular microbial communities can be responsible for the success or failure of the process. Insights into processes such as community succession are necessary to understand and predict future changes in community structure in natural and managed ecosystems. Often, bioreactor systems are operated for extended periods of time with the assumption that the microbial community structure remains unchanged. However, studies performed with a functionally stable methanogenic reactor revealed that microbial communities can be very dynamic, suggesting that stable function is not always correlated with stable community structure (17). Alternatively, perturbation studies with 1,1-dichlorethene, a very toxic compound (5, 15), have suggested that microbial communities can be quite resilient since the community structure recovered to its original state (19). Hence, there is a need to study the dynamics and stability of microbial communities with respect to their function.
In the study reported here, we compared the structures and functions of two reactor communities, one fed phenol only and the other fed phenol plus TCE for over 2 years of continuous operation. Our objectives were to determine the effect of long-term TCE application on community structure, succession, and reactor performance (performance was measured as rates of phenol and TCE degradation). We used sequencing batch reactors (SBRs) so that the phenol consumption phase was temporally separated from the TCE degradation phase to avoid competitive inhibition for the two substrates. Periodic changes in community structure and function were observed in the SBR that received phenol alone. These changes were not observed in the phenol-plus-TCE-fed reactor, in which long-term exposure to TCE selected for a community that was more stable in terms of structure and function.
|
|
|---|
1.1 liter) was removed before another 12-h cycle began. Suspended solids (200 ml) were removed every other day 1.5 h after phenol injection to maintain a mean cell residence time of 11 days.
TCE and phenol analysis.
Reactor performance was evaluated by calculation of the kinetics for phenol and TCE removal in batch assays. Samples (5 ml) from each reactor were taken 1.5 h after the end of the phenol addition step of the cycle and were transferred to three 20-ml glass vials that were crimp sealed with Teflon-coated butyl rubber stoppers. For TCE removal measurements, different amounts (10 to 50 µl) of TCE in water (
1,000 mg/liter) were injected into the vials, which were then incubated at 120 rpm in a rotary shaker at room temperature. Periodically, 0.1 ml of the headspace gas was removed with a 0.5 ml Pressure-Lok series A-2 gas syringe (Hamilton Company, Reno, Nev.) and injected into a Hewlett-Packard 5890A gas chromatograph equipped with a capillary column (DB624; 30 m by 0.53 mm [inside diameter]), a flame ionization detector, and an electron capture detector. The gas chromatograph was operated isothermally at 90°C with helium as the carrier gas (flow rate, 12 ml/min), and the injection port temperature was set at 250°C. The temperatures of the flame ionization detector and the electron capture detector were 250 and 350°C, respectively. Second-order rate constants for TCE removal, which normalized for biomass, were calculated from the TCE concentrations as described previously (12). Analyses of the headspace gas from the phenol-plus-TCE-fed reactor indicated that the injected TCE was completely degraded before the settling and decanting phases. For determinations of phenol removal rates, different amounts of phenol were injected into vials, which were then incubated as described above. Periodically, samples were removed with a syringe and filtered through 0.2-µm-pore-size nylon syringe filters (Wheaton). Filtrates (2 ml) were collected to measure the phenol concentrations with a high-performance liquid chromatograph equipped with a C18 column (Waters, Milford, Mass.) and a UV detector set to 235 nm. The solvent used was an acetonitrile-water mixture (60:40) at a flow rate of 1 ml/min. Analyses of samples from both reactors after phenol injection indicated that the phenol concentration decreased to zero after the first 2 h of the phenol recharge stage. Zero-order kinetic transformation coefficients for phenol removal were calculated from the phenol concentrations as described previously (12). The reactor biomass, expressed as dry weight, varied over the entire operational period. The means ± 95% confidence intervals for the phenol-fed and phenol-plus-TCE-fed reactors were 1,490 ± 66 and 1,641 ± 58 mg/liter, respectively. The dry weight of the phenol-fed reactor included approximately 10% extracellular polymeric substances, while the phenol-plus-TCE-fed reactor contained 45% extracellular polymeric substances (7).
Reactor sampling and DNA extraction.
Samples (
12 ml) of suspended solids were taken for community analysis at the same time that the kinetic assays were performed. The samples were centrifuged for 15 min at 6,000 x g and 4°C, the supernatant fractions were decanted, and the cell pellets were stored at 80°C. No major wall growth was observed in the reactors.
For nucleic acid extraction, cells were resuspended in TE buffer (10 mM Tris [pH 8.0], 1 mM EDTA [pH 8.0]) and ATL lysis buffer (QIAGEN, Chatsworth, Calif.) before they were mixed with a vortex mixer. Each cell suspension was subjected to two freeze-thaw cycles before it was digested with lysozyme (50 µl of a 100-mg/ml solution; Gibco BRL, Gaithersburg, Md.) and achromopeptidase (10 µl of a 25-mg/ml solution; Sigma, St. Louis, Mo.). The cell lysate was digested with RNase (10 mg/ml; Boehringer Mannheim, Indianapolis, Ind.), followed by proteinase K (45 µl of a 25-mg/ml solution; Gibco BRL) for 1 h at 60°C. Nucleic acids were purified by three phenol-chloroform-isoamyl alcohol extractions and were precipitated overnight at 20°C by addition of 0.1 volume of 3 M sodium acetate (pH 5.3) and 0.6 volume of 100% isopropanol. The DNA pellet was resuspended in 100 µl of sterile deionized water. The purity and quantity of the DNA were determined by UV spectrophotometry at 260 and 280 nm. The DNA was stored at 20°C.
PCR primers and conditions.
Amplification of bacterial 16S rRNA genes (rDNAs) from the community DNA was performed with primer 8F (5'-AGAGTTTGATCMTGGCTCAG-3', where M is A or C) (20), fluorescently labeled at the 5' end with the phosphoramidite dye 5-hexachlorofluorescein (Operon Inc., Alameda, Calif.), and primer 1392R (5'-ACGGGCGGTGTGTACA-3') (3). PCRs were carried out as described previously (9) with the following modifications: 20 ng of template DNA was added to each PCR mixture, and 25 cycles of amplification were performed. Amplification of the phenol hydroxylase genes from the community DNA was performed with primers FD2 (5'-GGCATCAARATYNNNGACTGG-3', where N is a mixture of A, C, G, and T; R is A or G; and Y is C or T) and Phe212 (5'-GTTGGTCAGCACGTACTCGAAGGAGAA-3') (42). These primers produced a 490-bp PCR product from the large subunit (
subunit) of the multicomponent phenol hydroxylase gene family. Primer FD2 was designed by comparative analysis of amino acid sequences of this protein from Pseudomonas sp. strain CF600, Pseudomonas putida strain H, P. putida P35X, Burkholderia cepacia G4, and Burkholderia sp. strain JS150. Primer specificity was tested by using a collection of isolates known to degrade phenol (7) and B. cepacia G4. The PCR products were sequenced and were found to be highly similar to phenol hydroxylase genes from known phenol degraders. PCR amplification was performed by using the thermal cycler and reaction mixture described previously (9), except that 25 pmol of each primer was used and five replicate reaction mixtures were prepared for each sample. The following touchdown thermal profile was used for amplification: an initial denaturation step of 5 min at 95°C, followed by 10 cycles of 94°C for 30 s, 30 s of annealing (in which the temperature started at 58°C and decreased 1°C every cycle), and 45 s at 72°C. In addition, 25 more cycles were carried out with an annealing temperature of 55°C. The quality and quantity of the PCR products were determined by electrophoresis of an aliquot (10 µl) of each PCR mixture on a 2.0% (wt/vol) agarose gel (Gibco BRL) in 1x Tris-acetate-EDTA.
PCR product purification, restriction digestion, and gel separation.
PCR product purification, restriction digestion, and gel separation for terminal restriction fragment length polymorphism (T-RFLP) analysis of 16S rDNAs were performed as described previously (9) with the HaeIII and CfoI restriction enzymes (Gibco BRL).
For the phenol hydroxylase genes, the PCR products from the five replicate reactions were combined and concentrated by roto-evaporation. The concentrated PCR products were separated in a 2.5% (wt/vol) agarose gel (Gibco BRL) in 1x Tris-acetate-EDTA, and the band of the correct size was excised from the gel and eluted in 60 µl of sterile filtered distilled water by using a QIAquick gel extraction kit (QIAGEN). The purified PCR products were concentrated by roto-evaporation and digested with 15 U of MspI and RsaI restriction enzymes combined (Gibco BRL). The resulting fragments were separated in a 4% (wt/vol) Metaphor agarose gel (FMC, Indianapolis, Ind.) in 1x Tris-borate-EDTA buffer at 140 V for 3.5 h at 4 °C and then stained with ethidium bromide (0.5 µg/ml). Gel images were analyzed by using GelCompar 4.5 (Applied Maths, Kortrijk, Belgium).
Analysis of T-RFLPs.
Electropherograms from each sample from both reactors were manually aligned with each other by using Genotyper 2.5 (Applied Biosystems Instruments, Foster City, Calif.). To avoid primer artifacts, fragments smaller than 35 bp were excluded from the analysis. The alignment resulted in a matrix in which different peaks were considered to be indicative of different ribotypes present in the community, and the peak heights were used as a measure of ribotype abundance. After alignment, samples for which the sum of the peak heights was less than 10,000 U were redigested and separated, since data below this threshold are unreliable for statistical analysis (8). Histograms with the abundance of each peak normalized by the total signal in each sample were used to compare the changes in community structure in each reactor. Similarity between samples from the same reactor or between samples from the two reactors was determined by correspondence analysis (CA) (25, 37) with the software package ADE-4 (http://biomserv.univ-lyon1.fr/ADE-4.html) (38). CA was chosen for the analysis of the T-RFLP profiles because this ordination method recovers the underlying structure of the data by analysis of the proportional changes in species abundance, providing a global view of the community structure without a severe effect due to noise. For CA, peak height was used as described above, and the data analysis was performed in a step-wise process. CA was initially performed on all the samples (from both reactors) to compare the two reactors (between- and within-treatment variability). Subsequently, each reactor data set was analyzed individually to evaluate within-treatment variability.
To evaluate the reproducibility of the T-RFLP methodology, nine replicate samples from each reactor (taken after 811 days of operation) and 32 lanes of an internal standard were analyzed by using two similarity indices. Similarity indices were used because they generate direct sample-to-sample comparisons, providing more precise estimates of the variability. The Sørensen similarity index takes into account the presence or absence of peaks, while the Steinhaus similarity index takes into account the abundance of each peak (25). Similarity matrixes were generated for each set of nine replicates and 32 lanes of the internal standard (Gene scan 2500; Applied Biosystems, Alameda, Calif.) with the SIMIL module of R-package 4.0 (http://www.fas.umontreal.ca/BIOL/legendre/) (25). Means, ±95% confidence intervals, and coefficients of variation were calculated from the similarity values.
|
|
|---|
![]() View larger version (46K): [in a new window] |
FIG. 1. Long-term functional responses of the phenol-fed and phenol-plus-TCE-fed reactors. (A) TCE transformation rates. (B) Phenol degradation rates. The arrows indicate the sampling days used for the community structure study.
|
In the phenol-plus-TCE-fed reactor, the phenol degradation rate started at an initial value of 0.5 mg of phenol mg (dry weight)1 h1 and then cycled in the early stage, and there were conspicuous decreases (0.19 and 0.22 mg of phenol mg [dry weight]1 h1) and increases (0.6 mg of phenol mg [dry weight]1 h1) (Fig. 1B). After day 145 the phenol degradation rate did not increase to values as high as those before this time, and the changes that occurred were less pronounced. No major wall growth was observed.
Community structure as determined by T-RFLP analysis.
To determine if the two reactor communities were similar after 30 days of operation with phenol and before TCE was added to one of them, T-RFLP analysis was performed with six different tetrameric restriction enzymes (Fig. 2). Each of the enzymes tested produced almost identical patterns for the two reactors, with only minor differences in the heights of some small peaks.
![]() View larger version (40K): [in a new window] |
FIG. 2. T-RFLP electropherograms of PCR-amplified 16S rDNA bacterial genes from replicate reactor communities before TCE application to the phenol-plus-TCE-fed reactor. Samples were digested with six different restriction enzymes individually.
|
To assess if the T-RFLP method is a suitable method for monitoring population dynamics, a reproducibility analysis was conducted (Table 1). The indices used measured similarity based on the presence or absence of a fragment (Sørensen index) and the abundance of each fragment (Steinhaus index). To have some reference for comparison, 32 lanes from the DNA marker GS2500 were analyzed. For the two indices the GS2500 standard resulted in the highest values, as expected (Sørensen index, 1.0; Steinhaus index, 0.931), and the phenol- and phenol-plus-TCE-fed reactor T-RFLP patterns showed good reproducibility (Table 1). The variability of the community measurements was lower than the variability in previous studies (9, 16, 33, 34), and the percentages of variation ranged from 3.3 to 6.2%. This is an important finding since in previous studies there were some disagreements concerning the reproducibility of the T-RFLP method (9, 33, 34), which were probably caused by the limited number of samples analyzed. The samples from the phenol-plus-TCE-fed reactor had a notable amount of polymer (7) that could have affected DNA extraction procedures and might account for the slightly higher variability in this reactor. Nevertheless, good reproducibility was observed for the nine replicate samples from each reactor, indicating that the samples were homogeneous and representative of the microbial community in each system.
|
View this table: [in a new window] |
TABLE 1. Reproducibility of the T-RFLP dataa
|
![]() View larger version (39K): [in a new window] |
FIG. 3. Changes in bacterial community structure over time. T-RFLP electropherograms were generated from a CfoI digest of the 16S rDNA PCR product. The shaded peaks indicate dynamic ribotypes.
|
![]() View larger version (24K): [in a new window] |
FIG. 4. Community structure and function relationships between reactors. (A) CA of T-RFLP data from reactor communities. Results from independent digests with CfoI and HaeIII were combined for the analysis. Circles and roman numbers indicate periods when the community structures of the two reactors were similar, as follows: cluster I, phenol-fed reactor on days 50 to 198, 567, and 605 to 811 and phenol-plus-TCE-fed reactor on days 50 to 104; cluster II, phenol-fed reactor on days 357 to 541 and phenol-plus-TCE-fed reactor on days 294 to 361, 399, and 414 to 811. T=0, day 0 of each reactor before TCE application to one of them. (B) TCE transformation rates for the two clusters identified by CA. Graph I, cluster I; graph II, cluster II. The dotted lines indicate the average TCE transformation rates (in liter per milligram per day) for the days analyzed. Dim, dimension; NM, not measured.
|
![]() View larger version (99K): [in a new window] |
FIG. 5. Relative abundance of bacterial terminal restriction fragments (T-RF) from reactor samples after digestion with CfoI. The numbers next to the symbols on the right indicate the sizes of the T-RFs (in base pairs).
|
![]() View larger version (19K): [in a new window] |
FIG. 6. CA of T-RFLP data for each reactor community individually. (A) Phenol-fed reactor; (B) phenol-plus-TCE-fed reactor. The numbers next to the symbols under each panel indicate the day periods. The arrows and numbers (1 to 3) in the graphs indicate the order of the changes in community structure. Results from independent digestions with CfoI and HaeIII were combined for the analysis. The circled areas indicate the different clusters. The mean ± 95% confidence interval for the TCE transformation rates for each of the different clusters is indicated next to the circle. Dim, dimension.
|
|
View this table: [in a new window] |
TABLE 2. TCE transformation rate comparisonsa m
|
To investigate changes in community structure at the functional level, phenol hydroxylase genes were amplified by PCR, digested with two restriction enzymes, and separated by electrophoresis (Fig. 7). Cluster analysis of the phenol hydroxylase fingerprints revealed that populations with the gene in the phenol-fed reactor exhibited periods that corresponded approximately to the periods found by T-RFLP analysis (Fig. 7). The fingerprints for days 0 to 198 of the phenol-fed reactor were different from the fingerprints for days 226 to 525, when the intensity of a
135-bp band increased and a band between 89 and 104 bp appeared. Furthermore, the fingerprints from days 630 to 811 were similar to profiles for days 0 to 198, which is consistent with the T-RFLP data (Fig. 5 and 6A). The phenol hydroxylase-bearing populations in the phenol-plus-TCE-fed reactor showed two periods of change (Fig. 7). During the first 50 days the community was stable, and then it changed, with the emergence of a band at
135 bp and the disappearance of two bands at
210 and
230 bp. This profile remained stable until the end of the experiment (day 811). Cluster analysis of the phenol hydroxylase fingerprints from both reactors confirmed that the two communities were similar during the initial days and the later period of the phenol-fed reactor, when the highest TCE transformation rates were observed. The communities were also similar during days 226 to 525 and days 127 to 811 of the phenol-fed reactor and the phenol-plus-TCE-fed reactor, respectively (Fig. 7), which were periods with low or intermediate TCE transformation rates.
![]() View larger version (87K): [in a new window] |
FIG. 7. Community restriction patterns of PCR-amplified phenol hydroxylase genes digested with RsaI and MspI. The number above each lane indicates the time (in days). The bars at the top indicate the samples with similar fingerprint patterns within and between reactors as determined by unweighted pair group with mathematical average cluster analysis. Lanes M contained markers. The numbers next to the arrows indicate the sizes of the fragments (in base pairs) in the marker lanes.
|
|
|
|---|
Community analysis by T-RFLP of 16S rDNAs revealed that there was correspondence between community structure and function in the reactors. Changes in TCE transformation rates for the phenol-fed reactor occurred periodically and corresponded to changes in community structure. CA of all the fragments obtained with two enzymes (HaeIII and CfoI) showed that different periods had particular community structures associated with them, and when TCE transformation rates were high, the community structures were similar (Fig. 5 and 6A). Fragment 367 was present during days 0 to 198 and 605 to 729 in the phenol-fed reactor, when the highest TCE transformation rates were observed. In contrast, fragments 206 and 95 were present during days 226 to 567, when the TCE transformation rates were intermediate or lower. Interestingly, a similar pattern was observed in the phenol-plus-TCE-fed reactor, in which fragment 367 was present when the transformation rates were high (days 0 to 104) and was replaced by fragments 206 and 95 when the transformation rates were low (127 to 811 days). Changes in the distribution of phenol-degrading populations can result from adaptations to phenol degradation or interactions between community members (e.g., predation by protozoans, which were present in both reactors). Mutations in the promoter region of the phenol degradation genes (11), repression by organic acids (4, 31), and activities of unknown gene repressors (6) are events that are possible in a mixed community and could affect the fitness of phenol-degrading populations, causing changes in community structure.
While the T-RFLP data suggest that there was a community structure-function relationship, it is possible that the changes observed at the 16S rDNA level were not related to changes in the functional genes. Analysis of the genes encoding phenol hydroxylase (an enzyme often involved in TCE cooxidation) revealed changes in community structure that corresponded to changes in the phenol hydroxylase fingerprints. Different phenol hydroxylase genotypes were present when the reactors had higher transformation rates (days 0 to 198 and 630 to 811 in the phenol-fed reactor and days 0 to 50 in the phenol-plus-TCE-fed reactor), suggesting that some phenol hydroxylases were more efficient in TCE cooxidation (Fig. 7). Furthermore, some phenol hydroxylase genotypes disappeared shortly after TCE was injected into one of the reactors, suggesting that genotypes that were efficient in TCE cometabolism were also sensitive to TCE toxicity.
Understanding community succession has been one of the most challenging topics in contemporary ecology. Current ecological theories of succession predict contrasting outcomes of the process. Succession, as defined by Clements (13), is an ordered unidirectional process of species replacement that culminates in a stable community (i.e., a climax community). More recent theories, however, propose that during succession, disturbances are frequent and communities never reach a climax stage. Our findings suggest that microbial community succession was continuous in the phenol-fed reactor, while a relatively stable community was formed in the phenol-plus-TCE-fed reactor. Furthermore, CA of the phenol-plus-TCE-fed reactor community revealed that most of the changes in the community occurred during the first 100 days (Fig. 6B) and then the community stabilized for 700 days, suggesting that the community reached a climax stage, as proposed by Clements (13). In contrast, the community in the phenol-fed reactor changed periodically over the 2-year period and never reached a stable phase. Interestingly, the community structure of the phenol-fed reactor during the initial days (days 0 to 198) was similar to the community structure at the end of the experiment (days 605 to 811). The similarity between community structures was determined by observing terminal restriction fragments present in the periods that were apparently displaced during the intermediate days (days 226 to 567). This is unusual in ecology since lost species usually do not return to a community. In our case, species were probably never entirely displaced from the community because they could remain as a refugia and their levels were below the detection limit of our methodologies (17, 27). Nevertheless, it can be hypothesized from the changes in community structure that some microbial communities are dynamic and never reach a climax stage, but rather are shifting mosaics in which changes are unpredictable, as observed in some plant communities (29), and only limited temporal stability occurs.
A remaining question is why the phenol-plus-TCE-fed community was more stable than the phenol-fed community. Cook (14) reviewed modern successional theory and concluded that one of the points of commonality in the new theories is that disturbances are frequent and they likely affect the outcome of succession. It is possible that the changes in community structure seen in the phenol-fed reactor were a response to random disturbances that occurred as part of the normal reactor operation. However, because of the way that the reactors were designed and operated, random disturbance should have affected both communities. Then why was one of our communities more stable? A possible answer to this question lies in the effects of TCE on the community. Toxicity studies with pure cultures have suggested that the epoxide generated from TCE cometabolism damages cellular components, rendering cells inactive (41). Furthermore, recent studies on B. cepacia G4 suggest that this strain secretes toxic soluble intermediates from the degradation of TCE, which affect cells that do not degrade TCE (43). Therefore, it is reasonable to hypothesize that TCE application selected for a community that was more resistant to TCE toxicity, and depending on the protection mechanism, the community could have been more resistant to environmental disturbances. For instance, an increase in extracellular polymeric substance accumulation was observed in the phenol-plus-TCE-fed reactor (7), a phenomenon that resulted in the formation of star-like flocs, which to some extent resembled biofilms. Biofilm communities are very stable and resilient to perturbations, suggesting that an analogous mechanism could have conferred increased stability to the TCE-fed community.
The results of this study suggest that the stability of microbial communities can be dependent upon random disturbances and the selective pressure applied to the ecosystem. More than 2 years of phenol and TCE feeding resulted in a community with reduced but stable TCE degradation rates and a relatively stable community structure. Conversely, long-term feeding of phenol resulted in a different community with higher TCE transformation rates but with continuous and unpredictable changes in community structure and function. These findings suggest that microbial communities in nature in similar niches could have contrasting dynamics depending upon particular physical, chemical, or biological conditions imposed.
This work was funded under grant R-817873-02-0 from the U.S. Environmental Protection Agency and under grant R-819605 from the Office of Research and Development, U.S. Environmental Protection Agency, to the Great Lakes and Mid-Atlantic Hazardous Substance Research Center and by NSF grant DEB 9120006 to the Center for Microbial Ecology.
Present address: Department of Civil and Environmental Engineering, Stanford University, Stanford, CA 94305. ![]()
|
|
|---|
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»