Source Attribution of Food-Borne Zoonoses in New Zealand: A Modified Hald Model

A Bayesian approach was developed by Hald et al .^{( 1 )} to estimate the contribution of different food sources to the burden of human salmonellosis in Denmark. This article describes the development of several modifications that can be used to adapt the model to different countries and pathogens. Our modified Hald model has several advantages over the original approach, which include the introduction of uncertainty in the estimates of source prevalence and an improved strategy for identifiability. We have applied our modified model to the two major food-borne zoonoses in New Zealand, namely, campylobacteriosis and salmonellosis. Major challenges were the data quality for salmonellosis and the inclusion of environmental sources of campylobacteriosis. We conclude that by modifying the Hald model we have improved its identifiability, made it more applicable to countries with less intensive surveillance, and feasible for other pathogens, in particular with respect to the inclusion of nonfood sources. The wider application and better understanding of this approach is of particular importance due to the value of the model for decision making and risk management.     

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.     

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.     

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.     

Quantitative Risk Assessment from Farm to Fork and Beyond: A Global Bayesian Approach Concerning Food-Borne Diseases

A novel approach to the quantitative assessment of food-borne risks is proposed. The basic idea is to use Bayesian techniques in two distinct steps: first by constructing a stochastic core model via a Bayesian network based on expert knowledge, and second, using the data available to improve this knowledge. Unlike the Monte Carlo simulation approach as commonly used in quantitative assessment of food-borne risks where data sets are used independently in each module, our consistent procedure incorporates information conveyed by data throughout the chain. It allows "back-calculation" in the food chain model, together with the use of data obtained "downstream" in the food chain. Moreover, the expert knowledge is introduced more simply and consistently than with classical statistical methods. Other advantages of this approach include the clear framework of an iterative learning process, considerable flexibility enabling the use of heterogeneous data, and a justified method to explore the effects of variability and uncertainty. As an illustration, we present an estimation of the probability of contracting a campylobacteriosis as a result of broiler contamination, from the standpoint of quantitative risk assessment. Although the model thus constructed is oversimplified, it clarifies the principles and properties of the method proposed, which demonstrates its ability to deal with quite complex situations and provides a useful basis for further discussions with different experts in the food chain.     

Modeling of interventions for reducing external Enterobacteriaceae contamination of broiler carcasses during processing

This article presents a mathematical model for the Enterobacteriaceae count on the surface of broiler chicken during slaughter and how it may be affected by different processing technologies. The model is based on a model originally developed for Campylobacter and has been adapted for Enterobacteriaceae using a Bayesian updating approach and hitherto unpublished data gathered from German abattoirs. The slaughter process in the model consists of five stages: input, scalding, defeathering, evisceration, washing, and chilling. The impact of various processing technologies along the broiler processing line on the Enterobacteriaceae count on the carcasses’ surface has been determined from literature data. The model is implemented in the software R and equipped with a graphical user interface which allows interactively to choose among different processing technologies for each stage along the processing line. Based on the choice of processing technologies the model estimates the Enterobacteriaceae count on the surface of each broiler chicken at each stage of processing. This result is then compared to a so-called baseline model which simulates a processing line with a fixed set of processing technologies. The model calculations showed how even very effective removal of bacteria on the exterior of the carcass in a previous step will be undone by the cross-contamination with leaked feces, if feces contain high concentrations of bacteria.     

A Practical Framework for the Construction of a Biotracing Model: Application to Salmonella in the Pork Slaughter Chain

A novel purpose of the use of mathematical models in quantitative microbial risk assessment (QMRA) is to identify the sources of microbial contamination in a food chain (i.e., biotracing). In this article we propose a framework for the construction of a biotracing model, eventually to be used in industrial food production chains where discrete numbers of products are processed that may be contaminated by a multitude of sources. The framework consists of steps in which a Monte Carlo model, simulating sequential events in the chain following a modular process risk modeling (MPRM) approach, is converted to a Bayesian belief network (BBN). The resulting model provides a probabilistic quantification of concentrations of a pathogen throughout a production chain. A BBN allows for updating the parameters of the model based on observational data, and global parameter sensitivity analysis is readily performed in a BBN. Moreover, a BBN enables “backward reasoning” when downstream data are available and is therefore a natural framework for answering biotracing questions. The proposed framework is illustrated with a biotracing model of Salmonella in the pork slaughter chain, based on a recently published Monte Carlo simulation model. This model, implemented as a BBN, describes the dynamics of Salmonella in a Dutch slaughterhouse and enables finding the source of contamination of specific carcasses at the end of the chain.     

A Bayesian Evidence Synthesis for Estimating Campylobacteriosis Prevalence

