Effectiveness and Efficiency of Controlling Campylobacter on Broiler Chicken Meat

Campylobacter bacteria are an important cause of foodborne infections. We estimated the potential costs and benefits of a large number of possible interventions to decrease human exposure to Campylobacter by consumption of chicken meat, which accounts for 20-40% of all cases of human campylobacteriosis in the Netherlands. For this purpose, a farm-to-fork risk assessment model was combined with economic analysis and epidemiological data. Reduction of contamination at broiler farms could be efficient in theory. However, it is unclear which hygienic measures need to be taken and the costs can be very high. The experimental treatment of colonized broiler flocks with bacteriophages has proven to be effective and could also be cost efficient, if confirmed in practice. Since a major decrease of infections at the broiler farm is not expected in the short term, additional measures in the processing plant were also considered. At this moment, guaranteed Campylobacter-free chicken meat at the retail level is not realistic. The most promising interventions in the processing plant are limiting fecal leakage during processing and separation of contaminated and noncontaminated flocks (scheduling), followed by decontamination of the contaminated flock. New (faster and more sensitive) test methods to detect Campylobacter colonization in broilers flocks are a prerequisite for successful scheduling scenarios. Other methods to decrease the contamination of meat of colonized flocks such as freezing and heat treatment are more expensive and/or less effective than chemical decontamination.     

Development of a Quantitative Microbial Risk Assessment for Human Salmonellosis Through Household Consumption of Fresh Minced Pork Meat in Belgium

A quantitative microbial risk assessment (QMRA) according to the Codex Alimentarius Principles is conducted to evaluate the risk of human salmonellosis through household consumption of fresh minced pork meat in Belgium. The quantitative exposure assessment is carried out by building a modular risk model, called the METZOON-model, which covers the pork production from farm to fork. In the METZOON-model, the food production pathway is split up in six consecutive modules: (1) primary production, (2) transport and lairage, (3) slaughterhouse, (4) postprocessing, (5) distribution and storage, and (6) preparation and consumption. All the modules are developed to resemble as closely as possible the Belgian situation, making use of the available national data. Several statistical refinements and improved modeling techniques are proposed. The model produces highly realistic results. The baseline predicted number of annual salmonellosis cases is 20,513 ( SD 9061.45). The risk is estimated higher for the susceptible population (estimate  4.713 × 10⁻⁵; SD 1.466 × 10⁻⁵  ) compared to the normal population  (estimate 7.704 × 10⁻⁶; SD 5.414 × 10⁻⁶)  and is mainly due to undercooking and to a smaller extent to cross-contamination in the kitchen via cook's hands.     

A Poultry-Processing Model for Quantitative Microbiological Risk Assessment

A poultry-processing model for a quantitative microbiological risk assessment (QMRA) of campylobacter is presented, which can also be applied to other QMRAs involving poultry processing. The same basic model is applied in each consecutive stage of industrial processing. It describes the effects of inactivation and removal of the bacteria, and the dynamics of cross-contamination in terms of the transfer of campylobacter from the intestines to the carcass surface and the environment, from the carcasses to the environment, and from the environment to the carcasses. From the model it can be derived that, in general, the effect of inactivation and removal is dominant for those carcasses with high initial bacterial loads, and cross-contamination is dominant for those with low initial levels. In other QMRA poultry-processing models, the input-output relationship between the numbers of bacteria on the carcasses is usually assumed to be linear on a logarithmic scale. By including some basic mechanistics, it is shown that this may not be realistic. As nonlinear behavior may affect the predicted effects of risk mitigations; this finding is relevant for risk management. Good knowledge of the variability of bacterial loads on poultry entering the process is important. The common practice in microbiology to only present geometric mean of bacterial counts is insufficient: arithmetic mean are more suitable, in particular, to describe the effect of cross-contamination. The effects of logistic slaughter (scheduled processing) as a risk mitigation strategy are predicted to be small. Some additional complications in applying microbiological data obtained in processing plants are discussed.     

Cross-Contamination During Food Preparation: A Mechanistic Model Applied to Chicken-Borne Campylobacter

