Using Average Lifetime Dose Rate for Intermittent Exposures to Carcinogens

The effect of using the average dose rate over a lifetime as a representative measure of exposure to carcinogens is investigated by comparing the true theoretical multistage intermittent-dosing lifetime low-dose excess risk to the theoretical multistage continuous-dosing lifetime risk corresponding to the average lifetime dose rate. It is concluded that low-dose risk estimates based on the average lifetime dose rate may overestimate the true risk by several orders of magnitude, but that they never underestimate the true risk by more than a factor of k/r, where k is the total number of stages in the multistage model and r is the number of stages that are dose-related.     

Using the Biological Two-Stage Model to Assess Risk from Short-Term Exposures

Two assumptions used in risk assessment are investigated: (1) the assumption of fraction of lifetime dose rate assumes that the risk from a fractional lifetime exposure at a given dose rate is equal to the risk from full lifetime exposure at that same fraction of the given dose rate; (2) the assumption of fraction of lifetime risk assumes that the risk from a fractional lifetime exposure at a given dose rate is equal to that same fraction of the risk from full lifetime exposure at the same dose rate. These two assumptions are equivalent when risk is a linear function of dose. Thus both can be thought of as generalizations of the assumption that cancer risk is proportional to the total accumulated lifetime dose (or average daily dose), which is often made to assess the risk from short-term exposures. In this paper, the age-specific cumulative hazard functions are derived using the two-stage model developed by Moolgavkar, Venzon, and Knudson for situations when the exposure occurs during a single period or a single instant. The two assumptions described above are examined for three types of carcinogens, initiator, completer, and promoter, in the context of the model. For initiator and completer, these two assumptions are equivalent in the low-dose region; for a promoter, using the fraction of lifetime risk assumption is generally more conservative than that of the fraction of lifetime dose rate assumption. Tables are constructed to show that the use of either the fraction of lifetime dose rate assumption or the fraction lifetime risk assumption can both underestimate and overestimate the true risk for the three types of carcinogens.     

The Multistage Model with a Time-Dependent Dose Pattern: Applications to Carcinogenic Risk Assessment1

A cancer risk assessment methodology based upon the Armitage–Doll multistage model of cancer is applied to animal bioassay data. The method utilizes the exact time-dependent dose pattern used in a bioassay rather than some single measure of dose such as average dose rate or cumulative dose. The methodology can be used to predict risks from arbitrary exposure patterns including, for example, intermittent exposure and short-term exposure occurring at an arbitrary age. The methodology is illustrated by applying it to a National Cancer Institute bioassay of ethylene dibromide in which dose rates were modified several times during the course of the experiment.     

Dose-Response and Risk Assessment of Airborne Hexavalent Chromium and Lung Cancer Mortality

This study evaluates the dose-response relationship for inhalation exposure to hexavalent chromium [Cr(VI)] and lung cancer mortality for workers of a chromate production facility, and provides estimates of the carcinogenic potency. The data were analyzed using relative risk and additive risk dose-response models implemented with both Poisson and Cox regression. Potential confounding by birth cohort and smoking prevalence were also assessed. Lifetime cumulative exposure and highest monthly exposure were the dose metrics evaluated. The estimated lifetime additional risk of lung cancer mortality associated with 45 years of occupational exposure to 1 microg/m3 Cr(VI) (occupational exposure unit risk) was 0.00205 (90%CI: 0.00134, 0.00291) for the relative risk model and 0.00216 (90%CI: 0.00143, 0.00302) for the additive risk model assuming a linear dose response for cumulative exposure with a five-year lag. Extrapolating these findings to a continuous (e.g., environmental) exposure scenario yielded an environmental unit risk of 0.00978 (90%CI: 0.00640, 0.0138) for the relative risk model [e.g., a cancer slope factor of 34 (mg/kg-day)-1] and 0.0125 (90%CI: 0.00833, 0.0175) for the additive risk model. The relative risk model is preferred because it is more consistent with the expected trend for lung cancer risk with age. Based on statistical tests for exposure-related trend, there was no statistically significant increased lung cancer risk below lifetime cumulative occupational exposures of 1.0 mg-yr/m3, and no excess risk for workers whose highest average monthly exposure did not exceed the current Permissible Exposure Limit (52 microg/m3). It is acknowledged that this study had limited power to detect increases at these low exposure levels. These cancer potency estimates are comparable to those developed by U.S. regulatory agencies and should be useful for assessing the potential cancer hazard associated with inhaled Cr(VI).     

