Estimating Listeria monocytogenes Growth in Ready-to-Eat Chicken Salad Using a Challenge Test for Quantitative Microbial Risk Assessment

Currently, there is a growing preference for convenience food products, such as ready-to-eat (RTE) foods, associated with long refrigerated shelf-lives, not requiring a heat treatment prior to consumption. Because Listeria monocytogenes is able to grow at refrigeration temperatures, inconsistent temperatures during production, distribution, and at consumer's household may allow for the pathogen to thrive, reaching unsafe limits. L. monocytogenes is the causative agent of listeriosis, a rare but severe human illness, with high fatality rates, transmitted almost exclusively by food consumption. With the aim of assessing the quantitative microbial risk of L. monocytogenes in RTE chicken salads, a challenge test was performed. Salads were inoculated with a three-strain mixture of cold-adapted L. monocytogenes and stored at 4, 12, and 16 °C for eight days. Results revealed that the salad was able to support L. monocytogenes’ growth, even at refrigeration temperatures. The Baranyi primary model was fitted to microbiological data to estimate the pathogen's growth kinetic parameters. Temperature effect on the maximum specific growth rate (μ_max) was modeled using a square-root-type model. Storage temperature significantly influenced μ_max of L. monocytogenes (p < 0.05). These predicted growth models for L. monocytogenes were subsequently used to develop a quantitative microbial risk assessment, estimating a median number of 0.00008726 listeriosis cases per year linked to the consumption of these RTE salads. Sensitivity analysis considering different time–temperature scenarios indicated a very low median risk per portion (<−7 log), even if the assessed RTE chicken salad was kept in abuse storage conditions.     

Growth and Heat Resistance Kinetic Variation Among Various Isolates of Salmonella and its Application to Risk Assessment

The abilities of cells of a particular type of bacteria to leave lag phase and begin the process of dividing or surviving heat treatment can depend on the serotypes or strains of the bacteria. This article reports an investigation of serotype-specific differences in growth and heat resistance kinetics of clinical and food isolates of Salmonella. Growth kinetics at 19 degrees C and 37 degrees C were examined in brain heart infusion broth and heat resistance kinetics for 60 degrees C were examined in beef gravy using a submerged coil heating apparatus. Estimates of the parameters of the growth curves suggests a small between-serotype variance of the growth kinetics. However, for inactivation, the results suggest a significant between-serotype effect on the asymptotic D-values, with an estimated between-serotype CV of about 20%. In microbial risk assessment, predictive microbiology is used to estimate growth and inactivation of pathogens. Often the data used for estimating the growth or inactivation kinetics are based on measurements on a cocktail--a mixture of approximately equal proportions of several serotypes or strains of the pathogen being studied. The expected growth or inactivation rates derived from data using cocktails are biased, reflecting the characteristics of the fastest growing or most heat resistant serotype of the cocktail. In this article, an adjustment to decrease this possible bias in a risk assessment is offered. The article also presents discussion of the effect on estimating growth when stochastic assumptions are incorporated in the model. In particular, equations describing the variation of relative growth are derived, accounting for the stochastic variations of the division of cells. For small numbers of cells, the expected value of the relative growth is not an appropriate "representative" value for actual relative growths that might occur.     

Meta-Analysis of Food Safety Information Based on a Combination of a Relational Database and a Predictive Modeling Tool

The management of microbial risk in food products requires the ability to predict growth kinetics of pathogenic microorganisms in the event of contamination and growth initiation. Useful data for assessing these issues may be found in the literature or from experimental results. However, the large number and variety of data make further development difficult. Statistical techniques, such as meta-analysis, are then useful to realize synthesis of a set of distinct but similar experiences. Moreover, predictive modeling tools can be employed to complete the analysis and help the food safety manager to interpret the data. In this article, a protocol to perform a meta-analysis of the outcome of a relational database, associated with quantitative microbiology models, is presented. The methodology is illustrated with the effect of temperature on pathogenic Escherichia coli and Listeria monocytogenes, growing in culture medium, beef meat, and milk products. Using a database and predictive models, simulations of growth in a given product subjected to various temperature scenarios can be produced. It is then possible to compare food products for a given microorganism, according to its growth ability in these products, and to compare the behavior of bacteria in a given foodstuff. These results can assist decisions for a variety of questions on food safety.     

A Risk Assessment of Campylobacteriosis and Salmonellosis Linked to Chicken Meals Prepared in Households in Dakar,Senegal

