Benzene Risk Estimation Using Radiation Equivalent Coefficients

We estimated benzene risk using a novel framework of risk assessment that employed the measurement of radiation dose equivalents to benzene metabolites and a PBPK model. The highest risks for 1 μg/m³ and 3.2 mg/m³ life time exposure of benzene estimated with a linear regression were 5.4 × 10⁻⁷ and 1.3 × 10⁻³, respectively. Even though these estimates were based on in vitro chromosome aberration test data, they were about one-sixth to one-fourteenth that from other studies and represent a fairly good estimate by using radiation equivalent coefficient as an "internal standard."     

Improving cancer dose-response characterization by using physiologically based pharmacokinetic modeling: an analysis of pooled data for acrylonitrile-induced brain tumors to assess cancer potency in the rat.

Historically, U.S. regulators have derived cancer slope factors by using applied dose and tumor response data from a single key bioassay or by averaging the cancer slope factors of several key bioassays. Recent changes in U.S. Environmental Protection Agency (EPA) guidelines for cancer risk assessment have acknowledged the value of better use of mechanistic data and better dose-response characterization. However, agency guidelines may benefit from additional considerations presented in this paper. An exploratory study was conducted by using rat brain tumor data for acrylonitrile (AN) to investigate the use of physiologically based pharmacokinetic (PBPK) modeling along with pooling of dose-response data across routes of exposure as a means for improving carcinogen risk assessment methods. In this study, two contrasting assessments were conducted for AN-induced brain tumors in the rat on the basis of (1) the EPA's approach, the dose-response relationship was characterized by using administered dose/concentration for each of the key studies assessed individually; and (2) an analysis of the pooled data, the dose-response relationship was characterized by using PBPK-derived internal dose measures for a combined database of ten bioassays. The cancer potencies predicted for AN by the contrasting assessments are remarkably different (i.e., risk-specific doses differ by as much as two to four orders of magnitude), with the pooled data assessments yielding lower values. This result suggests that current carcinogen risk assessment practices overestimate AN cancer potency. This methodology should be equally applicable to other data-rich chemicals in identifying (1) a useful dose measure, (2) an appropriate dose-response model, (3) an acceptable point of departure, and (4) an appropriate method of extrapolation from the range of observation to the range of prediction when a chemical's mode of action remains uncertain.     

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.     

Addressing Nonlinearity in the Exposure-Response Relationship for a Genotoxic Carcinogen: Cancer Potency Estimates for Ethylene Oxide

Ethylene oxide (EO) has been identified as a carcinogen in laboratory animals. Although the precise mechanism of action is not known, tumors in animals exposed to EO are presumed to result from its genotoxicity. The overall weight of evidence for carcinogenicity from a large body of epidemiological data in the published literature remains limited. There is some evidence for an association between EO exposure and lympho/hematopoietic cancer mortality. Of these cancers, the evidence provided by two large cohorts with the longest follow-up is most consistent for leukemia. Together with what is known about human leukemia and EO at the molecular level, there is a body of evidence that supports a plausible mode of action for EO as a potential leukemogen. Based on a consideration of the mode of action, the events leading from EO exposure to the development of leukemia (and therefore risk) are expected to be proportional to the square of the dose. In support of this hypothesis, a quadratic dose-response model provided the best overall fit to the epidemiology data in the range of observation. Cancer dose-response assessments based on human and animal data are presented using three different assumptions for extrapolating to low doses: (1) risk is linearly proportionate to dose; (2) there is no appreciable risk at low doses (margin-of-exposure or reference dose approach); and (3) risk below the point of departure continues to be proportionate to the square of the dose. The weight of evidence for EO supports the use of a nonlinear assessment. Therefore, exposures to concentrations below 37 microg/m3 are not likely to pose an appreciable risk of leukemia in human populations. However, if quantitative estimates of risk at low doses are desired and the mode of action for EO is considered, these risks are best quantified using the quadratic estimates of cancer potency, which are approximately 3.2- to 32-fold lower, using alternative points of departure, than the linear estimates of cancer potency for EO. An approach is described for linking the selection of an appropriate point of departure to the confidence in the proposed mode of action. Despite high confidence in the proposed mode of action, a small linear component for the dose-response relationship at low concentrations cannot be ruled out conclusively. Accordingly, a unit risk value of 4.5 x 10(-8) (microg/m3)(-1) was derived for EO, with a range of unit risk values of 1.4 x 10(-8) to 1.4 x 10(-7) (microg/m3)(-1) reflecting the uncertainty associated with a theoretical linear term at low concentrations.     

