Original Research ARTICLE
From Culture-Medium-Based Models to Applications to Food: Predicting the Growth of B. cereus in Reconstituted Infant Formulae
- 1Department of Chemical and Food Engineering, Federal University of Santa Catarina, Florianópolis, Brazil
- 2Nestlé Research Center, Lausanne, Switzerland
- 3Institute of Nutrition, University of Debrecen, Debrecen, Hungary
Predictive models of the growth of foodborne organisms are commonly based on data generated in laboratory medium. It is a crucial question how to apply the predictions to realistic food scenarios. A simple approach is to assume that the bias factor, i.e., the ratio between the maximum specific growth rate in culture medium and the food in question is constant in the region of interest of the studied environmental variables. In this study, we investigate the validity of this assumption using two well-known link functions, the square-root and the natural logarithm, both having advantageous properties when modeling the variation of the maximum specific growth rate with temperature. The main difference between the two approaches appears in terms of the respective residuals as the temperature decreases to its minimum. The model organism was Bacillus cereus. Three strains (B594, B596, and F4810/72) were grown in Reconstituted Infant Formulae, while one of them (F4810/72) was grown also in culture medium to calculate the bias factor. Their growth parameters were estimated using viable count measurements at temperatures ranging from 12 to 25°C. We utilized the fact that, if the bias factor is independent of the temperature, then the minimum growth temperature parameter of the square-root model of Ratkowsky et al. (1982) is the same for culture medium and food. We concluded, supported also by mathematical analysis, that the Ratkowsky model works well but its rearrangement for the natural logarithm of the specific growth rate is more appropriate for practical regression. On the other hand, when analyzing mixed culture data, available in the ComBase database, we observed a trend different from the one generated by pure cultures. This suggests that the identity of the strains dominating the growth of mixed cultures depends on the temperature. Such analysis can increase the accuracy of predictive models, based on culture medium, to food scenarios, bringing significant saving for the food industry.
Bacillus cereus is a Gram positive, spore-forming, facultative anaerobic, rod-shaped pathogen (Kotiranta et al., 2000). Bacillus cereus strains are ubiquitous in the environment and can be isolated from soil, water and vegetables (Althayer and Sutherland, 2006; El-Arabi and Griffiths, 2013). They are commonly found in raw materials and processed foods, such as rice, milk and dairy products, meat and meat products, pasteurized liquid egg, ready-to-eat vegetables, and spices (Ceuppens et al., 2011). Bacillus cereus can be responsible for food poisoning illnesses in two ways: by heat labile, diarrhea-causing enterotoxins produced in the small intestine, and by heat stable cereulide toxin produced in the food before ingestion (Ceuppens et al., 2011), causing emetic syndromes. Emetic strains are, therefore, a concern to the food industry since it is not possible to eliminate the cereulide once performed in the food. Growth and toxin production must be strictly controlled, especially in food targeted to sensitive populations.
Bacillus cereus can endure ultrahigh-temperature (UHT) pasteurization, resist spray drying and survive in final products (Mcauley et al., 2014). Moreover, according to a review published by the European Food Safety Agency (EFSA, 2005), B. cereus strains are highly variable in terms of their tolerance to high temperatures and their ability to grow. This is mainly dependent on their phylogenetic group (Carlin et al., 2013). Mathematical modeling can be a valuable tool to assess and quantify this variability. It is widely accepted that temperature is the most important external environmental factor affecting microbial response. Among the available predictive models, the model of Ratkowsky et al. (1982) is commonly used to predict the maximum specific growth rate in the suboptimal region of temperature.
However, developing and validating a new model to predict microbial behavior during the manufacturing or the shelf life of a food commodity require extensive experimental work. It is a good practice to use readily available published data and models in the literature or in user-friendly predictive microbiology tools. For example, ComBase (http://www.combase.cc) provides culture-medium-based predictive models for a large collection of micro-organisms including B. cereus. To establish a “correction factor” that could be used to predict the behavior of the organism in food from culture-medium based models would be valuable for the food industry. To quantify the similarity between prediction and observation, the accuracy and bias factors, Af and Bf, respectively, of Ross (1996) is commonly used for practical applications. The indicators ln(Af) and ln(Bf) are certain average differences between the natural logarithm of the predicted and observed ln(μ) values of the organism in the studied range of environmental conditions, where μ denotes the maximum specific growth rate under a given condition. In the case of Af, the difference is meant as an absolute value, while in the case of Bf it is meant with its sign. Consequently, a bias factor Bf = 1 means that, in a studied region, on average, the model predictions are neither over-estimating nor under-estimating the growth rate compared to the observations. However, this could happen in such a way, too, that the predictions are underestimations in one part of the region while they are overestimations in the other part. It would be desirable that, for a given food matrix, the bias factor is independent of the environmental conditions, primarily of the temperature, at least in the normal physiological growth region of the organism. In this case, culture-medium-based predictions, available from public software tools such as ComBase, could be readily applied to the food in question. Since culture medium broths provide optimal substrate for the organism, the bias factor μfood/μbroth should normally be < 1.
The main objectives of this paper are (i) to provide a numerical analysis for the connection between bias and the two most frequently used transformations of the maximum specific growth rate parameter, the square-root and the logarithm functions; and (ii) to test whether the bias factor for B. cereus in Reconstituted Infant Formulae (RIF) can be considered constant, at least in a sub-optimal region of the temperature.
Materials and Methods
In laminar flow cabinet, infant formulae milk powder was weighted into sterile bottles, warm (~50°C) sterile water was aseptically added and then mixed to dissolve, according to manufacturer's instructions to obtain 50 ml of RIF samples.
Three emetic strains of Bacillus cereus were used in this study. B594 and B596 isolated from cereals and filed in the Nestlé culture collection and F4810/72 a reference strain from the DSMZ culture collection isolated during an outbreak investigation and also referred to as DSMZ4312 as reported in Carlin et al. (2013). Stock cultures were formed using subcultures of each strain supplemented with glycerol and stored at −80°C until used. For each strain, one tube of frozen stock culture was used to inoculate 10 ml of BHI (Sigma-Aldrich) and stored for 24 h at 30°C. Then 100 μl of this culture was put into another 10 ml of BHI and incubated for 18 h at 30°C. The subculture was then enumerated both on selective (Bacara, BioMérieux) and a non-selective (TSAyE, Sigma-Aldrich) media, diluted to a target level of 106 CFU/ml before applying a thermal stress during 25 s at 72°C. The plates were incubated for 24 h at 30°C. The stressed culture was also enumerated both on the selective and non-selective media to assess the stress intensity. This protocol allowed to mimic the processing conditions that influences the physiological state of naturally contaminating B. cereus cells.
Prior to inoculation, RIF bottles were equilibrated at the targeted temperatures (9, 12, 15, 18, 22, 25°C for F4810/72 strain and 9, 12, 15, 18, and 22°C for B594 and B596 strains). Appropriate dilution of the inoculum was then added to 50 ml of RIF to reach an initial concentration of 2.5 log CFU/ml. Inoculated bottles were sampled for viable counts on Bacara medium at appropriate sampling times to describe the different phases of the growth curves. Three independent replicates were performed for each experiment. For the reference strain F4810/72, additional experiments were performed in BHI following the same protocol to calculate the bias factor. All experiments were performed with pure culture.
Primary and Secondary Modeling
For each temperature, each curve was fitted by the primary model of Baranyi and Roberts (1994) using the DMFit Excel Add-in, downloadable from http://www.combase.cc/index.php/en/tools. As a second step, the effects of the environment, history and strain were studied and modeled on the maximum specific growth rate (μ), duration of lag phase (λ), and the natural logarithm of the maximum population density (ymax).
The maximum population density and the h0 = μ·λ parameters were taken as constant, as the simplest (zero-order) model, obtained via taking the multiplicative average of their estimates from the primary model.
The model of Ratkowsky et al. (1982),
in the reparameterized version
was fitted to the square root of the specific growth rates to describe the effect of temperature, where a and b are constants and Tmin = −b/a is a nominal minimum growth temperature, at which the extrapolated maximum specific growth rate would be zero.
We use the above well-established model in its second form, Equation (2), so (though it leads to non-linear regression), the important Tmin parameter and its standard error can directly be obtained. To carry out the non-linear fitting, the method of Levenberg–Marquardt (Press et al., 1986) was programmed in Visual Basic for Applications assigned to MS Excel.
We also tested the Ratkowsky model in a rearranged form, with the natural logarithm link function applied to the observed maximum specific growth rates:
A measure of the deviation between observed and predicted growth, called the “bias factor” was introduced by Ross (1996). As per definition, its natural logarithm, ln(Bf) is the average value between the observed and predicted ln(μ) values where μ denotes the maximum specific growth rate of the organism. It is of common sense to expect the conditions (here the temperature), under which the μ-values were generated, randomly distributed in the modeled region, in which case
where E denotes the “expected value” operator,. Since the μpred predictions produced by commonly used software packages are often based on experiments carried out in culture medium broth, while practical observations (μobs) refer to food, the above expectation can be translated to
In our case, the studied food matrix is RIF, for which a bias factor can be derived via the above formula from the growth rate in broth.
Primary and Secondary Modeling
Examples for fitting the primary model can be seen in Figure 1.
Figure 1. Growth of heat stressed B. cereus strain F4810/72, in RFI, at 9°C (diamond), 12°C (cross), 15°C (triangle), 18°C (square), 22°C (star), and 25°C (circle). The sigmoid model of Baranyi and Roberts (1994) was fitted (continuous lines) to the log count curves at each temperature.
The parameter estimates from the primary modeling are given in the Supplementary Information. Figure 2 shows all specific rate estimates for the three studied strains at the different temperatures in broth and in RIF.
Figure 2. Maximum specific growth rate estimates as a function of temperature for all analyzed strains.
The estimates for the b and Tmin parameters, when the Ratkowsky model was fitted to the maximum specific growth rates, are shown in Table 1. The Tmin-estimate for the strain B594 was significantly higher than the respective estimates for the other two strains.
Table 1. Estimated parameters and their standard errors for the square root model in RIF (3 strains) and BHI (one strain).
The Ratkowsky model claims two major benefits: the linear model structure for the “ vs. temperature” relationship, and the constant variance of the measured values, independently of the temperature. In our case, because of biological replicates (three), it was possible to study the relative deviations of the specific growth rates (standard deviation divided by the mean) within the replicate curve-triplets. These showed no correlation with the temperature (p = 0.65), see Figure 3.
Figure 3. Relative deviations (standard deviation divided by the mean) of replicate specific growth rate estimates, in RIF, at the studied temperatures. There is no correlation (p = 0.65) between the relative errors and the temperature.
This suggests that the natural logarithm could also be a suitable link function for the maximum specific growth rate when modeling its dependence on the temperature. This comes from the relationships:
where the approximation is accurate at least for one digit if the relative error ε is less than 0.3. For ε-values over 0.3, the approximation in Equation (8) would have worse than one digit accuracy. From this, it also follows that, since our average relative errors are less than 10%, the standard error of fit of the secondary model for ln(μ) will be ca 0.1 (or possibly higher, if the secondary model describes the “μ vs. T” relationship inaccurately).
It can be readily seen that if the ε random error in the Equations (6–8) is constant, independently of the temperature, then the same cannot hold for the square-root model and vice versa. Revisiting Equation (7), one can obtain, by first order approximation:
This means that, if the natural logarithm transformation makes the deviation of the observed specific rates constant, then the deviation generated by the square-root function should tend to be smaller with temperature decreasing to Tmin. That is, the absolute residuals should show a decreasing trend with the temperature (and, therefore, with the μ-values) - as indeed we will see it later. On the other hand, if the square-root was the correct transformation to stabilize the variance still the natural logarithm of the μ-values is regressed in the secondary model, then the residuals should show increasing trend as the temperature decreases to Tmin:
The maximum specific growth rate of the strain F4810/72 was measured in both RIF and BHI, providing a good opportunity to investigate the bias factor. Notice that, if this is independent of the temperature, then the secondary model for ln(μ) for the two media (Equation 3) differ only by an additive constant (Equation 5). This is equivalent to the assumption that the Tmin parameter is the same for the BHI and RIF. Our investigations showed that the Tmin of this same strain in BHI was 5.13 ± 1.12 which is not significantly different (p = 0.35) from the Tmin-value in RIF 5.40 ± 0.88 (Table 1).
Then we use the formulae
from which the ratio Bf = (bfood/bbroth)2 is constant, so the secondary models for ln(μ), for broth and food, should be parallel and differ from each other by the
constant additive term. The opposite direction of this conclusion can be proven similarly, therefore, the assumption that Tmin is the same for culture medium and food is equivalent to the one that the ratio between maximum specific growth rates in culture medium broth and food is constant, independently of the temperature. In our situation, we showed, by F-test, that the strain F4810/72 has the same Tmin = 5.26 value for BHI and RIF (p = 0.35). Therefore, their model can be written as
Substituting the coefficients above in Equation (14), it can be calculated that, for this strain, the ratio between the maximum specific growth rates in RIF and culture medium broth is Bf = (0.050/0.061)2 = 0.67. That is, this strain grows at one third slower in RIF compared to BHI.
Figure 4 demonstrates well the equivalence between the two assumptions: common Tmin-value for the square-root and parallel models for the logarithm link functions. The model (3) fitted to the ln(μ) values of the strain F4810/72 in BHI and RIF are parallel because they have similar Tmin estimates. As follows from Equation (9), the deviation from the parallel behavior would be apparent at low temperatures only.
Figure 4. Model of Ratkowsky et al. (1982) with two different link functions: (A) Square-root; and (B) natural logarithm. The model describes the effect of temperature on the maximum specific growth rate of the strain F4810/72 in RIF (circles) and in BHI (crosses). The (B) plot also shows data from Carlin et al. (2013) (stars) for the same strain in broth.
The shown equivalence is independent of the question whether the square-root or the logarithm transformation stabilizes the variance of the μ-values. According to the Equations (7) and (9), both cannot be valid at the same time. Comparing the absolute residuals for both the square-root and logarithm link function, on all the data, the Figure 5 emerges.
Figure 5. Absolute residuals vs. temperature for (A) the square-root and (B) natural logarithm link functions.
The residuals with the square-root link function show a decreasing trend as the temperature decreases to Tmin (p = 0.004), while in the case of logarithm link function, it does not show correlation with the temperature (p = 0.37). Therefore, based on our data, while the Ratkowsky model accurately describes how the maximum specific growth rate depends on temperature, the logarithm link function is more suitable to be applied to the observed maximum specific growth rates when regressing them against temperature. The difficulty is that this difference between the two link functions can be detected at low temperatures only, where it is not easy to keep the environment constant for the required long time to reach the stationary phase, therefore, the environmental effects (e.g., pH decrease in the medium) rather than biological ones (linked to strain variability for example) can dominate the variability of the observed maximum specific growth rates.
The paper of Carlin et al. (2013) gives an opportunity to compare the kinetic parameters of the reference strain F4810/72 in broth as shown in Figure 4. Fitting the square-root model to the 12–25°C data for the same strain in that paper, the estimated parameters were not different at 5% significance level (p = 0.12).
In the same way, we can validate our rate estimates by the ComBase Predictor available from http://www.combase.cc. In Figure 6, the square root values of our specific growth rates can be compared with results from ComBase Predictor, in broth and milk, assuming a bias factor of Bf = 0.67 for the food scenario.
Figure 6. Square root of maximum specific growth rates of B. cereus vs. temperature from different sources: Our B594, B596, and F4810/72 strains in RIF (open diamond, circle and triangle, resp.); Our data for F4810/72 in broth (plus sign); and ComBase Predictor results (continuous line). The dashed line below the ComBase Predictor curve is the prediction obtained by using the Bf = 0.67 bias factor.
The validation plot in Figure 6 is a convincing proof of the diversity of the kinetic behavior of B. cereus strains. The ComBase Predictor is based on growth curves generated by a cocktail of six strains (Sutherland et al., 1996). A plausible explanation for the seemingly unexpected non-linear behavior of the predictions is that different strains were the dominant ones at different temperatures, while the same parameter of a pure culture show a consistent linear pattern with temperature.
For quantitative validation, we made an extensive use of the Bias and Accuracy factors of Ross (1996). We point out here that while acknowledging the useful applicability of these indicators, their definition needs some refinement, in agreement with Baranyi et al. (1999). When the average of the ln(μfood) - ln(μpred) values is taken, it is implicitly assumed that the probability distribution of this difference is independent of the temperature and possibly other environmental factors (Gill and Phillips, 1985; Buchanan and Bagi, 1997; Neumeyer et al., 1997; Mellefont et al., 2003). The constant bias-factor is a reasonable assumption in case of the temperature, with the rationale that all affecting biochemical reactions speed up or slow dow when temperature changes. It is less obvious with other environmental factors, like pH or water activity; nonetheless the assumption provides a convenient approximation that is easy to build in predictive packages.
The assumption of the temperature-independent bias factor is equivalent to the existence of a minimum growth temperature that is the same for the model and for the food matrix on which the model is tested. Indeed, this latter condition has been assumed by quite a few authors (Miles et al., 1997; Carlin et al., 2013; Aryani et al., 2015, 2016), and was observed in our situation, too, when comparing the temperature-dependent maximum specific growth rates in RIF and culture medium. The ComBase database (http://www.combase.cc) also provides a good opportunity to check how much the temperature-independence of the Bias factor is tenable.
In conclusion, we agree with many authors (Bernaerts et al., 2000; Ross et al., 2003; Powell et al., 2015; Den Besten et al., 2017) that, at sub-optimal temperatures, the Ratkowsky model is a good representation of the effect of temperature on the maximum specific growth rate measured in a pure culture, in both laboratory medium and food. However, the constant variance assumption does not necessarily hold at low temperatures. Besides, we established that the minimum growth temperature can be taken as the same Tmin value for culture medium and food, therefore, the bias factor is, indeed, independent of the temperature. In mixed cultures, however none of the above holds, and more complex developments (data and mathematical considerations) are needed for an accurate model; see Baranyi et al. (2017), which is, in a sense, a continuation of this paper.
NBdS: experiments and writing up. JB Experimental design, data analysis and writing up. BAMC Consultation. ME conception, experimental design and writing up.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The experiments were carried out at the Institute of Food Research, Norwich, UK. The authors wish to express their thanks to Susie George and Amanda Demeter for their help in the laboratory. NBdS thankfully acknowledges the scholarship from the Brazilian Coordination for the Improvement of Higher Education Personnel (CAPES).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2017.01799/full#supplementary-material
Althayer, M., and Sutherland, A. D. (2006). Bacillus cereus is common in the environment but emetic toxin producing isolates are rare. J. Appl. Microbiol. 100, 7–14. doi: 10.1111/j.1365-2672.2005.02764.x
Aryani, D. C., Den Besten, H. M. W., Hazeleger, W. C., and Zwietering, M. H. (2015). Quantifying strain variability in modeling growth of Listeria monocytogenes. Int. J. Food Microbiol. 208, 19–29. doi: 10.1016/j.ijfoodmicro.2015.05.006
Aryani, D. C., Zwietering, M. H., and Den Besten, H. M. W. (2016). The effect of different matrices on the growth kinetics and heat resistance of Listeria monocytogenes and Lactobacillus plantarum. Int. J. Food Microbiol. 238, 326–337. doi: 10.1016/j.ijfoodmicro.2016.09.012
Baranyi, J., Buss da Silva, N., and Ellouze, M. (2017). Relationships between the parameters of temperature-dependent growth models for Bacillus cereus strains. Front. Microbiol. 8:1890. doi: 10.3389/fmicb.2017.01890
Bernaerts, K., Versyck, K. J., and Van Impe, J. F. (2000). On the design of optimal dynamic experiments for parameter estimation of a Ratkowsky-type growth kinetics at suboptimal temperatures. Int. J. Food Microbiol. 54, 27–38. doi: 10.1016/S0168-1605(99)00140-3
Carlin, F., Albagnac, C., Rida, A., Guinebretière, M. H., Couvert, O., and Nguyen-The, C. (2013). Variation of cardinal growth parameters and growth limits according to phylogenetic affiliation in the Bacillus cereus Group. consequences for risk assessment. Food Microbiol. 33, 69–76. doi: 10.1016/j.fm.2012.08.014
Ceuppens, S., Rajkovic, A., Heyndricks, M., Tsilia, V., Van De Wiele, T., Boon, N., et al. (2011). Regulation of toxin production by Bacillus cereus and its food safety implications. Crit. Rev. Microbiol. 37, 188–213. doi: 10.3109/1040841X.2011.558832
Den Besten, H. M. W., Aryani, D. C., Metselaar, K. I., and Zwietering, M. H. (2017). Microbial variability in growth and heat resistance of a pathogen and a spoiler: all variabilities are equal but some are more equal than others. Int. J. Food Microbiol. 240, 24–31. doi: 10.1016/j.ijfoodmicro.2016.04.025
Gill, C. O., and Phillips, D. M. (1985). The effect of media composition on the relationship between temperature and growth rate of Escherichia coli. Food Microbiol. 2, 285–290. doi: 10.1016/0740-0020(85)90010-3
Mcauley, C. M., McMillan, K., Moore, S. C., Fegan, N., and Fox, E. M. (2014). Prevalence and characterization of foodborne pathogens from Australian dairy farm environments. J. Dairy Sci. 97, 7402–7412. doi: 10.3168/jds.2014-8735
Mellefont, L. A., McMeekin, T. A., and Ross, T. (2003). Performance evaluation of a model describing the effects of temperature, water activity, pH and lactic acid concentration on the growth of Escherichia coli. Int. J. Food Microbiol. 82, 45–58. doi: 10.1016/S0168-1605(02)00253-2
Miles, D. W., Ross, T., Olley, J., and McMeekin, T. A. (1997). Development and evaluation of a predictive model for the effect of temperature and water activity on the growth rate of Vibrio parahaemolyticus. Int. J. Food Microbiol. 38, 133–142. doi: 10.1016/S0168-1605(97)00100-1
Neumeyer, K., Ross, T., Thomson, G., and McMeekin, T. A. (1997). Validation of a model describing the effects of temperature and water activity on the growth of psychrotrophic pseudomonads. Int. J. Food Microbiol. 38, 55–63. doi: 10.1016/S0168-1605(97)00090-1
Powell, S. M., Ratkowsky, D. A., and Tamplin, M. L. (2015). Predictive model for the growth of spoilage bacteria on modified atmosphere packaged Atlantic salmon produced in Australia. Food Microbiol. 47, 111–115. doi: 10.1016/j.fm.2014.12.001
Ross, T., Ratkowsky, D. A., Mellefont, L. A., and McMeekin, T. A. (2003). Modelling the effects of temperature, water activity, pH and lactic acid concentration on the growth rate of Escherichia coli. Int. J. Food Microbiol. 82, 33–43. doi: 10.1016/S0168-1605(02)00252-0
Sutherland, J. P., Aheme, A., and Beaumont, A. L. (1996). Preparation and validation of a growth model for Bacillus cereus: the effects of temperature, pH, sodium chloride and carbon dioxide. Int. J. Food Microbiol. 30, 359–372. doi: 10.1016/0168-1605(96)00962-2
Keywords: Bacillus cereus, predictive microbiology, bias factor, reconstituted infant formulae, Ratkowsky model
Citation: Buss da Silva N, Baranyi J, Carciofi BAM and Ellouze M (2017) From Culture-Medium-Based Models to Applications to Food: Predicting the Growth of B. cereus in Reconstituted Infant Formulae. Front. Microbiol. 8:1799. doi: 10.3389/fmicb.2017.01799
Received: 03 May 2017; Accepted: 05 September 2017;
Published: 21 September 2017.
Edited by:Jean-Christophe Augustin, Ecole Nationale Vétérinaire d'Alfort, France
Reviewed by:Olivier Couvert, Laboratoire Universitaire de Biodiversité et Ecologie Microbienne (LUBEM), France
Maria Gougouli, Perrotis College, American Farm School, Greece
Copyright © 2017 Buss da Silva, Baranyi, Carciofi and Ellouze. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mariem Ellouze, firstname.lastname@example.org