Model Uncertainty and Risk Estimation for Experimental Studies of Quantal Responses

Experimental animal studies often serve as the basis for predicting risk of adverse responses in humans exposed to occupational hazards. A statistical model is applied to exposure-response data and this fitted model may be used to obtain estimates of the exposure associated with a specified level of adverse response. Unfortunately, a number of different statistical models are candidates for fitting the data and may result in wide ranging estimates of risk. Bayesian model averaging (BMA) offers a strategy for addressing uncertainty in the selection of statistical models when generating risk estimates. This strategy is illustrated with two examples: applying the multistage model to cancer responses and a second example where different quantal models are fit to kidney lesion data. BMA provides excess risk estimates or benchmark dose estimates that reflects model uncertainty.     

Estimates of Lifetime-Absorbed Daily Doses from the Use of Personal-Care Products Containing Polyacrylamide: A Monte Carlo Analysis

Estimates of the lifetime-absorbed daily dose (LADD) of acrylamide resulting from use of representative personal-care products containing polyacrylamides have been developed. All of the parameters that determine the amount of acrylamide absorbed by an individual vary from one individual to another. Moreover, for some parameters there is uncertainty as to which is the correct or representative value from a range of values. Consequently, the parameters used in the estimation of the LADD of acrylamide from usage of a particular product type (e.g., deodorant, makeup, etc.) were represented by distributions evaluated using Monte Carlo analyses.((1-4)) From these data, distributions of values for key parameters, such as the amount of acrylamide in polyacrylamide, absorption fraction, etc., were defined and used to provide a distribution of LADDs for each personal-care product. The estimated total acrylamide LADD (across all products) for males and females at the median, mean, and 95th percentile of the distribution of individual LADD values were 4.7 x 10(-8), 2.3 x 10(-7), and 7.3 x 10(-7) mg/kg/day for females and 3.6 x 10(-8), 1.7 x 10(-7), and 5.4 x 10(-7) mg/kg/day for males. The ratio of the LADDs to risk-specific dose corresponding to a target risk level of 1 x 10(-5), the acceptable risk level for this investigation, derived using approaches typically used by the FDA, the USEPA, and proposed for use by the European Union (EU) were also calculated. All ratios were well below 1, indicating that all the extra lifetime cancer risk from the use of polyacrylamide-containing personal-care products, in the manner assumed in this assessment, are well below acceptable levels. Even if it were assumed that an individual used all of the products together, the estimated LADD would still provide a dose that was well below the acceptable risk levels.     

Carcinogenic Risk Assessment with Time-Dependent Exposure Patterns

Previous applications of carcinogenic risk assessment using mathematical models of carcinogenesis have focused largely on the case where the level of exposure remains constant over time. In many situations, however, the dose of the carcinogen varies with time. In this paper, we discuss both the classical Armitage-Doll multistage model and the Moolgavkar-Venzon-Knudson two-stage birth-death-mutation model with time-dependent dosing regimens. Bounds on the degree of underestimation of risk that can occur through the use of a simple time-weighted average dose are derived by means of comparison with an equivalent constant dose corresponding to the actual risk under the time-dependent dosing regimen.     

Leukemia Risk Associated with Benzene Exposure in the Pliofilm Cohort. II. Risk Estimates

