Previous Article | Next Article ![]()
Applied and Environmental Microbiology, October 2005, p. 6325-6334, Vol. 71, No. 10
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.10.6325-6334.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Cumhuriyet Universitesi, Cevre Muhendisligi Bolumu, 58140 Sivas, Turkey,1 School of Civil Engineering and Geosciences, University of Newcastle upon Tyne, Newcastle NE1 7RU, United Kingdom,2 Chemical Engineering, Department of Food Technology and Environmental Technology, University of Cadiz, 11510 Puerto Real, Cadiz, Spain3
Received 25 January 2005/ Accepted 17 May 2005
|
|
|---|
|
|
|---|
The number of individuals should be the ideal benchmark for quantitative studies. Individual counts can be converted to biomass, biovolume, or proportion of biomass, and results obtained by more indirect methods are typically compared to the number of cells per unit volume (15). Fluorescence in situ hybridization (FISH) represents the "gold standard" for quantification of specific bacterial cells in the environment, against which other methods should be compared. Classical (27) and immunological (20) methods are subject to methodological biases, while nonmicroscopic 16S rRNA-based methods (8, 34) or PCR-based methods (13, 14, 18, 19) deliver a proportion of total cell counts, copy number, or relative signal intensities rather than an absolute number of cells or biomass.
A quantitative method may be evaluated with respect to its precision and its accuracy. Wagner et al. (43) originally evaluated the accuracy of FISH counts of AOB by using cell specific oxidation rates, an approach previously used to show that most-probable-number-based methods underestimate AOB numbers (7, 41). Wagner et al. were able to show that the number of AOB detected could, in principle, account for the nitrification rates observed. However, cell-specific reaction rates are likely to be a crude method for corroborating a quantitative procedure, because the rate will vary with environmental conditions and possibly between taxa. For example, published cell-specific reaction rates in pure cultures of AOB vary by one and a half orders of magnitude (0.9 to 53 femtomoles/cell/hour) (7, 24, 39). Cell-specific ammonia oxidation rates estimated in situ are equally variable, but lower, and range from 0.22 to 2.3 femtomoles/cell/hour (reported values of 2.3 femtomoles/cell/hour [10], 0.63 femtomoles/cell/hour [17], 0.22 femtomoles/cell/hour [43], and 0.25 to 0.97 femtomoles/cell/hour [38]). It is impossible to know if the disparity between the rates measured in pure culture studies and rates estimated from in situ measurements is due to overestimation of the AOB community size in situ or to differences in environmental conditions (rates are likely to be a function of temperature, oxygen and ammonia concentrations, AOB taxa present, and the three-dimensional structure of biofilms or flocs). This critique is not new. When Knowles et al. (22a) first proposed the concept of estimating AOB numbers from cell-specific rates, in 1965, they believed that observed uptake rates could be normalized against known maximum specific uptake rates determined in culture. Writing in 1979, Belser (7) pointed out that this approach could be undermined by a discrepancy between the behaviors of pure cultures and AOB in the environment. Much of what we have learned about AOB in the intervening years would appear to confirm this suspicion.
The precision of AOB enumeration by FISH was not explicitly considered in the earliest literature. However, Schramm and colleagues (38) reported that the Shapiro-Wilks test (typically a test for a normal distribution) "showed an uneven distribution for all data," and they expressed dissatisfaction with the exceptionally large standard deviations. They concluded that their results were only best estimates correct to an order of magnitude. However, high standard deviations and an uneven distribution would be expected if the data were not normally distributed (for example, if the data had a log-normal distribution). Log-normal distributions are associated with entities which grow and die (40).
To overcome, the apparent imprecision of AOB cell counts by FISH, image analysis tools have been used to measure the fluorescence from AOB as a proportion of the fluorescence from the Bacteria overall, and then this is converted to cell numbers by reference to Escherichia coli cells seeded at a known concentration (11). This elegant seven-step procedure yielded a coefficient of variation of 20%. However, the accuracy of the method was been evaluated on the basis of cell-specific ammonia oxidation rates in a single wastewater treatment plant.
Some investigators have reported FISH to be an inferior quantification method. Konuma et al. (23) compared the use of FISH (using Nso190 S-F-bAOB-0189-a-A-19) immunofluorescence and dot blot methods to enumerate AOB. They reported that quantitative FISH in activated sludge was confounded by weak signals, nonspecific binding, and autofluorescence and did not recommend its use. Rittmann and coworkers (34) used FISH and slot blot techniques to quantify AOB in a variety of activated sludge plants and compared their empirical and theoretical biomass estimates. Oligonucleotide probes Nso1225 (S-F-bAOB-1224-a-A-20) and Eub338 (S-D-Bact-0338-a-A-18) were used to detect AOB and Bacteria, respectively. The ratio of AOB to Bacteria obtained by slot blot analysis agreed with theoretical estimates. However, ratios of AOB biomass (obtained by FISH) to mixed-liquor volatile suspended solids (MLVSS) did not agree with theoretical predictions. Ratios obtained by FISH were much lower than predicted. No satisfactory explanation has been offered for this discrepancy. Rittmann et al. (34) tentatively suggested that the majority of AOB were not readily detectable by FISH, possibly because most of the biomass occurred as single cells rather than in microcolonies or failed to hybridize, perhaps due to permeabilization problems.
Rittman et al. (34) put forward an interesting methodology for assessing the amount of AOB that should be present in a system. In essence, they estimated the net production of AOB biomass from reduced ammonia by using the yield to convert ammonia consumed into biomass and simple mass balance concepts to account for losses. They used this technique to compare theoretical and measured biomasses in a number of plants. Thus, if measured and theoretical biomasses corresponded perfectly a plot of theoretical versus measured biomass would have a slope of 1, an r2 of 100, and an intercept of 0. Variation in the intercept would imply a systematic disagreement between the model of Rittman et al. and measured values. Low r2 values or a slope other than 1 would imply site-specific disagreement between theory and measurement. Unfortunately, in the original work the regression line was forced through the origin, and so the ability to assess systematic errors was lost. The accuracy of other published FISH studies cannot be retrospectively evaluated using the approach of Rittmann et al., because they typically focus on a single wastewater treatment plant and do not report the process variables required by the model.
We report a simple method for quantification of AOB by FISH. We explicitly consider the distribution of AOB microcolony sizes, which allows for the estimation of the proportion of AOB not detectable by FISH. The quantity of AOB per unit volume can be expressed as cells per unit volume or biomass. Using our methodology, we show that FISH and process-based estimates of AOB population size in several full- and lab-scale reactors are compatible and that cell-specific ammonia oxidation rates are very variable.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Typical operational parameters of the wastewater treatment plants
|
Sampling.
Grab samples of mixed liquor were preserved immediately in ethanol (sample/ethanol ration, 50:50 [vol/vol]), transported to the laboratory at 4°C, and stored at 20°C until analysis. Prior to FISH analysis, samples were fixed with 4% paraformaldehyde as described by Amann et al. (2). It has been reported in the literature (21, 26) that ethanol fixation can cause sufficient lysis of some gram-negative cells to affect the apparent proportional abundance of certain taxa. However, we did not observe lysis of AOB microcolonies when they were stored for up to 4 weeks.
Oligonucleotide probes.
Oligonucleotide probe nomenclature was based on the Oligonucleotide Probe Database protocol (1). Probes were labeled with fluorescein isothiocyanate, tetramethyl rhodamine isothiocyanate, or indocarbocyanine and were obtained commercially (Genosys, Cambridge, United Kingdom, or ThermoHybaid, Ulm, Germany). A calibration of cell numbers against microcolony dimensions was undertaken using Nso190 (S-F-bAOB-0189-a-A-19), Nsm156 (S-G-Nsm-0155-a-A-19), and Nsv443 (S-G-Nsp-0443-a-S-19) (29) and Nsm641 (S-*-Nsm-0641-a-A-23), a probe specific for the dominant AOB 16S rRNA gene sequence recovered from a laboratory-scale reactor (6). The probes used for counts were Nso1225 (S-F-bAOB-1224-a-A-20), Nso190 (S-F-bAOB-0189-a-A-19), Nsm641 (S-*-Nsm-0641-a-A-23), and Eub338 (S-D-Bact-0338-a-A-18). Negative control analyses using probe nonEub338 (S-D-Bact-0338-a-S-18) were conducted for all samples. It should be noted that Nso1225 has a single-base mismatch with the 16S rRNA of Nitrosococcus mobilis, which is common in reactors treating saline wastewaters, and studies of such systems should use a mixture of Nso1225 and probe NEU (22).
In situ hybridization.
All hybridizations were carried out as follows. Two hundred to 250 microliters of fixed activated sludge samples was placed in a 0.5-ml microcentrifuge tube. The samples were dehydrated in 60%, 80%, and 99.8% ethanol by successive suspension in 1 ml of the appropriate ethanol solution and centrifugation at 3,000 rpm in a microcentrifuge (MSE Microcentaur; MSE UK). After the final dehydration step, the supernatant was discarded and the pellet was resuspended in 38 µl of hybridization buffer (0.9 M NaCl, 20 mM Tris HCl, 0.01% sodium dodecyl sulfate [SDS], X% formamide, where X is the amount of formamide optimal for each probe [Table 2 ]), and 2 µl of labeled probe (50 ng/µl) was added. Negative control hybridizations were done without a probe and with probe nonEub338. Samples were hybridized overnight under the appropriate conditions (Table 2). After hybridization, samples were washed twice in a 0.5 ml washing buffer (20 mM Tris HCl, 0.01% SDS, 5 mM EDTA, and X NaCl, where X is the optimal concentration for each probe [Table 2]) for 15 min at the hybridization temperature, followed by a brief wash in 0.5 ml molecular biology grade water. The samples were centrifuged, and the pellet was resuspended in 10 to 100 µl of filtered, distilled, deionized water. Ten microliters of sample was spotted onto a gelatin-coated slide and allowed to air dry (2). A drop of Citifluor (Citifluor, Kent, United Kingdom) was added to the sample and a coverslip placed over the preparation. The edges of the coverslip were sealed using nail varnish, and prepared slides were stored in the dark at 4°C before viewing. Hybridizations to establish the relationship between cell numbers and microcolony dimensions were conducted in a similar manner except that a complex hybridization buffer [0.9 M NaCl, 50 mM sodium phosphate (pH 7.0), 5 mM EDTA, 0.1% SDS, 0.5 mg of poly(A) per ml, 10x Denhardt's solution] was used (2).
|
View this table: [in a new window] |
TABLE 2. Names, target positions, sequences, and specificities of oligonucleotide probes used during this study
|
![]() View larger version (12K): [in a new window] |
FIG. 1. Fluorescence conferred by probe Nso1225 to whole fixed cells of Nitrosomonas europaea and Ralstonia eutropha at different formamide concentrations. Optimization was done under the protocol of Daims et al. (11) (squares and dashed line), under the protocol used in this paper (diamonds and solid line), and to the nontarget species Ralstonia eutropha (triangles and dotted line). Error bars indicate standard deviations among individual cells in a sample.
|
Chemical and physical properties of mixed liquor.
All physical parameters (mixed-liquor suspended solids [MLSS] and MLVSS) and chemical parameters (ammonia and biological oxygen demand [BOD]) were determined using standard methods (3).
Statistical analysis.
Probability distributions, Anderson-Darling normality tests, analysis of variance, and multiple comparisons of means were undertaken with MINITAB v11 (Minitab Inc., State College, PA). Other statistical analyses were used as described by Sokal and Rohlf (40).
Calculation of area under normal curve.
It is possible to calculate the area under a normal distribution curve based on mean (µ) and standard deviation (
) values. However, if there are some data missing, both the mean and the standard deviation values are distorted and need to be corrected. This may be achieved by an iterative procedure described by Metcalfe (28). The procedure is as follows. The area of the unobserved proportion of the curve is calculated using the values for the "distorted" mean and standard deviation derived from the available data and the area under a standard normal distribution curve. The area calculated is then used to recalculate the values for the mean and standard deviation. The new mean and standard deviation are then used to recalculate the unobserved area under the curve. This procedure is iterated until the area of the recalculated mean and standard deviation reach a fixed value, and the corresponding unobserved area is determined using these values.
Sample size calculation.
The sample size required to achieve a particular power of discrimination was determined using a method we have described previously (12), which was itself derived from a protocol suggested by Sokal and Rohlf (40). Briefly, nested analysis of variance was used to establish that virtually all the observed variation occurred from field of view to field of view (as opposed to sample to sample), and the mean microcolony diameter was identified as the most important variable in the estimation of biomass or cell counts (see below). On this basis we were able to calculate the number of microcolonies that must be counted if we were to have a 95% chance of detecting a given difference in size, significant at the 95% level (Fig. 2A). For example, counting 42 microcolonies ensures that there is a 95% chance of detecting a size difference of 25%, significant at the 5% level. The number of fields of view that must be counted, then becomes a function of the number of microcolonies per field of view (Fig. 2B). To give an 80% chance of detecting a difference of one AOB microcolony per field of view between two samples at the 5% level of significance, it was calculated that a sample size of 46 fields of view was required.
![]() View larger version (12K): [in a new window] |
FIG. 2. A. Number of microcolonies that must be counted to ensure having a 95% chance of detecting a given difference in microcolony diameter significant at the 5% level. The number of microcolonies can be decreased by accepting a marginally lower chance of detecting a given difference. B. Numbers of fields of view (FOV) that must be counted to give an 80% chance of detecting a difference of one AOB microcolony per field of view between two samples at the 5% level of significance. It was calculated that a sample size of 46 fields of view was required.
|
The number of AOB microcolonies per unit volume was determined on the basis of the mean number of microcolonies per field of view, the area covered by the sample spot, the area of one field of view, and correction factors to take account of sample dilution and concentration steps, including the initial dilution in alcohol. The diameters of the microcolonies observed were measured directly and found to be log-normally distributed. Ellipsoidal microcolonies were accounted for by taking the arithmetic mean of the longest and shortest axes of the ellipse. Average microcolony volume was calculated by assuming that AOB microcolonies are spherical and using the geometric mean radius of AOB microcolonies in the equation 4/3
r3. The measurements of mean microcolony abundance and volume were combined to give the mean microcolony volume per unit volume of sample.
The relationship between microcolony volume and cell numbers was determined empirically to provide a calibration curve relating microcolony volume to cell numbers. This relationship was used to convert measurements of microcolony volume per unit volume of activated sludge mixed liquor to cell counts per unit volume.
In principle, cell counts can also be converted to biomass by calculation of cell volume and density (15). In practice, biomass estimates may be obtained with less error by measuring microcolony dimensions. AOB biomass was therefore determined by multiplying the mean volume of AOB microcolonies per unit volume of mixed liquor and the dry weight of cells per unit volume of an AOB microcolony (0.49 g/cm3 [15]). This figure was obtained by using a consensus value of 318 fg of carbon/µm3 of cellular biomass (15) to convert volume to biomass. Carbon accounts for about 50% of cellular biomass (15), and therefore the total density of a bacterial cell is about 636 fg (dry weight)/µm3 of cellular material, which is equivalent to 0.636 g/cm3. However, this value cannot be applied directly to microcolonies, because they contain void spaces. Our observations indicated that most microcolonies appeared to be perfectly packed (Fig. 3). The maximum theoretical packing efficiency for a three-dimensional object (perfect packing) is 77% (for a sphere) to 76% (for an ellipsoid) (36), i.e., there is a 23 to 24% void volume. By allowing for a 23% void space, the biomass density in an AOB microcolony was determined to be no greater than 0.490 g/cm3 (i.e., 77% of 0.636 g/cm3). This is a maximum possible density, because lower packing efficiency in the microcolonies would yield lower densities.
![]() View larger version (25K): [in a new window] |
FIG. 3. Single optical slice through a section of an activated sludge floc (from Wanlip), showing Nso1225-labeled microcolonies of a range of diameters, most of which are of above average (see Fig. 4) for this plant. Bar, 10 µm. Although the cells are very close together, they will not all appear to touch, as even perfectly packed spheroids make contact with adjacent particles at only one point on any given side.
|
![]() |
is hydraulic retention time,
x is biomass residence time, X is biomass, Y is yield, b is the endogenous respiration rate, and fd is the fraction of newly synthesized biomass that is degradable by endogenous decay. The subscripts v and aob represent total bacteria and AOB, respectively. Values for Yaob (0.34 kg volatile suspended solids/kg N), fd (0.8), and baob (0.15 day1) were taken from the literature (16, 34). Xv was measured directly and used to estimate sludge age as described previously (34). The ratio of active AOB to total MLVSS, Xaob/Xv, can be calculated using the equation above and employing known values for the parameters above and measured ammonia removal (
ammonia) and biomass concentration (Xv) (34) and was expressed as a percentage. This is a conservative estimate of the AOB biomass because it is based on ammonia alone; in reality, other forms of reduced nitrogen may become available for nitrification.
Comparison and combination of errors.
Errors were combined using standard equations for the combination of errors and coefficients of variation (30). The combined error for the multiplication of the microcolonies per unit volume of sample, cell count per unit volume, and volume of a microcolony would be (CV12 + CV22)1/2 + 3(CV32), where CV1, CV2, and CV3 represent the coefficients of variation for the microcolonies per unit volume, the slope of the curve relating cell number to volume, and the radius, respectively. The errors in the radius are multiplied by 3 because this value is cubed to calculate the volume; this provides the most conservative estimate of the error.
Cell-specific ammonia oxidation rates.
Cell-specific ammonia oxidation rates were calculated by the method of Daims et al. (11).
|
|
|---|
|
View this table: [in a new window] |
TABLE 3. BOD, ammonia, MLSS, and MLVSS in activated sludge samplesa
|
![]() View larger version (14K): [in a new window] |
FIG. 4. Probability distributions of microcolony diameter for raw measurements (lower panel) and log-transformed measurements (upper panel) for 89 microcolonies from a full-scale wastewater treatment plant (Wanlip), using Nso1225. The putative unobserved fractions are shown as the shaded area in the log-transformed data and were calculated to be less than ca. 3.5% of the total possible observations.
|
0.02%) of the overall biomass. Since the microcolony diameter and volume were log-normally distributed in all the plants observed, we conclude that virtually all the AOB biomass was detected by FISH (Fig. 4).
Estimating cell numbers.
Total AOB cell numbers may be determined by using optical sections obtained with a confocal microscope. It is evident that that there is a relationship between cell numbers and microcolony volume (Fig. 5). However, both cell numbers and microcolony volume were log-normally distributed and a Box-Cox analysis showed that a natural log transformation was appropriate for describing the relationship between microcolony volume and cell numbers. Because log-log transformations can be used to force linear relationships, both raw and transformed data are presented(Fig. 5). For the raw data, the r2 value was 0.89 and both the slope and the intercept were significant (P < 0.001). For the transformed data, the r2 value was 0.81 and both the slope (0.64 [standard error, 0.04; coefficient of variation, 7%]) and the intercept (2.1) were statistically significant (P < 0.001); the residuals were normally distributed, suggesting that the scatter observed is random. The intercept is greater than zero. This could imply some systematic error in the estimation of either cell numbers or microcolony volume. Removing the three obvious outliers (those outside the 95% prediction interval) in the transformed data set changed the r2, slope, and intercept to 0.90, 0.70 and 1.7, respectively.
![]() View larger version (19K): [in a new window] |
FIG. 5. Relationship between microcolony size and AOB cell numbers in activated sludge for the raw (A) and natural-log-transformed (B) data. For the raw data, the r2 value was 0.89 and both the slope (0.66) and the intercept (101) were statistically significant (P < 0.001). However, both sets of raw data were log-normally distributed, and a Box-Cox analysis showed that a natural log transformation was appropriate. For the natural-log-transformed data, the r2 value was 0.81 and both the slope (0.64) and the intercept (2.1) were statistically significant (P < 0.001). We recommend the use of the log-transformed data. CI, confidence interval; PI, prediction interval.
|
|
View this table: [in a new window] |
TABLE 4. AOB cells counts from full-scale reactors
|
![]() View larger version (16K): [in a new window] |
FIG. 6. Variation of the cell-specific ammonia oxidation rate with the concentration of AOB. The plant with the lowest rate is Wanlip, while the plant with the highest rate is Hydburn.
|
![]() View larger version (15K): [in a new window] |
FIG. 7. A. Relationship between theoretical and measured AOB fractions in full-scale activated sludge reactors in the United Kingdom and a bench-scale reactor. The regression line has a statistically significant (P = 0.003) slope of 1.28 (standard error of 0.20) and an intercept (2.8) that is statistically distinguishable from 0 (P = 0.022). The regression line explained 89% of the variation between the two estimates. Xaob/Xv is the proportion of the total biomass measured as MLVSS that is contributed by the AOB. Hydburn was identified as an outlier by using Dixon's test (P < 0.0.05). B. Relationship between theoretical and measured AOB fractions in full-scale activated sludge reactors in the United Kingdom and a bench-scale reactor, but with Hydburn removed. The slope is 1.23 (standard error of 0.16), r2 rises to 94%, and the intercept (2.47) is still significantly different from zero (although only marginally so) (P = 0.034). CI, confidence interval.
|
The combined site-specific error estimates for the number of cells per unit volume were 18% (Wanlip and Chorley), 19% (Preston), 31% (Stoke Bardolph), and 40% (Hydburn and bench-scale reactor). Not surprisingly the higher coefficients of variation were associated with lower cell counts. The number of cells per unit volume of mixed liquor was calculated by multiplying the number of microcolonies per unit volume of mixed liquor (coefficient or variation, < 5%), the mean microcolony volume based on measurement of the geometric mean diameter of the microcolony (mean coefficient of variation, 17%; range, 9 to 22%), and the number of cells per unit microcolony volume (from the slope in Fig. 5) (coefficient of variation, < 7%). Virtually all the error is attributable to the error in the microcolony mean radius, as this error is multiplied by 3 (as the radius is cubed).
The random error in the estimation of the biomass per unit volume is on the order of 39%. Again, much of the error is attributable to the error in the measurement of microcolony diameter, and the site-specific errors vary accordingly from 30% to 46% (30% for Wanlip and Chorley, 31% for Preston, 40% for Stoke Bardolph, and 46% for Hydburn and the bench-scale reactor). The biomass per unit volume was calculated by the multiplication of microcolonies per unit volume (coefficient of variation, < 5%), the mean microcolony volume based on measurement of the geometric mean diameter of the microcolony (mean coefficient of variation, 17%; range, 9 to 22%), and the estimate of AOB carbon per unit volume (coefficient of variation,
25%) (15). The errors due to the variation in the packing efficiency were not included in the analysis because the value employed in the calculations represented a fixed upper value, not a mean. However, since we know that the coefficient of variation in the number of cells per unit volume is about 11% in the size range used in biomass estimates, it appears that this is not a predominant source of error in biomass estimations.
Confocal microscopy versus conventional fluorescence microscopy.
To investigate whether image quality affects quantification, we counted Nitrosomonas spp. in a nitrifying activated sludge plant by quantitative FISH (with probe Nsm156), using CLSM and an epifluorescence microscope. In terms of both microcolony diameter and microcolony abundance, the values from epifluorescence microscopy were significantly lower than those obtained using CLSM (t test, P < 0.06). When these parameters were used to estimate the concentration of Nitrosomonas spp., the value obtained using epifluorescence microscopy was lower than that obtained using CLSM by about fourfold.
|
|
|---|
Belser's (7) critique of the use of cell-specific ammonia oxidation rates to corroborate quantification methods is well founded: it appears that these rates vary by approximately 3 orders of magnitude. This variation is not caused by over- or underestimation of AOB numbers in different plants, because (i) we are able to find more than 95% of the detectable AOB in a plant, and (ii) the variation in AOB biomass and cell numbers between plants is consistent with known variation in the ammonia consumed and plant characteristics. Thus, cell-specific ammonia oxidation rates emerge as a significant process variable. For example, the proportion of the biomass contributed by AOB in Wanlip (6%) is entirely consistent with theoretical estimates (7%), but the cell-specific ammonia oxidation rate observed here (0.03 femtomoles/cell/hour) is an order of magnitude lower than that previously observed in situ (0.22 femtomoles/cell/hour) (43). Thus, the estimate of cell numbers for this plant could have been dismissed as an overestimate. Conversely, we cannot confidently assume that data falling within the published range of cell-specific rates imply accuracy in quantification.
The description of the distribution of the cell numbers is a crucial part of the successful use of FISH to quantify AOB or indeed any other microbial community. Many authors have expressed dissatisfaction with the lack of precision obtained when counting even large numbers of cells by using FISH (26, 38). This has led to the belief that meaningful precision cannot be obtained with manual counting procedures, because of the difficulty of obtaining a sufficiently large sample size. In this and previous studies (12) we have determined the underlying distribution of the data, which has shown that log transformation of the data is required to permit the use of parametric statistical methods and modest sample sizes (without the use of image analysis). Calculation of arithmetic means from nontransformed data which are log-normally distributed (which is commonplace) not only will give large standard deviations but also will provide an erroneous mean value. For example, the data set describing the number of cells per microcolony (Fig. 5) was examined using the Box-Cox method and transformed using natural logarithms to give a back-transformed mean of 232 cells/microcolony and a standard deviation of 2. Using untransformed data, the arithmetic mean of the same data set was 300 cells/microcolony with a standard deviation of 241. Thus, by recognizing the underlying distribution, it is possible to reduce the variance in the data and make valid statistical comparisons using modest sample sizes in conjunction with powerful parametric statistics. A further advantage of examining the distribution of the AOB microcolony size is that we are able to establish the fraction of AOB that are detectable and demonstrate that those which occur as single cells or very small microcolonies represent a small proportion of the total AOB. We are thus able to disprove the tentative hypothesis of Rittmann et al. that discrepancies between modeled AOB biomass levels and those measured by FISH were due to single cells (34).
If quantitative methods are to improve, it is vital that we assess the nature and cause of our errors. A plot of theoretical versus experimental estimates provides a relatively plausible basis for such an assessment. Perfect agreement between a perfect theory and a perfect form of measurement would give us a slope of 1, an intercept of 0, and an r2 value of 100%.
In our comparison of measured AOB abundance and the abundance predicted from Rittman's model (34), the slope of the curve was statistically indistinguishable from 1, which also suggests good agreement between theory and measurement over a wide range of scales, ammonia concentrations, and sludge ages.
The r2 values suggest that between 6% (if Hydburn is an outlier) and 11% (if Hydburn is not an outlier) of the observed variation between plants cannot be accounted for by our method. The former value implies only modest room for improvement, so the Hydburn data are worthy of further consideration, having far fewer AOB than theoretically predicted. This plant was close to failure shortly before the measurements were taken, and thus it is possible that the AOB community was not at equilibrium or not limited by ammonia, that other nitrification processes (e.g., heterotrophic nitrification) were at work, or that the organisms were so stressed that the yield was very low. Any of these possibilities would mean that the abundance of AOB was not described by the theory, which assumes an ammonia-limited system at equilibrium. Alternatively, the probe employed might not have detected all the AOB present in the Hydburn plant. The probe Nso1225 has a mismatch with the 16S rRNA of N. mobilis and so could underestimate the AOB community if this organism is abundant in the treatment plant. Organisms related to N. mobilis have been typically associated with wastewater treatment plants treating saline wastes (22, 37).
The value of the intercept on the y axis of the plot of measured versus predicted AOB biomass tells us whether the theoretical estimates systematically overestimated (intercept of <0) or underestimated (intercept of >0) the amount of measured AOB biomass. The method we present appears to systematically overestimate the amount of measured biomass (or vice versa). This systematic error must, in part at least, represent the assimilation of ammonia by the non-AOB biomass, an important sink of ammonia that is neglected in the original calculations by Rittmann et al. However, other simplifying assumptions include an assumed yield, measuring the removal of ammonia rather than total reduced nitrogen, the use of consensus data on the density of individual bacterial cells (15), and that the AOB were packed in microcolonies with the maximum possible efficiency (the validity of this last assumption is weakest in the largest microcolonies). We wish to draw particular attention to the assumed yield. The AOB biomass estimates vary in a linear manner with this value. Therefore, using a lower yield would also have brought the data closer to an intercept of 0. It is interesting that the yield employed is the maximum theoretical value (35). An authoritative review has suggested that AOB yields are nearly always close to this maximum in all autotrophic nitrifying organisms (32). However, the available empirical values represent data from a limited number of taxa under good laboratory conditions, and yield could vary with environmental conditions. AOB in "real life" might obtain slightly less than the theoretical maximum yield used in our theoretical estimates (especially if subject to stress). The advent of trustworthy molecular tools for the quantification of AOB should allow engineers to ascertain AOB yields and incorporate them into design and management strategies for wastewater treatment plants and to empirically relate them to environmental conditions.
Why did Rittmann et al., using the same approach, fail to find agreement between the theoretical and FISH-based estimates of AOB biomass (34)? This is probably because a conventional fluorescence microscope was used and the assumed value for the density of AOB cells was low (0.1 g/cm3). It appears that a CLSM, or a microscope of equivalent performance, is required for meaningful quantitative FISH in activated sludge and probably other complex three-dimensional environments.
The precision of our manual cell counts (coefficients of variation, 17 to 40%; mean, 27%) is as good, or better, than those previously reported (Schramm et al. [manual counting], 3 to 50% [38]; Wagner et al. [manual counting], 19% [44]; Daims et al. [quantification using image analysis], 20% [11]). These are minimum estimates of the errors, since some sources of error are unreported in those studies. Importantly, our method applies only to organisms that reliably form microcolonies. Thus, although the method of Daims et al. (11) is more complex than our own, their quantification strategy may be more widely applicable. Most of the random variation in our method appears to be attributable to variation in the measurement of the diameter of the microcolony. Thus, if desired, greater precision can be achieved by improving measurements of AOB microcolony diameter. Improving the precision of other elements in the method will yield only limited improvement in the precision of the method.
It would seem that, in general, cell counts will be preferable to biomass estimates, at least until the coefficient of variation in the estimation of the conversion factors can be reduced. The high r2 values for estimated and measured biomasses suggest that our precision is perhaps better than we suggest, probably because actual biovolume-to-carbon ratios are relatively constant even though the estimation of the exact value of such ratios is subject to error.
Quantification is a strategically important aspect of microbial ecology. However, there is a world of difference between a number and the correct number. Important insights and practical applications will be gained if we can improve accuracy and precision in our methods.
One such insight might be that cell-specific ammonia oxidation rates vary by 2 to 3 orders of magnitude. Thus, AOB in some plants may be working 1,000 times harder than AOB in other plants. It could be significant that the highest cell-specific ammonia oxidation rate was seen in a treatment plant that was close to failure (Hydburn). It is also possible that the very low cell-specific ammonia oxidation rate in another plant means that this plant could be run more cheaply. We hypothesize that there is a threshold cell-specific ammonia oxidation rate below which stable performance may be expected. Operating a plant significantly below this threshold could incur needless extra aeration costs (typically the largest or second largest recurring expense in a treatment plant), and operating above this threshold may increase the risk of failure. Thus, the ability to count AOB (and indeed other functional groups) could help those operating biological treatment plants to more rationally balance costs and the risk of failure.
S.J.B. and G.C. thank Shell/NERC and Cumhuriyet University, respectively, for financial support.
|
|
|---|
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»