Previous Article | Next Article ![]()
Applied and Environmental Microbiology, April 2006, p. 2322-2330, Vol. 72, No. 4
0099-2240/06/$08.00+0 doi:10.1128/AEM.72.4.2322-2330.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Gerard Muyzer,2
Mark van Loosdrecht,2 and
Lutgarde Raskin1*
Department of Civil and Environmental Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois,1 Environmental Biotechnology Group, Department of Biotechnology, Delft University of Technology, Delft, The Netherlands2
Received 15 April 2005/ Accepted 13 January 2006
|
|
|---|
|
|
|---|
Because of the positive correlation between growth rates and rRNA levels (3, 7, 13), it is widely believed that the results obtained with rRNA-targeted hybridization techniques depend on both the abundance and the activity of microorganisms present in the biomass. In the case of an activated sludge wastewater treatment system, the average growth rate of the biomass is controlled by the solids retention time (SRT) imposed on the biomass by operation of the reactor. Therefore, the average growth rate is the same for all microbial groups present in a properly functioning activated sludge reactor at steady state. Consequently, it should be possible to calibrate the specific rRNA level of a population with laboratory chemostat studies and to extrapolate to activated sludge reactors based on the operational SRT. In contrast to the biomass in a laboratory-scale chemostat reactor, the biomass in most full-scale activated sludge wastewater treatment plants is subjected to temporally varying substrate concentrations due the plug flow or relatively limited mixing characteristics of the reactor or due to operation as a sequencing batch reactor. Such conditions lead to what is often described in the activated sludge literature as feast and famine cycles (2, 5, 31), which are defined as the intermittent presence of electron donors in the bulk liquid (feast period) and the absence of electron donors near the outlet of the reactor (famine period). Limited information concerning the effect of rapidly changing substrate concentrations on the dynamics of the rRNA pool is available. In the few studies in which this question was examined, the workers observed that bacterial cells growing in continuously fed reactors accumulated more RNA when the feeding rate was varied cyclically than when the feeding rate was kept constant (22, 26). Although these studies were performed with relatively high dilution rates compared to the SRTs found in activated sludge systems, it is possible that bacteria growing in activated sludge systems accumulate more rRNA than the amount that would be expected if only the SRT was considered a measure of the growth rate. Thus, to properly interpret microbial population abundance and dynamics, there is a need to ascertain the effect of the feast and famine cycles typically found in activated sludge treatment systems on the rRNA level observed in bacterial cells.
In addition to the mixing characteristics of bioreactors, other factors can affect the dynamics of the rRNA pool of microorganisms growing in activated sludge treatment systems. For example, the ability to modulate the growth rate by coutilization of several substrates (14) or the up-regulation of production of rRNA to counteract the presence of inhibitors of ribosome assembly (17) can potentially alter the size of the rRNA pool compared with the size that would be expected based on the SRT. In the current study, we focused on yet another factor, the capacity to uncouple the uptake of substrate and the formation of new biomass by accumulation of storage compounds, such as poly-ß-hydroxybutyrate (PHB). Determining the ability to accumulate storage compounds has been shown to be important for understanding the growth dynamics of activated sludge systems (2, 10, 31). Bacteria capable of accumulating storage compounds seem to rapidly take up substrate during the feast period and to use the substrate taken up that exceeds the capacity of new biomass production to synthesize storage compounds (8, 29, 31). During the famine period, they use the stored compounds to form new biomass and to meet maintenance requirements. It seems logical to hypothesize that bacteria that are not capable of accumulating storage compounds produce new biomass only during the feast period. Consequently, one could hypothesize that the impact of feast and famine cycles on the dynamics of the rRNA pool would be greatly reduced in bacteria that are capable of accumulating storage compounds.
The effects of feast and famine cycles and of storage metabolism on the dynamics of the rRNA pool were studied using two strains of Wautersia eutropha (formerly Ralstonia eutropha), a wild-type strain and a mutant that is not capable of accumulating PHB. The two strains were each grown in a chemostat and in a semibatch reactor (pulse fed every 3 h). The metabolic responses of the four cultures to pulse-wise addition of substrate were assayed to determine the impact of the ability to store PHB and the impact of the size of the rRNA pool on the metabolic response.
|
|
|---|
Culture conditions.
Microorganisms were cultured either in flasks or in reactors. The growth medium for the flask cultures contained (23) Na2HPO4 · 7H2O (7.9 g/liter), KH2PO4 (1.5 g/liter), NH4Cl (0.3 g/liter), MgSO4 · 7H2O (0.1 g/liter), CH3COONa · 7H2O (1.36 g/liter), and a trace element solution (33) (2 ml/liter). The flasks were incubated at 30°C and agitated on a shaking table at 300 rpm.
Each reactor consisted of a 3-liter chamber (Applikon, Schiedam, The Netherlands) with a 2-liter working volume and was equipped with three baffles and a central six-blade turbine rotating at 800 rpm. The reactors were fed medium containing (23) K2HPO4 (0.8g/liter), KH2PO4 (0.3 g/liter), NH4Cl (0.8 g/liter), MgSO4 · 7H2O (0.2 g/liter), CH3COONa · 7H2O (5.45 g/liter), silicone antifoaming agent (0.1 ml/liter; BDH, Poole, Dorset, England), and a trace element solution (33) (2 ml/liter).
Two reactor configurations were used: chemostat and semibatch. The same dilution rate (1 volume per day, corresponding to a hydraulic retention time of 24 h) was used for both configurations. For the semibatch reactor, one-eighth of the reactor volume (250 ml) was withdrawn (in approximately 3 min) every 3 h, and the reactor was refilled (in approximately 2 min) with 250 ml of fresh medium. The reactors were aerated with air so that the oxygen concentration was at least 75% of saturation for the chemostat and at least 50% of saturation for the semibatch reactor during the steady-state period (see below). The pH was controlled at 7.5 in both reactors by using solutions of 1 M HCl and 1 M NaOH. The influent line of the chemostat was equipped with an aerated flow breaker to minimize discontinuity in feed addition. Finally, the maximum liquid level in the reactors was controlled by a level controller (18), while the minimum level in the semibatch reactor was adjusted by adjusting the height of the withdrawing tube.
Experimental approach and analytical techniques.
Each reactor was monitored by measuring the dissolved oxygen (DO) concentration continuously and by measuring the optical density at 420 nm (OD420) and observing the culture with a phase-contrast microscope daily. Periodically, the purity of the culture was assessed by serial dilution plating on nutrient agar. At the same time, the biomass dry weight was determined by filtration through a 0.22-µm filter and drying the preparation to a constant weight in a microwave. We assumed that a reactor was at steady state when the DO concentration and the biomass level were constant (chemostat) or exhibited the same profiles for each cycle (semibatch reactor).
After steady state was reached, reactor operation was stopped, and each culture was tested by performing a pulse experiment equivalent to one cycle of the semibatch reactor (i.e., one-eighth of the reactor medium was replaced with new medium, and the culture was monitored for 180 min). During these experiments, 13 samples were obtained and were rapidly processed as follows. A portion of each sample was immediately filtered through a 0.22-µm filter and stored at 80°C to allow later analysis of acetate, soluble organic carbon (SOC), ammonium, nitrate, and nitrite contents; subsamples were immediately stored at 80°C for analysis of total organic carbon (TOC), protein, and RNA contents. In addition, a few drops of formaldehyde were added to subsamples kept for PHB analysis, which were stored at 80°C and freeze dried. Finally, the OD420 was determined, and the oxygen uptake rate (OUR) was determined by respirometry using a Jenway (Dunmow, Essex, England) model 9300 DO meter equipped with data acquisition software.
TOC and SOC concentrations were determined by determining the difference between the total and inorganic carbon concentrations using a Dohrman (Tekmar-Dohrmann, Cincinnati, OH) model DC-190 TOC analyzer. The particulate organic carbon (POC) concentration was determined by subtracting the SOC concentration from the TOC concentration. Acetate was measured by gas chromatography by direct column injection (injector temperature, 80°C), using an HP-INNOWAX (cross-linked polyethylene glycol) capillary column (Agilent Technology, Amstelveen, The Netherlands) and a flame ionization detector (temperature, 200°C). The carrier gas was helium, and the following program was used: 60°C for 3 min, increase at a rate of 5°C/min, and 150°C for 5 min (6). The ammonium concentration was determined spectrophotometrically at 623 nm by using the phenate method (standard method 4500-NH3.F) (1). The nitrite concentration was determined at 540 nm by using the colorimetric standard method (standard method 4500-NO2.B) (1). PHB was quantified as described by Smolders and coworkers (27) by first extracting the PHB from a freeze-dried cell pellet using a solution of concentrated HCl, isopropanol, and dichloromethane (1:4:5, vol/vol/vol) and then heating the preparation at 100°C for 2 h. The resulting organic phase was extracted with water, dried with Na2SO4, and analyzed by gas chromatography using a STABILWAX capillary column (Restek, Bellefonte, PA) at 200°C and a flame ionization detector (temperature, 240°C) with helium as the carrier gas. Benzoic acid was used as an internal standard, and the results were expressed as the mass proportion of the dried pellet (fPHB,MASS) (equation 1).
![]() | (1) |
Calculations.
The biomass was defined as the fraction of carbon that was structurally or enzymatically active in bacterial metabolism. Therefore, the highly variable PHB pool was not part of the biomass. The formula CH1.5O0.5 was assumed for PHB (formula weight of PHB [FWPHB], 21.5 g mol1), and the formula for the biomass was determined to be CH1.97O0.54N0.24 with a ash content of 7.3% (total formula weight of biomass with ash [FWH+ASH], 28.0 g mol1). The carbon concentrations for PHB (XPHB) (equation 2) and biomass (XH) (equation 3) were calculated by converting the PHB mass ratio (fPHB,MASS) to a carbon molar ratio using the formula weight of PHB and the total biomass and multiplying by the molar concentration of particulate carbon ([g POC liter1]/[12 g mol1]).
![]() | (2) |
![]() | (3) |
The continuous OUR profile of each semibatch reactor was calculated from the online DO measurements using equation 4. The oxygen rate transfer rate coefficient (kLa) was determined by minimizing the sum of squared differences between the calculated OUR and the OUR measured by respirometry.
![]() | (4) |
t is time between measurements, DO is the DO concentration, DOSAT is the saturation DO concentration, and DO(t) and DO(t +
t) are the two consecutive DO measurements. The conversion rates and the standard errors for acetate, ammonium, and PHB during the feast phase of the pulse experiments were determined by linear regression of concentrations versus time, assuming zero-order rates with respect to the concentrations of reactants. The conversion rates were reconciled to the elemental mass and degree-of-reduction balances (i.e., mass and electron balances were performed) to allow calculation of the biomass formation rate and to test for gross measurement errors or errors in the description of the overall reaction (data consistency). These calculations were performed by the technique of van der Heijden and coworkers (30), using the Macrobal software (12).
For the analysis of the famine phase, it was not possible to assume that the rates were zero order with respect to the concentrations of the reactants. In addition, the sampling frequency was much less than that during the feast phase. Consequently, the analysis had to be adapted to each culture. The analysis of the PHB(+) cultures started by describing the degradation rate for the PHB fraction (fPHB) by nonlinear regression of the PHB fraction over time, computed using the NLIN procedure of SAS/STAT (24). Although several equation were tested, the best fit was obtained with equation 5:
![]() | (5) |
is the minimum PHB fraction of the active biomass (29). After fitting, the PHB degradation rate at 110 min after the addition of acetate (approximately halfway through the famine phase) was estimated using equation 5. Based on the measured OUR at 110 min and the estimated PHB degradation rate, the rate of biomass formation was calculated based on mass and degree-of-reduction balances. For the PHB() cultures, it was assumed that the biomass degradation rate was balanced with the OUR. Although the rate of change of the biomass concentration could be calculated from other conversion rates using the elemental mass and degree-of-reduction balances, an insufficient number of independent conversion rates was determined to allow a test for data consistency. Instead, the data were tested for consistency with the metabolic model (see below).
Metabolic model.
A metabolic model was used in this study to analyze the stoichiometry of the biological reactions. This model has been described elsewhere (29) and was based on seven intracellular reactions describing (i) acetate uptake, (ii) PHB formation, (iii) PHB degradation, (iv) biomass monomer formation, (v) biomass monomer polymerization, (vi) acetyl-coenzyme A (acetyl-CoA) catabolism, and (vii) oxidative phosphorylation. The model assumed that 0.267 mol of CO2 was produced during the formation of 1 mol of biomass carbon monomers from acetyl-CoA (9), while the same reaction required 0.66 mol of ATP (28). The model further assumed that 1.5 mol of ATP was used to polymerize 1 mol of biomass monomers (32). Assuming that (i) the internal pools remained at steady state and (ii) PHB was not formed and degraded at the same time, the model yielded two linear stoichiometric equations (29), one for the feast phase (equation 6) and one for the famine phase (equation 7):
![]() | (6) |
![]() | (7) |
,
, and
are the maximum yield coefficients for biomass from acetate, for PHB from acetate, and for biomass from PHB, respectively, and mS and mPHB are the maintenance coefficients for growth on acetate and on PHB, respectively. The yield and maintenance coefficients are further defined as follows:
![]() | (8) |
![]() | (9) |
![]() | (10) |
![]() | (11) |
![]() | (12) |
) and the amount of ATP required for maintenance (mATP). The stoichiometric equation for the feast phase (equation 6) can be used to analyze the PHB() cultures by considering qPHB = 0. However, the stoichiometric equation for the famine phase (equation 7) cannot be used for the PHB() strain. By assuming that the amount of ATP required for maintenance was derived mainly from oxidative phosphorylation, the following stoichiometric equation was used to describe the famine phase of the mutant strain:
![]() | (13) |
and mATP.
![]() | (14) |
![]() | (15) |
![]() | (16) |
and mATP are the same for the two strains and remain constant for the two culture conditions, the four pulse experiments provide eight combinations of metabolic rates (four for the feast phase and four for the famine phase), these rates can be used to determine the two unknown model parameters (
and mATP) by linear regression of the right terms versus the left terms of the equations. In order to take into account the error associated with the two sides of the equations, major axis regression (16) was implemented in a Microsoft Excel worksheet to evaluate the values of
and mATP. Finally, the regression approach also provides a test for the assumption that there is constant stoichiometry among strains and culture conditions, as strong departures from linearity would suggest that the stoichiometry is not constant. |
|
|---|
![]() View larger version (31K): [in a new window] |
FIG. 1. Concentrations of acetate (mmol C/liter) ( ), PHB (mmol C/liter) ( ), ammonium (mmol/liter) ( ), and biomass (mmol C/liter) ( ) (right axis) during a pulse experiment. (a) PHB(+) in a semibatch reactor. (b) PHB() in a semibatch reactor. (c) PHB(+) in a chemostat reactor. (d) PHB() in a chemostat reactor. Note that the left y axes are broken between 12 and 21 mM.
|
Feast phase.
The rates of acetate, oxygen, and ammonium uptake and PHB formation were determined from measured concentrations. Although the biomass concentration clearly increased during the feast phase (Fig. 1), it was not possible to accurately determine the rate of biomass formation from these data or from the protein concentrations because the error for each measurement (approximately ±8% or ±2 mM C with respect to the biomass concentration) was large compared to the observed increase [approximately 3 mM C for the PHB(+) cultures and 6 mM C for the PHB() cultures]. Thus, the elemental mass and degree-of-reduction balances were used (i) to calculate the biomass formation rates, (ii) to reconcile the measured rates with the balances (Table 1), and (iii) to determine if the reconciled rates were different from the measured rates. The reconciled rates were not found to be significantly different from the measured rates (P > 0.25) (Table 1).
|
View this table: [in a new window] |
TABLE 1. Test conditions, characteristic times, specific conversion rates, and macromolecular compositions of the biomass during the feast phase
|
![]() View larger version (21K): [in a new window] |
FIG. 2. Optical density at 420 nm for semibatch reactors (a) and chemostat reactors (b) and specific OUR (c) calculated from the dissolved oxygen concentration for the semibatch reactors. , PHB(+) cultures; , PHB() cultures.
|
![]() View larger version (13K): [in a new window] |
FIG. 3. Specific growth rate during the feast phase as a function of the level of RNA in the biomass. and , PHB(+) cultures; and , PHB() cultures; solid symbols, semibatch reactor; open symbols, chemostat reactor.
|
![]() View larger version (13K): [in a new window] |
FIG. 4. Degradation of PHB in the PHB(+) cultures during the famine phase for the semibatch reactor ( ) and the chemostat reactor (). The lines indicate the fit obtained with equation 5.
|
Growth stoichiometry and metabolic model parameters.
The constancy of growth stoichiometry between cultures was tested by use of a metabolic model. The rates determined during the feast phase and during the famine phase were introduced into the linearized model equations (equations 14, 15, and 16), and the values for
and mATP were determined by major axis regression (Fig. 5). The values for
and mATP were found to be 2.88 ± 0.17 mol ATP · mol NADH 1 and 0.100 ± 0.042 mol ATP (mol C · h)1, respectively. The regression line satisfactorily described the data, suggesting that the metabolic model parameters for all cultures were similar (i.e., there was constant growth stoichiometry across cultures). Finally, the results shown in Fig. 5 suggest that the measured famine conversion rates were accurate and consistent with the elemental and degree-of-reduction balances if it was assumed that
and mATP are constant.
![]() View larger version (17K): [in a new window] |
FIG. 5. Major axis regression to determine the parameters of the metabolic model using equations 14, 15, and 16. The terms on the left of the equations were plotted on the y axis, while the multipliers of (terms on the right of the equations) were plotted on the x axis. For the regression line: = 2.88 ± 0.17 mol ATP (mol NADH)1, mATP = 0.100 ± 0.042 mol ATP (mol C · h)1), R2 = 0.98. and , PHB(+) cultures; and , PHB() cultures; solid symbols, semibatch reactor; open symbols, chemostat reactor.
|
|
|
|---|
The observed differences in RNA levels between cultures have direct implications for the interpretation of quantitative data for bacterial populations in activated sludge obtained by rRNA-targeted molecular techniques. The data show that reactor configurations and specific metabolic organizations affect the observed rRNA level for a given population, which may explain the wide range of cellular rRNA levels obtained for various populations in activated sludge (19, 20).
Conversion rates and growth stoichiometry.
Two important differences between cultures were observed for the conversion rates during the feast phase (Table 1): (i) the biomass formation rate of each culture was different and (ii) the acetate uptake rate of the chemostat-grown PHB() culture was approximately 33% lower than the rates for the other cultures. In general terms, these differences could be explained by differences in metabolic efficiency (biomass and PHB yields) or by differences in the kinetic abilities of the cellular systems. The possibility that the metabolic efficiencies of the cultures were different was tested using the metabolic model (equations 14, 15, and 16). Figure 5 shows that the metabolic efficiencies for cultures were similar as the combined conversion rates closely followed the expected linear relationship. Consequently, the differences in the conversion rates were more likely due to differences in the maximum instantaneous rates of limiting metabolic reactions.
The differences in the biomass formation rates were found to correlate well with the RNA levels (Fig. 3). Since most RNA is rRNA (3), measuring the total RNA also provides a measure of the level of ribosomes in cells. Therefore, it seems that the maximum instantaneous biomass formation rate during the feast phase is related to the level of ribosomes, which makes the pool of ribosomes likely to be limiting during this period. This is consistent with the near-maximal efficiency of ribosome activity immediately after a nutritional shift-up (4, 15), although bacterial cells typically accumulate more ribosomes than necessary to meet the protein synthesis need under balanced growth conditions (15).
The situation with the acetate uptake rate can be analyzed as follows. Based on the maximum growth rate of both strains (0.381 ± 0.007 h1, determined in four independent batch experiments [data not shown]) and the stoichiometric coefficients determined with the metabolic model (Table 2), a maximum acetate uptake rate for these cultures of 0.673 ± 0.057 mol C (mol C · h)1 was calculated (equation 6) (qPHB = 0). This maximum acetate uptake rate is close to the uptake rates observed during the feast period for the two PHB(+) cultures and the semibatch reactor-grown PHB() culture (Table 1). Thus, it appears that the acetate uptake pathway was fully induced in these cultures. This result is in agreement with the results of other studies that found no relationship between the steady-state growth rates of cultures and the acetate uptake rates after pulse-wise addition of substrate (21, 29). It is also likely that the acetate uptake pathway was fully induced in the chemostat-grown PHB() culture despite the lower measured acetate uptake rate. In the case of this culture, it is possible that the limitation of biomass formation (see above) may have restricted the turnover of acetyl-CoA, which in turn restricted the uptake of acetate because acetate must be charged on free CoA to be metabolized (36). This scenario is consistent with the results of a study of Escherichia coli that reported that there was a restriction in the acetate uptake rate due to a reduced acetyl-CoA turnover rate when there was nutritional shift-up (35). In this context, the PHB synthesis pathway may act as an acetyl-CoA sink to sustain high acetate uptake rates.
|
View this table: [in a new window] |
TABLE 2. Maximum yields and maintenance coefficients corresponding to the values for and mATP obtained
|
versus
x
) was quantified previously (2). In the current study, it was found that producing biomass from acetate via PHB reduced the biomass yield by 4%. Based on previously published data for
, the cost was determined to range from 4% to 10% (2). The data presented in this paper made it possible to quantify the benefit of accumulating PHB in feast and famine cycles for the conditions tested. It was found that the PHB(+) cultures and the semibatch reactor-grown PHB() culture had similar specific acetate uptake rates (Table 1) close to the maximum rates based on the maximum growth rate and yield. Therefore, the advantage of accumulating PHB does not seem to be the ability to develop a higher specific substrate uptake rate. However, only the rate for the chemostat-grown PHB() culture was unable to reach the maximum acetate uptake rate, indicating that repeated exposure to feast and famine cycles was necessary to develop this metabolic capacity. Thus, it seems that the ability of PHB-accumulating organisms to take up acetate is less sensitive to growth history. It is possible to quantify this benefit by realizing that the chemostat-grown PHB() culture took up only 50% of the available acetate in the time that it took the PHB(+) cultures to take up all the acetate, a percentage much higher than the 10% yield cost of PHB metabolism. If these results are extended to the ecology of activated sludge systems, they suggest that in an activated sludge system incoming populations that cannot accumulate PHB are not able to effectively compete upon introduction into the system. Therefore, the chance of becoming established is reduced.
Another competitive advantage of the PHB-accumulating bacteria may be their capacity to take up substrate faster after a period of famine. In the current experiments, a lag before the onset of acetate uptake at the beginning of the feast phase was observed for both PHB() cultures. As the length of the acetate uptake period was found to be between 30 and 35 min, a 5-min lag corresponds to uptake of an amount of acetate that is between 14% and 17% of the total amount of acetate available. This proportion of acetate is also greater than the yield cost of producing PHB (maximum, 10%).
In summary, it appears that the main advantage of bacteria that are capable of accumulating PHB and growing in activated sludge systems is their metabolic responsiveness. Accumulating a certain level of PHB allows the cells to take up acetate at the maximum uptake rate without delay and irrespective of the growth history. In the current experiments, these kinetic advantages corresponded to acetate uptake that was 14 to 50% greater than that in cultures that were not capable of accumulating PHB, values that are much higher than the metabolic cost of PHB storage.
This work was supported by grant BES 97-33826 from the National Science Foundation and by a Paul L. Bush Award (Water Environment Research Foundation). In addition, D.F. was supported by a graduate scholarship from the Fond Québécois de la Recherche sur la Nature et les Technologies, by a graduate travel award from the College of Engineering of the University of Illinois at Urbana-Champaign, and by a dissertation completion fellowship from the Graduate College of the University of Illinois at Urbana-Champaign.
Present address: Department of Microbiology and Immunology, University of British Columbia, Vancouver, BC, Canada. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»