The detailed work histories of the individual workers composing the Pliofilm cohort represent a unique resource for estimating the dose-respoonse for leukemia that may follow occupational exposure to benzene. In this paper, we report the results of analyzing the updated Pliofilm cohort using the proportional hazards model, a more sophisticated technique that uses more of the available exposure data than the conditional logistic model used by Rinsky et al. The more rigorously defined exposure estimates derived by Paustenbach et al. are consistent with those of Crump and Allen in giving estimates of the slope of the leukemogenic dose-response that are not as steep as the slope resulting from the exposure estimates of Rinsky et al. We consider estimates of 0.3-0.5 additional leukemia deaths per thousand workers with 45 ppm-years of cumulative benzene exposure to be the best estimates currently available of leukemia risk from occupational exposure to benzene. These risks were estimated in the proportional hazards model when the exposure estimates of Crump and Allen or of Paustenbach et al. were used to derive a cumulative concentration-by-time metric.     

Application of Geostatistics to Risk Assessment

Geostatistics offers two fundamental contributions to environmental contaminant exposure assessment: (1) a group of methods to quantitatively describe the spatial distribution of a pollutant and (2) the ability to improve estimates of the exposure point concentration by exploiting the geospatial information present in the data. The second contribution is particularly valuable when exposure estimates must be derived from small data sets, which is often the case in environmental risk assessment. This article addresses two topics related to the use of geostatistics in human and ecological risk assessments performed at hazardous waste sites: (1) the importance of assessing model assumptions when using geostatistics and (2) the use of geostatistics to improve estimates of the exposure point concentration (EPC) in the limited data scenario. The latter topic is approached here by comparing design-based estimators that are familiar to environmental risk assessors (e.g., Land's method) with geostatistics, a model-based estimator. In this report, we summarize the basics of spatial weighting of sample data, kriging, and geostatistical simulation. We then explore the two topics identified above in a case study, using soil lead concentration data from a Superfund site (a skeet and trap range). We also describe several areas where research is needed to advance the use of geostatistics in environmental risk assessment.     

Uncertainty and Variation in Indirect Exposure Assessments: An Analysis of Exposure to Tetrachlorodibenzo-p-Dioxin from a Beef Consumption Pathway

Indirect exposures to 2,3,7,8-tetrachlorodibenzo- p -dioxin (TCDD) and other toxic materials released in incinerator emissions have been identified as a significant concern for human health. As a result, regulatory agencies and researchers have developed specific approaches for evaluating exposures from indirect pathways. This paper presents a quantitative assessment of the effect of uncertainty and variation in exposure parameters on the resulting estimates of TCDD dose rates received by individuals indirectly exposed to incinerator emissions through the consumption of home-grown beef. The assessment uses a nested Monte Carlo model that separately characterizes uncertainty and variation in dose rate estimates. Uncertainty resulting from limited data on the fate and transport of TCDD are evaluated, and variations in estimated dose rates in the exposed population that result from location-specific parameters and individuals'behaviors are characterized. The analysis indicates that lifetime average daily dose rates for individuals living within 10 km of a hypothetical incinerator range over three orders of magnitude. In contrast, the uncertainty in the dose rate distribution appears to vary by less than one order of magnitude, based on the sources of uncertainty included in this analysis. Current guidance for predicting exposures from indirect exposure pathways was found to overestimate the intakes for typical and high-end individuals.     

A Probabilistic Arsenic Exposure Assessment for Children Who Contact Chromated Copper Arsenate (CCA)-Treated Playsets and Decks, Part 2: Sensitivity and Uncertainty Analyses

A probabilistic model (SHEDS-Wood) was developed to examine children's exposure and dose to chromated copper arsenate (CCA)-treated wood, as described in Part 1 of this two-part article. This Part 2 article discusses sensitivity and uncertainty analyses conducted to assess the key model inputs and areas of needed research for children's exposure to CCA-treated playsets and decks. The following types of analyses were conducted: (1) sensitivity analyses using a percentile scaling approach and multiple stepwise regression; and (2) uncertainty analyses using the bootstrap and two-stage Monte Carlo techniques. The five most important variables, based on both sensitivity and uncertainty analyses, were: wood surface residue-to-skin transfer efficiency; wood surface residue levels; fraction of hand surface area mouthed per mouthing event; average fraction of nonresidential outdoor time a child plays on/around CCA-treated public playsets; and frequency of hand washing. In general, there was a factor of 8 for the 5th and 95th percentiles and a factor of 4 for the 50th percentile in the uncertainty of predicted population dose estimates due to parameter uncertainty. Data were available for most of the key model inputs identified with sensitivity and uncertainty analyses; however, there were few or no data for some key inputs. To evaluate and improve the accuracy of model results, future measurement studies should obtain longitudinal time-activity diary information on children, spatial and temporal measurements of residue and soil concentrations on or near CCA-treated playsets and decks, and key exposure factors. Future studies should also address other sources of uncertainty in addition to parameter uncertainty, such as scenario and model uncertainty.     

Issues Related to Time Averaging of Exposure in Modeling Risks Associated with Intermittent Exposures to Lead

Typical exposures to lead often involve a mix of long-term exposures to relatively constant exposure levels (e.g., residential yard soil and indoor dust) and highly intermittent exposures at other locations (e.g., seasonal recreational visits to a park). These types of exposures can be expected to result in blood lead concentrations that vary on a temporal scale with the intermittent exposure pattern. Prediction of short-term (or seasonal) blood lead concentrations arising from highly variable intermittent exposures requires a model that can reliably simulate lead exposures and biokinetics on a temporal scale that matches that of the exposure events of interest. If exposure model averaging times (EMATs) of the model exceed the shortest exposure duration that characterizes the intermittent exposure, uncertainties will be introduced into risk estimates because the exposure concentration used as input to the model must be time averaged to account for the intermittent nature of the exposure. We have used simulation as a means of determining the potential magnitude of these uncertainties. Simulations using models having various EMATs have allowed exploration of the strengths and weaknesses of various approaches to time averaging of exposures and impact on risk estimates associated with intermittent exposures to lead in soil. The International Commission of Radiological Protection (ICRP) model of lead pharmacokinetics in humans simulates lead intakes that can vary in intensity over time spans as small as one day, allowing for the simulation of intermittent exposures to lead as a series of discrete daily exposure events. The ICRP model was used to compare the outcomes (blood lead concentration) of various time-averaging adjustments for approximating the time-averaged intake of lead associated with various intermittent exposure patterns. Results of these analyses suggest that standard approaches to time averaging (e.g., U.S. EPA) that estimate the long-term daily exposure concentration can, in some cases, result in substantial underprediction of short-term variations in blood lead concentrations when used in models that operate with EMATs exceeding the shortest exposure duration that characterizes the intermittent exposure. Alternative time-averaging approaches recommended for use in lead risk assessment more reliably predict short-term periodic (e.g., seasonal) elevations in blood lead concentration that might result from intermittent exposures. In general, risk estimates will be improved by simulation on shorter time scales that more closely approximate the actual temporal dynamics of the exposure.     

Consumption of Bovine Spongiform Encephalopathy (BSE) Contaminated Beef and the Risk of Variant Creutzfeldt‐Jakob Disease

To date, the variant Creutzfeldt‐Jakob disease (vCJD) risk assessments that have been performed have primarily focused on predicting future vCJD cases in the United Kingdom, which underwent a bovine spongiform encephalopathy (BSE) epidemic between 1980 and 1996. Surveillance of potential BSE cases was also used to assess vCJD risk, especially in other BSE‐prevalent EU countries. However, little is known about the vCJD risk for uninfected individuals who accidentally consume BSE‐contaminated meat products in or imported from a country with prevalent BSE. In this article, taking into account the biological mechanism of abnormal prion PrP^res aggregation in the brain, the probability of exposure, and the expected amount of ingested infectivity, we establish a stochastic mean exponential growth model of lifetime exposure through dietary intake. Given the findings that BSE agents behave similarly in humans and macaques, we obtained parameter estimates from experimental macaque data. We then estimated the accumulation of abnormal prions to assess lifetime risk of developing clinical signs of vCJD. Based on the observed number of vCJD cases and the estimated number of exposed individuals during the BSE epidemic period from 1980 to 1996 in the United Kingdom, an exposure threshold hypothesis is proposed. Given the age‐specific risk of infection, the hypothesis explains the observations very well from an extreme‐value distribution fitting of the estimated BSE infectivity exposure. The current BSE statistics in the United Kingdom are provided as an example.     

Exposure Reconstruction for the TCDD-Exposed NIOSH Cohort Using a Concentration- and Age-Dependent Model of Elimination

Recent studies demonstrating a concentration dependence of elimination of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) suggest that previous estimates of exposure for occupationally exposed cohorts may have underestimated actual exposure, resulting in a potential overestimate of the carcinogenic potency of TCDD in humans based on the mortality data for these cohorts. Using a database on U.S. chemical manufacturing workers potentially exposed to TCDD compiled by the National Institute for Occupational Safety and Health (NIOSH), we evaluated the impact of using a concentration- and age-dependent elimination model (CADM) (Aylward et al., 2005) on estimates of serum lipid area under the curve (AUC) for the NIOSH cohort. These data were used previously by Steenland et al. (2001) in combination with a first-order elimination model with an 8.7-year half-life to estimate cumulative serum lipid concentration (equivalent to AUC) for these workers for use in cancer dose-response assessment. Serum lipid TCDD measurements taken in 1988 for a subset of the cohort were combined with the NIOSH job exposure matrix and work histories to estimate dose rates per unit of exposure score. We evaluated the effect of choices in regression model (regression on untransformed vs. ln-transformed data and inclusion of a nonzero regression intercept) as well as the impact of choices of elimination models and parameters on estimated AUCs for the cohort. Central estimates for dose rate parameters derived from the serum-sampled subcohort were applied with the elimination models to time-specific exposure scores for the entire cohort to generate AUC estimates for all cohort members. Use of the CADM resulted in improved model fits to the serum sampling data compared to the first-order models. Dose rates varied by a factor of 50 among different combinations of elimination model, parameter sets, and regression models. Use of a CADM results in increases of up to five-fold in AUC estimates for the more highly exposed members of the cohort compared to estimates obtained using the first-order model with 8.7-year half-life. This degree of variation in the AUC estimates for this cohort would affect substantially the cancer potency estimates derived from the mortality data from this cohort. Such variability and uncertainty in the reconstructed serum lipid AUC estimates for this cohort, depending on elimination model, parameter set, and regression model, have not been described previously and are critical components in evaluating the dose-response data from the occupationally exposed populations.     

