Peak Exposures in Epidemiologic Studies and Cancer Risks: Considerations for Regulatory Risk Assessment

We review approaches for characterizing “peak” exposures in epidemiologic studies and methods for incorporating peak exposure metrics in dose–response assessments that contribute to risk assessment. The focus was on potential etiologic relations between environmental chemical exposures and cancer risks. We searched the epidemiologic literature on environmental chemicals classified as carcinogens in which cancer risks were described in relation to “peak” exposures. These articles were evaluated to identify some of the challenges associated with defining and describing cancer risks in relation to peak exposures. We found that definitions of peak exposure varied considerably across studies. Of nine chemical agents included in our review of peak exposure, six had epidemiologic data used by the U.S. Environmental Protection Agency (US EPA) in dose–response assessments to derive inhalation unit risk values. These were benzene, formaldehyde, styrene, trichloroethylene, acrylonitrile, and ethylene oxide. All derived unit risks relied on cumulative exposure for dose–response estimation and none, to our knowledge, considered peak exposure metrics. This is not surprising, given the historical linear no‐threshold default model (generally based on cumulative exposure) used in regulatory risk assessments. With newly proposed US EPA rule language, fuller consideration of alternative exposure and dose–response metrics will be supported. “Peak” exposure has not been consistently defined and rarely has been evaluated in epidemiologic studies of cancer risks. We recommend developing uniform definitions of “peak” exposure to facilitate fuller evaluation of dose response for environmental chemicals and cancer risks, especially where mechanistic understanding indicates that the dose response is unlikely linear and that short‐term high‐intensity exposures increase risk.     

Quantifying Cancer Risk from Radiation

Complex statistical models fitted to data from studies of atomic bomb survivors are used to estimate the human health effects of ionizing radiation exposures. We describe and illustrate an approach to estimate population risks from ionizing radiation exposure that relaxes many assumptions about radiation‐related mortality. The approach draws on developments in methods for causal inference. The results offer a different way to quantify radiation's effects and show that conventional estimates of the population burden of excess cancer at high radiation doses are driven strongly by projecting outside the range of current data. Summary results obtained using the proposed approach are similar in magnitude to those obtained using conventional methods, although estimates of radiation‐related excess cancers differ for many age, sex, and dose groups. At low doses relevant to typical exposures, the strength of evidence in data is surprisingly weak. Statements regarding human health effects at low doses rely strongly on the use of modeling assumptions.     

Methodology for Determining the Appropriateness of a Linear Dose‐Response Function

Microbial food safety risk assessment models can often at times be simplified by eliminating the need to integrate a complex dose‐response relationship across a distribution of exposure doses. This is possible if exposure pathways lead to pathogens at exposure that consistently have a small probability of causing illness. In this situation, the probability of illness will follow an approximately linear function of dose. Consequently, the predicted probability of illness per serving across all exposures is linear with respect to the expected value of dose. The majority of dose‐response functions are approximately linear when the dose is low. Nevertheless, what constitutes “low” is dependent on the parameters of the dose‐response function for a particular pathogen. In this study, a method is proposed to determine an upper bound of the exposure distribution for which the use of a linear dose‐response function is acceptable. If this upper bound is substantially larger than the expected value of exposure doses, then a linear approximation for probability of illness is reasonable. If conditions are appropriate for using the linear dose‐response approximation, for example, the expected value for exposure doses is two to three logs₁₀ smaller than the upper bound of the linear portion of the dose‐response function, then predicting the risk‐reducing effectiveness of a proposed policy is trivial. Simple examples illustrate how this approximation can be used to inform policy decisions and improve an analyst's understanding of risk.     

Characterizing Risk for Cumulative Risk Assessments

In assessing environmental health risks, the risk characterization step synthesizes information gathered in evaluating exposures to stressors together with dose–response relationships, characteristics of the exposed population, and external environmental conditions. This article summarizes key steps of a cumulative risk assessment (CRA) followed by a discussion of considerations for characterizing cumulative risks. Cumulative risk characterizations differ considerably from single chemical‐ or single source‐based risk characterization. CRAs typically focus on a specific population instead of a pollutant or pollutant source and should include an evaluation of all relevant sources contributing to the exposures in the population and other factors that influence dose–response relationships. Second, CRAs may include influential environmental and population‐specific conditions, involving multiple chemical and nonchemical stressors. Third, a CRA could examine multiple health effects, reflecting joint toxicity and the potential for toxicological interactions. Fourth, the complexities often necessitate simplifying methods, including judgment‐based and semi‐quantitative indices that collapse disparate data into numerical scores. Fifth, because of the higher dimensionality and potentially large number of interactions, information needed to quantify risk is typically incomplete, necessitating an uncertainty analysis. Three approaches that could be used for characterizing risks in a CRA are presented: the multiroute hazard index, stressor grouping by exposure and toxicity, and indices for screening multiple factors and conditions. Other key roles of the risk characterization in CRAs are also described, mainly the translational aspect of including a characterization summary for lay readers (in addition to the technical analysis), and placing the results in the context of the likely risk‐based decisions.     

