|
|
|
Scientific
Publications - Work Done by Microbiology Reader
Journal of Microbiological Methods, 1999, Oct, 38(1-2), 137-146 Estimation of temperature dependent growth rate and lag time of Listeria monocytogenes by optical density measurementsJean-Christophe Augustina, Laurent Rossob and Vincent Carliera a Service d’Hygiène et Industrie des Denrées Alimentaires d’Origine Animale, Ecole Nationale Vétérinaire d’Alfort, 7 avenue du Général de Gaulle, F-94704 Maisons Alfort Cedex, France b Pôle de Compétence en Sécurité des Aliments, Groupe Danone, Z.I. Le Teinchurier, F-19100 Brive-la-Gaillarde, France Received 12 April 1999; revised 1 July 1999; accepted 8 July 1999. Available online 17 September 1999. ABSTRACT An automated turbidimetric system, Bioscreen C, was used to monitor growth of
ten strains of Listeria monocytogenes at different temperatures. Several
methods for estimation of maximum specific growth rate ( Author Keywords: Automated turbidimetry; Lag time;
Listeria monocytogenes; Maximum specific growth rate; Predictive
microbiology.
1. INTRODUCTION Development of models in predictive microbiology requires large quantities of data. These data are generally obtained from viable count (VC) growth curves. Accurate estimations of growth parameters (lag time and maximum specific growth rate) need a large number of viable count determinations (Bratchell et al., 1989) which are labour and time-consuming. The measurement of optical density (OD) with automated turbidimetric systems provides for liquid and translucent media a rapid and inexpensive method to monitor bacterial growth. Turbidimetry is then an attractive method which has recently been used to acquire growth parameters ( McClure; Young and Barbosa et al., 1994 Jones; B; Br; Stephens; Arino and Cheroutre). There are, however, some problems with this technique and different methods have been proposed to derive growth parameters from OD measurements. For estimation of growth rates, authors frequently recommend the use of classical growth models on OD data from a knowledge of the relationship between OD and VC. It is generally assumed that there is a proportional relation between OD and VC (Corman; Young and Barbosa et al., 1994 Dalgaard; B; B; Br; Arino and Cheroutre) but quadratic ( McClure et al., 1993) and cubic ( Stephens et al., 1997) calibration equations has been proposed as well. Growth models (e.g. logistic, Gompertz) are used on the whole growth curves ( Corman; McClure and Barbosa et al., 1994 Dalgaard; B; B; Stephens; Arino and Cheroutre) or only in the exponential growth phases ( Young and Foegeding, 1993; Br; Stephens and Arino). Lag times has been estimated B and Cheroutre and Hudson and Mott (1994) by applying the Gompertz equation to OD growth curves for which the inoculum level was sufficiently high to produce an initial OD above the detection threshold of the turbidimeters used. When inoculum size is such that the detection threshold is not reached, the initial cell count and the calibration equation must be known. McClure et al. (1993) used the logistic function constrained at the initial cell counts to fit OD growth curves while other Br; Stephens and Arino) estimated lag times from the intersection between the extrapolated straight line, the slope of which is the estimated growth rate, at a point of visible growth and the log of the initial cell count horizontal. Whatever the method used, it seems that the growth parameters estimated from OD and those estimated from VC are unequal. Thus, Dalgaard et al. (1994) showed that growth rates estimated from OD data were smaller than those obtained from VC data. Hudson and Mott (1994) showed with Pseudomonas fragi that lag times calculated from OD data were shorter than those derived from VC data. Similarly, McClure et al. (1993) supposed that their lag times and/or doubling times estimates were shorter in the system they had used. In the present work, lag times and growth rates of ten strains of Listeria monocytogenes were estimated at different temperatures using both OD measurements and VC data, a novel method involving OD data and a few VC data is proposed to accurately estimate growth parameters.
2. MATERIALS AND METHODS 2.1. StrainsTen strains of Listeria monocytogenes were used, Table 1 indicates their reference and their origin. The strains were maintained by monthly transfers on tryptone soya agar (Oxoid, Unipath, Ltd., Basingstoke, Hampshire, UK) slopes stored at 4°C.
2.2. Growth experimentsL. monocytogenes strains were subcultured onto tryptone soya (Oxoid)
plus 0.6% yeast extract (AES, Combourg, France) agar slopes at 37°C for 24 h and
then into tryptone soya (Oxoid) plus 0.6% yeast extract (AES) (TSYE) broth at
30°C for 24 h. Inocula were then diluted in TSYE broth to obtain concentrations
of approximately 108 and 106 cfu×ml−1, and a
series of five half-dilutions was prepared from this last dilution. Aliquots of
400 The automated turbidimeter, Bioscreen C, was used to monitor the growth of L. monocytogenes by reading OD at a wavelength of 600 nm at regular time intervals. OD growth kinetics were constructed by ploting the OD of suspensions minus the OD of the non-inoculated medium vs. the time of incubation. In between readings, wells containing the dilution of about 106 cfu×ml−1 were sacrified to perform VC of the suspensions. VC growth curves were generated with 13–19 points. 2.3. Variance stabilizing transformations for VC and OD dataStabilization of the variance of the VC data was done by using the usual
logarithmic transformation (Dalgaard et al., 1994). The Box–Cox method ( Box et
al., 1978) was used on two sets of 90 replicates at 4.0 and 28.7°C to find the
variance stabilising transformation, s, for the OD data. The log of the
standard deviation of replicates was plotted against the log of the average OD
and the slopes ( 2.4. Estimation of growth parameters from VC growth curvesGrowth curves were fitted using the logistic growth model with delay, i.e., with a breakpoint at the transition between the lag and the exponential phase (Kono; Baranyi and Rosso):
where x(t) is the bacterial concentration (cfu×ml−1)
at time t (h), x0 is the initial bacterial
concentration (cfu×ml−1), xmax is the maximum
bacterial concentration (cfu×ml−1), lag is the lag time (h) and
This function, called f, was fitted to the logarithm of VC. VC maximum specific growth rates and lag times were estimated at 4.0, 7.8, 11.6, 15.4, 24.0 and 35.4°C. 2.5. Calibration of optical densities against cell concentrationOD was calibrated against VC in a given range of OD values using equation:
where pn is a polynomial of order n. In order to know what model between cubic, quadratic or linear was the most appropriate to calibrate OD data, F tests were used to compare nested models of degree n and n−1 (Bates and Watts, 1988). Calibrations were done with OD and VC data obtained from growth of suspensions with initial cell number of approximately 106 cfu×ml−1. 2.6. Estimation of growth rates from OD growth curvesOD growth rates were estimated at 4.0, 6.0, 7.8, 9.7, 11.6, 13.5, 15.4, 19.2, 24.0, 28.7 and 35.4°C. 2.6.1. From suspensions with initial OD above the bioscreen C’s detection thresholdThe model used to fit these growth curves, called model GP, was the modified Gompertz equation, g (Zwietering et al., 1990), which gives with calibration equation and variance stabilizing transformation:
where
ODa is the OD from suspensions with initial OD above the detection threshold, and A is the logarithmic increase of bacterial population (ln(cfu×ml−1)). 2.6.2. From suspensions with initial OD below the bioscreen C’s detection thresholdFrom the suspension with an initial cell count of approximately 106 cfu×ml−1 and the five half-dilutions series, generation time (Tg) was defined as the time separating two successive curves when OD was in the range of 0.010–0.600. This technique is valid assuming there is no effect of inoculum size on growth parameters. Maximum specific growth rate was derived from the average generation time achieved with the six growth curves using:
In the second method, the whole OD growth curve was fitted using the second term of the f model (Eq. (1)) with calibration equation and by constraining the initial OD as being pn(x0), with x0 known. This model, called W, is described by the equation:
where ODb is the OD from suspensions with initial OD below the detection threshold. In the third technique an exponential model with delay, called E, was used only on the exponential growth phase. Its equation is:
2.7. Estimation of lag times from OD growth curvesLag times were estimated from suspensions with an initial OD above the Bioscreen C’s detection threshold with model GP (Eq. (3)), and from suspensions with an initial OD below the detection threshold with models W and E ((6) and (7)). A new method consisted in calculating lag times by using the f model
(Eq. (1)), knowing
This model was called M. 2.8. Variance stabilising transformations for growth parametersTo compare maximum specific growth rates and lag times achieved with the different techniques, we had to stabilize the variances of these parameters. We chose for maximum specific growth rates the square root transformation which seems the most appropriate (Ratkowsky and Alber Schaffner and Ratkowsky) and the logarithmic transformation for lag times ( Ratkowsky and Alber Zwietering and Delignette). 2.9. Validation of the method selected for the acquisition of growth parameters from OD dataValidation of the best method was done by predicting bacterial growth from estimated growth parameters and by comparing these predictions with observed VC growth curves of L. monocytogenes Scott A incubated in TSYE broth at 6°C. The validation was done with eight growth curves achieved from different pre-incubation conditions applied after the culture prepared in TSYE broth at 30°C for 24 h: pre-incubations at 12°C for 10, 22 and 32 h; pre-incubations at 37°C for 1, 3 and 5 h and pre-incubations at 4°C for 52 days and 10 months. 2.10. Model fitFits were performed by linear or non-linear regression using the least squares criterion (Box et al., 1978). Estimation of parameters was carried out by minimizing the sum of the squared residuals (SSR) where SSR is defined as follows:
where n is the number of data points. The minimum SSR values were computed with the REGRESS and NLINFIT subroutines
of
3. RESULTS AND DISCUSSION 3.1. Determination of the variance stabilizing transformation for OD dataThe slopes of the regression lines of the log of the standard deviation of
the 90 replicates against the log of the average OD gave
3.2. Calibration equation of OD against VCCalibration equations were tested in the OD range of 0.0–0.6 for the ten strains because there were breakpoints at higher OD exhibiting an OD remaining at its maximum or decreasing slighty even though VC were increasing. The best fits of the data were achieved with polynomials of order 3 (P<0.001). However, parameter values of calibration equations were significantly different from one temperature to another and these values were not linked to the temperature. Furthermore, by using the average values of parameter estimates at the different temperatures, the best fits of the whole data were obtained with polynomials of order 1. We then used, for calibration equations, polynomials of order 1 with constant term equal to 0 because confidence intervals of this parameter always included this value. Eq. (2) can then be written in the OD range of 0.0–0.6:
where the k-values range from 2.85×10−10 to 4.46×10−10 for the ten strains. 3.3. Growth rates3.3.1. Estimation of growth rates with the dilution series methodThe hypothesis of no effect of inoculum size was checked by observing an almost constant time between successive OD growth curves for seven ten-fold dilutions (Fig. 2).
For all strains, the maximum specific growth rates obtained with the half-dilution series method were equal to those obtained with VC growth curves. However, this technique required several OD growth curves to obtain confident estimations of generation times. Indeed, the time separating the two successive growth curves was relatively variable (data not shown) and six successive dilutions were used to obtain average generation times. We moved then towards methods requiring only one OD growth curve. 3.3.2. Estimation of growth rates with other methods on OD growth curvesAs s=ln and pn(x)=k×x in the OD range of 0.0–0.6, the model GP (Eq. (3)) became:
the model W (Eq. (5)) became:
and the model E (Eq. (6)) became:
The OD range corresponding to the exponential growth phase was estimated by plotting the observed growth rate dOD/dt against OD (Fig. 3). dOD/dt was estimated by ΔOD/Δt where ΔOD is the difference of OD between two successive readings and Δt is the time interval between two readings. We could consider that there were good linear correlations between dOD/dt and OD for OD-values below 0.1. We then applied Eq. (10) in the OD range of 0.0 –0.1.
In (9) and (10), x0 was the average of two initital VC and xmax was the average of the maximum bacterial concentrations estimated with VC fits. These three models gave growth rates strongly correlated with reference ones.
The average square roots of ratios (Fig. 4) for the models GP, W
and E were, respectively, 0.97±0.06 (SD), 1.05±0.08 (SD) and 0.97±0.04
(SD). The best model for which the spread was the smallest was then the model
E. As Dalgaard et al. (1994), we observed that this technique
under-estimated the maximum specific growth rates and that the growth rates
achieved must be divided by 0.972=0.94. This is not surprising
because OD is only detected when bacterial suspensions reached high cell
concentrations, and at those high concentrations, the specific growth rate is
significantly lower than
3.4. Lag timesLag times were estimated with models GP, W, E and M. With this last model, the maximum specific growth rates were the values obtained with model E divided by the calibration factor 0.94. Always for this model, lag times were the averages of two estimations obtained with two VC in exponential growth phase. Estimated lag times from OD growth curves with models GP, W and E were correlated with those obtained with VC data but numerous aberrant values were obtained with these models (Fig. 5). The average value of the ln of the ratio of lag times for the model GP was −0.56±0.35 (SD) indicating that lag times achieved with this method were almost the half of VC lag times. This phenomenon was observed by Hudson and Mott (1994) with Pseudomonas fragi who obtained an average ln of the ratio of −0.54±0.35 (SD). With models W and E, lag times achieved were on average almost the same that VC estimates but the ratios ranged from 0.1 to 2.7 (average ln of the ratios of, respectively, 0.16±0.56 and -0.10±0.29). This can be explained in part by the variability of the factor of proportionality between OD and VC. It was then not possible to obtain reliable estimations of lag times with only the OD growth curves. We then preferred to use some VC in the exponential growth phase to estimate the lag time knowing the maximum specific growth from the OD data (model M). This method gave ratios of lag times less spread (Fig. 5) and lag times achieved with this method were almost equal to lag times estimated with VC growth curves (average ln of the ratio of 0.05±0.12).
3.5. ValidationModels E (Eq. (10)) and M (Eq. (7)) were chosen to determine growth parameters with OD growth curves. To use model M, the initial cell count, the maximum bacterial concentration and viable counts in the exponential growth phase must be known. Maximum bacterial concentration was supposed to be equal to the average value achieved with the VC data regressions. Initial cell counts were estimated from two initial VC. Maximum specific growth rates were estimated by using model E on OD growth curves in the range of 0.0–0.1 and by dividing the growth rates obtained by 0.94. Lag times were the averages of two estimates obtained by using model M on two VC in exponential growth phase. Knowing all these parameters, predicted growth curves were constructed with the growth model f (Eq. (1)). Predicted growth curves were consistent with the observed VC for the eight experiments (Fig. 6).
4. CONCLUSION Turbidimetric measurements are an effective way to accurately estimate growth
rates of Listeria monocytogenes at different temperatures. It was shown
that a reliable estimation of
ACKNOWLIDGEMENTS This work was supported by the Association Vétérinaire d’Hygiène Alimentaire.
We would like to thank Agnès Brouillaud, Inès Giovannacci, Isabelle and André
Lebert and Jocelyne Rocourt who kindly provided the L. monocytogenes
strains, and Olivier Cerf, Cécile Lahellec and Pierre Pardon for their numerous
advice during the elaboration and the follow-up of this work.
REFERENCES Alber, S.A. and Schaffner, D.W., 1992. Evaluation of data transformations used with square root and Schoolfield models for predicting bacterial growth rate. Appl. Environ. Microbiol. 58, pp. 3337–3342. Arino, S., Andolfo, S., Leitz, R., Perrier, L. and Rosso, L., 1998. Détermination des températures cardinales de souches de levure par une méthode de gradient thermique dynamique. In: . Baranyi, J., Roberts, T.A. and McClure, P., 1993. A non-autonomous differential equation to model bacterial growth. Food Microbiol. 10, pp. 43–59. Barbosa, W.B., Cabedo, L., Wederquist, H.J., Sofos, J.N., Schmidt, G.R., 1994. Growth variation among species and strains of Listeria in culture broth. J. Food Prot. 57, 765–769, 775.. Bates, D.M. and Watts, D.G., 1988. . Nonlinear Regression Analysis and Its Applications Wiley, New York. Bégot, C., Desnier, I., Daudin, J.D., Labadie, J.C. and Lebert, A., 1996. Recommendations for calculating growth parameters by optical density measurements. J. Microbiol. Methods 25, pp. 225–232. Bégot, C., Lebert, I. and Lebert, A., 1997. Variability of the response of 66 Listeria monocytogenes and Listeria innocua strains to different growth conditions. Food Microbiol. 14, pp. 403–412. Box, G.E.P., Hunter, W.G. and Hunter, J.S., 1978. . Statistics For Experimenters Wiley, New York. Bratchell, N., Gibson, A.M., Truman, M., Kelly, T.M. and Roberts, T.A., 1989. Predicting microbial growth: the consequences of quantity of data. Int. J. Food Microbiol. 8, pp. 47–58. Bréand, S., Fardel, G., Flandrois, J.P., Rosso, L. and Tomassone, R., 1997. A model describing the relationship between lag time and mild temperature increase duration. Int. J. Food Microbiol. 38, pp. 157–167. Cheroutre-Vialette, M., Lebert, I., Hebraud, M., Labadie, J.C. and Lebert, A., 1998. Effects of pH or aw stress on growth of Listeria monocytogenes. Int. J. Food Microbiol. 42, pp. 71–77. Corman, A., Carret, G., Pavé, A., Flandrois, J.P. and Couix, C., 1986. Bacterial growth measurement using an automated system: mathematical modelling and analysis of growth kinetics. Ann. Inst. Pasteur/Microbiol. 137B, pp. 133–143. Dalgaard, P., Ross, T., Kamperman, L., Neumeyer, K. and McMeekin, T.A., 1994. Estimation of bacterial growth rates from turbidimetric and viable count data. Int. J. Food Microbiol. 23, pp. 391–404. Delignette-Muller, M.L., 1998. Relation between the generation time and the lag time of bacterial growth kinetics. Int. J. Food Microbiol. 43, pp. 97–104. Hudson, J.A. and Mott, S.J., 1994. Comparison of lag times obtained from optical density and viable count data for a strain of Pseudomonas fragi. J. Food Safety 14, pp. 329–339. Jones, C.E., Shama, G., Andrew, P.W., Roberts, I.S. and Jones, D., 1995. Comparative growth of Listeria monocytogenes in defined media and demonstration of growth in continuous culture. J. Appl. Bacteriol. 78, pp. 66–70. Kono, T., 1968. Kinetics of microbial cell growth. Biotech. Bioeng. 10, pp. 105–131. McClure, P.J., Cole, M.B., Davies, K.W. and Anderson, W.A., 1993. The use of automated rurbidimetric data for the construction of kinetic models. J. Ind. Microbiol. 12, pp. 277–285. Ratkowsky, D.A., Ross, T., McMeekin, T.A. and Olley, J., 1991. Comparison of Arrhenius-type and Belehradek-type models for prediction of bacterial growth in foods. J. Appl. Bacteriol. 71, pp. 452–459. Ratkowsky, D.A., Ross, T., Macario, N., Dommett, T.W. and Kamperman, L., 1996. Choosing probability distributions for modelling generation time variability. J. Appl. Bacteriol. 80, pp. 131–137. Rosso, L., Bajard, S., Flandrois, J.P., Lahellec, C., Fournaud, J. and Veit, P., 1996. Differential growth of Listeria monocytogenes at 4 and 8°C: consequences for the shelf life of chilled products. J. Food Prot. 59, pp. 944–949. Schaffner, D.W., 1994. Application of a statistical bootstrapping technique to calculate growth rate variance for modelling psychotrophic pathogen growth. Int. J. Food Microbiol. 24, pp. 309–314. Stephens, P.J., Joynson, J.A., Davies, K.W., Holbrook, R., Lappin-Scott, H.M. and Humphrey, T.J., 1997. The use of an automated growth analyser to measure recovery times of single heat-injured Salmonella cells. J. Appl. Microbiol. 83, pp. 445–455. Young, K.M. and Foegeding, P.M., 1993. Acetic, lactic and citric acids and pH inhibition of Listeria monocytogenes Scott A and the effect on intracellular pH. J. Appl. Bacteriol. 74, pp. 515–520. Zwietering, M.H., Jongenburger, I., Rombouts, F.M. and van’t Riet, K., 1990. Modelling of the bacterial growth curve. Appl. Environ. Microbiol. 56, pp. 1875–1881. Zwietering, M.H., Cuppers, H.G.A.M., de Wit, J.C. and van’t Riet, K., 1994. Evaluation of data transformations and validation of a model for the effect of temperature on bacterial growth. Appl. Environ. Microbiol. 60, pp. 195–203.
|
© 2005
Transgalactic Ltd (manufacturer of Bioscreen C software) |
Privacy Statement | P.O. Box
1393, 00101 Helsinki, Finland,
Last modified: May 25, 2005
| ||||||