We used a quantitative microbiological risk assessment model to describe the risk of Campylobacter and Salmonella infection linked to chicken meals prepared in households in Dakar, Senegal. The model uses data collected specifically for this study, such as the prevalence and level of bacteria on the neck skin of chickens bought in Dakar markets, time‐temperature profiles recorded from purchase to consumption, an observational survey of meal preparation in private kitchens, and detection and enumeration of pathogens on kitchenware and cooks’ hands. Thorough heating kills all bacteria present on chicken during cooking, but cross‐contamination of cooked chicken or ready‐to‐eat food prepared for the meal via kitchenware and cooks’ hands leads to a high expected frequency of pathogen ingestion. Additionally, significant growth of Salmonella is predicted during food storage at ambient temperature before and after meal preparation. These high exposures lead to a high estimated risk of campylobacteriosis and/or salmonellosis in Dakar households. The public health consequences could be amplified by the high level of antimicrobial resistance of Salmonella and Campylobacter observed in this setting. A significant decrease in the number of ingested bacteria and in the risk could be achieved through a reduction of the prevalence of chicken contamination at slaughter, and by the use of simple hygienic measures in the kitchen. There is an urgent need to reinforce the hygiene education of food handlers in Senegal.     

A Model of Hygiene Practices and Consumption Patterns in the Consumer Phase

A mathematical model is presented, which addresses individual hygiene practices during food preparation and consumption patterns in private homes. Further, the model links food preparers and consumers based on their relationship to household types. For different age and gender groups, the model estimates (i) the probability of ingesting a meal where precautions have not been taken to avoid the transfer of microorganisms from raw food to final meal (a risk meal), exemplified by the event that the cutting board was not washed during food preparation, and (ii) the probability of ingesting a risk meal in a private home, where chicken was the prepared food item (a chicken risk meal). Chicken was included in the model, as chickens are believed to be the major source of human exposure to the foodborne pathogen Campylobacter. Monte Carlo simulations showed that the probability of ingesting a risk meal was highest for young males (aged 18-29 years) and lowest for the elderly above 60 years of age. Children aged 0-4 years had a higher probability of ingesting a risk meal than children aged 5-17 years. This difference between age and gender groups was ascribed to the variations in the hygiene levels of food preparers. By including the probability of ingesting a chicken meal at home, simulations revealed that all age groups, except the group above 60 years of age, had approximately the same probability of ingesting a chicken risk meal, the probability of females being slightly higher than that of males. The simulated results show that the probability of ingesting a chicken risk meal at home does not only depend on the hygiene practices of the persons preparing the food, but also on the consumption patterns of consumers, and the relationship between people preparing and ingesting food. This finding supports the need of including information on consumer behavior and preparation hygiene in the consumer phase of exposure assessments.     

An Integrated Model of Infection Risk in a Health-Care Environment

Certain respiratory tract infections can be transmitted by hand-to-mucous-membrane contact, inhalation, and/or direct respiratory droplet spray. In a room occupied by a patient with such a transmissible infection, pathogens present on textile and nontextile surfaces, and pathogens present in the air, provide sources of exposure for an attending health-care worker (HCW); in addition, close contact with the patient when the latter coughs allows for droplet spray exposure. We present an integrated model of pertinent source-environment-receptor pathways, and represent physical elements in these pathways as "states" in a discrete-time Markov chain model. We estimate the rates of transfer at various steps in the pathways, and their relationship to the probability that a pathogen in one state has moved to another state by the end of a specified time interval. Given initial pathogen loads on textile and nontextile surfaces and in room air, we use the model to estimate the expected pathogen dose to a HCW's mucous membranes and respiratory tract. In turn, using a nonthreshold infectious dose model, we relate the expected dose to infection risk. The system is illustrated with a hypothetical but plausible scenario involving a viral pathogen emitted via coughing. We also use the model to show that a biocidal finish on textile surfaces has the potential to substantially reduce infection risk via the hand-to-mucous-membrane exposure pathway.     

Physiologically Based Pharmacokinetic Modeling of the Lactational Transfer of Methylmercury