Using Physiologically-Based Pharmacokinetic Modeling to Address Nonlinear Kinetics and Changes in Rodent Physiology and Metabolism Due to Aging and Adaptation in Deriving Reference Values for Propylene Glycol Methyl Ether and Propylene Glycol Methyl Ether Acetate

Reference values, including an oral reference dose (RfD) and an inhalation reference concentration (RfC), were derived for propylene glycol methyl ether (PGME), and an oral RfD was derived for its acetate (PGMEA). These values were based on transient sedation observed in F344 rats and B6C3F1 mice during a two‐year inhalation study. The dose‐response relationship for sedation was characterized using internal dose measures as predicted by a physiologically‐based pharmacokinetic (PBPK) model for PGME and its acetate. PBPK modeling was used to account for changes in rodent physiology and metabolism due to aging and adaptation, based on data collected during Weeks 1, 2, 26, 52, and 78 of a chronic inhalation study. The peak concentration of PGME in richly perfused tissues (i.e., brain) was selected as the most appropriate internal dose measure based on a consideration of the mode of action for sedation and similarities in tissue partitioning between brain and other richly perfused tissues. Internal doses (peak tissue concentrations of PGME) were designated as either no‐observed‐adverse‐effect levels (NOAELs) or lowest‐observed‐adverse‐effect levels (LOAELs) based on the presence or the absence of sedation at each time point, species, and sex in the two‐year study. Distributions of the NOAEL and LOAEL values expressed in terms of internal dose were characterized using an arithmetic mean and standard deviation, with the mean internal NOAEL serving as the basis for the reference values, which was then divided by appropriate uncertainty factors. Where data were permitting, chemical‐specific adjustment factors were derived to replace default uncertainty factor values of 10. Nonlinear kinetics, which was predicted by the model in all species at PGME concentrations exceeding 100 ppm, complicate interspecies, and low‐dose extrapolations. To address this complication, reference values were derived using two approaches that differ with respect to the order in which these extrapolations were performed: (1) default approach of interspecies extrapolation to determine the human equivalent concentration (PBPK modeling) followed by uncertainty factor application, and (2) uncertainty factor application followed by interspecies extrapolation (PBPK modeling). The resulting reference values for these two approaches are substantially different, with values from the latter approach being seven‐fold higher than those from the former approach. Such a striking difference between the two approaches reveals an underlying issue that has received little attention in the literature regarding the application of uncertainty factors and interspecies extrapolations to compounds where saturable kinetics occur in the range of the NOAEL. Until such discussions have taken place, reference values based on the former approach are recommended for risk assessments involving human exposures to PGME and PGMEA.     

Symmetry, Identifiability, and Prediction Uncertainties in Multistage Clonal Expansion (MSCE) Models of Carcinogenesis

Many models of exposure-related carcinogenesis, including traditional linearized multistage models and more recent two-stage clonal expansion (TSCE) models, belong to a family of models in which cells progress between successive stages-possibly undergoing proliferation at some stages-at rates that may depend (usually linearly) on biologically effective doses. Biologically effective doses, in turn, may depend nonlinearly on administered doses, due to PBPK nonlinearities. This article provides an exact mathematical analysis of the expected number of cells in the last ("malignant") stage of such a "multistage clonal expansion" (MSCE) model as a function of dose rate and age. The solution displays symmetries such that several distinct sets of parameter values provide identical fits to all epidemiological data, make identical predictions about the effects on risk of changes in exposure levels or timing, and yet make significantly different predictions about the effects on risk of changes in the composition of exposure that affect the pharmacodynamic dose-response relation. Several different predictions for the effects of such an intervention (such as reducing carcinogenic constituents of an exposure) that acts on only one or a few stages of the carcinogenic process may be equally consistent with all preintervention epidemiological data. This is an example of nonunique identifiability of model parameters and predictions from data. The new results on nonunique model identifiability presented here show that the effects of an intervention on changing age-specific cancer risks in an MSCE model can be either large or small, but that which is the case cannot be predicted from preintervention epidemiological data and knowledge of biological effects of the intervention alone. Rather, biological data that identify which rate parameters hold for which specific stages are required to obtain unambiguous predictions. From epidemiological data alone, only a set of equally likely alternative predictions can be made for the effects on risk of such interventions.     