We develop a model for bacterial cross-contamination during food preparation in the domestic kitchen and apply this to the case of Campylobacter-contaminated chicken breast. Building blocks of the model are the routines performed during food preparation, with their associated probabilities of bacterial transfer between food items and kitchen utensils. The model is used in a quantitative microbiological risk assessment (QMRA) of Campylobacter in the Netherlands. Using parameter values from the literature and performing elementary sensitivity analyses, we show that cross-contamination can contribute significantly to the risk of Campylobacter infection and find that cleaning frequency of kitchen utensils and thoroughness of rinsing of raw food items after preparation has more impact on cross-contamination than previously emphasized. Furthermore, we argue that especially more behavioral data on hygiene during food preparation is needed for a comprehensive Campylobacter risk assessment.     

The BSE Risk of Processing Meat and Bone Meal in Nonruminant Feed: A Quantitative Assessment for the Netherlands

The total ban on use of meat and bone meal (MBM) in livestock feed has been very successful in reducing bovine spongiform encephalopathy (BSE) spread, but also implies a waste of high-quality proteins resulting in economic and ecological loss. Now that the BSE epidemic is fading out, a partial lifting of the MBM ban might be considered. The objective of this study was to assess the BSE risk for the Netherlands if MBM derived from animals fit for human consumption, i.e., category 3 MBM, would be used in nonruminant feed. A stochastic simulation model was constructed that calculates (1) the probability that infectivity of undetected BSE-infected cows ends up with calves and (2) the quantity of infectivity ( Q_inf ) consumed by calves in case of such an incident. Three pathways were considered via which infectivity can reach cattle: (1) cross-contamination in the feed mill, (2) cross-contamination on the primary farm, and (3) pasture contamination. Model calculations indicate that the overall probability that infectivity ends up with calves is 3.2%. In most such incidents the Q_inf is extremely small (median = 6.5 × 10⁻¹² ID₅₀; mean = 1.8 × 10⁻⁴ ID₅₀), corresponding to an average probability of 1.3 × 10⁻⁴ that an incident results in ≥1 new BSE infections. Cross-contamination in the feed mill is the most risky pathway. Combining model results with Dutch BSE prevalence estimates for the coming years, it can be concluded that the BSE risk of using category 3 MBM derived from Dutch cattle in nonruminant feed is very low.     

A Comparative Exposure Assessment of Campylobacter in Ontario,Canada

To inform source attribution efforts, a comparative exposure assessment was developed to estimate the relative exposure to Campylobacter, the leading bacterial gastrointestinal disease in Canada, for 13 different transmission routes within Ontario, Canada, during the summer. Exposure was quantified with stochastic models at the population level, which incorporated measures of frequency, quantity ingested, prevalence, and concentration, using data from FoodNet Canada surveillance, the peer‐reviewed and gray literature, other Ontario data, and data that were specifically collected for this study. Models were run with @Risk software using Monte Carlo simulations. The mean number of cells of Campylobacter ingested per Ontarian per day during the summer, ranked from highest to lowest is as follows: household pets, chicken, living on a farm, raw milk, visiting a farm, recreational water, beef, drinking water, pork, vegetables, seafood, petting zoos, and fruits. The study results identify knowledge gaps for some transmission routes, and indicate that some transmission routes for Campylobacter are underestimated in the current literature, such as household pets and raw milk. Many data gaps were identified for future data collection consideration, especially for the concentration of Campylobacter in all transmission routes.     

Attribution of Human VTEC O157 Infection from Meat Products: A Quantitative Risk Assessment Approach

To address the risk posed to human health by the consumption of VTEC O157 within contaminated pork, lamb, and beef products within Great Britain, a quantitative risk assessment model has been developed. This model aims to simulate the prevalence and amount of VTEC O157 in different meat products at consumption within a single model framework by adapting previously developed models. The model is stochastic in nature, enabling both variability (natural variation between animals, carcasses, products) and uncertainty (lack of knowledge) about the input parameters to be modeled. Based on the model assumptions and data, it is concluded that the prevalence of VTEC O157 in meat products (joints and mince) at consumption is low (i.e., <0.04%). Beef products, particularly beef burgers, present the highest estimated risk with an estimated eight out of 100,000 servings on average resulting in human infection with VTEC O157.     

A Quantitative Microbiological Risk Assessment for Campylobacter in Petting Zoos

