Abstract

The selection of a primary model to describe microbial growth in predictive food microbiology often appears to be subjective. The objective of this research was to check the performance of different mathematical models in predicting growth parameters, both by absorbance and plate count methods. For this purpose, growth curves of three different microorganisms (Bacillus cereus, Listeria monocytogenes, and Escherichia coli) grown under the same conditions, but with different initial concentrations each, were analysed. When measuring the microbial growth of each microorganism by optical density, almost all models provided quite high goodness of fit () for all growth curves. The growth rate remained approximately constant for all growth curves of each microorganism, when considering one growth model, but differences were found among models. Three-phase linear model provided the lowest variation for growth rate values for all three microorganisms. Baranyi model gave a variation marginally higher, despite a much better overall fitting. When measuring the microbial growth by plate count, similar results were obtained. These results provide insight into predictive microbiology and will help food microbiologists and researchers to choose the proper primary growth predictive model.

1. Introduction

Predictive microbiology enables, through the use of mathematical models, estimating the behaviour of microorganisms under certain circumstances [1], based upon the premise that the responses of microorganisms to environmental factors are reproducible. The ability to predict both the growth of microorganisms, as affected by different environmental factors, and the survival of microorganisms as a result of preservative treatments is an important tool for evaluating the safety and shelf life of food products.

Before predictive microbiology can be applied to the food industry, mathematical models that adequately describe microbial behaviour are needed. There are a number of sigmoid equations and several models that have been used as growth functions. They all differ in “ease of use” and number of parameters in the equation. Some authors have compared the behaviour of different growth models, from different viewpoints, including mathematical measures of goodness of fit [2] and/or other statistical criteria [35]. The usual measures of goodness of fit for model comparison in previous studies were done by calculating the bias () and accuracy () indices as provided by Ross [6], the coefficient of determination (), the residual mean square error (RMSE), or the -test. Other authors [7, 8] have focussed on direct comparisons of particular growth parameters as predicted by various models.

These studies have reached different conclusions. Hence, there is significant disagreement in literature on which is the best-fitting model for predictive microbiology. The selection of a model in predictive food microbiology often appears to be subjective. Based on reports in the literature, Gompertz, Baranyi, Richards, logistic, and three-phase linear models are the most widely used [5, 911].

The growth curve has been mostly expressed in terms of microbial numbers (concentration of colony forming units), but also in terms of optical density as an indirect measurement. The measurement of absorbance is a rapid, nondestructive, inexpensive, and relatively easy-to-automate method to monitor bacterial growth, as compared to many other techniques and particularly when compared to classical viable count methods. When modelling optical density growth curves, the fitted parameters are different from the population growth parameters derived from viable counts. The rate of increase of the optical density does not express the maximum specific growth rate and the detection time is not equivalent to the lag time, unless the initial inoculum is greater than the detection limit. In spite of the limitations related to detection thresholds, correlation with the parameters derived from viable count growth curves, and inability to model growth in turbid liquid foods and in solid food matrixes, numerous techniques and mathematical growth models have been used in recent years for estimation of growth rates and lag times from absorbance data [5, 7, 8, 1214]. In the opinion of Dalgaard and Koutsoumanis [7], absorbance techniques should be limited to conditions where high cell densities are reached, such as those resembling the growth of spoilage bacteria in foods. Even assuming the limitations of absorbance to build growth curves, it may be useful, if not to obtain very precise growth kinetic parameters, at least to compare the growth of different cultures or of the same cultures but in different conditions.

Work modelling the behaviour of bacteria in foods has shown that the lag phase is more difficult to predict than is the specific growth rate [15], mainly because of the influence in lag time of the physiological state of individual bacterial cells and, to a minor extent, of the inoculum size. The physiological state of the cells is affected by their previous growth environment and by exposure to stress conditions, which can extend the lag time considerably and also increase individual cell lag time variability [1618]. However, microorganisms with a similar precultural history exposed to the same favourable growth conditions should be in a similar optimum physiological state and, thus, its effect on lag time variability is negligible [19]. Regarding the inoculum size, Baranyi and Pin [20] showed that as the cell number in the inoculum decreases, the population lag increases by an amount that depends on the distribution of individual lag times and the maximum specific growth rate. Augustin et al. [21] showed that the inoculum level effect can be explained by an increasing variability in individual cell lag time when stress factors become more stringent, and Baranyi and Pin [20] found that, under optimum growth conditions, this effect would only be expected at inoculum levels below about 102-103 cells, because the impact of variability among a small population of cells can become more important on lag time [2, 22].

