A Bayesian Semiparametric Model for Radiation Dose‐Response Estimation

In evaluating the risk of exposure to health hazards, characterizing the dose‐response relationship and estimating acceptable exposure levels are the primary goals. In analyses of health risks associated with exposure to ionizing radiation, while there is a clear agreement that moderate to high radiation doses cause harmful effects in humans, little has been known about the possible biological effects at low doses, for example, below 0.1 Gy, which is the dose range relevant to most radiation exposures of concern today. A conventional approach to radiation dose‐response estimation based on simple parametric forms, such as the linear nonthreshold model, can be misleading in evaluating the risk and, in particular, its uncertainty at low doses. As an alternative approach, we consider a Bayesian semiparametric model that has a connected piece‐wise‐linear dose‐response function with prior distributions having an autoregressive structure among the random slope coefficients defined over closely spaced dose categories. With a simulation study and application to analysis of cancer incidence data among Japanese atomic bomb survivors, we show that this approach can produce smooth and flexible dose‐response estimation while reasonably handling the risk uncertainty at low doses and elsewhere. With relatively few assumptions and modeling options to be made by the analyst, the method can be particularly useful in assessing risks associated with low‐dose radiation exposures.     

Improving Causal Inferences in Risk Analysis

Recent headlines and scientific articles projecting significant human health benefits from changes in exposures too often depend on unvalidated subjective expert judgments and modeling assumptions, especially about the causal interpretation of statistical associations. Some of these assessments are demonstrably biased toward false positives and inflated effects estimates. More objective, data‐driven methods of causal analysis are available to risk analysts. These can help to reduce bias and increase the credibility and realism of health effects risk assessments and causal claims. For example, quasi‐experimental designs and analysis allow alternative (noncausal) explanations for associations to be tested, and refuted if appropriate. Panel data studies examine empirical relations between changes in hypothesized causes and effects. Intervention and change‐point analyses identify effects (e.g., significant changes in health effects time series) and estimate their sizes. Granger causality tests, conditional independence tests, and counterfactual causality models test whether a hypothesized cause helps to predict its presumed effects, and quantify exposure‐specific contributions to response rates in differently exposed groups, even in the presence of confounders. Causal graph models let causal mechanistic hypotheses be tested and refined using biomarker data. These methods can potentially revolutionize the study of exposure‐induced health effects, helping to overcome pervasive false‐positive biases and move the health risk assessment scientific community toward more accurate assessments of the impacts of exposures and interventions on public health.     

Chronic Health Risks from Aggregate Exposures to Ionizing Radiation and Chemicals: Scientific Basis for an Assessment Framework

Very little quantitative analysis is currently available on the cumulative effects of exposure to multiple hazardous agents that have either similar or different mechanisms of action. Over the past several years, efforts have been made to develop the methodologies for risk assessment of chemical mixtures, but mixed exposures to two or more dissimilar agents such as radiation and one or more chemical agents have not yet been addressed in any substantive way. This article reviews the current understanding of the health risks arising from mixed exposures to ionizing radiation and specific chemicals. Specifically discussed is how mixed radiation/chemical exposures, when evaluated in aggregation, were linked to chronic health endpoints such as cancer and intermediate health outcomes such as chromosomal aberrations. Also considered is the extent to which the current practices are consistent with the scientific understanding of the health risks associated with mixed-agent exposures. From this the discussion moves to the research needs for assessing the cumulative health risks from aggregate exposures to ionizing radiation and chemicals. The evaluation indicates that essentially no guidance has been provided for conducting risk assessment for two agents with different mechanisms of action (i.e., energy deposition from ionizing radiation versus DNA interactions with chemicals) but similar biological endpoints (i.e., chromosomal aberrations, mutations, and cancer). The literature review also reveals the problems caused by the absence of both the basic science and an appropriate evaluation framework for the combined effects of mixed-agent exposures. This makes it difficult to determine whether there is truly no interaction or somehow the interaction is masked by the scale of effect observation or inappropriate dose-response assumptions.     

Air Toxics and Health Risks in California: The Public Health Implications of Outdoor Concentrations