Stakeholders making decisions in public health and world trade need improved estimations of the burden‐of‐illness of foodborne infectious diseases. In this article, we propose a Bayesian meta‐analysis or more precisely a Bayesian evidence synthesis to assess the burden‐of‐illness of campylobacteriosis in France. Using this case study, we investigate campylobacteriosis prevalence, as well as the probabilities of different events that guide the disease pathway, by (i) employing a Bayesian approach on French and foreign human studies (from active surveillance systems, laboratory surveys, physician surveys, epidemiological surveys, and so on) through the chain of events that occur during an episode of illness and (ii) including expert knowledge about this chain of events. We split the target population using an exhaustive and exclusive partition based on health status and the level of disease investigation. We assume an approximate multinomial model over this population partition. Thereby, each observed data set related to the partition brings information on the parameters of the multinomial model, improving burden‐of‐illness parameter estimates that can be deduced from the parameters of the basic multinomial model. This multinomial model serves as a core model to perform a Bayesian evidence synthesis. Expert knowledge is introduced by way of pseudo‐data. The result is a global estimation of the burden‐of‐illness parameters with their accompanying uncertainty.     

Model Averaging in Microbial Risk Assessment Using Fractional Polynomials

The alleviation of food-borne diseases caused by microbial pathogen remains a great concern in order to ensure the well-being of the general public. The relation between the ingested dose of organisms and the associated infection risk can be studied using dose-response models. Traditionally, a model selected according to a goodness-of-fit criterion has been used for making inferences. In this article, we propose a modified set of fractional polynomials as competitive dose-response models in risk assessment. The article not only shows instances where it is not obvious to single out one best model but also illustrates that model averaging can best circumvent this dilemma. The set of candidate models is chosen based on biological plausibility and rationale and the risk at a dose common to all these models estimated using the selected models and by averaging over all models using Akaike's weights. In addition to including parameter estimation inaccuracy, like in the case of a single selected model, model averaging accounts for the uncertainty arising from other competitive models. This leads to a better and more honest estimation of standard errors and construction of confidence intervals for risk estimates. The approach is illustrated for risk estimation at low dose levels based on Salmonella typhi and Campylobacter jejuni data sets in humans. Simulation studies indicate that model averaging has reduced bias, better precision, and also attains coverage probabilities that are closer to the 95% nominal level compared to best-fitting models according to Akaike information criterion.     

Food Safety in the Domestic Environment: The Effect of Consumer Risk Information on Human Disease Risks

The improvement of food safety in the domestic environment requires a transdisciplinary approach, involving interaction between both the social and natural sciences. This approach is applied in a study on risks associated with Campylobacter on broiler meat. First, some web-based information interventions were designed and tested on participant motivation and intentions to cook more safely. Based on these self-reported measures, the intervention supported by the emotion "disgust" was selected as the most promising information intervention. Its effect on microbial cross-contamination was tested by recruiting a set of participants who prepared a salad with chicken breast fillet carrying a known amount of tracer bacteria. The amount of tracer that could be recovered from the salad revealed the transfer and survival of Campylobacter and was used as a measure of hygiene. This was introduced into an existing risk model on Campylobacter in the Netherlands to assess the effect of the information intervention both at the level of exposure and the level of human disease risk. We showed that the information intervention supported by the emotion "disgust" alone had no measurable effect on the health risk. However, when a behavioral cue was embedded within the instruction for the salad preparation, the risk decreased sharply. It is shown that a transdisciplinary approach, involving research on risk perception, microbiology, and risk assessment, is successful in evaluating the efficacy of an information intervention in terms of human health risks. The approach offers a novel tool for science-based risk management in the area of food safety.     

Uncertainty Evaluation of Human Risk Analysis (HRA) of Chemicals by Multiple Exposure Routes

The application of an ISO standard procedure (Guide to the Expression of Uncertainty in Measurement (GUM)) is here discussed to quantify uncertainty in human risk estimation under chronic exposure to hazardous chemical compounds. The procedure was previously applied to a simple model; in this article a much more complex model is used, i.e., multiple compound and multiple exposure pathways. Risk was evaluated using the usual methodologies: the deterministic reasonable maximum exposure (RME) and the statistical Monte Carlo method. In both cases, the procedures to evaluate uncertainty on risk values are detailed. Uncertainties were evaluated by different methodologies to account for the peculiarity of information about the single variable. The GUM procedure enables the ranking of variables by their contribution to uncertainty; it provides a criterion for choosing variables for deeper analysis. The obtained results show that the application of GUM procedure is easy and straightforward to quantify uncertainty and variability of risk estimation. Health risk estimation is based on literature data on a water table contaminated by three volatile organic compounds. Daily intake was considered by either ingestion of water or inhalation during showering. The results indicate one of the substances as the main contaminant, and give a criterion to identify the key component on which the treatment selection may be performed and the treatment process may be designed in order to reduce risk.     