Harnessing the Theoretical Foundations of the Exponential and Beta‐Poisson Dose‐Response Models to Quantify Parameter Uncertainty Using Markov Chain Monte Carlo

Dose‐response models are the essential link between exposure assessment and computed risk values in quantitative microbial risk assessment, yet the uncertainty that is inherent to computed risks because the dose‐response model parameters are estimated using limited epidemiological data is rarely quantified. Second‐order risk characterization approaches incorporating uncertainty in dose‐response model parameters can provide more complete information to decisionmakers by separating variability and uncertainty to quantify the uncertainty in computed risks. Therefore, the objective of this work is to develop procedures to sample from posterior distributions describing uncertainty in the parameters of exponential and beta‐Poisson dose‐response models using Bayes's theorem and Markov Chain Monte Carlo (in OpenBUGS). The theoretical origins of the beta‐Poisson dose‐response model are used to identify a decomposed version of the model that enables Bayesian analysis without the need to evaluate Kummer confluent hypergeometric functions. Herein, it is also established that the beta distribution in the beta‐Poisson dose‐response model cannot address variation among individual pathogens, criteria to validate use of the conventional approximation to the beta‐Poisson model are proposed, and simple algorithms to evaluate actual beta‐Poisson probabilities of infection are investigated. The developed MCMC procedures are applied to analysis of a case study data set, and it is demonstrated that an important region of the posterior distribution of the beta‐Poisson dose‐response model parameters is attributable to the absence of low‐dose data. This region includes beta‐Poisson models for which the conventional approximation is especially invalid and in which many beta distributions have an extreme shape with questionable plausibility.     

Dose‐Response Modeling for Inhalational Anthrax in Rabbits Following Single or Multiple Exposures

There is a need to advance our ability to characterize the risk of inhalational anthrax following a low‐dose exposure. The exposure scenario most often considered is a single exposure that occurs during an attack. However, long‐term daily low‐dose exposures also represent a realistic exposure scenario, such as what may be encountered by people occupying areas for longer periods. Given this, the objective of the current work was to model two rabbit inhalational anthrax dose‐response data sets. One data set was from single exposures to aerosolized Bacillus anthracis Ames spores. The second data set exposed rabbits repeatedly to aerosols of B. anthracis Ames spores. For the multiple exposure data the cumulative dose (i.e., the sum of the individual daily doses) was used for the model. Lethality was the response for both. Modeling was performed using Benchmark Dose Software evaluating six models: logprobit, loglogistic, Weibull, exponential, gamma, and dichotomous‐Hill. All models produced acceptable fits to either data set. The exponential model was identified as the best fitting model for both data sets. Statistical tests suggested there was no significant difference between the single exposure exponential model results and the multiple exposure exponential model results, which suggests the risk of disease is similar between the two data sets. The dose expected to cause 10% lethality was 15,600 inhaled spores and 18,200 inhaled spores for the single exposure and multiple exposure exponential dose‐response model, respectively, and the 95% lower confidence intervals were 9,800 inhaled spores and 9,200 inhaled spores, respectively.     

Model Uncertainty and Model Averaging in the Estimation of Infectious Doses for Microbial Pathogens

Food‐borne infection is caused by intake of foods or beverages contaminated with microbial pathogens. Dose‐response modeling is used to estimate exposure levels of pathogens associated with specific risks of infection or illness. When a single dose‐response model is used and confidence limits on infectious doses are calculated, only data uncertainty is captured. We propose a method to estimate the lower confidence limit on an infectious dose by including model uncertainty and separating it from data uncertainty. The infectious dose is estimated by a weighted average of effective dose estimates from a set of dose‐response models via a Kullback information criterion. The confidence interval for the infectious dose is constructed by the delta method, where data uncertainty is addressed by a bootstrap method. To evaluate the actual coverage probabilities of the lower confidence limit, a Monte Carlo simulation study is conducted under sublinear, linear, and superlinear dose‐response shapes that can be commonly found in real data sets. Our model‐averaging method achieves coverage close to nominal in almost all cases, thus providing a useful and efficient tool for accurate calculation of lower confidence limits on infectious doses.     

