Quantitative Risk Assessment of Listeria monocytogenes in French Cold-Smoked Salmon: II. Risk Characterization

A model for the assessment of exposure to Listeria monocytogenes from cold-smoked salmon consumption in France was presented in the first of this pair of articles (Pouillot et al ., 2007, Risk Analysis, 27:683–700). In the present study, the exposure model output was combined with an internationally accepted hazard characterization model, adapted to the French situation, to assess the risk of invasive listeriosis from cold-smoked salmon consumption in France in a second-order Monte Carlo simulation framework. The annual number of cases of invasive listeriosis due to cold-smoked salmon consumption in France is estimated to be 307, with a very large credible interval ([10; 12,453]), reflecting data uncertainty. This uncertainty is mainly associated with the dose-response model. Despite the significant uncertainty associated with the predictions, this model provides a scientific base for risk managers and food business operators to manage the risk linked to cold-smoked salmon contaminated with L. monocytogenes. Under the modeling assumptions, risk would be efficiently reduced through a decrease in the prevalence of L. monocytogenes or better control of the last steps of the cold chain (shorter and/or colder storage during the consumer step), whereas reduction of the initial contamination levels of the contaminated products and improvement in the first steps of the cold chain do not seem to be promising strategies. An attempt to apply the recent risk-based concept of FSO (food safety objective) on this example underlines the ambiguity in practical implementation of the risk management metrics and the need for further elaboration on these concepts.     

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.     

Modeling Logistic Performance in Quantitative Microbial Risk Assessment

In quantitative microbial risk assessment (QMRA), food safety in the food chain is modeled and simulated. In general, prevalences, concentrations, and numbers of microorganisms in media are investigated in the different steps from farm to fork. The underlying rates and conditions (such as storage times, temperatures, gas conditions, and their distributions) are determined. However, the logistic chain with its queues (storages, shelves) and mechanisms for ordering products is usually not taken into account. As a consequence, storage times—mutually dependent in successive steps in the chain—cannot be described adequately. This may have a great impact on the tails of risk distributions. Because food safety risks are generally very small, it is crucial to model the tails of (underlying) distributions as accurately as possible. Logistic performance can be modeled by describing the underlying planning and scheduling mechanisms in discrete-event modeling. This is common practice in operations research, specifically in supply chain management. In this article, we present the application of discrete-event modeling in the context of a QMRA for  Listeria monocytogenes  in fresh-cut iceberg lettuce. We show the potential value of discrete-event modeling in QMRA by calculating logistic interventions (modifications in the logistic chain) and determining their significance with respect to food safety.     

Risk Assessment of Listeriosis Linked to the Consumption of Two Soft Cheeses Made from Raw Milk: Camembert of Normandy and Brie of Meaux

This article reports a quantitative risk assessment of human listeriosis linked to the consumption of soft cheeses made from raw milk. Risk assessment was based on data purposefully acquired inclusively over the period 2000-2001 for two French cheeses, namely: Camembert of Normandy and Brie of Meaux. Estimated Listeria monocytogenes concentration in raw milk was on average 0.8 and 0.3 cells/L, respectively, in Normandy and Brie regions. A Monte Carlo simulation was used to account for the time-temperature history of the milk and cheeses from farm to table. It was assumed that cell progeny did not spread within the solid cheese matrix (as they would be free to do in liquid broth). Interaction between pH and temperature was accounted for in the growth model. The simulated proportion of servings with no L. monocytogenes cell was 88% for Brie and 82% for Camembert. The 99th percentile of L. monocytogenes cell numbers in servings of 27 g of cheese was 131 for Brie and 77 for Camembert at the time of consumption, corresponding respectively to three and five cells of L. monocytogenes per gram. The expected number of severe listeriosis cases would be < or =10(-3) and < or =2.5 x 10(-3) per year for 17 million servings of Brie of Meaux and 480 million servings of Camembert of Normandy, respectively.     

A Simulation Study of Quantitative Risk Assessment for Bivariate Continuous Outcomes