Hence, a sufficiently large population of microorganisms exposed to exactly the same favourable growth conditions and with a similar precultural history behave in a similar way; that is, they should show the same growth parameters, growth rate, and lag phase duration.

The objective of this research was to check the performance of different mathematical models in predicting growth parameters, by both absorbance and plate count methods. For this purpose, growth curves of three different microorganisms (Bacillus cereus, Listeria monocytogenes, and Escherichia coli) grown, each species, under the same conditions, but with different initial concentrations, were analysed.

2. Materials and Methods

2.1. Microorganisms

Bacillus cereus, Listeria monocytogenes, and Escherichia coli were chosen as representative microorganisms for spore-forming, Gram-positive, and Gram-negative bacteria, respectively.

B. cereus INRA-AVTZ415 was kindly provided by the Institut National de la Recherche Agronomique (INRA, Avignon, France). L. monocytogenes and E. coli type strains (CECT 4031 and CECT 515, resp.) were provided by the Spanish Type Culture Collection (CECT).

To inoculate the growth media, B. cereus vegetative cells were grown at 30°C in brain heart infusion broth (BHI; Scharlau Chemie S.A., Barcelona, Spain), until the stationary phase of growth was reached. L. monocytogenes vegetative cells were grown at 37°C in tryptic soy broth (TSB; Scharlau Chemie) supplemented (w/v) with 0.6% yeast extract (YE; Scharlau Chemie), until the stationary phase of growth was reached. E. coli vegetative cells were grown at 37°C in TSB+YE acidified to pH 5 with ClH (Panreac Química, Barcelona, Spain) until the stationary phase of growth was reached. These growth conditions were chosen as favorable for these microorganisms.

2.2. Optical Density Growth Curves

100-well microtitre plates were filled with 400 μL of the growth media (BHI for B. cereus, TSB+YE for L. monocytogenes, and pH 5 TSB+YE for E. coli) and were inoculated with the microorganisms and incubated in a Bioscreen C analyzer (Oy Growth Curves Ab Ltd., Helsinki, Finland). Initial concentrations in the growth media were 101, 102, 103, 104, 105, and 106 CFU mL−1 for B. cereus, 102, 103, 104, 105, 106, and 107 CFU mL−1 for L. monocytogenes, and 102, 104, and 106 CFU mL−1 for E. coli. In order to avoid variability derived from differences in the physiological state of different cultures, all the growth curves from each bacterium were obtained from a single bacterial culture. Growth media were incubated at 30°C for B. cereus and at 37°C for L. monocytogenes and E. coli. At 20 min intervals, the optical density (OD) of the samples using a wideband filter (420–580 nm) was measured.

For each combination of microorganism and initial concentration, 25–30 repetitions were performed, except for B. cereus 101 and 102 CFU mL−1 initial concentrations, where 13 and 21 repetitions were performed, and for L. monocytogenes 103, 105, and 107 CFU mL−1 initial concentrations, where 5 repetitions each were performed. Growth curves were obtained by plotting the OD against the exposure time. A total of 345 individual growth curves were generated by absorbance measurements.

2.3. Plate Count Growth Curves

50 mL flasks of the growth media were inoculated with the microorganisms and incubated with agitation at 500 rpm. Growth media and incubation temperatures were the same as those used for optical density growth curves. Initial concentrations in the growth media were 101, 103, and 105 CFU mL−1 for B. cereus and 102, 104, and 106 CFU mL−1 for L. monocytogenes and E. coli. At preset time intervals, samples were taken, properly diluted in buffered peptone water (BPW, Scharlau Chemie), and incubated in BHI agar (BHIA, Scharlau Chemie) for 24 h at 30°C for B. cereus and in tryptic soy agar (TSA, Scharlau Chemie) + YE for 24 h at 37°C for L. monocytogenes and E. coli. Growth curves were performed in duplicate.

Growth curves were obtained by plotting log CFU mL−1 against the exposure time.

2.4. Mathematical Models

Analyses of the growth curves were performed using five primary growth models. These growth models were based on either linear (derived from the Monod model) or nonlinear (Gompertz, Logistic, Richards, and Baranyi) equations (Table 1) and reparameterized to reflect microbial growth parameters as derived by Zwietering et al. [3].