Of the 188 hazardous air pollutants (HAPs) listed in the Clean Air Act, only a handful have information on human health effects, derived primarily from animal and occupational studies. Lack of consistent monitoring data on ambient air toxics makes it difficult to assess the extent of low-level, chronic, ambient exposures to HAPs that could affect human health, and limits attempts to prioritize and evaluate policy initiatives for emissions reduction. Modeled outdoor HAP concentration estimates from the U.S. Environmental Protection Agency's Cumulative Exposure Project were used to characterize the extent of the air toxics problem in California for the base year of 1990. These air toxics concentration estimates were used with chronic toxicity data to estimate cancer and noncancer hazards for individual HAPs and the risks posed by multiple pollutants. Although hazardous air pollutants are ubiquitous in the environment, potential cancer and noncancer health hazards posed by ambient exposures are geographically concentrated in three urbanized areas and in a few rural counties. This analysis estimated a median excess individual cancer risk of 2.7E-4 for all air toxics concentrations and 8600 excess lifetime cancer cases, 70% of which were attributable to four pollutants: polycyclic organic matter, 1,3 butadiene, formaldehyde, and benzene. For noncancer effects, the analysis estimated a total hazard index representing the combined effect of all HAPs considered. Each pollutant contributes to the index a ratio of estimated concentration to reference concentration. The median value of the index across census tracts was 17, due primarily to acrolein and chromium concentration estimates. On average, HAP concentrations and cancer and noncancer health risks originate mostly from area and mobile source emissions, although there are several locations in the state where point sources account for a large portion of estimated concentrations and health risks. Risk estimates from this study can provide guidance for prioritizing research, monitoring, and regulatory intervention activities to reduce potential hazards to the general population. Improved ambient monitoring efforts can help clarify uncertainties inherent in this analysis.     

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.     

Norovirus Dose–Response: Are Currently Available Data Informative Enough to Determine How Susceptible Humans Are to Infection from a Single Virus?

Two forms of single‐hit infection dose‐response models have previously been developed to assess available data from human feeding trials and estimate the norovirus dose‐response relationship. The mechanistic interpretations of these models include strong assumptions that warrant reconsideration: the first study includes an implicit assumption that there is no immunity to Norwalk virus among the specific study population, while the recent second study includes assumptions that such immunity could exist and that the nonimmune have no defensive barriers to prevent infection from exposure to just one virus. Both models addressed unmeasured virus aggregation in administered doses. In this work, the available data are reanalyzed using a generalization of the first model to explore these previous assumptions. It was hypothesized that concurrent estimation of an unmeasured degree of virus aggregation and important dose‐response parameters could lead to structural nonidentifiability of the model (i.e., that a diverse range of alternative mechanistic interpretations yield the same optimal fit), and this is demonstrated using the profile likelihood approach and by algebraic proof. It is also demonstrated that omission of an immunity parameter can artificially inflate the estimated degree of aggregation and falsely suggest high susceptibility among the nonimmune. The currently available data support the assumption of immunity within the specific study population, but provide only weak information about the degree of aggregation and susceptibility among the nonimmune. The probability of infection at low and moderate doses may be much lower than previously asserted, but more data from strategically designed dose‐response experiments are needed to provide adequate information.     

Bayesian Stacked Parametric Survival with Frailty Components and Interval‐Censored Failure Times: An Application to Food Allergy Risk

To better understand the risk of exposure to food allergens, food challenge studies are designed to slowly increase the dose of an allergen delivered to allergic individuals until an objective reaction occurs. These dose‐to‐failure studies are used to determine acceptable intake levels and are analyzed using parametric failure time models. Though these models can provide estimates of the survival curve and risk, their parametric form may misrepresent the survival function for doses of interest. Different models that describe the data similarly may produce different dose‐to‐failure estimates. Motivated by predictive inference, we developed a Bayesian approach to combine survival estimates based on posterior predictive stacking, where the weights are formed to maximize posterior predictive accuracy. The approach defines a model space that is much larger than traditional parametric failure time modeling approaches. In our case, we use the approach to include random effects accounting for frailty components. The methodology is investigated in simulation, and is used to estimate allergic population eliciting doses for multiple food allergens.     

Integrated Uncertainty Analysis for Ambient Pollutant Health Risk Assessment: A Case Study of Ozone Mortality Risk