The developmental neurotoxicity of methylmercury (MeHg) in humans has been described following catastrophic events in Minamata Bay, Japan and in Iraq, and following the exposure to lower doses elsewhere in the world. The most common route of MeHg exposure in humans is through the intake of contaminated food, especially fish. Although the precautions against the ingestion of potentially contaminated food during pregnancy are well recognized, precautions against the ingestion of MeHg during lactation are not so uniformly recognized. However, the continued development of the central nervous system during the early postnatal period serves to prolong the period during which this critical system is susceptible to the toxic insult of MeHg. Because no direct method is available to quantitatively assess the lactational transfer of MeHg to humans, a computer-aided simulation method was developed. An available gestational physiologically based pharmacokinetic model was refined and expanded to include parameters and algorithms specific for the elimination of MeHg in breast milk. The predictions of the completed model were compared with experimental data obtained from rodents, and the model parameters were allometrically scaled to humans. Finally, the model was validated by comparing its predictions against the available clinical data for MeHg distribution and elimination in mothers and their nursing infants. This model incorporated current and previous maternal exposures to MeHg to accurately predict the kinetics of MeHg excretion in breast milk and the daily intake by the nursing infant. This model may be used to quantify MeHg intake by the nursing infant, under different rates of maternal MeHg ingestion.     

Scheduling to Improve Field Service Quality*

This research addresses the problem of scheduling technicians to travel from customer site to customer site to perform emergency maintenance on office machines, computers, robots, telecommunications equipment, medical equipment, heating/cooling equipment, household appliances, and other equipment. We call this the Traveling Technician Problem (TTP). In its simplest form, the TTP is a multiserver, sequence-dependent, tardiness minimization problem. This research frames the TTP as a service quality maximization problem in which service quality is defined in terms of mean tardiness, mean technician phone response time, mean promise time, and mean response time. Tardiness is defined with respect to contractually guaranteed response times. Industry practice is to use dispatching rules to assign service calls to technicians. This research proposes scheduling procedures to maximize field service quality in a dynamic environment. A simulation experiment was used to compare three dispatching rules and three scheduling procedures for the TTP. The scheduling procedures dominated the dispatching rules on all four service quality measures. The proposed scheduling procedures hold promise for improving service quality in a wide variety of field service organizations and in other scheduling environments as well.     

Stochastically Modeling Listeria Monocytogenes Growth in Farm Tank Milk

This article presents a Listeria monocytogenes growth model in milk at the farm bulk tank stage. The main objective was to judge the feasibility and value to risk assessors of introducing a complex model, including a complete thermal model, within a microbial quantitative risk assessment scheme. Predictive microbiology models are used under varying temperature conditions to predict bacterial growth. Input distributions are estimated based on data in the literature, when it is available. If not, reasonable assumptions are made for the considered context. Previously published results based on a Bayesian analysis of growth parameters are used. A Monte Carlo simulation that forecasts bacterial growth is the focus of this study. Three scenarios that take account of the variability and uncertainty of growth parameters are compared. The effect of a sophisticated thermal model taking account of continuous variations in milk temperature was tested by comparison with a simplified model where milk temperature was considered as constant. Limited multiplication of bacteria within the farm bulk tank was modeled. The two principal factors influencing bacterial growth were found to be tank thermostat regulation and bacterial population growth parameters. The dilution phenomenon due to the introduction of new milk was the main factor affecting the final bacterial concentration. The results show that a model that assumes constant environmental conditions at an average temperature should be acceptable for this process. This work may constitute a first step toward exposure assessment for L. monocytogenes in milk. In addition, this partly conceptual work provides guidelines for other risk assessments where continuous variation of a parameter needs to be taken into account.     

Validation of a Stochastic Discrete Event Model Predicting Virus Concentration on Nurse Hands

Understanding healthcare viral disease transmission and the effect of infection control interventions will inform current and future infection control protocols. In this study, a model was developed to predict virus concentration on nurses’ hands using data from a bacteriophage tracer study conducted in Tucson, Arizona, in an urgent care facility. Surfaces were swabbed 2 hours, 3.5 hours, and 6 hours postseeding to measure virus spread over time. To estimate the full viral load that would have been present on hands without sampling, virus concentrations were summed across time points for 3.5‐ and 6‐hour measurements. A stochastic discrete event model was developed to predict virus concentrations on nurses’ hands, given a distribution of virus concentrations on surfaces and expected frequencies of hand‐to‐surface and orifice contacts and handwashing. Box plots and statistical hypothesis testing were used to compare the model‐predicted and experimentally measured virus concentrations on nurses’ hands. The model was validated with the experimental bacteriophage tracer data because the distribution for model‐predicted virus concentrations on hands captured all observed value ranges, and interquartile ranges for model and experimental values overlapped for all comparison time points. Wilcoxon rank sum tests showed no significant differences in distributions of model‐predicted and experimentally measured virus concentrations on hands. However, limitations in the tracer study indicate that more data are needed to instill more confidence in this validation. Next model development steps include addressing viral concentrations that would be found naturally in healthcare environments and measuring the risk reductions predicted for various infection control interventions.     