Model Averaging Using the Kullback Information Criterion in Estimating Effective Doses for Microbial Infection and Illness

Since the National Food Safety Initiative of 1997, risk assessment has been an important issue in food safety areas. Microbial risk assessment is a systematic process for describing and quantifying a potential to cause adverse health effects associated with exposure to microorganisms. Various dose-response models for estimating microbial risks have been investigated. We have considered four two-parameter models and four three-parameter models in order to evaluate variability among the models for microbial risk assessment using infectivity and illness data from studies with human volunteers exposed to a variety of microbial pathogens. Model variability is measured in terms of estimated ED01s and ED10s, with the view that these effective dose levels correspond to the lower and upper limits of the 1% to 10% risk range generally recommended for establishing benchmark doses in risk assessment. Parameters of the statistical models are estimated using the maximum likelihood method. In this article a weighted average of effective dose estimates from eight two- and three-parameter dose-response models, with weights determined by the Kullback information criterion, is proposed to address model uncertainties in microbial risk assessment. The proposed procedures for incorporating model uncertainties and making inferences are illustrated with human infection/illness dose-response data sets.     

Comparing Experimental Designs for Benchmark Dose Calculations for Continuous Endpoints

The BMD (benchmark dose) method that is used in risk assessment of chemical compounds was introduced by Crump (1984) and is based on dose-response modeling. To take uncertainty in the data and model fitting into account, the lower confidence bound of the BMD estimate (BMDL) is suggested to be used as a point of departure in health risk assessments. In this article, we study how to design optimum experiments for applying the BMD method for continuous data. We exemplify our approach by considering the class of Hill models. The main aim is to study whether an increased number of dose groups and at the same time a decreased number of animals in each dose group improves conditions for estimating the benchmark dose. Since Hill models are nonlinear, the optimum design depends on the values of the unknown parameters. That is why we consider Bayesian designs and assume that the parameter vector has a prior distribution. A natural design criterion is to minimize the expected variance of the BMD estimator. We present an example where we calculate the value of the design criterion for several designs and try to find out how the number of dose groups, the number of animals in the dose groups, and the choice of doses affects this value for different Hill curves. It follows from our calculations that to avoid the risk of unfavorable dose placements, it is good to use designs with more than four dose groups. We can also conclude that any additional information about the expected dose-response curve, e.g., information obtained from studies made in the past, should be taken into account when planning a study because it can improve the design.     

A Calibrated Human PBPK Model for Benzene Inhalation with Urinary Bladder and Bone Marrow Compartments

A physiologically‐based pharmacokinetic (PBPK) model of benzene inhalation based on a recent mouse model was adapted to include bone marrow (target organ) and urinary bladder compartments. Empirical data on human liver microsomal protein levels and linked CYP2E1 activities were incorporated into the model, and metabolite‐specific conversion rate parameters were estimated by fitting to human biomonitoring data and adjusting for background levels of urinary metabolites. Human studies of benzene levels in blood and breath, and phenol levels in urine were used to validate the rate of human conversion of benzene to benzene oxide, and urinary benzene metabolites from Chinese benzene worker populations provided model validation for rates of human conversion of benzene to muconic acid (MA) and phenylmercapturic acid (PMA), phenol (PH), catechol (CA), hydroquinone (HQ), and benzenetriol (BT). The calibrated human model reveals that while liver microsomal protein and CYP2E1 activities are lower on average in humans compared to mice, the mouse also shows far lower rates of benzene conversion to MA and PMA, and far higher conversion of benzene to BO/PH, and of BO/PH to CA, HQ, and BT. The model also differed substantially from existing human PBPK models with respect to several metabolic rate parameters of importance to interpreting benzene metabolism and health risks in human populations associated with bone marrow doses. The model provides a new methodological paradigm focused on integrating linked human liver metabolism data and calibration using biomonitoring data, thus allowing for model uncertainty analysis and more rigorous validation.     