The U.S. Environmental Protection Agency (EPA) uses health risk assessment to help inform its decisions in setting national ambient air quality standards (NAAQS). EPA's standard approach is to make epidemiologically‐based risk estimates based on a single statistical model selected from the scientific literature, called the “core” model. The uncertainty presented for “core” risk estimates reflects only the statistical uncertainty associated with that one model's concentration‐response function parameter estimate(s). However, epidemiologically‐based risk estimates are also subject to “model uncertainty,” which is a lack of knowledge about which of many plausible model specifications and data sets best reflects the true relationship between health and ambient pollutant concentrations. In 2002, a National Academies of Sciences (NAS) committee recommended that model uncertainty be integrated into EPA's standard risk analysis approach. This article discusses how model uncertainty can be taken into account with an integrated uncertainty analysis (IUA) of health risk estimates. It provides an illustrative numerical example based on risk of premature death from respiratory mortality due to long‐term exposures to ambient ozone, which is a health risk considered in the 2015 ozone NAAQS decision. This example demonstrates that use of IUA to quantitatively incorporate key model uncertainties into risk estimates produces a substantially altered understanding of the potential public health gain of a NAAQS policy decision, and that IUA can also produce more helpful insights to guide that decision, such as evidence of decreasing incremental health gains from progressive tightening of a NAAQS.     

Bounding Analysis of Drinking Water Health Risks from a Spill of Hydraulic Fracturing Flowback Water

A bounding risk assessment is presented that evaluates possible human health risk from a hypothetical scenario involving a 10,000‐gallon release of flowback water from horizontal fracturing of Marcellus Shale. The water is assumed to be spilled on the ground, infiltrates into groundwater that is a source of drinking water, and an adult and child located downgradient drink the groundwater. Key uncertainties in estimating risk are given explicit quantitative treatment using Monte Carlo analysis. Chemicals that contribute significantly to estimated health risks are identified, as are key uncertainties and variables to which risk estimates are sensitive. The results show that hypothetical exposure via drinking water impacted by chemicals in Marcellus Shale flowback water, assumed to be spilled onto the ground surface, results in predicted bounds between 10⁻¹⁰ and 10⁻⁶ (for both adult and child receptors) for excess lifetime cancer risk. Cumulative hazard indices (HI_CUMULATIVE) resulting from these hypothetical exposures have predicted bounds (5th to 95th percentile) between 0.02 and 35 for assumed adult receptors and 0.1 and 146 for assumed child receptors. Predicted health risks are dominated by noncancer endpoints related to ingestion of barium and lithium in impacted groundwater. Hazard indices above unity are largely related to exposure to lithium. Salinity taste thresholds are likely to be exceeded before drinking water exposures result in adverse health effects. The findings provide focus for policy discussions concerning flowback water risk management. They also indicate ways to improve the ability to estimate health risks from drinking water impacted by a flowback water spill (i.e., reducing uncertainty).     

Quantitative Assessment of Human MRSA Risks from Swine

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

An Approach for Modeling Noncancer Dose Responses with an Emphasis on Uncertainty

This paper presents an approach for characterizing the probability of adverse effects occurring in a population exposed to dose rates in excess of the Reference Dose (RfD). The approach uses a linear threshold (hockey stick) model of response and is based on the current system of uncertainty factors used in setting RfDs. The approach requires generally available toxicological estimates such as No-Observed-Adverse-Effect Levels (NOAELs) or Benchmark Doses and doses at which adverse effects are observed in 50% of the test animals (ED₅₀s). In this approach, Monte Carlo analysis is used to characterize the uncertainty in the dose response slope based on the range and magnitude of the key sources of uncertainty in setting protective doses. The method does not require information on the shape of the dose response curve for specific chemicals, but is amenable to the inclusion of such data. The approach is applied to four compounds to produce estimates of response rates for dose rates greater than the RfD     

Validation of a Novel Air Toxic Risk Model with Air Monitoring