Physiologically-Based Pharmacokinetic Modeling of a Mixture of Toluene and Xylene in Humans

A physiologically-based pharmacokinetic (PBPK) model for a mixture of toluene (TOL) and xylene (XYL), developed and validated in the rat, was used to predict the uptake and disposition kinetics of TOL/XYL mixture in humans. This was accomplished by substituting the rat physiological parameters and the blood:air partition coefficient with those of humans, scaling the maximal velocity for hepatic metabolism on the basis of body weight^0.75, and keeping all other model parameters species-invariant. The human TOL/XYL mixture PBPK model, developed based on the quantitative biochemical mechanism of interaction elucidated in the rat (i.e., competitive metabolic inhibition), simulated adequately the kinetics of TOL and XYL during combined exposures in humans. The simulations with this PBPK model indicate that an eight hour co-exposure to concentrations that remain within the current threshold limit values of TOL (50 ppm) and XYL (100 ppm) would not result in significant pharmacokinetic interferences, thus implying that data on biological monitoring of worker exposure to these solvents would be unaffected during co-exposures.     

Modeling Individual Cell Lag Time Distributions for Listeria monocytogenes

The food industry faces two paradoxical demands: on the one hand, foods need to be microbiologically safe for consumption and on the other hand, consumers want fresh, minimally processed foods. To meet these demands, more insight into the mechanisms of microbial growth is needed, which includes, among others, the microbial lag phase. This is the time needed by bacterial cells to adapt to a new environment (for example, after food product contamination) before starting an exponential growth regime. Since food products are often contaminated with low amounts of pathogenic microorganisms, it is important to know the distribution of these individual cell lag times to make accurate predictions concerning food safety. More precisely, cells with the shortest lag times (i.e., appearing in the left tail of the distribution) are largely decisive for the outgrowth of the population. In this study, an integrated modeling approach is proposed and applied to an existing data set of individual cell lag time measurements of Listeria monocytogenes. In a first step, a logistic modeling approach is applied to predict the fraction of zero-lag cells (which start growing immediately) as a function of temperature, pH, and water activity. For the nonzero-lag cells, the mean and variance of the lag time distribution are modeled with a hyperbolic-type model structure. This mean and variance allow identification of the parameters of a two-parameter Weibull distribution, representing the nonzero-lag cell lag time distribution. The integration of the developed models allows prediction of a global distribution of individual cell lag times for any combination of environmental conditions in the interpolation domain of the original temperature, pH, and water activity settings. The global fitting quality of the model is quantified using several measures indicating that the model gives accurate predictions, erring slightly on the fail-safe side when predicting the shortest lag times.     

Combining Quantitative Risk Assessment of Human Health,Food Waste,and Energy Consumption: The Next Step in the Development of the Food Cold Chain?

The preservation of perishable food via refrigeration in the supply chain is essential to extend shelf life and provide consumers with safe food. However, electricity consumed in refrigeration processes has an economical and an environmental impact. This study focuses on the cold chain of cooked ham, including transport, cold room in supermarket, display cabinet, transport by consumer, and domestic refrigerator, and aims to predict the risk for human health associated with Listeria monocytogenes, the amount of food wasted due to the growth of spoilage bacteria, and the electrical consumption to maintain product temperature through the cold chain. A set of eight intervention actions were tested to evaluate their impact on the three criteria. Results show that the modification of the thermostat of the domestic refrigerator has a high impact on food safety and food waste and a limited impact on the electrical consumption. Inversely, the modification of the airflow rate in the display cabinet has a high impact on electrical consumption and a limited impact on food safety and food waste. A cost–benefit analysis approach and two multicriteria decision analysis methods were used to rank the intervention actions. These three methodologies show that setting the thermostat of the domestic refrigerator to 4 °C presents the best compromise between the three criteria. The impact of decisionmaker preferences (criteria weight) and limitations of these three approaches are discussed. The approaches proposed by this study may be useful in decision making to evaluate global impact of intervention actions in issues involving conflicting outputs.     

Development of the PEARLS model (Particulate Exposure from Ambient to Regional Lung by Subgroup) and use of Monte Carlo simulation to predict internal exposure to 2.5 in Toronto.