The significance of petting zoos for transmission of Campylobacter to humans and the effect of interventions were estimated. A stochastic QMRA model simulating a child or adult visiting a Dutch petting zoo was built. The model describes the transmission of Campylobacter in animal feces from the various animal species, fences, and the playground to ingestion by visitors through touching these so‐called carriers and subsequently touching their lips. Extensive field and laboratory research was done to fulfill data needs. Fecal contamination on all carriers was measured by swabbing in 10 petting zoos, using Escherichia coli as an indicator. Carrier‐hand and hand‐lip touching frequencies were estimated by, in total, 13 days of observations of visitors by two observers at two petting zoos. The transmission from carrier to hand and from hand to lip by touching was measured using preapplied cow feces to which E. coli WG5 was added as an indicator. Via a Beta‐Poisson dose‐response function, the number of Campylobacter cases for the whole of the Netherlands (16 million population) in a year was estimated at 187 and 52 for children and adults, respectively, so 239 in total. This is significantly lower than previous QMRA results on chicken fillet and drinking water consumption. Scenarios of 90% reduction of the contamination (meant to mimic cleaning) of all fences and just goat fences reduces the number of cases by 82% and 75%, respectively. The model can easily be adapted for other fecally transmitted pathogens.     

Risk Assessment of Human Toxoplasmosis Associated with the Consumption of Pork Meat in Italy

Toxoplasmosis is a cosmopolitan disease and has a broad range of hosts, including humans and several wild and domestic animals. The human infection is mostly acquired through the consumption of contaminated food and pork meat has been recognized as one of the major sources of transmission. There are, however, certain fundamental differences between countries; therefore, the present study specifically aims to evaluate the exposure of the Italian population to Toxoplasma gondii through the ingestion of several types of pork meat products habitually consumed in Italy and to estimate the annual number of human infections within two subgroups of the population. A quantitative risk assessment model was built for this reason and was enriched with new elements in comparison to other similar risk assessments in order to enhance its accuracy. Sensitivity analysis and two alternative scenarios were implemented to identify the factors that have the highest impact on risk and to simulate different plausible conditions, respectively. The estimated overall average number of new infections per year among adults is 12,513 and 92 for pregnant women. The baseline model showed that almost all these infections are associated with the consumption of fresh meat cuts and preparations (mean risk of infection varied between 4.5 × 10⁻⁵ and 5.5 × 10⁻⁵) and only a small percentage is due to fermented sausages/salami. On the contrary, salt‐cured meat products seem to pose minor risk but further investigations are needed to clarify still unclear aspects. Among all the considered variables, cooking temperature and bradyzoites’ concentration in muscle impacted most the risk.     

A Linear Model for Managing the Risk of Antimicrobial Resistance Originating in Food Animals

A linear population risk model used by the U.S. Food and Drug Administration (FDA) Center for Veterinary Medicine (CVM) estimates the risk of human cases of campylobacteriosis caused by fluoroquinolone-resistant Campylobacter. Among the cases of campylobacteriosis attributed to domestically produced chicken, the fluoroquinolone resistance is assumed to result from the use of fluoroquinolones in poultry in the United States. Properties of the linear population risk model are contrasted with those of a farm-to-fork model commonly used for microbial risk assessments. The utility of the linear population model for the purpose for which it was used by CVM is discussed.     

An Ecological Model for Quantitative Risk Assessment for Schistosomiasis: The Case of a Patchy Environment in the Coastal Tropical Area of Northeastern Brazil

