Introduction
Campylobacter jejuni is a microaerobic bacterium that can survive with the status of viable but non-culturable, and thus it is difficult to detect (Gomes et al., 2018; Hsieh et al., 2018). In Korea, there were no reported Campylobacter outbreaks until 2002, but the number of annual outbreaks has continued to increase since then (MFDS, 2018).
The symptoms of C. jejuni foodborne illness are fever, diarrhea, abdominal pain, and severe diseases such as Guillain–Barré syndrome and Reiter's syndrome (Lopes et al., 2018; Schmutz et al., 2017). C. jejuni is also the causative pathogen of traveler's diarrhea (Hsieh et al., 2018). To treat C. jejuni foodborne illness, antibiotics are used, but C. jejuni has become a multi-drug resistant bacterium that is resistant to several antibiotics (Lee et al., 2017; Zhong et al., 2016). Thus, the risk of C. jejuni should be assessed to manage and ensure the food safety.
Microbial risk assessment (MRA) is conducted to support scientific evidences to control food safety (Nauta et al., 2009), and it evaluates the probability for the risk of foodborne pathogens in food. MRA is comprised of four steps: hazard identification, exposure assessment, hazard characterization, and risk characterization (CAC, 2015). The MRA scheme is very important for the analysis of the risk. As each food has various distribution environments such as time and temperature, these should be reflected in the MRA scheme.
The ground meat products in Korea might be the hamburger patty, meatball, and cutlet. The sales of ground meat products ranked fourth among processed meat products in Korea in 2016 (MFDS, 2017). The possibility for C. jejuni contamination in ground meat products has been suggested. However, the risk of C. jejuni in ground meat products has not been evaluated. Therefore, the objective of this study was to evaluate the risk of foodborne illness caused by C. jejuni in ground meat products in Korea.
Materials and Methods
Prevalence and initial contamination of Campylobacter jejuni in ground meat products
To identify the prevalence of C. jejuni, three types of ground meat products (patty, meatball, and cutlet) were collected from supermarkets and foodservices in Korea, and 224 samples were analyzed in this study. They were made from only one meat (pork, chicken, or duck meat) or from more than two different meat (pork, beef, and chicken). Each sample was cut to obtain a 25-g sample, and the samples were placed in sterilized sample bags (3M™, St. Paul, MN, USA) and pummeled for 60 s with 50 mL of 0.1% buffered peptone water (BPW; Becton, Dickinson and Company, Sparks, MD, USA) using a BagMixer (Interscience, St. Nom, France). The 1-mL aliquots of the homogenates were plated on three modified CCDA-Preston agar plates (mCCDA; Oxoid Ltd., Basingstoke, Hampshire, UK), and the plates were then incubated microaerobically with CampyGen (Oxoid Ltd.) in a tightly–sealed container at 42℃ for 48 h. After incubation, one typical Campylobacter colony was streaked on two Colombia agar (bioMérieux, Marcy–l'Étoile, France) plates, and each plate was incubated microaerobically or aerobically at 42℃ for 48 h. Only one colony that grew on the microaerobically–incubated plate was further analyzed by colony PCR for the identification of C. jejuni, using specific primers [5′–CAA ATA AAG TTA GAG GTA GAA TGT–3′ (forward) and 5′–CCA TAA GCA CTA GCT AGC TGA T–3′ (reverse); Wang et al., 1992].
The initial contamination of C. jejuni in ground meat products was calculated using prevalence data, to estimate the probability of illness caused by C. jejuni intake. The C. jejuni prevalence data (PR) collected from this study were substituted into the Beta distribution (α, β), where α was the number of positive samples plus one, and β was one plus the number of positive samples, subtracted from the number of total samples (Vose, 1998). The equation [–Ln(1–PR)/weight], presented by Sanaa et al. (2004), was used to calculate the initial contamination of C. jejuni in ground meat products.
Development of a predictive model
Determination of a model food to develop the predictive model
C. jejuni strains ATCC 33560 and NCTC 11168 stored at –70℃ were streaked on Colombia agar and incubated microaerobically at 42℃ for 48 h. Each colony was streaked on fresh Colombia agar plates using a sterilized swab, and the plates were incubated microaerobically at 42℃ for 48 h. Five milliliters of phosphate-buffered saline (PBS, pH 7.4; 0.2 g of KH2PO4, 1.5 g of Na2HPO4, 8.0 g of NaCl, and 0.2 g of KCl in 1 L of distilled water) was poured over the colonies, and the colonies were harvested by scraping with a glass rod. The harvested cells were transferred to 50 mL conical tubes and centrifuged (1,912×g, 15 min, and 4℃). The pellets were washed twice with PBS and re-suspended in PBS. The OD values of the cell suspensions at 600 nm were adjusted to 2.0 to obtain 5–6 Log CFU/mL of C. jejuni. Two strains were mixed, and the mixture was then used as the inoculum.
Three ground meat products (patty, meatball, and cutlet) were purchased from markets and cut to obtain 25 g portions. Samples were transferred to sample bags, and 100 μL of the inocula were inoculated into each sample. The samples were rubbed 30 times to smear the inoculum, and they were stored at 4℃ or 10℃ for up to 5 d. The samples were analyzed microbiologically at appropriate time intervals. Into each sample bag, 50 mL of 0.1% BPW was poured and pummeled for 60 s using a BagMixer. After pummeling, the homogenates were diluted serially using 9 mL of 0.1% BPW, and 100 μL of the diluents were spread–plated on mCCDA. The mCCDA plates were incubated at 42℃ for 48 h under the microaerobic condition.
Development of a predictive model
To develop a predictive model to describe the behavior of C. jejuni, 100–μL aliquots of inocula were inoculated on sample surfaces in sample bags and rubbed 30 times. The inoculated samples were left for 15 min. The samples were then placed under an aerobic condition (same condition as a commercial product on the market) and stored at 4℃, 10℃, 15℃, 25℃, or 30℃ for up to 15 days. To enumerate the C. jejuni cell counts, samples were placed in sample bags, containing 50 mL of 0.1% BPW and pummeled for 60 s. The homogenates were diluted with 9 mL of 0.1% BPW, and 0.1 mL of the diluents were spread-plated on mCCDA with a flame-sterilized rod. The plates were microaerobically incubated at 42℃ for 48 h, and the colonies were then manually counted. C. jejuni cell counts were fitted by the Weibull model.
Log(N ) = Log(N0) — (time/δ )ρ
Where, N is the number of cells at a specific time, N0 is the initial number of cells, ρ is the curve shape, and δ is the time of the first decimal reduction. To analyze the effect of the storage temperature on kinetic parameters, the Davey model was fitted to the δ and ρ values.
Y = a + (b/T ) + (c/T2)
Where, Y is the δ or ρ value, a, b, and c are constants, and T is the temperature. For validation of the developed predictive model, additional experiments were performed at 12℃ and 23℃. The observed data for 12℃ and 23℃ were compared with the predicted data obtained from the developed predictive model. Differences between the observed and predicted data were quantified by the root mean square error (RMSE).
\(R M S E=\sqrt{\sum(\text { Observed value }-\text { Predicted value })^{2} / \mathrm{n}}\)
Temperature and time data
The temperatures and times for the storage, display, and transport of ground meat products were collected through communication with a person who works in the industry and from literatures (Jung, 2011; Kim, 2002). The collected temperature and time data for each step were analyzed with @RISK (Palisade Corporation, Ithaca, NY, USA) to obtain appropriate probabilistic distributions.
Consumption data for ground meat products
Consumption data of ground meat products were obtained through face-to-face interview by Gallup Korea (Seoul, Korea) from May 29 to July 8, 2015. The interviewee were about 1,500 adults (≥19 years old) in Korea (excluding Jeju Island). Of the 1,500 people, 743 (49.5%) were male and 757 (50.5%) were female, and they were classified into five groups (19–29 years, 30s, 40s, 50s, and ≥60s) according to their ages. In the questionnaire, the approximate size of the food was presented to help the respondents for estimating food amount they consumed. The sampling was extracted by square root proportional allocation and proportional quota sampling method to guarantee the representativeness and accuracy of the interview. More than 20% of the people were called and confirmed to verify if the interview was conducted properly. Collected data (95% confidence interval) were computerized by the SPSS program (IBM Software, North Carolina, USA) to obtain raw data. The computerized data were analyzed by @RISK to find an appropriate probabilistic distribution for consumption of ground meat products.
Dose-response model and risk characterization
To evaluate the dose-response of C. jejuni after intake of the pathogens, we searched a dose–response models. The MRA scheme was prepared, as shown in Fig. 1. Especially, a simulation model was prepared with the initial contamination level data for C. jejuni, predictive models, probabilistic distributions for the temperature and time, the probabilistic distribution for the consumption data, and a dose–response model in Microsoft Excel (version 2007, Microsoft Corporation, Seattle, WA, USA). A simulation was then performed by @RISK for 10,000 iterations to calculate the risk of C. jejuni foodborne illness from the intake of ground meat products.
Fig. 1. Scheme of Campylobacter jejuni risk assessment in ground meat products.
Results and Discussion
Prevalence and initial contamination level of Campylobacter jejuni
To investigate the prevalence of C. jejuni in ground meat products, the three ground meat products were selected, and 96 patties, 73 meatballs, and 55 cutlets were analyzed. Of the 224 ground meat products, there was no C. jejuni in any of the samples. C. jejuni is usually detected in fresh meats, such as chicken, duck, turkey, beef, and pork (Hussain et al., 2007; Korsak et al., 2015; Lee et al., 2017; Narvaez-Bravo et al., 2017; Zhao et al., 2001). Soultos et al. (2015) also showed that there was no C. jejuni in 105 Greek meat products, called as souvlaki. Additionally, meat products (pork chops, fermented sausage, and beef/chicken wieners) were negative for C. jejuni (Bohaychuk et al., 2006). These previous studies suggest that C. jejuni does not usually exist in processed meat products, which is consistent with our data. To preserve the food, some hurdles (oxygen, water activity, and other competing bacteria) are involved, and C. jejuni is sensitive to the hurdles (Duffy et al., 2001; Leistner, 2000; Soultos et al., 2015; Zhao et al., 2001). Thus, C. jejuni might not be detected in ground meat products. Although there were no C. jejuni–positive samples, the initial contamination level of C. jejuni in ground meat products was estimated with the Beta distribution suggested by Vose (1998), in which the α and β parameters were 1 and 225, respectively. The distribution was substituted in the equation presented by Sanaa et al. (2004), and then the mean initial contamination level of C. jejuni in ground meat products was –4.0±0.6 Log CFU/g (Fig. 2).
Fig. 2. Probability density for initial contamination levels of Campylobacter jejuni in ground meat product.
Predictive model
C. jejuni cells were inoculated in ground meat products (patty, meatball, and cutlet) and stored at 4℃ or 10℃ for 5 days. The C. jejuni cell counts for all of the samples decreased, but the C. jejuni cell count decrease was slower for the cutlet than for the patty and meatball (data not shown). Hence, the cutlet was selected as a model food, considering it as the worst case scenario. To develop a predictive model, C. jejuni–inoculated cutlet samples were stored at 4℃, 10℃, 15℃, 25℃, or 30℃, and the C. jejuni cell counts in the cutlets decreased as the time and temperature increased. The survival data were fitted by the Weibull model, and then, the kinetic parameters (δ and ρ) were calculated (Table 1). The mean of the δ value, derived from the developed model, was 94.2 h at 4℃, and it decreased as the temperature increased (Table 1). These results indicate that Campylobacter can survive longer at low temperatures than at high temperatures. The ρ values were 1.2 at 4℃, 10℃, and 15℃ (concave shape) and 0.5 at 25℃ and 30℃ (convex shape). The Davey model was fitted to the δ and ρ values to develop the secondary model. The secondary models are presented in Fig. 3. Since the R2 values were relatively high, the secondary models were suitable to describe the correlation between the temperature and the δ and ρ values. The validation result showed that the RMSE value was 0.589 Log CFU/g. This result indicates that the developed model was appropriate to describe the kinetic behavior of Campylobacter in ground meat products.
Table 1. δ and ρ values (mean±SD) calculated by the Weibull model for Campylobacter jejuni survival in cutlet as a model food for ground meat products during storage at 4℃, 10℃, 15℃, 25℃, and 30℃
1) Required time for the first decimal reduction.
2) Shape of curve.
Fig. 3. Secondary model of Campylobacter jejuni in cutlet as a function of storage temperature. Symbol, observed value; line, fitted line with the Davey model.
Temperature and time during distribution
Ground meat products are usually exposed to a freezing environment during storage and display at markets. Kim (2002) reported that the mean temperature of the shelf for frozen cutlets was –18.8℃ with –23.6℃ as the minimum temperature and –14℃ as the maximum temperature. Thus, the Pert distribution (–23.6, –18.8, –14) was applied for the storage and display temperature at the market to obtain the probabilistic distribution (Table 2). Additionally, Jung (2011) presented that the minimum, mean, and maximum temperatures during transport from the market to the home were 10℃, 18℃, and 25℃, respectively. Thus, the Pert distribution (10, 18, 25) was applied to obtain the probabilistic distribution (Table 2). The times for the storage and display of ground meat products were investigated by communication with merchants at markets. Thus, the Pert distribution (0, 72, 120) for market storage and the Pert distribution (0, 36, 168) for market display were applied to obtain the probabilistic distributions (Table 2). Jung (2011) showed that the minimum and maximum times for transport were 0.325 h and 1.643 h, respectively. The mean time for transport was calculated as 0.984 h from the minimum and maximum times. Thus, the Pert distribution (0.325, 0.984, 1.643) was used for the transport time to obtain the probabilistic distribution (Table 2).
Table 2. Simulation model and formulas in Microsoft Excel spreadsheet used to estimate the risk of Campylobacter jejuni in ground meat products with @RISK
1)Personal communication with merchants at market.
Consumption amount and frequency
The collected data from the face–to–face interviews were analyzed by @RISK, and the analysis suggested that the best distribution for the data was the Exponential distribution. The probabilistic distribution showed that the mean value for the consumption amount of ground meat products was 28.3 g per day (Fig. 4). The consumption frequency of ground meat products was 6.98% per day (Table 2).
Fig. 4. Probabilistic distribution for consumption amount of ground meat product fitted by @RISK.
Risk characterization
The dose–response model suggested by Medema et al. (1996), Nauta et al. (2007), and Teunis and Havelaar (2000) was chosen as appropriate models for C. jejuni. The probability of illness (risk) from C. jejuni intake was calculated by multiplying the probability of infection (Pinf) by the probability of illness given infection (Pill|inf) (Table 2). The simulation model, including the initial contamination level of C. jejuni, predictive models with probabilistic distributions for the temperature and time, the probabilistic distribution for the consumption data, and the dose-response model, is presented in Table 2. According to the result of the simulation, the C. jejuni cell counts decreased gradually with distribution from the market to the home. It means that the risk of C. jejuni gradually decreases during distribution. Finally, the probability of illness caused by C. jejuni intake from ground meat products at home was 5.68×10–10 per person per day, which can be considered a very low risk (Table 3).
Table 3. Probability of Campylobacter jejuni foodborne illness per person per day caused by ground meat products consumption
In conclusion, the developed quantitative MRA model was useful for estimating the risk of C. jejuni in ground meat products in Korea. Although the concern for C. jejuni foodborne illness from the intake of ground meat products has increased, the results from this study indicate that the risk of C. jejuni foodborne illness caused by the consumption of ground meat products is very low in Korea.
Conflict of Interest
No potential conflict of interest relevant to this article was reported.
Acknowledgements
This research was supported by a grant (14162MFDS077) from the Ministry of Food and Drug Safety in 2016.
Author Contributions
Conceptualization: Yoon Y. Data curation: Yoon Y, Yoon K-S, Seo K. Formal analysis: Lee J, Ha J. Methodology: Kim S, Choi Y, Oh H. Software: Lee H, Lee S. Investigation: Lee Y, Kim Y. Writing - original draft: Lee J. Writing - review & editing: Lee J, Lee H, Lee S, Kim S, Ha J, Choi Y, Oh H, Lee Y, Kim Y, Yoon K-S, Seo K, Yoon Y.
Ethics Approval
This article does not require IRB/IACUC approval because there are no human and animal participants.
참고문헌
- Bohaychuk VM, Gensler GE, King RK, Manninen KI, Sorensen O, Wu JT, Stiles ME, McMullen LM. 2006. Occurrence of pathogens in raw and ready-to-eat meat and poultry products collected from the retail marketplace in Edmonton, Alberta, Canada. J Food Prot 69:2176-2182. https://doi.org/10.4315/0362-028X-69.9.2176
- Codex Alimentarius Commission [CAC]. 2015. Procedural manual. 24th ed. World Health Organization/Food and Agricultural Organization of the United Nations, Rome, Italy.
- Davey KR. 1991. Applicability of the Davey (linear Arrhenius) predictive model to the lag phase of microbial growth. J Appl Bacteriol Banner 70:253-257. https://doi.org/10.1111/j.1365-2672.1991.tb02933.x
- Duffy EA, Belk KE, Sofos JN, Bellinger GR, Pape A, Smith GC. 2001. Extent of microbial contamination in United States pork retail products. J Food Prot 64:172-178. https://doi.org/10.4315/0362-028X-64.2.172
- Gomes CN, Passaglia J, Vilela FP, da Silva FMHSP, Duque SS, Falcao JP. 2018. High survival rates of Campylobacter coli under different stress conditions suggest that more rigorous food control measures might be needed in Brazil. Food Microbiol 73:327-333. https://doi.org/10.1016/j.fm.2018.02.014
- Havelaar AH, de Wit MAS, van Koningsveld R. 2000. Health burden in the Netherlands (1990-1995) due to infections with thermophilic Campylobacter species. National Institute of Public Health and the Environment [RIVM], Bilthoven, Netheland. Report No.: 284550004.
- Hsieh YH, Simpson S, Kerdahi K, Sulaiman IM. 2018. A comparative evaluation study of growth conditions for culturing the isolates of Campylobacter spp. Curr Microbiol 75:71-78. https://doi.org/10.1007/s00284-017-1351-6
- Hussain I, Mahmood MS, Akhtar M, Khan A. 2007. Prevalence of Campylobacter species in meat, milk and other food commodities in Pakistan. Food Microbiol 24:219-222. https://doi.org/10.1016/j.fm.2006.06.001
- Jung HK. 2011. Consumer survey and hazard analysis for the improvement of food hygiene and safety in purchase. M.S. Thesis, Korea Univ., Seoul, Korea.
- Kim JG. 2002. A survey of the temperature control of refrigerators and freezers in retail food shops. Korean J Environ Health Soc 28:161-171.
- Korsak D, Mackiw E, Rozynek E, Zylowska M. 2015. Prevalence of Campylobacter spp. In retail chicken, turkey, pork, and beef meat in Poland between 2009 and 2013. J Food Prot 78:1024-1028. https://doi.org/10.4315/0362-028X.JFP-14-353
- Lee J, Jeong J, Lee H, Ha J, Kim S, Choi Y, Oh H, Seo K, Yoon Y, Lee S. 2017. Antibiotic susceptibility, genetic diversity, and the presence of toxin producing genes in Campylobacter isolates from poultry. Int J Environ Res Public Health 14:E1400. https://doi.org/10.3390/ijerph14111400
- Leistner L. 2000. Basic aspects of food preservation by hurdle technology. Int J Food Microbiol 55:181-186. https://doi.org/10.1016/S0168-1605(00)00161-6
- Lopes GV, Landgraf M, Destro MT. 2018. Occurrence of Campylobacter in raw chicken and beef from retail outlets in Sao Paulo, Brazil. J Food Saf 38:e12442. https://doi.org/10.1111/jfs.12442
- Mafart P, Couvert O, Gaillard S, Leguerinel I. 2002. On calculating sterility in thermal preservation methods: Application of the Weibull frequency distribution model. Int J Food Microbiol 72:107-113. https://doi.org/10.1016/S0168-1605(01)00624-9
- Medema GJ, Teunis PFM, Havelaar AH, Haas CN. 1996. Assessment of the dose-response relationship of Campylobacter jejuni. Int J Food Microbiol 30:101-111. https://doi.org/10.1016/0168-1605(96)00994-4
- Ministry of Food and Drug Safety [MFDS]. 2017. Establishment of statistics for production reports of animal products in 2016. MDFS, Korea. p 24.
- Ministry of Food and Drug Safety [MFDS]. Food poisoning statistics. Available from: http://www.foodsafetykorea.go.kr/portal/healthyfoodlife/foodPoisoningStat.do?menu_no=519&menu_grp=MENU_GRP02. Accessed at Mar 12, 2018.
- Narvaez-Bravo C, Taboada EN, Mutschall SK, Aslam M. 2017. Epidemiology of antimicrobial resistant Campylobacter spp. isolated from retail meats in Canada. Int J Food Microbiol 253:43-47. https://doi.org/10.1016/j.ijfoodmicro.2017.04.019
- Nauta M, Hill A, Rosenquist H, Brynestad S, Fetsch A, van der Logt P, Fazil A, Christensen B, Katsma E, Borck B, Havelaar A. 2009. A comparison of risk assessments on Campylobacter in broiler meat. Int J Food Microbiol 129:107-123. https://doi.org/10.1016/j.ijfoodmicro.2008.12.001
- Nauta MJ, Jacobs-Reitsma WF, Havelaar AH. 2007. A risk assessment model for Campylobacter in broiler meat. Risk Anal 27:845-861. https://doi.org/10.1111/j.1539-6924.2006.00834.x
- Sanaa M, Coroller L, Cerf O. 2004. Risk assessment of listeriosis linked to the consumption of two soft cheeses made from raw milk: Camembert of Normandy and Brie of Meaux. Risk Anal 24:389-399. https://doi.org/10.1111/j.0272-4332.2004.00440.x
- Schmutz C, Mausezahl D, Bless PJ, Hatz C, Schwenkglenks M, Urbinello D. 2017. Estimating healthcare costs of acute gastroenteritis and human campylobacteriosis in Switzerland. Epidemiol Infect 145:627-641. https://doi.org/10.1017/S0950268816001618
- Soultos N, Iossifidou E, Psomas E, Tzikas Z, Lazou T. 2015. Incidence and antimicrobial susceptibility of Salmonella, Listeria, and Campylobacter spp. in raw "souvlaki" marketed in Thessaloniki (Northern Greece). J Hell Vet Med Soc 66:161-168.
- Teunis PF, Havelaar AH. 2000. The beta Poisson dose-response model is not a single-hit model. Risk Anal 20:513-520. https://doi.org/10.1111/0272-4332.204048
- Vose DJ. 1998. The application of quantitative risk assessment to microbial food safety. J Food Prot 61:640-648. https://doi.org/10.4315/0362-028X-61.5.640
- Wang RF, Salvic MF, CaoWW. 1992. A rapid PCR method for direct detection of low numbers of Campylobacter jejuni. J Rapid Methods Autom Microbiol 1:101-108. https://doi.org/10.1111/j.1745-4581.1992.tb00074.x
- Zhao C, Ge B, De Villena J, Sudler R, Yeh E, Zhao S, White DG, Wagner D, Meng J. 2001. Prevalence of Campylobacter spp., Escherichia coli, and Salmonella serovars in retail chicken, turkey, pork, and beef from the Greater Washington, DC, area. Appl Environ Microbiol 67:5431-5436. https://doi.org/10.1128/AEM.67.12.5431-5436.2001
- Zhong X, Wu Q, Zhang J, Shen S. 2016. Prevalence, genetic diversity and antimicrobial susceptibility of Campylobacter jejuni isolated from retail food in China. Food Control 62:10-15. https://doi.org/10.1016/j.foodcont.2015.09.032
피인용 문헌
- 우렁쉥이에 대한 병원성 비브리오균 정량적 미생물 위해평가 vol.35, pp.1, 2020, https://doi.org/10.13103/jfhs.2020.35.1.51
- Antimicrobial Effects of EcoCal® and GF Bactostop® Formulated in Emulsified Sausages against Lactic Acid Bacteria vol.35, pp.3, 2019, https://doi.org/10.13103/jfhs.2020.35.3.279
- Quantitative microbial risk assessment for Clostridium perfringens foodborne illness following consumption of kimchi in South Korea vol.29, pp.8, 2020, https://doi.org/10.1007/s10068-020-00754-2
- Quantitative microbial risk assessment of Vibrio parahaemolyticus foodborne illness of sea squirt (Halocynthia roretzi) in South Korea vol.24, pp.2, 2021, https://doi.org/10.47853/fas.2021.e8
- Contamination of Clostridium perfringens in soy sauce, and quantitative microbial risk assessment for C. perfringens through soy sauce consumption vol.9, pp.4, 2019, https://doi.org/10.1002/fsn3.2182
- Growth of Salmonella in napa cabbage kimchi during fermentation vol.28, pp.4, 2019, https://doi.org/10.11002/kjfp.2021.28.4.532