Three modeling systems were used to estimate human health risks from air pollution: two versions of MNRiskS (for Minnesota Risk Screening), and the USEPA National Air Toxics Assessment (NATA). MNRiskS is a unique cumulative risk modeling system used to assess risks from multiple air toxics, sources, and pathways on a local to a state‐wide scale. In addition, ambient outdoor air monitoring data were available for estimation of risks and comparison with the modeled estimates of air concentrations. Highest air concentrations and estimated risks were generally found in the Minneapolis‐St. Paul metropolitan area and lowest risks in undeveloped rural areas. Emissions from mobile and area (nonpoint) sources created greater estimated risks than emissions from point sources. Highest cancer risks were via ingestion pathway exposures to dioxins and related compounds. Diesel particles, acrolein, and formaldehyde created the highest estimated inhalation health impacts. Model‐estimated air concentrations were generally highest for NATA and lowest for the AERMOD version of MNRiskS. This validation study showed reasonable agreement between available measurements and model predictions, although results varied among pollutants, and predictions were often lower than measurements. The results increased confidence in identifying pollutants, pathways, geographic areas, sources, and receptors of potential concern, and thus provide a basis for informing pollution reduction strategies and focusing efforts on specific pollutants (diesel particles, acrolein, and formaldehyde), geographic areas (urban centers), and source categories (nonpoint sources). The results heighten concerns about risks from food chain exposures to dioxins and PAHs. Risk estimates were sensitive to variations in methodologies for treating emissions, dispersion, deposition, exposure, and toxicity.     

Human Dose–Response Data for Francisella tularensis and a Dose‐ and Time‐Dependent Mathematical Model of Early‐Phase Fever Associated with Tularemia After Inhalation Exposure

Military health risk assessors, medical planners, operational planners, and defense system developers require knowledge of human responses to doses of biothreat agents to support force health protection and chemical, biological, radiological, nuclear (CBRN) defense missions. This article reviews extensive data from 118 human volunteers administered aerosols of the bacterial agent Francisella tularensis , strain Schu S4, which causes tularemia. The data set includes incidence of early‐phase febrile illness following administration of well‐characterized inhaled doses of F. tularensis . Supplemental data on human body temperature profiles over time available from de‐identified case reports is also presented. A unified, logically consistent model of early‐phase febrile illness is described as a lognormal dose–response function for febrile illness linked with a stochastic time profile of fever. Three parameters are estimated from the human data to describe the time profile: incubation period or onset time for fever; rise time of fever; and near‐maximum body temperature. Inhaled dose‐dependence and variability are characterized for each of the three parameters. These parameters enable a stochastic model for the response of an exposed population through incorporation of individual‐by‐individual variability by drawing random samples from the statistical distributions of these three parameters for each individual. This model provides risk assessors and medical decisionmakers reliable representations of the predicted health impacts of early‐phase febrile illness for as long as one week after aerosol exposures of human populations to F. tularensis .     

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.     

Development of a Screening Approach to Interpret Human Biomonitoring Data on Volatile Organic Compounds: Reverse Dosimetry on Biomonitoring Data for Trichloroethylene

A screening approach is developed for volatile organic compounds (VOCs) to estimate exposures that correspond to levels measured in fluids and/or tissues in human biomonitoring studies. The approach makes use of a generic physiologically-based pharmacokinetic (PBPK) model coupled with exposure pattern characterization, Monte Carlo analysis, and quantitative structure property relationships (QSPRs). QSPRs are used for VOCs with minimal data to develop chemical-specific parameters needed for the PBPK model. The PBPK model is capable of simulating VOC kinetics following multiple routes of exposure, such as oral exposure via water ingestion and inhalation exposure during shower events. Using published human biomonitoring data of trichloroethylene (TCE), the generic model is evaluated to determine how well it estimates TCE concentrations in blood based on the known drinking water concentrations. In addition, Monte Carlo analysis is conducted to characterize the impact of the following factors: (1) uncertainties in the QSPR-estimated chemical-specific parameters; (2) variability in physiological parameters; and (3) variability in exposure patterns. The results indicate that uncertainty in chemical-specific parameters makes only a minor contribution to the overall variability and uncertainty in the predicted TCE concentrations in blood. The model is used in a reverse dosimetry approach to derive estimates of TCE concentrations in drinking water based on given measurements of TCE in blood, for comparison to the U.S. EPA's Maximum Contaminant Level in drinking water. This example demonstrates how a reverse dosimetry approach can be used to facilitate interpretation of human biomonitoring data in a health risk context by deriving external exposures that are consistent with a biomonitoring data set, thereby permitting comparison with health-based exposure guidelines.     