Curve fitting of three-phase linear, Gompertz, logistic, and Richards models was done using the curve-fitting tool of Matlab 7.0 (Math Works, Natick, USA) with which 95% confidence limits (CL) for growth parameters and , RMSE, and sum of square error (SSE) of fit were calculated.

The curve fitting of Baranyi’s equation was done using DMFit 2.0 program and the model of Baranyi and Roberts [23] as kindly provided by Dr. József Baranyi. This program provided standard error for each growth parameter and and standard error of fit. With these values, 95% CL of growth parameters and RMSE of fit were calculated.

Analysis of variance, medians, and quartiles for box and whisker plots were calculated using StatGraphics (StatPoint Technologies, Warrenton, USA). values were always lower than 0.05.

3. Results and Discussion

3.1. Optical Density Growth Curves

Figure 1 shows the optical density growth curves plotted with the average OD values at each sampling time of B. cereus INRA-AVTZ 415 at 30°C in BHI (a), L. monocytogenes CECT 4031 at 37°C in TSB+YE (b), and E. coli CECT 515 at 37°C in pH 5 TSB+YE (c) starting at 102, 104, and 106 CFU mL−1. These average growth curves correspond to 25–30 individual growth curves each. The slopes of all the growth curves corresponding to a microorganism were parallel in the exponential growth phase; that is, the growth rates were similar, as they should correspond to different cultures of the same microorganism growing exactly in the same conditions. However, a progressive decrease of lag phase duration was observed as the initial concentration increased. This decrease in lag time can be easily explained because the culture spends less time in reaching a concentration where the absorbance increases. Actually, several authors have recently used simple equations to derive real lag time and provided the initial concentration of microorganisms and the observed lag time from OD measurements [2426]. The relatively low standard deviation values obtained, especially for growth curves starting at 104 and 106, are indicative of the repetitivity of the growth curves. It can be noticed that error bars are bigger as the initial concentration decreases, as a consequence of the increased variability among a small population of cells [2].

Table 2 shows the growth parameters given by three-phase linear, Gompertz, logistic, Richards, and Baranyi models for the growth curves of B. cereus plotted with the average OD values shown in Figure 1(a). All the models tested provided the values that could be expected for growth parameters of the three growth curves selected of this microorganism. Figure 1(a) shows that all three average growth curves of B. cereus were approximately parallel in their exponential growth phases. Although the curve starting at 102 CFU mL−1 had a slightly lower slope, only logistic and Richards models showed significant differences between this and the other two curves of B. cereus starting at higher initial concentrations (Table 2). When comparing the growth rate values given by the different models, three-phase linear model gave the lowest values (0.27–0.30 OD units h−1), followed by Baranyi model (0.28–0.31 OD units h−1), logistic and Gompertz (0.31–0.35 OD units h−1), and Richards model (0.30–0.42 OD units h−1). Also, all models provided shorter lag phases at higher initial concentrations, as shown in Figure 1(a). Initial absorbance () had values between 0.071 and 0.180 for all growth curves and models. The increase in absorbance from initial to final optical density () had values between 0.975 and 1.107 for all growth curves and models. Hence all models seemed to perform adequately, providing expected growth parameter values for these growth curves. values were higher than 0.995 for all growth curves and all models, except for three-phase linear model, which gave values as low as 0.989, hence being the model which provided worst fit to the data, as also it could be expected from a model consisting of three straight lines. All models gave RMSE values lower than 0.035.

Table 3 shows the growth parameters given by the five models for the growth curves of L. monocytogenes shown in Figure 1(b). Similar results were obtained, although in this case the only model which did not provide significant differences for the three growth rate values for the three initial concentrations tested was the three-phase linear. Again, the same order was obtained, three-phase linear model giving the lowest values (0.15–0.17 OD units h−1), followed by Baranyi model (0.15–0.20 OD units h−1), logistic and Gompertz (0.17–0.22 OD units h−1), and Richards model giving the highest values (0.18–0.32 OD units h−1). Also, all models provided shorter lag phases at higher initial concentrations. values were between 0.119 and 0.149 and values between 0.712 and 0.784. Again, values were higher than 0.990 for all growth curves and all models, except for three-phase linear model, which gave values as low as 0.930. Three-phase linear model also gave RMSE values as high as 0.061.