Air pollution is a current and growing concern for Canadians, and there is evidence that ambient levels that meet current exposure standards may be associated with mortality and morbidity in Toronto, Canada. Evaluating exposure is an important step in understanding the relationship between particulate matter (PM) exposure and health outcomes. This report describes the PEARLS model (Particulate Exposure from Ambient to Regional Lung by Subgroup), which predicts exposure distributions for 11 age-gender population subgroups in Toronto to PM2.5 (PM with a median aerodynamic diameter of 2.5 microm or less) using Monte Carlo simulation techniques. The model uses physiological and activity pattern characteristics of each subgroup to determine region-specific lung exposure to PM2.5, which is defined as the mass of PM2.5 deposited per unit time to each of five lung regions (two extrathoracic, bronchial, bronchiolar, and alveolar). The modeling results predict that children, toddlers, and infants have the broadest distributions of exposure, and the greatest chance of experiencing extreme exposures in the alveolar region of the lung. Importance analysis indicates that the most influential model variables are air exchange rate into indoor environments, time spent outdoors, and time spent at high activity levels. Additionally, a "critical point" was defined and introduced to the PEARLS to investigate the effects of possible threshold-pathogenic phenomena on subgroup exposure patterns. The analysis indicates that the subgroups initially predicted to be most highly exposed were likely to have the highest proportion of their population exposed above the critical point. Substantial exposures above the critical point were predicted in all subgroups for ambient concentrations of PM2.5 commonly observed in Toronto after continuous exposure of 24 hours or more.     

Consumer Phase Risk Assessment for Listeria monocytogenes in Deli Meats

The foodborne disease risk associated with the pathogen Listeria monocytogenes has been the subject of recent efforts in quantitative microbial risk assessment. Building upon one of these efforts undertaken jointly by the U.S. Food and Drug Administration and the U.S. Department of Agriculture (USDA), the purpose of this work was to expand on the consumer phase of the risk assessment to focus on handling practices in the home. One-dimensional Monte Carlo simulation was used to model variability in growth and cross-contamination of L. monocytogenes during food storage and preparation of deli meats. Simulations approximated that 0.3% of the servings were contaminated with >10(4) CFU/g of L. monocytogenes at the time of consumption. The estimated mean risk associated with the consumption of deli meats for the intermediate-age population was approximately 7 deaths per 10(11) servings. Food handling in homes increased the estimated mean mortality by 10(6)-fold. Of all the home food-handling practices modeled, inadequate storage, particularly refrigeration temperatures, provided the greatest contribution to increased risk. The impact of cross-contamination in the home was considerably less. Adherence to USDA Food Safety and Inspection Service recommendations for consumer handling of ready-to-eat foods substantially reduces the risk of listeriosis.     

Practice-Specific Risk Perceptions and Self-Reported Food Safety Practices

The relationship between risk perception and risk avoidance is typically analyzed using self-reported measures. However, in domains such as driving or food handling, the validity of responses about usual behavior is threatened because people think about the situations in which they are self-aware, such as when they encounter a hazard. Indeed, researchers have often noted a divergence between what people say about their behavior and how they actually behave. Thus, in order to draw conclusions about risk perceptions and risk avoidance from survey data, it is important to identify particular cognitive elements, such as those measured by questions about risk and safety knowledge, risk perceptions, or information search behavior, which may be effective antecedents of self-reported safety behavior. It is also important to identify and correct for potential sources of bias that may exist in the data. The authors analyze the Food and Drug Administration's 1998 Food Safety Survey to determine whether there are consistent cognitive antecedents for three types of safe food practices: preparation, eating, and cooling of foods. An assessment of measurement biases shows that endogeneity of food choices affects reports of food preparation. In addition, response bias affects reports of cooling practices as evidenced by its relation to knowledge and information search, a pattern of cognitive effects unique to cooling practices. After correcting for these biases, results show that practice-specific risk perceptions are the primary cognitive antecedents of safe food behavior, which has implications for the design of effective education messages about food safety.     

A Physiologically Based Pharmacokinetic Assessment of Tetrachloroethylene in Groundwater for a Bathing and Showering Determination

A two-step methodology is described to make a health-based determination for the bathing and showering use of the water from a private well contaminated with volatile organic chemicals. The chemical perchloroethylene (PERC) is utilized to illustrate the approach. First, a chemical-specific exposure model is used to predict the concentration of PERC in the shower air, shower water, and in the air above the bathtub. Second, a physiologically based pharmacokinetic (PBPK) model is used to predict the concentration of PERC delivered to the target tissue, the brain, since the focus is on neurological endpoints. The simulation exercise includes concurrent dermal and inhalation routes of exposure. A reference target tissue level (RTTL) in the brain is estimated using the PBPK model. A hazard index based on this benchmark guideline is used to make a regulatory determination for bathing and showering use of the contaminated water.     

