The application of non-default uncertainty factors in the U.S. EPA's Integrated Risk Information System (IRIS). Part I: UF(L), UF(S), and "other uncertainty factors"

Subchronic Exposure to Low-Level Lanthanum, Cerium, and Yttrium Mixtures Altered Cell Cycle and Increased Oxidative Stress Pathways in Human LO-2 Hepatocytes but Did Not Cause Malignant Transformation

Human exposures to rare earth elements are increasing with expanded use in aerospace, precision instruments, and new energy batteries, materials, and fertilizers. Individually these elements have low toxicity, although few investigations have examined the health effects of longer-term mixture exposures. We used the LO-2 cell line to examine the effects of graded exposures to lanthanum, cerium, and yttrium (LCY) mixtures at 1-, 100-, and 1000-fold their human background levels (0.31 μg/L La, 0.25 μg/L Ce, and 0.12 μg/L Y) on cell cycle, oxidative stress, and nuclear factor erythroid-2-related factor (NRF2) pathway biomarkers, assessing responses every 10 passages up to 100 passages. Cell migration, concanavalin A, malignant transformation, and tumorigenesis in nude mice were also examined. Mixed LCY exposures activated oxidative stress and the NRF2 pathway by the 30th passage and increased the proportion of cells in the S phase and cell cycle-specific biomarkers by the 40th passage. LCY exposures did not cause malignant transformation of hepatocytes or induced tumorigenesis in nude mice but enhanced cell proliferation, migration, and agglutination. Importantly, LCY mixtures with longer-term exposure activated the NRF2 pathway and altered the hepatocyte cell cycle at doses far below those used in previous toxicological studies. The consequences of LCY mixtures for public health merit further study.

Quantitative Characterization of Population-Wide Tissue- and Metabolite-Specific Variability in Perchloroethylene Toxicokinetics in Male Mice

Quantification of interindividual variability is a continuing challenge in risk assessment, particularly for compounds with complex metabolism and multi-organ toxicity. Toxicokinetic variability for perchloroethylene (perc) was previously characterized across 3 mouse strains and in 1 mouse strain with various degrees of liver steatosis. To further characterize the role of genetic variability in toxicokinetics of perc, we applied Bayesian population physiologically based pharmacokinetic (PBPK) modeling to the data on perc and metabolites in blood/plasma and tissues of male mice from 45 inbred strains from the Collaborative Cross (CC) mouse population. After identifying the most influential PBPK parameters based on global sensitivity analysis, we fit the model with a hierarchical Bayesian population analysis using Markov chain Monte Carlo simulation. We found that the data from 3 commonly used strains were not representative of the full range of variability in perc and metabolite blood/plasma and tissue concentrations across the CC population. Using interstrain variability as a surrogate for human interindividual variability, we calculated dose-dependent, chemical-, and tissue-specific toxicokinetic variability factors (TKVFs) as candidate science-based replacements for the default uncertainty factor for human toxicokinetic variability of 100.5. We found that toxicokinetic variability factors for glutathione conjugation metabolites of perc showed the greatest variability, often exceeding the default, whereas those for oxidative metabolites and perc itself were generally less than the default. Overall, we demonstrate how a combination of a population-based mouse model such as the CC with Bayesian population PBPK modeling can reduce uncertainty in human toxicokinetic variability and increase accuracy and precision in quantitative risk assessment.

Quantifying Uncertainty in Ecotoxicological Risk Assessment: MUST, a Modular Uncertainty Scoring Tool

Whether conducting a risk, hazard, or alternatives assessment, one invariably struggles with the task of reconciling multiple available values of toxicological thresholds into a single outcome. When combining multiple pieces of evidence from many different sources, it is important to consider the role of data uncertainty. Uncertainty is inherent to all scientific data. However, in toxicological assessments, controversies and uncertainties are typically understated; they lack methodological transparency; or they poorly integrate qualitative and quantitative sources of information. Similarly, in model development, data curation is rarely performed with sufficient rigor, particularly when applying big data statistics. To overcome the hurdles of a decision process that must reconcile divergent data, we developed an uncertainty scoring tool that can be trained to reproduce specific decision-making paradigms and ensure consistency in the practitioner's judgment across complex scenarios. While designed to aid with ecotoxicological assessments and predictive model development, the tool's applicability extends to any decision-making process that calls for synthesis of incongruent data. Here, we highlight the development process, as well as demonstrate the method's utility in several prototypical ecotoxicological case studies.

Incorporation of the glutathione conjugation pathway in an updated physiologically-based pharmacokinetic model for perchloroethylene in mice

BACKGROUND

