This study aimed to investigate the growth of Listeria monocytogenes in rice balls and to conduct its microbial risk assessment based on the Korean dietary pattern. Each tuna or ham rice ball was mixed with mayonnaise, soy sauce, or gochujang, a Korean traditional fermented red peeper paste, which was artificially contaminated with L. monocytogenes and then stored at 7°C–25°C to assess bacterial growth. Growth data were analyzed using three primary models (the Huang, Baranyi, and Gompertz models), and the growth pattern was found to fit well to the Baranyi model based on the following five statistical criteria: root mean square error (0.38–0.56), Akaike’s information criterion (−51.55–−26.99), coefficient of determination (0.72–0.97), bias factor (0.97–1.01), and accuracy factor (1.06–1.18). The effects of temperature on bacterial growth rate and lag time were evaluated using the square root model. The minimum growth temperature for L. monocytogenes in tuna or ham rice balls was the lowest when they were mixed with mayonnaise (−9.44°C or −15.37°C, respectively). Risk assessment using FDA-iRISK showed that tuna or ham rice balls mixed with gochujang exhibited the highest microbial risk among all the rice balls tested, regardless of the storage temperature. Tuna or ham rice balls mixed with gochujang had the highest disability-adjusted life years per year (0.015) followed by ham rice balls mixed with soy sauce (0.011–0.015) or mayonnaise (0.006–0.015) and then tuna rice balls mixed with soy sauce (0.006–0.008) or mayonnaise (<0.001). In conclusion, our results, determined using predictive growth models, allow the assessment of potential risk ranking associated with the consumption of rice balls contaminated with L. monocytogenes based on the number of illnesses experienced per serving and the disease burden.
1. Introduction
Listeria monocytogenes is widely distributed in the natural environment, grows at a wide range of temperatures from 0°C to 45°C, and is resistant to a variety of environmental conditions such as high salinity (10%–12% NaCl) [1], low water activity (minimum 0.83) [2], and low pH (minimum 4.4) [3]. According to the Centers for Disease Control and Prevention, L. monocytogenes is accountable for approximately 1,591 cases of listeriosis, 1,455 hospitalizations, and 255 deaths annually in the United States [4]. In general, compared with healthy individuals, the risk of fatal invasive listeriosis is high in vulnerable populations, including pregnant women, newborns, elderly, and immunocompromised patients, resulting in a death rate of 20%–30% [5, 6]. Reportedly, listeriosis is primarily associated with the consumption of refrigerated ready-to-eat (RTE) foods as these foods are not pasteurized [5, 7]. Thus, controlling the contamination of RTE foods with L. monocytogenes is crucial in the food industry.
In Korea, rice balls and kimbab are popular convenience foods and are consumed at a high rate of 54.4% among convenience food consumption as a substitute for meals [8, 9]. Both are made by combining cooked rice with side ingredients such as ham, tuna, meat, cheese, or vegetables. However, these are distinguished in their shape and the number of side ingredients. Korean rice balls are made from rice and fillings such as a few minced ingredients mixed with dressings formed into oval or triangular shapes. Kimbab is a seaweed rice roll containing several ingredients. According to the Institute of Food Technology [10], convenience foods such as rice balls, kimbab, and sushi that contain cooked rice and various other ingredients are considered potentially hazardous foods and require control of storage time and temperature during their shelf-life. Such foods have been reported to cause many cases of foodborne diseases linked to enteric bacteria in Korea [11]. Foodborne bacteria were detected in 21 out of 77 commercial kimbab as a result of hygiene inspection of convenience foods [12]. There are no reported cases on foodborne disease caused from L. monocytogenes in rice balls and kimbab in Korea. However, its outbreak possibility is high due to the rapid growth of HMR (home meal replacement) and CMR (convenient meal replacement) markets.
Several researchers have developed predictive growth models and quantitative microbial risk assessments (QMRAs) for Staphylococcus aureus in kimbab [13, 14]. L. monocytogenes is another significant pathogen present in refrigerated convenience foods such as kimbab, rice balls, and sushi owing to its ability to grow at low temperatures; however, risk assessment studies on this pathogen have been relatively limited. Rice balls, which contain a single main ingredient (usually tuna, ham, or sauce) other than rice, slightly differ from kimbab, which contains a number of ingredients. In Japan, the contamination risk of RTE foods containing raw minced tuna as a main ingredient with L. monocytogenes was 5.7%–12.1% [15]; moreover, a survey revealed that the contamination risk of sausages and ham with L. monocytogenes was ranked in top 10 on the necessity of priority management targeting Korean livestock coops [16]. Thus, in this study, the risk of L. monocytogenes was assessed in tuna or ham rice balls containing one of three popular sauces (mayonnaise, soy sauce, or gochujang) each as the main ingredient.
Predictive growth models are mathematical expressions that can analyze and predict the growth and survival of microorganisms for assessing food safety [17]. They can be useful in decision-making with respect to critical control points (CCP) in hazard analysis and CCP (HACCP) systems and microbial risk assessment for preventing foodborne illnesses and food spoilage [18]. The growth pattern of L. monocytogenes in foods is dependent on nutrients, water activity, pH, the presence of antagonists and synergists, redox potentials in the food matrix, and the environment in which food is stored [19]. The Integrated Pathogen Modeling Program, 2013 (IPMP, 2013) of the United States Department of Agriculture used in this study is one of several predictive software tools that models the growth of microorganisms. QMRA is a scientific method for assessing the risk of foodborne illness after the consumption of foods contaminated with certain pathogenic microorganisms [10, 20]. It includes four steps: hazard identification, exposure assessment, hazard characterization, and risk characterization. Therefore, QMRA can provide food safety standards for foodborne pathogenic microorganisms by estimating the risk of illness due to a certain pathogen in a certain population.
Food and Drug Administration- (FDA-) iRISK, which was developed by the US FDA, is an interactive, web-based, comparative risk assessment tool. It can calculate and determine risk ranks using mathematical logic and Monte Carlo simulations and using data on the number of diseases and then translate the result into public health metrics that facilitate the comparisons of public health burdens [21]. In this study, we used predictive modeling to investigate the growth characteristics of L. monocytogenes in various types of rice balls stored at different temperatures, and the predicted growth parameters obtained from a suitable model were used in risk assessment. QMRA of L. monocytogenes in various rice balls was then conducted using FDA-iRISK to compare, evaluate, and rank risk depending on different environmental conditions.
2. Materials and Methods2.1. Preparation of Rice Balls
Salad mixtures of tuna or ham were first prepared by mixing canned minced tuna (Dongwon F&B Co., Ltd., Changwon, Korea) with mayonnaise (Ottogi, Eumseong, Korea), soy sauce (Sempio, Incheon, Korea), or gochujang (CJ Cheiljedang Co., Ltd., Nonsan, Korea) in the ratios of 3 : 1, 15 : 1, or 10 : 1, respectively, or cubed ham (1 cm pieces; CJ Cheiljedang Co., Ltd., Seoul, Korea) with the same sauces in the ratios (weight : weight) of 5 : 1, 12 : 1, or 12 : 1, respectively. The bacterial cell number of initial rice balls was 2.26 log CFU/g before inoculation of L. monocytogenes. An overnight culture of L. monocytogenes (KCTC 3710) was inoculated into each salad mixture at approximately 3 log CFU/gram of rice balls. Each salad mixture was wrapped with cooked rice (Chucheong variety, Nonghyup, Anseong, Korea; ratio, 1 : 4), and each ball was sealed in a sterile plastic bag and stored at 7°C, 15°C, or25°C. L. monocytogenes was enumerated at each sampling time point by incubation on TSA (BD Difco, Sparks, MD) at 37°C for 24 h. All microbial analyses were conducted for each rice ball in four repetitions.
2.2. Primary Model
The growth rate (µmax and specific growth rate (SGR); log CFU/h) and lag time (LT; h) of L. monocytogenes were determined using three primary models (the Huang, Baranyi, and Gompertz models) using IPMP (IPMP, 2013) software.
The Huang model is described by the following equation:(1)Yt=Y0+Ymax−lneY0+eYmax−eY0e−µmaxBt,Bt=t+14ln1+e−4t−λ1+e4λ,where Y0, Ymax, and Y(t) are the natural logarithms of the initial bacterial population count, maximum count, and the count at time t, respectively; μmax is the maximum SGR; and λ is the lag-phase duration.
The Baranyi model is described by the following equation:(2)Yt=Y0+µmaxAt−ln1+expµmaxAt−1expYmax−Y0,At=t+1µmaxlnexp−µmaxt+exp−h0−exp−µmaxt−h0,where Y0, Ymax, Y(t), and µmax represent the same quantities as those in equation (1) and h0 is the physiological state of the microorganism under consideration. In the Baranyi model, the lag phase corresponds to h0/µmax, and h0 can be different and unique for each curve. In the Baranyi model, the initial growth curves observed at constant temperature conditions were analyzed to obtain each h0 value. After analyzing all the growth curves, the h0 values were averaged to produce a global value. All growth curves were then reanalyzed using the average value of h0 to determine µmax at each incubation temperature.
The Gompertz model is described by the following equation:(3)Yt=Y0+Ymax−Y0exp−expµmaxeYmax−Y0λ−t+1,where Y0, Ymax, Y(t), μmax, and λ represent the same quantities as those in equation (1).
2.3. Secondary Model
Temperature-dependent variations in the maximum SGR and LT were analyzed using square root and polynomial models, respectively, using Microsoft Excel 2016 (Microsoft Corp, Redmond, Washington). The equations used are shown in equations (4) to (6):(4)lnSGR:a+bT+cT2,(5)lnLT:a+bT+cT2,(6)Tmin:µmax=aT−Tmin,where SGR is in log CFU/h and LT is in hours. T is the temperature (ºC); Tmin is the intercept between the slope and the temperature axis; and a, b, and c are regression constants.
2.4. Model Validation
Model acceptability for L. monocytogenes in rice balls was evaluated based on the root mean square error (RMSE), coefficient of determination (R2), Akaike’s information criterion (AIC), bias factor, and accuracy factor.
RMSE is the difference between observed and predicted data. For RMSE, values close to zero are desirable as they indicate that the predicted values are very close to the observed values.
R2 is calculated using the following equation:(7)R2=1−∑ei2∑yi−y′2,where ei is the error of the predictive value, yi_is the predictive data, and ý is the average of predictive data.
AIC is calculated using the following equation:(8)AICc=−2lnlikelihood+2K+2KK+1n−K−1,where n is the number of observations and K is the number of parameters in the model. Bf and Af were calculated using equations (9) and (10), respectively:(9)Bf=10∑logpred/obs/n,(10)Af=10∑logpred/obs/n,where obs, pred, and n correspond to the observed value, predicted value, and repetition number of the observed data, respectively. Values near to 1 indicate the suitability of the model, whereas the values <0.7 or >1.5 indicate the unsuitability of the model.
2.5. Risk Assessment for Rice Balls Using FDA-iRISK
Results of risk assessment of all tested rice balls were compared using FDA-iRISK (version 4.0, FDA, US).
3. Results and Discussion3.1. Growth of L. monocytogenes in Tuna or Ham Rice Balls Mixed with Three Sauces
The initial number of L. monocytogenes cells in all rice balls was approximately 3.60 ± 0.36 log CFU/g; growth curves of L. monocytogenes in six rice balls at different storage temperatures are shown in Figure 1. Bacterial growth in the food matrix is dependent on several factors such as the presence of other microbiota, nutrients, pH, temperature, and water activity [22, 23]. Lower temperatures retarded bacterial growth in all rice balls. Distinctly different growth pattern in LTM (tuna rice balls mixed with mayonnaise; Figure 1(a)) at 7°C from other samples was observed. Mixing mayonnaise negatively affected bacterial growth in both tuna and ham rice balls. The growth rate of L. monocytogenes in LTM was slower than that in other samples, and its maximal growth did not exceed 6.5 log CFU/g. Additionally, mixing mayonnaise in ham rice balls resulted in the lowest maximum growth (<7 log CFU/g) of LHM (ham rice balls mixed with mayonnaise). The low growth in mayonnaise may be due to the presence of egg-derived antibacterial substances such as lysozyme and ovotransferrin, although these are present at very low levels [24, 25]. In general, L. monocytogenes had a higher growth rate in ham rice balls than in tuna. The maximum population of L. monocytogenes at 25°C was 4.5 log CFU/g in LTM, 6.88 log CFU/g in LTS, 7.38 log CFU/g in LTG, 6.84 log CFU/g in LHM, 8.00 log CFU/g in LHS, and 8.18 log CFU/g in LHG. The pH values of all rice balls decreased (from 5.38 to 5.12 in LTM; from 5.40 to 5.09 in LTS; from 5.45 to 5.07 in LTG; from 5.63 to 5.52 in LHM; from 5.55 to 5.47 in LHS; and from 5.67 to 5.59 in LHG) during storage at 7°C for 137 h. pH of tuna rice balls was lower than ham rice balls.
Growth of L. monocytogenes in rice balls by fitting the Baranyi model: tuna rice balls mixed with mayonnaise (a), soy sauce (b), and gochujang (c) and ham rice balls mixed with mayonnaise (d), soy sauce (e), and gochujang (f). Symbols, experimental data; line, fitted model.
3.2. Primary Modeling of L. monocytogenes Growth in Rice Balls
Many researchers agree that no single growth model can produce a consistent goodness of fit for all experimental data [26–29]; therefore, most data are assessed using several primary models such as Gompertz, Baranyi, and Huang models (Supplementary Table 1). The microbial growth models such as Gompertz, Baranyi, and Huang models are basically similar. These have been used to calculate microbial growth parameters including lag time, specific growth rate, and maximum population density to estimate the microbial growth in food [29–31]. The results for microbial growth parameters are dependent on which primary model is employed due to the differences of the equations produced from each primary model. The Gompertz model is an empirical equation, the Baranyi model is semitheoretical, and the Huang equation is a theoretical model developed. The results for growth parameters sometimes show that these models were not significantly different [31]. Therefore, they were equally capable of explaining microbial growth in food. In our study, statistical values including R2, RMSE, and AIC from all the models were similar for the predictive growth of L. monocytogenes in rice balls. In general, the Baranyi model was a better model for predicting L. monocytogenes in rice balls than the other models as it exhibited the lowest RMSE and AIC and the highest R2. However, the Huang model was suitable for LTS at 15°C and LHG at 7°C, whereas the Gompertz model was acceptable for LTG at 25°C and LHG at 25°C. Even in the same food matrix, the model suitability differed depending on temperature [27]. When compared with the average values for each temperature, the Baranyi model produced the smallest average values for RMSE, MSE, and AIC and the largest R2 in most rice balls except in those containing gochujang. LTG and LHG were most accurately modeled using the Gompertz model, but there was no significant difference between the three models. Therefore, we concluded that the Baranyi model was the best primary model for determining the growth of L. monocytogenes in rice balls at all temperatures, and so, this model was used to predict the growth pattern of L. monocytogenes in rice balls (Figure 1). The experimental observations coincided well with the predicted values, which also predicted the lag and exponential phases.
3.3. Secondary Models for L. monocytogenes and Validation
SGR and LT predicted by the Baranyi model were used in a polynomial equation to fit to secondary models. These models were developed to describe the effect of temperature on the SGR and LT, and these equations estimated the SGR and LT of L. monocytogenes as a function of temperature (Table 1). Based on SGRs and LTs obtained, mayonnaise was determined as the sauce that most interfered with the growth of L. monocytogenes followed by soy sauce and gochujang. As the storage temperature increased, SGR increased (0.02–0.08 at 7°C, 0.06–0.20 at 15°C, and 0.15–0.38 at 25°C). As storage temperature increased, LT decreased (42.55–278.53 h at 7°C, 3.71–93.65 h at 15°C, and 6.96–14.41 h at 25°C) under most growth conditions. These data demonstrate that the growth of the microorganisms is dependent on storage temperature. A short LT was observed for the bacteria in rice balls mixed with gochujang, indicating easy growth. Temperature has been reported to be an important environmental variable affecting microbial growth [32, 33].
Secondary models developed for the specific growth rate (SGR, log CFU/g/h) and lag time (LT, h) of L. monocytogenes in rice balls and statistical validation factors calculated for Baranyi models.
Code
Equation
R2
Tmin (°C)
Storage temperature
7°C
15°C
25°C
LTM1)
SGR = −0.0109 + 0.0127 × T + 0.0001 × T2
0.99
−9.43
0.08 ± 0.004
0.20 ± 0.003
0.36 ± 0.02
LT = 540.2198–43.118 × T + 0.8826 × T2
0.99
278.53 ± 15.34
93.62 ± 1.30
14.41 ± 0.59
LTS
SGR = –0.0155 + 0.0039 × T + 0.0001 × T2
0.97
−3.02
0.02 ± 0.00
0.06 ± 0.01
0.15 ± 0.02
LT = 87.2021–5.5066 × T + 0.0133 × T2
0.99
52.07 ± 0.00
21.63 ± 3.56
−6.96 ± 0.75
LTG
SGR = 0.0889–0.0189 × T + 0.0012 × T2
0.99
2.67
0.02 ± 0.002
0.08 ± 0.004
0.38 ± 0.01
LT = 1.7491 + 6.3248 × T−0.24796 × T2
0.98
30.75 ± 3.77
40.38 ± 1.98
4.91 ± 0.14
LHM
SGR = −0.0721 + 0.0217 × T−0.0004 × T2
0.95
−15.40
0.06 ± 0.001
0.16 ± 0.008
0.21 ± 0.03
LT = 530.4548–47.236 × T + 1.0475 × T2
0.95
250.00 ± 5.35
58.05 ± 2.68
4.35 ± 0.48
LHS
SGR = 0.0325–0.0034 × T + 0.0004 × T2
0.98
−3.34
0.03 ± 0.001
0.07 ± 0.003
0.19 ± 0.01
LT = 287.397–30.0968 × T + 0.7449 × T2
0.99
115.06 ± 5.14
3.71 ± 0.39
0.55 ± 0.04
LHG
SGR = 0.0418–0.0049 × T + 0.0004 × T2
0.99
−2.73
0.03 ± 0.002
0.07 ± 0.006
0.20 ± 0.007
LT = 87.4372–7.5618 × T +0.1591 × T2
0.99
42.55 ± 2.53
9.75 ± 0.87
−2.19 ± 0.07
1) LTM, LTS, and LTG are tuna rice balls mixed with mayonnaise, soy sauce, and gochujang contaminated by L. monocytogenes; LHM, LHS, and LHG are ham rice balls mixed with mayonnaise, soy sauce, and gochujang contaminated by L. monocytogenes.
R2 was estimated to evaluate the secondary models; higher R2 reflected more accurate model predictions. In the Baranyi model, which was the most accurate among the three primary models, R2 for SGR and LT model equations was >0.95, indicating good fit. Various developed models have been shown to be capable of predicting the growth of L. monocytogenes under the certain conditions [34]. Tmin is the minimum temperature that facilitates bacterial growth estimated by the traditional Ratkowsky square model; in this study, we identified slight differences in Tmin between tuna and ham rice balls, with ham rice balls exhibiting a lower value than tuna rice balls. Particularly, tuna or ham rice balls mixed with mayonnaise exhibited the lowest Tmin (−9.43°C and −15.40°C, respectively) followed by those mixed with soy sauce and gochujang. Tmin of L. monocytogenes in rice balls is much lower than other reports [34–37]. Tmin of L. monocytogenes grown in culture media under optimal conditions was reported from 0.5°C to 5°C [38]. Tmin for pathogenic bacteria is dependent on the food type owing to differences in nutrients, pH, water activity, salt content, and the presence of antagonists or synergists. The growth of L. monocytogenes at low temperatures has also been reported by Devlieghere et al. [36] and Franz et al. [37], who showed that the bacterium can grow at −4.6°C and −3.54°C in green leafy vegetables and cooked ham substrates, respectively. Tmin of L. monocytogenes is substantially lower than that of other pathogenic bacteria such as Staphylococcus aureus in kimbab (2.5°C) [39] and E. coli in beef meat (1.5°C to 4.7°C) [40].
3.4. Risk Assessment of Rice Balls Using FDA-iRISK
Table 2 presents the model inputs for risk scenarios of L. monocytogenes in rice balls. Due to the absence of the reference literature for rice balls, reports on kimbab were used for determining the initial prevalence of L. monocytogenes [41]; among the 30 kimbab samples from local markets in Korea, 83.3% samples were contaminated with pathogenic bacteria including L. monocytogenes (6.7%). Therefore, 0.067 was used as the input data for the initial prevalence. The consumption model included grams per eating occasion and annual eating occasions for kimbab [42, 45]. The initial unit mass was the average serving size of commercial rice balls sold in convenience stores. Data from the Baranyi model were included in FDA-iRISK during the process stage and were designated as “increased” in software. The dose-response model and disease burden were based on a study by Pouillot et al. [43] and Noordhout et al. [44]. The DALY (disability-adjusted life years) is one of the health impact metrics that integrates information on the severity and duration of illness to estimate disease burden [46, 47]. DALY is measured as the sum of years of life lost by premature mortality and years lived with disability. It considers both the probability of experiencing an illness, injury, or even death and the impact of the associated health effect. The WHO estimated that listeriosis resulted in 23,150 illness and 172,823 DALY worldwide in 2010 [48]. Thus, we calculated the DALY per case as 7.46 following DALY divided by illness and listeriosis cases.
Model inputs for food scenarios of L. monocytogenes in rice balls in FDA-iRISK 4.0.
Korea Health Statistics: Korea National Health and Nutrition Examination Survey (2016)
Dose-response model
Healthy people
7.90E − 12
Pouillot et al. [43]
Pregnant women
2.01E − 09
Elderly (more than 65 years old)
1.49E − 10
Disease burden
Healthy people
0.015 DALY per case
Noordhout et al. [44]
Pregnant women and elderly
12.91 DALY per case
1) LTM, LTS, and LTG are tuna rice balls mixed with mayonnaise, soy sauce, and gochujang contaminated by L. monocytogenes; LHM, LHS, and LHG are ham rice balls mixed with mayonnaise, soy sauce, and gochujang contaminated by L. monocytogenes. 2) Input data for the initial concentration input: minimum and maximum concentrations at 7°C/minimum and maximum concentrations at 15°C/minimum and maximum concentrations at 25°C. 3) Process stage input: minimum and maximum concentrations at 7°C/minimum and maximum concentrations at 15°C/minimum and maximum concentrations at 25°C.
Table 3 shows the predictions obtained from FDA-iRISK, including total number of illnesses, mean risk of illness, and total DALY per year for L. monocytogenes in tuna or ham rice balls mixed with the three different sauces and stored at different temperatures. In this study, the risk of illness from L. monocytogenes in rice balls varied with the main ingredients (i.e., tuna, ham, and the type of sauces) and storage temperature. Their risk rates ranged from 0.000003 to 0.041 based on the mean risk of illness and from 0.00000093 to 0.015 based on the total DALY per year. The highest risk, regardless of the storage temperature, occurred when rice balls were mixed with gochujang. As the storage temperature decreased, the risk of illness from L. monocytogenes did not decrease. The risks from LHS and LHM were also highest at 7°C. Despite the low storage temperature at 7°C, the risk of L. monocytogenes from rice balls mixed with mayonnaise or soy sauce was higher than that at 25°C. Variations in the microbiota population within a rice ball at room temperature may explain the interrupted growth of L. monocytogenes under nonsterilized condition compared with that of L. monocytogenes alone at a refrigeration temperature [49]. The growth of L. monocytogenes in ham rice balls was higher than that in tuna rice balls, even under refrigeration. In general, the use of mayonnaise as a sauce reduced the estimated risk from L. monocytogenes. Tuna or ham rice balls mixed with gochujang had the highest value for the mean risk of illness (4.1 cases per hundred servings) regardless of temperature. Ranking of the other products was the same as for total illness, with LTM having the lowest mean risk of illness (3 cases per million servings) stored at 7°C. However, this was still a much higher risk when compared with listeriosis caused by the consumption of cooked meat (0.000629 cases per million servings) and hot-smoked fish (0.000000989 cases per million servings), as reported by the EFSA [7]. These higher risks may have been due to the large serving size and intake frequency of rice balls in Korea. The total DALY per year, i.e., the annual health burden, of rice balls mixed with gochujang was 0.015 at all storage temperatures, which was the highest identified in this study. LHM exhibited the lowest total DALY (0.006) among ham rice balls. Hong, Kim, and Yoon [50] reported that the DALY of various hams contaminated with Campylobacter jejuni ranged from 0.000208 to 0.000346. LTM exhibited the lowest total DALY among all tested rice balls. The susceptibility to listeriosis varies within different population groups; moreover, the risk of illness from consuming rice balls contaminated with L. monocytogenes is much higher in pregnant (0.20–20.76) and elderly groups (0.02–20.76; Supplementary Table 2) than in the general population.
FDA-iRISK output: risk ranking of six rice balls stored at 7, 15, and 25°C.
Code
Temperature (°C)
Total illnesses2)
Mean risk of illness3)
Total DALY per year
LTM1)
7
0.0001b
0.000003
0.00000093
15
0.075
0.003
0.001
25
0.0005
0.00002
0.004
LTS
7
0.432
0.018
0.006
15
0.531
0.022
0.008
25
0.404
0.017
0.006
LTG
7
0.978
0.041
0.015
15
0.978
0.041
0.015
25
0.978
0.041
0.015
LHM
7
0.978
0.041
0.015
15
0.831
0.035
0.012
25
0.392
0.016
0.006
LHS
7
0.978
0.041
0.015
15
0.796
0.033
0.012
25
0.712
0.030
0.011
LHG
7
0.978
0.041
0.015
15
0.978
0.041
0.015
25
0.978
0.041
0.015
1) LTM, LTS, and LTG are tuna rice balls mixed with mayonnaise, soy sauce, and gochujang contaminated by L. monocytogenes; LHM, LHS, and LHG are ham rice balls mixed with mayonnaise, soy sauce, and gochujang contaminated by L. monocytogenes. 2) Total illness is the total number of illnesses generated from the scenario. The total lifespan in the scenario is 70 years (duration of exposure). 3) Mean risk of illness represents the average probability of disease based on one serving.
4. Conclusion
In this study, the growth of L. monocytogenes in various rice balls was found to differ depending on the main food ingredients, such as tuna or ham, and sauces, such as mayonnaise, soy sauce, and gochujang. Additionally, the Baranyi model was identified as the optimal model to describe L. monocytogenes growth in rice balls and can be useful for conducting future risk assessment. The relative risk of illness from six rice balls was compared using FDA-iRISK based on the dose-response model, consumption, and DALY. Overall, ham rice balls presented a higher risk of listeriosis compared with tuna rice balls. L. monocytogenes in rice balls containing mayonnaise had the longest LT, possibly because of less favorable growth conditions such as low pH and the presence of lysozyme and ovotransferrin, which inhibit bacterial growth. Thus, rice balls containing mayonnaise exhibited a lower risk of illness than those containing soy sauce or gochujang. Because the intake of rice balls by Koreans is very high, the risk of listeriosis from rice balls is high. Thus, consistently managing food handling is necessary to control L. monocytogenes because it can grow even at low temperatures. These results may be useful in future studies evaluating the growth characteristics of L. monocytogenes in rice balls.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Supplementary Materials
Supplementary Table 1: metrics of fit for models of L. monocytogenes growth in rice balls, fitted to experimental data (n = 4). Supplementary Table 2: risk ranking of six rice balls stored at 7, 15, and 25°C for susceptible population.
GarbuttJ.19972ndLondon, UKEdwardArnoldJayJ. M.20036thGaithursburg, ML, USAAspen Publishers, Inc.MontvilleT. J.MatthewsK. R.20053rdWashington, DC, USAASM PressScallanE.HoekstraR. M.AnguloF. J.Foodborne illness acquired in the United States-major pathogens201117171510.3201/eid1701.p111012-s2.0-78650979311WHORisk assessment of listeria monocytogenes in ready-to-eat foods2004Geneva, SwitzerlandWHOhttp://www.fao.org/3/a-y5394e.pdf/McLauchlinJ.MitchellR. T.SmerdonW. J.JewellK.Listeria monocytogenes and listeriosis: a review of hazard characterisation for use in microbiological risk assessment of foods2004921153310.1016/s0168-1605(03)00326-x2-s2.0-1642345296EFSA (European Food Safety Authority)Panel on Biological Hazards. Listeria monocytogenes contamination of ready-to-eat foods and the risk for human health in the EU2018165134PaeM.Dietary habits and perception toward food additives according to the frequency of consumption of convenience food at convenience stores among university students in cheongju201621214015110.5720/kjcn.2016.21.2.140KimS.-J.BuS. Y.ChoiM.-K.Preference and the frequency of processed food intake according to the type of residence of college students in Korea201520318819610.5720/kjcn.2015.20.3.188Institute of Food Technologists2001Chicago, IL, USAInstitute of Food TechnologistsKorea Food and Drug Administration2018Seoul, South Koreahttp://www.foodsafetykorea.go.kr/portal/safefoodlife/food/foodRvlv/foodRvlv.do/BaeY. M.LeeS. Y.Effect of microwave treatment and packaging methods on extending the shelf-life of RTE rice balls at room temperature201026165170BahkG.-J.HongC.-H.OhD.-H.HaS.-D.ParkK.-H.ToddE. C. D.Modeling the level of contamination of Staphylococcus aureus in ready-to-eat kimbab in Korea20066961340134610.4315/0362-028x-69.6.13402-s2.0-33745085932LeeH.LeeG.YoonE.Computation of maximum edible time using monitoring data of Staphylococcus aureus in kimbap and food MicroModel2004194954MiyaS.TakahashiH.IshikawaT.FujiiT.KimuraB.Risk of Listeria monocytogenes contamination of raw ready-to-eat seafood products available at retail outlets in Japan201076103383338610.1128/aem.01456-092-s2.0-77952245904HongS.-H.ParkN.-Y.JoH.-J.Risk ranking determination of combination of foodborne pathogens and livestock or livestock products201530111210.13103/jfhs.2015.30.1.1HongC. H.SimW. C.ChunS. J.Predictive growth model of native isolated Listeria monocytogenes on raw pork as a function of temperature and time200537850855FakruddinM.MazumderR. M.MannanK. S. B.Predictive microbiology: modeling microbial responses in food201140121131VermeulenA.SmigicN.RajkovicA.Performance of a growth-no growth model for listeria monocytogenes developed for mayonnaise-based salads: influence of strain variability, food matrix, inoculation level, and presence of sorbic and benzoic acid20077092118212610.4315/0362-028x-70.9.21182-s2.0-34748820232CODEX 19981998Orlando, FL, USACodex Committee on Food Hygiene. Report of the thirty-first sessionChenY.DennisS. B.HartnettE.FDA-iRISK-A comparative risk assessment system for evaluating and ranking food-hazard pairs: case studies on microbial hazards201376337638510.4315/0362-028x.jfp-12-3722-s2.0-84875533730CapozziV.FioccoD.AmodioM. L.GalloneA.SpanoG.Bacterial stressors in minimally processed food20091073076310510.3390/ijms100730762-s2.0-68949098194LebertI.Robles-OlveraV.LebertA.Application of polynomial models to predict growth of mixed cultures of Pseudomonas spp. and Listeria in meat2000611273910.1016/s0168-1605(00)00359-72-s2.0-0034307806MoonH. J.LimJ. G.YoonK. S.Comparative study of change in Salmonella enteritidis and Salmonella typhimurium populations in egg white and yolk201631534234810.13103/jfhs.2016.31.5.342SmittleR. B.Microbiological safety of mayonnaise, salad dressings, and sauces produced in the United States: a review20006381144115310.4315/0362-028x-63.8.11442-s2.0-0033865519BatyF.Delignette-MullerM.-L.Estimating the bacterial lag time: which model, which precision?200491326127710.1016/j.ijfoodmicro.2003.07.0022-s2.0-1242293076YeK.WangH.ZhangX.JiangY.XuX.ZhouG.Development and validation of a molecular predictive model to describe the growth of listeria monocytogenes in vacuum-packaged chilled pork201332124625410.1016/j.foodcont.2012.11.0172-s2.0-84871803981PlaM.-L.OltraS.EstebanM.-D.AndreuS.PalopA.Comparison of primary models to predict microbial growth by the plate count and absorbance methods2015201511410.1155/2015/3650252-s2.0-84945925339BuchananR. L.WhitingR. C.DamertW. C.When is simple good enough: a comparison of the gompertz, baranyi, and three-phase linear models for fitting bacterial growth curves199714431332610.1006/fmic.1997.01252-s2.0-0030745758LiM. Y.SunX. M.ZhaoG. M.Comparison of mathematical models of lactic acid bacteria growth in vacuum-packaged raw beef stored at different temperatures2013784M600M60410.1111/j.1750-3841.2012.02995.x2-s2.0-84875938451LiM.NiuH.ZhaoG.Analysis of mathematical models of Pseudomonas spp. growth in pallet-package pork stored at different temperatures201393485586410.1016/j.meatsci.2012.11.0482-s2.0-84872336646ZulianiV.LebertI.AugustinJ.-C.GarryP.VendeuvreJ.-L.LebertA.Modelling the behaviour of Listeria monocytogenes in ground pork as a function of pH, water activity, nature and concentration of organic acid salts2007103353655010.1111/j.1365-2672.2007.03283.x2-s2.0-34547902847Sant’AnaA. S.LandgrafM.DestroM. T.FrancoB. D. G. M.Prevalence and counts of Salmonella spp. in minimally processed vegetables in Sao Paulo, Brazil20112812351237LuoK.HongS.-S.WangJ.Development of predictive models for the growth kinetics of Listeria monocytogenes on fresh pork under different storage temperatures201578592192610.4315/0362-028x.jfp-14-4282-s2.0-84929094367SorrellsK. M.EniglD. C.HatfieldJ. R.Effect of pH, acidulant, time, and temperature on the growth and survival of Listeria monocytogenes198952857157310.4315/0362-028x-52.8.571DevlieghereF.GeeraerdA. H.VersyckK. J.VandewaetereB.Van ImpeJ.DebevereJ.Growth of Listeria monocytogenes in modified atmosphere packed cooked meat products: a predictive model2001181536610.1006/fmic.2000.03782-s2.0-0035121942FranzE.TrompS. O.RijgersbergH.Van Der Fels-KlerxH. J.Quantitative microbial risk assessment for Escherichia coli O157:H7, Salmonella, and Listeria monocytogenes in leafy green vegetables consumed at salad bars201073227428510.4315/0362-028x-73.2.2742-s2.0-77649088900TienungoonS.RatkowskyD. A.McmeekinT. A.RossT.Growth limits of listeria monocytogenesas a function of temperature, pH, NaCl, and lactic acid200066114979498710.1128/aem.66.11.4979-4987.20002-s2.0-0033759097DingT.RahmanS. M. E.OhD.-H.Inhibitory effects of low concentration electrolyzed water and other sanitizers against foodborne pathogens on oyster mushroom201122231832210.1016/j.foodcont.2010.07.0302-s2.0-78649510830HuangL.Growth Kinetics of Escherichia coli O157:H7 in mechanically-tenderized beef20101401404810.1016/j.ijfoodmicro.2010.02.0132-s2.0-77952290749ChoK. M.KambirandaD. M.KimS. W.Simultaneous detection of food-borne pathogenic bacteria in ready-to-eat kimbab using multiplex PCR method20081712401245ShinG.-N.KimY.-R.KimM.-H.Nutritional evaluation of convenience meals in convenience stores near the universities201722537538610.5720/kjcn.2017.22.5.375PouillotR.HoelzerK.ChenY.DennisS. B.Listeria monocytogenesDose response revisited-incorporating adjustments for variability in strain virulence and host susceptibility20153519010810.1111/risa.122352-s2.0-84922630441NoordhoutC. M.DevleesschauwerB.AnguloF.The global burden of listeriosis: a systematic review and meta-analysis20141410731082Korea centers for disease control and prevention. Korea health statistics 2016: korea national health and nutrition examination survey (KNHANES VII-1). http://knhanes.cdc.go.kr/HavelaarA. H.HaagsmaJ. A.MangenM.-J. J.Disease burden of foodborne pathogens in The Netherlands, 20092012156323123810.1016/j.ijfoodmicro.2012.03.0292-s2.0-84861001014KemmerenJ. M.MangenM.-J. J.Van DuynhovenY. T. H. P.HavelaarA. H.Priority Setting Of Foodborne Pathogens-Disease Burden And Costs Of Selected Enteric Pathogens2006Veinna, AustriaRijksinstituut voor Volksgezondheid en Milieu (RIVM)NoordhoutC. M.DevleesschauwerB.AnguloF.The global burden of listeriosis: a systematic review and meta-analysis20141410731082LeeY. J.JungB. S.YoonH. J.KimK.-T.PaikH.-D.LeeJ.-Y.Predictive model for the growth kinetics of Listeria monocytogenes in raw pork meat as a function of temperature201444162110.1016/j.foodcont.2014.03.0242-s2.0-84925659208HongS.KimH.YoonK.Survival and risk comparison of Campylobacter jejuni on various processed meat products201613658010.3390/ijerph130605802-s2.0-84973526664