Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2008, p. 286-293, Vol. 74, No. 1
0099-2240/08/$08.00+0 doi:10.1128/AEM.01006-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Laboratory of Microbial Ecology and Technology (LabMET), Ghent University, Coupure Links 653, B-9000 Gent, Belgium
Received 4 May 2007/ Accepted 23 October 2007
|
|
|---|
t(week)] parameter and (ii) the Pareto-Lorenz curve distribution pattern. During the whole sampling period, it was observed that neither of the reactor types maintained a static microbial community and that the SBR evolved more gradually than the MBR, particularly with respect to AOB (i.e., average weekly community changes of 12.6% ± 5.2% for the SBR and 24.6% ± 14.3% for the MBR). Based on the Pareto-Lorenz curves, it was observed that only a small group of AOB species played a numerically dominant role in the nitritation of both reactors, and this was true especially for the MBR. The remaining less dominant species were speculated to constitute a reserve of AOB which can proliferate to replace the dominant species. The value of these parameters in terms of tools to assist the operation of activated-sludge systems is discussed. |
|
|---|
Curtis and Sloan (8) postulated that the relationship between structure and function in a microbial community can be evaluated only through an understanding of the source of diversity from which the community is drawn, named the global reservoir or metacommunity. Up to now, diversity has been mainly measured statically, based on indices such as Shannon and Simpson indices. By adding Lorenz curves, a nonambiguous graphical presentation of species evenness can be provided (31, 35, 42), investigating interspecies abundance ratios and thus the internal community structure. Besides the importance of the community diversity and composition at one time point, the dynamics of microbial communities can play an important role in the functionality of a system (49). Nevertheless, due to infrequent sampling, the dynamics of functionally stable microbial communities has not yet been quantified in terms of a rate of change on a short-term basis (weeks).
Nitrification is the succession of the conversion of ammonium to nitrite (nitritation) by the ammonia-oxidizing bacteria (AOB) and the subsequent oxidation to nitrate (nitratation) by nitrite-oxidizing bacteria (NOB) (14). Ammonia oxidation is often a rate-limiting step in nitrogen removal in wastewater treatment plants (WWTP), as the AOB are influenced by various environmental factors such as light, dissolved oxygen, pH, and substrate concentrations (2, 19, 47, 48) and easily inhibited by various organic compounds or heavy metals (24, 26, 34). Curtis and Sloan (8) postulated that autotrophic ammonia oxidation in WWTP is governed by a small metacommunity, leading to easy replicate sampling of the microbially diverse community, which is thought to be stable. Therefore, the AOB form an ideal model group to study functional and microbial stability.
This paper describes a monitoring study of the physicochemical parameters, the floc structures, and the nitrifying bacterial communities of a sequential batch reactor (SBR) and a membrane bioreactor (MBR) operated in parallel. Both reactors were inoculated with the same sludge from a municipal WWTP and received the same ammonium-enriched water with equal volumetric loading rates. Other environmental parameters were kept the same as well. It was investigated whether the functional stability of the nitrification process corresponded with a static bacterial community. Hereto, the bacterial reactor dynamics were characterized by a rate of change as a function of time,
t(week), and the internal community structure was evaluated by Pareto-Lorenz evenness distribution curves.
|
|
|---|
Physicochemical parameters.
The nitrite and nitrate contents of the effluent were analyzed by ion chromatography as previously described (49). Ammonium (colorimetric), suspended solids (SS), and VSS were analyzed by standard methods (18).
Floc size distribution.
Floc size distribution was measured with a CIS-100 apparatus (Ankersmid, Edegem, Belgium), which combines a laser channel performing size measurements based on the time-of-transition principle with an image-analyzing channel allowing characterization of the dynamic size and shape. Equipment operation was performed according to the procedure described by Govoreanu et al. (17).
FISH analysis.
For FISH analysis, activated-sludge samples were fixed according to the protocol of Amann et al. (1). FISH assays were performed with fluorescence-labeled rRNA-targeted oligonucleotide probes by use of the method of Biesterfeld et al. (3). For the AOB species, the Cy3-labeled probe NSO1225 (36) was applied, and for the Nitrospira spp., the Alexa488-labeled probe S-G-Ntspa-0662-a-A-18 (9) was applied. Hybridizations were performed using the probe mixture and a hybridization buffer with a formamide concentration of 35%. Hybridized samples were visualized with a Nikon Eclipse TE300 epifluorescence microscope equipped with a Bio-Rad Radiance 2000 confocal system. Further image processing was performed using ImageJ 1.36b software freely available at http://rsb.info.nih.gov/ij.
Nucleic acid extraction.
Mixed liquor samples of 50 ml were centrifuged at 3,000 x g for 2 min, and supernatants were discarded. Nucleic acids were extracted from 0.5 g of the centrifuged sludge samples by using a low-pH hot phenol extraction procedure as previously described (44). The DNA concentration was measured spectrophotometrically at 260 nm (NanoDrop ND-1000 spectrophotometer; Isogen Life Science, Sint-Pieters-Leeuw, Belgium) and set to 5 ng µl–1 using DNase- and RNase-free filter-sterilized water (Sigma-Aldrich Chemie, Steinheim, Germany).
PCR amplification.
PCR was performed with primers targeting all bacteria (P338f with GC-clamp and P518r) to obtain DNA amplicons for further analysis of the total bacterial community by denaturing gradient gel electrophoresis (DGGE) (39). To study the AOB community, a nested PCR approach was used. The first round made use of the primers CTO189AB, CTO189C, and CTO653r to amplify the AOB specifically (27). The second round was performed with total bacterial primers (P338f with GC-clamp and P518r). PCRs were performed according to the protocol of Boon et al. (4). In all PCR runs, a positive control (1 µl of purified DNA from strain Nitrosomonas europaea ATCC 19718) and a negative control (1 µl of DNase- and RNase-free filter-sterilized water [Sigma-Aldrich Chemie, Steinheim, Germany]) were incorporated. After each PCR run, the size of the amplicon was verified by running it next to a low-range DNA Massruler (Fermentas, Burlington, ON, Canada) on a 1% agarose gel.
qRT-PCR.
Quantitative real-time PCR (qRT-PCR) was performed on an ABI Prism SDS 7000 (PE Applied Biosystems, Nieuwerkerk a/d Ijssel, The Netherlands). To evaluate the variation of the qRT-PCR technique, all amplification reactions were carried out in triplicate with SYBR green PCR master mix (Applied Biosystems). The following primer sets were used: (i) for the AOB, amoA-1F and amoA-2R (41); (ii) for the Nitrobacter spp., FGPS872 and FGPS1269 (10); and (iii) for the Nitrospira spp., NSR1113F and NSR 1264R (11). An oligonucleotide concentration of 150 nM and a DNA template volume of 5 µl were added to 20 µl of PCR master mix in MicroAmp optical 96-well reaction plates (Applied Biosystems). The real-time PCR protocol for all primer sets was as follows: 2 min at 50°C (for enzyme activation, as recommended by the manufacturer), 10 min at 95°C, and 40 cycles consisting of 1 min at 95°C, 1 min at the respective annealing temperatures, and 1 min at 60°C. The respective annealing temperatures were 54°C for the amoA-1F and amoA-2R primer set, 50°C for the FGPS872 and FGPS1269 primer set, and 65°C for the NSR1113F and NSR 1264R primer set. For the amoA, FGPS, and NSR standard curves, the respective standard sequences originated from Nitrosomonas europaea ATCC 19718, Nitrobacter winogradskyi ATCC 25391T, and Nitrospira moscoviensis NCMIB 13793T. In all experiments, appropriate negative controls containing no template DNA were subjected to the same procedure to exclude or detect any possible contamination or carryover. For the primer sets applied in this research, the respective amplification specificities for sludge samples were previously demonstrated by cloning (15). For each primer set, a standard curve was constructed by diluting a PCR-amplified and cloned standard sequence. Furthermore, melting curves were constructed by constant fluorescent measurement during a final 60°C-to-95°C heating at 0.1°C s–1 and routinely checked to confirm the purity of the amplified products. For the recalculation from gene copy numbers to numbers of cells, it was assumed that there were two copies of the amoA gene per AOB cell and one copy of the 16S rRNA gene per NOB cell (11).
DGGE analysis.
A Bio-Rad DGene system (Hercules, CA) was used for running 8% (wt/vol) polyacrylamide DGGE gels with a denaturing gradient ranging from 45 to 60% (100% denaturing contains 7 mol liter–1 urea and 40% formamide) as described previously (4). The desired DGGE bands were cut out, cloned with the pCR 2.1-TOPO cloning kit (Invitrogen, Carlsbad, CA) according to the manual instructions, and sequenced by IIT Biotech-Bioservice (Bielefeld, Germany). DNA sequence analysis was performed using the BLAST server of the National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov) with the BLAST algorithm and specifically with the BLASTN program.
The obtained DGGE patterns were subsequently processed using BioNumerics software version 2.0 (Applied Maths, Sint-Martens-Latem, Belgium). A matrix of similarities for the densitometric curves of the band patterns was calculated based on the Pearson product-moment correlation coefficients and used to perform moving-window analysis (49) by plotting the correlation between day x and day x – 7. This way, each data point in the graph is in itself a weekly comparison. The
t(week) values were calculated as the averages and standard deviations for the respective moving-window curve data points subtracted from the 100% similarity value. The higher the changes between the DGGE profiles of day x and day x – 7, the lower the moving-window curve data point and the higher the
t(week) values will be.
In order to graphically represent the evenness of the AOB communities, Lorenz distribution curves (31) were set up based on the DGGE profiles as previously described (35). Briefly, for each DGGE lane, the respective bands are ranked from high to low based on their intensities. Consecutively, the cumulative proportions of species are used as the x axis, and the y axis is represented by their respective cumulative proportions of intensities. Mathematically, this yields a convex curve. The more the Lorenz curve deviates from the theoretical perfect evenness line (i.e., the 45° diagonal), the less evenness can be observed in the structure of the studied community. A lower evenness means that a smaller fraction of the different species is present in dominant numbers. In this study, the Lorenz curves were also evaluated based on the Pareto principle (40). According to this principle, the cumulative y axis value (in casu the proportion of intensities) corresponding with the 20% level on the x axis (in casu the cumulative proportion of species) is evaluated. In economics, from which the Pareto principle originates, it was observed that a mere 20% of a country's population possesses 80% of its wealth (40). We refer to Lorenz curves which were evaluated on the 20% level of the x axis as "Pareto-Lorenz curves."
Reproducibility of the DGGE method.
The low-pH hot phenol nucleic acid extraction procedure was successfully applied previously to detect different kinds of ammonia oxidizers and both Nitrobacter and Nitrospira (nitrite-oxidizing) species (38). To test the reproducibility of the sampling, extraction, PCR, and DGGE, triplicate sludge samples were analyzed. Similarity values above 97% (data not shown) were noted between DGGE samples originating from the same reactor on the same day. This demonstrated the high reproducibility of the molecular methods.
Nucleotide sequence accession numbers.
Sequences from the DGGE analysis described above were deposited in the GenBank database under accession numbers EF446325 to EF446336.
|
|
|---|
For both reactors, the distribution patterns of the floc diameters were established at the end of the sampling period (see Fig. S2 in the supplemental material). For the SBR, the floc diameter was restricted to 300 µm, and 50% of the flocs had diameters above 165 µm. The MBR floc diameter maximum was 40 µm, and 50% of the flocs were smaller than 13 µm. Based on the steepness of the floc distribution curves, it was concluded that the SBR flocs were more of the same diameter (between 100 µm and 300 µm) and that only 10% of the flocs had a diameter below 100 µm. On the other hand, all flocs of the MBR ranged between 1 and 40 µm. To substantiate this observation, light microscope images were taken for both reactors on day 84 (Fig. 1A and B). These images supported the activated-sludge floc distribution measurements by indicating large and compact flocs in a clear suspension for the SBR, while the MBR's activated sludge consisted of dispersed microflocs.
![]() View larger version (91K): [in a new window] |
FIG. 1. Light microscopic images (magnification, x400) of the SBR (A) and MBR (B) flocs on day 84, and confocal microscopic images (magnification, x400) of FISH analysis of the SBR (C) and MBR (D) flocs focused on the AOB (red) and Nitrospira (NOB, green) communities on day 84. Scale bars, 50 µm.
|
![]() View larger version (9K): [in a new window] |
FIG. 2. Number of AOB (triangles) and Nitrospira (squares) for the SBR (full line) and the MBR (dashed line) as measured by qRT-PCR. Standard deviations (n = 3) are indicated but are very small.
|
Total bacterial DGGE analysis.
The DGGE analysis of all bacteria for both reactors showed more than 45 bands, indicating a very large number of species (see Fig. S3a and b in the supplemental material). Weekly deviations were established with the moving-window analysis (Fig. 3A). The SBR demonstrated weekly shifts of between 10 and 20% during the first month and changes of up to 10% during the last 2 months. However, the MBR shifted up to 20% during the whole sampling period. Average values for these weekly deviations were calculated and expressed as rates of change [
t(week)]. It was concluded that the total bacterial community changed more rapidly in the MBR [
t(week), 13.4% ± 4.8%] than in the SBR [
t(week), 9.6% ± 4.9%].
![]() View larger version (18K): [in a new window] |
FIG. 3. Moving-window analysis based on DGGE profiles for the SBR (full line) and MBR (dashed line) of the total bacteria (A) and the AOB subgroup (B). Each data point in the graph is in itself a weekly comparison, as it represents the correlation between the samples of day x and day x – 7.
|
|
View this table: [in a new window] |
TABLE 1. The 12 AOB DGGE bands with corresponding closest matches via the NCBI database and indication of similarity percentages and sources of origina
|
t(week) values were calculated. The SBR had a
t(week) value of 12.6% ± 5.2% and the MBR a
t(week) value of 24.6% ± 14.3% for their respective AOB subgroups. To asses the interspecies abundance ratios of the AOB, Pareto-Lorenz curve distribution patterns of the AOB DGGE profiles were plotted based on the numbers of bands and their intensities (Fig. 4). Over time, it was observed for the SBR reactor that 20% of the bands (number based) corresponded with 35 to 45% (on average 40%) of the cumulative band intensities. For the MBR reactor, 20% of the bands corresponded with 40 to 70% (on average 55%) of the cumulative band intensities. Thus, for both reactor types it was observed that overall only a small group of AOB species was numerically dominant and that this dominance was even more pronounced for the MBR.
![]() View larger version (13K): [in a new window] |
FIG. 4. Pareto-Lorenz distribution curves based on AOB PCR-DGGE analysis on days 1, 28, 56, and 84 of the SBR (A) and the MBR (B). The dashed vertical line at the 0.2 x axis level is plotted to evaluate the range of the Pareto values.
|
|
|
|---|
The SBR maintained complete nitrification during the entire experimental run. The MBR showed prolonged small nitrite effluent concentrations (0 to 2 mg NO2–-N liter–1) starting from around day 40, indicating incomplete nitratation. This lack of tight coupling between nitritation and nitratation could be linked with a decrease in the MBR sludge floc size (Fig. 1; also see Fig. S2 in the supplemental material). Both floc size distribution and light microscopic analyses indicated floc diameters in the SBR more than 10 times larger than those in the studied MBR, which had a microfloc structure. The lack of good floc structuring in the MBR resulted in sludge bulking, as supported by an SVI of 205 ml (g SS)–1 on day 84. Previously, floc size was reported to play an important role in maintaining complete nitrification, and MBR flocs are known to be smaller than flocs of a conventional activated-sludge WWTP (33), probably due to the lack of a settling and selection phase. In this research, confocal FISH analysis demonstrated the structural similarity of the AOB and Nitrospira clusters, i.e., distinct clusters of AOB and NOB in the flocs, for both SBR and MBR reactor sludges (Fig. 1C and D). Conclusively, it was demonstrated that flocs smaller than 40 µm in diameter can still contain distinct clusters of AOB and NOB but that larger floc sizes favor nitritation-nitratation tight coupling.
In this research, fairly comparable amounts of AOB and NOB were detected in both reactor types. Per mg VSS, a constant level of ca. 6 to 7 log units of AOB (the temporary 2- to 3-log unit drop was not taken into account) and 8 to 9 log units of Nitrospira cells were observed (Fig. 2). With qRT-PCR, AOB are reported to range in number between 5 and 8 log units and Nitrospira cells between 4 and 8 log units per mg VSS in WWTP samples (20, 21, 30). Comparing AOB (also amoA gene based) and Nitrospira cell numbers in the same WWTP samples, Harms et al. (21) reported ca. 10 times higher Nitrospira levels, whereas in the study of Layton et al. (30), the Nitrospira levels ranged from 100 times higher to 2 times lower than the AOB levels. In this research, it was observed that (i) even at startup, AOB values were evaluated as rather high and Nitrospira levels as very high, (ii) no further rise of either subgroup was noted though selective growth substrates were administered, and (iii) internal community dynamics were present. Therefore, this could indicate the existence of a maximal support capacity of AOB and NOB for this municipal WWTP sludge when a volumetric loading rate of 100 mg N liter–1 day–1 was administered; moreover, this finding was independent of the applied reactor configuration.
A static community composition was not observed in the SBR or the MBR. Nevertheless, the stability of the nitrification reaction was maintained. The rate of change of the bacterial communities was visualized by moving-window analysis based on the DGGE profiles (Fig. 3). The moving-window analysis technique has already shown to be useful for monitoring WWTP community shifts (49). In this research, the rate of change was quantified by introducing the "change over time" or "
t value," reported as an average value with its standard deviation. To enable comparison of
t values, the time window is indicated between brackets. Here,
t(week) values were applied to characterize short-term changes. The SBR was observed to change more gradually than the MBR both for the total bacterial (SBR, 9.6% ± 4.9%; MBR, 13.4% ± 4.8%) and for the AOB (SBR, 12.6% ± 5.2%; MBR, 24.6% ± 14.3%) levels, probably due to their different floc characteristics. The operational value of this parameter needs to be corroborated by studies on other and larger reactor types. However, it is of interest that a pragmatic analysis of the DGGE data as performed by moving-window analysis can visualize large differences in microbial community dynamics between two reactors inoculated with the same sludge.
To interpret the interspecies abundance, Pareto-Lorenz evenness distribution curves were plotted. These curves were based on the numbers of bands and the band intensities of the AOB DGGE profiles as described by Mertens et al. (35). It should be noted that the introduction of bias due to preferential amplification is considered as a possible disadvantage when applying conventional PCR (45). However, according to more-recent studies, this bias of preferential amplification may be overestimated, and the band intensity of 16S rRNA genes in a DGGE gel may correspond at least semiquantitatively with the abundance of the corresponding species (22, 23). Nevertheless, extra care should be taken when using nested PCR with primers having specificity issues, as nonspecific amplification might influence the intensity of specific bands. Therefore, it was decided in this research to work with the relative instead of the absolute band intensities. Over time, it was observed for the SBR reactor that 20% of the bands (number based) corresponded with on average 40% of the cumulative intensity of the bands. For the MBR reactor, 20% of the bands corresponded with on average 55% of the cumulative band intensity. This internal structuring of the highly dynamic AOB community indicated that at each moment only a small group of species played a numerically dominant role and that this dominance was even more pronounced for the MBR. The remaining less dominant DGGE bands are speculated to constitute a reserve of Nitrosomonas species which can proliferate to replace the dominant species. Just as for the rate of change, for this parameter as well more data need to be collected for other and larger reactors in order to substantiate the value for practicing this approach in interpreting the microbial ecology of activated-sludge communities.
Up to now, shifts in bacterial compositions have been reported on one hand as a result of changing environmental factors influencing the investigated processes (6, 46); on the other hand (i.e., under constant environmental factors), the sampling occurred infrequently (12), without visualizing short-term shifts. The DGGE monitoring of this study also reported dynamics of the relative interspecies abundance, without measurable environmental changes and based on weekly sampling. The latter observation was constantly noted, as visualized by the DGGE fingerprints and moving-window analysis. This correlates with findings of a study on anaerobic reactors performed earlier (12) and of research on activated-sludge SBRs (25), but our research examined the shifts on a short-term basis (2 to 3 weeks). Furthermore, it should be noted that this study was performed on 25-liter lab-scale reactors in order to control for most of the environmental factors influencing the nitrification process. Natural ecosystems as well as full-scale sewage treatment plants have larger dimensions and are often confronted with newly arriving allochthonous organisms. As it still remains unknown whether these effects lead to a more stable or more dynamic community structure, extrapolation of the observed results to natural ecosystems was not intended in this research phase.
The difficulty of having several repeats during such long-term experiments is a common problem when studying wastewater treatment systems. Often, single reactors are set up and compared over longer time periods (7, 16, 28, 29, 38, 50). In another study, we ran the reactor setup in triplicate for more than 100 days, and both functional and microbial community structure dynamics were highly reproducible (unpublished data). These observations are in agreement with other studies showing minor differences between appropriately repeated lab-scale setups (5, 25, 32, 43).
As observed in this research, floc structure is an important parameter to obtain complete nitrification. However, stability of the microbial community is not essential to maintain the functional stability of the nitrification process. Moving-window analysis and the derived rate of change values were demonstrated to be valuable tools for monitoring microbial community dynamics. Finally, Pareto-Lorenz curve distribution patterns facilitated the visualization of the evenness of the bacterial communities. This way, it was demonstrated that only a small group of AOB species played a numerically dominant role in the studied nitrifying reactors. The remaining less dominant species were speculated to constitute a reserve of AOB which can proliferate to replace the dominant species.
We thank Greet Van de Velde, Aleksandra Ziembiñska, Tijs Decuypere, Veerle De Schepper, and Winnok De Vos for their practical assistance in the completion of this research. Furthermore, Peter Aelterman, Tom Defoirdt, and Joke Geets are gratefully appreciated for their suggestions to improve the manuscript.
Published ahead of print on 2 November 2007. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»