A Compartmental Model for the Prediction of Breath Concentration and Absorbed Dose of Chloroform After Exposure While Showering

In order to predict the exhaled breath concentration of chloroform in individuals exposed to chloroform while showering, an existing physiologically based pharmacokinetic (PB-PK) model was modified to include a multicompartment, PB-PK model for the skin and a completely mixed shower exposure model. The PB-PK model of the skin included the stratum corneum as the principal resistance to absorption and a viable epidermis which is in dynamic equilibrium with the skin microcirculation. This model was calibrated with measured exhaled breath concentrations of chloroform in individuals exposed while showering with and without dermal absorption. The calibration effort indicated that the expected value of skin-blood partitioning coefficient would be 1.2 when the degree of transfer of chloroform from shower water into shower air was 61%. The stratum corneum permeability coefficient for chloroform was estimated to be within the range of 0.16-0.36 cm/hr and the expected value was 0.2 cm/hr. The estimated ratio of the dermally and inhaled absorbed doses ranged between 0.6 and 2.2 and the expected value was 0.75. These results indicate that for the purposes of risk assessment for dermal exposure to chloroform, a simple steady-state model can be used to predict the degree of dermal absorption and that a reasonable value of skin permeability coefficient for chloroform used in this model would be 0.2 cm/hr. Further research should be conducted to compare the elimination of chloroform via exhaled breath when different exposure routes are being compared. The model results from this study suggest that multiple measurements of exhaled breath concentrations after exposure may be necessary when making comparisons of breath concentrations that involve different exposure routes.     

Modeling Rabbit Responses to Single and Multiple Aerosol Exposures of Bacillus anthracis Spores

Survival models are developed to predict response and time‐to‐response for mortality in rabbits following exposures to single or multiple aerosol doses of Bacillus anthracis spores. Hazard function models were developed for a multiple‐dose data set to predict the probability of death through specifying functions of dose response and the time between exposure and the time‐to‐death (TTD). Among the models developed, the best‐fitting survival model (baseline model) is an exponential dose–response model with a Weibull TTD distribution. Alternative models assessed use different underlying dose–response functions and use the assumption that, in a multiple‐dose scenario, earlier doses affect the hazard functions of each subsequent dose. In addition, published mechanistic models are analyzed and compared with models developed in this article. None of the alternative models that were assessed provided a statistically significant improvement in fit over the baseline model. The general approach utilizes simple empirical data analysis to develop parsimonious models with limited reliance on mechanistic assumptions. The baseline model predicts TTDs consistent with reported results from three independent high‐dose rabbit data sets. More accurate survival models depend upon future development of dose–response data sets specifically designed to assess potential multiple‐dose effects on response and time‐to‐response. The process used in this article to develop the best‐fitting survival model for exposure of rabbits to multiple aerosol doses of B. anthracis spores should have broad applicability to other host–pathogen systems and dosing schedules because the empirical modeling approach is based upon pathogen‐specific empirically‐derived parameters.     

Topics in Microbial Risk Assessment: Dynamic Flow Tree Process

Microbial risk assessment is emerging as a new discipline in risk assessment. A systematic approach to microbial risk assessment is presented that employs data analysis for developing parsimonious models and accounts formally for the variability and uncertainty of model inputs using analysis of variance and Monte Carlo simulation. The purpose of the paper is to raise and examine issues in conducting microbial risk assessments. The enteric pathogen Escherichia coli O157:H7 was selected as an example for this study due to its significance to public health. The framework for our work is consistent with the risk assessment components described by the National Research Council in 1983 (hazard identification; exposure assessment; dose-response assessment; and risk characterization). Exposure assessment focuses on hamburgers, cooked a range of temperatures from rare to well done, the latter typical for fast food restaurants. Features of the model include predictive microbiology components that account for random stochastic growth and death of organisms in hamburger. For dose-response modeling, Shigella data from human feeding studies were used as a surrogate for E. coli O157:H7. Risks were calculated using a threshold model and an alternative nonthreshold model. The 95% probability intervals for risk of illness for product cooked to a given internal temperature spanned five orders of magnitude for these models. The existence of even a small threshold has a dramatic impact on the estimated risk.