Bayesian Temporal Source Attribution of Foodborne Zoonoses: Campylobacter in Finland and Norway

Statistical source attribution approaches of food‐related zoonoses can generally be based on reported diagnosed human cases and surveillance results from different food sources or reservoirs of bacteria. The attribution model, or probabilistic classifier, can thus be based on the (sub)typing information enabling comparison between human infections and samples derived from source surveillance. Having time series of both data allows analyzing temporal patterns over time providing a repeated natural experiment. A Bayesian approach combining both sources of information over a long time series is presented in the case of Campylobacter in Finland and Norway. The full model is transparently presented and derived from the Bayes theorem. Previous statistical source attribution approaches are here advanced (1) by explicit modeling of the cases not associated with any of the sources under surveillance over time, (2) by modeling uncertain prevalence in a food source by bacteria type over time, and (3) by implementing formal model fit assessment using posterior predictive discrepancy functions. Large proportion of all campylobacteriosis can be attributed to broiler, but considerable uncertainty remains over time. The source attribution is inherently incomplete if only the sources under surveillance are included in the model. All statistical source attribution approaches should include a model fit assessment for judgment of model performance with respect to relevant quantities of interest. It is especially relevant when the model aims at a synthesis of several incomplete information sources under significant uncertainty of explanatory variables.     

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.     

Fuzzy-Trace Theory, Risk Communication, and Product Labeling in Sexually Transmitted Diseases

Health care professionals are a major source of risk communications, but their estimation of risks may be compromised by systematic biases. We examined fuzzy-trace theory's predictions of professionals' biases in risk estimation for sexually transmitted infections (STIs) linked to: knowledge deficits (producing underestimation of STI risk, re-infection, and gender differences), gist-based mental representation of risk categories (producing overestimation of condom effectiveness for psychologically atypical but prevalent infections), retrieval failure for risk knowledge (producing greater risk underestimation when STIs are not specified), and processing interference involving combining risk estimates (producing biases in post-test estimation of infection, regardless of knowledge). One-hundred-seventy-four subjects (experts attending a national workshop, physicians, other health care professionals, and students) estimated the risk of teenagers contracting STIs, re-infection rates for males and females, and condom effectiveness in reducing infection risk. Retrieval was manipulated by asking estimation questions in two formats, a specific format that "unpacked" the STI category (infection types) and a global format that did not provide specific cues. Requesting estimates of infection risk after relevant knowledge was directly provided, isolating processing effects, assessed processing biases. As predicted, all groups of professionals underestimated the risk of STI transmission, re-infection, and gender differences, and overestimated the effectiveness of condoms, relative to published estimates. However, when questions provided better retrieval supports (specified format), estimation bias decreased. All groups of professionals also suffered from predicted processing biases. Although knowledge deficits contribute to estimation biases, the research showed that biases are also linked to fuzzy representations, retrieval failures, and processing errors Hence, interventions that are designed to improve risk perception among professionals must incorporate more than knowledge dissemination. They should also provide support for information representation, effective retrieval, and accurate processing.     

Importance of Uncertainty and Variability to Predicted Risks from Trophic Transfer of PCBs in Dredged Sediments

Biomagnification of organochlorine and other persistent organic contaminants by higher trophic level organisms represents one of the most significant sources of uncertainty and variability in evaluating potential risks associated with disposal of dredged materials. While it is important to distinguish between population variability (e.g., true population heterogeneity in fish weight, and lipid content) and uncertainty (e.g., measurement error), they can be operationally difficult to define separately in probabilistic estimates of human health and ecological risk. We propose a disaggregation of uncertain and variable parameters based on: (1) availability of supporting data; (2) the specific management and regulatory context (in this case, of the U.S. Army Corps of Engineers/U.S. Environmental Protection Agency tiered approach to dredged material management); and (3) professional judgment and experience in conducting probabilistic risk assessments. We describe and quantitatively evaluate several sources of uncertainty and variability in estimating risk to human health from trophic transfer of polychlorinated biphenyls (PCBs) using a case study of sediments obtained from the New York-New Jersey Harbor and being evaluated for disposal at an open water off-shore disposal site within the northeast region. The estimates of PCB concentrations in fish and dietary doses of PCBs to humans ingesting fish are expressed as distributions of values, of which the arithmetic mean or mode represents a particular fractile. The distribution of risk values is obtained using a food chain biomagnification model developed by Gobas by specifying distributions for input parameters disaggregated to represent either uncertainty or variability. Only those sources of uncertainty that could be quantified were included in the analysis. Results for several different two-dimensional Latin Hypercube analyses are provided to evaluate the influence of the uncertain versus variable disaggregation of model parameters. The analysis suggests that variability in human exposure parameters is greater than the uncertainty bounds on any particular fractile, given the described assumptions.     

