Considerable effort has been devoted to modelling the growth of pathogenic and spoilage bacteria in food. This is referred to as predictive microbiology. Mathematical equations are used to describe the effect of environment, for example temperature, or in more complex models temperature, pH, available water (Aw) and other factors that affect microbial behaviour.
The advantages of modelling have been described in the article on "Modelling in Food Technology" including reducing the costs and time required in determining the safe shelf life of new products or in undertaking pathogen challenge testing. However, caution and scientific expertise are required e.g. the Food Safety Authority of Ireland (2012) has cautioned the food industry on the use of predictive models.
The purpose of this article and supporting material is to illustrate how published research in predictive microbiology can be used in practice. This is an area in which Food Science and Food Technology undergraduates sometimes find difficulty. Part of the difficulty may exist because the steps in the calculations involved are not usually presented.
This article describes a model for the growth of salmonella on cut tomatoes and a calculator where you can enter the initial numbers of salmonella, the incubation temperature and the incubation time to obtain a prediction of final numbers.
Modelling microbial growth
Because of the critical importance of temperature early work attempted to modify the Arrhenius Law to describe microbial growth but this was either unsuccessful or the relationships derived were generally too complex for routine use. In a classic paper, Australian workers (Ratkowsky et al., 1982) proposed a relatively simple, two–factor empirical equation (equation 1) to describe the influence of temperature on microbial growth up to the maximum growth temperature of an organism, Tmax. This is often called the square root model.
Equation 1. √r= b (T-To).
r is the growth rate constant and denotes the number of generations of the organism produced in unit time, normally reciprocal hours, r is determined from regression analysis of growth curves of the organism concerned over selected temperatures up to Tmax.
T is temperature that can be recorded in either Celsius or kelvin.
The constant b is the slope of the regression line of √r versus temperature.
To is the temperature at which the regression line meets the temperature axis at √r =0. It is also a constant for a particular organism. Ratkowsky et al. (1982) defined To as “a conceptual temperature of no metabolic significance”.
Knowing the growth rate constant of a microorganism and the equation describing its relationship with initial numbers and temperature up to Tmax it is possible to model its growth at a particular temperature and after a designated time.
Let’s take a model from the literature to illustrate how to predict the numbers of a pathogen on a food after a specified time.
A number of Salmonella species can grow rapidly on fresh cut tomatoes and as a consequence there have been numerous incidents of food poisoning in the US and Europe from the consumption of infected tomatoes. Refer to Valdez et. al. (2012) for a review.
Pan and Schaffner (2010) have published a Ratkowsky model to describe salmonella growth as a function of storage temperature (equation 2) using several Salmonella species that are all human pathogens on cut tomatoes at temperatures ranging from 10°C to 35°C.
Equation 2. √Growth Rate = 0.026T- 0.1065
We will use this equation to predict the final number of salmonella produced from an inoculum of 10 cells / gram after 5 hours incubation at 22°C.
Using equation 2, the square root of the growth rate at 22 °C = 0.026 x 22 -0.1065
The growth rate = 0.46552
= 0.2167 generations h-1
Next the generation, or doubling time, is calculated. This is done by taking the reciprocal of the growth rate :
= 4.6147 hours/generation at 22°C.
The numbers of generations that have taken place in 5 hours are next calculated.
= 1.08345 generations
Next we will calculate the number of salmonella that have been produced after 3.4 generations.
= 10 x 2 1.08345
= 21.2 cells in 5 hours.
How accurate is this model?
The authors have indicated that the r2 value of the regression equation used to generate the specific growth rate was 0.9398. They also attempted to validate their work using the USDA Pathogen Modelling Program (PMP) and ComBase.
They reported that the PMP under-predicted growth at low temperatures and over predicted growth at high temperatures. ComBase predicted consistently slower growth rates than were observed in tomatoes, but showed parallel increases in growth rate with increasing temperature. However, the models tested were not developed using tomatoes and the discrepancies are therefore not unexpected.
Using ComBase to validate the microbial growth model
To complete this tutorial we will compare the prediction of salmonella growth using the equation derived by Pan and Schaffner (2010) with data obtained from ComBase.
To use ComBase registration is required. Uses register OnLine at https://www.combase.cc . Registration is free and once registered, the user logs on using a password and their Email address.
In this example we select ComBase Predictor and the growth options from the navigation menu on the left are completed (Figure 2).
Numbers are entered as log10 values, so the Initial Level box is filled using 1. A value of 1 for Physical State was selected on the assumption that the cells added were fully active and in the early log phase. The pH (4.6) and Aw (0.997) values used reflect values reported for tomato puree (Jakobsen, 1983). Note the pH and Aw values for tomatoes vary, pH values within the range 4.2 to about 4.6 have been reported. Labuza (1984) has reported Aw values for fresh tomatoes as varying from 0.991-0.998 and a value of 0.993 for pulped tomatoes.
Using these settings I used ComBase to obtain the number of salmonella predicted after 5 hours incubation at the temperatures shown in Table 1.
Table 1. Comparisons of salmonella growth on cut tomatoes at various temperatures using the ComBase microbial growth predictor and the Pan and Schaffner model
Numbers* of Salmonella species, CFU/g, on cut tomatoes after 5 hours at 22°C.
Pan and Schaffner model
* Based on an initial inoculum of 10 cells.
The results in Table 2 show that Combase Predictor gives higher growth rates for salmonella than the Pan and Schaffner model as found by Pan and Schaffner (2010).
EU standard for salmonella in food.
The current European Standard for salmonella in foods requires their absence in a 25 gram sample (European Commission, 2005).
Use the Pan and Schaffner model to predict salmonella numbers on cut tomatoes
Before using this model do note the widely accepted microbiological standard that salmonella should be absent in 25 grams of food. This does not mean that salmonella is absent! It just means that in the sample or samples studied there were none detectable in 25 grams.
When using the calculator try inputting 1 cell in say 1 kg of tomatoes and see what happens. You can also use 0.1 and lower numbers to simulate say 1 cell in 10 kg of tomatoes or lower. So, you can explore the consequences of temperature abuse over extended periods of time on tomatoes with low levels of initial microbial contamination.
The calculator was written in the Active Server Pages (ASP), also known as Classic ASP or ASP Classic, programming language. ASP is easy to learn and users can view the code and download the scripts for use either on their own PCs (you may need to install the free IIS or PWS components on your PC) or institution's Windows server. The coding in the scripts can easily be converted to PHP if required.
Mathematical models are very useful as discussed in the introduction, but some have certain limitations and often do not allow for the complex interactions that can occur in foods. Many models use broth or other media rather than foods as the modelling substrate and the results obtained require qualification and expert interpretation see for example Food Safety Authority of Ireland (2012).
I acknowledge helpful discussions with Professors Margaret Patterson and Donald Schaffner on predictive microbiology during the preparation of this tutorial. Any errors in this article are the responsibility of the author entirely.
European Commission. (2005). Commission Regulation (EC) no. 2073/2005 of 15 November 2005 on microbiological criteria for foodstuffs. Official J. Eur. Union L 338. 1–26.
Food Safety Authority of Ireland. (2012). Predictive Microbiology and Shelf-life. On Line. Available from: <https://www.fsai.ie/faq/shelf_life/predictive_microbiology.html>. Accessed: 28.04.2005.
Jakobsen, M. (1983). Filament hygrometer for water activity measurement: interlaboratory evaluation. Journal of the Association of Official Analytical Chemists. 66. 1106-1111.
Labuza,T.P. (1984). Moisture sorption: Practical Aspects of Isotherm measurement and use. St Paul, MN. AACC International Publishing.
Pan, W. and Schaffner, D.W. (2010). Modelling the Growth of Salmonella in Cut Red Round Tomatoes as a Function of Temperature. Journal of Food Protection. 8, 1408-1590.
Valadez, A. M., Schneider, K.R. and Danyluk, M.D. (2012). Outbreaks of Foodborne Diseases Associated with Tomatoes. Document number: FSHN12-08. Available from the EDIS website at https://edis.ifas.ufl.edu/fs192 .
On Line resources
Combase. Available from: https://www.combase.cc. Accessed: 30-03-2015.
Mullan, W.M.A. (2015). Calculator for predicting the growth of salmonella in cut tomatoes at 10°C to 35°C. [On-line]. Available from: https://www.dairyscience.info/newcalculators/sam.asp. Accessed: 30-03-2015.
Pathogen Modelling Program (PMP) Available from: https://pmp.errc.ars.usda.gov/PMPOnline.aspx. Accessed: 30-03-2015.
How to cite this article
Mullan, W.M.A. (2015). [On-line]. Available from: https://www.dairyscience.info/index.php/food-model/258-predict.html . Accessed: 2 April, 2020. Updated December 2019.