Table 4 shows the growth parameters given by the five models under study for the growth curves of E. coli shown in Figure 1(c). In this case, some unexpected results were obtained. Three-phase linear, Gompertz, logistic, and Baranyi models behaved as previously described for B. cereus and L. monocytogenes. Each model provided very similar values for the growth rates of the three growth curves, with only slight differences among them. However, only logistic and Baranyi models did not provide significant differences. For this microorganism, again, the three-phase linear model gave the lowest value for the growth rate (0.10–0.11 OD units h−1), but in this case together with Baranyi model, and again logistic and Gompertz models gave higher values (0.12–0.13 OD units h−1). However, Richards model gave unexpectedly low values for growth rates (<0.001 OD units h−1). All models provided shorter lag phases at higher initial concentrations. values were between 0.060 and 0.176 for all growth curves and models, except for the logistic model, which gave an exceptional low value of 0.000 (fixed at bound) for growth curve starting at 106 CFU mL−1. Logistic model gave values of SSE higher than 0.1 for all three growth curves of E. coli, which were almost twice the values obtained for any other growth curve shown in Table 2, 3, or 4. In this case all values were higher than 0.99 for all growth curves and all models, including the three-phase linear model.

After the results shown in Figure 1 and Tables 2, 3, and 4, the growth models selected for the rest of the investigation were three-phase linear, Gompertz, and Baranyi. The three selected models provided expected values for the growth parameters and fitted data properly, as shown by the statistics analysed, even in the case of the three-phase linear model, which was the one to provide worse goodness of fit, that is, lowest and highest RMSE values. Richards and logistic models were disregarded because they were not able to fit properly all these typical growth curves; that is, they gave abnormal values for growth parameters in some occasions.

Each model had a trend in providing higher or lower values, three-phase linear model giving consistently the lowest values, followed by Baranyi and Gompertz models, in this order, for both the growth rate and the lag phase. In a previous comparison of these three models, Buchanan et al. [2] already noticed and explained this effect on the basis of the nature of each model. These authors also highlighted the correlation existing between lag phase duration and specific growth rate values. This gives an explanation for the differences found in the values provided for growth rate and lag phase among these predictive models.

Hence, it seems that, depending on the predictive model chosen, values for growth rates and lag times will be consistently higher or lower. Then the question is which growth model can be considered the best for describing the true population growth and why. This can be answered with a deep analysis of the analytical aspects of the models, mechanistic elements, number of parameters, fitting properties, and so forth. However, such analyses have already been performed in the past [3], but researchers continue to use and compare different growth models [5, 12], probably because concluding results have not been reached yet.

In this context, our viewpoint is that the best performing model is the one which, for the growth parameters of different cultures of one microbial strain growing under exactly the same conditions and with the same precultural history, gives closer values.

With this purpose, an extensive analysis of growth curves was performed with the three selected growth models. Only growth rate was considered at this stage, since this parameter should be similar for all growth curves of the same microbial strain, even when starting from different initial inocula levels. A total of 345 individual growth curves were analysed, including those individual growth curves used to build the average OD growth curves shown in Figure 1.

Table 5 shows, as an example, growth rate values obtained for all E. coli individual growth curves with three-phase linear, Gompertz, and Baranyi growth models. With these data and those for the individual growth curves of B. cereus and L. monocytogenes, extensive statistical analyses were performed, including analyses of variance, medians, and quartiles. Analyses of variance showed that initial concentration did not influence growth rate and that significant differences were found among growth rate values given by the different models, as already pointed out. The statistical analyses also showed that, within models, growth rate data were not normally distributed, and significant differences were also found among median growth rate values.

The last step of this research was to analyse variation of growth rate values. Table 6 shows average and standard deviation for growth rate of all three microorganisms as obtained by three-phase linear, Gompertz, and Baranyi models, and Figure 2 shows box and whiskers plots for the growth rate of the three microorganisms. These plots are based on the median and withstand perturbations caused by outliers better than plots based on the average. Since growth rate data were not normally distributed, box and whiskers plots are more appropriate than average and standard deviation shown in Table 6.

The results shown in Figure 2 clearly show that Gompertz model has a higher degree of variation in the growth rate values than the three-phase linear and Baranyi models. Gompertz model also generated some outlier growth rate values which extended the “whiskers” several units in the case of B. cereus and L. monocytogenes. It is certainly possible to obtain correct values for these outliers by changing the initial values when performing the nonlinear regression, but the purpose of this research was not to optimize the models, but to test them in order to choose the model that performs best. Hence initial values were not changed. Three-phase linear model provided slightly less variation than Baranyi, although the goodness of fit of this model is considerably worse than that of Baranyi model.

3.2. Plate Count Growth Curves

