Previous Article | Next Article ![]()
Applied and Environmental Microbiology, March 2006, p. 2163-2169, Vol. 72, No. 3
0099-2240/06/$08.00+0 doi:10.1128/AEM.72.3.2163-2169.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Institute of Food Research, Norwich Research Park, Norwich NR4 7UA, United Kingdom
Received 29 July 2005/ Accepted 23 December 2005
|
|
|---|
|
|
|---|
Though the lag time of a growing bacterial population is a commonly used term in predictive microbiology, the relationship between the population lag time and the distribution of single-cell lag times is not straightforward. Even if ideal conditions are assumed, it is practically impossible to identify the distribution of the division times of single cells from the population growth (3) for simple numerical reasons. However, determining population lag from the distribution of the division times of single cells is not a problem. Therefore, models for the growth of single cells allow descriptions of bacterial growth at a deeper level.
Here, we will model the multiplication of cells by a stochastic birth process based on a set of random variables whose realizations are time dependent. Stochastic models are used to describe those processes that are not completely understood or are too complicated to be described deterministically (13). In our case, the first division times for a sufficiently large number of cells could be measured, but the complexity and the lack of knowledge about the intracellular processes during the lag were the main reasons for turning to stochastic models. The question to be answered is how accurately can we predict the growth of the bacterial culture, including the lag and the exponential phase, after we observe the first few division times of some single cells?
It is not easy to observe the division of single cells in such a way that statistically significant distributions can be derived from the observations. Elfwing et al. (7) described a novel automated method that enabled the user to observe the division times of a large number of individual cells during a relatively long interval. Using this technique together with viable-count growth curves measured by traditional plating, single-cell and population growth parameters were measured simultaneously. The measured distributions for the successive division intervals of the single cells were used to model the growth of the population as a stochastic process. The measured growth of the cell culture and that predicted by simulation were compared. The very good agreement justified the use of the simulation to study (i) the relationship between detection times and single-cell lag times and (ii) the effect of the inoculum size on the population lag time.
Detection times (measured by spectrophotometry) of single-cell-generated cultures are commonly used to estimate the lag times of single cells. The use of detection times is based on the assumption that, if the specific growth rate of the population and the detection level are known, the distribution of detection times can be identified with that of the lag times of the original single cells (8, 9, 11, 18, 23). However, this relationship could be seriously compromised by the variability of the first few generation times, when the cells are not necessarily in the exponential phase yet (15, 21). Our simulations, based on data produced by the flow chamber, will make it possible to study this question.
The effect of the inoculum size on the population lag times has also been studied by many microbiologists, who have arrived at mixed conclusions (6, 10, 17, 24). The effect of the inoculum size on the duration of the lag time has, so far, been reported for situations when either the inoculum was very low (1) or, because of stressing conditions, only a small proportion of the population was able to grow (1, 20, 22). The approach developed here could reliably model the population lag time with very low inoculum sizes.
|
|
|---|
Single-cell growth measurements.
The flow chamber was set up as described by Elfwing et al. (7). The system was initialized by flushing through sterile phosphate-buffered saline (0.2 g of KCl/liter, 0.2 g of KH2PO4/liter, 1.15 g of Na2PO4/liter, 8 g of NaCl/liter [pH 7.3]) at 0.7 ml/min, generating an average linear velocity of approximately 1.1 cm/s. One hundred microliters of the culture was inoculated onto the system. The flow was stopped for 15 min to allow cells to attach to the slide. Growth was initiated by flushing fresh medium (LB medium plus glucose) though the flow chamber at 0.7 ml/min. This flow was sufficient to feed the bacteria and remove one of the newborn daughter cells without removing significant numbers of cells already attached. The temperature in the flow chamber was ca. 25°C for the experiment with relatively long lag times and ca. 32°C for the experiment with relatively short lag times. The temperature was controlled by immersing both LB medium and phosphate-buffered saline in a water bath and measured by a thermocouple applied to the flow chamber outlet.
The cells were observed by dark-field microscopy (Zeiss Standard 25; Oberkochen, Germany) with a 10x objective (Zeiss, Oberkochen, Germany). A high resolution (1,040- by 1,392-pixel) charge-coupled-device camera (CoolSnap Pro cf; Roper Scientific, Trenton, NJ) was controlled by an image analysis program (Image Pro Plus; Media Cybernetics, Silver Spring, MD) that captured an image every 5 min. From the images, the sizes (in pixels) of the individual cells at the different times were recorded. The time of division was identified by the time a sudden drop in pixel size occurred; this sudden drop was attributed to the fact that, right after cell division, one daughter cell was removed by the flow of the fresh medium, while the other remained. The division times were recorded by using an in-house program written in Visual Basic. The time to the first division was assumed to consist of both the lag time and the first generation time. The subsequent generation times were defined as the time intervals between two successive divisions: ith generation time = time to the ith division time to the (i 1)th division (i > 1).
Population growth measurements.
Two batches of seven (0.5-liter) bottles containing 100 ml of LB medium were inoculated with different inoculum sizes (from 0.4 to 106 CFU/ml). One batch was incubated in a shaking water bath at 25°C and the other at 32°C. At each sampling time, bacterial counts were estimated by plating samples onto tryptone soy agar (Oxoid CM131).
Random number generation.
Random numbers had to be generated following the empirical distributions of the single-cell division/generation times measured in the flow chamber, as follows.
First, a histogram of the measured empirical distribution was created. The length of the bins in the histogram was estimated as indicated in reference 12. The length of the bins = L = 1.5 x (Q75 Q25) x n1/3, where Q75 is the third quartile or the value with the position 3 x (n + 1)/4 in the sorted vector of measurements (X1 ... Xn), Q25 is the first quartile or the value with the position (n + 1)/4 in the sorted vector of measurements, and n is the number of measurements.
Intervals were constructed starting from the median. The relative frequency for each interval is estimated as ni/n, where ni is the number of observations inside the ith interval.
Next, to generate random values from the histograms, two uniformly distributed random variables were generated: u1 between X1 and Xn and u2 between 0 and the greatest relative frequency. The value u1 is accepted as a value randomly generated from the empirical distribution only if 0 < fi < u2, with fi being the associated relative frequency of the interval to which u1 belongs.
Stochastic modeling of the population growth.
The distribution of the first division time of the first k generation times measured by image analysis for a single cells, are used to model the growth of the population. The value of k is usually around 3 or 4.
The assumptions made for this purpose are as follows. (i) The first division is the sum of the lag time and the first generation time. (ii) After the kth division, the daughter cells are in exponential growth phase, having the same distribution for their generation times as for the last (kth) observed one. (iii) The distributions of the second and successive generation times (up to the kth) are those measured by image analysis.
With the usual terminology for stochastic processes, the states of the system are represented as (Y1, Y2, ... , Yk), where Y1 is the number of cells that have not divided yet, Y2 is the number of cells at second generation time, and Yk is the number of cells at kth generation time. The initial state of the system, at the starting time, t = t0, is (N, 0, ... , 0), meaning that the N initial cells of the population have not started to divide yet.
To the N initial cells, the first division times are assigned by generating randomly N values, (F1... FN), following the empirical distribution of the first division time. The time t1 (the time at which the state of the system first changes) is that of the cell with the shortest division time: t1 = min(F1... FN).
At t = t1, N 1 cells have not divided yet, while the two daughter cells enter their second generation time. Thus, the state of the system is (N 1, 2, ... , 0). The total number of cells is N 1 + 2 = N + 1.
Two values, SG1 and SG2, are generated randomly from the empirical distribution of the second generation times for the two daughter cells. The next division time for these cells will be t1 + SG1 and t1 + SG2.
Now there are two possibilities: the time t2 of the second division in the population will be the minimum time to divide among the N 1 cells that have not divided yet (case 1), or the time t2 will be the division times of the two newborn cells (case 2): t2 = min(F1, ... FN 1, t1 + SG1, t1 + SG2).
In case 1, if t2 is the first division of a cell, then the state of the system will be (N 2, 4, ... , 0). The next division time of each newborn cell is the sum of t2 and a value randomly generated from the second generation time distribution. The time of the third division in the population will be t3 = min(F1, ... FN 2, t1 + SG1, t1 + SG2, t2 + SG3, t2 + SG4).
In case 2, if a cell from the second generation divides first, then the system goes into the state (N 1, 1, 2, ... , 0). The two newborn cells are at the third generation time. Two values will then be obtained following the third generation time distribution (TG1 and TG2). The time to the third division in the population will be t3 = min(F1, ... FN 1, t1 + SG1, t2 + TG1, t2 + TG2).
Continuing with case 2, the third division could happen to a cell that has not divided yet (case 2.1), a cell at second generation time (case 2.2), or a cell at third generation time (case 2.3).
In case 2.1, if the cell divides for its first time, the state of the system will be (N 2, 3, 2, ... , 0). Two values are randomly obtained from the second generation time distribution for the newborn cells, and the time of the fourth division in the population will be t4 = min(F1, ... FN 2, t1 + SG1, t2 + TG1, t2 + TG2, t3 + SG3, t3 + SG3).
In case 2.2, if the cell dividing is in its second generation time, the state of the system will be (N 1, 0, 4, ... , 0). To assign to the newborn cells the time for the next division, two values are randomly obtained from the third generation time distribution. The time of the fourth division in the population will be t4 = min(F1, ... FN 2, t2 + TG1, t2 + TG2, t3 + TG3, t3 + TG3).
In case 2.3, if the cell that divided was in its third generation time, the system goes into the state (N 1, 1, 1, 2, ... , 0). Two values are randomly obtained from the fourth generation time distribution for the newborn cells. If the fourth and/or the successive generation time distributions could not be measured, random values are generated from the third or last measured generation time distribution. The time of the fourth division in the population will be t4 = min(F1, ... FN 1, t1 + SG1, t2 + TG1, t3 + FG1, t3 + FG2).
The iteration can be continued in a similar manner until the population reaches a given number of cells. The system converges to the state (0, ... , NT), in which the total number of cells, NT, are in exponential phase and their generation times follow the kth generation time distribution.
Simulations.
The algorithm above was implemented in an in-house Excel add-in program written in Visual Basic. To validate the stochastic model, 100 simulated growth curves were generated for each experimentally measured growth curve.
To compare the distribution of the detection time with the distribution of the first division time, 100 growth curves were generated, starting from a single cell
To study the effect of the inoculum size on the lag time, 100 growth curves were generated for each of the inoculum sizes (1, 2, 4, 8, 16, and 100 cells).
Estimation of the parameters of population growth.
The maximum specific growth rate and the lag time were estimated by fitting the model of Baranyi and Roberts (5) to the growth curves.
Comparison of distributions.
The distributions of division times and the intervals between successive divisions were compared by
2 tests.
|
|
|---|
![]() View larger version (21K): [in a new window] |
FIG. 1. Distributions of the single-cell first division intervals measured with a flow chamber (7) by image analysis at 32°C (A) and 25°C (B). Rel. Freq., relative frequency; Standard dev., standard deviation.
|
|
View this table: [in a new window] |
TABLE 1. Lag time, maximum specific growth rates, µmax, and doubling times estimated from the measured and simulated growth curves with different inoculum sizes by fitting the model of Baranyi and Roberts (5)
|
![]() View larger version (24K): [in a new window] |
FIG. 2. Growth curves measured by viable counts (dots) and generated by simulation (continuous lines). The simulation was based on the distributions of the generation times of single cells measured by the image analysis-aided flow chamber technique. The agreement is very good at both 32 °C (A) and 25°C (B) for different inoculum sizes.
|
![]() View larger version (29K): [in a new window] |
FIG. 3. Distributions of the first division times (left panels) and the times to reach 107 cells (right panels) measured from simulated populations growing from one single cell. Simulations are based on the distributions for the single-cell growth parameters measured experimentally in the flow chamber at 32°C (A) and 25°C (B). Rel. Freq., relative frequency; Standard dev., standard deviation.
|
![]() View larger version (50K): [in a new window] |
FIG. 4. Distributions of the lag times estimated from 100 growth curves simulated for each inoculum size. Average values and standard deviations are shown next to each distribution. Growth curves were simulated with the single-cell measurements obtained at 32 °C (A) and 25°C (B) and fitted by the model of Baranyi and Roberts (5).
|
|
|
|---|
Recently, it has been found that the daughter cells inheriting the preexisting cellular extremes or old poles grow slower than those with newly formed poles (25). Stewart et al. reported that the successive generation times for cells inheriting the old pole increase but that they decrease for cells with newly formed poles. In our work, the division times were measured in cells that always inherited the old pole by which they were attached onto the flow chamber. However, the growth of these cells did not slow down in successive generation times but instead accelerated. This apparent contradiction can be explained by the different growth phases of the cells. While Stewart et al. (25) worked with cells growing exponentially, here the cells were in lag phase. As mentioned above, it has already been demonstrated that, after lag phase, the cells do not divide immediately at their maximum potential rate but reach it progressively (16, 19, 21). Hence, the aging process can be observed only after the cells reach this rate. Another explanation could be that the continuous supply of fresh medium in the flow chamber delayed the aging process described in reference 25.
Single-cell studies are often carried out by means of optical density methods (8, 9, 11, 18, 23) where approximately 1 cell per well is inoculated into a multiwell plate, and the turbidity is measured by an automated plate reader. From the detection time, it is a common practice to derive the single-cell lag time by a simple shift calculated from the population growth rate and the detection level. Some concerns have been raised (15, 21) about the perturbation of the relationship between lag and detection times by the variability of the first generation times, which are not even homogeneous. Here, we demonstrated that in the case of the single cells used in our work, the variance of the detection times would be equal to that of the first division times.
We did not see an inoculum effect on the lag time in the experimental growth curves, where the lowest inoculum used was ca. 40 cells per total volume (100 ml). For inoculum sizes of less than 40 cells, the first division and successive generation times measured by the flow chamber were used to simulate the growth of the population. The results showed that an inoculum effect on the lag time can be noticed only at very low inoculum sizes. However, this conclusion should be treated with caution, because the threshold number of cells to detect an inoculum effect on the duration of lag varies according to the distribution of the single-cell lag times. Indeed, differences in the lag time duration have been reported for greater inoculum sizes, but those cells were severely stressed by starvation (1, 10) or heating (24) or the growth conditions were very close to the growth-no growth boundary (20, 22).
The influence of the inoculum size on the population lag can be derived mainly from the distribution of the single-cell lag time. The population lag is shorter with a greater inoculum size. Any physiological activity during the lag time, such as the secretion of growth inducers (14), could contribute to the inoculum size effect, but that is out of the scope of this paper.
The successive divisions of each cell within a population were randomly generated from the histograms constructed by observing the division times of single cells in the flow chamber. We also tested a method in which gamma distributions were first fitted to those measured histograms, and then the single-cell division times were generated in the simulation to follow the fitted gamma distributions. As expected, the results with the two methods were similar. Hence, the developed stochastic process provides an accurate prediction of the growth of the bacterial population, regardless of whether an empirical or a theoretical distribution function is used to generate the single-cell division times. It is the lag time and the first few generation times of the original single cells which are of major importance for predicting the time by which low numbers of pathogenic bacteria grow above a dangerous level. Our modeling approach can therefore significantly improve predictive tools for quantitative microbial risk assessment purposes.
We thank Anders Elfwing and András Ballagi for help with the flow chamber setup and Susie George for help with the experiments.
|
|
|---|
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»