Uncertainty in Cancer Risk Estimates

Several existing databases compiled by Gold et al.^(1–3) for carcinogenesis bioassays are examined to obtain estimates of the reproducibility of cancer rates across experiments, strains, and rodent species. A measure of carcinogenic potency is given by the TD₅₀ (daily dose that causes a tumor type in 50% of the exposed animals that otherwise would not develop the tumor in a standard lifetime). The lognormal distribution can be used to model the uncertainty of the estimates of potency (TD₅₀) and the ratio of TD₅₀'s between two species. For near-replicate bioassays, approximately 95% of the TD₅₀'s are estimated to be within a factor of 4 of the mean. Between strains, about 95% of the TD₅₀'s are estimated to be within a factor of 11 of their mean, and the pure genetic component of variability is accounted for by a factor of 6.8. Between rats and mice, about 95% of the TD₅₀'s are estimated to be within a factor of 32 of the mean, while between humans and experimental animals the factor is 110 for 20 chemicals reported by Allen et al.⁽⁴⁾ The common practice of basing cancer risk estimates on the most sensitive rodent species-strain-sex and using interspecies dose scaling based on body surface area appears to overestimate cancer rates for these 20 human carcinogens by about one order of magnitude on the average. Hence, for chemicals where the dose-response is nearly linear below experimental doses, cancer risk estimates based on animal data are not necessarily conservative and may range from a factor of 10 too low for human carcinogens up to a factor of 1000 too high for approximately 95% of the chemicals tested to date. These limits may need to be modified for specific chemicals where additional mechanistic or pharmacokinetic information may suggest alterations or where particularly sensitive subpopu-lations may be exposed. Supralinearity could lead to anticonservative estimates of cancer risk. Underestimating cancer risk by a specific factor has a much larger impact on the actual number of cancer cases than overestimates of smaller risks by the same factor. This paper does not address the uncertainties in high to low dose extrapolation. If the dose-response is sufficiently nonlinear at low doses to produce cancer risks near zero, then low-dose risk estimates based on linear extrapolation are likely to overestimate risk and the limits of uncertainty cannot be established.     

Linking Indoor Air and Pharmacokinetic Models to Assess Tetrachloroethylene Risk

Physiologically based pharmacokinetic (PBPK) models describing the uptake, metabolism, and excretion of xenobiotic compounds are now proposed for use in regulatory health-risk assessments. In this study we investigate the extent of PCE metabolism arising from domestic respiratory exposure to tetrachloroethylene (PCE) from ground water, as predicted using a PBPK model. Indoor exposure patterns we use as input to the PBPK model are realistic ones generated from a three-compartment model describing volatilization of PCE from domestic water into household air. Values we use for the metabolic parameters of the PBPK model are estimated from data on urinary metabolites in workers exposed to PCE. It is shown that for respiratory PCE exposure due to typical levels of PCE in ground water, use of time-weighted average air concentrations with a steady-state PBPK model yields estimates of total metabolized PCE similar to those obtained using completely dynamic modeling, despite considerable uncertainty in key exposure- and metabolic-model parameters. These findings suggest that, for PCE, risk estimation taking pharmacokinetics into account may be accomplished using a simple analytic approach.     

Quantitative Cancer Risk Estimation for Formaldehyde