In order to double check the results obtained, several plate count growth curves of these same microorganisms were modelled with the five growth models. Figure 3 shows the growth curves plotted with the average log CFU mL−1 at each sampling time of B. cereus INRA-AVTZ 415 at 30°C in BHI (a), L. monocytogenes CECT 4031 at 37°C in TSB+YE (b), and E. coli CECT 515 at 37°C in pH 5 TSB+YE (c) starting at different initial concentration of microorganisms. As usual for plate count growth curves, these curves included less data per curve and did not show data points in all growth phases. In this way, models were forced into the common situation of scarce data points. Tables 7, 8, and 9 show the growth parameters given by three-phase linear, Gompertz, logistic, Richards, and Baranyi models for the growth curves of B. cereus, L. monocytogenes, and E. coli shown in Figures 3(a), 3(b), and 3(c), respectively.

Similar results to those obtained with absorbance data were obtained when modelling the data from plate count growth curves. Three-phase linear model gave the lowest values for growth rates and lag times, followed by Baranyi, Gompertz, and logistic, in that order, for all three microorganisms. For these plate count growth curves, Richards model was not able to provide satisfactory values for growth rate for any of the curves. It is worth noting that some of the growth curves did not reach the stationary growth phase, and both three-phase linear and Baranyi models did not provide values for maximum population levels in these cases, while all other models provided uncertain estimations for this parameter. Again, three-phase linear was the model which had the worst goodness of fit, with values as low as 0.92 and RMSE values as high as 0.49.

Table 10 shows average and standard deviation for growth rate of all three microorganisms as obtained by three-phase linear, Gompertz, logistic, and Baranyi models with plate count growth curves. Again, when comparing the similarities in the growth rate values (Table 10), three-phase linear was the model to give less variation (lower standard deviation) for the different growth curves.

Comparisons of the behaviour of different growth models reported in literature have reached different conclusions. Zwietering et al. [3] studied the growth of Lactobacillus plantarum in MRS medium at different temperatures and concluded that the Gompertz model was the best-fitting model. When these authors [3] extended the study to several microorganisms they reached similar conclusions. Buchanan et al. [2] reported that the three-phase linear model was more robust than the Gompertz and Baranyi models in terms of successfully fitting growth curve data. In their research [2] they fitted experimental data for E. coli O157:H7. Schepers et al. [27] found that the Richards model was the best growth model for Lactobacillus helveticus grown at different pH values, and Dalgaard and Koutsoumanis [7] also agreed that Richards model gave the best estimates for absorbance growth curves obtained with mixtures of different microbial strains isolated from spoiled seafood and incubated in different conditions to obtain a wide range of growth yields. López et al. [5] concluded, after a detailed statistical evaluation, that the Baranyi model showed the best behaviour for the growth curves studied. Baranyi and three-phase linear models showed the best fit for plate count data of Yersinia enterocolitica grown under different conditions of pH, temperature, and CO2. Richards model was the best-fitting optical density data of different bacterial and fungal species grown under different conditions. However, Mytilinaios et al. [12] found more recently that Baranyi was the most capable model to fit optical density data obtained for L. monocytogenes at different temperatures and pH values. The Weibull model also adequately described microbial growth [5]. Baty and Delignette-Muller [4] found that Baranyi was the best curve-fitting model for most curves, although they noted that the intermodel variability was frequently minor comparing to the imprecision of the parameter estimates, due to the low quantity and quality of the data used to build the growth curves. Actually, the low quality data used by these authors correspond to the datasets tabulated in the work by Buchanan et al. [2]. Pal et al. [28] also showed that Baranyi model provided the best fit for a majority of growth curves obtained for L. monocytogenes at low temperatures in liquid cultures, although there was no significant difference among all the primary growth models analysed. Maybe the use of so many microorganisms growing under such different growth conditions could explain, at least in part, the differences in the conclusions reached by authors regarding the best-fitting model.

4. Conclusions

Our results show that both Baranyi and three-phase linear models provide low variability for growth rate values when analysing similar growth curves, hence being the models of choice. Three-phase linear model gives the lowest variation for growth rates, while Baranyi gives a variation marginally higher, despite much better overall fitting.

These results provide insight into predictive microbiology and will help food microbiologists and researchers to choose the proper primary growth predictive model.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research was financially supported by the Ministerio de Economía y Competitividad of the Spanish Government and Fondo Europeo de Desarrollo Regional (FEDER) through Projects AGL2010-19775 and AGL2013-48993-C2-1-R.