Perchloroethylene (perc) induced target organ toxicity has been associated with tissue-specific metabolic pathways. Previous physiologically-based pharmacokinetic (PBPK) modeling of perc accurately predicted oxidative metabolites but suggested the need to better characterize glutathione (GSH) conjugation as well as toxicokinetic uncertainty and variability.

OBJECTIVES

We updated the previously published "harmonized" perc PBPK model in mice to better characterize GSH conjugation metabolism as well as the uncertainty and variability of perc toxicokinetics.

METHODS

The updated PBPK model includes expanded models for perc and its oxidative metabolite trichloroacetic acid (TCA), and physiologically-based sub-models for conjugative metabolites. Previously compiled mouse kinetic data in B6C3F1 and Swiss-Webster mice were augmented to include data from a recent study in male C57BL/6J mice that measured perc and metabolites in serum and multiple tissues. Hierarchical Bayesian population analysis using Markov chain Monte Carlo was conducted to characterize uncertainty and inter-strain variability in perc metabolism.

RESULTS

The updated model fit the data as well or better than the previously published "harmonized" PBPK model. Tissue dosimetry for both oxidative and conjugative metabolites was successfully predicted across the three strains of mice, with estimated residuals errors of 2-fold for majority of data. Inter-strain variability across three strains was evident for oxidative metabolism; GSH conjugation data were only available for one strain.

CONCLUSIONS

This updated PBPK model fills a critical data gap in quantitative risk assessment by predicting the internal dosimetry of perc and its oxidative and GSH conjugation metabolites and lays the groundwork for future studies to better characterize toxicokinetic variability.

Incorporating epigenetic data into the risk assessment process for the toxic metals arsenic, cadmium, chromium, lead, and mercury: strategies and challenges

Exposure to toxic metals poses a serious human health hazard based on ubiquitous environmental presence, the extent of exposure, and the toxicity and disease states associated with exposure. This global health issue warrants accurate and reliable models derived from the risk assessment process to predict disease risk in populations. There has been considerable interest recently in the impact of environmental toxicants such as toxic metals on the epigenome. Epigenetic modifications are alterations to an individual's genome without a change in the DNA sequence, and include, but are not limited to, three commonly studied alterations: DNA methylation, histone modification, and non-coding RNA expression. Given the role of epigenetic alterations in regulating gene and thus protein expression, there is the potential for the integration of toxic metal-induced epigenetic alterations as informative factors in the risk assessment process. In the present review, epigenetic alterations induced by five high priority toxic metals/metalloids are prioritized for analysis and their possible inclusion into the risk assessment process is discussed.

Physiologically based pharmacokinetic (PBPK) modeling of interstrain variability in trichloroethylene metabolism in the mouse

BACKGROUND

Quantitative estimation of toxicokinetic variability in the human population is a persistent challenge in risk assessment of environmental chemicals. Traditionally, interindividual differences in the population are accounted for by default assumptions or, in rare cases, are based on human toxicokinetic data.

OBJECTIVES

We evaluated the utility of genetically diverse mouse strains for estimating toxicokinetic population variability for risk assessment, using trichloroethylene (TCE) metabolism as a case study.

METHODS

We used data on oxidative and glutathione conjugation metabolism of TCE in 16 inbred and 1 hybrid mouse strains to calibrate and extend existing physiologically based pharmacokinetic (PBPK) models. We added one-compartment models for glutathione metabolites and a two-compartment model for dichloroacetic acid (DCA). We used a Bayesian population analysis of interstrain variability to quantify variability in TCE metabolism.

RESULTS

Concentration-time profiles for TCE metabolism to oxidative and glutathione conjugation metabolites varied across strains. Median predictions for the metabolic flux through oxidation were less variable (5-fold range) than that through glutathione conjugation (10-fold range). For oxidative metabolites, median predictions of trichloroacetic acid production were less variable (2-fold range) than DCA production (5-fold range), although the uncertainty bounds for DCA exceeded the predicted variability.

CONCLUSIONS

Population PBPK modeling of genetically diverse mouse strains can provide useful quantitative estimates of toxicokinetic population variability. When extrapolated to lower doses more relevant to environmental exposures, mouse population-derived variability estimates for TCE metabolism closely matched population variability estimates previously derived from human toxicokinetic studies with TCE, highlighting the utility of mouse interstrain metabolism studies for addressing toxicokinetic variability.

Analysis of worldwide Regulatory Guidance Values for less frequently regulated elemental surface soil contaminants