Of primary concern are irreversible effects, such as cancer induction, that formaldehyde exposure could have on human health. Dose-response data from human exposure situations would provide the most solid foundation for risk assessment, avoiding problematic extrapolations from the health effects seen in nonhuman species. However, epidemiologic studies of human formaldehyde exposure have provided little definitive information regarding dose-response. Reliance must consequently be placed on laboratory animal evidence. An impressive array of data points to significantly nonlinear relationships between rodent tumor incidence and administered dose, and between target tissue dose and administered dose (the latter for both rodents and Rhesus monkeys) following exposure to formaldehyde by inhalation. Disproportionately less formaldehyde binds covalently to the DNA of nasal respiratory epithelium at low than at high airborne concentrations. Use of this internal measure of delivered dose in analyses of rodent bioassay nasal tumor response yields multistage model estimates of low-dose risk, both point and upper bound, that are lower than equivalent estimates based upon airborne formaldehyde concentration. In addition, risk estimates obtained for Rhesus monkeys appear at least 10-fold lower than corresponding estimates for identically exposed Fischer-344 rats.     

Cancer Dose-Response Modeling of Epidemiological Data on Worker Exposures to Aldrin and Dieldrin

The paper applies classical statistical principles to yield new tools for risk assessment and makes new use of epidemiological data for human risk assessment. An extensive clinical and epidemiological study of workers engaged in the manufacturing and formulation of aldrin and dieldrin provides occupational hygiene and biological monitoring data on individual exposures over the years of employment and provides unusually accurate measures of individual lifetime average daily doses. In the cancer dose-response modeling, each worker is treated as a separate experimental unit with his own unique dose. Maximum likelihood estimates of added cancer risk are calculated for multistage, multistage-Weibull, and proportional hazards models. Distributional characterizations of added cancer risk are based on bootstrap and relative likelihood techniques. The cancer mortality data on these male workers suggest that low-dose exposures to aldrin and dieldrin do not significantly increase human cancer risk and may even decrease the human hazard rate for all types of cancer combined at low doses (e.g., 1 g/kg/day). The apparent hormetic effect in the best fitting dose-response models for this data set is statistically significant. The decrease in cancer risk at low doses of aldrin and dieldrin is in sharp contrast to the U.S. Environmental Protection Agency's upper bound on cancer potency based on mouse liver tumors. The EPA's upper bound implies that lifetime average daily doses of 0.0000625 and 0.00625 g/kg body weight/day would correspond to increased cancer risks of 0.000001 and 0.0001, respectively. However, the best estimate from the Pernis epidemiological data is that there is no increase in cancer risk in these workers at these doses or even at doses as large as 2 g/kg/day.     

Improving Cancer Dose–Response Characterization by Using Physiologically Based Pharmacokinetic Modeling: An Analysis of Pooled Data for Acrylonitrile-Induced Brain Tumors to Assess Cancer Potency in the Rat

Historically, U.S. regulators have derived cancer slope factors by using applied dose and tumor response data from a single key bioassay or by averaging the cancer slope factors of several key bioassays. Recent changes in U.S. Environmental Protection Agency (EPA) guidelines for cancer risk assessment have acknowledged the value of better use of mechanistic data and better dose–response characterization. However, agency guidelines may benefit from additional considerations presented in this paper. An exploratory study was conducted by using rat brain tumor data for acrylonitrile (AN) to investigate the use of physiologically based pharmacokinetic (PBPK) modeling along with pooling of dose–response data across routes of exposure as a means for improving carcinogen risk assessment methods. In this study, two contrasting assessments were conducted for AN-induced brain tumors in the rat on the basis of (1) the EPA's approach, the dose–response relationship was characterized by using administered dose/concentration for each of the key studies assessed individually; and (2) an analysis of the pooled data, the dose–response relationship was characterized by using PBPK-derived internal dose measures for a combined database of ten bioassays. The cancer potencies predicted for AN by the contrasting assessments are remarkably different (i.e., risk-specific doses differ by as much as two to four orders of magnitude), with the pooled data assessments yielding lower values. This result suggests that current carcinogen risk assessment practices overestimate AN cancer potency. This methodology should be equally applicable to other data-rich chemicals in identifying (1) a useful dose measure, (2) an appropriate dose–response model, (3) an acceptable point of departure, and (4) an appropriate method of extrapolation from the range of observation to the range of prediction when a chemical's mode of action remains uncertain.     