The neurotoxic effects of chemical agents are often investigated in controlled studies on rodents, with binary and continuous multiple endpoints routinely collected. One goal is to conduct quantitative risk assessment to determine safe dose levels. Yu and Catalano (2005) describe a method for quantitative risk assessment for bivariate continuous outcomes by extending a univariate method of percentile regression. The model is likelihood based and allows for separate dose‐response models for each outcome while accounting for the bivariate correlation. The approach to benchmark dose (BMD) estimation is analogous to that for quantal data without having to specify arbitrary cutoff values. In this article, we evaluate the behavior of the BMD relative to background rates, sample size, level of bivariate correlation, dose‐response trend, and distributional assumptions. Using simulations, we explore the effects of these factors on the resulting BMD and BMDL distributions. In addition, we illustrate our method with data from a neurotoxicity study of parathion exposure in rats.     

Risk Assessment of Listeria monocytogenes: Impact of Individual Cell Variability on the Exposure Assessment Step

Recently, the lag phase research in predictive microbiology is focusing more on the individual cell variability, especially for pathogenic microorganisms that typically occur in very low contamination levels, like Listeria monocytogenes. In this study, the effect of this individual cell lag phase variability was introduced in an exposure assessment study for L. monocytogenes in a liver paté. A basic framework was designed to estimate the contamination level of paté at the time of consumption, taking into account the frequency of contamination and the initial contamination levels of paté at retail. Growth was calculated on paté units of 150 g, comparing an individual-based approach with a classical population-based approach. The two different protocols were compared using simulations. If only the individual cell lag variability was taken into account, important differences were observed in cell density at the time of consumption between the individual-based approach and the classical approach, especially at low inoculum levels, resulting in high variability when using the individual-based approach. Although, when all variable factors were taken into account, no significant differences were observed between the different approaches, allowing the conclusion that the individual cell lag phase variability was overruled by the global variability of the exposure assessment framework. Even in more extreme conditions like a low inoculum level or a low water activity, no differences were created in cell density at the time of consumption between the individual-based approach and the classical approach. This means that the individual cell lag phase variability of L. monocytogenes has important consequences when studying specific growth cases, especially when the applied inoculum levels are low, but when performing more general exposure assessment studies, the variability between the individual cell lag phases is too limited to have a major impact on the total exposure assessment.     

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.     

A Semiparametric Approach to Risk Assessment for Quantitative Outcomes

Characterizing the dose-effect relationship and estimating acceptable exposure levels are the primary goals of quantitative risk assessment. A semiparametric approach is proposed for risk assessment with continuously measured or quantitative outcomes which has advantages over existing methods by requiring fewer assumptions. The approach is based on pairwise ranking between the response values in the control group and those in the exposed groups. The work generalizes the rank-based Wilcoxon-Mann-Whitney test, which for the two-group comparison is effectively a test of whether a response from the control group is different from (larger than) a response in an exposed group. We develop a regression framework that naturally extends this metric to model the dose effect in terms of a risk function. Parameters of the regression model can be estimated with standard software. However, inference requires an additional step to estimate the variance structure of the estimated parameters. An effective dose (ED) and associated lower confidence limit (LED) are easily calculated. The method is supported by a simulation study and is illustrated with a study on the effects of aconiazide. The method offers flexible modeling of the dose effect, and since it is rank-based, it is more resistant to outliers, nonconstant variance, and other departures from normality than previously described approaches.     

Quantitative Assessment of Human MRSA Risks from Swine

The public health community, news media, and members of the general public have expressed significant concern that methicillin‐resistant Staphylococcus aureus (MRSA) transmitted from pigs to humans may harm human health. Studies of the prevalence and dynamics of swine‐associated (ST398) MRSA have sampled MRSA at discrete points in the presumed causative chain leading from swine to human patients, including sampling bacteria from live pigs, retail meats, farm workers, and hospital patients. Nonzero prevalence is generally interpreted as indicating a potential human health hazard from MRSA infections, but quantitative assessments of resulting risks are not usually provided. This article integrates available data from several sources to construct a conservative (plausible upper bound) probability estimate for the actual human health harm (MRSA infections and fatalities) arising from ST398‐MRSA from pigs. The model provides plausible upper bounds of approximately one excess human infection per year among all U.S. pig farm workers, and one human infection per 31 years among the remaining total population of the United States. These results assume the possibility of transmission events not yet observed, so additional data collection may reduce these estimates further.     

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.     