Guidance values are used to regulate exposure to surface soil contamination. In the United States, element guidance values have been promulgated by at least 6 federal agencies, 46 states, and several regional, city, county, territorial, and autonomous Native American jurisdictions. Guidance values have also been promulgated in at least 74 other United Nations member states. A companion manuscript examined the values applied to the eight most frequently regulated elements (Pb, Cd, As, Ni, Cr, Hg, Cu, and Zn) each for which there are at least 300 guidance values. This manuscript extends analysis to the second tier of contaminants (Sb, Ba, Be, Co, F, Mn, Mo, Se, Ag, Tl, Sn, and V) each for which there are at least 100 guidance values. These values span from 3.4 (for Sn) to 6.1 (for Be) orders of magnitude. Their distributions resemble those of lognormal random variables, but also contain non-random value clusters. On average, with the exception of cobalt, the values used in the U.S. are higher (less conservative) than those used elsewhere. Only about 44% of U.S. values and 31% of all values fall within uncertainty bounds computed for the U.S. Environmental Protection Agency health risk model applied to the elements considered.

Prediction of acute mammalian toxicity using QSAR methods: a case study of sulfur mustard and its breakdown products

Predicting toxicity quantitatively, using Quantitative Structure Activity Relationships (QSAR), has matured over recent years to the point that the predictions can be used to help identify missing comparison values in a substance's database. In this manuscript we investigate using the lethal dose that kills fifty percent of a test population (LD₅₀) for determining relative toxicity of a number of substances. In general, the smaller the LD₅₀ value, the more toxic the chemical, and the larger the LD₅₀ value, the lower the toxicity. When systemic toxicity and other specific toxicity data are unavailable for the chemical(s) of interest, during emergency responses, LD₅₀ values may be employed to determine the relative toxicity of a series of chemicals. In the present study, a group of chemical warfare agents and their breakdown products have been evaluated using four available rat oral QSAR LD₅₀ models. The QSAR analysis shows that the breakdown products of Sulfur Mustard (HD) are predicted to be less toxic than the parent compound as well as other known breakdown products that have known toxicities. The QSAR estimated break down products LD₅₀ values ranged from 299 mg/kg to 5,764 mg/kg. This evaluation allows for the ranking and toxicity estimation of compounds for which little toxicity information existed; thus leading to better risk decision making in the field.

Worldwide regulatory guidance values for surface soil exposure to noncarcinogenic polycyclic aromatic hydrocarbons

Regulatory guidance values are used worldwide to control residential exposures to surface soil contamination. A total of 1791 values used in 39 nations to control exposures to eight polycyclic aromatic hydrocarbons (PAH) not normally considered to be carcinogenic (acenaphthene, acenaphthylene, anthracene, benzo(g,h,i)perylene, fluroanthene, fluorene, phenanthrene and pyrene) are examined. The guidance values vary over 5.8 to 8.3 orders of magnitude and have distributions that strongly resemble distributions of lognormal random variables. Where the U.S. Environmental Protection Agency (USEPA) has provided values, they fall near the upper end of the distributions and appear to exert influence on values from other U.S. jurisdictions. Approximately 84% of all U.S. values fall above a median PAH value. Uncertainty bounds approximated from the USEPA risk model encompass 28% of the regulatory values and these are predominantly from U.S. jurisdictions. An unusually high degree of toxicology data uncertainty for these PAH appears to be a significant factor influencing variability. The USEPA assigns its highest toxicology uncertainty factor value of 3000 to these components.

Prenatal developmental toxicity of decabromodiphenyl ethane in the rat and rabbit

BACKGROUND

The potential embryotoxic and teratogenic effects of decabromodiphenyl ethane (DBDPEthane; CASRN 84852-53-9) were evaluated in prenatal developmental studies using rats and rabbits and performed in accordance with international guidelines and Good Laboratory Practice standards. Preliminary dose-range-finding studies were conducted, which indicated doses up to 1,250 mg/kg-day were well tolerated by both rats and rabbits.

METHODS

For the developmental studies, animals were administered DBDPEthane via gavage at dosage levels of 0, 125, 400, or 1,250 mg/kg-day from gestation day (GD) 6 through 15 for rats and GDs 6 through 18 for rabbits. All female rats and rabbits were sacrificed on GD 20 or GD 29, respectively, and subjected to cesarean section. Fetuses were individually weighed, sexed, and examined for external, visceral and skeletal abnormalities.

RESULTS

No treatment-related mortality, abortions, or clinical signs of toxicity were observed during the study. Body weights, body weight gain, and food consumption were not affected by treatment. No significant internal abnormalities were observed in either species on necropsy. Cesarean section parameters were comparable between control and treated groups. No treatment-induced malformations or developmental variations occurred.

CONCLUSIONS

Based on these results, no evidence of maternal toxicity, developmental toxicity, or teratogenicity was observed in rats or rabbits treated with DBDPEthane at dosage levels up to 1,250 mg/kg-day.