Does Concern‐Driven Risk Management Provide a Viable Alternative to QRA?

This article discusses a concept of concern-driven risk management, in which qualitative expert judgments about whether concerns warrant specified risk management interventions are used in preference to quantitative risk assessment (QRA) to guide risk management decisions. Where QRA emphasizes formal quantitative assessment of the probable consequences caused by the recommended actions, and comparison to the probable consequences of alternatives, including the status quo, concern-driven risk management instead emphasizes perceived urgency or severity of the situation motivating recommended interventions. In many instances, especially those involving applications of the precautionary principle, no formal quantification or comparison of probable consequences for alternative decisions is seen as being necessary (or, perhaps, possible or desirable) prior to implementation of risk management measures. Such concern-driven risk management has been recommended by critics of QRA in several areas of applied risk management. Based on case studies and psychological literature on the empirical performance of judgment-based approaches to decision making under risk and uncertainty, we conclude that, although concern-driven risk management has several important potential political and psychological advantages over QRA, it is not clear that it performs better than (or as well as) QRA in identifying risk management interventions that successfully protect human health or achieve other desired consequences. Therefore, those who advocate replacing QRA with concern-driven alternatives, such as expert judgment and consensus decision processes, should assess whether their recommended alternatives truly outperform QRA, by the criterion of producing preferred consequences, before rejecting the QRA paradigm for practical applications.     

Bayesian Data Analysis of Severe Fatal Accident Risk in the Oil Chain

We analyze the risk of severe fatal accidents causing five or more fatalities and for nine different activities covering the entire oil chain. Included are exploration and extraction, transport by different modes, refining and final end use in power plants, heating or gas stations. The risks are quantified separately for OECD and non‐OECD countries and trends are calculated. Risk is analyzed by employing a Bayesian hierarchical model yielding analytical functions for both frequency (Poisson) and severity distributions (Generalized Pareto) as well as frequency trends. This approach addresses a key problem in risk estimation—namely the scarcity of data resulting in high uncertainties in particular for the risk of extreme events, where the risk is extrapolated beyond the historically most severe accidents. Bayesian data analysis allows the pooling of information from different data sets covering, for example, the different stages of the energy chains or different modes of transportation. In addition, it also inherently delivers a measure of uncertainty. This approach provides a framework, which comprehensively covers risk throughout the oil chain, allowing the allocation of risk in sustainability assessments. It also permits the progressive addition of new data to refine the risk estimates. Frequency, severity, and trends show substantial differences between the activities, emphasizing the need for detailed risk analysis.     

A Stochastic Assessment of the Public Health Risks of the Use of Macrolide Antibiotics in Food Animals

Campylobacteriosis is an important food-borne illness with more than a million U.S. cases annually. Antibiotic treatment is usually not required. However, erythromycin, a macrolide antibiotic, is recommended for the treatment of severe cases. Therefore, it is considered a critically important antibiotic and given special attention as to the risk that food animal use will lead to resistant infections and compromised human treatment. To assess this risk, we used a retrospective approach; estimating the number of campylobacteriosis cases caused by specific meat consumption utilizing the preventable fraction. We then determined the number of cases with macrolide resistance Campylobacter spp. based on a linear model relating the resistance fraction to on-farm macrolide use. In this article, we considered the uncertainties in the parameter estimates, utilized a more elaborate model of resistance development and separated C. coli and C. jejuni . There are no published data for the probability of compromised treatment outcome due to macrolide resistance. Therefore, our estimates of compromised treatment outcome were based on data for fluoroquinolone-resistant infections. The conservative results show the human health risks are extremely low. For example, the predicted risk of suboptimal human treatment of infection with C. coli from swine is only 1 in 82 million; with a 95% chance it could be as high as 1 in 49 million. Risks from C. jejuni in poultry or beef are even less. Reduced antibiotic use can adversely impact animal health. These low human risks should be weighed against the alternative risks.     

Use of Acute Toxicity to Estimate Carcinogenic Risk

Data on the effects of human exposure to carcinogens are limited, so that estimation of the risks of carcinogens must be obtained indirectly. Current risk estimates are generally based on lifetime animal bioassays which are expensive and which take more than two years to complete. We here show how data on acute toxicity can be used to make a preliminary estimate of carcinogenic risk and give an idea of the uncertainty in that risk estimate. The estimates obtained are biased upwards, and so are useful for setting interim standards and determining whether further study is worthwhile. A general scheme which incorporates the use of such estimates is outlined, and it is shown by example how adoption of the procedures suggested could have prevented regulatory hiatus in the past.     

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.