Quantitative Assessment of the Risk of Lung Cancer Associated with Occupational Exposure to Refractory Ceramic Fibers

We present the results of a quantitative assessment of the lung cancer risk associated with occupational exposure to refractory ceramic fibers (RCF). The primary sources of data for our risk assessment were two long-term oncogenicity studies in male Fischer rats conducted to assess the potential pathogenic effects associated with prolonged inhalation of RCF. An interesting feature of the data was the availability of the temporal profile of fiber burden in the lungs of experimental animals. Because of this information, we were able to conduct both exposure–response and dose–response analyses. Our risk assessment was conducted within the framework of a biologically based model for carcinogenesis, the two-stage clonal expansion model, which allows for the explicit incorporation of the concepts of initiation and promotion in the analyses. We found that a model positing that RCF was an initiator had the highest likelihood. We proposed an approach based on biological considerations for the extrapolation of risk to humans. This approach requires estimation of human lung burdens for specific exposure scenarios, which we did by using an extension of a model due to Yu. Our approach acknowledges that the risk associated with exposure to RCF depends on exposure to other lung carcinogens. We present estimates of risk in two populations: (1) a population of nonsmokers and (2) an occupational cohort of steelworkers not exposed to coke oven emissions, a mixed population that includes both smokers and nonsmokers.     

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.     

A Quantitative Microbial Risk Assessment Model for Legionnaires'' Disease: Animal Model Selection and Dose-Response Modeling

Legionnaires' disease (LD), first reported in 1976, is an atypical pneumonia caused by bacteria of the genus Legionella, and most frequently by L. pneumophila (Lp). Subsequent research on exposure to the organism employed various animal models, and with quantitative microbial risk assessment (QMRA) techniques, the animal model data may provide insights on human dose-response for LD. This article focuses on the rationale for selection of the guinea pig model, comparison of the dose-response model results, comparison of projected low-dose responses for guinea pigs, and risk estimates for humans. Based on both in vivo and in vitro comparisons, the guinea pig (Cavia porcellus) dose-response data were selected for modeling human risk. We completed dose-response modeling for the beta-Poisson (approximate and exact), exponential, probit, logistic, and Weibull models for Lp inhalation, mortality, and infection (end point elevated body temperature) in guinea pigs. For mechanistic reasons, including low-dose exposure probability, further work on human risk estimates for LD employed the exponential and beta-Poisson models. With an exposure of 10 colony-forming units (CFU) (retained dose), the QMRA model predicted a mild infection risk of 0.4 (as evaluated by seroprevalence) and a clinical severity LD case (e.g., hospitalization and supportive care) risk of 0.0009. The calculated rates based on estimated human exposures for outbreaks used for the QMRA model validation are within an order of magnitude of the reported LD rates. These validation results suggest the LD QMRA animal model selection, dose-response modeling, and extension to human risk projections were appropriate.     

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 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.     

Disparity in Quantitative Risk Assessment: A Review of Input Distributions

Monte Carlo simulations are commonplace in quantitative risk assessments (QRAs). Designed to propagate the variability and uncertainty associated with each individual exposure input parameter in a quantitative risk assessment, Monte Carlo methods statistically combine the individual parameter distributions to yield a single, overall distribution. Critical to such an assessment is the representativeness of each individual input distribution. The authors performed a literature review to collect and compare the distributions used in published QRAs for the parameters of body weight, food consumption, soil ingestion rates, breathing rates, and fluid intake. To provide a basis for comparison, all estimated exposure parameter distributions were evaluated with respect to four properties: consistency, accuracy, precision, and specificity. The results varied depending on the exposure parameter. Even where extensive, well-collected data exist, investigators used a variety of different distributional shapes to approximate these data. Where such data do not exist, investigators have collected their own data, often leading to substantial disparity in parameter estimates and subsequent choice of distribution. The present findings indicate that more attention must be paid to the data underlying these distributional choices. More emphasis should be placed on sensitivity analyses, quantifying the impact of assumptions, and on discussion of sources of variation as part of the presentation of any risk assessment results. If such practices and disclosures are followed, it is believed that Monte Carlo simulations can greatly enhance the accuracy and appropriateness of specific risk assessments. Without such disclosures, researchers will be increasing the size of the risk assessment "black box," a concern already raised by many critics of more traditional risk assessments.     