Benchmark Dose Analysis via Nonparametric Regression Modeling

Estimation of benchmark doses (BMDs) in quantitative risk assessment traditionally is based upon parametric dose‐response modeling. It is a well‐known concern, however, that if the chosen parametric model is uncertain and/or misspecified, inaccurate and possibly unsafe low‐dose inferences can result. We describe a nonparametric approach for estimating BMDs with quantal‐response data based on an isotonic regression method, and also study use of corresponding, nonparametric, bootstrap‐based confidence limits for the BMD. We explore the confidence limits’ small‐sample properties via a simulation study, and illustrate the calculations with an example from cancer risk assessment. It is seen that this nonparametric approach can provide a useful alternative for BMD estimation when faced with the problem of parametric model uncertainty.     

Building a Human Health Risk Assessment Ontology (RsO): A Proposed Framework

Over the last decade the health and environmental research communities have made significant progress in collecting and improving access to genomic, toxicology, exposure, health, and disease data useful to health risk assessment. One of the barriers to applying these growing volumes of information in fields such as risk assessment is the lack of informatics tools to organize, curate, and evaluate thousands of journal publications and hundreds of databases to provide new insights on relationships among exposure, hazard, and disease burden. Many fields are developing ontologies as a way of organizing and analyzing large amounts of complex information from multiple scientific disciplines. Ontologies include a vocabulary of terms and concepts with defined logical relationships to each other. Building from the recently published exposure ontology and other relevant health and environmental ontologies, this article proposes an ontology for health risk assessment (RsO) that provides a structural framework for organizing risk assessment information and methods. The RsO is anchored by eight major concepts that were either identified by exploratory curations of the risk literature or the exposure‐ontology working group as key for describing the risk assessment domain. These concepts are: (1) stressor, (2) receptor, (3) outcome, (4) exposure event, (5) dose‐response approach, (6) dose‐response metric, (7) uncertainty, and (8) measure of risk. We illustrate the utility of these concepts for the RsO with example curations of published risk assessments for ionizing radiation, arsenic in drinking water, and persistent pollutants in salmon.     

Annual Burden of Occupationally‐Acquired Influenza Infections in Hospitals and Emergency Departments in the United States

Infections among health‐care personnel (HCP) occur as a result of providing care to patients with infectious diseases, but surveillance is limited to a few diseases. The objective of this study is to determine the annual number of influenza infections acquired by HCP as a result of occupational exposures to influenza patients in hospitals and emergency departments (EDs) in the United States. A risk analysis approach was taken. A compartmental model was used to estimate the influenza dose received in a single exposure, and a dose–response function applied to calculate the probability of infection. A three‐step algorithm tabulated the total number of influenza infections based on: the total number of occupational exposures (tabulated in previous work), the total number of HCP with occupational exposures, and the probability of infection in an occupational exposure. Estimated influenza infections were highly dependent upon the dose–response function. Given current compliance with infection control precautions, we estimated 151,300 and 34,150 influenza infections annually with two dose–response functions (annual incidence proportions of 9.3% and 2.1%, respectively). Greater reductions in infectious were achieved by full compliance with vaccination and IC precautions than with patient isolation. The burden of occupationally‐acquired influenza among HCP in hospitals and EDs in the United States is not trivial, and can be reduced through improved compliance with vaccination and preventive measures, including engineering and administrative controls.     

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.     

Assessing human health response in life cycle assessment using ED10s and DALYs: part 2--Noncancer effects.