Modeling for Risk Assessment of Neurotoxic Effects

The regulation of noncancer toxicants, including neurotoxicants, has usually been based upon a reference dose (allowable daily intake). A reference dose is obtained by dividing a no-observed-effect level by uncertainty (safety) factors to account for intraspecies and interspecies sensitivities to a chemical. It is assumed that the risk at the reference dose is negligible, but no attempt generally is made to estimate the risk at the reference dose. A procedure is outlined that provides estimates of risk as a function of dose. The first step is to establish a mathematical relationship between a biological effect and the dose of a chemical. Knowledge of biological mechanisms and/or pharmacokinetics can assist in the choice of plausible mathematical models. The mathematical model provides estimates of average responses as a function of dose. Secondly, estimates of risk require selection of a distribution of individual responses about the average response given by the mathematical model. In the case of a normal or lognormal distribution, only an estimate of the standard deviation is needed. The third step is to define an adverse level for a response so that the probability (risk) of exceeding that level can be estimated as a function of dose. Because a firm response level often cannot be established at which adverse biological effects occur, it may be necessary to at least establish an abnormal response level that only a small proportion of individuals would exceed in an unexposed group. That is, if a normal range of responses can be established, then the probability (risk) of abnormal responses can be estimated. In order to illustrate this process, measures of the neurotransmitter serotonin and its metabolite 5-hydroxyindoleacetic acid in specific areas of the brain of rats and monkeys are analyzed after exposure to the neurotoxicant methylene-dioxymethamphetamine. These risk estimates are compared with risk estimates from the quantal approach in which animals are classified as either abnormal or not depending upon abnormal serotonin levels.     