We developed a stochastic model for quantitative risk assessment for the Schistosoma mansoni (SM) parasite, which causes an endemic disease of public concern. The model provides answers in a useful format for public health decisions, uses data and expert opinion, and can be applied to any landscape where the snail Biomphalaria glabrata is the main intermediate host (South and Central America, the Caribbean, and Africa). It incorporates several realistic and case‐specific features: stage‐structured parasite populations, periodic praziquantel (PZQ) drug treatment for humans, density dependence, extreme events (prolonged rainfall), site‐specific sanitation quality, environmental stochasticity, monthly rainfall variation, uncertainty in parameters, and spatial dynamics. We parameterize the model through a real‐world application in the district of Porto de Galinhas (PG), one of the main touristic destinations in Brazil, where previous studies identified four parasite populations within the metapopulation. The results provide a good approximation of the dynamics of the system and are in agreement with our field observations, i.e., the lack of basic infrastructure (sanitation level and health programs) makes PG a suitable habitat for the persistence and growth of a parasite metapopulation. We quantify the risk of SM metapopulation explosion and quasi‐extinction and the time to metapopulation explosion and quasi‐extinction. We evaluate the sensitivity of the results under varying scenarios of future periodic PZQ treatment (based on the Brazilian Ministry of Health's plan) and sanitation quality. We conclude that the plan might be useful to slow SM metapopulation growth but not to control it. Additional investments in better sanitation are necessary.     

Risk Assessment of Human Listeriosis from Semisoft Cheeses Made from Raw Sheep's Milk in Lazio and Tuscany (Italy)

Semisoft cheese made from raw sheep's milk is traditionally and economically important in southern Europe. However, raw milk cheese is also a known vehicle of human listeriosis and contamination of sheep cheese with Listeria monocytogenes has been reported. In the present study, we have developed and applied a quantitative risk assessment model, based on available evidence and challenge testing, to estimate risk of invasive listeriosis due to consumption of an artisanal sheep cheese made with raw milk collected from a single flock in central Italy. In the model, contamination of milk may originate from the farm environment or from mastitic animals, with potential growth of the pathogen in bulk milk and during cheese ripening. Based on the 48‐day challenge test of a local semisoft raw sheep's milk cheese we found limited growth only during the initial phase of ripening (24 hours) and no growth or limited decline during the following ripening period. In our simulation, in the baseline scenario, 2.2% of cheese servings are estimated to have at least 1 colony forming unit (CFU) per gram. Of these, 15.1% would be above the current E.U. limit of 100 CFU/g (5.2% would exceed 1,000 CFU/g). Risk of invasive listeriosis per random serving is estimated in the 10^?12 range (mean) for healthy adults, and in the 10^?10 range (mean) for vulnerable populations. When small flocks (10–36 animals) are combined with the presence of a sheep with undetected subclinical mastitis, risk of listeriosis increases and such flocks may represent a public health risk.     

Risk Ranking of Antimicrobial‐Resistant Hazards Found in Meat in Switzerland

Human exposure to bacteria resistant to antimicrobials and transfer of related genes is a complex issue and occurs, among other pathways, via meat consumption. In a context of limited resources, the prioritization of risk management activities is essential. Since the antimicrobial resistance (AMR) situation differs substantially between countries, prioritization should be country specific. The objective of this study was to develop a systematic and transparent framework to rank combinations of bacteria species resistant to selected antimicrobial classes found in meat, based on the risk they represent for public health in Switzerland. A risk assessment model from slaughter to consumption was developed following the Codex Alimentarius guidelines for risk analysis of foodborne AMR. Using data from the Swiss AMR monitoring program, 208 combinations of animal species/bacteria/antimicrobial classes were identified as relevant hazards. Exposure assessment and hazard characterization scores were developed and combined using multicriteria decision analysis. The effect of changing weights of scores was explored with sensitivity analysis. Attributing equal weights to each score, poultry‐associated combinations represented the highest risk. In particular, contamination with extended‐spectrum β‐lactamase/plasmidic AmpC‐producing Escherichia coli in poultry meat ranked high for both exposure and hazard characterization. Tetracycline‐ or macrolide‐resistant Enterococcus spp., as well as fluoroquinolone‐ or macrolide‐resistant Campylobacter jejuni , ranked among combinations with the highest risk. This study provides a basis for prioritizing future activities to mitigate the risk associated with foodborne AMR in Switzerland. A user‐friendly version of the model was provided to risk managers; it can easily be adjusted to the constantly evolving knowledge on AMR.     

Ecological Risk Assessment Conceptual Model Formulation for Nonindigenous Species

This article addresses the application of ecological risk assessment at the regional scale to the prediction of impacts due to invasive or nonindigenous species (NIS). The first section describes risk assessment, the decision-making process, and introduces regional risk assessment. A general conceptual model for the risk assessment of NIS is then presented based upon the regional risk assessment approach. Two diverse examples of the application of this approach are presented. The first example is based upon the dynamics of introduced plasmids into bacteria populations. The second example is the application risk assessment approach to the invasion of a coastal marine site of Cherry Point, Washington, USA by the European green crab. The lessons learned from the two examples demonstrate that assessment of the risks of invasion of NIS will have to incorporate not only the characteristics of the invasive species, but also the other stresses and impacts affecting the region of interest.     

A QMRA for the Transmission of ESBL‐Producing Escherichia coli and Campylobacter from Poultry Farms to Humans Through Flies

The public health significance of transmission of ESBL‐producing Escherichia coli and Campylobacter from poultry farms to humans through flies was investigated using a worst‐case risk model. Human exposure was modeled by the fraction of contaminated flies, the number of specific bacteria per fly, the number of flies leaving the poultry farm, and the number of positive poultry houses in the Netherlands. Simplified risk calculations for transmission through consumption of chicken fillet were used for comparison, in terms of the number of human exposures, the total human exposure, and, for Campylobacter only, the number of human cases of illness. Comparing estimates of the worst‐case risk of transmission through flies with estimates of the real risk of chicken fillet consumption, the number of human exposures to ESBL‐producing E. coli was higher for chicken fillet as compared with flies, but the total level of exposure was higher for flies. For Campylobacter, risk values were nearly consistently higher for transmission through flies than for chicken fillet consumption. This indicates that the public health risk of transmission of both ESBL‐producing E. coli and Campylobacter to humans through flies might be of importance. It justifies further modeling of transmission through flies for which additional data (fly emigration, human exposure) are required. Similar analyses of other environmental transmission routes from poultry farms are suggested to precede further investigations into flies.     

A Monte Carlo Risk Assessment Model for Acrylamide Formation in French Fries

The objective of this study is to estimate the likely human exposure to the group 2a carcinogen, acrylamide, from French fries by Irish consumers by developing a quantitative risk assessment model using Monte Carlo simulation techniques. Various stages in the French-fry-making process were modeled from initial potato harvest, storage, and processing procedures. The model was developed in Microsoft Excel with the @Risk add-on package. The model was run for 10,000 iterations using Latin hypercube sampling. The simulated mean acrylamide level in French fries was calculated to be 317 μg/kg. It was found that females are exposed to smaller levels of acrylamide than males (mean exposure of 0.20 μg/kg bw/day and 0.27 μg/kg bw/day, respectively). Although the carcinogenic potency of acrylamide is not well known, the simulated probability of exceeding the average chronic human dietary intake of 1 μg/kg bw/day (as suggested by WHO) was 0.054 and 0.029 for males and females, respectively. A sensitivity analysis highlighted the importance of the selection of appropriate cultivars with known low reducing sugar levels for French fry production. Strict control of cooking conditions (correlation coefficient of 0.42 and 0.35 for frying time and temperature, respectively) and blanching procedures (correlation coefficient −0.25) were also found to be important in ensuring minimal acrylamide formation.     

NUSAP Method for Evaluating the Data Quality in a Quantitative Microbial Risk Assessment Model for Salmonella in the Pork Production Chain

The numeral unit spread assessment pedigree (NUSAP) system was implemented to evaluate the quality of input parameters in a quantitative microbial risk assessment (QMRA) model for Salmonella spp. in minced pork meat. The input parameters were grouped according to four successive exposure pathways: (1) primary production (2) transport, holding, and slaughterhouse, (3) postprocessing, distribution, and storage, and (4) preparation and consumption. An inventory of 101 potential input parameters was used for building the QMRA model. The characteristics of each parameter were defined using a standardized procedure to assess (1) the source of information, (2) the sampling methodology and sample size, and (3) the distributional properties of the estimate. Each parameter was scored by a panel of experts using a pedigree matrix containing four criteria (proxy, empirical basis, method, and validation) to assess the quality, and this was graphically represented by means of kite diagrams. The parameters obtained significantly lower scores for the validation criterion as compared with the other criteria. Overall strengths of parameters related to the primary production module were significantly stronger compared to the other modules (the transport, holding, and slaughterhouse module, the processing, distribution, and storage module, and the preparation and consumption module). The pedigree assessment contributed to select 20 parameters, which were subsequently introduced in the QMRA model. The NUSAP methodology and kite diagrams are objective tools to discuss and visualize the quality of the parameters in a structured way. These two tools can be used in the selection procedure of input parameters for a QMRA, and can lead to a more transparent quality assurance in the QMRA.     

Quantitative Risk Assessment of Vibrio parahaemolyticus in Finfish: A Model of Raw Horse Mackerel Consumption in Japan

The aim of this study was to evaluate the effects of implemented control measures to reduce illness induced by Vibrio parahaemolyticus (V. parahaemolyticus) in horse mackerel (Trachurus japonicus), seafood that is commonly consumed raw in Japan. On the basis of currently available experimental and survey data, we constructed a quantitative risk model of V. parahaemolyticus in horse mackerel from harvest to consumption. In particular, the following factors were evaluated: bacterial growth at all stages, effects of washing the fish body and storage water, and bacterial transfer from the fish surface, gills, and intestine to fillets during preparation. New parameters of the beta‐Poisson dose‐response model were determined from all human feeding trials, some of which have been used for risk assessment by the U.S. Food and Drug Administration (USFDA). The probability of illness caused by V. parahaemolyticus was estimated using both the USFDA dose‐response parameters and our parameters for each selected pathway of scenario alternatives: washing whole fish at landing, storage in contaminated water, high temperature during transportation, and washing fish during preparation. The last scenario (washing fish during preparation) was the most effective for reducing the risk of illness by about a factor of 10 compared to no washing at this stage. Risk of illness increased by 50% by exposure to increased temperature during transportation, according to our assumptions of duration and temperature. The other two scenarios did not significantly affect risk. The choice of dose‐response parameters was not critical for evaluation of control measures.     

Campylobacter QMRA: A Bayesian Estimation of Prevalence and Concentration in Retail Foods Under Clustering and Heavy Censoring

A Bayesian statistical temporal‐prevalence‐concentration model (TPCM) was built to assess the prevalence and concentration of pathogenic campylobacter species in batches of fresh chicken and turkey meat at retail. The data set was collected from Finnish grocery stores in all the seasons of the year. Observations at low concentration levels are often censored due to the limit of determination of the microbiological methods. This model utilized the potential of Bayesian methods to borrow strength from related samples in order to perform under heavy censoring. In this extreme case the majority of the observed batch‐specific concentrations was below the limit of determination. The hierarchical structure was included in the model in order to take into account the within‐batch and between‐batch variability, which may have a significant impact on the sample outcome depending on the sampling plan. Temporal changes in the prevalence of campylobacter were modeled using a Markovian time series. The proposed model is adaptable for other pathogens if the same type of data set is available. The computation of the model was performed using OpenBUGS software.     

Quantitative Risk Assessment Relating to Adventitious Presence of Allergens in Food: A Probabilistic Model Applied to Peanut in Chocolate

Peanut allergy is a public health concern, owing to the high prevalence in France and the severity of the reactions. Despite peanut-containing product avoidance diets, a risk may exist due to the adventitious presence of peanut allergens in a wide range of food products. Peanut is not mentioned in their ingredients list, but precautionary labeling is often present. A method of quantifying the risk of allergic reactions following the consumption of such products is developed, taking the example of peanut in chocolate tablets. The occurrence of adventitious peanut proteins in chocolate and the dose-response relationship are estimated with a Bayesian approach using available published data. The consumption pattern is described by the French individual consumption survey INCA2. Risk simulations are performed using second-order Monte Carlo simulations, which separately propagates variability and uncertainty of the model input variables. Peanut allergens occur in approximately 36% of the chocolates, leading to a mean exposure level of 0.2 mg of peanut proteins per eating occasion. The estimated risk of reaction averages 0.57% per eating occasion for peanut-allergic adults. The 95% values of the risk stand between 0 and 3.61%, which illustrates the risk variability. The uncertainty, represented by the 95% credible intervals, is concentrated around these risk estimates. Children have similar results. The conclusion is that adventitious peanut allergens induce a risk of reaction for a part of the French peanut-allergic population. The method developed can be generalized to assess the risk due to the consumption of every foodstuff potentially contaminated by allergens.