In Part 1 of this article we developed an approach for the calculation of cancer effect measures for life cycle assessment (LCA). In this article, we propose and evaluate the method for the screening of noncancer toxicological health effects. This approach draws on the noncancer health risk assessment concept of benchmark dose, while noting important differences with regulatory applications in the objectives of an LCA study. We adopt the centraltendency estimate of the toxicological effect dose inducing a 10% response over background, ED10, to provide a consistent point of departure for default linear low-dose response estimates (betaED10). This explicit estimation of low-dose risks, while necessary in LCA, is in marked contrast to many traditional procedures for noncancer assessments. For pragmatic reasons, mechanistic thresholds and nonlinear low-dose response curves were not implemented in the presented framework. In essence, for the comparative needs of LCA, we propose that one initially screens alternative activities or products on the degree to which the associated chemical emissions erode their margins of exposure, which may or may not be manifested as increases in disease incidence. We illustrate the method here by deriving the betaED10 slope factors from bioassay data for 12 chemicals and outline some of the possibilities for extrapolation from other more readily available measures, such as the no observable adverse effect levels (NOAEL), avoiding uncertainty factors that lead to inconsistent degrees of conservatism from chemical to chemical. These extrapolations facilitated the initial calculation of slope factors for an additional 403 compounds; ranging from 10(-6) to 10(3) (risk per mg/kg-day dose). The potential consequences of the effects are taken into account in a preliminary approach by combining the betaED10 with the severity measure disability adjusted life years (DALY), providing a screening-level estimate of the potential consequences associated with exposures, integrated over time and space, to a given mass of chemical released into the environment for use in LCA.     

Consideration of Background Exposures in the Management of Hazardous Waste Sites: A New Approach to Risk Assessment

The current approach to health risk assessment of toxic waste sites in the U.S. may lead to considerable expenditure of resources without any meaningful reduction in population exposure. Risk assessment methods used generally ignore background exposures and consider only incremental risk estimates for maximally exposed individuals. Such risk estimates do not address true public health risks to which background exposures also contribute. The purpose of this paper is to recommend a new approach to risk assessment and risk management concerning toxic waste sites. Under this new approach, which we have called public health risk assessment, chemical substances would be classified into a level of concern based on the potential health risks associated with typical national and regional background exposures. Site assessment would then be based on the level of concern for the particular pollutants involved and the potential contribution of site contaminants to typical background human exposures. While various problems can be foreseen with this approach, the key advantage is that resources would be allocated to reduce the most important sources of human exposure, and site remediation decisions could be simplified by focussing on exposure assessment rather than questionable risk extrapolations.     

Fractional Poisson—A Simple Dose‐Response Model for Human Norovirus

This study utilizes old and new Norovirus (NoV) human challenge data to model the dose‐response relationship for human NoV infection. The combined data set is used to update estimates from a previously published beta‐Poisson dose‐response model that includes parameters for virus aggregation and for a beta‐distribution that describes variable susceptibility among hosts. The quality of the beta‐Poisson model is examined and a simpler model is proposed. The new model (fractional Poisson) characterizes hosts as either perfectly susceptible or perfectly immune, requiring a single parameter (the fraction of perfectly susceptible hosts) in place of the two‐parameter beta‐distribution. A second parameter is included to account for virus aggregation in the same fashion as it is added to the beta‐Poisson model. Infection probability is simply the product of the probability of nonzero exposure (at least one virus or aggregate is ingested) and the fraction of susceptible hosts. The model is computationally simple and appears to be well suited to the data from the NoV human challenge studies. The model's deviance is similar to that of the beta‐Poisson, but with one parameter, rather than two. As a result, the Akaike information criterion favors the fractional Poisson over the beta‐Poisson model. At low, environmentally relevant exposure levels (<100), estimation error is small for the fractional Poisson model; however, caution is advised because no subjects were challenged at such a low dose. New low‐dose data would be of great value to further clarify the NoV dose‐response relationship and to support improved risk assessment for environmentally relevant exposures.     

Constructing Scientific Authorities: Issue Framing of Chlorinated Disinfection Byproducts in Public Health

The practice of chlorine disinfection of drinking water to reduce microbial risks provides substantial benefits to public health. However, increasing concern around potential risks of cancer associated with exposure to chlorinated disinfection byproducts confuses this issue. This article examines the science agenda regarding chlorinated disinfection byproducts (CDBP) and cancer in Canada and the United States, focusing on the social construction of scientific knowledge claims and evidence. Data for this analysis were obtained from published documents as well as from in-depth interviews with epidemiologists and toxicologists centrally involved with the issue in both countries. Results of the analysis suggest that toxicological scientists want to close the door on the "chloroform issue" due to increasing evidence that chloroform is safe at low doses, because epidemiological scientists can no longer move forward the cancer science until significant improvements can be made in assessing human exposures, and because the scientific foci of research on DBP have shifted accordingly. Further, a distinction emerges in terms of how scientific uncertainties are interpreted when they cross-cut disciplines in the context of human health risk assessment. We suggest this tension reflects a balance of how uncertainty and authorities are managed in a mandated science-policy domain. Sufficient evidence was provided to keep the DBP issue on the regulatory agenda and to generate additional research, yet authorities and concomitant interpretations of uncertainty were contested. Such science generation and contestation inevitably influences complex risk assessment processes with respect to what water-related health risks are addressed and how.     