A Trichloroethylene Risk Assessment Using a Monte Carlo Analysis of Parameter Uncertainty in Conjunction with Physiologically-Based Pharmacokinetic Modeling

A Monte Carlo simulation is incorporated into a risk assessment for trichloroethylene (TCE) using physiologically-based pharmacokinetic (PBPK) modeling coupled with the linearized multistage model to derive human carcinogenic risk extrapolations. The Monte Carlo technique incorporates physiological parameter variability to produce a statistically derived range of risk estimates which quantifies specific uncertainties associated with PBPK risk assessment approaches. Both inhalation and ingestion exposure routes are addressed. Simulated exposure scenarios were consistent with those used by the Environmental Protection Agency (EPA) in their TCE risk assessment. Mean values of physiological parameters were gathered from the literature for both mice (carcinogenic bioassay subjects) and for humans. Realistic physiological value distributions were assumed using existing data on variability. Mouse cancer bioassay data were correlated to total TCE metabolized and area-under-the-curve (blood concentration) trichloroacetic acid (TCA) as determined by a mouse PBPK model. These internal dose metrics were used in a linearized multistage model analysis to determine dose metric values corresponding to 10^-6 lifetime excess cancer risk. Using a human PBPK model, these metabolized doses were then extrapolated to equivalent human exposures (inhalation and ingestion). The Monte Carlo iterations with varying mouse and human physiological parameters produced a range of human exposure concentrations producing a 10^-6 risk.     