Risk Characterization of Methyl tertiary Butyl Ether (MTBE) in Tap Water

Methyl tertiary butyl ether (MTBE) can enter surface water and groundwater through wet atmospheric deposition or as a result of fuel leaks and spills. About 30% of the U.S. population lives in areas where MTBE is in regular use. Ninety-five percent of this population is unlikely to be exposed to MTBE in tap water at concentrations exceeding 2 ppb, and most will be exposed to concentrations that are much lower and may be zero. About 5% of this population may be exposed to higher levels of MTBE in tap water, resulting from fuel tank leaks and spills into surface or groundwater used for potable water supplies. This paper describes the concentration ranges found and anticipated in surface and groundwater, and estimates the distribution of doses experienced by humans using water containing MTBE to drink, prepare food, and shower/bathe. The toxic properties (including potency) of MTBE when ingested, inhaled, and in contact with the skin are summarized. Using a range of human toxic potency values derived from animal studies, margins of exposure (MOE) associated with alternative chronic exposure scenarios are estimated to range from 1700 to 140,000. Maximum concentrations of MTBE in tap water anticipated not to cause adverse health effects are determined to range from 700 to 14,000 ppb. The results of this analysis demonstrate that no health risks are likely to be associated with chronic and subchronic human exposures to MTBE in tap water. Although some individuals may be exposed to very high concentrations of MTBE in tap water immediately following a localized spill, these exposures are likely to be brief in duration due to large-scale dilution and rapid volatilization of MTBE, the institution of emergency response and remediation measures to minimize human exposures, and the low taste and odor thresholds of MTBE which ensure that its presence in tap water is readily detected at concentrations well below the threshold for human injury.     

Potential Airborne Asbestos Exposure and Risk Associated with the Historical Use of Cosmetic Talcum Powder Products

Over time, concerns have been raised regarding the potential for human exposure and risk from asbestos in cosmetic‐talc–containing consumer products. In 1985, the U.S. Food and Drug Administration (FDA) conducted a risk assessment evaluating the potential inhalation asbestos exposure associated with the cosmetic talc consumer use scenario of powdering an infant during diapering, and found that risks were below levels associated with background asbestos exposures and risk. However, given the scope and age of the FDA's assessment, it was unknown whether the agency's conclusions remained relevant to current risk assessment practices, talc application scenarios, and exposure data. This analysis updates the previous FDA assessment by incorporating the current published exposure literature associated with consumer use of talcum powder and using the current U.S. Environmental Protection Agency's (EPA) nonoccupational asbestos risk assessment approach to estimate potential cumulative asbestos exposure and risk for four use scenarios: (1) infant exposure during diapering; (2) adult exposure from infant diapering; (3) adult exposure from face powdering; and (4) adult exposure from body powdering. The estimated range of cumulative asbestos exposure potential for all scenarios (assuming an asbestos content of 0.1%) ranged from 0.0000021 to 0.0096 f/cc‐yr and resulted in risk estimates that were within or below EPA's acceptable target risk levels. Consistent with the original FDA findings, exposure and corresponding health risk in this range were orders of magnitude below upper‐bound estimates of cumulative asbestos exposure and risk at ambient levels, which have not been associated with increased incidence of asbestos‐related disease.     

Diesel Emissions and Lung Cancer1

This paper uses two different methods to assess the potential risk of human lung cancer from exposure to diesel engine emissions. One method analyzes the best available epidemiological evidence on the lung cancer risks of persons exposed in their occupations to diesel engine emissions. The second conducts a comparative analysis of laboratory and epidemiological data on diesel engine emissions and two chemically related environmental exposures–coke oven emissions and roofing tar emissions. The estimates of potential risk derived from these two distinct methods are compared. The sources of uncertainty in each method are explicitly characterized. The value of these estimates for comparing the potential lung cancer risks from exposure to diesel engine emissions with other personal and societal risks are discussed. Also considered are the limitations of these results in predicting the possible excess incidence of lung cancer from ambient exposure to diesel emissions.     

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.