Health Risk Characterization of Chlorpyrifos Using Epidemiological Dose‐Response Data and Probabilistic Techniques: A Case Study with Rice Farmers in Vietnam

Various methods for risk characterization have been developed using probabilistic approaches. Data on Vietnamese farmers are available for the comparison of outcomes for risk characterization using different probabilistic methods. This article addresses the health risk characterization of chlorpyrifos using epidemiological dose‐response data and probabilistic techniques obtained from a case study with rice farmers in Vietnam. Urine samples were collected from farmers and analyzed for trichloropyridinol (TCP), which was converted into absorbed daily dose of chlorpyrifos. Adverse health response doses due to chlorpyrifos exposure were collected from epidemiological studies to develop dose‐adverse health response relationships. The health risk of chlorpyrifos was quantified using hazard quotient (HQ), Monte Carlo simulation (MCS), and overall risk probability (ORP) methods. With baseline (prior to pesticide spraying) and lifetime exposure levels (over a lifetime of pesticide spraying events), the HQ ranged from 0.06 to 7.1. The MCS method indicated less than 0.05% of the population would be affected while the ORP method indicated that less than 1.5% of the population would be adversely affected. With postapplication exposure levels, the HQ ranged from 1 to 32.5. The risk calculated by the MCS method was that 29% of the population would be affected, and the risk calculated by ORP method was 33%. The MCS and ORP methods have advantages in risk characterization due to use of the full distribution of data exposure as well as dose response, whereas HQ methods only used the exposure data distribution. These evaluations indicated that single‐event spraying is likely to have adverse effects on Vietnamese rice farmers.     

Bayesian Quantile Impairment Threshold Benchmark Dose Estimation for Continuous Endpoints

Quantitative risk assessment often begins with an estimate of the exposure or dose associated with a particular risk level from which exposure levels posing low risk to populations can be extrapolated. For continuous exposures, this value, the benchmark dose, is often defined by a specified increase (or decrease) from the median or mean response at no exposure. This method of calculating the benchmark dose does not take into account the response distribution and, consequently, cannot be interpreted based upon probability statements of the target population. We investigate quantile regression as an alternative to the use of the median or mean regression. By defining the dose–response quantile relationship and an impairment threshold, we specify a benchmark dose as the dose associated with a specified probability that the population will have a response equal to or more extreme than the specified impairment threshold. In addition, in an effort to minimize model uncertainty, we use Bayesian monotonic semiparametric regression to define the exposure–response quantile relationship, which gives the model flexibility to estimate the quantal dose–response function. We describe this methodology and apply it to both epidemiology and toxicology data.     

Ebola Virus Dose Response Model for Aerosolized Exposures: Insights from Primate Data

This study develops dose–response models for Ebolavirus using previously published data sets from the open literature. Two such articles were identified in which three different species of nonhuman primates were challenged by aerosolized Ebolavirus in order to study pathology and clinical disease progression. Dose groups were combined and pooled across each study in order to facilitate modeling. The endpoint of each experiment was death. The exponential and exact beta-Poisson models were fit to the data using maximum likelihood estimation. The exact beta-Poisson was deemed the recommended model because it more closely approximated the probability of response at low doses though both models provided a good fit. Although transmission is generally considered to be dominated by person-to-person contact, aerosolization is a possible route of exposure. If possible, this route of exposure could be particularly concerning for persons in occupational roles managing contaminated liquid wastes from patients being treated for Ebola infection and the wastewater community responsible for disinfection. Therefore, this study produces a necessary mathematical relationship between exposure dose and risk of death for the inhalation route of exposure that can support quantitative microbial risk assessment aimed at informing risk mitigation strategies including personal protection policies against occupational exposures.     

Regulating Under Uncertainty: Newsboy for Exposure Limits

Setting action levels or limits for health protection is complicated by uncertainty in the dose-response relation across a range of hazards and exposures. To address this issue, we consider the classic newsboy problem. The principles used to manage uncertainty for that case are applied to two stylized exposure examples, one for high dose and high dose rate radiation and the other for ammonia. Both incorporate expert judgment on uncertainty quantification in the dose-response relationship. The mathematical technique of probabilistic inversion also plays a key role. We propose a coupled approach, whereby scientists quantify the dose-response uncertainty using techniques such as structured expert judgment with performance weights and probabilistic inversion, and stakeholders quantify associated loss rates.