The Use of Multizone Models to Estimate an Airborne Chemical Contaminant Generation and Decay Profile: Occupational Exposures of Hairdressers to Vinyl Chloride in Hairspray During the 1960s and 1970s

Vinyl chloride (VC) was used as a propellant in a limited percentage of aerosol hairspray products in the United States from approximately 1967 to 1973. The question has arisen whether occupational exposures of hairdressers to VC-containing hairsprays in hair salons were sufficient to increase the risk for developing hepatic angiosarcoma (HAS). Transient two-zone and steady-state three-zone models were used to estimate the historical airborne concentration of VC for individual hairdressers using hairspray as well as estimated contributions from other hairdressers in the same salon. Concentrations of VC were modeled for small, medium, and large salons, as well as a representative home salon. Model inputs were determined using published literature, and variability in these inputs was also considered using Monte Carlo techniques. The 95th percentile for the daily time-weighted average exposure for small, medium, and large salons, assuming a market-share fraction of VC-containing hairspray use from the Monte Carlo analysis, was about 0.3 ppm, and for the home salon scenario was 0.1 ppm. The 95th percentile value for the cumulative lifetime exposure of the hairdressers was 2.8 ppm-years for the home salon scenario and 2.0 ppm-years for the small, medium, and large salon scenarios. If using the assumption that all hairsprays used in a salon contained VC, the 95th percentile of the theoretical lifetime cumulative dose was estimated to be 52–79 ppm-years. Estimated lifetime doses were all below the threshold dose for HAS of about 300 to 500 ppm-years reported in the published epidemiology literature.     

A Probabilistic Arsenic Exposure Assessment for Children Who Contact CCA-Treated Playsets and Decks, Part 1: Model Methodology, Variability Results, and Model Evaluation