The EPA Health Risk Assessment of Methylcyclopentadienyl Manganese Tricarbonyl (MMT)

This paper describes the U.S. Environmental Protection Agency's assessment of potential health risks associated with the possible widespread use of a manganese (Mn)-based fuel additive, methylcyclopentadienyl manganese tricarbonyl (MMT). This assessment was significant in several respects and may be instructive in identifying certain methodological issues of general relevance to risk assessment. A major feature of the inhalation health risk assessment was the derivation of Mn inhalation reference concentration (RfC) estimates using various statistical approaches, including benchmark dose and Bayesian analyses. The exposure assessment component used data from the Particle Total Exposure Assessment Methodology (PTEAM) study and other sources to estimate personal exposure levels of particulate Mn attributable to the permitted use of MMT in leaded gasoline in Riverside, CA, at the time of the PTEAM study; on this basis it was then possible to predict a distribution of possible future exposure levels associated with the use of MMT in all unleaded gasoline. Qualitative as well as quantitative aspects of the risk characterization are summarized, along with inherent uncertainties due to data limitations.     

Quantitative Assessment of Exposure to the Mycotoxin Ochratoxin A in Food

This article presents the methodology and the simulation results concerning the quantitative assessment of exposure to the fungus toxin named Ochratoxin A (OA) in food, in humans in France. We show that is possible to provide reliable calculations of exposure to OA with the conjugate means of a nonparametric-type method of simulation, a parametric-type method of simulation, and the use of bootstrap confidence intervals. In the context of the Monte Carlo simulation, the nonparametric method takes into account the consumptions and the contaminations in the simulations only via the raw data whereas the parametric method depends on the random samplings from distribution functions fitted to consumption and contamination data. Our conclusions are based on eight types of food only. Nevertheless, they are meaningful due to the major importance of these foodstuffs in human nourishment in France. This methodology can be applied whatever the food contaminant (pesticides, other mycotoxins, Cadmium, etc.) when data are available.     

Implications of Developmental Toxicity Study Design for Quantitative Risk Assessment

Standard experimental designs for conducting developmental toxicity studies typically include three- or four-dose levels in addition to a control group. Some researchers have suggested that designs with more exposure groups would improve dose-response characterization and risk estimation. Such proposals have not, however, been supported by the results of simulation studies, which instead back the use of fewer dose levels. This discrepancy is partly due to using a known dose–response pattern to generate data, making model choice obvious. While the carcinogenicity literature has explored implications of different study designs, little attention has been given to the role of design in developmental toxicity risk assessment (or noncancer toxicology in general). In this research, we explore the implications of various experimental designs for developmental toxicity by resampling data from a large study of 2,4,5-trichlorophenoxyacetic acid in mice. We compare the properties of benchmark dose (BMD) estimation for different design strategies by randomly selecting animals within particular dose groups from the entire 2,4,5-T database of over 77,000 birth outcomes to create smaller "pseudo-studies" that are representative of standard bioassay sample sizes. Our results show that experimental designs which include more dose levels have advantages in terms of risk characterization and estimation.     

A Risk Assessment for Occupational Acrylonitrile Exposure Using Epidemiology Data

The extensive data from the Blair et al.((1)) epidemiology study of occupational acrylonitrile exposure among 25460 workers in eight plants in the United States provide an excellent opportunity to update quantitative risk assessments for this widely used commodity chemical. We employ the semiparametric Cox relative risk (RR) regression model with a cumulative exposure metric to model cause-specific mortality from lung cancer and all other causes. The separately estimated cause-specific cumulative hazards are then combined to provide an overall estimate of age-specific mortality risk. Age-specific estimates of the additional risk of lung cancer mortality associated with several plausible occupational exposure scenarios are obtained. For age 70, these estimates are all markedly lower than those generated with the cancer potency estimate provided in the USEPA acrylonitrile risk assessment.((2)) This result is consistent with the failure of recent occupational studies to confirm elevated lung cancer mortality among acrylonitrile-exposed workers as was originally reported by O'Berg,((3)) and it calls attention to the importance of using high-quality epidemiology data in the risk assessment process.