Application of a Biologically-Based RFD Estimation Method to Tetrachlorodibenzo-P-Dioxin (TCDD) Mediated Immune Suppression and Enzyme Induction

The current methods for a reference dose (RfD) determination can be enhanced through the use of biologically-based dose-response analysis. Methods developed here utilizes information from tetrachlorodibenzo- p -dioxin (TCDD) to focus on noncancer endpoints, specifically TCDD mediated immune system alterations and enzyme induction. Dose-response analysis, using the Sigmoid-Emax (EMAX) function, is applied to multiple studies to determine consistency of response. Through the use of multiple studies and statistical comparison of parameter estimates, it was demonstrated that the slope estimates across studies were very consistent. This adds confidence to the subsequent effect dose estimates. This study also compares traditional methods of risk assessment such as the NOAEL/safety factor to a modified benchmark dose approach which is introduced here. Confidence in the estimation of an effect dose (ED₁₀) was improved through the use of multiple datasets. This is key to adding confidence to the benchmark dose estimates. In addition, the Sigmoid-Emax function when applied to dose-response data using nonlinear regression analysis provides a significantly improved fit to data increasing confidence in parameter estimates which subsequently improve effect dose estimates.     

Estimation of the Leukemia Risk in Human Populations Exposed to Benzene from Tobacco Smoke Using Epidemiological Data

Several epidemiological studies have demonstrated an association between occupational benzene exposure and increased leukemia risk, in particular acute myeloid leukemia (AML). However, there is still uncertainty as to the risk to the general population from exposure to lower environmental levels of benzene. To estimate the excess risk of leukemia from low‐dose benzene exposure, various methods for incorporating epidemiological data in quantitative risk assessment were utilized. Tobacco smoke was identified as one of the main potential sources of benzene exposure and was the focus of this exposure assessment, allowing further investigation of the role of benzene in smoking‐induced leukemia. Potency estimates for benzene were generated from individual occupational studies and meta‐analysis data, and an exposure assessment for two smoking subgroups (light and heavy smokers) carried out. Subsequently, various techniques, including life‐table analysis, were then used to evaluate both the excess lifetime risk and the contribution of benzene to smoking‐induced leukemia and AML. The excess lifetime risk for smokers was estimated at between two and six additional leukemia deaths in 10,000 and one to three additional AML deaths in 10,000. The contribution of benzene to smoking‐induced leukemia was estimated at between 9% and 24% (U_pperCL 14–31%). For AML this contribution was estimated as 11–30% (U_pperCL 22–60%). From the assessments carried out here, it appears there is an increased risk of leukemia from low‐level exposure to benzene and that benzene may contribute up to a third of smoking‐induced leukemia. Comparable results from using methods with varying degrees of complexity were generated.     

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.     

Reproducibility of the Dose-Response Curve of Steroid-Induced Cleft Palate in Mice

Pregnant CD-1 mice were exposed to cortisone acetate at doses ranging from 20 to 100 mg/kg/ day on days 10-13 by oral and intramuscular routes. Multiple replicate assays were conducted under identical conditions to assess the reproducibility of the dose–response curve for cleft palate. The data were fitted to the probit, logistic, multistage or Armitage-Doll, and Weibull dose-response model separately for each route of exposure. The curves were then tested for parallel slopes (probit and logistic models) or coincidence of model parameters (multistage and Weibull models). The 19 replicate experiments had a wide range of slope estimates, wider for the oral than for the intramuscular experiments. For all models and both routes of exposure the null hypothesis of equality of slopes was rejected at a significant level of p < 0.001. For the intramuscular group of replicates, rejection of slope equality could in part be explained by not maintaining a standard dosing regime. The rejection of equivalence of dose-response curves from replicate studies showed that it is difficult to reproduce dose-response data of a single study within the limits defined by the dose-response model. This has important consequences for quantitative risk assessment, public health measures, or development of mechanistic theories which are typically based on a single animal bioassay.     

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.     

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.