Concerns have been raised regarding the safety of young children who may contact arsenic residues while playing on and around chromated copper arsenate (CCA)-treated wood playsets and decks. Although CCA registrants voluntarily canceled the production of treated wood for residential use in 2003, the potential for exposure from existing structures and surrounding soil still poses concerns. The EPA's Office of Research and Development developed and applied the probabilistic Stochastic Human Exposure and Dose Simulation model for wood preservatives (SHEDS-Wood) to estimate children's absorbed dose of arsenic from CCA. Skin contact with, and nondietary ingestion of, arsenic in soil and wood residues were considered for the population of children in the United States who frequently contact CCA-treated wood playsets and decks. Model analyses were conducted to assess the range in population estimates and the impact of potential mitigation strategies such as the use of sealants and hand washing after play events. The results show predicted central values for lifetime annual average daily dose values for arsenic ranging from 10(-6) to 10(-5) mg/kg/day, with predicted 95th percentiles on the order of 10(-5) mg/kg/day. There were several orders of magnitude between lower and upper percentiles. Residue ingestion via hand-to-mouth contact was determined to be the most significant exposure route for most scenarios. Results of several alternative scenarios were similar to baseline results, except for the scenario with greatly reduced residue concentrations through hypothetical wood sealant applications; in this scenario, exposures were lower, and the soil ingestion route dominated. SHEDS-Wood estimates are typically consistent with, or within the range of, other CCA exposure models.     

An Application of the Holistochastic Human Exposure Methodology to Naturally Occurring Arsenic in Bangladesh Drinking Water

The occurrence of arsenic in drinking water is an issue of considerable interest. In the case of Bangladesh, arsenic concentrations have been closely monitored since the early 1990s through an extensive sampling network. The focus of the present work is methodological. In particular, we propose the application of a holistochastic framework of human exposure to study lifetime population damage due to arsenic exposure across Bangladesh. The Bayesian Maximum Entropy theory is an important component of this framework, which possesses solid theoretical foundations and offers powerful tools to assimilate a variety of knowledge bases (physical, epidemiologic, toxicokinetic, demographic, etc.) and uncertainty sources (soft data, measurement errors, etc.). The holistochastic exposure approach leads to physically meaningful and informative spatial maps of arsenic distribution in Bangladesh drinking water. Global indicators of the adverse health effects on the population are generated, and valuable insight is gained by blending information from different scientific disciplines. The numerical results indicate an increased lifetime bladder cancer probability for the Bangladesh population due to arsenic. The health effect estimates obtained and the associated uncertainty assessments are valuable tools for a broad spectrum of end-users.     

Impact of Acquired Immunity and Dose‐Dependent Probability of Illness on Quantitative Microbial Risk Assessment

Dose‐response models in microbial risk assessment consider two steps in the process ultimately leading to illness: from exposure to (asymptomatic) infection, and from infection to (symptomatic) illness. Most data and theoretical approaches are available for the exposure‐infection step; the infection‐illness step has received less attention. Furthermore, current microbial risk assessment models do not account for acquired immunity. These limitations may lead to biased risk estimates. We consider effects of both dose dependency of the conditional probability of illness given infection, and acquired immunity to risk estimates, and demonstrate their effects in a case study on exposure to Campylobacter jejuni. To account for acquired immunity in risk estimates, an inflation factor is proposed. The inflation factor depends on the relative rates of loss of protection over exposure. The conditional probability of illness given infection is based on a previously published model, accounting for the within‐host dynamics of illness. We find that at low (average) doses, the infection‐illness model has the greatest impact on risk estimates, whereas at higher (average) doses and/or increased exposure frequencies, the acquired immunity model has the greatest impact. The proposed models are strongly nonlinear, and reducing exposure is not expected to lead to a proportional decrease in risk and, under certain conditions, may even lead to an increase in risk. The impact of different dose‐response models on risk estimates is particularly pronounced when introducing heterogeneity in the population exposure distribution.