Previous Article | Next Article ![]()
Applied and Environmental Microbiology, April 2006, p. 2721-2729, Vol. 72, No. 4
0099-2240/06/$08.00+0 doi:10.1128/AEM.72.4.2721-2729.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Wageningen University, Laboratory of Food Microbiology, P.O. Box 8129, 6700 EV Wageningen, The Netherlands
Received 22 June 2005/ Accepted 30 January 2006
|
|
|---|
|
|
|---|
In order to prevent illness, product contamination at manufacture and/or during preparation and growth after reconstitution must be minimized by appropriate control measures. Mathematical models are useful tools for evaluating the effectiveness of control measures. Depending on the source and the history of contaminating bacterial cells, which influence their physiological state, and the suitability of the product to sustain their growth, i.e., the product's (intrinsic) conditions and the environmental (extrinsic) conditions, microbial cells will either start to grow immediately or show a distinct phase of no apparent growth (the lag phase). In the case of E. sakazakii cells, reconstituted infant formulae offer rich growth environments that allow immediate proliferation provided that the cells are in a sound physiological state, that the external conditions (mainly temperature) are favorable, and that there is sufficient time for growth. Should lag times be apparent before growth, this may be a result of an injury to the E. sakazakii cells, from which they may gradually recover, as is evidenced by the start of cell proliferation (16). Baranyi and Roberts (1) emphasized that the lag time is a period of adjustment to a new environment, during which only intracellular conditions change. Growth models can simulate growth after reconstitution, and the effects of key intrinsic or extrinsic conditions can be determined. To develop growth models, insight into parameters describing growth of the microorganism, such as lag time and specific growth rate, is required.
This study describes the effects of a number of preculturing conditions on key growth parameters for E. sakazakii growing in reconstituted (with sterile water) powdered infant formula. Furthermore, the effects of temperature on specific growth rate and lag time were quantified and compared with literature values. Viable counts were used to construct growth curves that were used to derive the key growth parameters by curve fitting with the modified Gompertz equation as the primary growth model (19). During the secondary modeling step, the square root Ratkowsky model (12) and the secondary Rosso model (13) were fitted to the estimates of the specific growth rates at various temperatures. Likewise, the lag time data were fitted with the logarithm of the inverse of the Ratkowsky model and the hyperbolic equation (20). The resulting parameters permit useful predictions of the growth of E. sakazakii in reconstituted infant formula and aid in the design of effective control measures to reduce exposure of susceptible consumers in both hospital and household settings.
|
|
|---|
Preparation of the bacterial suspension.
Strain ATCC 29544 was cultured by transferring 250 µl of the stock culture into 250 ml of BHI, followed by incubation at 37°C. Cells were incubated for 3, 6, 16, 24, and 72 h to obtain, respectively, exponential-phase, early-stationary-phase, mid-stationary-phase, stationary-phase, and late-stationary-phase cells. Strains MC10, MM9, and S94 were incubated at 37°C for 18 h to obtain mid-stationary-phase cells. In order to obtain enough cells in the lag phase, 1 ml of a stock culture of ATCC 29544 that had been maintained at 20°C was diluted in 2 ml of glycerol, whereupon the suspension was centrifuged (Hermle Z 231 M; B. Hermle GmbH & Co., Gosheim, Germany) for 10 min at 10,000 x g. Harvested cells were transferred into 30 ml of BHI and incubated for 1.5 h at 37°C.
In order to obtain cells for spiking of the dry infant formula, these BHI-grown cultures were centrifuged for 5 min at 20°C at 2,958 x g (Mistral 3000i; MSE, Leicester, United Kingdom). Cells were washed twice in 1% physiological salt solution and subsequently suspended in 30 ml of 1% physiological salt solution, for the lag-phase cells, or 250 ml of 1% physiological salt solution, resulting in a cell suspension of about 104 CFU/ml for lag-phase cells and between 105 and 107 CFU/ml for the other growth phases.
Spiking of the infant formula.
Infant formula was bought in local shops; the numbers of viable bacteria in the infant formula were sufficiently low to prevent them from influencing the growth of the spiked cells (data not shown). The obtained bacterial suspension was sprayed on commercial dry infant formula (1:50, wt/wt) with a perfume sprayer (designed by Gérard Brinard, DA Drogisterij, Leusden, The Netherlands). The final bacterial concentration at 3 to 4 days after spiking was 102 to 104 CFU/g of dry powder, and the infant formula maintained a water activity of <0.22. All growth experiments were performed within 10 days after spiking the infant formula.
Design.
The initial estimates of the lowest temperature supporting growth (Tmin), the highest temperature supporting growth (Tmax), and the specific growth rate (µopt) at the optimal growth temperature (Topt), shown in Table 1, were used to predict the specific growth rates at various temperatures. The lag time at each temperature was initially estimated by the k value, which is the product of the specific growth rate (µm) and lag time (
) and which is known to have a large variability but can be considered constant over a wide range of temperatures (18).
|
View this table: [in a new window] |
TABLE 1. Initial estimates of parameters for growth of E. sakazakii in reconstituted infant formula, based on published data
|
Growth experiments.
To prepare samples for growth experiments, 10-g portions of spiked infant formula were reconstituted in 100 ml of sterilized tap water. In experiments to determine the effects of various growth phases, bottles with strain ATCC 29544 were incubated as follows (the number of bottles is given in parentheses): 10°C (n = 9), 21°C (n = 11), 29°C (n = 14), 37°C (n = 14). Mid-stationary-phase cells of strain ATCC 29544 were additionally used to assess (in duplicate) growth in reconstituted infant formula at the following temperatures: 8, 14, 38, 39, 41, 43, 45, 46, 47, 48, 49, and 50°C. Mid-stationary-phase cells of strains MM9, MC10, and S94 were incubated only at 29 and 37°C.
Growth of E. sakazakii was measured at various time intervals, depending on the temperature of incubation; appropriate dilutions were made in peptone saline solution (NaCl [8.5 g/liter] and neutralized Bacteriological Peptone [1g/liter]; Oxoid, Basingstoke, England). Samples were surface plated onto tryptone soy agar (Oxoid, Basingstoke, England) with a spiral plater (Eddy Jet; IUL Instruments, I.K.S., Leerdam, The Netherlands). Inoculated plates were incubated for 20 to 24 h at 37°C before manual counting.
Data analysis.
For describing the evolution of the microbial count with time, a primary model was used. The three kinetic parameters, namely, lag time, specific growth rate, and maximum population density, were estimated by fitting with the modified Gompertz equation (equation 1):
![]() | (1) |
[h]), the specific growth rate (µm [h1]), and the dimensionless asymptotic value (A) (19) at the tested temperature for each growth curve. In order to obtain reliable estimates for the growth parameters, experimental growth curves had to meet the following requirements: (i) at least two data points should fall within the lag time, unless the lag time was shorter than 2 h; (ii) for data points within the exponential phase, there should be at least three data points over a range of 3 h and 3 logs; and (iii) three data points should fall within the stationary phase, at least 1 h apart. Growth curves that did not meet these requirements were excluded.
A Bélehrádek-type model (equation 2), also known as the (expanded) square root model of Ratkowsky (12), was used to describe the relation between the specific growth rate and the temperature. This model contains four parameters, of which two are easily interpretable, Tmin and Tmax. If Tmin < T < Tmax, then
![]() | (2) |
Tmin or T
Tmax, then µm = 0, where Tmin is the extrapolated minimum temperature (°C) at which the specific growth rate (µm [h1]) is zero, Tmax is the extrapolated maximum temperature at which µm is zero, and b (°C1 h0.5) and c (°C1) are the so-called Ratkowsky parameters (12).
The secondary growth model of Rosso et al. (13) (equation 3) was used as well to describe the effect of temperature on growth rate. This model contains all four interpretable parameters: µopt, Tmin, Tmax, and Topt. If Tmin < T < Tmax, then
![]() | 3 |
Tmin or T
Tmax, then µm = 0, where Tmin and Tmax are defined as in equation 2, Topt is the temperature (°C) at which the specific growth rate µm (h1) is optimal, and µopt is the µm at the optimal temperature.
The logarithm of the inverse of the secondary Ratkowsky model (equation 4) and the hyperbolic equation (equation 5) were used (20) to describe the lag time/temperature relation.
![]() | (4) |
![]() | (5) |
Statistical analysis.
In order to determine whether preculturing conditions have a significant effect on lag times and/or specific growth rates, the data obtained with strain ATCC 29544 at 10, 21, 29, and 37°C were subjected to univariate analysis of variance. Lag time data were log transformed and specific growth rate data were square root transformed, in order to obtain homogeneity of variance. A significance level of 5% was used. All data analyses were performed with SPSS (SPSS, release 11.5, for Microsoft Windows 95/98/NT/2000; SPSS Inc., Chicago, Ill.).
Fitting was done by minimizing the residual sum of squares (RSS) with both the solver function in Excel and Table Curve 2D for Windows, version 2.03, for verification.
Standard deviations were calculated with Excel and are reported as plus or minus the means.
|
|
|---|
![]() View larger version (10K): [in a new window] |
FIG. 1. (A) Predicted and measured growth curves of E. sakazakii ATCC 29544 at 21°C, precultured to mid-stationary phase. The different symbols indicate different replicate experiments. The solid line is the design growth curve at 21°C, calculated with the cardinal values shown in Table 1 under "initial estimate"; dotted lines represent fits of the modified Gompertz equation to each single experiment. (B) Predicted and measured growth curves of E. sakazakii ATCC 29544 at 37°C, precultured to mid-stationary phase. The different symbols indicate different replicate experiments. The solid line is the design growth curve at 37°C, calculated with the cardinal values shown in Table 1; dotted lines represent fits of the modified Gompertz equation to each single experiment.
|
.
) (Fig. 2A and B). With increasing temperature, µm strongly increased. Values for specific growth rates varied from 0.12 ± 0.04 h1 at 10°C to 2.29 ± 0.45 h1 at 37°C, and no apparent effect of the physiological state was found. In Fig. 2B, it is shown that the lag time decreased with increasing temperature and that there was also no apparent effect of the various physiological growth phases on the lag time.
![]() View larger version (8K): [in a new window] |
FIG. 2. (A) Square roots of specific growth rate data for various physiological growth phases of strain ATCC 29544 as a function of temperature. x, lag-phase cells; , exponential-phase cells; , early-stationary-phase cells; x , stationary-phase cells; , late-stationary-phase cells. (B) Log lag times for various physiological growth states of ATCC 29544 as a function of temperature. x, lag-phase cells; , exponential-phase cells; , early-stationary-phase cells; x , stationary-phase cells; , late-stationary-phase cells.
|
· µm) was also not significantly influenced by cell history (Table 2). |
View this table: [in a new window] |
TABLE 2. Statistical evaluation of univariate analyses of variance for the effects of the physiological growth phase of E. sakazakii strain ATCC 29544 cells on lag time, specific growth rate, and the product of both (k) at 10, 21, 29, and 37°C
|
.
![]() View larger version (10K): [in a new window] |
FIG. 3. Square roots of measured and fitted specific growth rates as a function of temperature. Shown are growth rates of ATCC 29544 precultures to various growth phases, as estimated by the fit of the modified Gompertz model to each individual growth curve , and growth rates for mid-stationary-phase cells of MM9 , S94 , and MC10 . Other symbols: +, growth rates published by Nazarowec-White and Farber (9); : growth rates published by Iversen et al. (5); dashed line, fit by the secondary growth model of Ratkowsky; dotted line, fit by the secondary growth model of Rosso (fitted to square root transformed data of the present study only).
|
![]() View larger version (9K): [in a new window] |
FIG. 4. Log lag times as a function of temperature. , ATCC 29544 precultures to various growth phases. Also shown are data for mid-stationary-phase cells of MM9 , S94 , and MC10 ; lag times at 10 and 23°C as published by Nazarowec-White and Farber (9) (+); and logarithmic transformed lag time data from the present study, modeled with the hyperbolic model (dashed line) and the reciprocal Ratkowsky model (solid line).
|
All of the estimated specific growth rates were combined and modeled with equations 2 and 3, as shown in Fig. 3, where the square root of the specific growth rate (
µm) is plotted as a function of temperature. The secondary growth parameters derived from optimal fits are shown in Tables 3 and 4. From Fig. 3 it is apparent that the differences between the fits of these models are smaller than the experimental variability. The RSS for the secondary Rosso model (equation 3) was 1.31, and the RSS for the square root Ratkowsky model (equation 2) was 1.28. Although its RSS was slightly higher, the secondary Rosso model was used for further evaluation because it consists of four parameters that all have biological meaning and can be interpreted as such. Transforming the square root µm data from Fig. 3 back to µm, it can be seen that the specific growth rate of E. sakazakii varied from 0.115 h1 (
µm = 0.339 h0.5) at 10°C to 1.113 h1 (
µm = 1.063 h0.5) at 46°C and 2.242 h1 (
µm = 1.498 h0.5) at 37°C. For comparison, the experimental specific growth data measured by Iversen et al. (5) and Nazarowec-White and Farber (9) have been transformed and included in Fig. 3.
|
View this table: [in a new window] |
TABLE 3. Parameter values for the effects of temperature on the specific growth rate (µm) of E. sakazakii in reconstituted infant formula resulting from fits by the Ratkowsky secondary growth model (equation 2)
|
|
View this table: [in a new window] |
TABLE 4. Parameter values for the effects of temperature on the specific growth rate (µm) of E. sakazakii in reconstituted infant formula resulting from fits by the Rosso secondary growth model (equation 3)
|
The effect of temperature on lag time was described by fitting the reciprocal square root relation (equation 4) and the hyperbolic model (equation 5) to the logarithmic transformation of the data. It is assumed that the Tmin and Tmax values are equal to the Tmin and Tmax values of equation 3, describing the specific growth rate. Parameters estimated with both models are given in Table 5. The logarithmic transformed lag times, as fitted by the reciprocal square root relation and the hyperbolic model, are represented in Fig. 4. From this graph it can be concluded that both models fit the data reasonably well. The RSS for the hyperbolic model was 2.11, and the RSS for the reciprocal square root relation was 2.29. Figure 5 shows the dimensionless parameter k, the product of the lag time and the specific growth rate (µm ·
), as a function of temperature. The k value between 8 and 47°C was 5.08 ± 3.37. At temperatures ranging from 20 to 46°C, values for k were between 0.82 and 11.6, with an average of 4.05 ± 1.92. Below 20°C, it increased to an average of 10.06 ± 4.49. Use of a k value of 5.08 to predict the lag time at any temperature and equation 3 to determine the specific growth rate resulted in an RSS of 16.7. Based on the RSS values for the models describing the lag time, it can be concluded that the hyperbolic model and reciprocal square root model best convey the experimental lag times. The reciprocal square root model is recommended, since it has the ability to increase the lag time at higher temperatures and it contains more interpretable parameters.
|
View this table: [in a new window] |
TABLE 5. Parameter values for the effects of temperature on the lag time ( ) of E. sakazakii in reconstituted infant formula
|
![]() View larger version (7K): [in a new window] |
FIG. 5. Parameter k (the product of and µm) as a function of temperature for each growth experiment with strain ATCC 29544 precultured to various growth phases and grown at temperatures from 8 to 47°C. x, lag-phase cells; , exponential-phase cells; , early-stationary-phase cells; x , stationary-phase cells; , late-stationary-phase cells. The dotted line represents the average value for the data from 20 to 46°C.
|
![]() View larger version (8K): [in a new window] |
FIG. 6. Initial inoculum of strain ATCC 29544 precultured to various physiological states and maximum population densities at various temperatures. , maximum population density; , initial number of cells. The dotted line represents the average value for the maximum population density data from 8 to 46°C.
|
|
|
|---|
Our specific growth rate data were compared to specific growth rates published by Nazarowec-White and Farber (9) and Iversen et al. (5), as shown in Fig. 3. In our study, the specific growth rates of E. sakazakii cells spiked into dry infant formula were comparable to those of Nazarowec-White and Farber (9), who reported similar specific growth rates. The specific growth rates reported by Iversen et al. (5) were consistently over 10% lower than ours. One cause for the difference with the latter study may be differences in growth conditions: in our experiments, spiked dry powdered infant formula samples reconstituted with sterile tap water were used as growth medium, whereas Iversen et al. measured growth in other media by inoculating them with an overnight tryptic soy broth culture. Another difference is the use by Iversen et al. of the rapid automated bacterial impedance technique, which measures growth only at much higher levels (>106 CFU/ml) than the plate count method used in our study (15). In the present study, the plate count technique was used and growth was measured in a more representative range and over multiple logs of bacterial counts; thus, the estimates for lag time and specific growth rate can be considered more accurate (2).
The experimental design, based on initial estimates of specific growth rate and lag time, allowed prediction of growth with enough accuracy to determine the sampling times and dilutions appropriate for measuring growth curves, such that our quality requirements were met for 95% of the experimental growth curves. The modified Gompertz model (equation 1), the secondary Rosso model (equation 3), and the reciprocal square root model (equation 4) successfully estimated the lag times and specific growth rates for the whole growth temperature range (Tables 4 and 5). The theoretical minimal and maximal growth temperatures, as fitted with both secondary growth models, were 2.19°C and 48.9°C (equation 2) and 3.60°C and 47.6°C (equation 3), respectively. In practice, however, growth of E. sakazakii strains has been observed between 5.5°C and 47°C (3).
The fact that small numbers of E. sakazakii cells have occurred in dry infant formulae (4, 8, 9), in combination with the relatively short lag times and high specific growth rates found in this study, underscores the need for careful preparation and use of dry infant formulae. A study by Pagotto et al., however, showed that 105 CFU/ml of certain E. sakazakii strains could be lethal to suckling mice after ingestion, though not all strains appeared to be pathogenic (10). As there is a lack of information about virulence factors, not all E. sakazakii strains necessarily need to be regarded as potential pathogens. Nevertheless, the results described in this paper (Tables 4 and 5) permit predictions of the growth of E. sakazakii in reconstituted infant formulae under conditions that closely mimic conditions in the actual food product in both hospital and household settings, and they will be useful in designing effective control measures.
Financial support by Nestlé Research Center, Lausanne, Switzerland, is gratefully acknowledged.
|
|
|---|
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»